  The baseline level of transcription, which is variable and difficult to
quantify, seriously complicates the normalization of comparative transcriptomic
data, but its biological importance remains unappreciated. We show that this
currently neglected ingredient is essential for controlling gene network
multistability and therefore cellular differentiation. Basal expression is
correlated to the degree of chromatin loosening measured by DNA accessibility,
and systematically leads to cellular dedifferentiation as assessed by
transcriptomic signatures, irrespective of the molecular and cellular tools
used. Modeling gene network motifs formally involved in developmental
bifurcations, reveals that the epigenetic landscapes of Waddington are
restructured by the level of non specific expression, such that the attractors
of progenitor and differentiated cells can be mutually exclusive. This
mechanism is universal and holds beyond the particular nature of the genes
involved, provided the multistable circuits are correctly described with
autonomous basal expression. These results explain the relationships long
established between gene expression noise, chromatin decondensation and
cellular dedifferentiation, and highlight how heterochromatin maintenance is
essential for preventing pathological cellular reprogramming, age-related
diseases and cancer.


  Summary: Coarse-grained normal mode analysis (NMA) is a fast computational
technique to study the dynamics of biomolecules. Here we present the
Najmanovich Research Group Toolkit for Elastic Networks (NRGTEN). NRGTEN is a
Python toolkit that implements four different NMA models in addition to popular
and novel metrics to benchmark and measure properties from these models.
Furthermore, the toolkit is available as a public Python package and is easily
extensible for the development or implementation of additional NMA models. The
inclusion of the ENCoM model (Elastic Network Contact Model) developed in our
group within NRGTEN is noteworthy, owing to its account for the specific
chemical nature of atomic interactions. This makes possible some unique
predictions of the effect of mutations, such as on stability (via changes in
vibrational entropy differences), on the transition probability between
different conformational states or on the flexibility profile of the whole
macromolecule/complex (to study allostery and signalling). In addition, all NMA
models can be used to generate conformational ensembles from a starting
structure to aid in protein-protein, protein-ligand or other docking studies
among applications. NRGTEN is freely available via a public Python package
which can be easily installed on any modern machine and includes a detailed
user guide hosted online. Availability and implementation:
https://github.com/gregorpatof/nrgten_package/ Contact:
rafael.najmanovich@umontreal.ca


  This paper introduces new methods to analyze the changing progression of
COVID-19 cases to deaths in different waves of the pandemic. First, an
algorithmic approach partitions each country or state's COVID-19 time series
into a first wave and subsequent period. Next, offsets between case and death
time series are learned for each country via a normalized inner product.
Combining these with additional calculations, we can determine which countries
have most substantially reduced the mortality rate of COVID-19. Finally, our
paper identifies similarities in the trajectories of cases and deaths for
European countries and U.S. states. Our analysis refines the popular conception
that the mortality rate has greatly decreased throughout Europe during its
second wave of COVID-19; instead, we demonstrate substantial heterogeneity
throughout Europe and the U.S. The Netherlands exhibited the largest reduction
of mortality, a factor of 16, followed by Denmark, France, Belgium, and other
Western European countries, greater than both Eastern European countries and
U.S. states. Some structural similarity is observed between Europe and the
United States, in which Northeastern states have been the most successful in
the country. Such analysis may help European countries learn from each other's
experiences and differing successes to develop the best policies to combat
COVID-19 as a collective unit.


  As the COVID-19 outbreak evolves around the world, the World Health
Organization (WHO) and its Member States have been heavily relying on staying
at home and lock down measures to control the spread of the virus. In the last
months, various signs showed that the COVID-19 curve was flattening, but even
the partial lifting of some containment measures (e.g., school closures and
telecommuting) appear to favor a second wave of the disease. The accurate
evaluation of possible countermeasures and their well-timed revocation are
therefore crucial to avoid future waves or reduce their duration. In this
paper, we analyze patient and route data of infected patients from January 20,
2020, to May 31, 2020, collected by the Korean Center for Disease Control &
Prevention (KCDC). This data analysis helps us to characterize patient mobility
patterns and then use this characterization to parameterize simulations to
evaluate different what-if scenarios. Although this is not a definitive model
of how COVID-19 spreads in a population, its usefulness and flexibility are
illustrated using real-world data for exploring virus spread under a variety of
circumstances.


  With growing exchanges of people and merchandise between countries, epidemics
have become an issue of increasing importance and huge amounts of data are
being collected every day. Hence, analyses that were usually run in personal
computers and desktops are no longer feasible. It is now common to run such
tasks in High-performance computing (HPC) environments and/or dedicated
systems. On the other hand we are often dealing in these analyses with graphs
and trees, and running algorithms to find patterns in such structures. Hence,
although graph oriented databases and processing systems can be of much help in
this setting, as far as we know there is no solution relying on these
technologies to address large scale phylogenetic analysis challenges. This
project aims to develop a modular framework for large scale phylogenetic
analysis that exploits such technologies, namely Neo4j. We address this
challenge by proposing and developing a framework which allows representing
large phylogenetic networks and trees, as well as ancillary data, that supports
queries on such data, and allows the deployment of algorithms for
inferring/detecting patterns and pre-computing visualizations, as a Neo4j
plugin. This framework is innovative and brings several advantages to the
phylogenetic analysis, such as by storing the phylogenetic trees will avoid
having to compute them again, and by using multilayer networks will make the
comparison between them more efficient and scalable. Experimental results
showcase that it can be very efficient in the mostly used operations and that
the supported algorithms comply with their time complexity.


  We consider a class of density-dependent branching processes which
generalises exponential, logistic and Gompertz growth. A population begins with
a single individual, grows exponentially initially, and then growth may slow
down as the population size moves towards a carrying capacity. At a time while
the population is still growing superlinearly, a fixed number of individuals
are sampled and their coalescent tree is drawn. Taking the sampling time and
carrying capacity simultaneously to infinity, we prove convergence of the
coalescent tree to a limiting tree which is in a sense universal over our class
of models.


  Many unicellular organisms allocate their key proteins asymmetrically between
the mother and daughter cells, especially in a stressed environment. A recent
theoretical model is able to predict when the asymmetry in segregation of key
proteins enhances the population fitness, extrapolating the solution at two
limits where the segregation is perfectly asymmetric (asymmetry $a$ = 1) and
when the asymmetry is small ($0 \leq a \ll 1$). We generalize the model by
introducing stochasticity and use a transport equation to obtain a
self-consistent equation for the population growth rate and the distribution of
the amount of key proteins. We provide two ways of solving the self-consistent
equation: numerically by updating the solution for the self-consistent equation
iteratively and analytically by expanding moments of the distribution. With
these more powerful tools, we can extend the previous model by Lin et al. to
include stochasticity to the segregation asymmetry. We show the stochastic
model is equivalent to the deterministic one with a modified effective
asymmetry parameter ($a_{\rm eff}$). We discuss the biological implication of
our models and compare with other theoretical models.


  The SIR-compartment model is among the simplest models that describe the
spread of a disease through a population. The model makes the unrealistic
assumption that the population through which the disease is spreading is
well-mixed. Although real populations have heterogeneities in contacts not
represented in the SIR model, it nevertheless well fits real U.S. state data at
multiple points throughout the pandemic. Here we demonstrate mathematically how
closely the simple continuous SIR model approximates a model which includes
heterogeneous contacts, and provide insight onto how one can interpret
parameters gleaned from regression in the context of heterogeneous dynamics.


  Self-supervised representation learning of biological sequence embeddings
alleviates computational resource constraints on downstream tasks while
circumventing expensive experimental label acquisition. However, existing
methods mostly borrow directly from large language models designed for NLP,
rather than with bioinformatics philosophies in mind. Recently, contrastive
mutual information maximization methods have achieved state-of-the-art
representations for ImageNet. In this perspective piece, we discuss how viewing
evolution as natural sequence augmentation and maximizing information across
phylogenetic "noisy channels" is a biologically and theoretically desirable
objective for pretraining encoders. We first provide a review of current
contrastive learning literature, then provide an illustrative example where we
show that contrastive learning using evolutionary augmentation can be used as a
representation learning objective which maximizes the mutual information
between biological sequences and their conserved function, and finally outline
rationale for this approach.


  Graphs representing real world systems may be studied from their underlying
community structure. A community in a network is an intuitive idea for which
there is no consensus on its objective mathematical definition. The most used
metric in order to detect communities is the modularity, though many
disadvantages of this parameter have already been noticed in the literature. In
this work, we present a new approach based on a different metric: the surprise.
Moreover, the biases of different community detection algorithms and benchmark
networks are thoroughly studied, identified and commented about.


  We analytically study the SEIR (Susceptible Exposed Infectious Removed)
epidemic model. The aim is to provide simple analytical expressions for the
peak and asymptotic values and their characteristic times of the populations
affected by the COVID-19 pandemic.


  Herein, we provide estimations for the effective reproduction number $R_e$
for the greater metropolitan area of Athens, Greece during the first wave of
the pandemic (February 26-May 15, 2020). For our calculations, we implemented,
in a comparative approach, the two most widely used methods for the estimation
of $R_e$, that by Wallinga and Teunis and by Cori et al. Data were retrieved
from the national database of SARS-CoV-2 infections in Greece. Our analysis
revealed that the expected value of Re dropped below 1 around March 15, shortly
after the suspension of the operation of educational institutions of all levels
nationwide on March 10, and the closing of all retail activities (cafes, bars,
museums, shopping centres, sports facilities and restaurants) on March 13. On
May 4, the date on which the gradual relaxation of the strict lockdown
commenced, the expected value of $R_e$ was slightly below 1, however with
relatively high levels of uncertainty due to the limited number of notified
cases during this period. Finally, we discuss the limitations and pitfalls of
the methods utilized for the estimation of the $R_e$, highlighting that the
results of such analyses should be considered only as indicative by policy
makers.


  Complex systems, such as life and languages, are governed by principles of
evolution. The analogy and comparison between biology and
linguistics\cite{alphafold2, RoseTTAFold, lang_virus, cell language, faculty1,
language of gene, Protein linguistics, dictionary, Grammar of pro_dom,
complexity, genomics_nlp, InterPro, language modeling, Protein language
modeling} provide a computational foundation for characterizing and analyzing
protein sequences, human corpora, and their evolution. However, no general
mathematical formula has been proposed so far to illuminate the origin of
quantitative hallmarks shared by life and language. Here we show several new
statistical relationships shared by proteins and words, which inspire us to
establish a general mechanism of evolution with explicit formulations that can
incorporate both old and new characteristics. We found natural selection can be
quantified via the entropic formulation by the principle of least effort to
determine the sequence variation that survives in evolution. Besides, the
origin of power law behavior and how changes in the environment stimulate the
emergence of new proteins and words can also be explained via the introduction
of function connection network. Our results demonstrate not only the
correspondence between genetics and linguistics over their different
hierarchies but also new fundamental physical properties for the evolution of
complex adaptive systems. We anticipate our statistical tests can function as
quantitative criteria to examine whether an evolution theory of sequence is
consistent with the regularity of real data. In the meantime, their
correspondence broadens the bridge to exchange existing knowledge, spurs new
interpretations, and opens Pandora's box to release several potentially
revolutionary challenges. For example, does linguistic arbitrariness conflict
with the dogma that structure determines function?


  The COVID-19 disease has forced countries to make a considerable
collaborative effort between scientists and governments to provide indicators
to suitable follow-up the pandemic's consequences. Mathematical modeling plays
a crucial role in quantifying indicators describing diverse aspects of the
pandemic. Consequently, this work aims to develop a clear, efficient, and
reproducible methodology for parameter optimization, whose implementation is
illustrated using data from three representative regions from Chile and a
suitable generalized SIR model together with a fitted positivity rate. Our
results reproduce the general trend of the infected's curve, distinguishing the
reported and real cases. Finally, our methodology is robust, and it allows us
to forecast a second outbreak of COVID-19 and the infection fatality rate of
COVID-19 qualitatively according to the reported dead cases.


  We study the univariate moment problem of piecewise-constant density
functions on the interval $[0,1]$ and its consequences for an inference problem
in population genetics. We show that, up to closure, any collection of $n$
moments is achieved by a step function with at most $n-1$ breakpoints and that
this bound is tight. We use this to show that any point in the $n$th
coalescence manifold in population genetics can be attained by a piecewise
constant population history with at most $n-2$ changes. Both the moment cones
and the coalescence manifold are projected spectrahedra and we describe the
problem of finding a nearest point on them as a semidefinite program.


  The present Health Crisis tests the response of modern science and medicine
to finding treatment for a new COVID-19 disease. The presentation on the world
stage of antivirals such as remdesivir, obeys to the continuous investigation
of biologically active molecules with multiple theoretical, computational and
experimental tools. Diseases such as COVID:19 remind us that research into
active ingredients for therapeutic purposes should cover all available sources,
such as plants. In the present work, in silico tools, specifically docking
study, were used to evaluate the binding and inhibition capacity of an
antiviral such as remdesivir on the NSP-12 protein of SARS-CoV, a polymerase
that is key in the replication of the SARS-COV virus. The results are then
compared with a docking analysis of two natural products (Alpha-Bisabolol and
betalain) with SARS-CoV protein, in order to find more candidates for COVID-19
virus replication inhibitors. in addition to increasing studies that help
explain the specific mechanisms of the SARs-CoV-2 virus, remembering that we
will have to live with the virus for an indefinite time from now on. Finally,
natural products such as betalains may have inhibitory effects of a small order
but in conjunction with other synergistic active ingredients they may increase
their inhibition effect on NSP-12 protein of SARS-CoV.


  In a single-state repeated game, zero-determinant strategies can unilaterally
force functions of the payoffs to take values in particular closed intervals.
When the explicit use of a determinant is absent from the analysis, they are
instead called autocratic. We extend their study to the setting of finite state
games with deterministic transitions. For a given game we find that the
endpoints of the intervals of enforceable values must satisfy fixed point
equations. From these extreme enforceable values we show it is always possible
to construct finite memory strategies to enforce a particular value. An
algorithm is presented which will approximate the enforceable values in each
state. Finally, we present formulas from which the exact solutions can be
recovered from the approximate solution.


  How large ecosystems can create and maintain the remarkable biodiversity we
see in nature is probably one of the biggest open questions in science,
attracting attention from different fields, from Theoretical Ecology to
Mathematics and Physics. In this context, modeling the stable coexistence of
species competing for limited resources is a particularly challenging task.
From a mathematical point of view, coexistence in competitive dynamics can be
achieved when dominance among species forms intransitive loops. However, these
relationships usually lead to species' relative abundances neutrally cycling
without converging to a stable equilibrium. Although in recent years several
mechanisms have been proposed, models able to explain species coexistence in
competitive communities are still limited. Here we identify locality in the
interactions as one of the simplest mechanisms leading to stable species
coexistence. We consider a simplified ecosystem where individuals of each
species lay on a spatial network and interactions are possible only between
nodes within a certain distance. Varying such distance allows to interpolate
between local and global competition. Our results demonstrate, within the scope
of our model, that species coexist reaching a stable equilibrium when two
conditions are met: individuals are embedded in space and can only interact
with other individuals within a short distance. On the contrary, when one of
these ingredients is missing, large oscillations and neutral cycles emerge.


  Network-based models of epidemic spread have become increasingly popular in
recent decades. Despite a rich foundation of such models, few low-dimensional
systems for modeling SIS-type diseases have been proposed that manage to
capture the complex dynamics induced by the network structure. We analyze one
recently introduced model and derive important epidemiological quantities for
the system. We derive the epidemic threshold and analyze the bifurcation that
occurs, and we use asymptotic techniques to derive an approximation for the
endemic equilibrium when it exists. We consider the sensitivity of this
approximation to network parameters, and the implications for disease control
measures are found to be in line with the results of existing studies.


  Two factors that are often ignored but could play a crucial role in the
progression of an infectious disease are the distributions of inherent
susceptibility ($\sigma_{inh}$) and external infectivity ($\iota_{ext}$), in a
given population. While the former is determined by the immunity of an
individual towards a disease, the latter depends on the duration of exposure to
the infection. We model the spatio-temporal propagation of a pandemic using a
generalized SIR (Susceptible-Infected-Removed) model by introducing the
susceptibility and infectivity distributions to understand their combined
effects, which appear to remain inadequately addressed till date. We consider
the coupling between $\sigma_{inh}$ and $\iota_{ext}$ through a new Critical
Infection Parameter (CIP) ($\gamma_c$). We find that the neglect of these
distributions, as in the naive SIR model, results in an overestimation of the
amount of infection in a population, which leads to incorrect (higher)
estimates of the infections required to achieve the herd immunity threshold.
Additionally, we include the effects of seeding of infection in a population by
long-range migration. We solve the resulting master equations by performing
Kinetic Monte Carlo Cellular Automata (KMC-CA) simulations. Importantly, our
simulations can reproduce the multiple infection peak scenario of a pandemic.
The latent interactions between disease migration and the distributions of
susceptibility and infectivity can render the progression a character vastly
different from the naive SIR model. In particular, inclusion of these
additional features renders the problem a character of a living percolating
system where the disease cluster survives by migrating from region to region.


  Motivated by COVID-19, we develop and analyze a simple stochastic model for a
disease spread in human population. We track how the number of infected and
critically ill people develops over time in order to estimate the demand that
is imposed on the hospital system. To keep this demand under control, we
consider a class of simple policies for slowing down and reopening the society
and we compare their efficiency in mitigating the spread of the virus from
several different points of view. We find that in order to avoid overwhelming
of the hospital system, a policy must impose a harsh lockdown or it must react
swiftly (or both). While reacting swiftly is universally beneficial, being
harsh pays off only when the country is patient about reopening and when the
neighboring countries coordinate their mitigation efforts. Our work highlights
the importance of acting decisively when closing down and the importance of
patience and coordination between neighboring countries when reopening.


  We have studied the dynamic evolution of the Covid-19 pandemic in Argentina.
The marked heterogeneity in population density and the very extensive geography
of the country becomes a challenge itself. Standard compartment models fail
when they are implemented in the Argentina case. We extended a previous
successful model to describe the geographical spread of the AH1N1 influenza
epidemic of 2009 in two essential ways: we added a stochastic local mobility
mechanism, and we introduced a new compartment in order to take into account
the isolation of infected asymptomatic detected people. Two fundamental
parameters drive the dynamics: the time elapsed between contagious and
isolation of infected individuals ($\alpha$) and the ratio of people isolated
over the total infected ones ($p$). The evolution is more sensitive to the
$p-$parameter. The model not only reproduces the real data but also predicts
the second wave before the former vanishes. This effect is intrinsic of
extensive countries with heterogeneous population density and interconnection.
The model presented here becomes a good predictor of the effects of public
policies as, for instance, the unavoidable vaccination campaigns starting at
present in the world.


  The mathematical interpretation of interventions for the mitigation of
epidemics and pandemics in the literature often involves finding the optimal
time to initiate an intervention and/or the use of infections to manage impact.
Whilst these methods may work in theory, in order to implement they may require
information which is likely not available whilst one is in the midst of an
epidemic, or they may require impeccable data about infection levels in the
community. In practice, testing and cases data is only as good as the policy of
implementation and the compliance of the individuals, which means that
understanding the levels of infections becomes difficult or complicated from
the data that is provided. In this paper, we aim to develop a different
approach to the mathematical modelling of interventions, not based on
optimality, but based on demand and capacity of local authorities who have to
deal with the epidemic on a day to day basis. In particular, we use data-driven
modelling to calibrate an Susceptible Exposed Infectious Recovered-Died
(SEIR-D) model to infer parameters that depict the dynamics of the epidemic in
a region of the UK. We use the calibrated parameters for forecasting scenarios
and understand, given a maximum capacity of hospital healthcare services, how
the timing of interventions, severity of interventions, and conditions for the
releasing of interventions affect the overall epidemic-picture.


  Enhancer-promoter interactions (EPIs) regulate the expression of specific
genes in cells, and EPIs are important for understanding gene regulation, cell
differentiation and disease mechanisms. EPI identification through the wet
experiments is costly and time-consuming, and computational methods are in
demand. In this paper, we propose a deep neural network-based method EPIHC
based on sequence-derived features and genomic features for the EPI prediction.
EPIHC extracts features from enhancer and promoter sequences respectively using
convolutional neural networks (CNN), and then design a communicative learning
module to captures the communicative information between enhancer and promoter
sequences. EPIHC also take the genomic features of enhancers and promoters into
account. At last, EPIHC combines sequence-derived features and genomic features
to predict EPIs. The computational experiments show that EPIHC outperforms the
existing state-of-the-art EPI prediction methods on the benchmark datasets and
chromosome-split datasets, and the study reveal that the communicative learning
module can bring explicit information about EPIs, which is ignore by CNN.
Moreover, we consider two strategies to improve performances of EPIHC in the
cross-cell line prediction, and experimental results show that EPIHC
constructed on training cell lines exhibit improved performances for the other
cell lines.


  Enzymes have been shown to diffuse faster in the presence of their reactants.
Recently, we revealed new insights into this process of enhanced diffusion
using single-particle tracking (SPT) with total internal reflection
fluorescence (TIRF) microscopy. We found that the mobility of individual
enzymes was enhanced three fold in the presence of the substrate, and the
motion remained Brownian. In this work, we compare different experimental
designs, as well as different data analysis approaches, for studying single
enzyme diffusion. We first tether enzymes directly on supported lipid bilayers
(SLBs) to constrain the diffusion of enzymes to two dimensions. This
experimental design recovers the 3-fold enhancement in enzyme diffusion in the
presence of the substrate, as we observed before. We also simplify our system
by replacing the bulky polymers used in the prior chamber design with a
SLB-coated surface and glycerol. Using this newly-designed SLB/glycerol
chamber, we compare two different analysis approaches for SPT: the mean-squared
displacement (MSD) analysis and the jump-length analysis. We find that the MSD
analysis requires high viscosity and large particles to accurately report the
diffusion coefficient, while jump-length analysis depends less on the viscosity
or size. Furthermore, the SLB-glycerol chamber fails to reproduce the enhanced
diffusion of enzymes because glycerol inhibits enzyme activity.


  For the reaction networks with one-dimensional stoichiometric subspaces, we
show the following results. (1) If the maximum number of positive steady states
is an even number N, then the maximum number of stable positive steady states
is N/2. (2) If the maximum number of positive steady states is an odd number N,
then we provide a condition on the network such that the maximum number of
stable positive steady states is (N-1)/2 if this condition is satisfied, and
this maximum number is (N+1)/2 otherwise.


  We develop and analyze a spatial temporal model of light driven
ecotoxicological processes, motivated by an aquatic predator-prey system of
algae and \textsl{Daphnia} subject to a contaminant. Population dynamics are
driven by light, which is periodic in time and varies with spatial depth. The
existence and uniqueness of spatial and temporal dependent periodic solutions
are shown and analytical functions of the solutions under parameter constraints
are presented. We conduct Turing stability analyses of solutions with respect
to perturbations of initial conditions. Given a perturbation to a periodic
equilibrium state, we show the system will return to this equilibrium state as
long as motility is fast enough and/or the reservoir depth is shallow enough.
Analytical results assume some Dirichlet boundary conditions that match the
periodic equilibrium state, however numerical simulations with more relaxed
boundary conditions capture similar periodic solutions. The work sheds light
onto spatially dependent population dynamics that are driven by periodic
forces, such as light levels.


  Mathematical models describing SARS-CoV-2 dynamics and the corresponding
immune responses in patients with COVID-19 can be critical to evaluate possible
clinical outcomes of antiviral treatments. In this work, based on the concept
of virus spreadability in the host, antiviral effectiveness thresholds are
determined to establish whether or not a treatment will be able to clear the
infection. In addition, the virus dynamic in the host -- including the
time-to-peak and the final monotonically decreasing behavior -- is chracterized
as a function of the treatment initial time. Simulation results, based on nine
real patient data, show the potential clinical benefits of a treatment
classification according to patient critical parameters. This study is aimed at
paving the way for the different antivirals being developed to tackle
SARS-CoV-2.


  As the year 2020 draws to an end, several new strains have been reported for
the SARS-CoV-2 coronavirus, the agent responsible for the COVID-19 pandemic
that has afflicted us all this past year. However, it is difficult to
comprehend the scale, in sequence space, geographical location and time, at
which SARS-CoV-2 mutates and evolves in its human hosts. To get an appreciation
for the rapid evolution of the coronavirus, we built interactive scalable
vector graphics maps that show daily nucleotide variations in genomes from the
six most populated continents compared to that of the initial, ground-zero
SARS-CoV-2 isolate sequenced at the beginning of the year. Availability:
Mutation time maps are available from https://bcgsc.github.io/SARS2/


  In genetic networks, information of relevance to the organism is represented
by the concentrations of transcription factor molecules. In order to extract
this information the cell must effectively "measure"' these concentrations, but
there are physical limits to the precision of these measurements. We explore
this trading between bits of precision in measuring concentration and bits of
relevant information that can be extracted, using the gap gene network in the
early fly embryo as an example. We argue that cells in the embryo can extract
all the available information about their position, but only if the
concentration measurements approach the physical limits to information
capacity. These limits necessitate the observed proliferation of enhancer
elements with sensitivities to combinations of transcription factors, but fine
tuning of the parameters of these multiple enhancers is not required.


  There is growing appreciation that gene function is connected to the dynamic
structure of the chromosome. Here we explore the interplay between
three-dimensional structure and transcriptional activity at the single cell
level. We show that inactive loci are spatially more compact than active ones,
and that within active loci the enhancer driving transcription is closest to
the promoter. On the other hand, even this shortest distance is too long to
support direct physical contact between the enhancer-promoter pair when the
locus is transcriptionally active. Artificial manipulation of genomic
separations between enhancers and the promoter produces changes in physical
distance and transcriptional activity, recapitulating the correlation seen in
wild-type embryos, but disruption of topological domain boundaries has no
effect. Our results suggest a complex interdependence between transcription and
the spatial organization of cis-regulatory elements.


  The exploration of selected single nucleotide polymorphisms (SNPs) to
identify genetic diversity between different sequencing population pools
(Pool-seq) is a fundamental task in genetic research. As underlying sequence
reads and their alignment are error-prone and univariate statistical solutions
only take individual positions of the genome into account, the identification
of selected SNPs remains a challenging process. Deep learning models like
convolutional neural networks (CNNs) are able to consider large input areas in
their decisions. We suggest an unsupervised pipeline to be independent of a
rarely known ground truth. We train a supervised discriminator CNN to
distinguish alignments from different populations and utilize the model for
unsupervised SNP calling by applying explainable artificial intelligence
methods. Our proposed multivariate method is based on two main assumptions: We
assume (i) that instances having a high predictive certainty of being
distinguishable are likely to contain genetic variants, and (ii) that selected
SNPs are located at regions with input features having the highest influence on
the model's decision process. We directly compare our method with statistical
results on two different Pool-seq datasets and show that our solution is able
to extend statistical results.


  Super-resolution microscopy has catalyzed valuable insights into the
sub-cellular, mechanistic details of many different biological processes across
a wide range of cell types. Fluorescence polarization spectroscopy tools have
also enabled important insights into cellular processes through identifying
orientational changes of biological molecules typically at an ensemble level.
Here, we combine these two biophysical methodologies in a single home-made
instrument to enable the simultaneous detection of orthogonal fluorescence
polarization signals from single fluorescent protein molecules used as common
reporters on the localization of proteins in cellular processes. These enable
measurement of spatial location to a super-resolved precision better than the
diffraction-limited optical resolution, as well as estimation of molecular
stoichiometry based on the brightness of individual fluorophores. In this
innovation we have adapted a millisecond timescale microscope used for
single-molecule detection to enable splitting of fluorescence polarization
emissions into two separate imaging channels for s- and p- polarization
signals, which are imaged onto separate halves of the same high sensitivity
back-illuminated CMOS camera detector. We applied this fluorescence
polarization super-resolved imaging modality to a range of test fluorescent
samples relevant to the study of biological processes, including purified
monomeric green fluorescent protein, single combed DNA molecules, and protein
assemblies and complexes from live Escherichia coli and Saccharomyces
cerevisiae cells. Our findings are qualitative but demonstrate promise in
showing how fluorescence polarization and super-resolved localization
microscopy can be combined on the same sample to enable simultaneous
measurements of polarization and stoichiometry of tracked molecular complexes,
as well as the translational diffusion coefficient.


  Bacteria live in environments that are continuously fluctuating and changing.
Exploiting any predictability of such fluctuations can lead to an increased
fitness. On longer timescales bacteria can "learn" the structure of these
fluctuations through evolution. However, on shorter timescales, inferring the
statistics of the environment and acting upon this information would need to be
accomplished by physiological mechanisms. Here, we use a model of metabolism to
show that a simple generalization of a common regulatory motif (end-product
inhibition) is sufficient both for learning continuous-valued features of the
statistical structure of the environment and for translating this information
into predictive behavior; moreover, it accomplishes these tasks near-optimally.
We discuss plausible genetic circuits that could instantiate the mechanism we
describe, including one similar to the architecture of two-component signaling,
and argue that the key ingredients required for such predictive behavior are
readily accessible to bacteria.


  Non-proportional hazards data are routinely encountered in randomized
clinical trials. In such cases, classic Cox proportional hazards model can
suffer from severe power loss, with difficulty in interpretation of the
estimated hazard ratio since the treatment effect varies over time. We propose
CauchyCP, an omnibus test of change-point Cox regression models, to overcome
both challenges while detecting signals of non-proportional hazards patterns.
Extensive simulation studies demonstrate that, compared to existing treatment
comparison tests under non-proportional hazards, the proposed CauchyCP test 1)
controls the type I error better at small $\alpha$ levels ($< 0.01$); 2)
increases the power of detecting time-varying effects; and 3) is more
computationally efficient. The superior performance of CauchyCP is further
illustrated using retrospective analyses of two randomized clinical trial
datasets and a pharmacogenetic biomarker study dataset. The R package
$\textit{CauchyCP}$ is publicly available on CRAN.


  A major strategy to prevent the spread of COVID-19 is the limiting of
in-person contacts. However, this is impractical or impossible for the many
disabled people who do not live in care facilities, but still require
caregivers. We seek to determine which interventions can prevent infections
among disabled people and their caregivers. We simulate transmission with a
model that includes susceptible, exposed, asymptomatic, symptomatically ill,
hospitalized, and removed individuals. The networks on which we simulate
disease spread incorporate heterogeneity in the risks of different types of
interactions, time-dependent lockdown and reopening measures, and contact
distributions for four different groups (caregivers, disabled people, essential
workers, and the general population). We find the probability of becoming
infected is largest for caregivers and second largest for disabled people. Our
analysis of network structure illustrates that caregivers have the largest
modal eigenvector centrality. We find that two interventions --
contact-limiting by all groups and mask-wearing by disabled people and
caregivers -- most reduce the cases among disabled people and caregivers. We
also test which group spreads COVID-19 most readily by seeding infections in a
subset of each group. We find caregivers are the most potent spreaders of
COVID-19, particularly to other caregivers and to disabled people. We test
where to use limited vaccine doses most effectively and find (1) vaccinating
caregivers better protects disabled people than vaccinating the general
population or essential workers and (2) vaccinating caregivers protects
disabled people about as much as vaccinating disabled people themselves. Our
results highlight the potential effectiveness of mask-wearing, contact-limiting
throughout society, and strategic vaccination for limiting the exposure of
disabled people and their caregivers to COVID-19.


  Many dynamical processes of complex systems can be understood as the dynamics
of a group of nodes interacting on a given network structure. However, finding
such interaction structure and node dynamics from time series of node
behaviours is tough. Conventional methods focus on either network structure
inference task or dynamics reconstruction problem, very few of them can work
well on both. This paper proposes a universal framework for reconstructing
network structure and node dynamics at the same time from observed time-series
data of nodes. We use a differentiable Bernoulli sampling process to generate a
candidate network structure, and use neural networks to simulate the node
dynamics based on the candidate network. We then adjust all the parameters with
a stochastic gradient descent algorithm to maximize the likelihood function
defined on the data. The experiments show that our model can recover various
network structures and node dynamics at the same time with high accuracy. It
can also work well on binary, discrete and continuous time-series data, and the
reconstruction results are robust against noise and missing information.


  Both FCM and PCM clustering methods have been widely applied to pattern
recognition and data clustering. Nevertheless, FCM is sensitive to noise and
PCM occasionally generates coincident clusters. PFCM is an extension of the PCM
model by combining FCM and PCM, but this method still suffers from the
weaknesses of PCM and FCM. In the current paper, the weaknesses of the PFCM
algorithm are corrected and the enhanced possibilistic fuzzy c-means (EPFCM)
clustering algorithm is presented. EPFCM can still be sensitive to noise.
Therefore, we propose an interval type-2 enhanced possibilistic fuzzy c-means
(IT2EPFCM) clustering method by utilizing two fuzzifiers $(m_1, m_2)$ for fuzzy
memberships and two fuzzifiers $({\theta}_1, {\theta}_2)$ for possibilistic
typicalities. Our computational results show the superiority of the proposed
approaches compared with several state-of-the-art techniques in the literature.
Finally, the proposed methods are implemented for analyzing microarray gene
expression data.


  Motivated by the recent outbreak of coronavirus (COVID-19), we propose a
stochastic model of epidemic temporal growth and mitigation based on a
time-modulated Hawkes process. The model is sufficiently rich to incorporate
specific characteristics of the novel coronavirus, to capture the impact of
undetected, asymptomatic and super-diffusive individuals, and especially to
take into account time-varying counter-measures and detection efforts. Yet, it
is simple enough to allow scalable and efficient computation of the temporal
evolution of the epidemic, and exploration of what-if scenarios. Compared to
traditional compartmental models, our approach allows a more faithful
description of virus specific features, such as distributions for the time
spent in stages, which is crucial when the time-scale of control (e.g.,
mobility restrictions) is comparable to the lifetime of a single infection. We
apply the model to the first and second wave of COVID-19 in Italy, shedding
light into several effects related to mobility restrictions introduced by the
government, and to the effectiveness of contact tracing and mass testing
performed by the national health service.


  Mina and Andersen, authors of the Perspectives in Science: COVID-19 Testing:
One Size Does Not Fit All have referred to results and adopted conclusions from
recently published governmental report Pavelka et al. "The effectiveness of
population wide, rapid antigen test based screening in reducing SARS-CoV-2
infection prevalence in Slovakia" without critical consideration, and rigorous
verification. We demonstrate that the authors refer to conclusions that are not
supported by experimental data. Further, there is a lack of objective,
independent information and studies regarding the widespread, public testing
program currently in force in the Slovak Republic. We offer an alternative
explanation of observed data as they have been provided by the Slovak Republic
government to fill this information gap. We also provide explanations and
conclusions that more accurately describe viral spread dynamics. Drawing from
available public data and our simple but rigorous analysis, we show that it is
not possible to make clear conclusions about any positive impact of the public
testing program in the Slovak Republic. In particular, it is not possible to
conclude that this testing program forces the curve down for the SARS-CoV-2
virus outbreak. We think that Pavelka et al. did not consider many fundamental
phenomena in their proposed computer simulations and data analysis - in
particular: the complexity of SARS-CoV-2 virus spread. In complex
spatio-temporal dynamical systems, small spatio-temporal fluctuations can
dramatically change the dynamics of virus spreading on large scales.


  As the large amount of sequencing data accumulated in past decades and it is
still accumulating, we need to handle the more and more sequencing data. As the
fast development of the computing technologies, we now can handle a large
amount of data by a reasonable of time using the neural network based model.
This tutorial will introduce the the mathematical model of the single cell
variational inference (scVI), which use the variational auto-encoder (building
on the neural networks) to learn the distribution of the data to gain insights.
It was written for beginners in the simple and intuitive way with many
deduction details to encourage more researchers into this field.


  Segmentation and genome annotation (SAGA) algorithms are widely used to
understand genome activity and gene regulation. These algorithms take as input
epigenomic datasets, such as chromatin immunoprecipitation-sequencing
(ChIP-seq) measurements of histone modifications or transcription factor
binding. They partition the genome and assign a label to each segment such that
positions with the same label exhibit similar patterns of input data. SAGA
algorithms discover categories of activity such as promoters, enhancers, or
parts of genes without prior knowledge of known genomic elements. In this
sense, they generally act in an unsupervised fashion like clustering
algorithms, but with the additional simultaneous function of segmenting the
genome. Here, we review the common methodological framework that underlies
these methods, review variants of and improvements upon this basic framework,
catalogue existing large-scale reference annotations, and discuss the outlook
for future work.


  In this paper, we investigate a novel 3-compartment model of HIV infection of
CD4$^+$ T-cells with a mass action term by including two versions: one baseline
ODE model and one delay-differential equation (DDE) model with a constant
discrete time delay. Similar to various endemic models, the dynamics within the
ODE model is fully determined by the basic reproduction term $R_0$. If $R_0<1$,
the disease-free (zero) equilibrium will be asymptotically stable and the
disease gradually dies out. On the other hand, if $R_0>1$, there exists a
positive equilibrium that is globally/orbitally asymptotically stable within
the interior of a predefined region. To present the incubation time of the
virus, a constant delay term $\tau$ is added, forming a DDE model. In this
model, this time delay (of the transmission between virus and healthy cells)
can destabilize the system, arising periodic solutions through Hopf
bifurcation. Finally, numerical simulations are conducted to illustrate and
verify the results.


  The COVID-19 crisis highlighted the importance of non-medical interventions,
such as testing and isolation of infected individuals, in the control of
epidemics. Here, we show how to minimize testing needs while maintaining the
number of infected individuals below a desired threshold. We find that the
optimal policy is adaptive, with testing rates that depend on the epidemic
state. Additionally, we show that such epidemic state is difficult to infer
with molecular tests alone, which are highly sensitive but have a short
detectability window. Instead, we propose the use of baseline serology testing,
which is less sensitive but detects past infections, for the purpose of state
estimation. Validation of such combined testing approach with a stochastic
model of epidemics shows significant cost savings compared to non-adaptive
testing strategies that are the current standard for COVID-19.


  The transmission of COVID-19 is dependent on social contacts, the rate of
which have varied during the pandemic due to mandated and voluntary social
distancing. Changes in transmission dynamics eventually affect hospital
admissions and we have used this connection in order to model and predict
regional hospital admissions in Sweden during the COVID-19 pandemic. We use an
SEIR-model for each region in Sweden in which the infectivity is assumed to
depend on mobility data in terms of public transport utilisation and mobile
phone usage. The results show that the model can capture the timing of the
first and beginning of the second wave of the pandemic. Further, we show that
for two major regions of Sweden models with public transport data outperform
models using mobile phone usage. The model assumes a three week delay from
disease transmission to hospitalisation which makes it possible to use current
mobility data to predict future admissions.


  Networks with a scale-free degree distribution are widely thought to promote
cooperation in various games. Herein, by studying the well-known prisoner's
dilemma game, we demonstrate that this need not necessarily be true. For the
very same degree sequence and degree distribution, we present a variety of
possible behaviour. We reassess the perceived importance of hubs in a network
towards the maintenance of cooperation. We also reevaluate the dependence of
cooperation on network clustering and assortativity.


  This paper explores the genotype-phenotype relationship. It outlines
conditions under which the dependence of a quantitative trait on the genome
might be predictable, based on measurement of a limited subset of genotypes. It
uses the theory of real-valued Boolean functions in a systematic way to
translate trait data into the Fourier domain. Important trait features, such as
the roughness of the trait landscape or the modularity of a trait have a simple
Fourier interpretation. Ruggedness at a gene location corresponds to high
sensitivity to mutation, while a modular organization of gene activity reduces
such sensitivity.
  Traits where rugged loci are rare will naturally compress gene data in the
Fourier domain, leading to a sparse representation of trait data, concentrated
in identifiable, low-level coefficients. This Fourier representation of a trait
organizes epistasis in a form which is isometric to the trait data. As Fourier
matrices are known to be maximally incoherent with the standard basis, this
permits employing compressive sensing techniques to work from data sets that
are relatively small -- sometimes even of polynomial size -- compared to the
exponentially large sets of possible genomes.
  This theory provides a theoretical underpinning for systematic use of Boolean
function machinery to dissect the dependency of a trait on the genome and
environment.


  Some prokaryotes possess CRISPR-Cas systems that provide adaptive immunity to
viruses guided by DNA segments called spacers acquired from invading phage.
However, the patchy incidence and limited memory breadth of CRISPR-Cas systems
suggest that their fitness benefits are offset by costs. Here, we propose that
cross-reactive CRISPR targeting can lead to heterologous autoimmunity, whereby
foreign spacers guide self-targeting in a spacer-length dependent fashion.
Balancing antiviral defense against autoimmunity predicts a scaling relation
between spacer length and CRISPR repertoire size. We find evidence for this
scaling through comparative analysis of sequenced prokaryotic genomes, and show
that this association also holds at the level of CRISPR types. In contrast, the
scaling is absent in strains with nonfunctional CRISPR loci. Finally, we
demonstrate that stochastic spacer loss can explain variations around the
scaling relation, even between strains of the same species. Our results suggest
that heterologous autoimmunity is a selective factor shaping the evolution of
CRISPR-Cas systems.


  A fractional compartmental mathematical model for the spread of the COVID-19
disease is proposed. Special focus has been done on the transmissibility of
super-spreaders individuals. Numerical simulations are shown for data of
Galicia, Spain, and Portugal. For each region, the order of the Caputo
derivative takes a different value, that is not close to one, showing the
relevance of considering fractional models.


  CRISPR-Cas is an adaptive immune mechanism that has been harnessed for a
variety of genetic engineering applications: the Cas9 protein recognises a
2-5nt DNA motif, known as the PAM, and a programmable crRNA binds a target DNA
sequence that is then cleaved. While off-target activity is undesirable, it
occurs because cross-reactivity was beneficial in the immune system on which
the machinery is based. Here, a stochastic model of the target recognition
reaction was derived to study the specificity of the innate immune mechanism in
bacteria. CRISPR systems with Cas9 proteins that recognised PAMs of varying
lengths were tested on self and phage DNA. The model showed that the energy
associated with PAM binding impacted mismatch tolerance, cleavage probability,
and cleavage time. Small PAMs allowed the CRISPR to balance catching mutant
phages, avoiding self-targeting, and quickly dissociating from critically
non-matching sequences. Additionally, the results revealed a lower tolerance to
mismatches in the PAM and a PAM-proximal region known as the seed, as seen in
experiments. This work illustrates the role that the Cas9 protein has in
dictating the specificity of DNA cleavage that can aid in preventing off-target
activity in biotechnology applications.


  The question of whether biological populations survive or are eventually
driven to extinction has long been examined using mathematical models. In this
work we study population survival or extinction using a stochastic, discrete
lattice-based random walk model where individuals undergo movement, birth and
death events. The discrete model is defined on a two-dimensional hexagonal
lattice with periodic boundary conditions. A key feature of the discrete model
is that crowding effects are introduced by specifying two different crowding
functions that govern how local agent density influences movement events and
birth/death events. The continuum limit description of the discrete model is a
nonlinear reaction-diffusion equation, and we focus on crowding functions that
lead to linear diffusion and a bistable source term that is often associated
with the strong Allee effect. Using both the discrete and continuum modelling
tools we explore the complicated relationship between the long-term survival or
extinction of the population and the initial spatial arrangement of the
population. In particular, we study different spatial arrangements of initial
distributions: (i) a well-mixed initial distribution where the initial density
is independent of position in the domain; (ii) a vertical strip initial
distribution where the initial density is independent of vertical position in
the domain; and, (iii) several forms of two-dimensional initial distributions
where the initial population is distributed in regions with different shapes.
Our results indicate that the shape of the initial spatial distribution of the
population affects extinction of bistable populations. All software required to
solve the discrete and continuum models used in this work are available on
GitHub.


  Inspired by the MaCzek Visual Basic program we provide an R package, RMaCzek,
that produces Czekanowski's diagram. Our package permits any seriation and
distance method the user provides. In this paper we focus on the OLO and
QAP_2SUM methods from the seriation package. We illustrate the possibilities of
our package with three anthropological studies, one socio-economical one and a
phylogenetically motivated simulation study.


  The evolution of epidemiological parameters, such as instantaneous
reproduction number Rt, is important for understanding the transmission
dynamics of infectious diseases. Current estimates of time-varying
epidemiological parameters often face problems such as lagging observations,
averaging inference, and improper quantification of uncertainties. To address
these problems, we propose a Bayesian data assimilation framework for
time-varying parameter estimation. Specifically, this framework is applied to
Rt estimation, resulting in the state-of-the-art DARt system. With DARt, time
misalignment caused by lagging observations is tackled by incorporating
observation delays into the joint inference of infections and Rt; the drawback
of averaging is overcome by instantaneously updating upon new observations and
developing a model selection mechanism that captures abrupt changes; the
uncertainty is quantified and reduced by employing Bayesian smoothing. We
validate the performance of DARt and demonstrate its power in revealing the
transmission dynamics of COVID-19. The proposed approach provides a promising
solution for accurate and timely estimating transmission dynamics from reported
data.


  We have learned to live with many potentially deadly viruses for which there
is no vaccine, no immunity, and no cure. We do not live in constant fear of
these viruses, instead, we have learned how to outsmart them and reduce the
harm they cause. A new mathematical model that combines the spread of diseases
that do not confer immunity together with the evolution of human behaviors
indicates that we may be able to fight new diseases with the same type of
strategy we use to fight viruses like HIV.


  Robustness to genetic or environmental disturbances is often considered as a
key property of living systems. Yet, in spite of being discussed since the
1950s, how robustness emerges from the complexity of genetic architectures and
how it evolves still remains unclear. In particular, whether or not robustness
is independent to various sources of perturbations conditions the range of
adaptive scenarios that can be considered. For instance, selection for
robustness to heritable mutations is likely to be modest and indirect, and its
evolution might result from indirect selection on a pleiotropically-related
character (e.g., homeostasis). Here, I propose to treat various robustness
measurements as quantitative characters, and study theoretically, by
individual-based simulations, their propensity to evolve independently. Based
on a simple evolutionary model of a gene regulatory network, I showed that five
measurements of the robustness of gene expression to genetic or non-genetic
disturbances were substantially correlated. Yet, robustness was mutationally
variable in several dimensions, and robustness components could evolve
differentially under direct selection pressure. Therefore, the fact that the
sensitivity of gene expression to mutations and environmental factors rely on
the same gene networks does not preclude distinct evolutionary histories of
robustness components. This article has been peer-reviewed and recommended by
Peer Community In Evolutionary Biology Posted: 30th March 2022
(https://doi.org/10.24072/pci.evolbiol.100138)


  Intelligence is often discussed in terms of neural networks in the cerebral
cortex, whose evolution has presumably been influenced by Darwinian selection.
Here we present molecular evidence that one of the many kinesin motors, Kif14,
has evolved to exhibit special features in its amino acid sequence that could
have evolved to improve neural networks. The improvement is quantified by
comparison of Kif14 sequences for 12 species. The special feature is level sets
of hydrophobic extrema in water wave profiles based on several hydropathic
scales. The most effective scale is a new one based on fractals, indicative of
approach of globular curvatures to self-organized criticality.


  A minimal model for gene expression, consisting of a switchable promoter
together with the resulting messenger RNA, is equivalent to a Poisson channel
with a binary Markovian input process. Determining its capacity is an
optimization problem with respect to two parameters: the average sojourn times
of the promoter's active (ON) and inactive (OFF) state. An expression for the
mutual information is found by solving the associated filtering problem
analytically on the level of distributions. For fixed peak power, three
bandwidth-like constraints are imposed by lower-bounding (i) the average
sojourn times (ii) the autocorrelation time and (iii) the average time until a
transition. OFF-favoring optima are found for all three constraints, as
commonly encountered for the Poisson channel. In addition, constraint (i)
exhibits a region that favors the ON state, and (iii) shows ON-favoring local
optima.


  The design of antibacterial-releasing coatings or wrapping materials with
controlled drug release capability is a promising strategy to minimise risks of
infection and medical device failure in vivo. Collagen fibres have been
employed as medical device building block, although they still fail to display
controlled release capability, competitive wet-state mechanical properties, and
retained triple helix organisation. We investigated this challenge by pursuing
a multiscale design approach integrating drug encapsulation, in-situ covalent
crosslinking and fibre spinning. By selecting ciprofloxacin (Cip) as a typical
antibacterial drug, wet spinning was selected as a triple helix-friendly route
towards Cip-encapsulated collagen fibres; whilst in situ crosslinking of
fibre-forming triple helices with 1,3 phenylenediacetic acid (Ph) was
hypothesised to yield Ph-Cip {\pi}-{\pi} stacking aromatic interactions and
enable controlled drug release. Higher tensile modulus and strength were
measured in Ph crosslinked fibres compared to state-of-the-art carbodiimide
crosslinked controls. Cip-encapsulated Ph-crosslinked fibres revealed decreased
elongation at break and significantly-enhanced drug retention in vitro with
respect to Cip-free variants and carbodiimide-crosslinked controls,
respectively. This multiscale manufacturing strategy provides new insight
aiming at wet spun collagen triple helices with nanoscale-regulated tensile
properties and drug release capability.


  SARS-CoV-2 is a global challenge due to its ability to spread much faster
than SARS-CoV, which was attributed to the mutations in the receptor binding
domain (RBD). These mutations enhanced the electrostatic interactions.
Recently, a new strain was reported in the UK that includes a mutation (N501Y)
in the RBD, that possibly increases the infection rate. Using Molecular
Dynamics simulations (MD) and Monte Carlo (MC) sampling, we showed that the
N501 mutation enhances the electrostatic interactions due to the formation of a
strong hydrogen bond between SARS-CoV-2-T500 and ACE2-D355 near the mutation
site. In addition, we observed that the electrostatic interactions between the
SARS-CoV-2 and ACE2 in the wild type and the mutant are dominated by
salt-bridges formed between SARS-CoV-2-K417 and ACE2-D30, SARS-CoV-2-K458,
ACE2-E23, and SARS-CoV-2-R403 and ACE2-E37. These interactions contributed more
than 40 % of the total binding energies.


  SARS-CoV-2 is what has caused the COVID-19 pandemic. Early viral infection is
mediated by the SARS-CoV-2 homo-trimeric Spike (S) protein with its receptor
binding domains (RBDs) in the receptor-accessible state. We performed molecular
dynamics simulation on the S protein with a focus on the function of its
N-terminal domains (NTDs). Our study reveals that the NTD acts as a "wedge" and
plays a crucial regulatory role in the conformational changes of the S protein.
The complete RBD structural transition is allowed only when the neighboring NTD
that typically prohibits the RBD's movements as a wedge detaches and swings
away. Based on this NTD "wedge" model, we propose that the NTD-RBD interface
should be a potential drug target.


  In unicellular organisms such as bacteria and in most viruses, mutations
mainly occur during reproduction. Thus, genotypes with a high birth rate should
have a higher mutation rate. However, standard models of asexual adaptation
such as the 'replicator-mutator equation' often neglect this effect. In this
study, we investigate the emergence of a positive dependence between the birth
rate and the mutation rate in models of asexual adaptation and the consequences
of this dependence. We show that it emerges naturally at the population scale,
based on a large population limit of a stochastic timecontinuous
individual-based model with elementary assumptions. We derive a
reaction-diffusion framework that describes the evolutionary trajectories and
steady states in the presence of this dependence. When this model is coupled
with a phenotype to fitness landscape with two optima, one for birth, the other
one for survival, a new trade-off arises in the population. Compared to the
standard approach with a constant mutation rate, the symmetry between birth and
survival is broken. Our analytical results and numerical simulations show that
the trajectories of mean phenotype, mean fitness and the stationary phenotype
distribution are in sharp contrast with those displayed for the standard model.
Here, we obtain trajectories of adaptation where the mean phenotype of the
population is initially attracted by the birth optimum, but eventually
converges to the survival optimum, following a hook-shaped curve which
illustrates the antagonistic effects of mutation on adaptation.


  In this work, we develop an inmate population model with a sentencing length
structure. The sentence length structure of new inmates represents the problem
data and can usually be estimated from the histograms corresponding to the
conviction times that are sentenced in a given population. We obtain a
transport equation, typically known as the McKendrick equation, the homogenous
version of which is included in population models with age structures. Using
this equation, we compute the inmate population and entry/exit rates in
equilibrium, which are the values to consider in the design of a penitentiary
system. With data from the Chilean penitentiary system, we illustrate how to
perform these computations. In classifying the inmate population into two
groups of sentence lengths (short and long), we incorporate the SIS
(susceptible-infected-susceptible) epidemiological model, which considers the
entry of infective individuals. We show that a failure to consider the
structure of the sentence lengths -- as is common in epidemiological models
developed for inmate populations -- for prevalences of new inmates below a
certain threshold induces an underestimation of the prevalence in the prison
population at steady state. The threshold depends on the basic reproduction
number associated with the nonstructured SIS model with no entry of new
inmates. We illustrate our findings with analytical and numerical examples for
different distributions of sentencing lengths.


  As the major target of many vaccines and neutralizing antibodies against
SARS-CoV-2, the spike (S) protein is observed to mutate over time. In this
paper, we present statistical approaches to tackle some challenges associated
with the analysis of S-protein data. We build a Bayesian hierarchical model to
study the temporal and spatial evolution of S-protein sequences, after grouping
the sequences into representative clusters. We then apply sampling methods to
investigate possible changes to the S-protein's 3-D structure as a result of
commonly observed mutations. While the increasing spread of D614G variants has
been noted in other research, our results also show that the co-occurring
mutations of D614G together with S477N or A222V may spread even more rapidly,
as quantified by our model estimates.


  Cell division times in microbial populations display significant
fluctuations, that impact the population growth rate in a non-trivial way. If
fluctuations are uncorrelated among different cells, the population growth rate
is predicted by the Euler-Lotka equation, which is a classic result in
mathematical biology. However, cell division times can be significantly
correlated, due to physical properties of cells that are passed through
generations. In this paper, we derive an equation remarkably similar to the
Euler-Lotka equation which is valid in the presence of correlations. Our exact
result is based on large deviation theory and does not require particularly
strong assumptions on the underlying dynamics. We apply our theory to a
phenomenological model of bacterial cell division in E.coli and to experimental
data. We find that the discrepancy between the growth rate predicted by the
Euler-Lotka equation and our generalized version is relatively small, but large
enough to be measurable by our approach.


  In this paper, we formulate the space-dependent variable-order fractional
master equation to model clustering of particles, organelles, inside living
cells. We find its solution in the long time limit describing non-uniform
distribution due to a space dependent fractional exponent. In the continuous
space limit, the solution of this fractional master equation is found to be
exactly the same as the space-dependent variable-order fractional diffusion
equation. In addition, we show that the clustering of lysosomes, an essential
organelle for healthy functioning of mammalian cells, exhibit space-dependent
fractional exponents. Furthermore, we demonstrate that the non-uniform
distribution of lysosomes in living cells is accurately described by the
asymptotic solution of the space-dependent variable-order fractional master
equation. Finally, Monte Carlo simulations of the fractional master equation
validate our analytical solution.


  Computational inference of dated evolutionary histories relies upon various
hypotheses about RNA, DNA, and protein sequence mutation rates. Using mutation
rates to infer these dated histories is referred to as molecular clock
assumption. Coalescent theory is a popular class of evolutionary models that
implements the molecular clock hypothesis to facilitate computational inference
of dated phylogenies. Cancer and virus evolution are two areas where these
methods are particularly important.
  Methodologically, phylogenetic inference methods require a tree space over
which the inference is performed, and geometry of this space plays an important
role in statistical and computational aspects of tree inference algorithms. It
has recently been shown that molecular clock, and hence coalescent, trees
possess a unique geometry, different from that of classical phylogenetic tree
spaces which do not model mutation rates.
  Here we introduce and study a space of discrete coalescent trees, that is, we
assume that time is discrete, which is inevitable in many computational
formalisations. We establish several geometrical properties of the space and
show how these properties impact various algorithms used in phylogenetic
analyses. Our tree space is a discretisation of a known time tree space, called
t-space, and hence our results can be used to approximate solutions to various
open problems in t-space. Our tree space is also a generalisation of another
known trees space, called the ranked nearest neighbour interchange space, hence
our advances in this paper imply new and generalise existing results about
ranked trees.


  Corynebacterium glutamicum is a Gram-positive bacterium found in soil where
the condition changes demand plasticity of the regulatory machinery. The study
of such machinery at the global scale has been challenged by the lack of data
integration. Here, we report three regulatory network models for C. glutamicum:
strong (3040 interactions) constructed solely with regulations previously
supported by directed experiments; all evidence (4665 interactions) containing
the strong network, regulations previously supported by non-directed
experiments, and protein-protein interactions with a direct effect on gene
transcription; and sRNA (5222 interactions) containing the all evidence network
and sRNA-mediated regulations. Compared to the previous version (2018), the
strong and all evidence networks increased by 75 and 1225 interactions,
respectively. We analyzed the system-level components of the three networks to
identify how they differ and compared their structures against those for the
networks of more than 40 species. The inclusion of the sRNAs regulations
changed the proportions of the system-level components and increased the number
of modules but decreased their size. The C. glutamicum regulatory structure
contrasted with other bacterial regulatory networks. Finally, we used the
strong networks of three model organisms to provide insights and future
directions of the C. glutamicum regulatory network characterization.


  Quantitatively predictive models of biomolecular circuits are important tools
for the design of synthetic biology and molecular communication circuits. The
information content of typical time-lapse single-cell data for the inference of
kinetic parameters is not only limited by measurement uncertainty and intrinsic
stochasticity, but also by the employed perturbations. Novel microfluidic
devices enable the synthesis of temporal chemical concentration profiles. The
informativeness of a perturbation can be quantified based on mutual
information. We propose an approximate method to perform optimal experimental
design of such perturbation profiles. To estimate the mutual information we
perform a multivariate log-normal approximation of the joint distribution over
parameters and observations and scan the design space using Metropolis-Hastings
sampling. The method is demonstrated by finding optimal perturbation sequences
for synthetic case studies on a gene expression model with varying reporter
characteristics.


  We present a simple enzymatic system that is capable of a biphasic response
under competitive inhibition. This is arguably the simplest system that can be
said to be hormetic


  In this paper, we propose a two-group SIR epidemic model to simulate the
outcome of stay-at-home policy and wearing face masks during the first COVID-19
epidemic wave in the United States. Based on our proposed model, we further use
the optimal control approach (with the objective of minimizing total deaths) to
find the optimal dynamical distribution of face masks between the healthcare
workers and the general public. It is not surprising that all the face masks
should be solely reserved for the healthcare workers if the supply is short.
However, when the supply is indeed sufficient, our numerical study indicates
that the general public should share a large portion of face masks at the
beginning of an epidemic wave to dramatically reduce the death toll. This
interesting result partially contradicts with the guideline advised by the US
Surgeon General and the Centers for Disease Control and Prevention (CDC) in
March 2020. The optimality of this sounding CDC guideline highly depends on the
supply level of face masks that changes frequently, and hence it should be
adjusted according to the supply of face masks.


  Enzymes and proteins are live driven biochemicals, which has a dramatic
impact over the environment, in which it is active. So, therefore, it is highly
looked-for to build such a robust and highly accurate automatic and
computational model to accurately predict enzymes nature. In this study, a
novel split amino acid composition model named piSAAC is proposed. In this
model, protein sequence is discretized in equal and balanced terminus to fully
evaluate the intrinsic correlation properties of the sequence. Several
state-of-the-art algorithms have been employed to evaluate the proposed model.
A 10-folds cross-validation evaluation is used for finding out the authenticity
and robust-ness of the model using different statistical measures e.g.
Accuracy, sensitivity, specificity, F-measure and area un-der ROC curve. The
experimental results show that, probabilistic neural network algorithm with
piSAAC feature extraction yields an accuracy of 98.01%, sensitivity of 97.12%,
specificity of 95.87%, f-measure of 0.9812and AUC 0.95812, over dataset S1,
accuracy of 97.85%, sensitivity of 97.54%, specificity of 96.24%, f-measure of
0.9774 and AUC 0.9803 over dataset S2. Evident from these excellent empirical
results, the proposed model would be a very useful tool for academic research
and drug designing related application areas.


  We present a model of the evolution of control systems in a genome under
environmental constraints. The model conceptually follows the Jacob and Monod
model of gene control. Genes contain control elements which respond to the
internal state of the cell as well as the environment to control expression of
a coding region. Control and coding regions evolve to maximize a fitness
function between expressed coding sequences and the environment. 118 runs of
the model run to an average of 1.4 x 10^6 `generations' each with a range of
starting parameters probed the conditions under which genomes evolved a
`default style' of control. Unexpectedly, the control logic that evolved was
not significantly correlated to the complexity of the environment. Genetic
logic was strongly correlated with genome complexity and with the fraction of
genes active in the cell at any one time. More complex genomes correlated with
the evolution of genetic controls in which genes were active (`default on'),
and a low fraction of genes being expressed correlated with a genetic logic in
which genes were biased to being inactive unless positively activated (`default
off' logic). We discuss how this might relate to the evolution of the complex
eukaryotic genome, which operates in a `default off' mode.


  Identifying groups that share common features among datasets through
clustering analysis is a typical problem in many fields of science,
particularly in post-omics and systems biology research. In respect of this,
quantifying how a measure can cluster or organize intrinsic groups is important
since currently there is no statistical evaluation of how ordered is, or how
much noise is embedded in the resulting clustered vector. Many of the
literature focuses on how well the clustering algorithm orders the data, with
several measures regarding external and internal statistical measures; but none
measure has been developed to statistically quantify the noise in an arranged
vector posterior a clustering algorithm, i.e., how much of the clustering is
due to randomness. Here, we present a quantitative methodology, based on
autocorrelation, to assess this problem.


  This paper presents an analysis of the publication of datasets collected via
Google Dataset Search, specialized in families of RNA viruses, whose
terminology was obtained from the National Cancer Institute (NCI) thesaurus
developed by the US Department of Health and Human Services. The objective is
to determine the scope and reuse capacity of the available data, determine the
number of datasets and their free access, the proportion in reusable download
formats, the main providers, their publication chronology, and to verify their
scientific provenance. On the other hand, we also define possible relationships
between the publication of datasets and the main pandemics that have occurred
during the last 10 years. The results obtained highlight that only 52% of the
datasets are related to scientific research, while an even smaller fraction
(15%) are reusable. There is also an upward trend in the publication of
datasets, especially related to the impact of the main epidemics, as clearly
confirmed for the Ebola virus, Zika, SARS-CoV, H1N1, H1N5, and especially the
SARS-CoV-2 coronavirus. Finally, it is observed that the search engine has not
yet implemented adequate methods for filtering and monitoring the datasets.
These results reveal some of the difficulties facing open science in the
dataset field.


  The COVID-19 epidemic is the last of a long list of pandemics that have
affected humankind in the last century. In this paper, we propose a novel
mathematical epidemiological model named SUIHTER from the names of the seven
compartments that it comprises: susceptible uninfected individuals (S),
undetected (both asymptomatic and symptomatic) infected (U), isolated (I),
hospitalized (H), threatened (T), extinct (E), and recovered (R). A suitable
parameter calibration that is based on the combined use of least squares method
and Markov Chain Monte Carlo (MCMC) method is proposed with the aim of
reproducing the past history of the epidemic in Italy, surfaced in late
February and still ongoing to date, and of validating SUIHTER in terms of its
predicting capabilities. A distinctive feature of the new model is that it
allows a one-to-one calibration strategy between the model compartments and the
data that are daily made available from the Italian Civil Protection. The new
model is then applied to the analysis of the Italian epidemic with emphasis on
the second outbreak emerged in Fall 2020. In particular, we show that the
epidemiological model SUIHTER can be suitably used in a predictive manner to
perform scenario analysis at national level.


  We derive the full kinetic equations describing the evolution of the
probability density distribution for a structured population such as cells
distributed according to their ages and sizes. The kinetic equations for such a
"sizer-timer" model incorporates both demographic and individual cell growth
rate stochasticities. Averages taken over the densities obeying the kinetic
equations can be used to generate a second order PDE that incorporates the
growth rate stochasticity. On the other hand, marginalizing over the densities
yields a modified birth-death process that shows how age and size influence
demographic stochasticity. Our kinetic framework is thus a more complete model
that subsumes both the deterministic PDE and birth-death master equation
representations for structured populations.


  Evolutionary game theory assumes that players replicate a highly scored
player's strategy through genetic inheritance. However, when learning occurs
culturally, it is often difficult to recognize someone's strategy just by
observing the behaviour. In this work, we consider players with memory-one
stochastic strategies in the iterated prisoner's dilemma, with an assumption
that they cannot directly access each other's strategy but only observe the
actual moves for a certain number of rounds. Based on the observation, the
observer has to infer the resident strategy in a Bayesian way and chooses his
or her own strategy accordingly. By examining the best-response relations, we
argue that players can escape from full defection into a cooperative
equilibrium supported by Win-Stay-Lose-Shift in a self-confirming manner,
provided that the cost of cooperation is low and the observational learning
supplies sufficiently large uncertainty.


  Sulfur mustard (SM), a chemical warfare agent, is a strong alkylating
compound that readily reacts with numerous biomolecules. The goal of the
present work was to define and validate new biomarkers of exposure to SM that
could be easily accessible in urine or plasma. Because investigations using SM
are prohibited by the Organization for the Prohibition of Chemical Weapons, we
worked with 2-chloroethyl ethyl sulfide (CEES), a monofunctional analog of SM.
We developed an ultra-high-pressure liquid chromatography - tandem mass
spectrometry approach (UHPLC-MS/MS) to the conjugate of CEES to glutathione and
two of its metabolites, the cysteine and the N-acetyl-cysteine conjugates. The
N7-guanine adduct of CEES (N7Gua-CEES) was also targeted. After synthesizing
the specific biomarkers, a solid phase extraction protocol and a UHPLC-MS/MS
method with isotopic dilution were optimized. We were able to quantify
N7Gua-CEES in the DNA of HaCaT keratinocytes and of explants of human skin
exposed to CEES. N7Gua-CEES was also detected in the culture medium of these
two models, together with the glutathione and the cysteine conjugates. In
contrast, the N-acetyl-cysteine conjugate was not detected. The method was then
applied to plasma from mice cutaneously exposed to CEES. All four markers could
be detected. Our present results thus validate both the analytical technique
and the biological relevance of new, easily quantifiable biomarkers of exposure
to CEES. Because CEES behaves very similarly to SM, the results are promising
for application to this toxic of interest.


  Metabolite structure identification has become the major bottleneck of the
mass spectrometry based metabolomics research. Till now, number of mass spectra
databases and search algorithms have been developed to address this issue.
However, two critical problems still exist: the low chemical component record
coverage in databases and significant MS/MS spectra variations related to
experiment equipment and parameter settings. In this work, we considered the
molecule fragment as basic building blocks of the metabolic components which
had relatively consistent signatures in MS/MS spectra. And from a bottom-up
point of view, we built a fragment centered database, MSFragDB, by reorganizing
the data from the Human Metabolome Database (HMDB) and developed an
intensity-free searching algorithm to search and rank the most relative
metabolite according to the users' input. We also proposed the concept of
fragment network, a graph structure that encoded the relationship between the
molecule fragments to find close motif that indicated a specific chemical
structure. Although based on the same dataset as the HMDB, validation results
implied that the MSFragDB had a higher hit ratio and furthermore, estimated
possible taxonomy that a query spectrum belongs to when the corresponding
chemical component was missing in the database. Aid by the Fragment Network,
the MSFragDB was also proved to be able to estimate the right structure while
the MS/MS spectrum suffers from the precursor-contamination. The strategy
proposed is general and can be adopted in existing databases. We believe
MSFragDB and Fragment Network can improve the performance of structure
identification with existing data. The beta version of the database is freely
available at www.xrzhanglab.com/msfragdb/.


  We developed a model and a software package for stochastic simulations of
transmission of COVID-19 and other similar infectious diseases, that takes into
account contact network structures and geographical distribution of population
density, detailed up to a level of location of individuals. Our analysis
framework includes a surrogate model optimization process for quick fitting of
the model's parameters to the observed epidemic curves for cases,
hospitalizations and deaths. This set of instruments (the model, the simulation
code, and the optimizer) is a useful tool for policymakers and epidemic
response teams who can use it to forecast epidemic development scenarios in
local environments (on the scale from towns to large countries) and design
optimal response strategies. The simulation code also includes a geospatial
visualization subsystem, presenting detailed views of epidemic scenarios
directly on population density maps. We used the developed framework to draw
predictions for COVID-19 spreading in the canton of Geneva, Switzerland.


  Profiling of whole transcriptomes has become a cornerstone of molecular
biology and an invaluable tool for the characterization of clinical phenotypes
and the identification of disease subtypes. Analyses of these data are becoming
ever more sophisticated as we move beyond simple comparisons to consider
networks of higher-order interactions and associations. Gene regulatory
networks model the regulatory relationships of transcription factors and genes
and have allowed the identification of differentially regulated processes in
disease systems. In this perspective we discuss gene targeting scores, which
measure changes in inferred regulatory network interactions, and their use in
identifying disease-relevant processes. In addition, we present an example
analysis or pancreatic ductal adenocarcinoma demonstrating the power of gene
targeting scores to identify differential processes between complex phenotypes;
processes which would have been missed by only performing differential
expression analysis. This example demonstrates that gene targeting scores are
an invaluable addition to gene expression analysis in the characterization of
diseases and other complex phenotypes.


  Mendelian randomization is the use of genetic variants to assess the
existence of a causal relationship between a risk factor and an outcome of
interest. Here, we focus on two-sample summary-data Mendelian randomization
analyses with many correlated variants from a single gene region, and
particularly on cis-Mendelian randomization studies which use protein
expression as a risk factor. Such studies must rely on a small, curated set of
variants from the studied region; using all variants in the region requires
inverting an ill-conditioned genetic correlation matrix and results in
numerically unstable causal effect estimates. We review methods for variable
selection and estimation in cis-Mendelian randomization with summary-level
data, ranging from stepwise pruning and conditional analysis to principal
components analysis, factor analysis and Bayesian variable selection. In a
simulation study, we show that the various methods have a comparable
performance in analyses with large sample sizes and strong genetic instruments.
However, when weak instrument bias is suspected, factor analysis and Bayesian
variable selection produce more reliable inferences than simple pruning
approaches, which are often used in practice. We conclude by examining two case
studies, assessing the effects of LDL-cholesterol and serum testosterone on
coronary heart disease risk using variants in the HMGCR and SHBG gene regions
respectively.


  In a few years, space telescopes will investigate our Galaxy to detect
evidence of life, mainly by observing rocky planets. In the last decade, the
observation of exoplanet atmospheres and the theoretical works on biosignature
gasses have experienced a considerable acceleration. The~most attractive
feature of the realm of exoplanets is that 40\% of M dwarfs host super-Earths
with a minimum mass between 1 and 30 Earth masses, orbital periods shorter than
50 days, and radii between those of the Earth and Neptune (1--3.8 R$_\oplus$).
Moreover, the recent finding of cyanobacteria able to use far-red (FR) light
for oxygenic photosynthesis due to the synthesis of chlorophylls $d$ and $f$,
extending in vivo light absorption up to 750\ nm, suggests the possibility of
exotic photosynthesis in planets around M dwarfs. Using innovative laboratory
instrumentation, we exposed different cyanobacteria to an M dwarf star
simulated irradiation, comparing their responses to those under solar and FR
simulated lights.~As expected, in FR light, only the cyanobacteria able to
synthesize chlorophyll $d$ and $f$ could grow. Surprisingly, all strains, both
able or unable to use FR light, grew and photosynthesized under the M dwarf
generated spectrum in a similar way to the solar light and much more
efficiently than under the FR one. Our findings highlight the importance of
simulating both the visible and FR light components of an M dwarf spectrum to
correctly evaluate the photosynthetic performances of oxygenic organisms
exposed under such an exotic light~condition.


  Increasing number in global COVID-19 cases demands for mathematical model to
analyze the interaction between the virus dynamics and the response of innate
and adaptive immunity. Here, based on the assumption of a weak and delayed
response of the innate and adaptive immunity in SARS-CoV-2 infection, we
constructed a mathematical model to describe the dynamic processes of immune
system. Integrating theoretical results with clinical COVID-19 patients' data,
we classified the COVID-19 development processes into three typical modes of
immune responses, correlated with the clinical classification of mild &
moderate, severe and critical patients. We found that the immune efficacy (the
ability of host to clear virus and kill infected cells) and the lymphocyte
supply (the abundance and pool of na\"ive T and B cell) play important roles in
the dynamic process and determine the clinical outcome, especially for the
severe and critical patients. Furthermore, we put forward possible treatment
strategies for the three typical modes of immune response. We hope our results
can help to understand the dynamical mechanism of the immune response against
SARS-CoV-2 infection, and to be useful for the treatment strategies and vaccine
design.


  We introduce and study the dynamics of an \emph{immortal} critical branching
process. In the classic, critical branching process, particles give birth to a
single offspring or die at the same rates. Even though the average population
is constant in time, the ultimate fate of the population is extinction. We
augment this branching process with immortality by positing that either: (a) a
single particle cannot die, or (b) there exists an immortal stem cell that
gives birth to ordinary cells that can subsequently undergo critical branching.
We discuss the new dynamical aspects of this immortal branching process.


  Collective migration -- the directed, coordinated motion of many
self-propelled agents -- is a fascinating emergent behavior exhibited by active
matter that has key functional implications for biological systems. Extensive
studies have elucidated the different ways in which this phenomenon may arise.
Nevertheless, how collective migration can persist when a population is
confronted with perturbations, which inevitably arise in complex settings, is
poorly understood. Here, by combining experiments and simulations, we describe
a mechanism by which collectively migrating populations smooth out large-scale
perturbations in their overall morphology, enabling their constituents to
continue to migrate together. We focus on the canonical example of chemotactic
migration of Escherichia coli, in which fronts of cells move via directed
motion, or chemotaxis, in response to a self-generated nutrient gradient. We
identify two distinct modes in which chemotaxis influences the morphology of
the population: cells in different locations along a front migrate at different
velocities due to spatial variations in (i) the local nutrient gradient and in
(ii) the ability of cells to sense and respond to the local nutrient gradient.
While the first mode is destabilizing, the second mode is stabilizing and
dominates, ultimately driving smoothing of the overall population and enabling
continued collective migration. This process is autonomous, arising without any
external intervention; instead, it is a population-scale consequence of the
manner in which individual cells transduce external signals. Our findings thus
provide insights to predict, and potentially control, the collective migration
and morphology of cell populations and diverse other forms of active matter.


  The COVID-19 pandemic has, worldwide and up to December 2020, caused over 1.7
million deaths, and put the world's most advanced healthcare systems under
heavy stress. In many countries, drastic restrictive measures adopted by
political authorities, such as national lockdowns, have not prevented the
outbreak of new pandemic's waves. In this article, we propose an integrated
detection-estimation-forecasting framework that, using publicly available data,
is designed to: (i) learn relevant features of the pandemic (e.g., the
infection rate); (ii) detect as quickly as possible the onset (or the
termination) of an exponential growth of the contagion; and (iii) reliably
forecast the pandemic evolution. The proposed solution is validated by
analyzing the COVID-19 second and third waves in the USA.


  The eastern oyster is a keystone species and ecosystem engineer. However,
restoration efforts of wild oysters are often unsuccessful, in that they do not
produce a robust population of oysters that are able to successfully reproduce.
Furthermore, the dynamics of wild oyster fertilization is not yet well
understood. Through conducting an experiment predicated on quantifying the
influence of elementary aspects of fertilization kinetics--sperm concentration,
gamete age, and success rate--we found that, as stochastic as the mating
process may seem, there are correlations which fundamentally serve as the
framework for assessing long-term sustainability, reef structure, and
hydrodynamic parameters in relation to fertilization. We then focused on
mathematically defining a procedure which simulated a concentration
distribution of a single sperm and egg release where there existed conditions
necessary for breeding to take place. We found a very significant impact of
both gamete age and sperm concentration on fertilization rate ($p < 0.0001$).
Our hydrodynamic model demonstrates that distance can also drastically
influence broadcast spawning. This could be used as a foundation for developing
a flexible model for wild oyster fertilization based on placement, initial
seawater conditions, and size of the starting population. The results of this
research could be implemented into a more user-friendly program which would
accept multiple variables as inputs and output the probability of fertilization
given arbitrary values. By accounting for environmental deviations, this
generalization would increase its compatibility with the public and actualize
the project's intended purpose: enhance the planning of oyster reef restoration
projects.