  Strategies for sustaining cooperation and preventing exploitation by selfish
agents in repeated games have mostly been restricted to Markovian strategies
where the response of an agent depends on the actions in the previous round.
Such strategies are characterized by lack of learning. However, learning from
accumulated evidence over time and using the evidence to dynamically update our
response is a key feature of living organisms. Bayesian inference provides a
framework for such evidence-based learning mechanisms. It is therefore
imperative to understand how strategies based on Bayesian learning fare in
repeated games with Markovian strategies. Here, we consider a scenario where
the Bayesian player uses the accumulated evidence of the opponent's actions
over several rounds to continuously update her belief about the reactive
opponent's strategy. The Bayesian player can then act on her inferred belief in
different ways. By studying repeated Prisoner's dilemma games with such
Bayesian inferential strategies, both in infinite and finite populations, we
identify the conditions under which such strategies can be evolutionarily
stable. We find that a Bayesian strategy that is less altruistic than the
inferred belief about the opponent's strategy can outperform a larger set of
reactive strategies, whereas one that is more generous than the inferred belief
is more successful when the benefit-to-cost ratio of mutual cooperation is
high. Our analysis reveals how learning the opponent's strategy through
Bayesian inference, as opposed to utility maximization, can be beneficial in
the long run, in preventing exploitation and eventual invasion by reactive
strategies.


  Small open reading frames are understudied as they have been historically
excluded from genome annotations. However, evidence for the functional
significance of small proteins in various cellular processes accumulates.
Proteins with less than 70 residues can also confer resistance to antimicrobial
compounds, including intracellularly-acting protein toxins, membrane-acting
antimicrobial peptides and various small-molecule antibiotics. Such herein
coined Small Antimicrobial Resistance Proteins (SARPs) have emerged on
evolutionary timescales or can be enriched from protein libraries using
laboratory evolution. Our review consolidates existing knowledge on SARPs and
highlights recent advancements in proteomics and genomics that reveal pervasive
translation of unannotated genetic regions into small proteins that show
features of known SARPs. The potential contribution of small proteins to
antimicrobial resistance is awaiting exploration.


  In this paper, we present a new approach for Covid-19 Pandemic spreading
simulation based on fuzzy multi agents. The agent parameters consider
distribution of the population according to age, and the index of
socio-economic fragility. Medical knowledge affirms that the COVID-19 main risk
factors are age and obesity. The worst medical situation is caused by the
combination of these two risk factors which in almost99% of cases finish in
ICU. The appearance of virus variants is another aspect parameter by our
simulation through a simplified modeling of the contagiousness. Using real data
from people from West Indies (Guadeloupe, F.W.I.), we modeled the infection
rate of the risk population, if neither vaccination nor barrier gestures are
respected. The results show that hospital capacities are exceeded, and the
number of deaths exceeds 2% of the infected population, which is close to the
reality.


  In the last decade, some algebraic tools have been successfully applied to
phylogenetic reconstruction. These tools are mainly based on the knowledge of
equations describing algebraic varieties associated to phylogenetic trees
evolving under Markov processes of molecular substitution, the so called
phylogenetic invariants. Although the theory involved allows to explicitly
obtain these equations for all equivariant models (which include some of the
most popular nucleotide substitution models), practical uses of these algebraic
tools have been restricted to the case of the general Markov model. Arguably,
one of the reasons for this restriction is that knowledge of linear
representation theory is required before making these equations explicit.
  With the aim of enlarging the practical uses of algebraic phylogenetics, in
this paper we prove that phylogenetic invariants for trees evolving under
equivariant models can be derived from phylogenetic invariants for the general
Markov model, without the need of representation theory. Our main result states
that the algebraic variety corresponding to a phylogenetic tree evolving under
an equivariant model is an irreducible component of the variety corresponding
to the same tree under the general Markov model cut with the linear space
defined by the model. We also prove that, for any equivariant model, those
phylogenetic invariants that are relevant for practical uses (e.g. tree
reconstruction) can be simply deduced from a single rank constraint on the
matrices obtained by flattening the joint distribution at the leaves of the
tree. This condition can be easily tested from singular values of the matrices
and extends our results from trees to phylogenetic networks.


  We study the response of a quantitative trait to exponential directional
selection in a finite haploid population at the genetic and the phenotypic
level. We assume an infinite sites model, in which the number of new mutations
per generation in the population follows a Poisson distribution (with mean
$\Theta$) and each mutation occurs at a new, previously monomorphic site.
Mutation effects are beneficial and drawn from a distribution. Sites are
unlinked and contribute additively to the trait. Assuming that selection is
stronger than random genetic drift, we model the initial phase of the dynamics
by a supercritical Galton-Watson process. This enables us to obtain
time-dependent results. We show that the copy-number distribution of the mutant
in generation n, conditioned on non-extinction until n, is described accurately
by the deterministic increase from an initial distribution with mean 1. This
distribution is related to the absolutely continuous part $W^+$ of the random
variable, typically denoted $W$, that characterizes the stochasticity
accumulating during the mutant's sweep. On this basis, we derive explicitly the
(approximate) time dependence of the mutant frequency distribution, of the
expected mean and variance of the trait and of the expected number of
segregating sites. Unexpectedly, we obtain highly accurate approximations for
all times, even for the quasi-stationary phase where we refine classical
results. In addition, we find that $\Theta$ is the main determinant of the
pattern of adaptation at the genetic level, i.e., whether the initial
allele-frequency dynamics are best described by sweep-like patterns at few loci
or small allele-frequency shifts at many. The selection strength determines
primarily the rate of adaptation. The accuracy of our results is tested by
comprehensive simulations in a Wright-Fisher framework.


  Peptides play a pivotal role in a wide range of biological activities through
participating in up to 40% protein-protein interactions in cellular processes.
They also demonstrate remarkable specificity and efficacy, making them
promising candidates for drug development. However, predicting peptide-protein
complexes by traditional computational approaches, such as Docking and
Molecular Dynamics simulations, still remains a challenge due to high
computational cost, flexible nature of peptides, and limited structural
information of peptide-protein complexes. In recent years, the surge of
available biological data has given rise to the development of an increasing
number of machine learning models for predicting peptide-protein interactions.
These models offer efficient solutions to address the challenges associated
with traditional computational approaches. Furthermore, they offer enhanced
accuracy, robustness, and interpretability in their predictive outcomes. This
review presents a comprehensive overview of machine learning and deep learning
models that have emerged in recent years for the prediction of peptide-protein
interactions.


  The most popular and universally predictive protein simulation models employ
all-atom molecular dynamics (MD), but they come at extreme computational cost.
The development of a universal, computationally efficient coarse-grained (CG)
model with similar prediction performance has been a long-standing challenge.
By combining recent deep learning methods with a large and diverse training set
of all-atom protein simulations, we here develop a bottom-up CG force field
with chemical transferability, which can be used for extrapolative molecular
dynamics on new sequences not used during model parametrization. We demonstrate
that the model successfully predicts folded structures, intermediates,
metastable folded and unfolded basins, and the fluctuations of intrinsically
disordered proteins while it is several orders of magnitude faster than an
all-atom model. This showcases the feasibility of a universal and
computationally efficient machine-learned CG model for proteins.


  Protein mutations can significantly influence protein solubility, which
results in altered protein functions and leads to various diseases. Despite of
tremendous effort, machine learning prediction of protein solubility changes
upon mutation remains a challenging task as indicated by the poor scores of
normalized Correct Prediction Ratio (CPR). Part of the challenge stems from the
fact that there is no three-dimensional (3D) structures for the wild-type and
mutant proteins. This work integrates persistent Laplacians and pre-trained
Transformer for the task. The Transformer, pretrained with hunderds of millions
of protein sequences, embeds wild-type and mutant sequences, while persistent
Laplacians track the topological invariant change and homotopic shape evolution
induced by mutations in 3D protein structures, which are rendered from
AlphaFold2. The resulting machine learning model was trained on an extensive
data set labeled with three solubility types. Our model outperforms all
existing predictive methods and improves the state-of-the-art up to 15%.


  Biological mixtures, such as the cellular cytoplasm, are composed of a large
number of different components. From this heterogeneity, ordered mesoscopic
structures emerge, such as liquid phases with controlled composition. These
structures compete with each other for the same components. This raises several
questions, such as what types of interactions allow the retrieval of multiple
ordered mesoscopic structures, and what are the physical limitations for the
retrieval of said structures. In this work, we develop an analytically
tractable model for liquids capable of retrieving states with target
compositions. We name this model the liquid Hopfield model in reference to
corresponding work in the theory of associative neural networks. By solving
this model, we show that non-linear repulsive interactions are necessary for
retrieval of target structures. We demonstrate that this is because liquid
mixtures at low temperatures tend to transition to phases with few components,
a phenomenon that we term localization. Taken together, our results demonstrate
a trade-off between retrieval and localization phenomena in liquid mixtures.


  Plastid genomes (plastomes) of angiosperms are of great interest among
biologists. High-throughput sequencing is making many such genomes accessible,
increasing the need for tools to perform rapid comparative analysis. This
exploratory analysis investigates whether the Pangenome Research Tool Kit
(PGR-TK) is suitable for analyzing plastomes. After determining the optimal
parameters for this tool on plastomes, we use it to compare sequences from each
of the genera - Magnolia, Solanum, Fragaria and Cotoneaster, as well as a
combined set from 20 rosid genera. PGR-TK recognizes large-scale plastome
structures, such as the inverted repeats, among combined sequences from distant
rosid families. If the plastid genomes are rotated to the same starting point,
it also correctly groups different species from the same genus together in a
generated cladogram. The visual approach of PGR-TK provides insights into
genome evolution without requiring gene annotations.


  We consider the problem of antibody sequence design given 3D structural
information. Building on previous work, we propose a fine-tuned inverse folding
model that is specifically optimised for antibody structures and outperforms
generic protein models on sequence recovery and structure robustness when
applied on antibodies, with notable improvement on the hypervariable CDR-H3
loop. We study the canonical conformations of complementarity-determining
regions and find improved encoding of these loops into known clusters. Finally,
we consider the applications of our model to drug discovery and binder design
and evaluate the quality of proposed sequences using physics-based methods.


  We propose a kinetic model for understanding the link between opinion
formation phenomena and epidemic dynamics. The recent pandemic has brought to
light that vaccine hesitancy can present different phases and temporal and
spatial variations, presumably due to the different social features of
individuals. The emergence of patterns in societal reactions permits to design
and predict the trends of a pandemic. This suggests that the problem of vaccine
hesitancy can be described in mathematical terms, by suitably coupling a
kinetic compartmental model for the spreading of an infectious disease with the
evolution of the personal opinion of individuals, in the presence of leaders.
The resulting model makes it possible to predict the collective compliance with
vaccination campaigns as the pandemic evolves and to highlight the best
strategy to set up for maximizing the vaccination coverage. We conduct
numerical investigations which confirm the ability of the model to describe
different phenomena related to the spread of an epidemic.


  Deep learning is emerging as an effective tool in drug discovery, with
potential applications in both predictive and generative models. Generative
Flow Networks (GFlowNets/GFNs) are a recently introduced method recognized for
the ability to generate diverse candidates, in particular in small molecule
generation tasks. In this work, we introduce double GFlowNets (DGFNs). Drawing
inspiration from reinforcement learning and Double Deep Q-Learning, we
introduce a target network used to sample trajectories, while updating the main
network with these sampled trajectories. Empirical results confirm that DGFNs
effectively enhance exploration in sparse reward domains and high-dimensional
state spaces, both challenging aspects of de-novo design in drug discovery.


  We present a simple mathematical model that captures the evolutionary
capabilities of a prebiotic compartment or protocell. In the model the
protocell contains an autocatalytic set whose chemical dynamics is coupled to
the growth-division dynamics of the compartment. Bistability in the dynamics of
the autocatalytic set results in a protocell that can exist with two distinct
growth rates. Stochasticity in chemical reactions plays the role of mutations
and causes transitions from one growth regime to another. We show that the
system exhibits `natural selection', where a `mutant' protocell in which the
autocatalytic set is active arises by chance in a population of inactive
protocells, and then takes over the population because of its higher growth
rate or `fitness'. The work integrates three levels of dynamics: intracellular
chemical, single protocell, and population (or ecosystem) of protocells..


  In today's world,the risk of emerging and re-emerging epidemics have
increased.The recent advancement in healthcare technology has made it possible
to predict an epidemic outbreak in a region.Early prediction of an epidemic
outbreak greatly helps the authorities to be prepared with the necessary
medications and logistics required to keep things in control. In this article,
we try to predict the epidemic outbreak (influenza, hepatitis and malaria) for
the state of New York, USA using machine and deep learning algorithms, and a
portal has been created for the same which can alert the authorities and health
care organizations of the region in case of an outbreak. The algorithm takes
historical data to predict the possible number of cases for 5 weeks into the
future. Non-clinical factors like google search trends,social media data and
weather data have also been used to predict the probability of an outbreak.


  Many crucial biological processes rely on networks of protein-protein
interactions. Predicting the effect of amino acid mutations on protein-protein
binding is vital in protein engineering and therapeutic discovery. However, the
scarcity of annotated experimental data on binding energy poses a significant
challenge for developing computational approaches, particularly deep
learning-based methods. In this work, we propose SidechainDiff, a
representation learning-based approach that leverages unlabelled experimental
protein structures. SidechainDiff utilizes a Riemannian diffusion model to
learn the generative process of side-chain conformations and can also give the
structural context representations of mutations on the protein-protein
interface. Leveraging the learned representations, we achieve state-of-the-art
performance in predicting the mutational effects on protein-protein binding.
Furthermore, SidechainDiff is the first diffusion-based generative model for
side-chains, distinguishing it from prior efforts that have predominantly
focused on generating protein backbone structures.


  With the rise of Transformers and Large Language Models (LLMs) in Chemistry
and Biology, new avenues for the design and understanding of therapeutics have
opened up to the scientific community. Protein sequences can be modeled as
language and can take advantage of recent advances in LLMs, specifically with
the abundance of our access to the protein sequence datasets. In this paper, we
developed the GPCR-BERT model for understanding the sequential design of G
Protein-Coupled Receptors (GPCRs). GPCRs are the target of over one-third of
FDA-approved pharmaceuticals. However, there is a lack of comprehensive
understanding regarding the relationship between amino acid sequence, ligand
selectivity, and conformational motifs (such as NPxxY, CWxP, E/DRY). By
utilizing the pre-trained protein model (Prot-Bert) and fine-tuning with
prediction tasks of variations in the motifs, we were able to shed light on
several relationships between residues in the binding pocket and some of the
conserved motifs. To achieve this, we took advantage of attention weights, and
hidden states of the model that are interpreted to extract the extent of
contributions of amino acids in dictating the type of masked ones. The
fine-tuned models demonstrated high accuracy in predicting hidden residues
within the motifs. In addition, the analysis of embedding was performed over 3D
structures to elucidate the higher-order interactions within the conformations
of the receptors.


  Cancer is a complex disease that is characterized by uncontrolled growth and
division of cells. It involves a complex interplay between genetic and
environmental factors that lead to the initiation and progression of tumors.
Recent advances in molecular dynamics simulations have revolutionized our
understanding of the molecular mechanisms underlying cancer initiation and
progression. Molecular dynamics simulations enable researchers to study the
behavior of biomolecules at an atomic level, providing insights into the
dynamics and interactions of proteins, nucleic acids, and other molecules
involved in cancer development. In this review paper, we provide an overview of
the latest advances in molecular dynamics simulations of cancer cells. We will
discuss the principles of molecular dynamics simulations and their applications
in cancer research. We also explore the role of molecular dynamics simulations
in understanding the interactions between cancer cells and their
microenvironment, including signaling pathways, proteinprotein interactions,
and other molecular processes involved in tumor initiation and progression. In
addition, we highlight the current challenges and opportunities in this field
and discuss the potential for developing more accurate and personalized
simulations. Overall, this review paper aims to provide a comprehensive
overview of the current state of molecular dynamics simulations in cancer
research, with a focus on the molecular mechanisms underlying cancer initiation
and progression.


  It has been recently shown that the exponential growth rate of a population
of bacterial cells starting from a single cell shows transient oscillations due
to early synchronized bursts of division. These oscillations are enhanced by
cell size regulation and contain information about single-cell growth
statistics. Here, we report a phase transition in these oscillations as a
function of growth rate variability. Below the transition point, these
oscillations become asymptotically deterministic and can be measured
experimentally, while above the transition point, the stochasticity in
population growth dominates the oscillations and masks all the information
about single cell growth statistics. The analytically calculated transition
point, which roughly corresponds to $13\%$ variability in single-cell growth
rate, falls within physiologically relevant parameters. Additionally, we show
that the oscillations can stochastically emerge even when the initial state
contains multiple cells with out-of-phase division cycles. We show that the
amplitude and the phase of these oscillations are stochastic and would vary
across repeated measurements with the same initial conditions. We provide
analytic expressions as well as numerical estimates for the typical oscillation
amplitude and the number of generations before the amplitude falls below a
given measurement threshold for E. coli in multiple growth conditions.


  Competition between parasitoids can reduce the success of pest control in
biological programs using two species as bio-control agents or when multiple
species exploit the same host crop. Parasitoid foraging behavior and the
ability to identify already parasitized hosts affect the efficacy of parasitoid
species as bio-agents to regulate pest insects. We evaluated the behavioural
changes of parasitoids according to the quality of hosts ({\it i.e.},
previously parasitised or not), and the characterisation of these transitions
over time via multi-state models. We evaluated the effects of previous
parasitism of the brown stinkbug {\it Euschistus heros} eggs on the parasitism
rate of the species {\it Trissolcus basalis} and {\it Telenomus podisi}. We
successively modelled the choice of eggs (with three possibilities: non
parasitised eggs, eggs previously parasitised by {\it T. podisi}, and eggs
previously parasitised by {\it T. basalis}) and the conditional behaviour given
the choice (walking, drumming, ovipositing or marking the chosen egg). We
consider multi-state models in two successive stages to calculate double
transition probabilities, and the statistical methodology is based on the
maximum likelihood procedure. Using the Cox model and assuming a stationary
process, we verified that the treatment effect was significant for the choice,
indicating that the two parasitoid species have different choice patterns. For
the second stage, i.e. behaviour given the choice, the results also showed the
influence of the species on the conditional behaviour, especially for
previously parasitised eggs. Specifically, {\it T.podisi} avoids intraspecific
competition and makes decisions faster than {\it T. basalis}. In this work, we
emphasise the methodological contribution with multi-state models, especially
in the context of double transitions.


  At what level does selective pressure effectively act? When considering the
reproductive dynamics of interacting and mutating agents, it has long been
debated whether selection is better understood by focusing on the individual or
if hierarchical selection emerges as a consequence of joint adaptation. Despite
longstanding efforts in theoretical ecology there is still no consensus on this
fundamental issue, most likely due to the difficulty in obtaining adequate data
spanning sufficient number of generations and the lack of adequate tools to
quantify the effect of hierarchical selection. Here we capitalise on recent
advances in information-theoretic data analysis to advance this state of
affairs by investigating the emergence of high-order structures -- such as
groups of species -- in the collective dynamics of the Tangled Nature model of
evolutionary ecology. Our results show that evolutionary dynamics can lead to
clusters of species that act as a selective group, that acquire
information-theoretic agency. Overall, our findings provide quantitative
evidence supporting the relevance of high-order structures in evolutionary
ecology, which can emerge even from relatively simple processes of adaptation
and selection.


  We consider a stochastic individual-based model of adaptive dynamics for an
asexually reproducing population with mutation, with linear birth and death
rates, as well as a density-dependent competition. To depict repeating changes
of the environment, all of these parameters vary over time as piecewise
constant and periodic functions, on an intermediate time-scale between those of
stabilization of the resident population (fast) and exponential growth of
mutants (slow). Studying the growth of emergent mutants and their invasion of
the resident population in the limit of small mutation rates for a
simultaneously diverging population size, we are able to determine their
effective growth rates. We describe this growth as a mesoscopic scaling-limit
of the orders of population sizes, where we observe an averaging effect of the
invasion fitness. Moreover, we prove a limit result for the sequence of
consecutive macroscopic resident traits that is similar to the so-called
trait-substitution-sequence.


  The ecological dynamics of interacting predator and prey populations can
display sustained oscillations, as for instance predicted by the
Rosenzweig-MacArthur predator-prey model. The presence of demographic
stochasticity, due to the finiteness of population sizes, alters the amplitude
and frequency of these oscillations. Here we present a method for
characterizing the effects of demographic stochasticity on the limit cycle
attractor of the Rosenzweig-MacArthur. We show that an angular Brownian motion
well describes the frequency oscillations. In the vicinity of the bifurcation
point, we obtain an analytical approximation for the angular diffusion
constant. This approximation accurately captures the effect of demographic
stochasticity across parameter values.


  Single-cell RNA sequencing (scRNAseq) has revolutionized our ability to
explore biological systems by enabling the study of gene expression at the
individual cell level. However, handling and analyzing this data often require
specialized expertise. In this contribution, we present scX, an R package built
on top of the Shiny framework, designed to simplify the analysis, exploration,
and visualization of single-cell experiments. scX offers straightforward access
to essential scRNAseq analyses, encompassing marker identification, gene
expression profiling, and differential gene expression analysis. Implemented as
a local web application with an intuitive graphical interface, scX allows users
to create customized, publication-ready plots. Additionally, it seamlessly
integrates with popular single-cell Seurat and SingleCellExperiment R objects,
facilitating the rapid processing and visualization of diverse datasets. In
summary, scX serves as a valuable tool for effortless exploration and sharing
of single-cell data, alleviating some of the complexities associated with
scRNAseq analysis.


  Protein engineering is an established method for tailoring enzymatic
reactivity. A commonly used method is directed evolution, where the mutagenesis
and natural selection process is mimicked and accelerated in the laboratory.
Here, we describe a reliable method for generating saturation mutagenesis
libraries by golden gate cloning in a broad host range plasmid containing the
pBBR1 replicon. The applicability is demonstrated by generating a mutant
library of the iron nitrogenase gene cluster (anfHDGK) of Rhodobacter
capsulatus, which is subsequently screened for the improved formation of
molecular hydrogen.


  In this paper, we consider a two-player two-strategy game with random payoffs
in a population subdivided into $d$ demes, each containing $N$ individuals at
the beginning of any given generation and experiencing local extinction and
recolonization with some fixed probability $m$ after reproduction and selection
among offspring. Within each deme, offspring engage in random pairwise
interactions, and the payoffs are assumed to have means and variances
proportional to the inverse of the population size. By verifying the conditions
given in Ethier and Nagylaki (1980) to approximate Markov chains with two time
scales, we establish that the discrete-time evolutionary dynamics with $Nd$
generations as unit of time converges to a continuous-time diffusion as
$d\rightarrow\infty$. The infinitesimal mean and variance of this diffusion are
expressed in terms of the population-scaled means and variances of the payoffs
besides identity-by-descent measures between offspring in the same deme in a
neutral population. We show that the probability for a strategy to fix in the
population starting from an initial frequency $(Nd)^{-1}$ generally increases
as the payoffs to that strategy exhibit less variability or the payoffs to the
other strategy more variability. As a result, differences in variability can
make this fixation probability for cooperation larger than the corresponding
one for defection. As the deme-scaled extinction rate $\nu=mN$ decreases for
$N$ large enough and $m$ small enough, creating a higher level of identity
among offspring within demes, the differences between the population-scaled
variances of the payoffs for interacting offspring of different types increases
this effect to a greater extent than the differences for interacting offspring
of the same type.


  In this paper, we study a finite population undergoing discrete,
nonoverlapping generations, that is structured into $D$ demes, each containing
$N$ individuals of two possible types, $A$ and $B$, whose viability
coefficients, $s_A$ and $s_B$, respectively, vary randomly from one generation
to the next. We assume that the means, variances and covariance of the
viability coefficients are inversely proportional to the number of demes $D$,
while higher-order moments are negligible in comparison to $1/D$. We use a
discrete-time Markov chain with two time scales to model the evolutionary
process, and we demonstrate that as the number of demes $D$ approaches
infinity, the accelerated Markov chain converges to a diffusion process for any
deme size $N\geq 2$. This diffusion process allows us to evaluate the fixation
probability of type $A$ following its introduction as a single mutant in a
population that was fixed for type $B$. We explore the impact of increasing the
variability in the viability coefficients on this fixation probability. At
least when $N$ is large enough, it is shown that increasing this variability
for type $B$ or decreasing it for type $A$ leads to an increase in the fixation
probability of a single $A$. The effect of the population-scaled variances,
$\sigma^2_A$ and $\sigma^2_B$, can even cancel the effects of the
population-scaled means, $\mu_A$ and $\mu_B$. We also show that the fixation
probability of a single $A$ increases as the deme-scaled migration rate
increases. Moreover, this probability is higher for type $A$ than for type $B$
if the population-scaled geometric mean is higher for type $A$ than for type
$B$, which means that $\mu_A-\sigma_A^2/2>\mu_B-\sigma_B^2/2$.


  MolecularWebXR is our new website for education, science communication and
scientific peer discussion in chemistry and biology built on WebXR. It
democratizes multi-user, inclusive virtual reality (VR) experiences that are
deeply immersive for users wearing high-end headsets, yet allow participation
by users with consumer devices such as smartphones, possibly inserted into
cardboard goggles for immersivity, or even computers or tablets. With no
installs as it is all web-served, MolecularWebXR enables multiple users to
simultaneously explore, communicate and discuss chemistry and biology concepts
in immersive 3D environments, manipulating objects with their bare hands,
either present in the same real space or scattered throughout the globe thanks
to built-in audio features. A series of preset rooms cover educational material
on chemistry and structural biology, and an empty room can be populated with
material prepared ad hoc using moleculARweb's VMD-based PDB2AR tool. We
verified ease of use and versatility by users aged 12-80 in entirely virtual
sessions or mixed real-virtual sessions at science outreach events, student
instruction, scientific collaborations, and conference lectures. MolecularWebXR
is available for free use without registration at https://molecularwebxr.org,
and a blog post version of this preprint with embedded videos is available at
https://go.epfl.ch/molecularwebxr-blog-post.


  The adaptive immune system engages in an arms race with evolving viruses,
trying to generate new responses to viral strains that continually move away
from the set of variants that have already elicited a functional immune
response. In previous work, it has been argued that this dynamical process can
lead to a propagating pulse of ever-changing viral population and concomitant
immune response. Here, we introduce a new stochastic model of viral-host
coevolution and a deterministic approximation thereof with which to study these
dynamics. We show that there is indeed a possible pulse solution, but for a
large host population size, the pulse becomes unstable to the generation of new
infections in its wake, leading to an extended endemic infection pattern. This
time-dependent endemic pattern eventually reaches a fluctuating steady-state
when the infected population size reaches a large enough fraction of the total
number of hosts.


  This paper presents a modeling framework for tracking the spawning migration
of the capelin, which is a fish species in the Barents Sea. The framework
combines an individual-based model (IBM) with artificial neural networks
(ANNs). The ANNs determine the direction of the fish's movement based on local
environmental information, while a genetic algorithm and fitness function
assess the suitability of the proposed directions. The framework's efficacy is
demonstrated by comparing the spatial distributions of modeled and empirical
potential spawners.
  The proposed model successfully replicates the southeastward movement of
capelin during their spawning migration, accurately capturing the distribution
of spawning fish over historical spawning sites along the eastern coast of
northern Norway.
  Furthermore, the paper compares three migration models: passive swimmers,
taxis movement based on temperature gradients, and restricted-area search,
along with our proposed approach. The results reveal that our approach
outperforms the other models in mimicking the migration pattern. Most spawning
stocks managed to reach the spawning sites, unlike the other models where water
currents played a significant role in pushing the fish away from the coast. The
temperature gradient detection model and restricted-area search model are found
to be inadequate for accurately simulating capelin spawning migration in the
Barents Sea due to complex oceanographic conditions.


  Recent progress has boosted the resolving power of optical microscopies to
spatial dimensions well below the diffraction limit. Yet, axial
super-resolution and axial single-molecule localisation typically require more
complicated implementations than their lateral counterparts. In the present
work, we propose a simple solution for axial metrology by providing a
multi-layered single-excitation, dual-emission test slide, in which axial
distance is colour-encoded. Our test slide combines on a standard microscope
coverslip substrate two flat, thin, uniform and brightly emitting fluorophore
layers, separated by a nanometric transparent spacer layer having a refractive
index close to a biological cell. The ensemble is sealed in an index-matched
protective polymer. As a proof-of-principle, we estimate the light confinement
resulting from evanescent-wave excitation in total internal reflection
fluorescence (TIRF) microscopy. Our test sample permits, even for the
non-expert user, a facile axial metrology at the sub-100-nm scale, a critical
requirement for axial super-resolution, as well as near-surface imaging,
spectroscopy and sensing.


  We analyse an extended SIRS epidemic model in which immunity at the
individual level wanes gradually at exponential rate, but where the waning rate
may differ between individuals, for instance as an effect of differences in
immune systems. The model also includes vaccination schemes aimed to reach and
maintain herd immunity. We consider both the informed situation where the
individual waning parameters are known, thus allowing selection of vaccinees
being based on both time since last vaccination as well as on the individual
waning rate, and the more likely uninformed situation where individual waning
parameters are unobserved, thus only allowing vaccination schemes to depend on
time since last vaccination. The optimal vaccination policies for both the
informed and uniformed heterogeneous situation are derived and compared with
the homogeneous waning model (meaning all individuals have the same immunity
waning rate), as well as to the classic SIRS model where immunity at the
individual level drops from complete immunity to complete susceptibility in one
leap. It is shown that the classic SIRS model requires least vaccines, followed
by the SIRS with homogeneous gradual waning, followed by the informed situation
for the model with heterogeneous gradual waning. The situation requiring most
vaccines for herd immunity is the most likely scenario, that immunity wanes
gradually with unobserved individual heterogeneity. For parameter values chosen
to mimic COVID-19 and assuming perfect initial immunity and cumulative immunity
of 12 months, the classic homogeneous SIRS epidemic suggests that vaccinating
individuals every 15 months is sufficient to reach and maintain herd immunity,
whereas the uninformed case for exponential waning with rate heterogeneity
corresponding to a coefficient of variation being 0.5, requires that
individuals instead need to be vaccinated every 4.4 months.


  We study the tipping point collective dynamics of an adaptive
susceptible-infected-susceptible (SIS) epidemiological network in a
data-driven, machine learning-assisted manner. We identify a
parameter-dependent effective stochastic differential equation (eSDE) in terms
of physically meaningful coarse mean-field variables through a deep-learning
ResNet architecture inspired by numerical stochastic integrators. We construct
an approximate effective bifurcation diagram based on the identified drift term
of the eSDE and contrast it with the mean-field SIS model bifurcation diagram.
We observe a subcritical Hopf bifurcation in the evolving network's effective
SIS dynamics, that causes the tipping point behavior; this takes the form of
large amplitude collective oscillations that spontaneously -- yet rarely --
arise from the neighborhood of a (noisy) stationary state. We study the
statistics of these rare events both through repeated brute force simulations
and by using established mathematical/computational tools exploiting the
right-hand-side of the identified SDE. We demonstrate that such a collective
SDE can also be identified (and the rare events computations also performed) in
terms of data-driven coarse observables, obtained here via manifold learning
techniques, in particular Diffusion Maps. The workflow of our study is
straightforwardly applicable to other complex dynamics problems exhibiting
tipping point dynamics.


  We analyse a non-local parabolic integro-differential equation modelling the
evolutionary dynamics of a phenotypically-structured population in a changing
environment. Such models arise in the study of species adapting to climate
change, or cancer adapting to therapy. Our results concern the long-time
behaviour, in the small mutation limit, of the model. The main novelty of our
work is that the time- and trait-dependent per capita growth rate is
characterised by having multiple (locally) optimal traits which shift at
possibly different velocities. When the velocities are the same we find that
the solution concentrates on a point-set which depends on the shifting speed.
These points can be thought of as "lagged optima" and the fitness value at each
lagged optima is the "lagged fitness". When the velocities of the optimal
traits are different and each has a potentially different optimal fitness, we
find that the solution in fact concentrates as a Dirac delta function on the
positive lagged optimum with maximum lagged fitness. Our results imply that in
populations undergoing competition in temporally changing environments, both
the true optimal fitness and the required rate of adaptation for each of the
diverging optimal traits contribute to the eventual dominance of one trait.


  Aligning EM density maps and fitting atomic models are essential steps in
single particle cryogenic electron microscopy (cryo-EM), with recent methods
leveraging various algorithms and machine learning tools. As aligning maps
remains challenging in the presence of a map that only partially fits the other
(e.g. one subunit), we here propose a new procedure, EMPOT (EM Partial
alignment with Optimal Transport), for partial alignment of 3D maps. EMPOT
first finds a coupling between 3D point-cloud representations, which is
associated with their so-called unbalanced Gromov Wasserstein divergence, and
second, uses this coupling to find an optimal rigid body transformation. Upon
running and benchmarking our method with experimental maps and structures, we
show that EMPOT outperforms standard methods for aligning subunits of a protein
complex and fitting atomic models to a density map, suggesting potential
applications of Partial Optimal Transport for improving Cryo-EM pipelines.


  Biodiversity is declining at an unprecedented rate, impacting ecosystem
services necessary to ensure food, water, and human health and well-being.
Understanding the distribution of species and their habitats is crucial for
conservation policy planning. However, traditional methods in ecology for
species distribution models (SDMs) generally focus either on narrow sets of
species or narrow geographical areas and there remain significant knowledge
gaps about the distribution of species. A major reason for this is the limited
availability of data traditionally used, due to the prohibitive amount of
effort and expertise required for traditional field monitoring. The wide
availability of remote sensing data and the growing adoption of citizen science
tools to collect species observations data at low cost offer an opportunity for
improving biodiversity monitoring and enabling the modelling of complex
ecosystems. We introduce a novel task for mapping bird species to their
habitats by predicting species encounter rates from satellite images, and
present SatBird, a satellite dataset of locations in the USA with labels
derived from presence-absence observation data from the citizen science
database eBird, considering summer (breeding) and winter seasons. We also
provide a dataset in Kenya representing low-data regimes. We additionally
provide environmental data and species range maps for each location. We
benchmark a set of baselines on our dataset, including SOTA models for remote
sensing tasks. SatBird opens up possibilities for scalably modelling properties
of ecosystems worldwide.


  Can a micron sized sack of interacting molecules autonomously learn an
internal model of a complex and fluctuating environment? We draw insights from
control theory, machine learning theory, chemical reaction network theory, and
statistical physics to develop a general architecture whereby a broad class of
chemical systems can autonomously learn complex distributions. Our construction
takes the form of a chemical implementation of machine learning's optimization
workhorse: gradient descent on the relative entropy cost function. We show how
this method can be applied to optimize any detailed balanced chemical reaction
network and that the construction is capable of using hidden units to learn
complex distributions. This result is then recast as a form of integral
feedback control. Finally, due to our use of an explicit physical model of
learning, we are able to derive thermodynamic costs and trade-offs associated
to this process.


  Obstructive sleep apnea (OSA) is a sleep respiratory disease characterized by
sleep snoring accompanied by apnea and daytime sleeplessness. It is a complex
disease, with the multifactorial etiology, and the pathology is incompletely
understood. Genioglossus (GG), the largest dilator of upper airway, whose
fatigue is strongly correlated to onset of OSA. This brief review was to
investigate the pathogenesis of OSA targeting on GG from different risk factors
as gender, obesity, and aging, and the molecular mechanism of GG injury in OSA
pathogenesis. We hope to find the targeted molecular mechanism on GG in OSA
treatment.


  Patterns in a semiarid ecosystem are important because they directly and
indirectly affect ecological processes, biodiversity, and ecosystem resilience.
Understanding the causes and effects of these patterns is critical for
long-term land surface management and conservation efforts in semiarid regions,
which are especially sensitive to climate change and human-caused disturbances.
It is known that there is a regular connection between the vegetation and the
living species in a habitat since some animals evolved to live in a semiarid
ecosystem and rely on plants for food. In this work, we have constructed a
coupled mathematical model to connect the water resource, vegetation and living
organisms and have investigated how the soil surface affects the resulting
patterns for the long term. This study contributes to a better understanding of
ecological patterns and processes in semiarid environments by shedding light on
the complex interaction mechanisms that depend on the structure of semiarid
ecosystems. The findings provide further critical insight into the influence of
efforts for improving ecosystem resilience and adjusting to the challenges
posed by climate change and human activities.


  Switching interacting particle systems studied in probability theory are the
stochastic processes of hopping particles on a lattice consisting of slow
particles and fast particles, where the switching between these two types of
particles happens randomly with a given transition rate. In the present paper,
we show that such stochastic processes of many particles are useful in modeling
group behaviors of ants. Recently the situation-depending switching between two
distinct types of primarily relied cues for ants in selecting foraging paths
has been experimentally studied by the research group of the last author of
this paper. We propose a discrete-time interacting random-walk model on a
square lattice, in which two kinds of hopping rules are included. We report the
numerical simulation results which exhibit the global changes in selected
homing paths from trailing paths of the `pheromone road' to
nearly-optimal-paths depending on the switching parameters. By introducing two
kinds of order parameters characterizing the switching-parameter dependence of
homing-duration distributions, we discuss such global changes as phase
transitions realized in path selections of ants. Critical phenomena associated
with the continuous phase transitions are also studied.


  Random matrix theory allows for the deduction of stability criteria for
complex systems using only a summary knowledge of the statistics of the
interactions between components. As such, results like the well-known
elliptical law are applicable in a myriad of different contexts. However, it is
often assumed that all components of the complex system in question are
statistically equivalent, which is unrealistic in many applications. Here, we
introduce the concept of a finely-structured random matrix. These are random
matrices with element-specific statistics, which can be used to model systems
in which the individual components are statistically distinct. By supposing
that the degree of `fine structure' in the matrix is small, we arrive at a
succinct `modified' elliptical law. We demonstrate the direct applicability of
our results to the niche and cascade models in theoretical ecology, as well as
a model of a neural network, and a directed network with arbitrary degree
distribution. The simple closed form of our central results allow us to draw
broad qualitative conclusions about the effect of fine structure on stability.


  Metagenomics, the study of genome sequences of diverse organisms cohabiting
in a shared environment, has experienced significant advancements across
various medical and biological fields. Metagenomic analysis is crucial, for
instance, in clinical applications such as infectious disease screening and the
diagnosis and early detection of diseases such as cancer. A key task in
metagenomics is to determine the species present in a sample and their relative
abundances. Currently, the field is dominated by either alignment-based tools,
which offer high accuracy but are computationally expensive, or alignment-free
tools, which are fast but lack the needed accuracy for many applications. In
response to this dichotomy, we introduce MetaFast, a tool based on heuristics,
to achieve a fundamental improvement in accuracy-runtime tradeoff over existing
methods. MetaFast delivers accuracy comparable to the alignment-based and
highly accurate tool Metalign but with significantly enhanced efficiency. In
MetaFast, we accelerate memory-frugal reference database indexing and
filtering. We further employ heuristics to accelerate read mapping. Our
evaluation demonstrates that MetaFast achieves a 4x speedup over Metalign
without compromising accuracy. MetaFast is publicly available at
https://github.com/CMU-SAFARI/MetaFast.


  DNA is not only the genetic material of life, but also a favorable material
for a new computing model. Various research works based on DNA computing have
been carried out in recent years. DNA sequence design is the foundation of such
research. The sequence quality directly affects the universality, robustness,
and stability of DNA computing. How to design DNA sequences depends on the
biological properties and target requirements, which is a typical combinatorial
optimization problem. In this paper, in order to design DNA sequences with
high-quality, we propose a novel meta-heuristic evolutionary algorithm, termed
the static virus spread algorithm (SVS). Through this algorithm, we focus on
the constraints of universal DNA sequence design and produce a large number of
DNA sequences with non-complementarity and small difference in melting
temperature as the objectives, and fully considering the balanced proportion of
the four bases. The computer simulation and polyacrylamide gel electrophoresis
experiments show that the high-quality DNA sequences designed by this algorithm
are effective, which is expected to provide a convenient tool for sequence
preparation before DNA biochemical operations.


  While molecular pre-training has shown great potential in enhancing drug
discovery, the lack of a solid physical interpretation in current methods
raises concerns about whether the learned representation truly captures the
underlying explanatory factors in observed data, ultimately resulting in
limited generalization and robustness. Although denoising methods offer a
physical interpretation, their accuracy is often compromised by ad-hoc noise
design, leading to inaccurate learned force fields. To address this limitation,
this paper proposes a new method for molecular pre-training, called sliced
denoising (SliDe), which is based on the classical mechanical intramolecular
potential theory. SliDe utilizes a novel noise strategy that perturbs bond
lengths, angles, and torsion angles to achieve better sampling over
conformations. Additionally, it introduces a random slicing approach that
circumvents the computationally expensive calculation of the Jacobian matrix,
which is otherwise essential for estimating the force field. By aligning with
physical principles, SliDe shows a 42\% improvement in the accuracy of
estimated force fields compared to current state-of-the-art denoising methods,
and thus outperforms traditional baselines on various molecular property
prediction tasks.


  A mathematical development of the Walsh transform, Walsh basis, and Walsh
coefficients is given. The author was prompted to write this by a wish to give
a unified treatment of epistatic coordinates as they are used in evolutionary
biology. At the end of the article, opinions are expressed regarding the
usefulness of these concepts for the practical researcher.


  We present and analyze an actively controlled
Susceptible-Infected-Susceptible (actSIS) model of interconnected populations
to study how risk aversion strategies, such as social distancing, affect
network epidemics. A population using a risk aversion strategy reduces its
contact rate with other populations when it perceives an increase in infection
risk. The network actSIS model relies on two distinct networks. One is a
physical contact network that defines which populations come into contact with
which other populations and thus how infection spreads. The other is a
communication network, such as an online social network, that defines which
populations observe the infection level of which other populations and thus how
information spreads. We prove that the model, with these two networks and
populations using risk aversion strategies, exhibits a transcritical
bifurcation in which an endemic equilibrium emerges. For regular graphs, we
prove that the endemic infection level is uniform across populations and
reduced by the risk aversion strategy, relative to the network SIS endemic
level. We show that when communication is sufficiently sparse, this initially
stable equilibrium loses stability in a secondary bifurcation. Simulations show
that a new stable solution emerges with nonuniform infection levels.


  Investigations of infectious disease outbreaks often focus on identifying
place- and context-dependent factors responsible for emergence and spread,
resulting in phenomenological narratives ill-suited to developing generalizable
predictive and preventive measures. We contend that case-control hypothesis
testing is a more powerful framework for epidemiological investigation. The
approach, widely used in medical research, involves identifying
counterfactuals, with case-control comparisons drawn to test hypotheses about
the conditions that manifest outbreaks. Here we outline the merits of applying
a case-control framework as epidemiological study design. We first describe a
framework for iterative multidisciplinary interrogation to discover minimally
sufficient sets of factors that can lead to disease outbreaks. We then lay out
how case-control comparisons can respectively center on pathogen(s), factor(s),
or landscape(s) with vignettes focusing on pathogen transmission. Finally, we
consider how adopting case-control approaches can promote evidence-based
decision making for responding to and preventing outbreaks.


  We study an individual-based stochastic epidemic model in which infected
individuals become susceptible again following each infection (generalized SIS
model). Specifically, after each infection, the infectivity is a random
function of the time elapsed since the infection, and each recovered individual
loses immunity gradually (equivalently, becomes gradually susceptible) after
some time according to a random susceptibility function. The epidemic dynamics
is described by the average infectivity and susceptibility processes in the
population together with the numbers of infected and susceptible/uninfected
individuals. In \cite{forien-Zotsa2022stochastic}, a functional law of large
numbers (FLLN) is proved as the population size goes to infinity, and
asymptotic endemic behaviors are also studied. In this paper, we prove a
functional central limit theorem (FCLT) for the stochastic fluctuations of the
epidemic dynamics around the FLLN limit. The FCLT limit for the aggregate
infectivity and susceptibility processes is given by a system of stochastic
non-linear integral equation driven by a two-dimensional Gaussian process.


  We demonstrate the fabrication of sharp nanopillars of high aspect ratio onto
specialized atomic force microscopy (AFM) microcantilevers and their use for
high-speed AFM of DNA and nucleoproteins in liquid. The fabrication technique
uses localized charged-particle-induced deposition with either a focused beam
of helium ions or electrons in a helium ion microscope (HIM) or scanning
electron microscope (SEM). This approach enables customized growth onto
delicate substrates with nanometer-scale placement precision and in-situ
imaging of the final tip structures using the HIM or SEM. Tip radii of <10 nm
are obtained and the underlying microcantilever remains intact. Instead of the
more commonly used organic precursors employed for bio-AFM applications, we use
an organometallic precursor (tungsten hexacarbonyl) resulting in
tungsten-containing tips. Transmission electron microscopy reveals a thin layer
of carbon on the tips. Consequently, the interaction of the new tips with
biological specimens is likely very similar to that of standard carbonaceous
tips, with the added benefit of robustness. A further advantage of the
organometallic tips is that compared to carbonaceous tips they better withstand
UV-ozone cleaning treatments to remove residual organic contaminants between
experiments. Such contaminants are unavoidable during the scanning of soft
biomolecules in liquid. Our tips can also be grown onto the blunted tips of
previously-used cantilevers, thus providing a means to recycle specialized
cantilevers and restore their performance to the original manufacturer
specifications. Finally, a focused helium ion beam milling technique to reduce
the tip radii and thus further improve lateral spatial resolution in the AFM
scans is demonstrated.


  This paper presents the Ensemble Nucleotide Byte-level Encoder-Decoder
(ENBED) foundation model, analyzing DNA sequences at byte-level precision with
an encoder-decoder Transformer architecture. ENBED uses a sub-quadratic
implementation of attention to develop an efficient model capable of
sequence-to-sequence transformations, generalizing previous genomic models with
encoder-only or decoder-only architectures. We use Masked Language Modeling to
pre-train the foundation model using reference genome sequences and apply it in
the following downstream tasks: (1) identification of enhancers, promotors and
splice sites, (2) identification of biological function annotations of genomic
sequences, (3) recognition of sequences containing base call mismatches and
insertion/deletion errors, an advantage over tokenization schemes involving
multiple base pairs, which lose the ability to analyze with byte-level
precision, and (4) generating mutations of the Influenza virus using the
encoder-decoder architecture and validating them against real-world
observations. In each of these tasks, we demonstrate significant improvement as
compared to the existing state-of-the-art results.


  Predicting patient features from single-cell data can unveil cellular states
implicated in health and disease. Linear models and average cell type
expressions are typically favored for this task for their efficiency and
robustness, but they overlook the rich cell heterogeneity inherent in
single-cell data. To address this gap, we introduce GMIL, a framework
integrating Generalized Linear Mixed Models (GLMM) and Multiple Instance
Learning (MIL), upholding the advantages of linear models while modeling
cell-state heterogeneity. By leveraging predefined cell embeddings, GMIL
enhances computational efficiency and aligns with recent advancements in
single-cell representation learning. Our empirical results reveal that GMIL
outperforms existing MIL models in single-cell datasets, uncovering new
associations and elucidating biological mechanisms across different domains.


  This paper introduces an efficient approach to reduce the computational cost
of simulating collective behaviors, such as fish schooling, using
Individual-Based Models (IBMs). The proposed technique employs adaptive and
dynamic load-balancing domain partitioning, which utilizes unsupervised
machine-learning models to cluster a large number of simulated individuals into
sub-schools based on their spatial-temporal locations. It also utilizes Voronoi
tessellations to construct non-overlapping simulation subdomains. This approach
minimizes agent-to-agent communication and balances the load both spatially and
temporally, ultimately resulting in reduced computational complexity.
  Experimental simulations demonstrate that this partitioning approach
outperforms the standard regular grid-based domain decomposition, achieving a
reduction in computational cost while maintaining spatial and temporal load
balance. The approach presented in this paper has the potential to be applied
to other collective behavior simulations requiring large-scale simulations with
a substantial number of individuals.


  Despite the breakthroughs in biomarker discovery facilitated by differential
gene analysis, challenges remain, particularly at the single-cell level.
Traditional methodologies heavily rely on user-supplied cell annotations,
focusing on individually expressed data, often neglecting the critical
interactions between biological conditions, such as healthy versus diseased
states. In response, here we introduce scBeacon, an innovative framework built
upon a deep contrastive siamese network. scBeacon pioneers an unsupervised
approach, adeptly identifying matched cell populations across varied
conditions, enabling a refined differential gene analysis. By utilizing a
VQ-VAE framework, a contrastive siamese network, and a greedy iterative
strategy, scBeacon effectively pinpoints differential genes that hold potential
as key biomarkers. Comprehensive evaluations on a diverse array of datasets
validate scBeacon's superiority over existing single-cell differential gene
analysis tools. Its precision and adaptability underscore its significant role
in enhancing diagnostic accuracy in biomarker discovery. With the emphasis on
the importance of biomarkers in diagnosis, scBeacon is positioned to be a
pivotal asset in the evolution of personalized medicine and targeted
treatments.


  The embedding problem of Markov matrices in Markov semigroups is a classic
problem that regained a lot of impetus and activities through recent needs in
phylogeny and population genetics. Here, we give an account for dimensions
$d\leqslant 4$, including a complete and simplified treatment of the case
$d=3$, and derive the results in a systematic fashion, with an eye on the
potential applications.
  Further, we reconsider the setup of the corresponding problem for
time-inhomogeneous Markov chains, which is needed for real-world applications
because transition rates need not be constant over time. Additional cases of
this more general embedding occur for any $d\geqslant 3$. We review the known
case of $d=3$ and describe the setting for future work on $d=4$.


  The review is devoted to analysis of mathematical models used for describing
epidemic processes. A main focus is done on the models that are based on
partial differential equations (PDEs), especially those that were developed and
used for the COVID-19 pandemic modelling. Our attention is paid preferable to
the studies in which not only results of numerical simulations are presented
but analytical results as well. In particular, travelling fronts (waves), exact
solutions, estimation of key epidemic parameters of the epidemic models with
governing PDEs (typically reaction-diffusion equations) are discussed. The
review may serve as a valuable source for researchers and practitioners in the
field of mathematical modelling in epidemiology.


  The replicator equation in evolutionary game theory describes the change in a
population's behaviors over time given suitable incentives. It arises when
individuals make decisions using a simple learning process - imitation. A
recent emerging framework builds upon this standard model by incorporating
game-environment feedback, in which the population's actions affect a shared
environment, and in turn, the changing environment shapes incentives for future
behaviors. In this paper, we investigate game-environment feedback when
individuals instead use a boundedly rational learning rule known as logit
learning. We characterize the resulting system's complete set of fixed points
and their local stability properties, and how the level of rationality
determines overall environmental outcomes in comparison to imitative learning
rules. We identify a large parameter space for which logit learning exhibits a
wide range of dynamics as the rationality parameter is increased from low to
high. Notably, we identify a bifurcation point at which the system exhibits
stable limit cycles. When the population is highly rational, the limit cycle
collapses and a tragedy of the commons becomes stable.


  Motivation: Curation of literature in life sciences is a growing challenge.
The continued increase in the rate of publication, coupled with the relatively
fixed number of curators worldwide presents a major challenge to developers of
biomedical knowledgebases. Very few knowledgebases have resources to scale to
the whole relevant literature and all have to prioritise their efforts.
  Results: In this work, we take a first step to alleviating the lack of
curator time in RNA science by generating summaries of literature for
non-coding RNAs using large language models (LLMs). We demonstrate that
high-quality, factually accurate summaries with accurate references can be
automatically generated from the literature using a commercial LLM and a chain
of prompts and checks. Manual assessment was carried out for a subset of
summaries, with the majority being rated extremely high quality. We also
applied the most commonly used automated evaluation approaches, finding that
they do not correlate with human assessment. Finally, we apply our tool to a
selection of over 4,600 ncRNAs and make the generated summaries available via
the RNAcentral resource. We conclude that automated literature summarization is
feasible with the current generation of LLMs, provided careful prompting and
automated checking are applied.
  Availability: Code used to produce these summaries can be found here:
https://github.com/RNAcentral/litscan-summarization and the dataset of contexts
and summaries can be found here:
https://huggingface.co/datasets/RNAcentral/litsumm-v1. Summaries are also
displayed on the RNA report pages in RNAcentral (https://rnacentral.org/)


  Synthetic biologists and molecular programmers design novel nucleic acid
reactions, with many potential applications. Good visualization tools are
needed to help domain experts make sense of the complex outputs of folding
pathway simulations of such reactions. Here we present ViDa, a new approach for
visualizing DNA reaction folding trajectories over the energy landscape of
secondary structures. We integrate a deep graph embedding model with common
dimensionality reduction approaches, to map high-dimensional data onto 2D
Euclidean space. We assess ViDa on two well-studied and contrasting DNA
hybridization reactions. Our preliminary results suggest that ViDa's
visualization successfully separates trajectories with different folding
mechanisms, thereby providing useful insight to users, and is a big improvement
over the current state-of-the-art in DNA kinetics visualization.


  Single-cell RNA sequencing (scRNA-seq) is important to transcriptomic
analysis of gene expression. Recently, deep learning has facilitated the
analysis of high-dimensional single-cell data. Unfortunately, deep learning
models may leak sensitive information about users. As a result, Differential
Privacy (DP) is increasingly used to protect privacy. However, existing DP
methods usually perturb whole neural networks to achieve differential privacy,
and hence result in great performance overheads. To address this challenge, in
this paper, we take advantage of the uniqueness of the autoencoder that it
outputs only the dimension-reduced vector in the middle of the network, and
design a Differentially Private Deep Contrastive Autoencoder Network (DP-DCAN)
by partial network perturbation for single-cell clustering. Since only partial
network is added with noise, the performance improvement is obvious and
twofold: one part of network is trained with less noise due to a bigger privacy
budget, and the other part is trained without any noise. Experimental results
of six datasets have verified that DP-DCAN is superior to the traditional DP
scheme with whole network perturbation. Moreover, DP-DCAN demonstrates strong
robustness to adversarial attacks. The code is available at
https://github.com/LFD-byte/DP-DCAN.


  Visualization tools can help synthetic biologists and molecular programmers
understand the complex reactive pathways of nucleic acid reactions, which can
be designed for many potential applications and can be modelled using a
continuous-time Markov chain (CTMC). Here we present ViDa, a new visualization
approach for DNA reaction trajectories that uses a 2D embedding of the
secondary structure state space underlying the CTMC model. To this end, we
integrate a scattering transform of the secondary structure adjacency, a
variational autoencoder, and a nonlinear dimensionality reduction method. We
augment the training loss with domain-specific supervised terms that capture
both thermodynamic and kinetic features. We assess ViDa on two well-studied DNA
hybridization reactions. Our results demonstrate that the domain-specific
features lead to significant quality improvements over the state-of-the-art in
DNA state space visualization, successfully separating different folding
pathways and thus providing useful insights into dominant reaction mechanisms.


  Clinical variant classification of pathogenic versus benign genetic variants
remains a pivotal challenge in clinical genetics. Recently, the proposition of
protein language models has improved the generic variant effect prediction
(VEP) accuracy via weakly-supervised or unsupervised training. However, these
VEPs are not disease-specific, limiting their adaptation at point-of-care. To
address this problem, we propose a disease-specific \textsc{pro}tein language
model for variant \textsc{path}ogenicity, termed ProPath, to capture the
pseudo-log-likelihood ratio in rare missense variants through a siamese
network. We evaluate the performance of ProPath against pre-trained language
models, using clinical variant sets in inherited cardiomyopathies and
arrhythmias that were not seen during training. Our results demonstrate that
ProPath surpasses the pre-trained ESM1b with an over $5\%$ improvement in AUC
across both datasets. Furthermore, our model achieved the highest performances
across all baselines for both datasets. Thus, our ProPath offers a potent
disease-specific variant effect prediction, particularly valuable for disease
associations and clinical applicability.


  Spatial agent-based models are increasingly used to investigate the evolution
of solid tumours subject to localised cell-cell interactions and
microenvironmental heterogeneity. Here we present a non-technical step by step
guide to developing such a model from first principles, aimed at both aspiring
modellers and other biologists and oncologists who wish to understand the
assumptions and limitations of this approach. Stressing the importance of
tailoring the model structure to that of the biological system, we describe
methods of increasing complexity, from the basic Eden growth model up to
off-lattice simulations with diffusible factors. We examine choices that
unavoidably arise in model design, such as implementation, parameterisation,
visualisation, and reproducibility. Each topic is illustrated with examples
drawn from recent research studies and state of the art modelling platforms. We
emphasise the benefits of simpler models that aim to match the complexity of
the phenomena of interest, rather than that of the entire biological system.


  Most sequence sketching methods work by selecting specific $k$-mers from
sequences so that the similarity between two sequences can be estimated using
only the sketches. Estimating sequence similarity is much faster using sketches
than using sequence alignment, hence sketching methods are used to reduce the
computational requirements of computational biology software packages.
Applications using sketches often rely on properties of the $k$-mer selection
procedure to ensure that using a sketch does not degrade the quality of the
results compared with using sequence alignment. In particular the window
guarantee ensures that no long region of the sequence goes unrepresented in the
sketch.
  A sketching method with a window guarantee corresponds to a Decycling Set,
aka an unavoidable sets of $k$-mers. Any long enough sequence must contain a
$k$-mer from any decycling set (hence, it is unavoidable). Conversely, a
decycling set defines a sketching method by selecting the $k$-mers from the
set. Although current methods use one of a small number of sketching method
families, the space of decycling sets is much larger, and largely unexplored.
Finding decycling sets with desirable characteristics is a promising approach
to discovering new sketching methods with improved performance (e.g., with
small window guarantee).
  The Minimum Decycling Sets (MDSs) are of particular interest because of their
small size. Only two algorithms, by Mykkeltveit and Champarnaud, are known to
generate two particular MDSs, although there is a vast number of alternative
MDSs. We provide a simple method that allows one to explore the space of MDSs
and to find sets optimized for desirable properties. We give evidence that the
Mykkeltveit sets are close to optimal regarding one particular property, the
remaining path length.


  Infectious zoonotic disease emergence, through spillover events, is of global
concern and has the potential to cause significant harm to society, as recently
demonstrated by COVID-19. More than 70% of the 400 infectious diseases that
emerged in the past five decades have a zoonotic origin, including all recent
pandemics. There have been several approaches used to predict the risk of
spillover through some of the known or suspected infectious disease emergence
drivers, largely using correlative approaches. Here, we predict the spatial
distribution of spillover risk by approximating general transmission through
animal and human interactions. These mass action interactions are approximated
through the multiplication of the spatial distribution of zoonotic viral
diversity and human population density. Although our results indicate higher
risk in regions along the equator and in Southeast Asia where both viral
diversity and human population density are high, it should be noted that this
is primarily a conceptual exercise. We compared our spillover risk map to key
factors, including the model inputs of zoonotic viral diversity estimate map,
human population density map, and the spatial distribution of species richness.
Despite the limitations of this approach, this viral spillover map is a step
towards developing a more comprehensive spillover risk prediction system to
inform global monitoring.


  In this study, we employed a Raman spectroscopic analysis of the amide I band
of polymerized samples of tubulin exposed to pulsed low-intensity NIR radiation
(810 nm, 10 Hz, 22.5 J/cm$^{2}$ dose). Peaks in the Raman fingerprint region
(300$\unicode{x2013}$1900 cm$^{-1}$), in particular, in the amide I band
(1600$\unicode{x2013}$1700 cm$^{-1}$), can be used to quantify the percentage
of protein secondary structures. Under this band, hidden signals of
$\mathrm{C}$=$\mathrm{O}$ stretching, belonging to different structures, are
superimposed$\unicode{x2014}$producing a complex signal as a result. An
accurate decomposition of the amide I band is therefore required for the
reliable analysis of the conformation of proteins, which we achieved through a
straightforward method employing a Voigt profile. This approach was validated
through secondary structure analyses of unexposed control samples, for which
comparisons with other values available in the literature could be conducted.
Subsequently, using this validated method, we present novel findings of
statistically significant alterations in the secondary structures of
NIR-exposed tubulin, characterized by a notable decrease in alpha-helix content
and a concurrent increase in beta-sheets compared to the control samples. The
alpha-helix to beta-sheet transition suggests that PBM reduces microtubule
stability and introduces dynamism to allow for the remodeling and,
consequently, refreshing of microtubule structures. This newly discovered
mechanism could have implications for reducing the risks associated with brain
aging, including neurodegenerative diseases like Alzheimer's disease.


  Individual-based models of contagious processes are useful for predicting
epidemic trajectories and informing intervention strategies. In such models,
the incorporation of contact network information can capture the non-randomness
and heterogeneity of realistic contact dynamics. In this paper, we consider
Bayesian inference on the spreading parameters of an SIR contagion on a known,
static network, where information regarding individual disease status is known
only from a series of tests (positive or negative disease status). When the
contagion model is complex or information such as infection and removal times
is missing, the posterior distribution can be difficult to sample from.
Previous work has considered the use of Approximate Bayesian Computation (ABC),
which allows for simulation-based Bayesian inference on complex models.
However, ABC methods usually require the user to select reasonable summary
statistics. Here, we consider an inference scheme based on the Mixture Density
Network compressed ABC (MDN-ABC), which minimizes the expected posterior
entropy in order to learn informative summary statistics. This allows us to
conduct Bayesian inference on the parameters of a partially observed contagious
process while also circumventing the need for manual summary statistic
selection. This methodology can be extended to incorporate additional
simulation complexities, including behavioral change after positive tests or
false test results.


  Microbial communities are diverse biological systems that include taxa from
across multiple kingdoms of life. Notably, interactions between bacteria and
fungi play a significant role in determining community structure. However,
these statistical associations across kingdoms are more difficult to infer than
intra-kingdom associations due to the nature of the data involved using
standard network inference techniques. We quantify the challenges of
cross-kingdom network inference from both a theoretical and practical viewpoint
using synthetic and real-world microbiome data. We detail the theoretical issue
presented by combining compositional data-sets drawn from the same environment,
e.g. 16S and ITS sequencing of a single set of samples, and survey common
network inference techniques for their ability to handle this error. We then
test these techniques for the accuracy and usefulness of their intra- and
inter-kingdom associations by inferring networks from a set of simulated
samples for which a ground-truth set of associations is known. We show that
while two methods mitigate the error of cross-kingdom inference, there is
little difference between techniques for key practical applications including
identification of strong correlations and identification of possible keystone
taxa (i.e. hub nodes in the network). Furthermore, we identify a signature of
the error caused transkingdom network inference and demonstrate that it appears
in networks constructed using real-world environmental microbiome data.


  Models are often employed to integrate knowledge about epidemics across
scales and simulate disease dynamics. While these approaches have played a
central role in studying the mechanics underlying epidemics, we lack ways to
reliably predict how the relationship between virulence (the harm to hosts
caused by an infection) and transmission will evolve in certain virus-host
contexts. In this study, we invoke evolutionary invasion analysis -- a method
used to identify the evolution of uninvadable strategies in dynamical systems
-- to examine how the virulence-transmission dichotomy can evolve in models of
virus infections defined by different natural histories. We reveal that
peculiar ecologies drive different evolved relationships between virulence and
transmission. Specifically, we discover patterns of virulence evolution between
epidemics of various kinds (SARS-CoV-2 and hepatitis C virus) and that varying
definitions of virulence alter our predictions for how viruses will evolve. We
discuss the findings in light of contemporary conversations in the public
health sector around the possibility of predicting virus evolution and in more
extensive theoretical discussions involving virulence evolution in emerging
infectious diseases.


  In recent years, the scientific community has become increasingly interested
on peptides with non-canonical amino acids due to their superior stability and
resistance to proteolytic degradation. These peptides present promising
modifications to biological, pharmacological, and physiochemical attributes in
both endogenous and engineered peptides. Notwithstanding their considerable
advantages, the scientific community exhibits a conspicuous absence of an
effective pre-trained model adept at distilling feature representations from
such complex peptide sequences. We herein propose PepLand, a novel pre-training
architecture for representation and property analysis of peptides spanning both
canonical and non-canonical amino acids. In essence, PepLand leverages a
comprehensive multi-view heterogeneous graph neural network tailored to unveil
the subtle structural representations of peptides. Empirical validations
underscore PepLand's effectiveness across an array of peptide property
predictions, encompassing protein-protein interactions, permeability,
solubility, and synthesizability. The rigorous evaluation confirms PepLand's
unparalleled capability in capturing salient synthetic peptide features,
thereby laying a robust foundation for transformative advances in
peptide-centric research domains. We have made all the source code utilized in
this study publicly accessible via GitHub at
https://github.com/zhangruochi/pepland


  Aquaporins (AQPs) are recognized as transmembrane water channels that
facilitate selective water permeation through their monomeric pores. Among the
AQP family, AQP6 has a unique characteristic as an anion channel, which is
allosterically controlled by pH conditions and is eliminated by a single amino
acid mutation. However, the molecular mechanism of anion permeation through
AQP6 remains unclear. Using molecular dynamics simulations in the presence of a
transmembrane voltage utilizing an ion concentration gradient, we show that
chloride ions permeate through the pore corresponding to the central axis of
the AQP6 homotetramer. Under low pH conditions, a subtle opening of the
hydrophobic selective filter (SF), located near the extracellular part of the
central pore, becomes wetted and enables anion permeation. Our simulations also
indicate that a single mutation (N63G) in human AQP6, located at the central
pore, significantly reduces anion conduction, consistent with experimental
data. Moreover, we demonstrate the pH-sensing mechanism in which the
protonation of H184 and H189 under low pH conditions allosterically triggers
the gating of the SF region. These results suggest a unique pH-dependent
allosteric anion permeation mechanism in AQP6 and could clarify the role of the
central pore in some of the AQP tetramers.


  A novel mathematical model for the hypothalamic pituitary adrenal (HPA) axis
is proposed to comprehend the oscillations observed in hormone concentration
and potential dysfunction within the HPA axis in stressful situation. This
model integrates impact of hippocampal receptors on the secretion of
corticotropin-releasing hormone (CRH), an additional signaling pathway
involving Arginine Vasopressin (AVP) for the production and secretion of
adrenocorticotropic hormone (ACTH), the inclusion of a daylight-related
function for modeling circadian rhythms, and a short negative feedback loop
from the pituitary to the hypothalamus in a minimal mechanistic model of the
HPA axis. This expansion allows us to estimate model parameters that led to a
significant reduction in the mean absolute percent error, thereby enhancing the
model's predictive accuracy with a demonstration of a strong fit to the
validation dataset. Through sensitivity and correlation analyses, our study
shows the parameters that exert the most significant influence on the dynamics
of cortisol within the system. This study reveals intricate interdependencies
within the model and among the various components of the HPA axis. Despite
current advancements in the model to comprehend the dynamics of the HPA Axis,
these models are still lacking in addressing major factors such as the impact
of genetic and epigenetics, the role of the amygdala in processing and
transmitting stress signals, and other minor impacting factors.


  Simple and low-cost biosensing solutions are suitable for point-of-care
applications aiming to overcome the gap between scientific concepts and
technological production. To compete with sensitivity and selectivity of golden
standards, such as liquid chromatography, the functionalization of biosensors
is continuously optimized to enhance the signal and improve their performance,
often leading to complex chemical assay development. In this research, the
efforts are made on optimizing the methodology for electrochemical reduction of
graphene oxide to produce thin film-modified gold electrodes. Under the
employed specific conditions, 20 cycles of cyclic voltammetry (CV) are shown to
be optimal for superior electrical activation of graphene oxide into
electrochemically reduced graphene oxide (ERGO). This platform is further used
to develop a matrix metalloproteinase 2 (MMP-2) biosensor, where specific
anti-MMP2 aptamers are utilized as a biorecognition element. MMP-2 is a protein
which is typically overexpressed in tumor tissues, with important roles in
tumor invasion, metastasis as well as in tumor angiogenesis. Based on
impedimetric measurements, we were able to detect as low as 3.32 pg/mL of MMP-2
in PBS with a dynamic range of 10 pg/mL - 10 ng/mL. Besides high specificity,
ERGO-based aptasensor showed a potential of reuse due to demonstrated
successful signal restoration after experimental detection of MMP-2.


  Many imaging techniques for biological systems -- like fixation of cells
coupled with fluorescence microscopy -- provide sharp spatial resolution in
reporting locations of individuals at a single moment in time but also destroy
the dynamics they intend to capture. These snapshot observations contain no
information about individual trajectories, but still encode information about
movement and demographic dynamics, especially when combined with a
well-motivated biophysical model. The relationship between spatially evolving
populations and single-moment representations of their collective locations is
well-established with partial differential equations (PDEs) and their inverse
problems. However, experimental data is commonly a set of locations whose
number is insufficient to approximate a continuous-in-space PDE solution. Here,
motivated by popular subcellular imaging data of gene expression, we embrace
the stochastic nature of the data and investigate the mathematical foundations
of parametrically inferring demographic rates from snapshots of particles
undergoing birth, diffusion, and death in a nuclear or cellular domain. Toward
inference, we rigorously derive a connection between individual particle paths
and their presentation as a Poisson spatial process. Using this framework, we
investigate the properties of the resulting inverse problem and study factors
that affect quality of inference. One pervasive feature of this experimental
regime is the presence of cell-to-cell heterogeneity. Rather than being a
hindrance, we show that cell-to-cell geometric heterogeneity can increase the
quality of inference on dynamics for certain parameter regimes. Altogether, the
results serve as a basis for more detailed investigations of subcellular
spatial patterns of RNA molecules and other stochastically evolving populations
that can only be observed for single instants in their time evolution.


  Advanced omics technologies and facilities generate a wealth of valuable data
daily; however, the data often lacks the essential metadata required for
researchers to find and search them effectively. The lack of metadata poses a
significant challenge in the utilization of these datasets. Machine
learning-based metadata extraction techniques have emerged as a potentially
viable approach to automatically annotating scientific datasets with the
metadata necessary for enabling effective search. Text labeling, usually
performed manually, plays a crucial role in validating machine-extracted
metadata. However, manual labeling is time-consuming; thus, there is an need to
develop automated text labeling techniques in order to accelerate the process
of scientific innovation. This need is particularly urgent in fields such as
environmental genomics and microbiome science, which have historically received
less attention in terms of metadata curation and creation of gold-standard text
mining datasets.
  In this paper, we present two novel automated text labeling approaches for
the validation of ML-generated metadata for unlabeled texts, with specific
applications in environmental genomics. Our techniques show the potential of
two new ways to leverage existing information about the unlabeled texts and the
scientific domain. The first technique exploits relationships between different
types of data sources related to the same research study, such as publications
and proposals. The second technique takes advantage of domain-specific
controlled vocabularies or ontologies. In this paper, we detail applying these
approaches for ML-generated metadata validation. Our results show that the
proposed label assignment approaches can generate both generic and
highly-specific text labels for the unlabeled texts, with up to 44% of the
labels matching with those suggested by a ML keyword extraction algorithm.


  Phosphatidylserine (PS) is a negatively-charged glycerophospholipid found
mainly in the plasma membrane (PM) and in the late secretory/endocytic
compartments, where it regulates cellular activity and can mediate apoptosis.
Export of PS from the endoplasmic reticulum, its site of synthesis, to other
compartments, and its transbilayer asymmetry must therefore be precisely
regulated. We review recent findings on non-vesicular transport of PS by lipid
transfer proteins (LTPs) at membrane contact sites, on PS flip-flop between
membrane leaflets by flippases and scramblases, and on PS nano-clustering at
the PM. We also discuss emerging data on cooperation between scramblases and
LTPs, how perturbation of PS distribution can lead to disease, and the specific
role of PS in viral infection.


  Disease gene prioritization assigns scores to genes or proteins according to
their likely relevance for a given disease based on a provided set of seed
genes. Here, we describe a new algorithm for disease gene prioritization based
on continuous-time quantum walks using the adjacency matrix of a
protein-protein interaction (PPI) network. Our algorithm can be seen as a
quantum version of a previous method known as the diffusion kernel, but,
importantly, has higher performance in predicting disease genes, and also
permits the encoding of seed node self-loops into the underlying Hamiltonian,
which offers yet another boost in performance. We demonstrate the success of
our proposed method by comparing it to several well-known gene prioritization
methods on three disease sets, across seven different PPI networks. In order to
compare these methods, we use cross-validation and examine the mean reciprocal
ranks and recall values. We further validate our method by performing an
enrichment analysis of the predicted genes for coronary artery disease. We also
investigate the impact of adding self-loops to the seeds, and argue that they
allow the quantum walker to remain more local to low-degree seed nodes.


  Proteins populate a manifold in the high-dimensional sequence space whose
geometrical structure guides their natural evolution. Leveraging
recently-developed structure prediction tools based on transformer models, we
first examine the protein sequence landscape as defined by the folding score
function. This landscape shares characteristics with optimization challenges
encountered in machine learning and constraint satisfaction problems. Our
analysis reveals that natural proteins predominantly reside in wide, flat
minima within this energy landscape. To investigate further, we employ
statistical mechanics algorithms specifically designed to explore regions with
high local entropy in relatively flat landscapes. Our findings indicate that
these specialized algorithms can identify valleys with higher entropy compared
to those found using traditional methods such as Monte Carlo Markov Chains. In
a proof-of-concept case, we find that these highly entropic minima exhibit
significant similarities to natural sequences, especially in critical key sites
and local entropy. Additionally, evaluations through Molecular Dynamics
suggests that the stability of these sequences closely resembles that of
natural proteins. Our tool combines advancements in machine learning and
statistical physics, providing new insights into the exploration of sequence
landscapes where wide, flat minima coexist alongside a majority of narrower
minima.


  Computer-assisted methods have emerged as valuable tools for retrosynthesis
analysis. However, quantifying the plausibility of generated retrosynthesis
routes remains a challenging task. We introduce Retro-BLEU, a statistical
metric adapted from the well-established BLEU score in machine translation, to
evaluate the plausibility of retrosynthesis routes based on reaction template
sequences analysis. We demonstrate the effectiveness of Retro-BLEU by applying
it to a diverse set of retrosynthesis routes generated by state-of-the-art
algorithms and compare the performance with other evaluation metrics. The
results show that Retro-BLEU is capable of differentiating between plausible
and implausible routes. Furthermore, we provide insights into the strengths and
weaknesses of Retro-BLEU, paving the way for future developments and
improvements in this field.


  To end the pandemic caused by the coronavirus disease 2019 (COVID-19), many
countries have implemented large scale vaccinations for their nationals since
late 2020. Most of the vaccines have been allocated to wealthy countries while
few have been sent to countries with lower average income. Unequal distribution
of vaccines not only deepens the gap between the rich and poor but also
elongates the pandemic, as evidenced by emergent pathogenic variants of the
virus. We study the relationship between the pandemic ending time and vaccine
allocation time between two countries with mutual migration of peoples. Using
two coupled Susceptible-Infected-Recovered-Deceased (SIRD) model with
vaccinations, we show that sharing vaccines in time is beneficial to both
countries with and without the threat of virus variants. This shows that in the
matter of vaccine distribution the self-interest and altruism are not mutually
exclusive.


  Identifying molecular signatures from complex disease patients with
underlying symptomatic similarities is a significant challenge in analysis of
high dimensional multi-omics data. Topological data analysis (TDA) provides a
way of extracting such information from the geometric structure of the data and
identify multiway higher-order relationships. Here, we propose an application
of Harmonic persistent homology which overcomes the limitations of ambiguous
assignment of the topological information to the original elements in a
representative topological cycle from the data. When applied to multi-omics
data, this leads to the discovery of hidden patterns highlighting the
relationships between different omic profiles, while allowing for common tasks
in multiomics analyses such as disease subtyping, and most importantly
biomarker identification for similar latent biological pathways that are
associated with complex diseases. Our experiments on multiple cancer data show
that harmonic persistent homology and TDA can be very useful in dissecting
muti-omics data and identify biomarkers while detecting representative cycles
of the data which also predicts disease subtypes.