  This paper develops an individual-based stochastic network SIR model for the
empirical analysis of the Covid-19 pandemic. It derives moment conditions for
the number of infected and active cases for single as well as multigroup
epidemic models. These moment conditions are used to investigate the
identification and estimation of the transmission rates. The paper then
proposes a method that jointly estimates the transmission rate and the
magnitude of under-reporting of infected cases. Empirical evidence on six
European countries matches the simulated outcomes once the under-reporting of
infected cases is addressed. It is estimated that the number of actual cases
could be between 4 to 10 times higher than the reported numbers in October 2020
and declined to 2 to 3 times in April 2021. The calibrated models are used in
the counterfactual analyses of the impact of social distancing and vaccination
on the epidemic evolution, and the timing of early interventions in the UK and
Germany.


  In this article we compare two multi-patch models for the spread of Zika
virus based on an SIRUV model. When the commuting between patches is ceased we
expect that all the patches follow the dynamics of the single patch model. We
show in an example that the effective population size should be used rather
than the population size of the respective patch.


  Inferring functional relationships within complex networks from static
snapshots of a subset of variables is a ubiquitous problem in science. For
example, a key challenge of systems biology is to translate cellular
heterogeneity data obtained from single-cell sequencing or flow-cytometry
experiments into regulatory dynamics. We show how static population snapshots
of co-variability can be exploited to rigorously infer properties of gene
expression dynamics when gene expression reporters probe their upstream
dynamics on separate time-scales. This can be experimentally exploited in
dual-reporter experiments with fluorescent proteins of unequal maturation
times, thus turning an experimental bug into an analysis feature. We derive
correlation conditions that detect the presence of closed-loop feedback
regulation in gene regulatory networks. Furthermore, we show how genes with
cell-cycle dependent transcription rates can be identified from the variability
of co-regulated fluorescent proteins. Similar correlation constraints might
prove useful in other areas of science in which static correlation snapshots
are used to infer causal connections between dynamically interacting
components.


  Branching processes are widely used to model evolutionary and population
dynamics as well as the spread of infectious diseases. To characterize the
dynamics of their growth or spread, the basic reproduction number $R_0$ has
received considerable attention. In the context of infectious diseases, it is
usually defined as the expected number of secondary cases produced by an
infectious case in a completely susceptible population. Typically $R_0>1$
indicates that an outbreak is expected to continue and to grow exponentially,
while $R_0<1$ usually indicates that an outbreak is expected to terminate after
some time.
  In this work, we show that fluctuations of the dynamics in time can lead to a
continuation of outbreaks even when the expected number of secondary cases from
a single case is below $1$. Such fluctuations are usually neglected in
modelling of infectious diseases by a set of ordinary differential equations,
such as the classic SIR model. We showcase three examples: 1) extinction
following an Ornstein-Uhlenbeck process, 2) extinction switching randomly
between two values and 3) mixing of two populations with different $R_0$
values. We corroborate our analytical findings with computer simulations.


  Mathematical models have come to play a key role in global pandemic
preparedness and outbreak response: helping to plan for disease burden,
hospital capacity, and inform non-pharmaceutical interventions. Such models
have played a pivotal role in the COVID-19 pandemic, with transmission models
and, by consequence, modelers guiding global, national, and local responses to
SARS-CoV-2. However, these models have systematically failed to account for the
social and structural factors which lead to socioeconomic, racial, and
geographic health disparities. Why do epidemiologic models of emerging
infections ignore known structural drivers of disparate health outcomes? What
have been the consequences of this limited framework? What should be done to
develop a more holistic approach to modeling-based decision-making during
pandemics? In this Perspective, we evaluate potential historical and political
explanations for the exclusion of drivers of disparity in infectious disease
models for emerging infections, which have often been characterized as "equal
opportunity infectors" despite ample evidence to the contrary. We look to
examples from other disease systems (HIV, STIs) as a potential blueprint for
how social connections, environmental, and structural factors can be integrated
into a coherent, rigorous, and interpretable modeling framework. We conclude by
outlining principles to guide modeling of emerging infections in ways that
represent the causes of inequity in infection as central rather than peripheral
mechanisms.


  Microbes are essentially yet convolutedly linked with human lives on the
earth. They critically interfere in different physiological processes and thus
influence overall health status. Studying microbial species is used to be
constrained to those that can be cultured in the lab. But it excluded a huge
portion of the microbiome that could not survive on lab conditions. In the past
few years, the culture-independent metagenomic sequencing enabled us to explore
the complex microbial community coexisting within and on us. Metagenomics has
equipped us with new avenues of investigating the microbiome, from studying a
single species to a complex community in a dynamic ecosystem. Thus, identifying
the involved microbes and their genomes becomes one of the core tasks in
metagenomic sequencing. Metagenome-assembled genomes are groups of contigs with
similar sequence characteristics from de novo assembly and could represent the
microbial genomes from metagenomic sequencing. In this paper, we reviewed a
spectrum of tools for producing and annotating metagenome-assembled genomes
from metagenomic sequencing data and discussed their technical and biological
perspectives.


  Phylostratigraphy suggests that new genes are continually born de novo from
non-genic sequences, and the genes that persist found new lineages,
contributing to the adaptive evolution of organisms. While recent evidence
supports the view that de novo gene birth is frequent and widespread, the
mechanisms underlying this process are yet to be discovered. Here we
hypothesize and examine a potential general mechanism of gene birth driven by
the accumulation of beneficial mutations at non-genic loci. To demonstrate this
possibility, we model this mechanism within the boundaries set by current
knowledge on mutation effects. Estimates from this analysis are in line with
observations of recurrent and extensive gene birth in genomics studies. Thus,
we propose that, rather than being inactive and silent, non-genic regions are
likely to be dynamic storehouses of potential genes.


  Epilepsy is a disorder of the nervous system that can affect people of any
age group. With roughly 50 million people worldwide diagnosed with the
disorder, it is one of the most common neurological disorders. The EEG is an
indispensable tool for diagnosis of epileptic seizures in an ideal case, as
brain waves from an epileptic person will present distinct abnormalities.
However, in real world situations there will often be biological and electrical
noise interference, as well as the issue of a multichannel signal, which
introduce a great challenge for seizure detection. For this study, the Temple
University Hospital (TUH) EEG Seizure Corpus dataset was used. This paper
proposes a novel channel selection method which isolates different frequency
ranges within five channels. This is based upon the frequencies at which normal
brain waveforms exhibit. A one second window was selected, with a 0.5 second
overlap. Wavelet signal denoising was performed using Daubechies 4 wavelet
decomposition, thresholding was applied using minimax soft thresholding
criteria. Filter banking was used to localise frequency ranges from five
specific channels. Statistical features were then derived from the outputs.
After performing bagged tree classification using 500 learners, a test accuracy
of 0.82 was achieved.


  Current questions in ecology revolve around instabilities in the dynamics on
spatial networks and particularly the effect of node heterogeneity. We extend
the Master Stability Function formalism to inhomogeneous biregular networks
having two types of spatial nodes. Notably, this class of systems also allows
the investigation of certain types of dynamics on higher-order networks.
Combined with the Generalized Modelling approach to study the linear stability
of steady states, this is a powerful tool to numerically asses the stability of
large ensembles of systems. We analyze the stability of ecological
metacommunities with two distinct types of habitats analytically and
numerically in order to identify several sets of conditions under which the
dynamics can become stabilized by dispersal. Our analytical approach allows
general insights into stabilizing and destabilizing effects in metapopulations.
Specifically, we show that maladaptive dispersal may be stable under certain
conditions.


  The Covid-19 pandemic has taken millions of lives, demonstrating the tragedy
and disruption of respiratory diseases, and how difficult they can be to
manage. However, there is still significant debate in the scientific community
as to which transmission pathways are most significant and how settings and
behaviour affect risk of infection, which all have implications for which
mitigation strategies are most effective. This study presents a general model
to estimate the rate of viral transfer between individuals, objects, and the
air. The risk of infection to individuals in a setting is then computed
considering the behaviour and interactions of individuals between themselves
and the environment in the setting, survival times of the virus on different
surface types and in the air, and mitigating interventions (ventilation, hand
disinfection, surface cleaning, etc.). The model includes discrete events such
as touch events, individuals entering/leaving the setting, and cleaning events.
  We demonstrate the model capabilities on three case studies to quantify and
understand the relative risk associated with the different transmission
pathways and the effectiveness of mitigation strategies in different settings.
The results show the importance of considering all transmission pathways and
their interactions, with each scenario displaying different dominant pathways
depending on the setting and behaviours of individuals therein. The flexible
model, which is freely available, can be used to quickly simulate the spread of
any respiratory virus via the modelled transmission pathways and the efficacy
of potential mitigation strategies in any enclosed setting by making reasonable
assumptions regarding the behaviour of its occupants. It is hoped that the
model can be used to inform sensible decision-making regarding viral infection
mitigations that are targeted to specific settings and pathogens.


  We prove that nested canalizing functions are the minimum-sensitivity Boolean
functions for any given activity ratio and we characterize the sensitivity
boundary which has a nontrivial fractal structure. We further observe, on an
extensive database of regulatory functions curated from the literature, that
this bound severely constrains the robustness of biological networks. Our
findings suggest that the accumulation near the "edge of chaos" in these
systems is a natural consequence of a drive towards maximum stability while
maintaining plasticity in transcriptional activity.


  Autocatalysis underlies the ability of chemical and biochemical systems to
replicate. Recently, Blokhuis et al. gave a stoechiometric definition of
autocatalysis for reaction networks, stating the existence of a combination of
reactions such that the balance for all autocatalytic species is strictly
positive, and investigated minimal autocatalytic networks, called {\em
autocatalytic cores}. By contrast, spontaneous autocatalysis -- namely,
exponential amplification of all species internal to a reaction network,
starting from a diluted regime, i.e. low concentrations -- is a dynamical
property.
  We introduce here a topological condition (Top) for autocatalysis, namely:
restricting the reaction network description to highly diluted species, we
assume existence of a strongly connected component possessing at least one
reaction with multiple products (including multiple copies of a single
species). We find this condition to be necessary and sufficient for
stoechiometric autocatalysis. When degradation reactions have small enough
rates, the topological condition further ensures dynamical autocatalysis,
characterized by a strictly positive Lyapunov exponent giving the instantaneous
exponential growth rate of the system.
  The proof is generally based on the study of auxiliary Markov chains. We
provide as examples general autocatalytic cores of Type I and Type III in the
typology of Blokhuis et al. In a companion article, Lyapunov exponents and the
behavior in the growth regime are studied quantitatively beyond the present
diluted regime .


  We study the dynamics about equilibria of an infinite dimension
coagulation-fragmentation-death model for the silicosis disease mechanism
introduced recently by da Costa, Drmota, and Grinfeld [Modelling silicosis:
structure of equilibria, Euro. J. Appl. Math., 31 (6), (2020) 950-967] in the
case where the rate coefficients are piecewise constant.


  A principle of evolutionary adaptation is applied to the Lotka--Volterra
models, in particular to the food webs. We present a relatively simple
computational algorithm of optimization with respect to a given criterion. This
algorithm boils down to a sequence of easy to solve linear programming
problems. As a criterion for the optimization we use the total weighted
population size of the given community and an ecological fitness, which is an
analogue of the potential energy in physics. We show by computational
experiments that it is almost always possible to substantially increase the
total weighed population size for an especially simple food web -- food chain;
we also show that food chains are evolutionary unstable under the given
optimization criteria and, if allowed, evolve into more complicated structures
of food webs.


  Prediction of protein-ligand binding affinity is a major goal in drug
discovery. Generally, free energy gap is calculated between two states (e.g.,
ligand binding and unbinding). The energy gap implicitly includes the effects
of changes in protein dynamics induced by the binding ligand. However, the
relationship between protein dynamics and binding affinity remains unclear.
Here, we propose a novel method that represents protein behavioral change upon
ligand binding with a simple feature that can be used to predict protein-ligand
affinity. From unbiased molecular simulation data, an unsupervised deep
learning method measures the differences in protein dynamics at a
ligand-binding site depending on the bound ligands. A dimension-reduction
method extracts a dynamic feature that is strongly correlated to the binding
affinities. Moreover, the residues that play important roles in protein-ligand
interactions are specified based on their contribution to the differences.
These results indicate the potential for dynamics-based drug discovery.


  Multilabel learning is an important topic in machine learning research.
Evaluating models in multilabel settings requires specific cross validation
methods designed for multilabel data. In this article, we show that the most
widely used cross validation split quality measures do not behave adequately
with multilabel data that has strong class imbalance. We present improved
measures and an algorithm, optisplit, for optimising cross validations splits.
We present an extensive comparison of various types of cross validation methods
in which we show that optisplit produces more even cross validation splits than
the existing methods and that it is fast enough to be used on big Gene Ontology
(GO) datasets.


  We review epidemiological models for the propagation of the COVID-19 pandemic
during the early months of the outbreak: from February to May 2020. The aim is
to propose a methodological review that highlights the following
characteristics: (i) the epidemic propagation models, (ii) the modeling of
intervention strategies, (iii) the models and estimation procedures of the
epidemic parameters and (iv) the characteristics of the data used. We finally
selected 80 articles from open access databases based on criteria such as the
theoretical background, the reproducibility, the incorporation of interventions
strategies, etc. It mainly resulted to phenomenological, compartmental and
individual-level models. A digital companion including an online sheet, a
Kibana interface and a markdown document is proposed. Finally, this work
provides an opportunity to witness how the scientific community reacted to this
unique situation.


  In this study, we present an epidemic-controlled SIRD model with two types of
control strategies: mask wear and screening. The aim of this study is to
minimize the number of deceased keeping a minimal cost of mask advertising and
screening. The model is proved to be well-posed and to have an invariant
region. Also, a thorough study of the dynamics is effected and the basic
reproduction number is used to study the stability of the steady states of the
model. As for the optimal control analysis, the existence of an optimal control
was checked. Then its characterization was carried out using Pontryagin's
minimum principle. Numerical simulations are conducted after that with
different values of maximal screening for comparison. The findings of the
optimal control analysis and numerical simulations both reveal that the optimal
pair of strategies contribute enormously in lowering the number of infected and
dead individuals. Although zero infection is not achieved in the population,
this study implies that carrying an optimal approach constitutes a major step
in controlling the spread of the disease to the barest minimum that can buy
time for middle and low-income countries to carry on with their vaccination
strategies.


  We introduce a novel compartmental model accounting for the effects of
vaccine efficacy, deployment rates and timing of initiation of deployment. We
simulate different scenarios and initial conditions, and we find that higher
abundancy and rate of deployment of low efficacy vaccines lowers the cumulative
number of deaths in comparison to slower deployment of high efficacy vaccines.
We also forecast that, at the same daily deployment rate, the earlier
introduction of vaccination schemes with lower efficacy would also lower the
number of deaths with respect to a delayed introduction of high efficacy
vaccines, which can however, still achieve lower numbers of infections and
better herd immunity.


  The last decade has seen increasing theoretical and empirical support for
beneficial reversals of dominance, which reduce the costs of genetic trade-offs
and enable selection to maintain genetic variation for fitness. Beneficial
dominance reversals are characterized by the beneficial allele for a given
context (e.g. habitat, developmental stage, trait, or sex) being dominant in
that context but recessive where deleterious. This context-dependence at least
partially mitigates the fitness consequence of heterozygotes carrying one
non-beneficial allele for their context and can result in balancing selection
that maintains alternative alleles. Dominance reversals are theoretically
plausible and are supported by mounting empirical evidence. Here we highlight
the importance of beneficial dominance reversals as a mechanism for the
mitigation of genetic conflict and review the theory and empirical evidence for
them. We identify some areas in need of further research and development and
outline three methods that leverage signatures of dominance reversal toward
identifying antagonistic genetic variation. There is ample scope for the
development of new empirical methods as well as reanalysis of existing data
through the lens of dominance reversals. A greater focus on this topic will
expand our understanding of the mechanisms that resolve genetic conflict and
whether they maintain genetic variation.


  We consider biochemical systems associated with a generalised class of Petri
nets with possibly negative token numbers. We show that the existence of a
structural polyhedral Lyapunov function for the biochemical system is
equivalent to the boundedness of the associated Petri net evolution or,
equivalently, to the finiteness of the number of states reachable from each
initial condition. For networks that do not admit a polyhedral Lyapunov
function, we investigate whether it is possible to enforce polyhedral
structural stability by applying a strong negative feedback on some pinned
nodes: in terms of the Petri net, this is equivalent to turning pinned nodes
into black holes that clear any positive or negative incoming token. If such
nodes are chosen so that the transformed Petri net has bounded discrete
trajectories, then there exists a stabilising pinning control: the biochemical
network becomes Lyapunov stable if a sufficiently strong local negative
feedback is applied to the pinned nodes. These results allow us to structurally
identify the critical nodes to be locally controlled so as to ensure the
stability of the whole network.


  Compartmental models have long served as important tools in mathematical
epidemiology, with their usefulness highlighted by the recent COVID-19
pandemic. However, most of the classical models fail to account for certain
features of this disease and others like it, such as the ability of exposed
individuals to recover without becoming infectious, or the possibility that
asymptomatic individuals can indeed transmit the disease but at a lesser rate
than the symptomatic.
  In the first part of this paper we propose two new compartmental
epidemiological models and study their equilibria, obtaining an endemic
threshold theorem for the first model. In the second part of the paper, we
treat the second model as an affine control system with two controls:
vaccination and mitigation. We show that this system is static feedback
linearizable, present some simulations, and investigate of an optimal control
version of the problem. We conclude with some open problems and ideas for
future research.


  A tiny bacterial flagellar motor rotates in both counter-clockwise (CCW) and
clockwise (CW) rotational directions. The most important measurable
characteristic of the flagellar motor is its torque versus angular speed
relationship in CCW or CW modes, which is found to be non-symmetrical with each
other, and still, such a phenomenon is not clearly understood.Here, we explain
this asymmetry through a mechanistic model based on the detailed torque
analysis for the rotation of the motor and the revolutionary as well as
spinning motion of the filament and bead. We find out that the asymmetry
results from the conformational changing of the hook due to rotational
switching, rather than any non-symmetric changes in the potential of mean force
generated by the stator-rotor interactions. In CCW mode, when the hook remains
bend and flexible, the revolution motion predominates and the restoring torque
in this motion, originated due to drag, governs the shape of the torque-speed
curve. However, in CW mode, spinning motion dominates as the hook becomes
straight and rigid, and the linear torque-speed relation arises due to the
restoring torque for the drag corresponding to this motion. Our study indicates
the significant role of the hook's conformational change upon the biological
functions of the motor and paves the way for further experimental exploration
on the structural origin of such asymmetry.


  The occurrence and distributions of wildlife populations and communities are
shifting as a result of global changes. To evaluate whether these shifts are
negatively impacting biodiversity processes, it is critical to monitor the
status, trends, and effects of environmental variables on entire communities.
However, modeling the dynamics of multiple species simultaneously can require
large amounts of diverse data, and few modeling approaches exist to
simultaneously provide species and community level inferences. We present an
"integrated community occupancy model" (ICOM) that unites principles of data
integration and hierarchical community modeling in a single framework to
provide inferences on species-specific and community occurrence dynamics using
multiple data sources. We use simulations to compare the ICOM to previously
developed hierarchical community occupancy models and single species integrated
distribution models. We then apply our model to assess the occurrence and
biodiversity dynamics of foliage-gleaning birds in the White Mountain National
Forest in the northeastern USA from 2010-2018 using three independent data
sources. Simulations reveal that integrating multiple data sources in the ICOM
increased precision and accuracy of species and community level inferences
compared to single data source models, although benefits of integration were
dependent on data source quality (e.g., amount of replication). Compared to
single species models, the ICOM yielded more precise species-level estimates.
Within our case study, the ICOM had the highest out-of-sample predictive
performance compared to single species models and models that used only a
subset of the three data sources. The ICOM offers an attractive approach to
estimate species and biodiversity dynamics, which is additionally valuable to
inform management objectives of both individual species and their broader
communities.


  Spatiotemporal processes have the potential to be one of the most influential
factors governing how fisheries targeting sedentary species respond to
harvesting. Despite this, management strategy evaluation often fails to account
for space or does so at low resolutions due to compute constraints. In this
paper, a multi-species spatiotemporal size-structured operating model for
sedentary species is presented. The model combines a spatially continuous
Gaussian Markov Random Field model of the population dynamics with an areal
harvesting model that supports preferential targeting and site selection
constraints (e.g., economic constraints). This approach is very compute
efficient, which makes it feasible to simulate realistic fisher dynamics and
catch data at true spatial scale (e.g., the swept area of a dredge). The New
Zealand surfclam fishery was used as a case study to demonstrate the
versatility of the operating model and to showcase the simulation of localized
depletion, which was manifest in the generation of realistic
catch-per-unit-effort data that were uncorrelated with the trends in population
abundance. The model is available as part of the open-source R package
spatialSim.


  Herewith we discuss a network model of the ferroptosis avascular and vascular
tumor growth based on our previous proposed framework. Chiefly, ferroptosis
should be viewed as a first order phase transition characterized by a
supercritical Andronov Hopf bifurcation, with the emergence of limit cycle. The
increase of the population of the oxidized PUFA fragments, take as the control
parameter, involves an inverse Feigenbaum, (a cascade of saddle foci
Shilnikov's bifurcations) scenario, which results in the stabilization of the
dynamics and in a decrease of complexity.


  As the second wave in India mitigates, COVID-19 has now infected about 29
million patients countrywide, leading to more than 350 thousand people dead. As
the infections surged, the strain on the medical infrastructure in the country
became apparent. While the country vaccinates its population, opening up the
economy may lead to an increase in infection rates. In this scenario, it is
essential to effectively utilize the limited hospital resources by an informed
patient triaging system based on clinical parameters. Here, we present two
interpretable machine learning models predicting the clinical outcomes,
severity, and mortality, of the patients based on routine non-invasive
surveillance of blood parameters from one of the largest cohorts of Indian
patients at the day of admission. Patient severity and mortality prediction
models achieved 86.3% and 88.06% accuracy, respectively, with an AUC-ROC of
0.91 and 0.92. We have integrated both the models in a user-friendly web app
calculator, https://triage-COVID-19.herokuapp.com/, to showcase the potential
deployment of such efforts at scale.


  The effective reproduction number is a key figure to monitor the course of
the COVID-19 pandemic. In this study we consider a retrospective modelling
approach for estimating the effective reproduction number based on death counts
during the first year of the pandemic in Germany. The proposed Bayesian
hierarchical model incorporates splines to estimate reproduction numbers
flexibly over time while adjusting for varying effective infection fatality
rates. The approach also provides estimates of dark figures regarding
undetected infections over time. Results for Germany illustrate that estimated
reproduction numbers based on death counts are often similar to classical
estimates based on confirmed cases. However, considering death counts proves to
be more robust against shifts in testing policies: during the second wave of
infections, classical estimation of the reproduction number suggests a
flattening/ decreasing trend of infections following the "lockdown light" in
November 2020, while our results indicate that true numbers of infections
continued to rise until the "second lockdown" in December 2020. This
observation is associated with more stringent testing criteria introduced
concurrently with the "lockdown light", which is reflected in subsequently
increasing dark figures of infections estimated by our model. These findings
illustrate that the retrospective viewpoint can provide additional insights
regarding the course of the pandemic. In light of progressive vaccinations,
shifting the focus from modelling confirmed cases to reported deaths with the
possibility to incorporate effective infection fatality rates might be of
increasing relevance for the future surveillance of the pandemic.


  The paradox of the plankton highlights the apparent contradiction between
Gause's law of competitive exclusion and the observed diversity of
phytoplankton. It is well known that phytoplankton dynamics depend heavily on
two main resources: light and nutrients. Here we treat light as a continuum of
resources rather than a single resource by considering the visible light
spectrum. We propose a spatially explicit reaction-diffusion-advection model to
explore under what circumstance coexistence is possible from mathematical and
biological perspectives. Furthermore, we provide biological context as to when
coexistence is expected based on the degree of niche differentiation within the
light spectrum and overall turbidity of the water.


  We are concerned with polynomial ordinary differential systems that arise
from modelling chemical reaction networks. For such systems, which may be of
high dimension and may depend on many parameters, it is frequently of interest
to obtain a reduction of dimension in certain parameter ranges. Singular
perturbation theory, as initiated by Tikhonov and Fenichel, provides a path
toward such reductions. In the present paper we discuss parameter values that
lead to singular perturbation reductions (so-called Tikhonov-Fenichel parameter
values, or TFPVs). An algorithmic approach is known, but it is feasible for
small dimensions only. Here we characterize conditions for classes of reaction
networks for which TFPVs arise by turning off reactions (by setting rate
parameters to zero), or by removing certain species (which relates to the
classical quasi-steady state approach to model reduction). In particular, we
obtain definitive results for the class of complex balanced reaction networks
(of deficiency zero) and first order reaction networks.


  We compute the allele frequencies of the alpha (B.1.1.7), beta (B.1.351) and
delta (B.167.2) variants of SARS-CoV-2 from almost two million genome sequences
on the GISAID repository. We find that the frequencies of a majority of the
defining mutations in alpha rose towards the end of 2020 but drifted apart
during spring 2021, a similar pattern being followed by delta during summer of
2021. For beta we find a more complex scenario with frequencies of some
mutations rising and some remaining close to zero. Our results point to that
what is generally reported as single variants is in fact a collection of
variants with different genetic characteristics. For all three variants we
further find some alleles with a clearly deviating time series.


  Phylogenetic networks are a generalization of evolutionary or phylogenetic
trees that are commonly used to represent the evolution of species which cross
with one another. A special type of phylogenetic network is an {\em
$X$-cactus}, which is essentially a cactus graph in which all vertices with
degree less than three are labelled by at least one element from a set $X$ of
species. In this paper, we present a way to {\em encode} $X$-cactuses in terms
of certain collections of partitions of $X$ that naturally arise from
$X$-cactuses. Using this encoding, we also introduce a partial order on the set
of $X$-cactuses (up to isomorphism), and derive some structural properties of
the resulting partially ordered set. This includes an analysis of some
properties of its least upper and greatest lower bounds. Our results not only
extend some fundamental properties of phylogenetic trees to $X$-cactuses, but
also provides a new approach to solving topical problems in phylogenetic
network theory such as deriving consensus networks.


  Dual graphs have been applied to model RNA secondary structures with
pseudoknots, or intertwined base pairs. In previous works, a linear-time
algorithm was introduced to partition dual graphs into maximally connected
components called blocks and determine whether each block contains a pseudoknot
or not. As pseudoknots can not be contained into two different blocks, this
characterization allow us to efficiently isolate smaller RNA fragments and
classify them as pseudoknotted or pseudoknot-free regions, while keeping these
sub-structures intact. Moreover we have extended the partitioning algorithm by
classifying a pseudoknot as either recursive or non-recursive in order to
continue with our research in the development of a library of building blocks
for RNA design by fragment assembly. In this paper we present a methodology
that uses our previous results and classify pseudoknots into the classical
H,K,L, and M types, based upon a novel representation of RNA secondary
structures as dual directed graphs (i.e., digraphs). This classification would
help the systematic analysis of RNA structure and prediction as for example the
implementation of more accurate folding algorithms.


  Consider a projector matrix P, representing the first order reduced density
matrix in a basis of orthonormal atom-centric basis functions. A mathematical
question arises, and that is, how to break P into its natural component kernel
projector matrices, while preserving N-representability of P. The answer relies
upon 2- projector triple products, P'jPP'j. The triple product solutions,
applicable within the quantum crystallography of large molecules, are
determined by a new form of the Clinton equations, which - in their original
form - have long been used to ensure N-representability of density matrices
consistent with X-ray diffraction scattering factors.


  An insufficient supply of effective SARS-CoV-2 vaccine in most countries
demands an effective vaccination strategy to minimize the damage caused by the
disease. Currently, many countries vaccinate their population in descending
order of age (i.e. descending order of fatality rate) to minimize the deaths
caused by the disease; however, the effectiveness of this strategy needs to be
quantitatively assessed. We employ the susceptible-infected-recovered-dead
(SIRD) model to investigate various vaccination strategies. We constructed a
metapopulation model with heterogeneous contact and fatality rates and
investigated the effectiveness of vaccination strategies to reduce epidemic
mortality. We found that the fatality-based strategy, which is currently
employed in many countries, is more effective when the contagion rate is high
and vaccine supply is low, but the contact-based method outperforms the
fatality-based strategy when there is a sufficiently high supply of the
vaccine. We identified a discontinuous transition of the optimal vaccination
strategy and path-dependency analogous to hysteresis. This transition and
path-dependency imply that combining the fatality-based and contact-based
strategies is ineffective in reducing the number of deaths. Furthermore, we
demonstrate that such phenomena occur in real-world epidemic diseases, such as
tuberculosis and COVID-19. We also show that the conclusions of this research
are valid even when the complex epidemic stages, efficacy of the vaccine, and
reinfection are considered.


  The SI model is the most basic of all compartmental models used to describe
the spreading of information through a population. Despite its apparent
simplicity, the analytic solution of this model on networks is still lacking.
We address this problem here, using a novel formulation inspired by the
mathematical treatment of many-body quantum systems. This allows us to organize
the time-dependent expectation values for the state of individual nodes in
terms of contributions from subgraphs of the network. We compute these
contributions systematically and find a set of symmetry relations among
subgraphs of differing topologies. We use our novel approach to compute the
spreading of information on three different sample networks. The exact
solution, which matches with Monte-Carlo simulations, visibly departs from the
mean-field results.


  The analysis of differential gene expression from RNA-Seq data has become a
standard for several research areas mainly involving bioinformatics. The steps
for the computational analysis of these data include many data types and file
formats, and a wide variety of computational tools that can be applied alone or
together as pipelines. This paper presents a review of differential expression
analysis pipeline, addressing its steps and the respective objectives, the
principal methods available in each step and their properties, bringing an
overview in an organized way in this context. In particular, this review aims
to address mainly the aspects involved in the differentially expressed gene
(DEG) analysis from RNA sequencing data (RNA-Seq), considering the
computational methods and its properties. In addition, a timeline of the
evolution of computational methods for DEG is presented and discussed, as well
as the relationships existing between the main computational tools are
presented by an interaction network. A discussion on the challenges and gaps in
DEG analysis is also highlighted in this review.


  Understanding how to effectively control an epidemic spreading on a network
is a problem of paramount importance for the scientific community. The ongoing
COVID-19 pandemic has highlighted the need for policies that mitigate the
spread, without relying on pharmaceutical interventions, that is, without the
medical assurance of the recovery process. These policies typically entail
lockdowns and mobility restrictions, having thus nonnegligible socio-economic
consequences for the population. In this paper, we focus on the problem of
finding the optimum policies that "flatten the epidemic curve" while limiting
the negative consequences for the society, and formulate it as a nonlinear
control problem over a finite prediction horizon. We utilize the model
predictive control theory to design a strategy to effectively control the
disease, balancing safety and normalcy. An explicit formalization of the
control scheme is provided for the susceptible--infected--susceptible epidemic
model over a network. Its performance and flexibility are demonstrated by means
of numerical simulations.


  Biocatalysis is a promising approach to sustainably synthesize
pharmaceuticals, complex natural products, and commodity chemicals at scale.
However, the adoption of biocatalysis is limited by our ability to select
enzymes that will catalyze their natural chemical transformation on non-natural
substrates. While machine learning and in silico directed evolution are
well-posed for this predictive modeling challenge, efforts to date have
primarily aimed to increase activity against a single known substrate, rather
than to identify enzymes capable of acting on new substrates of interest. To
address this need, we curate 6 different high-quality enzyme family screens
from the literature that each measure multiple enzymes against multiple
substrates. We compare machine learning-based compound-protein interaction
(CPI) modeling approaches from the literature used for predicting drug-target
interactions. Surprisingly, comparing these interaction-based models against
collections of independent (single task) enzyme-only or substrate-only models
reveals that current CPI approaches are incapable of learning interactions
between compounds and proteins in the current family level data regime. We
further validate this observation by demonstrating that our no-interaction
baseline can outperform CPI-based models from the literature used to guide the
discovery of kinase inhibitors. Given the high performance of non-interaction
based models, we introduce a new structure-based strategy for pooling residue
representations across a protein sequence. Altogether, this work motivates a
principled path forward in order to build and evaluate meaningful predictive
models for biocatalysis and other drug discovery applications.


  Boltzmann machines are energy-based models that have been shown to provide an
accurate statistical description of domains of evolutionary-related protein and
RNA families. They are parametrized in terms of local biases accounting for
residue conservation, and pairwise terms to model epistatic coevolution between
residues. From the model parameters, it is possible to extract an accurate
prediction of the three-dimensional contact map of the target domain. More
recently, the accuracy of these models has been also assessed in terms of their
ability in predicting mutational effects and generating in silico functional
sequences. Our adaptive implementation of Boltzmann machine learning, adabmDCA,
can be generally applied to both protein and RNA families and accomplishes
several learning set-ups, depending on the complexity of the input data and on
the user requirements. The code is fully available at
https://github.com/anna-pa-m/adabmDCA. As an example, we have performed the
learning of three Boltzmann machines modeling the Kunitz and Beta-lactamase2
protein domains and TPP-riboswitch RNA domain. The models learned by adabmDCA
are comparable to those obtained by state-of-the-art techniques for this task,
in terms of the quality of the inferred contact map as well as of the
synthetically generated sequences. In addition, the code implements both
equilibrium and out-of-equilibrium learning, which allows for an accurate and
lossless training when the equilibrium one is prohibitive in terms of
computational time, and allows for pruning irrelevant parameters using an
information-based criterion.


  RNA 3D architectures are stabilized by sophisticated networks of
(non-canonical) base pair interactions, which can be conveniently encoded as
multi-relational graphs and efficiently exploited by graph theoretical
approaches and recent progresses in machine learning techniques. RNAglib is a
library that eases the use of this representation, by providing clean data,
methods to load it in machine learning pipelines and graph-based deep learning
models suited for this representation. RNAglib also offers other utilities to
model RNA with 2.5D graphs, such as drawing tools, comparison functions or
baseline performances on RNA applications. The method and data is distributed
as a fully documented pip package.
  Availability: https://rnaglib.cs.mcgill.ca


  Deep neural networks such as AlphaFold and RoseTTAFold predict remarkably
accurate structures of proteins compared to other algorithmic approaches. It is
known that biologically small perturbations in the protein sequence do not lead
to drastic changes in the protein structure. In this paper, we demonstrate that
RoseTTAFold does not exhibit such a robustness despite its high accuracy, and
biologically small perturbations for some input sequences result in radically
different predicted protein structures. This raises the challenge of detecting
when these predicted protein structures cannot be trusted. We define the
robustness measure for the predicted structure of a protein sequence to be the
inverse of the root-mean-square distance (RMSD) in the predicted structure and
the structure of its adversarially perturbed sequence. We use adversarial
attack methods to create adversarial protein sequences, and show that the RMSD
in the predicted protein structure ranges from 0.119\r{A} to 34.162\r{A} when
the adversarial perturbations are bounded by 20 units in the BLOSUM62 distance.
This demonstrates very high variance in the robustness measure of the predicted
structures. We show that the magnitude of the correlation (0.917) between our
robustness measure and the RMSD between the predicted structure and the ground
truth is high, that is, the predictions with low robustness measure cannot be
trusted. This is the first paper demonstrating the susceptibility of
RoseTTAFold to adversarial attacks.


  The recent global surge in COVID-19 infections has been fueled by new
SARS-CoV-2 variants, namely Alpha, Beta, Gamma, Delta, etc. The molecular
mechanism underlying such surge is elusive due to 4,653 non-degenerate
mutations on the spike protein, which is the target of most COVID-19 vaccines.
The understanding of the molecular mechanism of transmission and evolution is a
prerequisite to foresee the trend of emerging vaccine-breakthrough variants and
the design of mutation-proof vaccines and monoclonal antibodies. We integrate
the genotyping of 1,489,884 SARS-CoV-2 genomes isolates, 130 human antibodies,
tens of thousands of mutational data points, topological data analysis, and
deep learning to reveal SARS-CoV-2 evolution mechanism and forecast emerging
vaccine-escape variants. We show that infectivity-strengthening and
antibody-disruptive co-mutations on the S protein RBD can quantitatively
explain the infectivity and virulence of all prevailing variants. We
demonstrate that Lambda is as infectious as Delta but is more
vaccine-resistant. We analyze emerging vaccine-breakthrough co-mutations in 20
countries, including the United Kingdom, the United States, Denmark, Brazil,
and Germany, etc. We envision that natural selection through infectivity will
continue to be the main mechanism for viral evolution among unvaccinated
populations, while antibody disruptive co-mutations will fuel the future growth
of vaccine-breakthrough variants among fully vaccinated populations. Finally,
we have identified the co-mutations that have the great likelihood of becoming
dominant: [A411S, L452R, T478K], [L452R, T478K, N501Y], [V401L, L452R, T478K],
[K417N, L452R, T478K], [L452R, T478K, E484K, N501Y], and [P384L, K417N, E484K,
N501Y]. We predict they, particularly the last four, will break through
existing vaccines. We foresee an urgent need to develop new vaccines that
target these co-mutations.


  How to strategically allocate the available vaccines is a crucial issue for
pandemic control. In this work, we propose a mathematical framework for optimal
stabilizing vaccine allocation, where our goal is to send the infections to
zero as soon as possible with a fixed number of vaccine doses. This framework
allows us to efficiently compute the optimal vaccine allocation policy for
general epidemic spread models including SIS/SIR/SEIR and a new model of
COVID-19 transmissions. By fitting the real data in New York State to our
framework, we found that the optimal stabilizing vaccine allocation policy
suggests offering vaccines priority to locations where there are more
susceptible people and where the residents spend longer time outside the home.
Besides, we found that offering vaccines priority to young adults (20-29) and
middle-age adults (20-44) can minimize the cumulative infected cases and the
death cases. Moreover, we compared our method with five age-stratified
strategies in \cite{bubar2021model} based on their epidemics model. We also
found it's better to offer vaccine priorities to young people to curb the
disease and minimize the deaths when the basic reproduction number $R_0$ is
moderately above one, which describes the most world during COVID-19. Such
phenomenon has been ignored in \cite{bubar2021model}.


  In cellular reprogramming, almost all epigenetic memories of differentiated
cells are erased by the overexpression of few genes, regaining pluripotency,
potentiality for differentiation. Considering the interplay between oscillatory
gene expression and slower epigenetic modifications, such reprogramming is
perceived as an unintuitive, global attraction to the unstable manifold of a
saddle, which represents pluripotency. The universality of this scheme is
confirmed by the repressilator model, and by gene regulatory networks randomly
generated and those extracted from embryonic stem cells.


  Routine single-sample haplotype-resolved assembly remains an unresolved
problem. Here we describe a new algorithm that combines PacBio HiFi reads and
Hi-C chromatin interaction data to produce a haplotype-resolved assembly
without the sequencing of parents. Applied to human and other vertebrate
samples, our algorithm consistently outperforms existing single-sample assembly
pipelines and generates assemblies of comparable quality to the best
pedigree-based assemblies.


  Scaffold based drug discovery (SBDD) is a technique for drug discovery which
pins chemical scaffolds as the framework of design. Scaffolds, or molecular
frameworks, organize the design of compounds into local neighborhoods. We
formalize scaffold based drug discovery into a network design. Utilizing
docking data from SARS-CoV-2 virtual screening studies and JAK2 kinase assay
data, we showcase how a scaffold based conception of chemical space is
intuitive for design. Lastly, we highlight the utility of scaffold based
networks for chemical space as a potential solution to the intractable
enumeration problem of chemical space by working inductively on local
neighborhoods.


  With the rapid global spread of COVID-19, more and more data related to this
virus is becoming available, including genomic sequence data. The total number
of genomic sequences that are publicly available on platforms such as GISAID is
currently several million, and is increasing with every day. The availability
of such \emph{Big Data} creates a new opportunity for researchers to study this
virus in detail. This is particularly important with all of the dynamics of the
COVID-19 variants which emerge and circulate. This rich data source will give
us insights on the best ways to perform genomic surveillance for this and
future pandemic threats, with the ultimate goal of mitigating or eliminating
such threats. Analyzing and processing the several million genomic sequences is
a challenging task. Although traditional methods for sequence classification
are proven to be effective, they are not designed to deal with these specific
types of genomic sequences. Moreover, most of the existing methods also face
the issue of scalability. Previous studies which were tailored to coronavirus
genomic data proposed to use spike sequences (corresponding to a subsequence of
the genome), rather than using the complete genomic sequence, to perform
different machine learning (ML) tasks such as classification and clustering.
However, those methods suffer from scalability issues. In this paper, we
propose an approach called Spike2Vec, an efficient and scalable feature vector
representation for each spike sequence that can be used for downstream ML
tasks. Through experiments, we show that Spike2Vec is not only scalable on
several million spike sequences, but also outperforms the baseline models in
terms of prediction accuracy, F1 score, etc.


  Complex pattern formation is an essential characteristic of plants and their
ageing, growth, and evolution. Perception of these patterns is an intrinsic
nature of plant-dependent animals for coexistence. Areca nut consisting of
complex patterns is considered to be addictive for humans and has increased
adverse health effects. However, no critical study is performed on the complex
pattern of the areca nut. A large number of areca nuts has been studied since
2017 to develop a low-cost tool for the LMICs to categorise areca nuts. We
present the first finding to identify similarities among complex networks of
differently aged areca nuts by investigating the internal patterns of randomly
chosen nuts from the same age group. We developed a smartphone camera-based
high-resolution measurement with comprehensive biophysical mathematics and a
quantum mechanical concept called density of states (DOS). We found that the
DOS can provide a unique coefficient to represent age and ageing together using
a single number. The average of these single numbers for less aged nuts and
highly aged nuts are 4.9 and 3.8, respectively. If fruit looks aged from its
external morphology as well as internal morphology, our method identifies the
intrinsic similarities among the ageing networks without implementing any
computationally expensive search algorithm. We show clear evidence of the
diversity of ageing from relative and absolute colour vision perspectives. We
have also conducted further analyses of local DOS, Fourier decomposition,
correlation study, spectral decomposition, and structural similarity index.


  In epidemic or pandemic situations, resources for testing the infection
status of individuals may be scarce. Although group testing can help to
significantly increase testing capabilities, the (repeated) testing of entire
populations can exceed the resources of any country. We thus propose an
extension of the theory of group testing that takes into account the fact that
definitely specifying the infection status of each individual is impossible.
Our theory builds on assigning to each individual an infection status
(healthy/infected), as well as an associated cost function for erroneous
assignments. This cost function is versatile, e.g., it could take into account
that false negative assignments are worse than false positive assignments and
that false assignments in critical areas, such as health care workers, are more
severe than in the general population. Based on this model, we study the
optimal use of a limited number of tests to minimize the expected cost. More
specifically, we utilize information-theoretic methods to give a lower bound on
the expected cost and describe simple strategies that can significantly reduce
the expected cost over currently known strategies. A detailed example is
provided to illustrate our theory.


  As the global need for large-scale data storage is rising exponentially,
existing storage technologies are approaching their theoretical and functional
limits in terms of density and energy consumption, making DNA based storage a
potential solution for the future of data storage. Several studies introduced
DNA based storage systems with high information density (petabytes/gram).
However, DNA synthesis and sequencing technologies yield erroneous outputs.
Algorithmic approaches for correcting these errors depend on reading multiple
copies of each sequence and result in excessive reading costs. The
unprecedented success of Transformers as a deep learning architecture for
language modeling has led to its repurposing for solving a variety of tasks
across various domains. In this work, we propose a novel approach for
single-read reconstruction using an encoder-decoder Transformer architecture
for DNA based data storage. We address the error correction process as a
self-supervised sequence-to-sequence task and use synthetic noise injection to
train the model using only the decoded reads. Our approach exploits the
inherent redundancy of each decoded file to learn its underlying structure. To
demonstrate our proposed approach, we encode text, image and code-script files
to DNA, produce errors with high-fidelity error simulator, and reconstruct the
original files from the noisy reads. Our model achieves lower error rates when
reconstructing the original data from a single read of each DNA strand compared
to state-of-the-art algorithms using 2-3 copies. This is the first
demonstration of using deep learning models for single-read reconstruction in
DNA based storage which allows for the reduction of the overall cost of the
process. We show that this approach is applicable for various domains and can
be generalized to new domains as well.


  The waning of immunity after recovery or vaccination is a major factor
accounting for the severity and prolonged duration of an array of epidemics,
ranging from COVID-19 to diphtheria and pertussis. To study the effectiveness
of different immunity level-based vaccination schemes in mitigating the impact
of waning immunity, we construct epidemiological models that mimic the latter's
effect. The total susceptible population is divided into an arbitrarily large
number of discrete compartments with varying levels of disease immunity. We
then vaccinate various compartments within this framework, comparing the value
of $R_0$ and the equilibria locations for our systems to determine an optimal
immunization scheme under natural constraints. Relying on perturbative
analysis, we establish a number of results concerning the location, existence,
and uniqueness of the system's endemic equilibria, as well as results on
disease-free equilibria. In addition, we numerically simulate the dynamics
associated with our model in the case of pertussis in Canada, fitting our model
to available time-series data. Our analytical results are applicable to a wide
range of systems composed of arbitrarily many ODEs.


  "Personal responsibility", one of the basic principles of modern law,
requires one to be responsible for what he did. However, personal
responsibility is far from the only norm ruling human interactions, especially
in social and economic activities. In many collective communities such as among
enterprise colleagues and family members, one's personal interests are often
bound to others' -- once one member breaks the rule, a group of people have to
bear the punishment or sanction. Such a mechanism is termed "joint liability".
Although many real-world cases have demonstrated that joint liability helps to
maintain collective collaboration, a deep and systematic theoretical analysis
on how and when joint liability promotes cooperation is lacking. Here we use
evolutionary game theory to model an interacting system with joint liability,
where one's losing credit could deteriorate the reputation of the whole group.
We provide the analytical condition to predict when cooperation evolves in the
presence of joint liability, which is verified by simulations. We also
analytically prove that joint liability can greatly promote cooperation. Our
work stresses that joint liability is of great significance in promoting the
current economic propensity.


  We consider the dynamics of a three-species system incorporating the Allee
Effect, focussing on its influence on the emergence of extreme events in the
system. First we find that under Allee effect the regular periodic dynamics
changes to chaotic. Further, we find that the system exhibits unbounded growth
in the vegetation population after a critical value of the Allee parameter. The
most significant finding is the observation of a critical Allee parameter
beyond which the probability of obtaining extreme events becomes non-zero for
all three population densities. Though the emergence of extreme events in the
predator population is not affected much by the Allee effect, the prey
population shows a sharp increase in the probability of obtaining extreme
events after a threshold value of the Allee parameter, and the vegetation
population also yields extreme events for sufficiently strong Allee effect.
Lastly we consider the influence of additive noise on extreme events. First, we
find that noise tames the unbounded vegetation growth induced by Allee effect.
More interestingly, we demonstrate that stochasticity drastically diminishes
the probability of extreme events in all three populations. In fact for
sufficiently high noise, we do not observe any more extreme events in the
system. This suggests that noise can mitigate extreme events, and has
potentially important bearing on the observability of extreme events in
naturally occurring systems.


  We present a detailed set-based analysis of the well-known SIR and SEIR
epidemic models subjected to hard caps on the proportion of infective
individuals, and bounds on the allowable intervention strategies, such as
social distancing, quarantining and vaccination. We describe the admissible and
maximal robust positively invariant (MRPI) sets of these two models via the
theory of barriers. We show how the sets may be used in the management of
epidemics, for both perfect and imperfect/uncertain models, detailing how
intervention strategies may be specified such that the hard infection cap is
never breached, regardless of the basic reproduction number. The results are
clarified with detailed examples.


  Motivation: Omics data, such as transcriptomics or phosphoproteomics, are
broadly used to get a snap-shot of the molecular status of cells. In
particular, changes in omics can be used to estimate the activity of pathways,
transcription factors and kinases based on known regulated targets, that we
call footprints. Then the molecular paths driving these activities can be
estimated using causal reasoning on large signaling networks. Results: We have
developed FUNKI, a FUNctional toolKIt for footprint analysis. It provides a
user-friendly interface for an easy and fast analysis of several omics data,
either from bulk or single-cell experiments. FUNKI also features different
options to visualise the results and run post-analyses, and is mirrored as a
scripted version in R. Availability: FUNKI is a free and open-source
application built on R and Shiny, available in GitHub at
https://github.com/saezlab/ShinyFUNKI under GNU v3.0 license and accessible
also in https://saezlab.shinyapps.io/funki/ Contact: pub.saez@uni-heidelberg.de
Supplementary information: We provide data examples within the app, as well as
extensive information about the different variables to select, the results, and
the different plots in the help page.


  We revisit the spatial ${\lambda}$-Fleming-Viot process introduced in [1].
Particularly, we are interested in the time $T_0$ to the most recent common
ancestor for two lineages. We distinguish between the case where the process
acts on the entire two-dimensional plane, and on a finite rectangle. Utilizing
a differential equation linking $T_0$ with the physical distance between the
lineages, we arrive at simple and reasonably accurate approximation schemes for
both cases. Furthermore, our analysis enables us to address the question of
whether the genealogical process of the model "comes down from infinity", which
has been partly answered before in [2].


  Methods for time series prediction and classification of gene regulatory
networks (GRNs) from gene expression data have been treated separately so far.
The recent emergence of attention-based recurrent neural networks (RNN) models
boosted the interpretability of RNN parameters, making them appealing for the
understanding of gene interactions. In this work, we generated synthetic time
series gene expression data from a range of archetypal GRNs and we relied on a
dual attention RNN to predict the gene temporal dynamics. We show that the
prediction is extremely accurate for GRNs with different architectures. Next,
we focused on the attention mechanism of the RNN and, using tools from graph
theory, we found that its graph properties allow to hierarchically distinguish
different architectures of the GRN. We show that the GRNs respond differently
to the addition of noise in the prediction by the RNN and we relate the noise
response to the analysis of the attention mechanism. In conclusion, this work
provides a a way to understand and exploit the attention mechanism of RNN and
it paves the way to RNN-based methods for time series prediction and inference
of GRNs from gene expression data.


  Despite many studies on the transmission mechanism of the Severe acute
respiratory syndrome coronavirus 2 (SARS-CoV-2), it remains still challenging
to efficiently reduce mortality. In this work, we apply a two-population
Susceptible-Infected-Removed (SIR) model to investigate the COVID-19 spreading
when contacts between elderly and non-elderly individuals are reduced due to
the high mortality risk of elderly people. We discover that the reduction of
connections between two populations can delay the death curve but cannot well
reduce the final mortality. We propose a merged SIR model, which advises
elderly individuals to interact less with their non-elderly connections at the
initial stage but interact more with their non-elderly relationships later, to
reduce the final mortality. Finally, immunizing elderly hub individuals can
also significantly decrease mortality.


  The deer tick, $\textit{Ixodes scapularis}$, is a vector for numerous human
diseases, including Lyme disease, anaplasmosis, and babesiosis. Concern is
rising in the US and abroad as the population and range of this species grow
and new diseases emerge. Herein I consider the potential for control of
$\textit{I. scapularis}$ using the Sterile Insect Technique (SIT), which acts
by reducing net fertility through release of sterile males. I construct a
population model with density-dependent and -independent growth, migration, and
an Allee effect (decline of the population when it is small), and use this
model to simulate sterile tick release in both single- and multi-patch
frameworks. I test two key concerns with implementing $\textit{I. scapularis}$
SIT: that the ticks' lengthy life course could make control take too long and
that low migration might mean sterile males need thorough manual dispersal to
all parts of the control area. Results suggest that typical $\textit{I.
scapularis}$ SIT programs will take about eight years, a prediction near the
normal range for the technique, but that thorough distribution of sterile ticks
over the control area is indeed critical, increasing expense substantially by
necessitating aerial release. With particularly high rearing costs also
expected for $\textit{I. scapularis}$, the latter finding suggests that
cost-effectiveness improvements to aerial release may be a prerequisite to
$\textit{I. scapularis}$ SIT.


  Since Press and Dyson's ingenious discovery of ZD (zero-determinant) strategy
in the repeated Prisoner's Dilemma game, several studies have confirmed the
existence of ZD strategy in repeated multiplayer social dilemmas. However, few
researches study the evolutionary performance of multiplayer ZD strategies,
especially from a theoretical perspective. Here, we use a newly proposed
state-clustering method to theoretically analyze the evolutionary dynamics of
two representative ZD strategies: generous ZD strategies and extortionate ZD
strategies. Apart from the competitions between the two strategies and some
classical strategies, we consider two new settings for multiplayer ZD
strategies: competitions in the whole ZD strategy space and competitions in the
space of all memory-1 strategies. Besides, we investigate the influence of
level of generosity and extortion on the evolutionary dynamics of generous and
extortionate ZD, which was commonly ignored in previous studies. Theoretical
results show players with limited generosity are at an advantageous place and
extortioners extorting more severely hold their ground more readily. Our
results may provide new insights into better understanding the evolutionary
dynamics of ZD strategies in repeated multiplayer games.


  We study the impact on the epidemiological dynamics of a class of restrictive
measures that are aimed at reducing the number of contacts of individuals who
have a higher risk of being infected with a transmittable disease. Such
measures are currently either implemented or at least discussed in numerous
countries worldwide to ward off a potential new wave of COVID-19 across Europe.
They come in the form of Health Passes (HP), which grant full access to public
life only to individuals with a certificate that proves that they have either
been fully vaccinated, have recovered from a previous infection or have
recently tested negative to SARS-Cov-19 . We develop both a compartmental model
as well as an epidemic Renormalisation Group approach, which is capable of
describing the dynamics over a longer period of time, notably an entire
epidemiological wave. Introducing different versions of HPs in this model, we
are capable of providing quantitative estimates on the effectiveness of the
underlying measures as a function of the fraction of the population that is
vaccinated and the vaccination rate. We apply our models to the latest COVID-19
wave in several European countries, notably Germany and Austria, which validate
our theoretical findings.


  There are many mathematical models of biochemical cell signaling pathways
that contain a large number of elements (species and reactions). This is
sometimes a big issue for identifying critical model elements and describing
the model dynamics. Thus, techniques of model reduction can be used as a
mathematical tool in order to minimize the number of variables and parameters.
In this thesis, we review some well-known methods of model reduction for cell
signaling pathways. We have also developed some approaches that provide us a
great step forward in model reduction. The techniques are quasi steady state
approximation (QSSA), quasi equilibrium approximation (QEA), lumping of species
and entropy production analysis. They are applied on protein translation
pathways with microRNA mechanisms, chemical reaction networks, extracellular
signal regulated kinase (ERK) pathways, NFkB signal transduction pathways,
elongation factors EFTu and EFTs signaling pathways and Dihydrofolate reductase
(DHFR) pathways. The main aim of this thesis is to reduce the complex cell
signaling pathway models. This provides one a better understanding of the
dynamics of such models and gives an accurate approximate solution. Results
show that there is a good agreement between the original models and the
simplified models.


  Chemotaxis-driven invasions have been proposed across a broad spectrum of
biological processes, from cancer to ecology. The influential system of
equations introduced by Keller and Segel has proven a popular choice in the
modelling of such phenomena, but in its original form restricts to a
homogeneous population. To account for the possibility of phenotypic
heterogeneity, we extend to the case of a population continuously structured
across space, time and phenotype, where the latter determines variation in
chemotactic responsiveness, proliferation rate, and the level of chemical
environment modulation. The extended model considered here comprises a
non-local partial differential equation for the local phenotype distribution of
cells which is coupled, through an integral term, with a differential equation
for the concentration of an attractant, which is sensed and degraded by the
cells. In the framework of this model, we concentrate on a
chemotaxis/proliferation trade-off scenario, where the cell phenotypes span a
spectrum of states from highly-chemotactic but minimally-proliferative to
minimally-chemotactic but highly-proliferative. Using a combination of
numerical simulation and formal asymptotic analysis, we explore the properties
of travelling-wave solutions. The results of our study demonstrate how
incorporating phenotypic heterogeneity may lead to a highly-structured wave
profile, where cells in different phenotypic states dominate different spatial
positions across the invading wave, and clarify how the phenotypic structuring
of the wave can be shaped by trade-offs between chemotaxis and proliferation.


  A decomposition of a chemical reaction network (CRN) is produced by
partitioning its set of reactions. The partition induces networks, called
subnetworks, that are "smaller" than the given CRN which, at this point, can be
called parent network. A complex is called a common complex if it occurs in at
least two subnetworks in a decomposition. A decomposition is said to be
incidence independent if the image of the incidence map of the parent network
is the direct sum of the images of the subnetworks' incidence maps. It has been
recently discovered that the complex balanced equilibria of the parent network
and its subnetworks are fundamentally connected in an incidence independent
decomposition. In this paper, we utilized the set of common complexes and a
developed criterion to investigate decomposition's incidence independence
properties. A framework was also developed to analyze decomposition classes
with similar structure and incidence independence properties. We identified
decomposition classes that can be characterized by their sets of common
complexes and studied their incidence independence. Some of these decomposition
classes occur in some biological and chemical models. Finally, a sufficient
condition was obtained for the complex balancing of some power law kinetic
(PLK) systems with incidence independent and complex balanced decompositions.
This condition led to a generalization of the Defficiency Zero Theorem for some
PLK systems.


  Coarse-graining is a powerful tool for extending the reach of dynamic models
of proteins and other biological macromolecules. Topological coarse-graining,
in which biomolecules or sets thereof are represented via graph structures, is
a particularly useful way of obtaining highly compressed representations of
molecular structure, and simulations operating via such representations can
achieve substantial computational savings. A drawback of coarse-graining,
however, is the loss of atomistic detail - an effect that is especially acute
for topological representations such as protein structure networks (PSNs).
Here, we introduce an approach based on a combination of machine learning and
physically-guided refinement for inferring atomic coordinates from PSNs. This
"neural upscaling" procedure exploits the constraints implied by PSNs on
possible configurations, as well as differences in the likelihood of observing
different configurations with the same PSN. Using a 1 $\mu$s atomistic
molecular dynamics trajectory of A$\beta_{1-40}$, we show that neural upscaling
is able to effectively recapitulate detailed structural information for
intrinsically disordered proteins, being particularly successful in recovering
features such as transient secondary structure. These results suggest that
scalable network-based models for protein structure and dynamics may be used in
settings where atomistic detail is desired, with upscaling employed to impute
atomic coordinates from PSNs.


  Previous studies of the flexibilities of ancestral proteins suggests that
proteins evolve their function by altering their native state ensemble. Here we
propose a more direct method of visualizing this by measuring the changes in
the vibrational density of states (VDOS) of proteins as they evolve. Through
analysis of VDOS profiles of ancestral and extant proteins we observe that
$\beta$-lactamase and thioredoxins evolve by altering their density of states
in the terahertz region. Particularly, the shift in VDOS profiles between
ancestral and extant proteins suggests that nature utilize dynamic allostery
for functional evolution. Moreover, we also show that VDOS profile of
individual position can be used to describe the flexibility changes,
particularly those without any amino acid substitution.


  The multidomain protein gelsolin (GSN) is composed of six homologous modules,
sequentially named G1 to G6. Single point substitutions in this protein are
responsible for AGel amyloidosis, a hereditary disease characterized by
progressive corneal lattice dystrophy, cutis laxa, and polyneuropathy. Several
different amyloidogenic variants of GSN have been identified over the years,
but only the most common D187N/Y mutants, in G2, have been thoroughly
characterized, and the underlying functional mechanistic link between mutation,
altered protein structure, susceptibility to aberrant furin cleavage and
aggregative potential resolved. Little is known about the recently identified
mutations A551P, E553K and M517R hosted at the interface between G4 and G5,
whose aggregation process likely follows an alternative pathway. We demonstrate
that these three substitutions impair temperature and pressure stability of GSN
but do not increase its susceptibility to furin cleavage, the first event of
the canonical aggregation pathway. The variants are also characterized by a
higher tendency to aggregate in the unproteolysed forms and show a higher
proteotoxicity in a C. elegans-based assay. Structural studies point to a
destabilization of the interface between G4 and G5 due to three different
structural determinants: beta-strand breaking, steric hindrance and/or charge
repulsion, all implying the impairment of interdomain contacts. All available
evidence suggests that the rearrangement of the protein global architecture
triggers a furin-independent aggregation of the protein, supporting the
existence of a non-canonical pathway of gelsolin amyloidosis pathogenesis.


  The COVID-19 pandemic, caused by the coronavirus SARS-CoV-2, has led to a
wide range of non-pharmaceutical interventions being implemented around the
world to curb transmission. However, the economic and social costs of some of
these measures, especially lockdowns, has been high. An alternative and widely
discussed public health strategy for the COVID-19 pandemic would have been to
'shield' those most vulnerable to COVID-19 (minimising their contacts with
others), while allowing infection to spread among lower risk individuals with
the aim of reaching herd immunity. Here we retrospectively explore the
effectiveness of this strategy using a stochastic SEIR framework, showing that
even under the unrealistic assumption of perfect shielding, hospitals would
have been rapidly overwhelmed with many avoidable deaths among lower risk
individuals. Crucially, even a small (20%) reduction in the effectiveness of
shielding would have likely led to a large increase (>150%) in the number of
deaths compared to perfect shielding. Our findings demonstrate that shielding
the vulnerable while allowing infections to spread among the wider population
would not have been a viable public health strategy for COVID-19 and is
unlikely to be effective for future pandemics.


  Arid zones contain a diverse set of microbes capable of survival under dry
conditions, some of which can form relationships with plants under drought
stress conditions to improve plant health. We studied squash (Cucurbita pepo
L.) root microbiome under historically arid and humid sites, both in situ and
performing a common garden experiment. Plants were grown in soils from sites
with different drought levels, using in situ collected soils as the microbial
source. We described and analyzed bacterial diversity by 16S rRNA gene
sequencing (N=48) from the soil, rhizosphere, and endosphere. Proteobacteria
were the most abundant phylum present in humid and arid samples, while
Actinobacteriota abundance was higher in arid ones. The Beta-diversity analyses
showed split microbiomes between arid and humid microbiomes, and aridity and
soil pH levels could explain it. These differences between humid and arid
microbiomes were maintained in the common garden experiment, showing that it is
possible to transplant in situ diversity to the greenhouse. We detected a total
of 1009 bacterial genera; 199 exclusively associated with roots under arid
conditions. With shotgun metagenomic sequencing of rhizospheres (N=6), we
identified 2969 protein families in the squash core metagenome and found an
increased number of exclusively protein families from arid (924) than humid
samples (158). We found arid conditions enriched genes involved in protein
degradation and folding, oxidative stress, compatible solute synthesis, and ion
pumps associated with osmotic regulation. Plant phenotyping allowed us to
correlate bacterial communities with plant growth. Our study revealed that it
is possible to evaluate microbiome diversity ex-situ and identify critical
species and genes involved in plant-microbe interactions in historically arid
locations.


  The cataclysmic contagion based calamity -- Covid-19 has shown us a clear
need for a comprehensive community based strategy that overcomes the sheer
complexity of controlling it and the caveats of current methods. In this
regard, as seen in earlier epidemics, testing has always been an integral part
of containment policy. But one has to consider the optimality of a testing
scheme based on the resultant disease spread in the community and not based on
purely increasing testing efficiency. Therefore, taking a decision is no easy
feat and must consider the community utility constrained by its priorities,
budget, risks, collateral and abilities which can be encoded into the
optimization of the strategy. We thus propose a simple pooling strategy that is
easy to customize and practical to implement, unlike other complex and
computationally intensive methods.


  Bayesian inference is a popular and widely-used approach to infer phylogenies
(evolutionary trees). However, despite decades of widespread application, it
remains difficult to judge how well a given Bayesian Markov chain Monte Carlo
(MCMC) run explores the space of phylogenetic trees. In this paper, we
investigate the Monte Carlo error of phylogenies, focusing on high-dimensional
summaries of the posterior distribution, including variability in estimated
edge/branch (known in phylogenetics as "split") probabilities and tree
probabilities, and variability in the estimated summary tree. Specifically, we
ask if there is any measure of effective sample size (ESS) applicable to
phylogenetic trees which is capable of capturing the Monte Carlo error of these
three summary measures. We find that there are some ESS measures capable of
capturing the error inherent in using MCMC samples to approximate the posterior
distributions on phylogenies. We term these tree ESS measures, and identify a
set of three which are useful in practice for assessing the Monte Carlo error.
Lastly, we present visualization tools that can improve comparisons between
multiple independent MCMC runs by accounting for the Monte Carlo error present
in each chain. Our results indicate that common post-MCMC workflows are
insufficient to capture the inherent Monte Carlo error of the tree, and
highlight the need for both within-chain mixing and between-chain convergence
assessments.


  We consider a population constituted by two types of individuals; each of
them can produce offspring in two different islands (as a particular case the
islands can be interpreted as active or dormant individuals). We model the
evolution of the population of each type using a two-type Feller diffusion with
immigration, and we study the frequency of one of the types, in each island,
when the total population size in each island is forced to be constant at a
dense set of times. This leads to the solution of a SDE which we call the
asymmetric two-island frequency process. We derive properties of this process
and obtain a large population limit when the total size of each island tends to
infinity. Additionally, we compute the fluctuations of the process around its
deterministic limit. We establish conditions under which the asymmetric
two-island frequency process has a moment dual. The dual is a continuous-time
two-dimensional Markov chain that can be interpreted in terms of mutation,
branching, pairwise branching, coalescence, and a novel mixed
selection-migration term. Also, we conduct a stability analysis of the limiting
deterministic dynamical system and present some numerical results to study
fixation and a new form of balancing selection. When restricting to the
seedbank model, we observe that some combinations of the parameters lead to
balancing selection. Besides finding yet another way in which genetic
reservoirs increase the genetic variability, we find that if a population that
sustains a seedbank competes with one that does not, the seed producers will
have a selective advantage if they reproduce faster, but will not have a
selective disadvantage if they reproduce slower: their worst case scenario is
balancing selection.


  To contain the spread of the virus and stop the overcrowding of hospitalized
patients, the coronavirus pandemic crippled healthcare facilities, mandating
lockdowns and promoting remote work. As a result, telehealth has become
increasingly popular for offering low-risk care to patients. However, the
difficulty of preventing the next potential waves of infection has increased by
constant virus mutation into new forms and a general lack of test kits,
particularly in developing nations. In this research, a unique cloud-based
application for the early identification of individuals who may have COVID-19
infection is proposed. The application provides five modes of diagnosis from
possible symptoms (f1), cough sound (f2), specific blood biomarkers (f3), Raman
spectral data of blood specimens (f4), and ECG signal paper-based image (f5).
When a user selects an option and enters the information, the data is sent to
the cloud server. The deployed machine learning (ML) and deep learning (DL)
models classify the data in real time and inform the user of the likelihood of
COVID-19 infection. Our deployed models can classify with an accuracy of 100%,
99.80%, 99.55%, 95.65%, and 77.59% from f3, f4, f5, f2, and f1 respectively.
Moreover, the sensitivity for f2, f3, and f4 is 100%, which indicates the
correct identification of COVID positive patients. This is significant in
limiting the spread of the virus. Additionally, another ML model, as seen to
offer 92% accuracy serves to identify patients who, out of a large group of
patients admitted to the hospital cohort, need immediate critical care support
by estimating the mortality risk of patients from blood parameters. The
instantaneous multimodal nature of our technique offers multiplex and accurate
diagnostic methods, highlighting the effectiveness of telehealth as a simple,
widely available, and low-cost diagnostic solution, even for future pandemics.


  Proteins perform critical processes in all living systems: converting solar
energy into chemical energy, replicating DNA, as the basis of highly performant
materials, sensing and much more. While an incredible range of functionality
has been sampled in nature, it accounts for a tiny fraction of the possible
protein universe. If we could tap into this pool of unexplored protein
structures, we could search for novel proteins with useful properties that we
could apply to tackle the environmental and medical challenges facing humanity.
This is the purpose of protein design.
  Sequence design is an important aspect of protein design, and many successful
methods to do this have been developed. Recently, deep-learning methods that
frame it as a classification problem have emerged as a powerful approach.
Beyond their reported improvement in performance, their primary advantage over
physics-based methods is that the computational burden is shifted from the user
to the developers, thereby increasing accessibility to the design method.
Despite this trend, the tools for assessment and comparison of such models
remain quite generic. The goal of this paper is to both address the timely
problem of evaluation and to shine a spotlight, within the Machine Learning
community, on specific assessment criteria that will accelerate impact.
  We present a carefully curated benchmark set of proteins and propose a number
of standard tests to assess the performance of deep learning based methods. Our
robust benchmark provides biological insight into the behaviour of design
methods, which is essential for evaluating their performance and utility. We
compare five existing models with two novel models for sequence prediction.
Finally, we test the designs produced by these models with AlphaFold2, a
state-of-the-art structure-prediction algorithm, to determine if they are
likely to fold into the intended 3D shapes.


  The information content of symbolic sequences (such as nucleic- or amino acid
sequences, but also neuronal firings or strings of letters) can be calculated
from an ensemble of such sequences, but because information cannot be assigned
to single sequences, we cannot correlate information to other observables
attached to the sequence. Here we show that an information score obtained from
multivariate (multiple-variable) correlations within sequences of a "training"
ensemble can be used to predict observables of out-of-sample sequences with an
accuracy that scales with the complexity of correlations, showing that
functional information emerges from a hierarchy of multi-variable correlations.


  Node embedding is a powerful approach for representing the structural role of
each node in a graph. $\textit{Node2vec}$ is a widely used method for node
embedding that works by exploring the local neighborhoods via biased random
walks on the graph. However, $\textit{node2vec}$ does not consider edge weights
when computing walk biases. This intrinsic limitation prevents
$\textit{node2vec}$ from leveraging all the information in weighted graphs and,
in turn, limits its application to many real-world networks that are weighted
and dense. Here, we naturally extend $\textit{node2vec}$ to
$\textit{node2vec+}$ in a way that accounts for edge weights when calculating
walk biases, but which reduces to $\textit{node2vec}$ in the cases of
unweighted graphs or unbiased walks. We empirically show that
$\textit{node2vec+}$ is more robust to additive noise than $\textit{node2vec}$
in weighted graphs using two synthetic datasets. We also demonstrate that
$\textit{node2vec+}$ significantly outperforms $\textit{node2vec}$ on a
commonly benchmarked multi-label dataset (Wikipedia). Furthermore, we test
$\textit{node2vec+}$ against GCN and GraphSAGE using various challenging gene
classification tasks on two protein-protein interaction networks. Despite some
clear advantages of GCN and GraphSAGE, they show comparable performance with
$\textit{node2vec+}$. Finally, $\textit{node2vec+}$ can be used as a general
approach for generating biased random walks, benefiting all existing methods
built on top of $\textit{node2vec}$. $\textit{Node2vec+}$ is implemented as
part of $\texttt{PecanPy}$, which is available at
https://github.com/krishnanlab/PecanPy .


  The mechanism of SARS-CoV-2 evolution and transmission is elusive and its
understanding, a prerequisite to forecast emerging variants, is of paramount
importance. SARS-CoV-2 evolution is driven by the mechanisms at molecular and
organism scales and regulated by the transmission pathways at the population
scale. In this review, we show that infectivity-based natural selection was
discovered as the mechanism for SARS-CoV-2 evolution and transmission in July
2020. In April 2021, we proved beyond all doubt that such a natural selection
via infectivity-based transmission pathway remained the sole mechanism for
SARS-CoV-2 evolution. However, we reveal that antibody-disruptive co-mutations
[Y449S, N501Y] debuted as a new vaccine-resistant transmission pathway of viral
evolution in highly vaccinated populations a few months ago. Over one year ago,
we foresaw that mutations spike protein RBD residues, 452 and 501, would "have
high chances to mutate into significantly more infectious COVID-19 strains".
Mutations on these residues underpin prevailing SARS-CoV-2 variants Alpha,
Beta, Gamma, Delta, Epsilon, Theta, Kappa, Lambda, and Mu at present and are
expected to be vital to emerging variants. We anticipate that viral evolution
will combine RBD co-mutations at these two sites, creating future variants that
are tens of times more infectious than the original SARS-CoV-2. Additionally,
two complementary transmission pathways of viral evolution: infectivity and
vaccine-resistant, will prolong our battle with COVID-19 for years. We predict
that RBD co-mutation [A411S, L452R, T478K], [L452R, T478K, N501Y], [L452R,
T478K, E484K, N501Y], [K417N, L452R, T478K], and [P384L, K417N, E484K, N501Y]
will have high chances to grow into dominating variants due to their high
infectivity and/or strong ability to break through current vaccines, calling
for the development of new vaccines and antibody therapies.


  Classical swine fever (CSF) in Ecuador is prevalent since 1940, pig farming
represents an important economic and cultural sector. Recently, the National
Veterinary Service (NVS) has implemented individual identification of pigs,
movement control and mandatory vaccination against CSF, looking for a future
eradication. Our aim was to characterise the pig premises according to risk
criteria, analyse the effect of random and targeted strategies to control CSF
and consider the temporal development of the network. We used social network
analysis (SNA), SIRS (susceptible, infected, recovered, susceptible) network
modelling and temporal network analysis. The data set contained 751,003
shipments and 6 million pigs from 2017 to 2019. 165,593 premises were involved:
144,118 farms, 138 industrials, 21,337 traders, and 51 markets. On annual
average, 124,976 premises (75%) received or sent one movement with 1.5 pigs, in
contrast, 166 (0.01%) with 1,372 movements and 11,607 pigs. Simulations
resulted in CSF mean prevalence of 29.93%; Targeted selection strategy reduced
the prevalence to 3.3%, while 24% with random selection. Selection of high-risk
premises in every province was the best strategy using available surveillance
infrastructure. Notably, selecting 10 traders/markets reduced the CSF
prevalence to 4%, evidencing their prime influence over the network. Temporal
analysis showed an overestimation of 38% (causal fidelity) in the number of
transmission paths; The steps to cross the network were 4.3 (average path
length), but take approximately 233 days. In conclusion, surveillance
strategies applied by the NVS could be more efficient to find cases, reduce the
spread of diseases and enable the implementation of risk-based surveillance. To
focus the efforts on target selection of high-risk premises, special attention
should be given to markets/traders which proved similar disease spread
potential.


  Ecological and evolutionary dynamics have been historically regarded as
unfolding at broadly separated timescales. However, these two types of
processes are nowadays well documented to much more tightly than traditionally
assumed, especially in communities of microorganisms. With this motivation in
mind, here we scrutinize recent experimental results showing evidence of rapid
evolution of tolerance by lag in bacterial populations that are periodically
exposed to antibiotic stress in laboratory conditions. In particular, the
distribution of single-cell lag times evolves its average value to
approximately fit the antibiotic-exposure time. Moreover, the distribution
develops right skewed heavy tails, revealing the presence of individuals with
anomalously large lag times. Here, we develop a parsimonious individual-based
model mimicking the actual demographic processes of the experimental setup.
Individuals are characterized by a single phenotypic trait: their intrinsic lag
time, which is transmitted with variation to the progeny. The model (in a
version in which the amplitude of phenotypic variations grows with the
parent\'s lag time) is able to reproduce quite well the key empirical
observations. Furthermore, we develop a general mathematical framework allowing
us to describe with good accuracy the properties of the stochastic model by
means of a macroscopic equation, which generalizes the Crow Kimura equation in
population genetics. From a broader perspective, this work represents a
benchmark for the mathematical framework designed to tackle much more general
eco-evolutionary problems, thus paving the road to further research avenues.


  The potential waning of the vaccination immunity to COVID-19 could pose
threats to public health, as it is tenable that the timing of such waning would
synchronize with the near-complete restoration of normalcy. Should also testing
be relaxed, we might witness a resurgent COVID-19 wave in winter 2021/2022. In
response to this risk, an additional vaccine dose, the booster shot, is being
administered worldwide. In a projected study with an outlook of six months, we
explore the interplay between the rate at which boosters are distributed and
the extent to which testing practices are implemented, using a highly granular
agent-based model tuned on a medium-sized U.S. town. Theoretical projections
indicate that the administration of boosters at the rate at which the vaccine
is currently administered could yield a severe resurgence of the pandemic.
Projections suggest that the peak levels of mid spring 2021 in the vaccination
rate may prevent such a scenario to occur, although exact agreement between
observations and projections should not be expected due to continuously
evolving nature of the pandemics. Our study highlights the importance of
testing, especially to detect asymptomatic individuals in the near future, as
the release of the booster reaches full speed.