  A variety of missense mutations and a stop mutation in the gene coding for
transmembrane protein 240 (TMEM240) have been reported to be the causative
mutations of spinocerebellar ataxia 21 (SCA21). We aimed to investigate the
expression of TMEM240 protein in mouse brain at the tissue, cellular, and
subcellular levels. Immunofluorescence labeling showed TMEM240 to be expressed
in various areas of the brain, with the highest levels in the hippocampus,
isocortex, and cerebellum. In the cerebellum, TMEM240 was detected in the deep
nuclei and the cerebellar cortex. The protein was expressed in all three layers
of the cortex and various cerebellar neurons. TMEM240 was localized to
climbing, mossy, and parallel fiber afferents projecting to Purkinje cells, as
shown by coimmunostaining with VGLUT1 and VGLUT2. Co-immunostaining with
synaptophysin, post-synaptic fractionation, and confirmatory electron
microscopy showed TMEM240 to be localized to the post-synaptic side of synapses
near the Purkinje-cell soma. Similar results were obtained in human cerebellar
sections. These data suggest that TMEM240 may be involved in the organization
of the cerebellar network, particularly in synaptic inputs converging on
Purkinje cells. This study is the first to describe TMEM240 expression in the
normal mouse brain.


  Homozygous mutations in the progranulin gene (GRN) are associated with
neuronal ceroid lipofuscinosis 11 (CLN11), a rare lysosomal-storage disorder
characterized by cerebellar ataxia, seizures, retinitis pigmentosa, and
cognitive disorders, usually beginning between 13 and 25 years of age. This is
a rare condition, previously reported in only four families. In contrast,
heterozygous GRN mutations are a major cause of frontotemporal dementia
associated with neuronal cytoplasmic TDP-43 inclusions. We identified
homozygous GRN mutations in six new patients. The phenotypic spectrum is much
broader than previously reported, with two remarkably distinct presentations,
depending on the age of onset. A childhood/juvenile form is characterized by
classical CLN11 symptoms at an early age at onset. Unexpectedly, other
homozygous patients presented a distinct delayed phenotype of frontotemporal
dementia and parkinsonism after 50 years; none had epilepsy or cerebellar
ataxia. Another major finding of this study is that all GRN mutations may not
have the same impact on progranulin protein synthesis. A hypomorphic effect of
some mutations is supported by the presence of residual levels of plasma
progranulin and low levels of normal transcript detected in one case with a
homozygous splice-site mutation and late onset frontotemporal dementia. This is
a new critical finding that must be considered in therapeutic trials based on
replacement strategies. The first neuropathological study in a homozygous
carrier provides new insights into the pathological mechanisms of the disease.
Hallmarks of neuronal ceroid lipofuscinosis were present. The absence of TDP-43
cytoplasmic inclusions markedly differs from observations of heterozygous
mutations, suggesting a pathological shift between lysosomal and TDP-43
pathologies depending on the mono or bi-allelic status. An intriguing
observation was the loss of normal TDP-43 staining in the nucleus of some
neurons, which could be the first stage of the TDP-43 pathological process
preceding the formation of typical cytoplasmic inclusions. Finally, this study
has important implications for genetic counselling and molecular diagnosis.
Semi-dominant inheritance of GRN mutations implies that specific genetic
counseling should be delivered to children and parents of CLN11 patients, as
they are heterozygous carriers with a high risk of developing dementia. More
broadly, this study illustrates the fact that genetic variants can lead to
different phenotypes according to their mono- or bi-allelic state, which is a
challenge for genetic diagnosis.


  Urothelial cell carcinoma (UCC) is the ninth most common cancer that accounts
for 4.7% of all the new cancer cases globally. UCC development and progression
are due to complex and stochastic genetic programmes. To study the cascades of
molecular events underlying the poor prognosis that may lead to limited
treatment options for advanced disease and resistance to conventional therapies
in UCC, transcriptomics technology (RNA-Seq), a method of analysing the RNA
content of a sample using modern high-throughput sequencing platforms has been
employed. Here we review the principles of RNA-Seq technology and summarize
recent studies on human bladder cancer that employed this technique to unravel
the pathogenesis of the disease, identify biomarkers, discover pathways and
classify the disease state. We list the commonly used computational platforms
and software that are publicly available for RNA-Seq analysis. Moreover, we
discussed the future perspectives for RNA-Seq studies on bladder cancer and
recommend the application of new technology called single cell sequencing
(scRNA-Seq) to further understand the disease. Keywords: Transcriptome
profiling, RNA-sequencing, genomics, bioinformatics, bladder cancer


  Background: Progressive supranuclear palsy (PSP) is a rare neurodegenerative
disease causing parkinsonian symptoms. Altered DNA methylation of the
microtubule-associated protein tau gene correlates with the expression changes
in Alzheimer's disease and Parkinson's disease brains. However, few studies
examine the sequences beyond the constitutive promoter.Objectives: Because
activating different microtubule associated protein tau gene control regions
via methylation might regulate the differential tau expression constituting the
specific signatures of individual tauopathies, we compared methylation of a
candidate promoter, intron 0.Methods: We assessed DNA methylation in the brains
of patients with different tauopathies (35 Alzheimer's disease, 10 corticobasal
degeneration, and 18 PSP) and 19 controls by intron 0 pyrosequencing. We also
evaluated methylation in an independent cohort of 11 PSP cases and 12 controls.
Frontal (affected by tau pathology) and occipital (unaffected) cortices were
analyzed.Results: In the initial samples, one CpG island site in intron 0
(CpG1) showed significant hypomethylation in PSP-affected frontal cortices when
compared with controls (p = 0.022). Such hypomethylation was observed in
replicate samples, but not in occipital cortices or other tauopathies. PSP and
control samples (combining the initial and replicate samples) remained
significantly different after adjustment for potential confounding factors
(age, H1/H1 diplotype; p = 0.0005). PSP-affected tissues exhibited
microtubule-associated protein tau RNA hyperexpression when compared with
controls (p = 0.004), although no correlation with CpG1 methylation was
observed.Conclusions: This exploratory study suggests that regions other than
the constitutive promoter may be involved in microtubule-associated protein tau
gene regulation in tauopathies and that intron 0 hypomethylation may be a
specific epigenetic signature of PSP. These preliminary findings require
confirmation.


  BACKGROUND: Spinocerebellar ataxia types 19 and 22 (SCA19/22) are rare
conditions in which relatively isolated cerebellar involvement is frequently
associated with cognitive impairment. Here, we report on new clinical features
and provide details of the cognitive profile in two SCA19/22 families.METHODS:
Two families displaying an autosomal-dominant form of cerebellar ataxia
underwent clinical examinations and genetic testing.RESULTS: In addition to the
classical clinical features of SCA, a wide spectrum of cognitive disorders
(including visuospatial impairments) was observed. Eight patients had mild
Parkinsonism, and five had epilepsy. Genetic testing showed that the KCND3
mutation (c.679_681delTTC, p.F227del) was present in both families.CONCLUSIONS:
Our findings broaden the phenotypic spectrum of SCA19/22, and suggest that
KCND3 should be included in the list of candidate genes for epilepsy,
Parkinsonism and cognitive impairment.


  This paper describes a mathematical model for the spread of a virus through
an isolated population of a given size. The model uses three, color-coded
components, called molecules (red for infected and still contagious; green for
infected, but no longer contagious; and blue for uninfected). In retrospect,
the model turns out to be a digital analogue for the well-known SIR model of
Kermac and McKendrick (1927). In our RGB model, the number of accumulated
infections goes through three phases, beginning at a very low level, then
changing to a transition ramp of rapid growth, and ending in a plateau of final
values. Consequently, the differential change or growth rate begins at 0, rises
to a peak corresponding to the maximum slope of the transition ramp, and then
falls back to 0. The properties of these time variations, including the slope,
duration, and height of the transition ramp, and the width and height of the
infection rate, depend on a single parameter - the time that a red molecule is
contagious divided by the average time between collisions of the molecules.
Various temporal milestones, including the starting time of the transition
ramp, the time that the accumulating number of infections obtains its maximum
slope, and the location of the peak of the infection rate depend on the size of
the population in addition to the contagious lifetime ratio. Explicit formulas
for these quantities are derived and summarized. Finally, Appendix E has been
added to describe the effect of vaccinations.


  In this paper, we propose a simplified bidimensional Wolbachia infestation
model in a population of Aedes aegypti mosquitoes, preserving the main features
associated with the biology of this species that can be found in
higher-dimensional models. Namely, our model represents the maternal
transmission of the Wolbachia symbiont, expresses the reproductive phenotype of
cytoplasmic incompatibility, accounts for different fecundities and mortalities
of infected and wild insects, and exhibits the bistable nature leading to the
so-called principle of competitive exclusion. Since Wolbachia-based biocontrol
is now accepted as an ecologically friendly and potentially cost-effective
method for prevention and control of dengue and other arboviral infections, it
is essential to have reduced models with the main biological characteristics of
Aedes aegypti in the presence of Wolbachia-carriers because such models help to
simplify the mathematical analysis for determining appropriate biocontrol
strategies. Using tools borrowed from monotone dynamical system theory, in the
proposed model, we prove the existence of an invariant threshold manifold that
allows us to provide practical recommendations for performing single and
periodic releases of Wolbachia-carrying mosquitoes, seeking the eventual
elimination of wild insects that are capable of transmitting infections to
humans. We illustrate these findings with numerical simulations using parameter
values corresponding to the wMelPop strain of Wolbachia that is considered the
best virus blocker but induces fitness loss in its carriers.


  This paper develops an easily-implementable version of Page's CUSUM
quickest-detection test, designed to work in certain composite hypothesis
scenarios with time-varying data statistics. The decision statistic can be cast
in a recursive form and is particularly suited for on-line analysis. By
back-testing our approach on publicly-available COVID-19 data we find reliable
early warning of infection flare-ups, in fact sufficiently early that the tool
may be of use to decision-makers on the timing of restrictive measures that may
in the future need to be taken.


  Due to structural and functional abnormalities or genetic variations and
mutations, there may be dysfunctional molecules within an intracellular
signaling network that do not allow the network to correctly regulate its
output molecules, such as transcription factors. This disruption in signaling
interrupts normal cellular functions and may eventually develop some
pathological conditions. In this paper, computation capacity of signaling
networks is introduced as a fundamental limit on signaling capability and
performance of such networks. The computation capacity measures the maximum
number of computable inputs, that is, the maximum number of input values for
which the correct functional output values can be recovered from the erroneous
network outputs, when the network contains some dysfunctional molecules. This
contrasts with the conventional communication capacity that measures instead
the maximum number of input values that can be correctly distinguished based on
the erroneous network outputs.
  The computation capacity is higher than the communication capacity, if the
network response function is not a one-to-one function of the input signals. By
explicitly incorporating the effect of signaling errors that result in the
network dysfunction, the computation capacity provides more information about
the network and its malfunction. Two examples of signaling networks are studied
here, one regulating caspase3 and another regulating NFkB, for which
computation and communication capacities are analyzed. Higher computation
capacities are observed for both networks. One biological implication of this
finding is that signaling networks may have more capacity than that specified
by the conventional communication capacity metric. The effect of feedback is
also studied. In summary, this paper reports findings on a new fundamental
feature of the signaling capability of cell signaling networks.


  Transients are fundamental to ecological systems with significant
implications to management, conservation, and biological control. We uncover a
type of transient synchronization behavior in spatial ecological networks whose
local dynamics are of the chaotic, predator-prey type. In the parameter regime
where there is phase synchronization among all the patches, complete
synchronization (i.e., synchronization in both phase and amplitude) can arise
in certain pairs of patches as determined by the network symmetry - henceforth
the phenomenon of "synchronization within synchronization." Distinct patterns
of complete synchronization coexist but, due to intrinsic instability or noise,
each pattern is a transient and there is random, intermittent switching among
the patterns in the course of time evolution. The probability distribution of
the transient time is found to follow an algebraic scaling law with a divergent
average transient lifetime. Based on symmetry considerations, we develop a
stability analysis to understand these phenomena. The general principle of
symmetry can also be exploited to explain previously discovered,
counterintuitive synchronization behaviors in ecological networks.


  How can we learn a dynamical system to make forecasts, when some variables
are unobserved? For instance, in COVID-19, we want to forecast the number of
infected and death cases but we do not know the count of susceptible and
exposed people. While mechanics compartment models are widely used in epidemic
modeling, data-driven models are emerging for disease forecasting. We first
formalize the learning of physics-based models as AutoODE, which leverages
automatic differentiation to estimate the model parameters. Through a benchmark
study on COVID-19 forecasting, we notice that physics-based mechanistic models
significantly outperform deep learning. Our method obtains a 57.4% reduction in
mean absolute errors for 7-day ahead COVID-19 forecasting compared with the
best deep learning competitor. Such performance differences highlight the
generalization problem in dynamical system learning due to distribution shift.
We identify two scenarios where distribution shift can occur: changes in data
domain and changes in parameter domain (system dynamics). Through systematic
experiments on several dynamical systems, we found that deep learning models
fail to forecast well under both scenarios. While much research on distribution
shift has focused on changes in the data domain, our work calls attention to
rethink generalization for learning dynamical systems.


  We propose a $\tau$-leaping simulation algorithm for stochastic systems
subject to fast environmental changes. Similar to conventional $\tau$-leaping
the algorithm proceeds in discrete time steps, but as a principal addition it
captures environmental noise beyond the adiabatic limit. The key idea is to
treat the input rates for the $\tau$-leaping as (clipped) Gaussian random
variables with first and second moments constructed from the environmental
process. In this way, each step of the algorithm retains environmental
stochasticity to sub-leading order in the time scale separation between system
and environment. We test the algorithm on several toy examples with discrete
and continuous environmental states, and find good performance in the regime of
fast environmental dynamics. At the same time, the algorithm requires
significantly less computing time than full simulations of the combined system
and environment. In this context we also discuss several methods for the
simulation of stochastic population dynamics in time-varying environments with
continuous states.


  A method is proposed to select the suitable sets of potential parameters for
a one-dimensional mesoscopic Hamiltonian model, first introduced to describe
the DNA melting transition and later extended to investigate thermodynamic and
dynamical properties of nucleic acids. The DNA base pair fluctuations are
considered as time dependent trajectories whose initial condition sets the no
crossing constraint enforced in the path integral for the first-passage
probability. Performing the path integration at room temperature, relations are
established between the cutoff on the amplitude of the base pair fluctuations
and the model parameters. In particular, it is shown that the non-linear
stacking parameter should be $\sim 1$. The formalism here developed may be
applied to compute the lifetime of open base pairs in three-dimensional helical
models for DNA molecules.


  Autosomal dominant cerebellar ataxia corresponds to a clinically and
genetically heterogeneous group of neurodegenerative disorders that primarily
affect the cerebellum. Here, we report the identification of the causative gene
in spinocerebellar ataxia 21, an autosomal-dominant disorder previously mapped
to chromosome 7p21.3-p15.1. This ataxia was firstly characterized in a large
French family with slowly progressive cerebellar ataxia, accompanied by severe
cognitive impairment and mental retardation in two young children. Following
the recruitment of 12 additional young family members, linkage analysis enabled
us to definitively map the disease locus to chromosome 1p36.33-p36.32. The
causative mutation, (c.509C4T/p.P170L) in the transmembrane protein gene
TMEM240, was identified by whole exome sequencing and then was confirmed by
Sanger sequencing and co-segregation analyses. Index cases from 368 French
families with autosomal-dominant cerebellar ataxia were also screened for
mutations. In seven cases, we identified a range of missense mutations
(c.509C4T/p.P170L, c.239C4T/p.T80M, c.346C4T/p.R116C, c.445G4A/p.E149K,
c.511C4T/p.R171W), and a stop mutation (c.489C4G/p.Y163*) in the same gene.
TMEM240 is a small, strongly conserved transmembrane protein of unknown
function present in cerebellum and brain. Spinocerebellar ataxia 21 may be a
particular early-onset disease associated with severe cognitive impairment.


  Mathematical and computational models can assist in gaining an understanding
of cell behavior at many levels of organization. Here, we review models in the
literature that focus on eukaryotic cell motility at 3 size scales:
intracellular signaling that regulates cell shape and movement, single cell
motility, and collective cell behavior from a few cells to tissues. We survey
recent literature to summarize distinct computational methods (phase-field,
polygonal, Cellular Potts, and spherical cells). We discuss models that bridge
between levels of organization, and describe levels of detail, both biochemical
and geometric, included in the models. We also highlight links between models
and experiments. We find that models that span the 3 levels are still in the
minority.


  Whirling and swerving, a bacterium is swimming in a test tube, foraging for
food. On the surface of a vibrating bath, a droplet starts walking. A certain
similarity, but mostly dissimilarity, between the physical memory that emerges
in Couder's droplet experiments and the biological memory of the bacterium is
noted. It serves as a starting point for a short perspective and speculation on
the multilevel, loopy memory of living matter.


  Stratifying cancer patients based on their gene expression levels allows
improving diagnosis, survival analysis and treatment planning. However, such
data is extremely highly dimensional as it contains expression values for over
20000 genes per patient, and the number of samples in the datasets is low. To
deal with such settings, we propose to incorporate prior biological knowledge
about genes from ontologies into the machine learning system for the task of
patient classification given their gene expression data. We use ontology
embeddings that capture the semantic similarities between the genes to direct a
Graph Convolutional Network, and therefore sparsify the network connections. We
show this approach provides an advantage for predicting clinical targets from
high-dimensional low-sample data.


  Single-particle cryo-electron microscopy (cryo-EM) is an emerging imaging
modality capable of visualizing proteins and macro-molecular complexes at
near-atomic resolution. The low electron-doses used to prevent sample radiation
damage, result in images where the power of the noise is 100 times greater than
the power of the signal. To overcome the low-SNRs, hundreds of thousands of
particle projections acquired over several days of data collection are averaged
in 3D to determine the structure of interest. Meanwhile, recent image
super-resolution (SR) techniques based on neural networks have shown state of
the art performance on natural images. Building on these advances, we present a
multiple-image SR algorithm based on deep internal learning designed
specifically to work under low-SNR conditions. Our approach leverages the
internal image statistics of cryo-EM movies and does not require training on
ground-truth data. When applied to a single-particle dataset of apoferritin, we
show that the resolution of 3D structures obtained from SR micrographs can
surpass the limits imposed by the imaging system. Our results indicate that the
combination of low magnification imaging with image SR has the potential to
accelerate cryo-EM data collection without sacrificing resolution.


  Identifying subgroups and properties of cancer biopsy samples is a crucial
step towards obtaining precise diagnoses and being able to perform personalized
treatment of cancer patients. Recent data collections provide a comprehensive
characterization of cancer cell data, including genetic data on copy number
alterations (CNAs). We explore the potential to capture information contained
in cancer genomic information using a novel topology-based approach that
encodes each cancer sample as a persistence diagram of topological features,
i.e., high-dimensional voids represented in the data. We find that this
technique has the potential to extract meaningful low-dimensional
representations in cancer somatic genetic data and demonstrate the viability of
some applications on finding substructures in cancer data as well as comparing
similarity of cancer types.


  An analysis of the COVID-19 epidemic is proposed on the basis of the epiMOX
dashboard (publicly accessible at https://www.epimox.polimi.it) that deals with
data of the epidemic trends and outbreaks in Italy from late February 2020. Our
analysis provides an immediate appreciation of the past epidemic development,
together with its current trends by fostering a deeper interpretation of
available data through several critical epidemic indicators. In addition, we
complement the epiMOX dashboard with a predictive tool based on an
epidemiological compartmental model, named SUIHTER, for the forecast on the
near future epidemic evolution.


  Boltzmann machines (BM) are widely used as generative models. For example,
pairwise Potts models (PM), which are instances of the BM class, provide
accurate statistical models of families of evolutionarily related protein
sequences. Their parameters are the local fields, which describe site-specific
patterns of amino-acid conservation, and the two-site couplings, which mirror
the coevolution between pairs of sites. This coevolution reflects structural
and functional constraints acting on protein sequences during evolution. The
most conservative choice to describe the coevolution signal is to include all
possible two-site couplings into the PM. This choice, typical of what is known
as Direct Coupling Analysis, has been successful for predicting residue
contacts in the three-dimensional structure, mutational effects, and in
generating new functional sequences. However, the resulting PM suffers from
important over-fitting effects: many couplings are small, noisy and hardly
interpretable; the PM is close to a critical point, meaning that it is highly
sensitive to small parameter perturbations. In this work, we introduce a
general parameter-reduction procedure for BMs, via a controlled iterative
decimation of the less statistically significant couplings, identified by an
information-based criterion that selects either weak or statistically
unsupported couplings. For several protein families, our procedure allows one
to remove more than $90\%$ of the PM couplings, while preserving the predictive
and generative properties of the original dense PM, and the resulting model is
far away from criticality, hence more robust to noise.


  High-throughput sequencing can detect tens of thousands of genes in parallel,
providing opportunities for improving the diagnostic accuracy of multiple
diseases including sepsis, which is an aggressive inflammatory response to
infection that can cause organ failure and death. Early screening of sepsis is
essential in clinic, but no effective diagnostic biomarkers are available yet.
Here, we present a novel method, Recurrent Logistic Regression, to identify
diagnostic biomarkers for sepsis from the blood transcriptome data. A panel
including five immune-related genes, LRRN3, IL2RB, FCER1A, TLR5, and S100A12,
are determined as diagnostic biomarkers (LIFTS) for sepsis. LIFTS discriminates
patients with sepsis from normal controls in high accuracy (AUROC = 0.9959 on
average; IC = [0.9722-1.0]) on nine validation cohorts across three independent
platforms, which outperforms existing markers. Our analysis determined an
accurate prediction model and reproducible transcriptome biomarkers that can
lay a foundation for clinical diagnostic tests and biological mechanistic
studies.


  The traditional long-term solutions for epidemic control involve eradication
or population immunity. Here, we analytically derive the existence of a third
viable solution: a stable equilibrium at low case numbers, where
test-trace-and-isolate policies partially compensate for local spreading
events, and only moderate restrictions remain necessary. In this equilibrium,
daily cases stabilize around ten new infections per million people or less.
However, stability is endangered if restrictions are relaxed or case numbers
grow too high. The latter destabilization marks a tipping point beyond which
the spread self-accelerates. We show that a lockdown can reestablish control
and that recurring lockdowns are not necessary given sustained, moderate
contact reduction. We illustrate how this strategy profits from vaccination and
helps mitigate variants of concern. This strategy reduces cumulative cases (and
fatalities) 4x more than strategies that only avoid hospital collapse. In the
long term, immunization, large-scale testing, and international coordination
will further facilitate control.


  Containment of epidemic outbreaks entails great societal and economic costs.
Cost-effective containment strategies rely on efficiently identifying infected
individuals, making the best possible use of the available testing resources.
Therefore, quickly identifying the optimal testing strategy is of critical
importance. Here, we demonstrate that machine learning can be used to identify
which individuals are most beneficial to test, automatically and dynamically
adapting the testing strategy to the characteristics of the disease outbreak.
Specifically, we simulate an outbreak using the archetypal
susceptible-infectious-recovered (SIR) model and we use data about the first
confirmed cases to train a neural network that learns to make predictions about
the rest of the population. Using these prediction, we manage to contain the
outbreak more effectively and more quickly than with standard approaches.
Furthermore, we demonstrate how this method can be used also when there is a
possibility of reinfection (SIRS model) to efficiently eradicate an endemic
disease.


  Consider a mathematical model of evolutionary adaptation of fitness landscape
and mutation matrix as a reaction to population changes. As a basis, we use an
open quasispecies model, which is modified to include explicit death flow. We
assume that evolutionary parameters of mutation and selection processes vary in
a way to maximize the mean fitness of the system. From this standpoint,
Fisher's theorem of natural selection is being rethought and discussed. Another
assumption is that system dynamics has two significant timescales. According to
our central hypothesis, major evolutionary transitions happen in the
steady-state of the corresponding dynamical system, so the evolutionary time is
much slower than the one of internal dynamics. For the specific cases of
quasispecies systems, we show how our premises form the fitness landscape
adaptation process.


  People tend to align their use of language to the linguistic behaviour of
their own ingroup and to simultaneously diverge from the language use of
outgroups. This paper proposes to model this phenomenon of sociolinguistic
identity maintenance as an evolutionary game in which individuals play the
field and the dynamics are supplied by a multi-population extension of the
replicator-mutator equation. Using linearization, the stabilities of all
dynamic equilibria of the game in its fully symmetric two-population special
case are found. The model is then applied to an empirical test case from
adolescent sociolinguistic behaviour. It is found that the empirically attested
population state corresponds to one of a number of stable equilibria of the
game under an independently plausible value of a parameter controlling the rate
of linguistic mutations. An asymmetric three-population extension of the game,
explored with numerical solution methods, furthermore predicts to which
specific equilibrium the system converges.


  Many systems may switch to an undesired state due to internal failures or
external perturbations, of which critical transitions toward degraded ecosystem
states are a prominent example. Resilience restoration focuses on the ability
of spatially-extended systems and the required time to recover to their desired
states under stochastic environmental conditions. While mean-field approaches
may guide recovery strategies by indicating the conditions needed to
destabilize undesired states, these approaches are not accurately capturing the
transition process toward the desired state of spatially-extended systems in
stochastic environments. The difficulty is rooted in the lack of mathematical
tools to analyze systems with high dimensionality, nonlinearity, and stochastic
effects. We bridge this gap by developing new mathematical tools that employ
nucleation theory in spatially-embedded systems to advance resilience
restoration. We examine our approach on systems following mutualistic dynamics
and diffusion models, finding that systems may exhibit single-cluster or
multi-cluster phases depending on their sizes and noise strengths, and also
construct a new scaling law governing the restoration time for arbitrary system
size and noise strength in two-dimensional systems. This approach is not
limited to ecosystems and has applications in various dynamical systems, from
biology to infrastructural systems.


  Double-stranded DNA (dsDNA) has been established as an efficient medium for
charge migration, bringing it to the forefront of the field of molecular
electronics as well as biological research. The charge migration rate is
controlled by the electronic couplings between the two nucleobases of DNA/RNA.
These electronic couplings strongly depend on the intermolecular geometry and
orientation. Estimating these electronic couplings for all the possible
relative geometries of molecules using the computationally demanding
first-principles calculations requires a lot of time as well as computation
resources. In this article, we present a Machine Learning (ML) based model to
calculate the electronic coupling between any two bases of dsDNA/dsRNA of any
length and sequence and bypass the computationally expensive first-principles
calculations. Using the Coulomb matrix representation which encodes the atomic
identities and coordinates of the DNA base pairs to prepare the input dataset,
we train a feedforward neural network model. Our NN model can predict the
electronic couplings between dsDNA base pairs with any structural orientation
with a MAE of less than 0.014 eV. We further use the NN predicted electronic
coupling values to compute the dsDNA/dsRNA conductance.


  The idea that structural disorder might be a novel mechanism of protein
interaction is widespread in the Literature, although the number of
statistically significant structural studies supporting this is surprisingly
low. At variance with previous works, our conclusions rely exclusively on a
large-scale analysis of all the 134337 X-ray crystallographic structures of the
Protein Data Bank averaged over clusters of almost identical protein sequences.
In this work, we explore the complexity of the organization of all the
interaction interfaces observed when a protein lies in alternative complexes,
showing that interfaces progressively add up in a hierarchical way. We further
investigate the connection of this complexity with different measures of
structural disorder: the standard missing residues and a new definition, called
"soft disorder", that covers all the flexible and structurally amorphous
residues of a protein. We show evidences that both the interaction interfaces
and the soft disordered regions tend to involve roughly the same amino-acids of
the protein, and preliminary results suggesting that soft disorder spots those
surface regions where new interfaces are progressively accommodated by complex
formation. Our results suggest that disordered regions not only carry crucial
information about the location of alternative interfaces within complexes, but
also of the order of the assembly. We verify these hypotheses in several
examples. We finally compare our measures of disorder with several disorder
predictors, showing that these latter are optimized to predict the residues
that are missing in all the alternative structures of a protein, and they are
not able to catch the progressive evolution of the disordered regions upon
complex formation. Yet, the predicted residues, if not missing, tend to be
characterized as soft disordered.


  The problem of estimating the structure of a graph from observed data is of
growing interest in the context of high-throughput genomic data, and
single-cell RNA sequencing in particular. These, however, are challenging
applications, since the data consist of high-dimensional counts with high
variance and over-abundance of zeros. Here, we present a general framework for
learning the structure of a graph from single-cell RNA-seq data, based on the
zero-inflated negative binomial distribution. We demonstrate with simulations
that our approach is able to retrieve the structure of a graph in a variety of
settings and we show the utility of the approach on real data.


  As of November 2020, the number of COVID-19 cases is increasing rapidly in
many countries. In Europe, the virus spread slowed considerably in the late
spring due to strict lockdown, but a second wave of the pandemic grew
throughout the fall. In this study, we first reconstruct the time evolution of
the effective reproduction numbers ${\cal R}(t)$ for each country by
integrating the equations of the classic SIR model. We cluster countries based
on the estimated ${\cal R}(t)$ through a suitable time series dissimilarity.
The result suggests that simple dynamical mechanisms determine how countries
respond to changes in COVID-19 case counts. Inspired by these results, we
extend the SIR model to include a social response to explain the number $X(t)$
of new confirmed daily cases. As a first-order model, we assume that the social
response is on the form $d_t{\cal R}=-\nu (X-X^*)$, where $X^*$ is a threshold
for response. The response rate $\nu$ depends on whether $X^*$ is below or
above this threshold, on three parameters $\nu_1,\;\nu_2,\,\nu_3,$, and on $t$.
When $X<X^*$, $\nu= \nu_1$, describes the effect of relaxed intervention when
the incidence rate is low. When $X>X^*$, $\nu=\nu_2\exp{(-\nu_3t)}$, models the
impact of interventions when incidence rate is high. The parameter $\nu_3$
represents the fatigue, i.e., the reduced effect of intervention as time
passes. The proposed model reproduces typical evolving patterns of COVID-19
epidemic waves observed in many countries. Estimating the parameters
$\nu_1,\,\nu_2,\,\nu_3$ and initial conditions, such as ${\cal R}_0$, for
different countries helps to identify important dynamics in their social
responses. One conclusion is that the leading cause of the strong second wave
in Europe in the fall of 2020 was not the relaxation of interventions during
the summer, but rather the general fatigue to interventions developing in the
fall.


  This technical report describes the rationale and technical details for the
dynamic causal modelling of mitigated epidemiological outcomes based upon a
variety of timeseries data. It details the structure of the underlying
convolution or generative model (at the time of writing on 6-Nov-20). This
report is intended for use as a reference that accompanies the predictions in
following dashboard: https://www.fil.ion.ucl.ac.uk/spm/covid-19/dashboard


  SARS-CoV-2 is a new human coronavirus (CoV), which emerged in China in late
2019 and is responsible for the global COVID-19 pandemic that caused more than
59 million infections and 1.4 million deaths in 11 months. Understanding the
origin of this virus is an important issue and it is necessary to determine the
mechanisms of its dissemination in order to contain future epidemics. Based on
phylogenetic inferences, sequence analysis and structure-function relationships
of coronavirus proteins, informed by the knowledge currently available on the
virus, we discuss the different scenarios evoked to account for the origin -
natural or synthetic - of the virus. The data currently available is not
sufficient to firmly assert whether SARS-CoV2 results from a zoonotic emergence
or from an accidental escape of a laboratory strain. This question needs to be
solved because it has important consequences on the evaluation of risk/benefit
balance of our interaction with ecosystems, the intensive breeding of wild and
domestic animals, as well as some lab practices and on scientific policy and
biosafety regulations. Regardless of its origin, studying the evolution of the
molecular mechanisms involved in the emergence of pandemic viruses is essential
to develop therapeutic and vaccine strategies and to prevent future zoonoses.
This article is a translation and update of a French article published in
M{\'e}decine/Sciences, Aug/Sept 2020 (http://doi.org/10.1051/medsci/2020123).


  Genome assembly is a fundamental problem in Bioinformatics, requiring to
reconstruct a source genome from an assembly graph built from a set of reads
(short strings sequenced from the genome). A notion of genome assembly solution
is that of an arc-covering walk of the graph. Since assembly graphs admit many
solutions, the goal is to find what is definitely present in all solutions, or
what is safe. Most practical assemblers are based on heuristics having at their
core unitigs, namely paths whose internal nodes have unit in-degree and
out-degree, and which are clearly safe. The long-standing open problem of
finding all the safe parts of the solutions was recently solved [RECOMB 2016]
yielding a 60% increase in contig length. This safe and complete genome
assembly algorithm was followed by other works improving the time bounds, as
well as extending the results for different notions of assembly solution. But
it remained open whether one can be complete also for models of genome assembly
of practical applicability.
  In this paper we present a universal framework for obtaining safe and
complete algorithms which unify the previous results, while also allowing for
easy generalisations to assembly problems including many practical aspects.
This is based on a novel graph structure, called the hydrostructure of a walk,
which highlights the reachability properties of the graph from the perspective
of the walk. The hydrostructure allows for simple characterisations of the
existing safe walks, and of their new practical versions. Almost all of our
characterisations are directly adaptable to optimal verification algorithms,
and simple enumeration algorithms. Most of these algorithms are also improved
to optimality using an incremental computation procedure and a previous optimal
algorithm of a specific model.


  Pandemics, like the 1918 Spanish Influenza and COVID-19, spread through
regions of the World in subsequent waves. There is, however, no consensus on
the origin of this pattern, which may originate from human behaviour rather
than from the virus diffusion itself. Time-honoured models of the SIR type or
others based on complex networks describe well the exponential spread of the
disease, but cannot naturally accommodate the wave pattern. Nevertheless,
understanding this time-structure is of paramount importance in designing
effective prevention measures. Here we propose a consistent picture of the wave
pattern based on the epidemic Renormalisation Group (eRG) framework, which is
guided by the global symmetries of the system under time rescaling. We show
that the rate of spreading of the disease can be interpreted as a time-dilation
symmetry, while the final stage of an epidemic episode corresponds to reaching
a time scale-invariant state. We find that the endemic period between two waves
is a sign of instability in the system, associated to near-breaking of the time
scale-invariance. This phenomenon can be described in terms of an eRG model
featuring complex fixed points. Our results demonstrate that the key to control
the arrival of the next wave of a pandemic is in the strolling period in
between waves, i.e. when the number of infections grow linearly. Thus, limiting
the virus diffusion in this period is the most effective way to prevent or
delay the arrival of the next wave. In this work we establish a new guiding
principle for the formulation of mid-term governmental strategies to curb
pandemics and avoid recurrent waves of infections, deleterious in terms of
human life loss and economic damage.


  Genome-wide association studies (GWAS) require accurate cohort phenotyping,
but expert labeling can be costly, time-intensive, and variable. Here we
develop a machine learning (ML) model to predict glaucomatous optic nerve head
features from color fundus photographs. We used the model to predict vertical
cup-to-disc ratio (VCDR), a diagnostic parameter and cardinal endophenotype for
glaucoma, in 65,680 Europeans in the UK Biobank (UKB). A GWAS of ML-based VCDR
identified 299 independent genome-wide significant (GWS; $P\leq5\times10^{-8}$)
hits in 156 loci. The ML-based GWAS replicated 62 of 65 GWS loci from a recent
VCDR GWAS in the UKB for which two ophthalmologists manually labeled images for
67,040 Europeans. The ML-based GWAS also identified 92 novel loci,
significantly expanding our understanding of the genetic etiologies of glaucoma
and VCDR. Pathway analyses support the biological significance of the novel
hits to VCDR, with select loci near genes involved in neuronal and synaptic
biology or known to cause severe Mendelian ophthalmic disease. Finally, the
ML-based GWAS results significantly improve polygenic prediction of VCDR and
primary open-angle glaucoma in the independent EPIC-Norfolk cohort.


  Progress in real-time, simultaneous in vivo detection of multiple
neurotransmitters will help accelerate advances in neuroscience research. The
need for development of probes capable of stable electrochemical detection of
rapid neurotransmitter fluctuations with high sensitivity and selectivity and
sub-second temporal resolution has, therefore, become compelling. Additionally,
a higher spatial resolution multi-channel capability is required to capture the
complex neurotransmission dynamics across different brain regions. These
research needs have inspired the introduction of glassy carbon (GC)
microelectrode arrays on flexible polymer substrates through carbon MEMS
(C-MEMS) microfabrication process followed by a novel pattern transfer
technique. These implantable GC microelectrodes offer unique advantages in
electrochemical detection of electroactive neurotransmitters through the
presence of active carboxyl, carbonyl, and hydroxyl functional groups. In
addition, they offer fast electron transfer kinetics, capacitive
electrochemical behavior, and wide electrochemical window. Here, we combine the
use of these GC microelectrodes with the fast scan cyclic voltammetry (FSCV)
technique to optimize the co-detection of dopamine and serotonin in vitro and
in vivo. We demonstrate that using optimized FSCV triangular waveform at scan
rates lower than 700 V/s and holding and switching at potentials of 0.4 and 1V
respectively, it is possible to discriminate voltage reduction and oxidation
peaks of serotonin and dopamine, with serotonin contributing distinct multiple
oxidation peaks. Taken together, our results present a compelling case for a
carbon-based MEA platform rich with active functional groups that allows for
repeatable and stable detection of electroactive multiple neurotransmitters at
concentrations as low as 10 nM


  Proteins constitute a large group of macromolecules with a multitude of
functions for all living organisms. Proteins achieve this by adopting distinct
three-dimensional structures encoded by the sequence of their constituent amino
acids in one or more polypeptides. In this paper, the statistical modelling of
the protein backbone torsion angles is considered. Two new distributions are
proposed for toroidal data by applying the M\"obius transformation to the
bivariate von Mises distribution. Marginal and conditional distributions in
addition to sine-skewed versions of the proposed models are also developed.
Three big data sets consisting of bivariate information about protein domains
are analysed to illustrate the strength of the flexible proposed models.
Finally, a simulation study is done to evaluate the obtained maximum likelihood
estimates and also to find the best method of generating samples from the
proposed models to use as the proposal distributions in the Markov Chain Monte
Carlo sampling method for predicting the 3D structure of proteins.


  The vulnerability to trawling on a species level is determined by a species'
individual combination of biological traits that is related to the ecosystem
functions. Trait-based indices of physical resistance RI and reproductive
potential RPI were developed and combined into an overall vulnerability index
on a species level, the RRI or Resistance and Reproductive Potential Index. The
indices are used to explore how resistance and reproductive potentia change
over a trawling gradient. The RRI allows for dividing the benthic community
into groups expressing different levels of vulnerability that can be linked to
ecosystem functions. The RRI index opens up the possibility of scenario
modelling by simulating the extinction or loss of vulnerable species and its
effects on functions. The validity of the trait-based RRI index was explored by
comparing individual species' RRI scores to empirically observed responses over
a trawling gradient based on a previously published data. RRI score and
observed responses (regression slopes) were significantly correlated providing
support for the rationality of the approach taken. Data analyses evidenced
increases of resistance and resilience indices over the trawling gradient,
demonstrating that communities lost vulnerable species with increasing
trawling. When exploring the effects of trawling on the bioturbation, we found
it to be disproportionately affected though the loss of vulnerable species. The
proposed indices provide new insights into the link of species vulnerability
and function.


  In recent years, cell-based assays have been frequently used in molecular
interaction analysis. Cell-based assays complement traditional biochemical and
biophysical methods, as they allow for molecular interaction analysis, mode of
action studies and even drug screening processes to be performed under
physiologically relevant conditions. In most cellular assays, biomolecules are
usually labeled to achieve specificity. In order to overcome some of the
drawbacks associated with label-based assays, we have recently introduced
cell-based molography as a biosensor for the analysis of specific molecular
interactions involving native membrane receptors in living cells. Here, we
expand this assay to cytosolic protein-protein interactions. First, we created
a biomimetic membrane receptor by tethering one cytosolic interaction partner
to the plasma membrane. The artificial construct is then coherently arranged
into a two-dimensional pattern within the cytosol of living cells. Thanks to
the molographic sensor, the specific interactions between the coherently
arranged protein and its endogenous interaction partners become visible in
real-time without the use of a fluorescent label. This method turns out to be
an important extension of cell-based molography because it expands the range of
interactions that can be analyzed by molography to those in the cytosol of
living cells.


  Recently, the perfect phylogeny model with persistent characters has
attracted great attention in the literature. It is based on the assumption that
complex traits or characters can only be gained once and lost once in the
course of evolution. Here, we consider a generalization of this model, namely
Dollo parsimony, that allows for multiple character losses. More precisely, we
take a combinatorial perspective on the notion of Dollo-$k$ characters, i.e.
traits that are gained at most once and lost precisely $k$ times throughout
evolution. We first introduce an algorithm based on the notion of spanning
subtrees for finding a Dollo-$k$ labeling for a given character and a given
tree in linear time. We then compare persistent characters (consisting of the
union of Dollo-0 and Dollo-1 characters) and general Dollo-$k$ characters.
While it is known that there is a strong connection between Fitch parsimony and
persistent characters, we show that Dollo parsimony and Fitch parsimony are in
general very different. Moreover, while it is known that there is a direct
relationship between the number of persistent characters and the Sackin index
of a tree, a popular index of tree balance, we show that this relationship does
not generalize to Dollo-$k$ characters. In fact, determining the number of
Dollo-$k$ characters for a given tree is much more involved than counting
persistent characters, and we end this manuscript by introducing a recursive
approach for the former. This approach leads to a polynomial time algorithm for
counting the number of Dollo-$k$ characters, and both this algorithm as well as
the algorithm for computing Dollo-$k$ labelings are publicly available in the
Babel package for BEAST 2.


  The importance of a strict quarantine has been widely debated during the
COVID-19 epidemic even from the purely epidemiological point of view. One
argument against strict lockdown measures is that once the strict quarantine is
lifted, the epidemic comes back, and so the cumulative number of infected
individuals during the entire epidemic will stay the same. We consider an SIR
model on a network and follow the disease dynamics, modeling the phases of
quarantine by changing the node degree distribution. We show that the system
reaches different steady states based on the history: the outcome of the
epidemic is path-dependent despite the same final node degree distribution. The
results indicate that two-phase route to the final node degree distribution (a
strict phase followed by a soft phase) are always better than one phase (the
same soft one) unless all the individuals have the same number of connections
at the end (the same degree); in the latter case, the overall number of
infected is indeed history-independent. The modeling also suggests that the
optimal procedure of lifting the quarantine consists of releasing nodes in the
order of their degree - highest first.


  A formal analogy of fluctuating diffusivity to thermodynamics is discussed
for messenger RNA molecules fluorescently fused to a protein in living cells.
Regarding the average value of the fluctuating diffusivity of such RNA-protein
particles as the analog of the internal energy, the analogs of the quantity of
heat and work are identified. The Clausius-like inequality is shown to hold for
the entropy associated with diffusivity fluctuations, which plays a role
analogous to the thermodynamic entropy, and the analog of the quantity of heat.
The change of the statistical fluctuation distribution is also examined from a
geometric perspective. The present discussions may contribute to a deeper
understanding of the fluctuating diffusivity in view of the laws of
thermodynamics.


  We present the phylogenetic quartet reconstruction method SAQ (Semi-algebraic
quartet reconstruction). SAQ is consistent with the most general Markov model
of nucleotide substitution and, in particular, it allows for rate heterogeneity
across lineages. Based on the algebraic and semi-algebraic description of
distributions that arise from the general Markov model on a quartet, the method
outputs normalized weights for the three trivalent quartets (which can be used
as input of quartet-base methods). We show that SAQ is a highly competitive
method that outperforms most of the well known reconstruction methods on data
simulated under the general Markov model on 4-taxon trees. Moreover, it also
achieves a high performance on data that violates the underlying assumptions.


  For fast development of COVID-19, it is only feasible to use drugs (off label
use) or approved natural products that are already registered or been assessed
for safety in previous human trials. These agents can be quickly assessed in
COVID-19 patients, as their safety and pharmacokinetics should already be well
understood. Computational methods offer promise for rapidly screening such
products for potential SARS-CoV-2 activity by predicting and ranking the
affinities of these compounds for specific virus protein targets. The
RNA-dependent RNA polymerase (RdRP) is a promising target for SARS-CoV-2 drug
development given it has no human homologs making RdRP inhibitors potentially
safer, with fewer off-target effects that drugs targeting other viral proteins.
We combined robust Vina docking on RdRP with molecular dynamic (MD) simulation
of the top 80 identified drug candidates to yield a list of the most promising
RdRP inhibitors. Literature reviews revealed that many of the predicted
inhibitors had been shown to have activity in in vitro assays or had been
predicted by other groups to have activity. The novel hits revealed by our
screen can now be conveniently tested for activity in RdRP inhibition assays
and if conformed testing for antiviral activity invitro before being tested in
human trials


  There are often multiple diseases with cross immunity competing for
vaccination resources. Here we investigate the optimal vaccination program in a
two-layer Susceptible-Infected-Removed (SIR) model, where two diseases with
cross immunity spread in the same population, and vaccines for both diseases
are available. We identify three scenarios of the optimal vaccination program,
which prevents the outbreaks of both diseases at the minimum cost. We
analytically derive a criterion to specify the optimal program based on the
costs for different vaccines.


  We study a spatial, one-shot prisoner's dilemma (PD) model in which selection
operates on both an organism's behavioral strategy (cooperate or defect) and
its choice of when to implement that strategy across a set of discrete time
slots. Cooperators evolve to fixation regularly in the model when we add time
slots to lattices and small-world networks, and their portion of the population
grows, albeit slowly, when organisms interact in a scale-free network. This
selection for cooperators occurs across a wide variety of time slots and it
does so even when a crucial condition for the evolution of cooperation on
graphs is violated--namely, when the ratio of benefits to costs in the PD does
not exceed the number of spatially-adjacent organisms.


  The bipolar organization of the microtubule-based mitotic spindle is
essential for the faithful segregation of chromosomes in cell division. Despite
our extensive knowledge of genes and proteins, the physical mechanism of how
the ensemble of microtubules can assemble into a proper bipolar shape remains
elusive. Here, we study the pathways of spindle self-organization using
cell-free Xenopus egg extracts and computer-based automated shape analysis. Our
microscopy assay allows us to simultaneously record the growth of hundreds of
spindles in the bulk cytoplasm and systematically analyze the shape of each
structure over the course of self-organization. We find that spindles that are
maturing into a bipolar shape take a route that is distinct from those ending
up with faulty structures, such as those of a tripolar shape. Moreover, matured
structures are highly stable with little occasions of transformation between
different shape phenotypes. Visualizing the movement of microtubules further
reveals a fraction of microtubules that assemble between extra poles and push
the poles apart, suggesting the presence of active extensile force that
prevents pole coalescence. Together, we propose that a proper control over the
magnitude and location of the extensile, pole-pushing force is crucial for
establishing spindle bipolarity while preventing multipolarity.


  Certain subterranean environments of Earth have naturally accumulated
long-lived radionuclides, such as 238U, 232Th and 40K, near the presence of
liquid water. In these natural radioactive environments (NRE), water radiolysis
can produce chemical species of biological importance, such as H2. Although the
proposal of radioactive decay as an alternative source of energy for living
systems has existed for more than thirty years, this hypothesis gained strength
after the recent discovery of a peculiar ecosystem in a gold mine in South
Africa, whose existence is dependent on chemical species produced by water
radiolysis. In this work, we calculate the chemical disequilibrium generated
locally by water radiolysis due gamma radiation and analyse the possible
contribution of this disequilibrium for the emergence of life, considering
conditions of early Earth and having as reference the alkaline hydrothermal
vent (AHV) theory. Results from our kinetic model points out the similarities
between the conditions caused by water radiolysis and those found on alkaline
hydrothermal systems. Our model produces a steady increase of pH with time,
which favours the precipitation of minerals with catalytic activity for
protometabolism, as well as a natural electrochemical gradient in this aqueous
environment. In conclusion, we described a possible free-energy conversion
mechanism that could be a requisite for emergence of life in Hadean Earth.


  Despite substantial progress in non-equilibrium physics, steady-state (s.s.)
probabilities remain intractable to analysis. For a Markov process, s.s.
probabilities can be expressed in terms of transition rates using the
Matrix-Tree theorem (MTT) in the graph-based linear framework. The MTT reveals
that, away from equilibrium, s.s. probabilities become globally dependent on
all rates, with expressions growing exponentially in the system size. This
overwhelming complexity and lack of thermodynamic interpretation have greatly
impeded analysis. Here, we show that s.s. probabilities are proportional to the
average of $\exp(-S(P))$, where $S(P)$ is the entropy generated along minimal
paths, $P$, in the graph, and the average is taken over a probability
distribution on spanning trees. Assuming Arrhenius rates, this "arboreal"
distribution becomes Boltzmann-like, with the energy of a tree being its total
edge barrier energy. This reformulation offers a thermodynamic interpretation
that smoothly generalises equilibrium statistical mechanics and reorganises the
expression complexity: the number of distinct minimal-path entropies depends on
the entropy production index, a new graph-theoretic measure of non-equilibrium
complexity, not on graph size. We demonstrate the power of this reformulation
by extending Hopfield's analysis of discrimination by kinetic proofreading to
any graph with index 1. We derive a general formula for the error ratio and use
it to show that local energy dissipation can yield optimal discrimination
through global synergy.


  In a world of hardening borders, nations may deprive themselves of enjoying
the benefits of cooperative immigrants. Here, we analyze the effect of
efficient cooperative immigrants on a population playing public goods games. We
considered a population structured on a square lattice with individuals playing
public goods games with their neighbors. The demographics are determined by
stochastic birth, death, and migration. The strategies spread through imitation
dynamics. Our model shows that cooperation among natives can emerge due to
social contagion of good role-model agents that can provide better quality
public goods. Only a small fraction of efficient cooperators, among immigrants,
is enough to trigger cooperation across the native population. We see that
native cooperation achieves its peak at moderate values of immigration rate.
Such efficient immigrant cooperators act as nucleation centers for the growth
of cooperative clusters, that eventually dominate defection.


  The processes and mechanisms underlying the origin and maintenance of
biological diversity have long been of central importance in ecology and
evolution. The competitive exclusion principle states that the number of
coexisting species is limited by the number of resources, or by the species'
similarity in resource use. Natural systems such as the extreme diversity of
unicellular life in the oceans provide counter examples. It is known that
mathematical models incorporating population fluctuations can lead to
violations of the exclusion principle. Here we use simple eco-evolutionary
models to show that a certain type of population dynamics, boom-bust dynamics,
can allow for the evolution of much larger amounts of diversity than would be
expected with stable equilibrium dynamics. Boom-bust dynamics are characterized
by long periods of almost exponential growth (boom) and a subsequent population
crash due to competition (bust). When such ecological dynamics are incorporated
into an evolutionary model that allows for adaptive diversification in
continuous phenotype spaces, desynchronization of the boom-bust cycles of
coexisting species can lead to the maintenance of high levels of diversity.


  Across the tree of life, populations have evolved the capacity to contend
with suboptimal conditions by engaging in dormancy, whereby individuals enter a
reversible state of reduced metabolic activity. The resulting seed banks are
complex, storing information and imparting memory that gives rise to
multi-scale structures and networks spanning collections of cells to entire
ecosystems. We outline the fundamental attributes and emergent phenomena
associated with dormancy and seed banks, with the vision for a unifying and
mathematically based framework that can address problems in the life sciences,
ranging from global change to cancer biology.


  Honey bees make decisions regarding foraging and nest-site selection in
groups ranging from hundreds to thousands of individuals. To effectively make
these decisions bees need to communicate within a spatially distributed group.
However, the spatiotemporal dynamics of honey bee communication have been
mostly overlooked in models of collective decisions, focusing primarily on mean
field models of opinion dynamics. We analyze how the spatial properties of the
nest or hive, and the movement of individuals with different belief states
(uncommitted or committed) therein affect the rate of information transmission
using spatially-extended models of collective decision-making within a hive.
Honeybees waggle-dance to recruit conspecifics with an intensity that is a
threshold nonlinear function of the waggler concentration. Our models range
from treating the hive as a chain of discrete patches to a continuous line
(long narrow hive). The combination of population-thresholded recruitment and
compartmentalized populations generates tradeoffs between rapid information
propagation with strong population dispersal and recruitment failures resulting
from excessive population diffusion and also creates an effective colony-level
signal-detection mechanism whereby recruitment to low quality objectives is
blocked.


  For protein sequence datasets, unlabeled data has greatly outpaced labeled
data due to the high cost of wet-lab characterization. Recent deep-learning
approaches to protein prediction have shown that pre-training on unlabeled data
can yield useful representations for downstream tasks. However, the optimal
pre-training strategy remains an open question. Instead of strictly borrowing
from natural language processing (NLP) in the form of masked or autoregressive
language modeling, we introduce a new pre-training task: directly predicting
protein profiles derived from multiple sequence alignments. Using a set of
five, standardized downstream tasks for protein models, we demonstrate that our
pre-training task along with a multi-task objective outperforms masked language
modeling alone on all five tasks. Our results suggest that protein sequence
models may benefit from leveraging biologically-inspired inductive biases that
go beyond existing language modeling techniques in NLP.


  During their lifetimes, individuals in populations pass through different
states, and the notion of an occupancy time describes the amount of time an
individual spends in a given set of states. Questions related to this idea were
studied in a recent paper by Roth and Caswell for cases where the environmental
conditions are constant. However, it is truly important to consider the case
where environments are changing randomly or in directional way through time, so
the transition probabilities between different states change over time,
motivating the use of time-dependent stage-structured models.
  Using absorbing inhomogenous Markov chains and the discrete-time
McKendrick--von F{\"o}rster equation, we derive explicit formulas for the
occupancy time, its expectation, and its higher-order moments for
stage-structured models with time-dependent transition rates. We apply our
approach to study a time-dependent model of the Southern Fulmar, and obtain
insights into how the number of breeding attempts depends on external
conditions that vary through time.


  Regulatory T cells (Tregs) play a crucial role in mediating immune response.
Yet an algorithmic understanding of the role of Tregs in adaptive immunity
remains lacking. Here, we present a biophysically realistic model of Treg
mediated self-tolerance in which Tregs bind to self-antigens and locally
inhibit the proliferation of nearby activated T cells. By exploiting a duality
between ecological dynamics and constrained optimization, we show that Tregs
tile the potential antigen space while simultaneously minimizing the overlap
between Treg activation profiles. We find that for sufficiently high Treg
diversity, Treg mediated self-tolerance is robust to fluctuations in
self-antigen concentrations but lowering the Treg diversity results in a sharp
transition -- related to the Gardner transition in perceptrons -- to a regime
where changes in self-antigen concentrations can result in an auto-immune
response. We propose a novel experimental test of this transition in
immune-deficient mice and discuss potential implications for autoimmune
diseases.


  The study of the origins of crops is of interest from both a fundamental
evolutionary understanding viewpoint, and from an applied agricultural
technology perspective. The date palm (Phoenix dactylifera L.) is the iconic
fruit crop of hot and arid regions of North Africa and the Middle East,
producing sugar-rich fruits, known as dates. There are many different cultivars
each with distinctive fruit traits, and there are many wild Phoenix species
too, which in total form a complex of related species. The understanding of
plant domestication involves multiple disciplines, including phylogeography,
population genetics and archaeology. In the past decade, they have prompted new
discoveries on the evolutionary history of date palm, but a complete
understanding of its origins remains to be elucidated, along with the genetic
architecture of its domestication syndrome. In this chapter, we review the
current state of the art regarding the origins of the domesticated date palm.
We first discuss whether date palms are domesticated, and highlight how they
diverge from their wild Phoenix relatives. We then outline patterns in the
population genetic and archaeobotanical data, and review different models for
the origins of domesticated date palms by highlighting sources of evidence that
are either consistent or inconsistent with each model. We then review the
process of date palm domestication, and emphasize the human activities that
have prompted its domestication. We particularly focus on the evolution of
fruit traits.


  Variational Bayesian phylogenetic inference (VBPI) provides a promising
general variational framework for efficient estimation of phylogenetic
posteriors. However, the current diagonal Lognormal branch length approximation
would significantly restrict the quality of the approximating distributions. In
this paper, we propose a new type of VBPI, VBPI-NF, as a first step to empower
phylogenetic posterior estimation with deep learning techniques. By handling
the non-Euclidean branch length space of phylogenetic models with carefully
designed permutation equivariant transformations, VBPI-NF uses normalizing
flows to provide a rich family of flexible branch length distributions that
generalize across different tree topologies. We show that VBPI-NF significantly
improves upon the vanilla VBPI on a benchmark of challenging real data Bayesian
phylogenetic inference problems. Further investigation also reveals that the
structured parameterization in those permutation equivariant transformations
can provide additional amortization benefit.


  Zhou et al. reported the discovery of RmYN02, a strain closely related to
SARS-CoV-2, which is claimed to contain a natural PAA amino acid insertion at
the S1/S2 junction of the spike protein at the same position of the PRRA
insertion that has created a polybasic furin cleavage site in SARS-CoV-2. The
authors support with their findings the theory that the furin cleavage site
insertion present in SARS-CoV-2 is natural. Because no nucleotide alignment
with closely related strains of the region coding for the supposed insertion is
provided by Zhou et al., we have applied several alignment algorithms to search
for the most parsimonious alignments. We conclude that RmYN02 does not contain
an insertion at the S1/S2 junction when compared to its closest relatives at
the nucleotide level, but rather a 6-nucleotide deletion and that the claimed
PAA insertion is more likely to be the result of mutations. A close examination
of RmYN02 sequencing records and assembly methods is wishful. In conclusion,
SARS-CoV-2, with its 12-nucleotide insertion at the S1/S2 junction remains
unique among its sarbecovirus relatives.


  A great variety of complex physical, natural and artificial systems are
governed by statistical distributions, which often follow a standard
exponential function in the bulk, while their tail obeys the Pareto power law.
The recently introduced $\kappa$-statistics framework predicts distribution
functions with this feature. A growing number of applications in different
fields of investigation are beginning to prove the relevance and effectiveness
of $\kappa$-statistics in fitting empirical data. In this paper, we use
$\kappa$-statistics to formulate a statistical approach for epidemiological
analysis. We validate the theoretical results by fitting the derived
$\kappa$-Weibull distributions with data from the plague pandemic of 1417 in
Florence as well as data from the COVID-19 pandemic in China over the entire
cycle that concludes in April 16, 2020. As further validation of the proposed
approach we present a more systematic analysis of COVID-19 data from countries
such as Germany, Italy, Spain and United Kingdom, obtaining very good agreement
between theoretical predictions and empirical observations. For these countries
we also study the entire first cycle of the pandemic which extends until the
end of July 2020. The fact that both the data of the Florence plague and those
of the Covid-19 pandemic are successfully described by the same theoretical
model, even though the two events are caused by different diseases and they are
separated by more than 600 years, is evidence that the $\kappa$-Weibull model
has universal features.


  Compound-protein pairs dominate FDA-approved drug-target pairs and the
prediction of compound-protein affinity and contact (CPAC) could help
accelerate drug discovery. In this study we consider proteins as multi-modal
data including 1D amino-acid sequences and (sequence-predicted) 2D residue-pair
contact maps. We empirically evaluate the embeddings of the two single
modalities in their accuracy and generalizability of CPAC prediction (i.e.
structure-free interpretable compound-protein affinity prediction). And we
rationalize their performances in both challenges of embedding individual
modalities and learning generalizable embedding-label relationship. We further
propose two models involving cross-modality protein embedding and establish
that the one with cross interaction (thus capturing correlations among
modalities) outperforms SOTAs and our single modality models in affinity,
contact, and binding-site predictions for proteins never seen in the training
set.


  In the present article, we investigate the effects of dormancy on an abstract
population genetic level. We first provide a short review of seed bank models
in population genetics, and the role of dormancy for the interplay of
evolutionary forces in general, before we discuss two recent paradigmatic
models, referring to spontaneous resp. simultaneous switching of individuals
between the active and the dormant state. We show that both mechanisms give
rise to non-trivial mathematical objects, namely the (continuous) seed bank
diffusion and the seed bank diffusion with jumps, as well as their dual
processes, the seed bank coalescent and the seed bank coalescent with
simultaneous switching.


  The spread of infectious diseases is a major challenge in our contemporary
world, especially after the recent outbreak of Coronavirus disease 2019
(COVID-19). The quarantine strategy is one of the important intervention
measures to control the spread of an epidemic by greatly minimizing the
likelihood of contact between infected and susceptible individuals. In this
study, we analyze the impact of various stochastic disturbances on the epidemic
dynamics during the quarantine period. For this purpose, we present an SIQS
epidemic model that incorporates the stochastic transmission and the L\'evy
noise in order to simulate both small and massive perturbations. Under
appropriate conditions, some interesting asymptotic properties are proved,
namely: ergodicity, persistence in the mean, and extinction of the disease. The
theoretical results show that the dynamics of the perturbed model are
determined by parameters that are closely related to the stochastic noises. Our
work improves many existing studies in the field of mathematical epidemiology
and provides new techniques to predict and analyze the dynamic behavior of
epidemics.


  We outline recent developments in artificial intelligence (AI) and machine
learning (ML) techniques for integrative structural biology of intrinsically
disordered proteins (IDP) ensembles. IDPs challenge the traditional protein
structure-function paradigm by adapting their conformations in response to
specific binding partners leading them to mediate diverse, and often complex
cellular functions such as biological signaling, self organization and
compartmentalization. Obtaining mechanistic insights into their function can
therefore be challenging for traditional structural determination techniques.
Often, scientists have to rely on piecemeal evidence drawn from diverse
experimental techniques to characterize their functional mechanisms. Multiscale
simulations can help bridge critical knowledge gaps about IDP structure
function relationships - however, these techniques also face challenges in
resolving emergent phenomena within IDP conformational ensembles. We posit that
scalable statistical inference techniques can effectively integrate information
gleaned from multiple experimental techniques as well as from simulations, thus
providing access to atomistic details of these emergent phenomena.


  The SARS-CoV-2 pandemic has created a global race for a cure. One approach
focuses on designing a novel variant of the human angiotensin-converting enzyme
2 (ACE2) that binds more tightly to the SARS-CoV-2 spike protein and diverts it
from human cells. Here we formulate a novel protein design framework as a
reinforcement learning problem. We generate new designs efficiently through the
combination of a fast, biologically-grounded reward function and sequential
action-space formulation. The use of Policy Gradients reduces the compute
budget needed to reach consistent, high-quality designs by at least an order of
magnitude compared to standard methods. Complexes designed by this method have
been validated by molecular dynamics simulations, confirming their increased
stability. This suggests that combining leading protein design methods with
modern deep reinforcement learning is a viable path for discovering a Covid-19
cure and may accelerate design of peptide-based therapeutics for other
diseases.


  Large-scale testing is considered key to assess the state of the current
COVID-19 pandemic. Yet, the link between the reported case numbers and the true
state of the pandemic remains elusive. We develop mathematical models based on
competing hypotheses regarding this link, thereby providing different
prevalence estimates based on case numbers, and validate them by predicting
SARS-CoV-2-attributed death rate trajectories. Assuming that individuals were
tested based solely on a predefined risk of being infectious implies the
absolute case numbers reflect the prevalence, but turned out to be a poor
predictor, consistently overestimating growth rates at the beginning of two
COVID-19 epidemic waves. In contrast, assuming that testing capacity is fully
exploited performs better. This leads to using the percent-positive rate as a
more robust indicator of epidemic dynamics, however we find it is subject to a
saturation phenomenon that needs to be accounted for as the number of tests
becomes larger.


  Non-pharmaceutical interventions (NPIs) have played a crucial role in
controlling the spread of COVID-19. Nevertheless, NPI efficacy varies
enormously between and within countries, mainly because of population and
behavioural heterogeneity. In this work, we adapted a multi-group SEIRA model
to study the spreading dynamics of COVID-19 in Chile, representing
geographically separated regions of the country by different groups. We use
national mobilization statistics to estimate the connectivity between regions
and data from governmental repositories to obtain COVID-19 spreading and death
rates in each region. We then assessed the effectiveness of different NPIs by
studying the temporal evolution of the reproduction number Rt. Analyzing
data-driven and model-based estimates of Rt, we found a strong coupling of
different regions, highlighting the necessity of organized and coordinated
actions to control the spread of SARS-CoV-2. Finally, we evaluated different
scenarios to forecast the evolution of COVID-19 in the most densely populated
regions, finding that the early lifting of restriction probably will lead to
novel outbreaks.


  Graph neural networks for molecular property prediction are frequently
underspecified by data and fail to generalise to new scaffolds at test time. A
potential solution is Bayesian learning, which can capture our uncertainty in
the model parameters. This study benchmarks a set of Bayesian methods applied
to a directed MPNN, using the QM9 regression dataset. We find that capturing
uncertainty in both readout and message passing parameters yields enhanced
predictive accuracy, calibration, and performance on a downstream molecular
search task.


  This book brings new mathematical rigour to the ongoing vigorous debate on
how to quantify biological diversity. The question "what is diversity?" has
surprising mathematical depth, and breadth too: this book involves parts of
mathematics ranging from information theory, functional equations and
probability theory to category theory, geometric measure theory and number
theory. It applies the power of the axiomatic method to a biological problem of
pressing concern, but the new concepts and theorems are also motivated from a
purely mathematical perspective.
  The main narrative thread requires no more than an undergraduate course in
analysis. No familiarity with entropy or diversity is assumed.


  The 2019 novel coronavirus (SARS-CoV-2) pandemic has resulted in more than a
million deaths, high morbidities, and economic distress worldwide. There is an
urgent need to identify medications that would treat and prevent novel diseases
like the 2019 coronavirus disease (COVID-19). Drug repurposing is a promising
strategy to discover new medical indications of the existing approved drugs due
to several advantages in terms of the costs, safety factors, and quick results
compared to new drug design and discovery. In this work, we explore
computational data-driven methods for drug repurposing and propose a dedicated
graph neural network (GNN) based drug repurposing model, called Dr-COVID.
Although we analyze the predicted drugs in detail for COVID-19, the model is
generic and can be used for any novel diseases. We construct a four-layered
heterogeneous graph to model the complex interactions between drugs, diseases,
genes, and anatomies. We pose drug repurposing as a link prediction problem.
Specifically, we design an encoder based on the scalable inceptive graph neural
network (SIGN) to generate embeddings for all the nodes in the four-layered
graph and propose a quadratic norm scorer as a decoder to predict treatment for
a disease. We provide a detailed analysis of the 150 potential drugs (such as
Dexamethasone, Ivermectin) predicted by Dr-COVID for COVID-19 from different
pharmacological classes (e.g., corticosteroids, antivirals, antiparasitic). Out
of these 150 drugs, 46 drugs are currently in clinical trials. Dr-COVID is
evaluated in terms of its prediction performance and its ability to rank the
known treatment drugs for diseases as high as possible. For a majority of the
diseases, Dr-COVID ranks the actual treatment drug in the top 15.


  A Boolean network (BN) with $n$ components is a discrete dynamical system
described by the successive iterations of a function $f:\{0,1\}^n \to
\{0,1\}^n$. This model finds applications in biology, where fixed points play a
central role. For example, in genetic regulations, they correspond to cell
phenotypes. In this context, experiments reveal the existence of positive or
negative influences among components: component $i$ has a positive (resp.
negative) influence on component $j$ meaning that $j$ tends to mimic (resp.
negate) $i$. The digraph of influences is called signed interaction digraph
(SID), and one SID may correspond to a large number of BNs (which is, in
average, doubly exponential according to $n$). The present work opens a new
perspective on the well-established study of fixed points in BNs. When
biologists discover the SID of a BN they do not know, they may ask: given that
SID, can it correspond to a BN having at least/at most $k$ fixed points?
Depending on the input, we prove that these problems are in $\textrm{P}$ or
complete for $\textrm{NP}$, $\textrm{NP}^{\textrm{NP}}$,
$\textrm{NP}^{\textrm{#P}}$ or $\textrm{NEXPTIME}$. In particular, we prove
that it is $\textrm{NP}$-complete (resp. $\textrm{NEXPTIME}$-complete) to
decide if a given SID can correspond to a BN having at least two fixed points
(resp. no fixed point).


  Epidemic spreading can be suppressed by the introduction of containment
measures such as social distancing and lock downs. Yet, when such measures are
relaxed, new epidemic waves and infection cycles may occur. Here we explore
this issue in compartmentalized epidemic models on graphs in presence of a
feedback between the infection state of the population and the structure of its
social network for the case of discontinuous control. We show that in random
graphs the effect of containment measures is simply captured by a
renormalization of the effective infection rate that accounts for the change in
the branching ratio of the network. In our simple setting, a piece-wise
mean-field approximations can be used to derive analytical formulae for the
number of epidemic waves and their length. A variant of the model with
imperfect information is used to model data of the recent covid-19 epidemics in
the Basque Country and Lombardy, where we estimate the extent of social network
disruption during lock downs and characterize the dynamical attractors.


  We consider the transmission dynamics of COVID-19 which is characterized by
two distinct features. One is the existence of asymptomatic carriers which is a
hidden variable in the problem. The other is the issue of latency which means
that among the symptomatic carriers there could be a fraction whose symptoms
develop after a couple of days. Our modelling is restricted to what we call the
Phase -1 of the disease. During this phase the disease sets in and the number
of infected people starts growing fast ( the number of new cases per day keeps
growing on an average ) and then it slows down ( the number of new cases per
day starts decreasing ) with the number of new cases decreasing to about one
tenth of its peak value or even smaller). We define Phase-1 to be over when the
daily cases start rising once again. We write down a Henon-like map to take
various effects into account for this first phase.


  How we choose to represent our data has a fundamental impact on our ability
to subsequently extract information from them. Machine learning promises to
automatically determine efficient representations from large unstructured
datasets, such as those arising in biology. However, empirical evidence
suggests that seemingly minor changes to these machine learning models yield
drastically different data representations that result in different biological
interpretations of data. This begs the question of what even constitutes the
most meaningful representation. Here, we approach this question for
representations of protein sequences, which have received considerable
attention in the recent literature. We explore two key contexts in which
representations naturally arise: transfer learning and interpretable learning.
In the first context, we demonstrate that several contemporary practices yield
suboptimal performance, and in the latter we demonstrate that taking
representation geometry into account significantly improves interpretability
and lets the models reveal biological information that is otherwise obscured.


  We obtain general inequalities constraining the difference between the
average of an arbitrary function of a phenotypic trait, which includes the
fitness landscape of the trait itself, in the presence or in the absence of
natural selection. These inequalities imply bounds on the strength of
selection, which can be measured from the statistics of trait values and
divisions along lineages. The upper bound is related to recent generalizations
of linear response relations in Stochastic Thermodynamics, and shares common
features with Fisher's fundamental theorem of natural selection, and with its
generalization by Price, although they define different measures of selection.
The lower bound follows from recent improvements on Jensen's inequality, and
both bounds depend on the variability of the fitness landscape. We illustrate
our results using numerical simulations of growing cell colonies and with
experimental data of time-lapse microscopy experiments of bacteria cell
colonies.


  Between-species coevolution, and in particular antagonistic host-parasite
coevolution, is a major process shaping within-species diversity. In this paper
we investigate the role of various stochastic processes affecting the outcome
of the deterministic coevolutionary models. Specifically, we assess 1) the
impact of genetic drift and mutation on the maintenance of polymorphism at the
interacting loci, and 2) the change in neutral allele frequencies across the
genome of both coevolving species due to co-demographic population size
changes. We find that genetic drift decreases the likelihood to observe classic
balancing selection signatures, and that for most realistic values of the
coevolutionary parameters, balancing selection signatures cannot be seen at the
host loci. Further, we reveal that contrary to classic expectations, fast
changes in parasite population size due to eco-evo feedbacks can be tracked by
the allelic site-frequency spectrum measured at several time points. Changes in
host population size are, however, less pronounced and thus not observable.
Finally, we also review several understudied stochastic processes occurring in
host-parasite coevolution which are of importance to predict maintenance of
polymorphism at the underlying loci and the genome-wide nucleotide diversity of
host and parasite populations.


  The COVID-19 pandemic has challenged authorities at different levels of
government administration around the globe. When faced with diseases of this
severity, it is useful for the authorities to have prediction tools to estimate
in advance the impact on the health system and the human, material, and
economic resources that will be necessary. In this paper, we construct an
extended Susceptible-Exposed-Infected-Recovered model that incorporates the
social structure of Mar del Plata, the $4^\circ$ most inhabited city in
Argentina and head of the Municipality of General Pueyrred\'on. Moreover, we
consider detailed partitions of infected individuals according to the illness
severity, as well as data of local health resources, to bring these predictions
closer to the local reality. Tuning the corresponding epidemic parameters for
COVID-19, we study an alternating quarantine strategy, in which a part of the
population can circulate without restrictions at any time, while the rest is
equally divided into two groups and goes on successive periods of normal
activity and lockdown, each one with a duration of $\tau$ days. Besides, we
implement a random testing strategy over the population. We found that $\tau =
7$ is a good choice for the quarantine strategy since it matches with the
weekly cycle as it reduces the infected population. Focusing on the health
system, projecting from the situation as of September 30, we foresee a
difficulty to avoid saturation of ICU, given the extremely low levels of
mobility that would be required. In the worst case, our model estimates that
four thousand deaths would occur, of which 30\% could be avoided with proper
medical attention. Nonetheless, we found that aggressive testing would allow an
increase in the percentage of people that can circulate without restrictions,
being the equipment required to deal with the additional critical patients
relatively low.


  Temperature sensing is a ubiquitous cell behavior, but the fundamental limits
to the precision of temperature sensing are poorly understood. Unlike in
chemical concentration sensing, the precision of temperature sensing is not
limited by extrinsic fluctuations in the temperature field itself. Instead, we
find that precision is limited by the intrinsic copy number, turnover, and
binding kinetics of temperature-sensitive proteins. Developing a model based on
the canonical TlpA protein, we find that a cell can estimate temperature to
within 2%. We compare this prediction with in vivo data on temperature sensing
in bacteria.


  The aim of this paper is two-fold. First, we propose a new computational
method to investigate the particularities of evolution. Second, we apply this
method to a model of gene regulatory networks (GRNs) and explore the evolution
of mutational robustness and bistability. Living systems have developed their
functions through evolutionary processes. To understand the particularities of
this process theoretically, evolutionary simulation (ES) alone is insufficient
because the outcomes of ES depend on evolutionary pathways. We need a reference
system for comparison. An appropriate reference system for this purpose is an
ensemble of the randomly sampled genotypes. However, generating high-fitness
genotypes by simple random sampling is difficult because such genotypes are
rare. In this study, we used the multicanonical Monte Carlo method developed in
statistical physics to construct a reference ensemble of GRNs and compared it
with the outcomes of ES. We obtained the following results. First, mutational
robustness was significantly higher in ES than in the reference ensemble at the
same fitness level. Second, the emergence of a new phenotype, bistability, was
delayed in evolution. Third, the bistable group of GRNs contains many
mutationally fragile GRNs compared with those in the non-bistable group. This
suggests that the delayed emergence of bistability is a consequence of the
mutation-selection mechanism.


  Abrupt transitions leading to algal blooms are quite well known in aquatic
ecosystems and have important implications for the environment. These ecosystem
shifts have been largely attributed to nutrient dynamics and food web
interactions. Contamination with heavy metals such as copper can modulate such
ecological interactions which in turn may impact ecosystem functioning.
Motivated by this, we explored the effect of copper enrichment on such regime
shifts in planktonic systems. We integrated copper contamination to a minimal
phytoplankton-zooplankton model which is known to demonstrate abrupt
transitions between ecosystem states. Our results suggest that both the toxic
and deficient concentration of copper in water bodies can lead to regime shift
to an algal dominated alternative stable state. Further, interaction with fish
density can also lead to collapse of population cycles thus leading to algal
domination in the intermediate copper ranges. Environmental stochasticity may
result in state transition much prior to the tipping point and there is a
significant loss in the bimodality on increasing intensity and redness of
noise. Finally, the impending state shifts due to contamination cannot be
predicted by the generic early warning indicators unless the transition is
close enough. Overall the study provides fresh impetus to explore regime shifts
in ecosystems under the influence of anthropogenic changes like chemical
contamination.


  Less than 1% of protein sequences are structurally and functionally
annotated. Natural Language Processing (NLP) community has recently embraced
self-supervised learning as a powerful approach to learn representations from
unlabeled text, in large part due to the attention-based context-aware
Transformer models. In this work we present a modification to the RoBERTa model
by inputting during pre-training a mixture of binding and non-binding protein
sequences (from STRING database). However, the sequence pairs have no label to
indicate their binding status, as the model relies solely on Masked Language
Modeling (MLM) objective during pre-training. After fine-tuning, such approach
surpasses models trained on single protein sequences for protein-protein
binding prediction, TCR-epitope binding prediction, cellular-localization and
remote homology classification tasks. We suggest that the Transformer's
attention mechanism contributes to protein binding site discovery. Furthermore,
we compress protein sequences by 64% with the Byte Pair Encoding (BPE)
vocabulary consisting of 10K subwords, each around 3-4 amino acids long.
Finally, to expand the model input space to even larger proteins and
multi-protein assemblies, we pre-train Longformer models that support 2,048
tokens. Further work in token-level classification for secondary structure
prediction is needed. Code available at:
https://github.com/PaccMann/paccmann_proteomics