  The breakdown of trusted sources of information is probably one of the most
serious problems today, since in the absence of a common ground, it will be
impossible to address the problems that trouble our contemporary world. The
Covid-19 pandemic is just a recent situation where the lack of agreed stances
has led to failure and hopelessness. In fact, disinformation surrounding the
Covid-19 has been a distinctive feature of this pandemic since its very
beginning and has hampered what is perhaps the most important initiative to
prevent the spread of the coronavirus, viz., an effective communication between
scientifically-minded health authorities and the general public. To investigate
how disinformation threatens epistemic security, here we propose and solve
analytically an evolutionary game-theoretic model where the individuals must
accurately estimate some property of their hazardous environment. They can
either explore the environment or copy the estimate from another individual,
who may display a distorted version of its estimate. We find that the
exploration-only strategy is optimal when the environment is relatively safe
and the individuals are not reliable. In this doomsday scenario, disinformation
erodes trust and suppresses the ability of the individuals to share information
with one another.


  Competitive systems can exhibit both hierarchical (transitive) and cyclic
(intransitive) structures. Despite theoretical interest in cyclic competition,
which offers richer dynamics, and occupies a larger subset of the space of
possible competitive systems, most real-world systems are predominantly
transitive. Why? Here, we introduce a generic mechanism which promotes
transitivity, even when there is ample room for cyclicity. Consider a
competitive system where outcomes are mediated by competitor attributes via a
performance function. We demonstrate that, if competitive outcomes depend
smoothly on competitor attributes, then similar competitors compete
transitively. We quantify the rate of convergence to transitivity given the
similarity of the competitors and the smoothness of the performance function.
Thus, we prove the adage regarding apples and oranges. Similar objects admit
well ordered comparisons. Diverse objects may not. To test that theory, we run
a series of evolution experiments designed to mimic genetic training
algorithms. We consider a series of canonical bimatrix games and an ensemble of
random performance functions that demonstrate the generality of our mechanism,
even when faced with highly cyclic games. We vary the training parameters
controlling the evolution process, and the shape parameters controlling the
performance function, to evaluate the robustness of our results. These
experiments illustrate that, if competitors evolve to optimize performance,
then their traits may converge, leading to transitivity.


  Accurate ADMET (an abbreviation for "absorption, distribution, metabolism,
excretion, and toxicity") predictions can efficiently screen out undesirable
drug candidates in the early stage of drug discovery. In recent years, multiple
comprehensive ADMET systems that adopt advanced machine learning models have
been developed, providing services to estimate multiple endpoints. However,
those ADMET systems usually suffer from weak extrapolation ability. First, due
to the lack of labelled data for each endpoint, typical machine learning models
perform frail for the molecules with unobserved scaffolds. Second, most systems
only provide fixed built-in endpoints and cannot be customised to satisfy
various research requirements. To this end, we develop a robust and endpoint
extensible ADMET system, HelixADMET (H-ADMET). H-ADMET incorporates the concept
of self-supervised learning to produce a robust pre-trained model. The model is
then fine-tuned with a multi-task and multi-stage framework to transfer
knowledge between ADMET endpoints, auxiliary tasks, and self-supervised tasks.
Our results demonstrate that H-ADMET achieves an overall improvement of 4%,
compared with existing ADMET systems on comparable endpoints. Additionally, the
pre-trained model provided by H-ADMET can be fine-tuned to generate new and
customised ADMET endpoints, meeting various demands of drug research and
development requirements.


  Molecular dynamics (MD) simulations allow atomistic insights into chemical
and biological processes. Accurate MD simulations require computationally
demanding quantum-mechanical calculations, being practically limited to short
timescales and few atoms. For larger systems, efficient, but much less reliable
empirical force fields are used. Recently, machine learned force fields (MLFFs)
emerged as an alternative means to execute MD simulations, offering similar
accuracy as ab initio methods at orders-of-magnitude speedup. Until now, MLFFs
mainly capture short-range interactions in small molecules or periodic
materials, due to the increased complexity of constructing models and obtaining
reliable reference data for large molecules, where long-ranged many-body
effects become important. This work proposes a general approach to constructing
accurate MLFFs for large-scale molecular simulations (GEMS) by training on
"bottom-up" and "top-down" molecular fragments of varying size, from which the
relevant physicochemical interactions can be learned. GEMS is applied to study
the dynamics of alanine-based peptides and the 46-residue protein crambin in
aqueous solution, allowing nanosecond-scale MD simulations of >25k atoms at
essentially ab initio quality. Our findings suggest that structural motifs in
peptides and proteins are more flexible than previously thought, indicating
that simulations at ab initio accuracy might be necessary to understand dynamic
biomolecular processes such as protein (mis)folding, drug-protein binding, or
allosteric regulation.


  Stock and flow diagrams are widely used in epidemiology to model the dynamics
of populations. Although tools already exist for building these diagrams and
simulating the systems they describe, we have created a new package called
StockFlow, part of the AlgebraicJulia ecosystem, which uses ideas from category
theory to overcome notable limitations of existing software. Compositionality
is provided by the theory of decorated cospans: stock and flow diagrams can be
composed to form larger ones in an intuitive way formalized by the operad of
undirected wiring diagrams. Our approach also cleanly separates the syntax of
stock and flow diagrams from the semantics they can be assigned. We consider
semantics in ordinary differential equations, although others are possible. As
an example, we explain code in StockFlow that implements a simplified version
of a COVID-19 model used in Canada.


  A steadily growing computational power is employed to perform molecular
dynamics simulations of biological macromolecules, which represents at the same
time an immense opportunity and a formidable challenge. In fact, large amounts
of data are produced, from which useful, synthetic, and intelligible
information has to be extracted to make the crucial step from knowing to
understanding. Here we tackled the problem of coarsening the conformational
space sampled by proteins in the course of molecular dynamics simulations. We
applied different schemes to cluster the frames of a dataset of protein
simulations; we then employed an information-theoretical framework, based on
the notion of resolution and relevance, to gauge how well the various
clustering methods accomplish this simplification of the configurational space.
Our approach allowed us to identify the level of resolution that optimally
balances simplicity and informativeness; furthermore, we found that the most
physically accurate clustering procedures are those that induce an ultrametric
structure of the low-resolution space, consistently with the hypothesis that
the protein conformational landscape has a self-similar organisation. The
proposed strategy is general and its applicability extends beyond that of
computational biophysics, making it a valuable tool to extract useful
information from large datasets.


  Protein Structure Prediction (PSP) is an unsolved problem in the field of
computational biology. The problem of protein structure prediction is about
predicting the native conformation of a protein, while its sequence of amino
acids is known. Regarding processing limitations of current computer systems,
all-atom simulations for proteins are typically unpractical; several reduced
models of proteins have been proposed. Additionally, due to intrinsic hardness
of calculations even in reduced models, many computational methods mainly based
on artificial intelligence have been proposed to solve the problem. Agent-based
modeling is a relatively new method for modeling systems composed of
interacting items. In this paper we proposed a new approach for protein
structure prediction by using agent-based modeling (ABM) in two dimensional
hydrophobic-hydrophilic model. We broke the whole process of protein structure
prediction into two steps: the first step, which was introduced in our previous
paper, is about biasing the linear sequence to gain a primary energy, and the
next step, which will be explained in this paper, is about using ABM with a
predefined set of rules, to find the best conformation in the least possible
amount of time and steps. This method was implemented in NETLOGO. We have
tested this algorithm on several benchmark sequences ranging from 20 to 50-mers
in two dimensional Hydrophobic-Hydrophilic lattice models. Comparing to the
result of the other algorithms, our method is capable of finding the best known
conformations in a significantly shorter time. A major problem in PSP
simulation is that as the sequence length increases the time consumed to
predict a valid structure will exponentially increase. In contrast, by using
MAS2HP the effect of increase in sequence length on spent time has changed from
exponentially to linear.


  In this paper, based on the epidemiological microparasite model, a
parasite-host model is established by considering the fear effect of
susceptible individuals on infectors. We explored the pattern formation with
the help of numerical simulation, and analyzed the effects of fear effect,
infected host mortality, population diffusion rate and reducing reproduction
ability of infected hosts on population activities in different degrees.
Theoretically, we give the general conditions for the stability of the model
under non-diffusion and considering the Turing instability caused by diffusion.
Our results indicate how fear affects the distribution of the uninfected and
infected hosts in the habitat and quantify the influence of the fear factor on
the spatiotemporal pattern of the population. In addition, we analyze the
influence of natural death rate, reproduction ability of infected hosts, and
diffusion level of uninfected (infected) hosts on the spatiotemporal pattern,
respectively. The results present that the growth of pattern induced by
intensified fear effect follows the certain rule: cold spots $\rightarrow$ cold
spots-stripes $\rightarrow$ cold stripes $\rightarrow$ hot stripes
$\rightarrow$ hot spots-stripes $\rightarrow$ hot spots. Interestingly, the
natural mortality and fear effect take the opposite effect on the growth order
of the pattern. From the perspective of biological significance, we find that
the degree of fear effect can reshape the distribution of population to meet
the previous rule.


  Diverse T and B cell repertoires play an important role in mounting effective
immune responses against a wide range of pathogens and malignant cells. The
number of unique T and B cell clones is characterized by T and B cell receptors
(TCRs and BCRs), respectively. Although receptor sequences are generated
probabilistically by recombination processes, clinical studies found a high
degree of sharing of TCRs and BCRs among different individuals. In this work,
we use a general probabilistic model for T/B cell receptor clone abundances to
define "publicness" or "privateness" and information-theoretic measures for
comparing the frequency of sampled sequences observed across different
individuals. We derive mathematical formulae to quantify the mean and the
variances of clone richness and overlap. Our results can be used to evaluate
the effect of different sampling protocols on abundances of clones within an
individual as well as the commonality of clones across individuals. Using
synthetic and empirical TCR amino acid sequence data, we perform simulations to
study expected clonal commonalities across multiple individuals. Based on our
formulae, we compare these simulated results with the analytically predicted
mean and variances of the repertoire overlap. Complementing the results on
simulated repertoires, we derive explicit expressions for the richness and its
uncertainty for specific, single-parameter truncated power-law probability
distributions. Finally, the information loss associated with grouping together
certain receptor sequences, as is done in spectratyping, is also evaluated. Our
approach can be, in principle, applied under more general and mechanistically
realistic clone generation models.


  Emerging marine infectious diseases pose a substantial threat to marine
ecosystems and the conservation of their biodiversity. Compartmental models of
epidemic transmission in marine sessile organisms, available only recently, are
based on non-spatial descriptions in which space is homogenised and parasite
mobility is not explicitly accounted for. However, in realistic scenarios
epidemic transmission is conditioned by the spatial distribution of hosts and
the parasites mobility patterns, calling for a explicit description of space.
In this work we develop a spatially-explicit individual-based model to study
disease transmission by waterborne parasites in sessile marine populations. We
investigate the impact of spatial disease transmission through extensive
numerical simulations and theoretical analysis. Specifically, the effects of
parasite mobility into the epidemic threshold and the temporal progression of
the epidemic are assessed. We show that larger values of pathogen mobility
imply more severe epidemics, as the number of infections increases, and shorter
time-scales to extinction. An analytical expression for the basic reproduction
number of the spatial model is derived as function of the non-spatial
counterpart, which characterises a transition between a disease-free and a
propagation phase, in which the disease propagates over a large fraction of the
system.


  This work focuses on the formation of mating aggregates in zooplankton. In
particular, sexual encounters are behaviourally supported by males actively
swimming in search for females, and approaching them for mating once they are
found. While the random search leads to a diffusive flux of individuals, the
approaching for encounter supports attraction. Thus, we ask whether these
competing mechanisms of diffusion and attraction can support aggregation and
lead to the formation of mating clusters. To answer our question we formulate a
model in which particles performing random walks can briefly make contact with
other particles if they are found within a particular distance from each other.
Our analysis shows that this model supports clustering in a way analogous to
the process of colloid aggregation. Following that, we analyze a dataset of 3D
trajectories of swimming copepods and show that the results compare well with
our model. These results support the hypothesis that pair-interactions promote
mating aggregates in zooplankton and are sufficient to overcome the diffusive
nature of their mate searching behavior. Our results are useful for
understanding small-scale clustering of zooplankton, which is crucial for
predicting encounter rates and reproduction rates in the ocean.


  The COVID-19 pandemic created enormous public health and socioeconomic
challenges. The health effects of vaccination and non-pharmaceutical
interventions (NPIs) were often contrasted with significant social and economic
costs. We describe a general framework aimed to derive adaptive cost-effective
interventions, adequate for both recent and emerging pandemic threats. We also
quantify the net health benefits and propose a reinforcement learning approach
to optimise adaptive NPIs. The approach utilises an agent-based model
simulating pandemic responses in Australia, and accounts for a heterogeneous
population with variable levels of compliance fluctuating over time and across
individuals. Our analysis shows that a significant net health benefit may be
attained by adaptive NPIs formed by partial social distancing measures, coupled
with moderate levels of the society's willingness to pay for health gains
(health losses averted). We demonstrate that a socially acceptable balance
between health effects and incurred economic costs is achievable over a long
term, despite possible early setbacks.


  In this work we present a general framework for the modeling of the
transmission dynamics of macroparasites which do not reproduce within the host
like Ascaris lumbricoides, Trichuris trichiura, Necator americanus y
Ancylostoma duodenale. The basic models are derived from general probabilistic
models for the parasite density-dependent mating probability. Here we
considered the particular, and common case, of a negative binomial distribution
for the number of parasites in hosts. We find the basic reproductive number and
we show that the system exhibit a saddle-node bifurcation at some value of the
basic reproduction number. We also found the equilibria and basic reproduction
number of a model for the more general case of heteregeneous host populations.


  The mechanisms by which a protein's 3D structure can be determined based on
its amino acid sequence have long been one of the key mysteries of biophysics.
Often simplistic models, such as those derived from geometric constraints,
capture bulk real-world 3D protein-protein properties well. One approach is
using protein contact maps to better understand proteins' properties. Here, we
investigate the emergent behaviour of contact maps for different geometrically
constrained models and real-world protein systems. We derive an analytical
approximation for the distribution of model amino acid distances, $s$, by means
of a mean-field approach. This approximation is then validated for simulations
using a 2D and 3D random interaction model, as well as from contact maps of
real-world protein data. Using data from the RCSB Protein Data Bank (PDB) and
AlphaFold~2 database, the analytical approximation is fitted to protein chain
lengths of $L\approx100$, $L\approx200$, and $L\approx300$. While a universal
scaling behaviour for protein chains of different lengths could not be deduced,
we present evidence that the amino acid distance distributions can be
attributed to geometric constraints of protein chains in bulk and amino acid
sequences only play a secondary role.


  The Traill island model of Gilg et al. (2003) is a landmark attempt at
mechanistic modelling of the cyclic population dynamics of rodents, focusing on
a high Arctic community. It models the dynamics of one prey, the collared
lemming, and four predators : the stoat, the Arctic fox, the long-tailed skua
and the snowy owl. In the present short note, we first summarize how the model
works in light of theory on seasonally forced predator-prey systems, with a
focus on the temporal dynamics of predation rates. We show notably how the
impact of generalist predation, which is able here to initiate population
declines, differs slightly from that of generalist predation in other
mechanistic models of rodent-mustelid interactions such as Turchin & Hanski
(1997). We then provide a low-dimensional approximation with a single
generalist predator compartment that mimics the essential features of the
Traill island model: cycle periodicity, amplitude, shape, as well as
generalist-induced declines. This simpler model should be broadly applicable to
model other lemming populations that predominantly grow under the snow during
the winter period. Matlab computer codes for Gilg et al. (2003), its
two-dimensional approximation, as well as alternative lemming population
dynamics models are provided.


  We consider the dynamic infection spread model that is based on the discrete
SIR model which assumes infections to be spread over time via infected and
non-isolated individuals. In our system, the main objective is not to minimize
the number of required tests to identify every infection, but instead, to
utilize the available, given testing capacity $T$ at each time instance to
efficiently control the infection spread. We introduce and study a novel
performance metric, which we coin as $\epsilon$-disease control time. This
metric can be used to measure how fast a given algorithm can control the spread
of a disease. We characterize the performance of dynamic individual testing
algorithm and introduce a novel dynamic SAFFRON based group testing algorithm.
We present theoretical results and implement the proposed algorithms to compare
their performances.


  Directed evolution is a versatile technique in protein engineering that
mimics the process of natural selection by iteratively alternating between
mutagenesis and screening in order to search for sequences that optimize a
given property of interest, such as catalytic activity and binding affinity to
a specified target. However, the space of possible proteins is too large to
search exhaustively in the laboratory, and functional proteins are scarce in
the vast sequence space. Machine learning (ML) approaches can accelerate
directed evolution by learning to map protein sequences to functions without
building a detailed model of the underlying physics, chemistry and biological
pathways. Despite the great potentials held by these ML methods, they encounter
severe challenges in identifying the most suitable sequences for a targeted
function. These failures can be attributed to the common practice of adopting a
high-dimensional feature representation for protein sequences and inefficient
search methods. To address these issues, we propose an efficient, experimental
design-oriented closed-loop optimization framework for protein directed
evolution, termed ODBO, which employs a combination of novel low-dimensional
protein encoding strategy and Bayesian optimization enhanced with search space
prescreening via outlier detection. We further design an initial sample
selection strategy to minimize the number of experimental samples for training
ML models. We conduct and report four protein directed evolution experiments
that substantiate the capability of the proposed framework for finding of the
variants with properties of interest. We expect the ODBO framework to greatly
reduce the experimental cost and time cost of directed evolution, and can be
further generalized as a powerful tool for adaptive experimental design in a
broader context.


  Movements in imagery captivate the human eye and imagination. They are also
of interest in variety of scientific disciplines that study spatiotemporal
dynamics. Popular methods for quantifying movement in imagery include particle
image velocimetry and digital image correlation. Both methods are widely
applied in engineering and materials science, but less applied in other
disciplines. This paper describes an implementation of a basic digital image
correlation algorithm in R open source software as well as an extension
designed to quantify persistent movement velocities in sequences of three or
more images. Algorithms are applied in the novel arena of landscape ecology to
quantify population movement and to produce vector fields for easy
visualization of complex movement patterns across space. Functions to
facilitate analyses are available in the ICvectorfields software package. These
methods and functions are likely to produce novel insights in theoretical and
landscape ecology because they facilitate visualization and comparison of
theoretical and observed data in complex and heterogeneous environments.


  We present two methods to reveal protein-ligand unbinding mechanisms in
biased unbinding simulations by clustering trajectories into ensembles
representing unbinding paths. The first approach is based on a contact
principal component analysis for reducing the dimensionality of the input data,
followed by identification of unbinding paths and training a machine learning
model for trajectory clustering. The second approach clusters trajectories
according to their pairwise mean Euclidean distance employing the neighbor-net
algorithm, which takes into account input data bias in the distances set and is
superior to dendrogram construction. Finally, we describe a more complex case
where the reaction coordinate relevant for path identification is a single
intra-ligand hydrogen bond, highlighting the challenges involved in unbinding
path reaction coordinate detection.


  Deep Learning (DL) is a surprisingly successful branch of machine learning.
The success of DL is usually explained by focusing analysis on a particular
recent algorithm and its traits. Instead, we propose that an explanation of the
success of DL must look at the population of all algorithms in the field and
how they have evolved over time. We argue that cultural evolution is a useful
framework to explain the success of DL. In analogy to biology, we use
`development' to mean the process converting the pseudocode or text description
of an algorithm into a fully trained model. This includes writing the
programming code, compiling and running the program, and training the model. If
all parts of the process don't align well then the resultant model will be
useless (if the code runs at all!). This is a constraint. A core component of
evolutionary developmental biology is the concept of deconstraints -- these are
modification to the developmental process that avoid complete failure by
automatically accommodating changes in other components. We suggest that many
important innovations in DL, from neural networks themselves to hyperparameter
optimization and AutoGrad, can be seen as developmental deconstraints. These
deconstraints can be very helpful to both the particular algorithm in how it
handles challenges in implementation and the overall field of DL in how easy it
is for new ideas to be generated. We highlight how our perspective can both
advance DL and lead to new insights for evolutionary biology.


  Protein complexes are macromolecules essential to the functioning and
well-being of all living organisms. As the structure of a protein complex, in
particular its region of interaction between multiple protein subunits (i.e.,
chains), has a notable influence on the biological function of the complex,
computational methods that can quickly and effectively be used to refine and
assess the quality of a protein complex's 3D structure can directly be used
within a drug discovery pipeline to accelerate the development of new
therapeutics and improve the efficacy of future vaccines. In this work, we
introduce the Equivariant Graph Refiner (EGR), a novel E(3)-equivariant graph
neural network (GNN) for multi-task structure refinement and assessment of
protein complexes. Our experiments on new, diverse protein complex datasets,
all of which we make publicly available in this work, demonstrate the
state-of-the-art effectiveness of EGR for atomistic refinement and assessment
of protein complexes and outline directions for future work in the field. In
doing so, we establish a baseline for future studies in macromolecular
refinement and structure analysis.


  A structural genetic model incorporating a modern understanding of the genome
and common practice in genome-wide association studies is derived
mathematically. The model shows the Haldane map distance as a direct
consequence of the structure of the genome. An expression for genetic
correlation is derived under the model and compared to data resulting from the
BSB mouse cross. A correlation test plot is introduced for this comparison and
shows the close agreement of the model and empirical results. Noteworthy
departures in this plot indicate regions which warrant further investigation.


  Massive molecular simulations of drug-target proteins have been used as a
tool to understand disease mechanism and develop therapeutics. This work
focuses on learning a generative neural network on a structural ensemble of a
drug-target protein, e.g. SARS-CoV-2 Spike protein, obtained from
computationally expensive molecular simulations. Model tasks involve
characterizing the distinct structural fluctuations of the protein bound to
various drug molecules, as well as efficient generation of protein
conformations that can serve as an complement of a molecular simulation engine.
Specifically, we present a geometric autoencoder framework to learn separate
latent space encodings of the intrinsic and extrinsic geometries of the protein
structure. For this purpose, the proposed Protein Geometric AutoEncoder
(ProGAE) model is trained on the protein contact map and the orientation of the
backbone bonds of the protein. Using ProGAE latent embeddings, we reconstruct
and generate the conformational ensemble of a protein at or near the
experimental resolution, while gaining better interpretability and
controllability in term of protein structure generation from the learned latent
space. Additionally, ProGAE models are transferable to a different state of the
same protein or to a new protein of different size, where only the dense layer
decoding from the latent representation needs to be retrained. Results show
that our geometric learning-based method enjoys both accuracy and efficiency
for generating complex structural variations, charting the path toward scalable
and improved approaches for analyzing and enhancing high-cost simulations of
drug-target proteins.


  Proteins interact to form complexes to carry out essential biological
functions. Computational methods have been developed to predict the structures
of protein complexes. However, an important challenge in protein complex
structure prediction is to estimate the quality of predicted protein complex
structures without any knowledge of the corresponding native structures. Such
estimations can then be used to select high-quality predicted complex
structures to facilitate biomedical research such as protein function analysis
and drug discovery. We challenge this significant task with DProQ, which
introduces a gated neighborhood-modulating Graph Transformer (GGT) designed to
predict the quality of 3D protein complex structures. Notably, we incorporate
node and edge gates within a novel Graph Transformer framework to control
information flow during graph message passing. We train and evaluate DProQ on
four newly-developed datasets that we make publicly available in this work. Our
rigorous experiments demonstrate that DProQ achieves state-of-the-art
performance in ranking protein complex structures.


  Phylogenetic tree reconciliation is employed for the examination of
coevolution between host and symbiont species. An important concern is the
requirement for dependable cost values when selecting event-based parsimonious
reconciliation. Although certain approaches deduce event probabilities unique
to each pair of host and symbiont trees, which can subsequently be converted
into cost values, a significant limitation lies in their inability to model the
invasion of diverse host species by the same symbiont species (termed as a
spread event), which is believed to occur in symbiotic relationships. Invasions
lead to the observation of multiple associations between symbionts and their
hosts (indicating that a symbiont is no longer exclusive to a single host),
which are incompatible with the existing methods of coevolution. We present
AmoCoala, an enhanced version of the tool Coala, that provides a more realistic
estimation of cophylogeny event probabilities for a given pair of host and
symbiont trees, even in the presence of spread events. We expand the classical
4-event coevolutionary model to include 2 additional spread events (vertical
and horizontal spreads) that lead to multiple associations. By incorporating
spread events, our reconciliation model enables a more accurate consideration
of multiple associations. This improvement enhances the precision of estimated
cost sets, paving the way to a more reliable reconciliation of host and
symbiont trees. Our results showcase that AmoCoala produces biologically
plausible reconciliation scenarios, further emphasizing its effectiveness. The
software is accessible at https://github.com/sinaimeri/AmoCoala


  Spatially resolved transcriptomics (ST) measures gene expression along with
the spatial coordinates of the measurements. The analysis of ST data involves
significant computation complexity. In this work, we propose gene expression
dimensionality reduction algorithm that retains spatial structure. We combine
the wavelet transformation with matrix factorization to select
spatially-varying genes. We extract a low-dimensional representation of these
genes. We consider Empirical Bayes setting, imposing regularization through the
prior distribution of factor genes. Additionally, We provide visualization of
extracted representation genes capturing the global spatial pattern. We
illustrate the performance of our methods by spatial structure recovery and
gene expression reconstruction in simulation. In real data experiments, our
method identifies spatial structure of gene factors and outperforms regular
decomposition regarding reconstruction error. We found the connection between
the fluctuation of gene patterns and wavelet technique, providing smoother
visualization. We develop the package and share the workflow generating
reproducible quantitative results and gene visualization. The package is
available at https://github.com/OliverXUZY/waveST.


  Gene regulatory network (GRN) refers to the complex network formed by
regulatory interactions between genes in living cells. In this paper, we
consider inferring GRNs in single cells based on single cell RNA sequencing
(scRNA-seq) data. In scRNA-seq, single cells are often profiled from mixed
populations and their cell identities are unknown. A common practice for single
cell GRN analysis is to first cluster the cells and infer GRNs for every
cluster separately. However, this two-step procedure ignores uncertainty in the
clustering step and thus could lead to inaccurate estimation of the networks.
To address this problem, we propose to model scRNA-seq by the mixture
multivariate Poisson log-normal (MPLN) distribution. The precision matrices of
the MPLN are the GRNs of different cell types and can be jointly estimated by
maximizing MPLN's lasso-penalized log-likelihood. We show that the MPLN model
is identifiable and the resulting penalized log-likelihood estimator is
consistent. To avoid the intractable optimization of the MPLN's log-likelihood,
we develop an algorithm called VMPLN based on the variational inference method.
Comprehensive simulation and real scRNA-seq data analyses reveal that VMPLN
performs better than the state-of-the-art single cell GRN methods.


  We study the dynamic evolution of COVID-19 cased by the Omicron variant via a
fractional susceptible-exposedinfected-removed (SEIR) model. Preliminary data
suggest that the symptoms of Omicron infection are not prominent and the
transmission is therefore more concealed, which causes a relatively slow
increase in the detected cases of the new infected at the beginning of the
pandemic. To characterize the specific dynamics, the Caputo-Hadamard fractional
derivative is adopted to refined the classical SEIR model. Based on the
reported data, we infer the fractional order, timedependent parameters, as well
as unobserved dynamics of the fractional SEIR model via fractional
physics-informed neural networks (fPINNs). Then, we make short-time predictions
using the learned fractional SEIR model.


  Protein function is inherently linked to its localization within the cell,
and fluorescent microscopy data is an indispensable resource for learning
representations of proteins. Despite major developments in molecular
representation learning, extracting functional information from biological
images remains a non-trivial computational task. Current state-of-the-art
approaches use autoencoder models to learn high-quality features by
reconstructing images. However, such methods are prone to capturing noise and
imaging artifacts. In this work, we revisit deep learning models used for
classifying major subcellular localizations, and evaluate representations
extracted from their final layers. We show that simple convolutional networks
trained on localization classification can learn protein representations that
encapsulate diverse functional information, and significantly outperform
autoencoder-based models. We also propose a robust evaluation strategy to
assess quality of protein representations across different scales of biological
function.


  Aedes albopictus has been newly recorded in Agdal district at Rabat in
Morocco. The establishment of this invasive mosquito could affect the public
health by causing serious epidemics despite of its high nuisance in urban and
sub-urban areas. Vector control is mainly based on environmental management but
chemical insecticides can be used to reduce adult mosquito densities during
peak periods. However, the level of susceptibility of this mosquito to
insecticides has not been studied yet in Morocco. This paper reports the
results of the first study conducted to monitor the insecticide resistance of
adult and larva Ae. albopictus to the insecticides currently used in the vector
control. The study was carried out during May-June 2018 at Rabat from the
north-west of the country. Adult susceptibility tests were performed following
the WHO test procedures. One organochlorate (DDT 4%), one pyrethroids
(cyfluthrin 0.15%), one carbamate (bendiocarb 0.1%) and one organophosphate
(fenithrothion 1%) were tested at diagnostic doses (DD). The five-fold DD of
bendiocarb were also used to yield information on the intensity of mosquito
adult resistance. The results of the performed susceptibility bioassay showed
that the vector is susceptible to cyfluthrin and resistant to DDT, bendiocarb
and fenithrothion. Larval bioassays to temephos were conducted according to WHO
standard operating protocol to establish the dose-mortality relationship and
deduct the LC50 and LC90 then resistance ratios. We show that larval
populations of Ae. albopictus are still sensible to this insecticide. This
information could help policy-makers to plan insecticide resistance management.


  Developing accurate and reliable methods to estimate vaccine protection is a
key goal in immunology and public health. While several statistical methods
have been proposed, their potential inaccuracy in capturing fast intra-seasonal
waning of vaccine-induced protection needs to be rigorously investigated. To
compare statistical methods for vaccine effectiveness (VE) estimation, we
generated simulated data using a multiscale agent-based model of an epidemic
with an acute viral infection and differing extents of VE waning. We extended
the previously proposed framework for VE measures based on the observational
data richness to assess changes of vaccine-induced protection with time. While
VE measures based on hard-to-collect information (e.g. exact timing of
exposures) were accurate, usually VE studies rely on time-to-infection data and
the Cox proportional hazard model. We found that its extension utilizing scaled
Schoenfeld residuals, previously proposed for capturing VE waning, was
unreliable in capturing both the degree of waning and its functional form and
identified the mathematical factors contributing to this unreliability. We
showed that partitioning time and including a time-vaccine interaction term in
the Cox model significantly improved estimation of VE waning, even in the case
of dramatic, rapid waning. We also proposed how to optimize the partitioning
scheme. Using simulated data, we compared different measures of VE for
capturing the intra-seasonal waning of vaccine-induced protection. We propose
an extension of the Cox model based on including a time-vaccine interaction
term with further optimization of partitioning time. These findings may guide
future analysis of VE waning in observational data.


  Predators and their foraging strategies often determine ecosystem structure
and function. Yet, the role of protozoan predators in microbial soil ecosystems
remains elusive despite the importance of these ecosystems to global
biogeochemical cycles. In particular, amoebae -- the most abundant soil
protozoan predators of bacteria -- remineralize soil nutrients and shape the
bacterial community. However, their foraging strategies and their role as
microbial ecosystem engineers remain unknown. Here we present a multi-scale
approach, connecting microscopic single-cell analysis and macroscopic whole
ecosystem dynamics, to expose a phylogenetically widespread foraging strategy,
in which an amoeba population spontaneously partitions between cells with fast,
polarized movement and cells with slow, unpolarized movement. Such
differentiated motion gives rise to efficient colony expansion and consumption
of the bacterial substrate. From these insights we construct a theoretical
model that predicts how disturbances to amoeba growth rate and movement disrupt
their predation efficiency. These disturbances correspond to distinct classes
of bacterial defenses, which allows us to experimentally validate our
predictions. All considered, our characterization of amoeba foraging identifies
amoeba mobility, and not amoeba growth, as the core determinant of predation
efficiency and a key target for bacterial defense systems.


  Identifying cell clusters is a critical step for single-cell transcriptomics
study. Despite the numerous clustering tools developed recently, the rapid
growth of scRNA-seq volumes prompts for a more (computationally) efficient
clustering method. Here, we introduce Secuer, a Scalable and Efficient speCtral
clUstERing algorithm for scRNA-seq data. By employing an anchor-based bipartite
graph representation algorithm, Secuer enjoys reduced runtime and memory usage
by orders of magnitude, especially for ultra-large datasets profiling over 1 or
even 10 million cells. Meanwhile, Secuer also achieves better or comparable
accuracy than competing methods in small and moderate benchmark datasets.
Furthermore, we showcase that Secuer can also serve as a building block for a
new consensus clustering method, Secuer-consensus, which again greatly improves
the runtime and scalability of state-of-the-art consensus clustering methods
while also maintaining the accuracy. Overall, Secuer is a versatile, accurate,
and scalable clustering framework suitable for small to ultra-large single-cell
clustering tasks.


  The slogan "nobody is safe until everybody is safe" is a dictum to raise
awareness that in an interconnected world, pandemics such as COVID-19, require
a global approach. Motivated by the ongoing COVID-19 pandemic, we model here
the spread of a virus in interconnected communities and explore different
vaccination scenarios, assuming that the efficacy of the vaccination wanes over
time. We start with susceptible populations and consider a
susceptible-vaccinated-infected-recovered model with unvaccinated ("Bronze"),
moderately vaccinated ("Silver") and very well vaccinated ("Gold") communities,
connected through different types of networks via a diffusive linear coupling
for local spreading. We show that when considering interactions in
"Bronze"-"Gold" and "Bronze"-"Silver" communities, the "Bronze" community is
driving an increase in infections in the "Silver" and "Gold" communities. This
shows a detrimental, unidirectional effect of non-vaccinated to vaccinated
communities. Regarding the interactions between "Gold", "Silver" and "Bronze"
communities in a network, we find that two factors play central role: the
coupling strength in the dynamics and network density. When considering the
spread of a virus in Barab\'asi-Albert networks, infections in "Silver" and
"Gold" communities are lower than in "Bronze" communities. We find that the
"Gold" communities are the best in keeping their infection levels low. However,
a small number of "Bronze" communities are enough to give rise to an increase
in infections in moderately and well-vaccinated communities. When studying the
spread of a virus in a dense Erd\H{o}s-R\'enyi, and sparse Watts-Strogatz and
Barab\'asi-Albert networks, the communities reach the disease-free state in the
dense Erd\H{o}s-R\'enyi networks, but not in the sparse Watts-Strogatz and
Barab\'asi-Albert networks. However, we also find that if all these networks...


  Atomically detailed simulations of RNA folding have proven very challenging
in view of the difficulties of developing realistic force fields and the
intrinsic computational complexity of sampling rare conformational transitions.
To tackle both these issues, we extend to RNA an enhanced path sampling method
previously successfully applied to proteins. In this scheme, the information
about the RNA's native structure is harnessed by a soft history-dependent
biasing force, which is added to the atomistic force field, thus promoting the
generation of productive folding trajectories. Here, we report on the results
of simulations in explicit solvent of RNA molecules from 20 to 47 nucleotides
long and increasing topological complexity. From a statistical analysis of the
folding pathways we infer that, differently from proteins, the underlying free
energy landscape is significantly frustrated, even for relatively small chains
with a simple topology. The folding mechanism and the thermodynamics are in
agreement with the available experiments and some of the existing
coarse-grained models. This scheme provides a fully microscopic
characterization of RNA folding, relating the kinetics and dynamics of the
transition to the chemistry of the chain and its solvent. Therefore, it
provides a transferable framework that sets the stage for future translational
applications.


  Modern molecular technologies have revolutionized our understanding of
bacterial epidemiology, but reported data across different settings remain
under-integrated in common theoretical frameworks. Pneumococcus serotype
co-colonization, caused by the polymorphic bacteria Streptococcus pneumoniae,
has been increasingly investigated in recent years. While the global genomic
diversity and serotype distribution of S. pneumoniae are well-characterized,
there is limited information on how co-colonization patterns vary globally,
critical for understanding bacterial evolution and dynamics. Gathering a rich
dataset of cross-sectional pneumococcal colonization studies in the literature,
we quantified patterns of transmission intensity and co-colonization prevalence
in children populations across 17 geographic locations. Fitting these data to
an SIS model with co-colonization under the assumption of similarity among
interacting strains, our analysis reveals strong patterns of negative
co-variation between transmission intensity ($R_0$) and susceptibility to
co-colonization ($k$). In support of the stress-gradient hypothesis in ecology
(SGH), pneumococcus serotypes appear to compete more in high-transmission
settings and less in low-transmission settings, a trade-off which ultimately
leads to a conserved ratio of single to co-colonization $\mu=1/(R_0-1)k$.
Within our mathematical model, such conservation suggests preservation of
'stability-diversity-complexity' regimes in multi-strain coexistence. We find
no major study differences in serotype composition, pointing to underlying
adaptation of the same set of serotypes across environments. Our work
highlights that understanding pneumococcus transmission patterns from global
epidemiological data can benefit from simple analytical approaches that account
for quasi-neutrality among strains, co-colonization, as well as variable
environmental adaptation.


  Stereotypes are generalized beliefs about groups of people, which are used to
make decisions and judgments about them. Although such heuristics can be useful
when decisions must be made quickly, or when information is lacking, they can
also serve as the basis for prejudice and discrimination. In this paper we
study the evolution of stereotypes through group reciprocity. We characterize
the warmth of a stereotype as the willingness to cooperate with an individual
based solely on the identity of the group they belong to. We show that when
stereotypes are coarse, such group reciprocity is less likely to evolve, and
stereotypes tend to be negative. We also show that, even when stereotypes are
broadly positive, individuals are often overly pessimistic about the
willingness of those they stereotype to cooperate. We then show that the
tendency for stereotyping itself to evolve is driven by the costs of cognition,
so that more people are stereotyped with greater coarseness as costs increase.
Finally we show that extrinsic "shocks", in which the benefits of cooperation
are suddenly reduced, can cause stereotype warmth and judgement bias to turn
sharply negative, consistent with the view that economic and other crises are
drivers of out-group animosity.


  Mathematical modelling of epidemiological and coevolutionary dynamics is
widely being used to improve disease management strategies of infectious
diseases. Many diseases present some form of intra-host quiescent stage, also
known as covert infection, while others exhibit dormant stages in the
environment. As quiescent/dormant stages can be resistant to drug, antibiotics,
fungicide treatments, it is of practical relevance to study the influence of
these two life-history traits on the coevolutionary dynamics. We develop first
a deterministic coevolutionary model with two parasite types infecting one host
type and study analytically the stability of the dynamical system. We
specifically derive a stability condition for a five-by-five system of
equations with quiescence. Second, we develop a stochastic version of the model
to study the influence of quiescence on stochasticity of the system dynamics.
We compute the steady state distribution of the parasite types which follows a
multivariate normal distribution. Furthermore, we obtain numerical solutions
for the covariance matrix of the system under symmetric and asymmetric
quiescence rates between parasite types. When parasite strains are identical,
quiescence increases the variance of the number of infected individuals at high
transmission rate and vice versa when the transmission rate is low. However,
when there is competition between parasite strains with different quiescent
rates, quiescence generates a moving average behaviour which dampen off
stochasticity and decreases the variance of the number of infected hosts. The
strain with the highest rate of entering quiescence determines the strength of
the moving average and the magnitude of reduction of stochasticity. Thus, it is
worth investigating simple models of multi-strain parasite under
quiescence/dormancy to improve disease management strategies.


  In this paper, we present a mathematical model for the interaction between
honey bees and mites. The dynamics of a mite-infested honey bee colony and the
evaluation of the commonly used mite-control strategies (traditional,
mechanical and chemical) are studied. The mite-free and mite reproduction
numbers are derived using the next generation operator approach. The
mathematical analysis of the model reveals that in the absence of mites, the
colony survives if the mite-free reproduction number is greater than unity
otherwise it goes extinct. Stability and sensitivity analyses of the model
reveal that the egg laying rate of the queen bee is key in regulating mite
reproduction. Adult bee grooming and hygienic behavior of worker bees have also
been shown to play a vital role in reducing parasitism. Using the
Volterra-Lyapunov stable matrix approach, the mite-infested equilibrium is
confirmed to be globally asymptotically stable when the mite-free reproduction
number is greater than unity and the mite reproduction number is greater than
unity. It is also shown that varroa mite control strategies that focus on
limiting mite-reproduction such as caging the queen bee and using a young queen
bee quickly reduces the mite reproduction number to a value less than unity
compared to those that are intended to kill the mites.


  Most evolutionary-oriented deep generative models do not explicitly consider
the underlying evolutionary dynamics of biological sequences as it is performed
within the Bayesian phylogenetic inference framework. In this study, we propose
a method for a deep variational Bayesian generative model (EvoVGM) that jointly
approximates the true posterior of local evolutionary parameters and generates
sequence alignments. Moreover, it is instantiated and tuned for continuous-time
Markov chain substitution models such as JC69, K80 and GTR. We train the model
via a low-variance stochastic estimator and a gradient ascent algorithm. Here,
we analyze the consistency and effectiveness of EvoVGM on synthetic sequence
alignments simulated with several evolutionary scenarios and different sizes.
Finally, we highlight the robustness of a fine-tuned EvoVGM model using a
sequence alignment of gene S of coronaviruses.


  Identification of human genes involved in the aging process is critical due
to the incidence of many diseases with age. A state-of-the-art approach for
this purpose infers a weighted dynamic aging-specific subnetwork by mapping
gene expression (GE) levels at different ages onto the protein-protein
interaction network (PPIN). Then, it analyzes this subnetwork in a supervised
manner by training a predictive model to learn how network topologies of known
aging- vs. non-aging-related genes change across ages. Finally, it uses the
trained model to predict novel aging-related genes. However, the best current
subnetwork resulting from this approach still yields suboptimal prediction
accuracy. This could be because it was inferred using outdated GE and PPIN
data. Here, we evaluate whether analyzing a weighted dynamic aging-specific
subnetwork inferred from newer GE and PPIN data improves prediction accuracy
upon analyzing the best current subnetwork inferred from outdated data.
Unexpectedly, we find that not to be the case. To understand this, we perform
aging-related pathway and Gene Ontology (GO) term enrichment analyses. We find
that the suboptimal prediction accuracy, regardless of which GE or PPIN data is
used, may be caused by the current knowledge about which genes are
aging-related being incomplete, or by the current methods for inferring or
analyzing an aging-specific subnetwork being unable to capture all of the
aging-related knowledge. These findings can potentially guide future directions
towards improving supervised prediction of aging-related genes via -omics data
integration.


  In multicellular systems, the single-cell behaviors should be coordinated
consistently with the overall population dynamics and functions. However, the
interrelation between single-cell rules and the population-level goal is still
elusive. In this work, we reveal that these two levels are naturally connected
via a gradient flow structure of the heterogeneous cellular population and that
biologically prevalent single-cell rules such as unidirectional type-switching
and hierarchical order in types emerge from this structure. We also demonstrate
the gradient flow structure in a standard model of the T-cell immune response.
This theoretical framework works as a basis for understanding multicellular
dynamics and functions.


  The extent of parallel evolution at the genotypic level is quantitatively
linked to the distribution of beneficial fitness effects (DBFE) of mutations.
The standard view, based on light-tailed distributions (i.e. distributions with
finite moments), is that the probability of parallel evolution in duplicate
populations is inversely proportional to the number of available mutations, and
moreover that the DBFE is sufficient to determine the probability when the
number of available mutations is large. Here we show that when the DBFE is
heavy-tailed, as found in several recent experiments, these expectations are
defied. The probability of parallel evolution decays anomalously slowly in the
number of mutations or even becomes independent of it, implying higher
repeatability of evolution. At the same time, the probability of parallel
evolution is non-self-averaging, that is, it does not converge to its mean
value even when a large number of mutations are involved. This behavior arises
because the evolutionary process is dominated by only a few mutations of high
weight. Consequently, the probability varies widely across systems with the
same DBFE. Contrary to the standard view, the DBFE is no longer sufficient to
determine the extent of parallel evolution, making it much less predictable. We
illustrate these ideas theoretically and through analysis of empirical data on
antibiotic resistance evolution.


  Elastic strain energy is stored and released from long, distal tendons such
as the Achilles during locomotion, reducing locomotor energy cost by minimising
muscle shortening distance and speed, and thus activation. However, numerous
additional, often unrecognised, advantages of long tendons may speculatively be
of greater evolutionary advantage, including the reduced limb inertia afforded
by shorter and lighter muscles (reducing proximal muscle force requirement);
reduced energy dissipation during the foot-ground collision; capacity to store
and reuse the muscle work done to dampen the vibrations triggered by
foot-ground collisions; and attenuation of work-induced muscle damage.
Cumulatively, these effects should reduce both neuromotor fatigue and sense of
locomotor effort, allowing humans to choose to move at faster speeds for
longer. As these benefits are greater at faster locomotor speeds, they are
consistent with the hypothesis that running gaits used by our ancestors exerted
substantial evolutionary pressure on Achilles tendon length.


  Predicted inter-chain residue-residue contacts can be used to build the
quaternary structure of protein complexes from scratch. However, only a small
number of methods have been developed to reconstruct protein quaternary
structures using predicted inter-chain contacts. Here, we present an
agent-based self-learning method based on deep reinforcement learning
(DRLComplex) to build protein complex structures using inter-chain contacts as
distance constraints. We rigorously tested DRLComplex on two standard datasets
of homodimeric and heterodimeric protein complexes (i.e., the CASP-CAPRI
homodimer and Std_32 heterodimer datasets) using both true and predicted
interchain contacts as inputs. Utilizing true contacts as input, DRLComplex
achieved high average TM-scores of 0.9895 and 0.9881 and a low average
interface RMSD (I_RMSD) of 0.2197 and 0.92 on the two datasets, respectively.
When predicted contacts are used, the method achieves TM-scores of 0.73 and
0.76 for homodimers and heterodimers, respectively. Our experiments find that
the accuracy of reconstructed quaternary structures depends on the accuracy of
the contact predictions. Compared to other optimization methods for
reconstructing quaternary structures from inter-chain contacts, DRLComplex
performs similar to an advanced gradient descent method and better than a
Markov Chain Monte Carlo simulation method and a simulated annealing-based
method, validating the effectiveness of DRLComplex for quaternary
reconstruction of protein complexes.


  Consider a branching process with a homogeneous reproduction law. Sampling a
single cell uniformly from the population at a time $T > 0$ and looking along
the sampled cell's ancestral lineage, we find that the reproduction law is
heterogeneous - the expected reproductive output of ancestral cells on the
lineage from time $0$ to time $T$ continuously increases. This `inspection
paradox' is due to sampling bias, that cells with a larger number of offspring
are more likely to have one of their descendants sampled by virtue of their
prolificity, and the bias's strength grows with the random population size
and/or the sampling time $T$. Our main result explicitly characterises the
evolution of reproduction rates and sizes along the sampled ancestral lineage
as a mixture of Poisson processes, which simplifies in special cases. The
ancestral bias helps to explain recently observed variation in mutation rates
along lineages of the developing human embryo.


  Sequence diversity is one of the major challenges in the design of
diagnostic, prophylactic and therapeutic interventions against viruses. Herein,
we present DiMA, a tool designed to facilitate the dissection of sequence
diversity dynamics for viruses. As a base, DiMA provides a quantitative
overview of sequence diversity by use of Shannon's entropy, applied via a
user-defined k-mer sliding window to an input alignment file. Distinctively,
the key feature is that DiMA interrogates diversity dynamics by dissecting each
k-mer position to various diversity motifs, defined based on the incidence of
distinct sequences. At a given position, an index is a predominant sequence,
while all the others are (total) variants to the index, sub-classified into the
major (most common) variant, minor variants (occurring more than once and of
frequency lower than the major), and the unique (singleton) variants. Moreover,
DiMA allows for metadata enrichment of the motifs. DiMA is big data ready and
provides an interactive output, depicting multiple facets of sequence
diversity, with download options. It enables comparative genome/proteome
diversity dynamics analyses, within and between sequences of viral species. The
web server is publicly available at https://dima.bezmialem.edu.tr.


  We investigate the dynamics of music creation style distributions to
understand cultural evolution involving advanced intelligence. Using
statistical modelling methods and several musical statistics extracted from
datasets of popular music created in Japan and the United States (the US), we
explored the dynamics of cluster structures and constructed a fitness-based
evolutionary model to analyze and predict the evolution of music creation style
distributions. We found that intra-cluster dynamics, such as the contraction of
a cluster and the shift of a cluster centre, as well as inter-cluster dynamics
represented by clusters' relative frequencies, often exhibit notable dynamical
modes that hold across the cultures and different musical aspects.
Additionally, we found that the evolutionary model incorporating these
dynamical modes is effective for predicting the future creation style
distributions, with the predictions of cluster frequencies and cluster
variances often having comparable contributions. Our results highlight the
relevance of intra-cluster dynamics in cultural evolution, which have often
been overlooked in previous studies.


  Our world is ambiguous and this is reflected in the data we use to train our
algorithms. This is particularly true when we try to model natural processes
where collected data is affected by noisy measurements and differences in
measurement techniques. Sometimes, the process itself is ambiguous, such as in
the case of RNA folding, where the same nucleotide sequence can fold into
different structures. This suggests that a predictive model should have similar
probabilistic characteristics to match the data it models. Therefore, we
propose a hierarchical latent distribution to enhance one of the most
successful deep learning models, the Transformer, to accommodate ambiguities
and data distributions. We show the benefits of our approach (1) on a synthetic
task that captures the ability to learn a hidden data distribution, (2) with
state-of-the-art results in RNA folding that reveal advantages on highly
ambiguous data, and (3) demonstrating its generative capabilities on
property-based molecule design by implicitly learning the underlying
distributions and outperforming existing work.


  The terrestrial water cycle links the soil and atmosphere moisture reservoirs
through four fluxes: precipitation, evaporation, runoff, and atmospheric
moisture convergence (net import of water vapor to balance runoff). Each of
these processes is essential for human and ecosystem well-being. Predicting how
the water cycle responds to changes in vegetation cover remains a challenge.
Recently, changes in plant transpiration across the Amazon basin were shown to
be associated disproportionately with changes in rainfall, suggesting that even
small declines in transpiration (e.g., from deforestation) would lead to much
larger declines in rainfall. Here, constraining these results by the law of
mass conservation, we show that in a sufficiently wet atmosphere, forest
transpiration can control atmospheric moisture convergence such that increased
transpiration enhances atmospheric moisture import and resulting water yield.
Conversely, in a sufficiently dry atmosphere increased transpiration reduces
atmospheric moisture convergence and water yield. This previously unrecognized
dichotomy can explain the otherwise mixed observations of how water yield
responds to re-greening, as we illustrate with examples from China's Loess
Plateau. Our analysis indicates that any additional precipitation recycling by
additional vegetation increases precipitation but decreases local water yield
and steady-state runoff. Therefore, in the drier regions and early stages of
ecological restoration, the role of vegetation can be confined to precipitation
recycling, while once a wetter stage is achieved, additional vegetation
enhances atmospheric moisture convergence and water yield. Evaluating the
transition between regimes, and recognizing the potential of vegetation for
enhancing moisture convergence, are crucial for characterizing the consequences
of deforestation and for motivating and guiding ecological restoration.


  The COVID-19 pandemic has been a great catastrophe that upended human lives
and caused millions of deaths all over the world. The rapid spread of the
virus, with its early-stage exponential growth and subsequent 'waves', caught
many medical professionals and decision-makers unprepared. Even though
epidemiological models have been known for almost a century (since the 'Spanish
Influenza' pandemic of 1918-20), the real-life spread of the SARS-CoV-2 virus
often confounded the modelers. While the general framework of epidemiological
models like SEIR (susceptible-exposed-infected-recovered) or SIR
(susceptible-exposed-infected) was not in question, the behavior of model
parameters turned out to be unpredictable and complicated. In particular, while
the 'basic' reproduction number, R0, can be considered a constant (for the
original SARS-CoV-2 virus, prior to the emergence of variants, R0 is between
2.5 and 3.0), the 'effective' reproduction number, R(t), was a complex function
of time, influenced by human behavior in response to the pandemic (e.g.,
masking, lockdowns, transition to remote work, etc.) To better understand these
phenomena, we model the first year of the pandemic (between February 2020 and
February 2021) for a number of localities (fifty US states, as well as several
countries) using a simple SIR model. We show that the evolution of the pandemic
can be described quite successfully by assuming that R(t) behaves in a
'viscoelastic' manner, as a sum of two or three 'damped oscillators' with
different natural frequencies and damping coefficients. These oscillators
likely correspond to different sub-populations having different reactions to
proposed mitigation measures. The proposed approach can offer future data
modelers new ways to fit the reproduction number evolution with time (as
compared to the purely data-driven approaches most prevalent today).


  This thesis is focused in the study of the evolution of the metabolic
repertoire of Streptomyces, which are renowned as proficient producers of
bioactive Natural Products (NPs). The main goal of my work was to contribute
into the understanding of the evolutionary mechanisms behind the evolution of
NP biosynthetic pathways. Specifically, the development of a bioinformatic
method that helps into the discovery of new NP biosynthetic pathways from
actinobacterial genome sequences with emphasis on members of the genus
Streptomyces. I developed this method using a comparative and functional
genomics perspective. My studies indicate that central metabolic enzymes were
expanded in a genus-specific manner in Actinobacteria, and that they have been
extensively recruited for the biosynthesis of NPs. Based in these observations,
I developed EvoMining, a bioinformatic pipeline for the identificatoon of novel
biosynthetic pathways in microbial genomes. Using EvoMining several new NP
biosynthetic pathways have been predicted in different members of the phylum
Actinobacteria, including the model organism S. lividans 66. To test this
approach, the genome sequence of this model strain was obtained, and its
analysis led to the discovery of an unprecedented system for peptide bond
formation, as well as a biosynthetic pathway for an arsenic-containing
metabolite. Moreover, this work also led to the identification of expansions on
a conserved metabolic node in the glycolytic pathway of Streptomyces. These
expansions occurred before the radiation of Streptomyces and are concomitant
with the evolution of their capability to produce NPs. Experimental analyses
indicate that this node evolved to mediate the interplay between central an NP
metabolism.


  As a necessary process in drug development, finding a drug compound that can
selectively bind to a specific protein is highly challenging and costly.
Drug-target affinity (DTA), which represents the strength of drug-target
interaction (DTI), has played an important role in the DTI prediction task over
the past decade. Although deep learning has been applied to DTA-related
research, existing solutions ignore fundamental correlations between molecular
substructures in molecular representation learning of drug compound
molecules/protein targets. Moreover, traditional methods lack the
interpretability of the DTA prediction process. This results in missing feature
information of intermolecular interactions, thereby affecting prediction
performance. Therefore, this paper proposes a DTA prediction method with
interactive learning and an autoencoder mechanism. The proposed model enhances
the corresponding ability to capture the feature information of a single
molecular sequence by the drug/protein molecular representation learning module
and supplements the information interaction between molecular sequence pairs by
the interactive information learning module. The DTA value prediction module
fuses the drug-target pair interaction information to output the predicted
value of DTA. Additionally, this paper theoretically proves that the proposed
method maximizes evidence lower bound (ELBO) for the joint distribution of the
DTA prediction model, which enhances the consistency of the probability
distribution between the actual value and the predicted value. The experimental
results confirm mutual transformer-drug target affinity (MT-DTA) achieves
better performance than other comparative methods.


  DNA pattern matching is essential for many widely used bioinformatics
applications. Disease diagnosis is one of these applications, since analyzing
changes in DNA sequences can increase our understanding of possible genetic
diseases. The remarkable growth in the size of DNA datasets has resulted in
challenges in discovering DNA patterns efficiently in terms of run time and
power consumption. In this paper, we propose an efficient hardware and software
codesign that determines the chance of the occurrence of repeat-expansion
diseases using DNA pattern matching. The proposed design parallelizes the DNA
pattern matching task using associative memory realized with analog
content-addressable memory and implements an algorithm that returns the maximum
number of consecutive occurrences of a specific pattern within a DNA sequence.
We fully implement all the required hardware circuits with PTM 45-nm
technology, and we evaluate the proposed architecture on a practical human DNA
dataset. The results show that our design is energy-efficient and significantly
accelerates the DNA pattern matching task compared to previous approaches
described in the literature.


  A coronavirus disease 2019 (COVID-19) model with quarantine and standard
incidence rate is first developed, then a novel analysis approach for finding
the ultimate lower bound of COVID-19 infectious individuals is proposed, which
means that the COVID-19 pandemic is uniformly persistent if the control
reproduction number $\mathcal{R}_{c}>1$. This approach can be applied to other
related biomathematical models, and some existing works can be improved by
using it. In addition, the COVID-19-free equilibrium $V^0$ is locally
asymptotically stable (LAS) if $\mathcal{R}_{c}<1$ and linearly stable if
$\mathcal{R}_{c}=1$, respectively; while $V^0$ is unstable if
$\mathcal{R}_{c}>1$.


  Learning from 3D protein structures has gained wide interest in protein
modeling and structural bioinformatics. Unfortunately, the number of available
structures is orders of magnitude lower than the training data sizes commonly
used in computer vision and machine learning. Moreover, this number is reduced
even further, when only annotated protein structures can be considered, making
the training of existing models difficult and prone to over-fitting. To address
this challenge, we introduce a new representation learning framework for 3D
protein structures. Our framework uses unsupervised contrastive learning to
learn meaningful representations of protein structures, making use of proteins
from the Protein Data Bank. We show, how these representations can be used to
solve a large variety of tasks, such as protein function prediction, protein
fold classification, structural similarity prediction, and protein-ligand
binding affinity prediction. Moreover, we show how fine-tuned networks,
pre-trained with our algorithm, lead to significantly improved task
performance, achieving new state-of-the-art results in many tasks.


  With the covid-19 pandemic still ongoing and an enormous amount of test data
available, the lessons learned over the last two years need to be developed to
a point where they can provide understanding for tackling new variants and
future diseases. The SIR-model commonly used to model disease spread, predicts
exponential initial growth, which helps establish the infectiousness of a
disease in the early days of an outbreak. Unfortunately, the exponential growth
becomes muddied by spatial, finite-size, and non-equilibrium effects in
realistic systems, and robust estimates that may be used in prediction and
description are still lacking. I here establish a second quantization framework
that allows introduction of arbitrarily complicated spatial behavior, and I
show that a simplified version of this model is in good agreement with both the
growth of different covid-19 variants in Denmark and analytical results from
the theory of branched polymers. Denmark is well-suited for comparison, because
the number of tests with variant information in early December 2021 is very
high, so the spread of a single variant can be followed. I expect this model to
build bridges between the epidemic modeling and solid state communities. The
long-term goal of the particular analysis in this paper is to establish priors
that allow better early estimates for the infectiousness of a new disease.


  Recent studies emphasized the necessity to identify key (human) biological
processes and pathways targeted by the Coronaviridae family of viruses,
especially SARS-CoV-2. COVID-19 caused up to 33-55\% death rates in COVID-19
patients with malignant neoplasms and Alzheimer's disease. Given this scenario,
we identified biological processes and pathways which are most likely affected
by COVID-19. The associations between various diseases and human genes known to
interact with viruses from Coronaviridae family were obtained from the IntAct
COVID-19 data set annotated with DisGeNET data. We constructed the disease-gene
network to identify genes that are involved in various comorbid diseased
states. Communities from the disease-gene network through Louvain method were
identified and functional enrichment through over-representation analysis
methodology was used to discover significant biological processes and pathways
shared between COVID-19 and other diseases. The IntAct COVID-19 data set
comprised of 828 human genes and 10,473 diseases that together constituted
nodes in the disease-gene network. Each of the 70,210 edges connects a human
gene with an associated disease. The top 10 genes linked to most number of
diseases were VEGFA, BCL2, CTNNB1, ALB, COX2, AGT, HLA-A, HMOX1, FGT2 and COMT.
The most vulnerable group of patients thus discovered had comorbid conditions
such as carcinomas, malignant neoplasms and Alzheimer's disease. Finally, we
identified 37 potentially useful biological processes and pathways for improved
therapies.


  The emergence of eusocial species is both very rare in evolutionary history
and results in remarkably successful species. By inverting an agent based
model, agent rules are discovered that display behaviors characteristic of
eusocial species as well as other behaviors that lead to unexpected population
dynamics. By holding the agents' genome constant across the colony and allowing
the agents' rules to evolve, the individual behaviors exhibit phenotypic
plasticity in response to environmental cues. The phenotypically driven
reduction of intrinsic growth rates and the emergence of non-reproducing
phenotypes both demonstrate selection pressure at the colony (system) level.
The efficiency of an evolved colony is shown to have a strong relationship to
the computational capacity of the agents. Various other emergent behaviors,
both eusocial and otherwise novel, are identified and discussed. A path forward
to more capable eusocial populations and inter-colony evolution is outlined.


  Many important problems involving molecular property prediction from 3D
structures have limited data, posing a generalization challenge for neural
networks. In this paper, we describe a pre-training technique based on
denoising that achieves a new state-of-the-art in molecular property prediction
by utilizing large datasets of 3D molecular structures at equilibrium to learn
meaningful representations for downstream tasks. Relying on the well-known link
between denoising autoencoders and score-matching, we show that the denoising
objective corresponds to learning a molecular force field -- arising from
approximating the Boltzmann distribution with a mixture of Gaussians --
directly from equilibrium structures. Our experiments demonstrate that using
this pre-training objective significantly improves performance on multiple
benchmarks, achieving a new state-of-the-art on the majority of targets in the
widely used QM9 dataset. Our analysis then provides practical insights into the
effects of different factors -- dataset sizes, model size and architecture, and
the choice of upstream and downstream datasets -- on pre-training.


  The computational drug repositioning aims to discover new uses for marketed
drugs, which can accelerate the drug development process and play an important
role in the existing drug discovery system. However, the number of validated
drug-disease associations is scarce compared to the number of drugs and
diseases in the real world. Too few labeled samples will make the
classification model unable to learn effective latent factors of drugs,
resulting in poor generalization performance. In this work, we propose a
multi-task self-supervised learning framework for computational drug
repositioning. The framework tackles label sparsity by learning a better drug
representation. Specifically, we take the drug-disease association prediction
problem as the main task, and the auxiliary task is to use data augmentation
strategies and contrast learning to mine the internal relationships of the
original drug features, so as to automatically learn a better drug
representation without supervised labels. And through joint training, it is
ensured that the auxiliary task can improve the prediction accuracy of the main
task. More precisely, the auxiliary task improves drug representation and
serving as additional regularization to improve generalization. Furthermore, we
design a multi-input decoding network to improve the reconstruction ability of
the autoencoder model. We evaluate our model using three real-world datasets.
The experimental results demonstrate the effectiveness of the multi-task
self-supervised learning framework, and its predictive ability is superior to
the state-of-the-art model.


  Deep-learning models based on whole-slide digital pathology images (WSIs)
become increasingly popular for predicting molecular biomarkers. Instance-based
models has been the mainstream strategy for predicting genetic alterations
using WSIs although bag-based models along with self-attention mechanism-based
algorithms have been proposed for other digital pathology applications. In this
paper, we proposed a novel Attention-based Multiple Instance Mutation Learning
(AMIML) model for predicting gene mutations. AMIML was comprised of successive
1-D convolutional layers, a decoder, and a residual weight connection to
facilitate further integration of a lightweight attention mechanism to detect
the most predictive image patches. Using data for 24 clinically relevant genes
from four cancer cohorts in The Cancer Genome Atlas (TCGA) studies (UCEC, BRCA,
GBM and KIRC), we compared AMIML with one popular instance-based model and four
recently published bag-based models (e.g., CHOWDER, HE2RNA, etc.). AMIML
demonstrated excellent robustness, not only outperforming all the five baseline
algorithms in the vast majority of the tested genes (17 out of 24), but also
providing near-best-performance for the other seven genes. Conversely, the
performance of the baseline published algorithms varied across different
cancers/genes. In addition, compared to the published models for genetic
alterations, AMIML provided a significant improvement for predicting a wide
range of genes (e.g., KMT2C, TP53, and SETD2 for KIRC; ERBB2, BRCA1, and BRCA2
for BRCA; JAK1, POLE, and MTOR for UCEC) as well as produced outstanding
predictive models for other clinically relevant gene mutations, which have not
been reported in the current literature. Furthermore, with the flexible and
interpretable attention-based MIL pooling mechanism, AMIML could further
zero-in and detect predictive image patches.


  In this paper, by resorting to classical methods of statistical mechanics, we
build a kinetic model able to reproduce the observed statistical weight
distribution of many diverse species. The kinetic description of the time
variations of the weight distribution is based on elementary interactions that
describe in a qualitative and quantitative way successive evolutionary updates,
and determine explicit equilibrium distributions. Numerical fittings on
mammalian eutherians of the order Chiroptera population illustrates the
effectiveness of the approach.


  In this work we analyze some models used to explain the origins of
intra-strand parity and strand compositional asymmetries in bacterial genomes.
Due to the particular way that these two features emerge in bacterial DNA, we
performed our analysis from the perspective that they are complementary
phenomena that should be addressed together. Although most of the models for
these features try to explain them as consequence of evolutionary mechanisms,
recently it was proposed that they could be `relics' of some primordial genome
that were conserved thorough out the genome evolution. We shall pay special
attention to the S-H model, which is, up to the date, the unique model proposed
as a possible explanation for intra-strand parity and strand compositional
asymmetries in primordial genomes as mere consequence of randomness under
chemical/physical constraints. In particular, we shall discuss possible
directions to test some of the hypotheses of the S-H model, and we will present
a possible formulation of the S-H model as an evolutive model too.


  Molecular surface representations have been advertised as a great tool to
study protein structure and functions, including protein-ligand binding
affinity modeling. However, the conventional surface-area-based methods fail to
deliver a competitive performance on the energy scoring tasks. The main reason
is the lack of crucial physical and chemical interactions encoded in the
molecular surface generations. We present novel molecular surface
representations embedded in different scales of the element interactive
manifolds featuring the dramatically dimensional reduction and accurately
physical and biological properties encoders. Those low-dimensional
surface-based descriptors are ready to be paired with any advanced machine
learning algorithms to explore the essential structure-activity relationships
that give rise to the element interactive surface area-based scoring functions
(EISA-score). The newly developed EISA-score has outperformed many
state-of-the-art models, including various well-established surface-related
representations, in standard PDBbind benchmarks.


  A quest to determine the complete sequence of a human DNA from telomere to
telomere started three decades ago and was finally completed in 2021. This
accomplishment was a result of a tremendous effort of numerous experts who
engineered various tools and performed laborious manual inspection to achieve
the first gapless genome sequence. However, such method can hardly be used as a
general approach to assemble different genomes, especially when the assembly
speed is critical given the large amount of data. In this work, we explore a
different approach to the central part of the genome assembly task that
consists of untangling a large assembly graph from which a genomic sequence
needs to be reconstructed. Our main motivation is to reduce human-engineered
heuristics and use deep learning to develop more generalizable reconstruction
techniques. Precisely, we introduce a new learning framework to train a graph
convolutional network to resolve assembly graphs by finding a correct path
through them. The training is supervised with a dataset generated from the
resolved CHM13 human sequence and tested on assembly graphs built using real
human PacBio HiFi reads. Experimental results show that a model, trained on
simulated graphs generated solely from a single chromosome, is able to
remarkably resolve all other chromosomes. Moreover, the model outperforms
hand-crafted heuristics from a state-of-the-art \textit{de novo} assembler on
the same graphs. Reconstructed chromosomes with graph networks are more
accurate on nucleotide level, report lower number of contigs, higher genome
reconstructed fraction and NG50/NGA50 assessment metrics.


  We derive the Hessian geometric structure of nonequilibrium chemical reaction
networks (CRN) on the flux and force spaces induced by the Legendre duality of
convex dissipation functions and characterize their dynamics as a generalized
flow. With this structure, we can extend theories of nonequilibrium systems
with quadratic dissipation functions to more general ones with nonquadratic
ones, which are pivotal for studying chemical reaction networks. By applying
generalized notions of orthogonality in Hessian geometry to chemical reaction
networks, we obtain two generalized decompositions of the entropy production
rate, each of which captures gradient-flow and minimum-dissipation aspects in
nonequilibrium dynamics.


  Background: Recent epidemic outbreaks such as the SARS-CoV-2 pandemic and the
mpox outbreak in 2022 have demonstrated the value of genomic sequencing data
for tracking the origin and spread of pathogens. Laboratories around the globe
generated new sequences at unprecedented speed and volume and bioinformaticians
developed new tools and dashboards to analyze this wealth of data. However, a
major challenge that remains is the lack of simple and efficient approaches for
accessing and processing sequencing data.
  Results: The Lightweight API for Sequences (LAPIS) facilitates rapid
retrieval and analysis of genomic sequencing data through a REST API. It
supports complex mutation- and metadata-based queries and can perform
aggregation operations on massive datasets. LAPIS is optimized for typical
questions relevant to genomic epidemiology. Using a newly-developed in-memory
database engine, it has a high speed and throughput: between 25 January and 4
February 2023, the SARS-CoV-2 instance of LAPIS, which contains 14.5 million
sequences, processed over 20 million requests with a mean response time of 411
ms and a median response time of 1 ms. LAPIS is the core engine behind our
dashboards on genspectrum.org and we currently maintain public LAPIS instances
for SARS-CoV-2 and mpox.
  Conclusions: Powered by an optimized database engine and available through a
web API, LAPIS enhances the accessibility of genomic sequencing data. It is
designed to serve as a common backend for dashboards and analyses with the
potential to be integrated into common database platforms such as GenBank.


  Circadian rhythms are biological rhythms of approximately 24 h that persist
even under constant conditions without environmental daily cues. The molecular
circadian clock machinery generates the physiological rhythms, which can be
transmitted into the downstream output system. Owing to the stochastic nature
of the biochemical reactions, the oscillation period of circadian rhythms
exhibited by individual organisms or cells is not constant on a daily basis
with variations as high as 10 % in terms of the coefficient of variation.
Although the fluctuations in circadian rhythm is measured through a reporter
such as bioluminescence or fluorescence, experimentally confirming whether the
fluctuations found in the reporter system are the same as those in the clock
itself is challenging. This study numerically and analytically investigated a
coupled system of a circadian clock and its output system, and then compared
the fluctuations in the oscillation period of the two systems. We found that
the amount of fluctuation in the output system is smaller than that in the
circadian clock when the degradation rate of the molecules responsible for the
output system was at typical values. The results obtained imply that the output
system improves the accuracy of the circadian rhythm without the need for any
special denoising processes.


  The concept of entropy in statistical physics is related to the existence of
irreversible macroscopic processes. In this work, we explore a recently
introduced entropy formula for a class of stochastic processes with more than
one absorbing state that is extensively used in population genetics models. We
will consider the Moran process as a paradigm for this class, and will extend
our discussion to other models outside this class. We will also discuss the
relation between non-extensive entropies in physics and epistasis (i.e., when
the effects of different alleles are not independent) and the role of
symmetries in population genetic models.


  Molecular conformer generation is a fundamental task in computational
chemistry. Several machine learning approaches have been developed, but none
have outperformed state-of-the-art cheminformatics methods. We propose
torsional diffusion, a novel diffusion framework that operates on the space of
torsion angles via a diffusion process on the hypertorus and an
extrinsic-to-intrinsic score model. On a standard benchmark of drug-like
molecules, torsional diffusion generates superior conformer ensembles compared
to machine learning and cheminformatics methods in terms of both RMSD and
chemical properties, and is orders of magnitude faster than previous
diffusion-based models. Moreover, our model provides exact likelihoods, which
we employ to build the first generalizable Boltzmann generator. Code is
available at https://github.com/gcorso/torsional-diffusion.


  Although previous infection and vaccination provide protection against
SARS-CoV-2 infection, both reinfection and breakthrough infection are possible
events whose occurrence would increase with time after first exposure to the
antigen and with the emergence of new variants of the virus. Periodic
vaccination could counteract this decline in protection. In the present work,
our aim was to develop and explore a model of SARS-CoV-2 spread with
vaccination, reinfection and breakthrough infection. A modified deterministic
SIS (Susceptible-Infected-Susceptible) model represented by a system of
differential equations was designed. As in any SIS model, the population was
divided into susceptible and infected individuals. But in our design,
susceptible individuals were, in turn, grouped into three consecutive
categories whose susceptibility increases with time after infection or
vaccination. The model was studied by means of computer simulations, which were
analysed qualitatively. The results obtained show that the prevalence, after
oscillating between peaks and valleys, reaches a plateau phase. Moreover, as
might be expected, the magnitude of the peaks and plateaus increases as the
infection rate rises, the vaccination rate decreases and the rate of decay of
protection conferred by vaccination or previous infection increases. Therefore,
the present study suggests that, at least under certain conditions, the spread
of SARS-CoV-2, although it could experience fluctuations, would finally evolve
into an endemic form, with a more or less stable prevalence that would depend
on the levels of infection and vaccination, and on the kinetics of
post-infection and post-vaccination protection. However, it should be kept in
mind that our development is a theoretical scheme with many limitations. For
this reason, its predictions should be considered with great care.


  Radiomics has shown a capability for different types of cancers such as
glioma to predict the clinical outcome. It can have a non-invasive means of
evaluating the immunotherapy response prior to treatment. However, the use of
deep convolutional neural networks (CNNs)-based radiomics requires large
training image sets. To avoid this problem, we investigate a new imaging
features that model distribution with a Gaussian mixture model (GMM) of learned
3D CNN features. Using these deep radiomic features (DRFs), we aim to predict
the immune marker status (low versus high) and overall survival for glioma
patients. We extract the DRFs by aggregating the activation maps of a
pre-trained 3D-CNN within labeled tumor regions of MRI scans that corresponded
immune markers of 151 patients. Our experiments are performed to assess the
relationship between the proposed DRFs, three immune cell markers (Macrophage
M1, Neutrophils and T Cells Follicular Helper), and measure their association
with overall survival. Using the random forest (RF) model, DRFs was able to
predict the immune marker status with area under the ROC curve (AUC) of 78.67,
83.93 and 75.67\% for Macrophage M1, Neutrophils and T Cells Follicular Helper,
respectively. Combined the immune markers with DRFs and clinical variables,
Kaplan-Meier estimator and Log-rank test achieved the most significant
difference between predicted groups of patients (short-term versus long-term
survival) with p\,=\,4.31$\times$10$^{-7}$ compared to p\,=\,0.03 for Immune
cell markers, p\,=\,0.07 for clinical variables , and
p\,=\,1.45$\times$10$^{-5}$ for DRFs. Our findings indicate that the proposed
features (DRFs) used in RF models may significantly consider prognosticating
patients with brain tumour prior to surgery through regularly acquired imaging
data.


  Food profiling is an essential step in any food monitoring system needed to
prevent health risks and potential frauds in the food industry. Significant
improvements in sequencing technologies are pushing food profiling to become
the main computational bottleneck. State-of-the-art profilers are unfortunately
too costly for food profiling.
  Our goal is to design a food profiler that solves the main limitations of
existing profilers, namely (1) working on massive data structures and (2)
incurring considerable data movement for a real-time monitoring system. To this
end, we propose Demeter, the first platform-independent framework for food
profiling. Demeter overcomes the first limitation through the use of
hyperdimensional computing (HDC) and efficiently performs the accurate
few-species classification required in food profiling. We overcome the second
limitation by using an in-memory hardware accelerator for Demeter (named
Acc-Demeter) based on memristor devices. Acc-Demeter actualizes several
domain-specific optimizations and exploits the inherent characteristics of
memristors to improve the overall performance and energy consumption of
Acc-Demeter.
  We compare Demeter's accuracy with other industrial food profilers using
detailed software modeling. We synthesize Acc-Demeter's required hardware using
UMC's 65nm library by considering an accurate PCM model based on silicon-based
prototypes. Our evaluations demonstrate that Acc-Demeter achieves a (1)
throughput improvement of 192x and 724x and (2) memory reduction of 36x and 33x
compared to Kraken2 and MetaCache (2 state-of-the-art profilers), respectively,
on typical food-related databases. Demeter maintains an acceptable profiling
accuracy (within 2% of existing tools) and incurs a very low area overhead.


  We consider a generic class of gene circuits affected by nonlinear extrinsic
noise. To address this nonlinearity we introduce a general perturbative
methodology based on assuming timescale separation between noise and genes
dynamics, with fluctuations exhibiting a large but finite correlation time. We
apply this methodology to the case of the toggle switch, and by considering
biologically relevant log-normal fluctuations, we find that the system exhibits
noise-induced transitions. The system becomes bimodal in regions of the
parameter space where it would be deterministically monostable. We show that by
including higher order corrections our methodology allows one to obtain correct
predictions for the occurrence of transitions even for not so large correlation
time of the fluctuations, overcoming thereby limitations of previous
theoretical approaches. Interestingly we find that at intermediate noise
intensities the noise-induced transition in the toggle switch affects one of
the genes involved, but not the other one.