  In this paper, we consider the replicator-mutator dynamics for pairwise
social dilemmas where the payoff entries are random variables. The randomness
is incorporated to take into account the uncertainty, which is inevitable in
practical applications and may arise from different sources such as lack of
data for measuring the outcomes, noisy and rapidly changing environments, as
well as unavoidable human estimate errors. We analytically and numerically
compute the probability that the replicator-mutator dynamics has a given number
of equilibria for four classes of pairwise social dilemmas (Prisoner's Dilemma,
Snow-Drift Game, Stag-Hunt Game and Harmony Game). As a result, we characterise
the qualitative behaviour of such probabilities as a function of the mutation
rate. Our results clearly show the influence of the mutation rate and the
uncertainty in the payoff matrix definition on the number of equilibria in
these games. Overall, our analysis has provided novel theoretical contributions
to the understanding of the impact of uncertainty on the behavioural diversity
in a complex dynamical system.


  We analyze the management of the Italian pandemic during the five identified
waves. We considered the following problems: (i) The composition of the CTS
("Scientific Technical Committee"), which was composed entirely of doctors,
mainly virologists, without mathematical epidemiologists, statisticians,
physicists, etc. In fact, a pandemic has a behavior described by mathematical,
stochastic and probabilistic criteria; (ii) Political interference in security
measures and media propaganda; (iii) The initial stages of the vaccination
campaign, ignoring the age factor, and (iv) The persistence of the pandemic due
to the population unvaccinated (anti-vax or "no-vax"), which amounts to about
six to seven million people, including 10% of anti-vax doctors.


  In the field of functional genomics, the analysis of gene expression profiles
through Machine and Deep Learning is increasingly providing meaningful insight
into a number of diseases. The paper proposes a novel algorithm to perform
Feature Selection on genomic-scale data, which exploits the reconstruction
capabilities of autoencoders and an ad-hoc defined Explainable Artificial
Intelligence-based score in order to select the most informative genes for
diagnosis, prognosis, and precision medicine. Results of the application on a
Chronic Lymphocytic Leukemia dataset evidence the effectiveness of the
algorithm, by identifying and suggesting a set of meaningful genes for further
medical investigation.


  Formation of base pairs between the nucleotides of an RNA sequence gives rise
to a complex and often highly branched RNA structure. While numerous studies
have demonstrated the functional importance of the high degree of RNA branching
-- for instance, for its spatial compactness or interaction with other
biological macromolecules -- RNA branching topology remains largely unexplored.
Here, we use the theory of randomly branching polymers to explore the scaling
properties of RNAs by mapping their secondary structures onto planar tree
graphs. Focusing on random RNA sequences of varying lengths, we determine the
two scaling exponents related to their topology of branching. Our results
indicate that ensembles of RNA secondary structures are characterized by
annealed random branching and scale similarly to self-avoiding trees in three
dimensions. We further show that the obtained scaling exponents are robust upon
changes in nucleotide composition, tree topology, and folding energy
parameters. Finally, in order to apply the theory of branching polymers to
biological RNAs, whose length cannot be arbitrarily varied, we demonstrate how
both scaling exponents can be obtained from the distributions of the related
topological quantities of individual RNA molecules with fixed length. In this
way, we establish a framework to study the branching properties of RNA and
compare them to other known classes of branched polymers. By understanding the
scaling properties of RNA related to its branching structure we aim to improve
our understanding of the underlying principles and open up the possibility to
design RNA sequences with desired topological properties.


  The vast amount of sequencing data presently available allow the scientific
community to explore a range of genetic variables that may drive and progress
cancer. A myriad of predictive tools has been proposed, allowing researchers
and clinicians to compare and prioritize driver genes and mutations and their
relative pathogenicity. However, there is little consensus on the computational
approach or a golden standard for comparison. Hence, benchmarking the different
tools depends highly on the input data, indicating that overfitting is still a
massive problem. One of the solutions is to limit the scope and usage of
specific tool. However, such limitations forces researchers to walk on a
tightrope between creating and using high-quality tools for a specific purpose
and describing the complex alterations driving cancer. While the knowledge of
cancer development increases every day, many bioinformatic pipelines rely on
single nucleotide variants or alterations in a vacuum without accounting for
cellular compartment, mutational burden, or disease progression. Even within
bioinformatics and computational cancer biology, the research fields work in
silos, risking overlooking potential synergies or breakthroughs. Here, we
provide an overview of databases and datasets for building or testing
predictive tools for discovery of cancer drivers. We introduce predictive tools
for driver genes, driver mutations, and the impact of these based on structural
analysis. Additionally, we suggest and recommend directions in the field to
avoid silo-research, moving in the direction of integrative frameworks.


  Stochastic processes on graphs can describe a great variety of phenomena
ranging from neural activity to epidemic spreading. While many existing methods
can accurately describe typical realizations of such processes, computing
properties of extremely rare events is a hard task. Particularly so in the case
of recurrent models, in which variables may return to a previously visited
state. Here, we build on the matrix product cavity method, extending it
fundamentally in two directions: first, we show how it can be applied to Markov
processes biased by arbitrary reweighting factors that concentrate most of the
probability mass on rare events. Second, we introduce an efficient scheme to
reduce the computational cost of a single node update from exponential to
polynomial in the node degree. Two applications are considered: inference of
infection probabilities from sparse observations within the SIRS epidemic
model, and the computation of both typical observables and large deviations of
several kinetic Ising models.


  Deep learning-based approaches for generating novel drug molecules with
specific properties have gained a lot of interest in the last years. Recent
studies have demonstrated promising performance for string-based generation of
novel molecules utilizing reinforcement learning. In this paper, we develop a
unified framework for using reinforcement learning for de novo drug design,
wherein we systematically study various on- and off-policy reinforcement
learning algorithms and replay buffers to learn an RNN-based policy to generate
novel molecules predicted to be active against the dopamine receptor DRD2. Our
findings suggest that it is advantageous to use at least both top-scoring and
low-scoring molecules for updating the policy when structural diversity is
essential. Using all generated molecules at an iteration seems to enhance
performance stability for on-policy algorithms. In addition, when replaying
high, intermediate, and low-scoring molecules, off-policy algorithms display
the potential of improving the structural diversity and number of active
molecules generated, but possibly at the cost of a longer exploration phase.
Our work provides an open-source framework enabling researchers to investigate
various reinforcement learning methods for de novo drug design.


  Motile organisms can form stable agglomerates such as cities or colonies. In
the outbreak of a highly contagious disease, the control of large-scale
epidemic spread depends on factors like the number and size of agglomerates,
travel rate between them, and disease recovery rate. While the emergence of
agglomerates permits early interventions, it also explains longer real
epidemics. In this work, we study the spread of susceptible-infected-recovered
epidemics in one-dimensional spatially-structured systems. By working in one
dimension, we mimic microorganisms in narrow channels and establish a necessary
foundation for future investigation in higher dimensions. We employ a model of
self-propelled particles which spontaneously form multiple clusters. As the
rate of stochastic reorientation decreases, clusters become larger and less
numerous. Besides examining the time evolution averaged over many epidemics, we
show how the final number of ever-infected individuals depends non-trivially on
single-individual parameters. In particular, the number of ever-infected
individuals first increases with the reorientation rate since particles escape
sooner from clusters and spread the disease. For higher reorientation rate,
travel between clusters becomes too diffusive and the clusters too small,
decreasing the number of ever-infected individuals.


  The cooperative management of rice terraces in Bali reveals an interesting
phenomenon that stems from the feedback loop between human decisions and the
ecosystem process. In particular, spatial patterning is observed, which is
heavily reliant on the farmer's decision to plant crops as well as the response
from the physical environment like pest damage and water shortage. A recent
study proposed an evolutionary game theoretic model to infer particular power
laws governing this spatial patterning along the Bali region. In this paper, we
show a snapshot of rice patches in Bali with colors to indicate the different
stages of rice growth


  Invasive non-Typhoidal Salmonella (iNTS) is one of the leading causes of
blood stream infections in Sub-Saharan Africa, especially among children. iNTS
can be difficult to diagnose, particularly in areas where malaria is endemic
and difficult to treat, partly because of the emergence of antibiotic
resistance. We developed a mathematical model to evaluate the impact of a
vaccine for iNTS in 49 countries of sub-Saharan Africa. Without vaccination we
estimate 9.2 million new iNTS cases among children below 5 years old in these
49 countries from 2022 to 2038, 6.2 million of which between 2028 and 2038. The
introduction of a 85% (95%) efficacy vaccine in 2028 would prevent 2.6 (2.9)
million of these new infections. We provide the country-specific impact of a
iNTS vaccine considering the different age structures and vaccine coverage
levels.


  This paper provides a first attempt to incorporate the massive discontinuous
changes in the spatio-temporal dynamics of epidemics. Namely, we propose an
extended class of epidemic models, governed by coupled stochastic semilinear
partial differential equations, driven by pure-jump L\'evy noise. Based on the
considered type of incidence functions, by virtue of semi-group theory, a
truncation technique and Banach fixed point theorem, we prove the existence and
pathwise uniqueness of mild solutions, depending continuously on the initial
datum. Moreover, by means of a regularization technique, based on the resolvent
operator, we acquire that mild solutions can be approximated by a suitable
converging sequence of strong solutions. With this result at hand, for positive
initial states, we derive the almost-sure positiveness of the obtained
solutions. Finally, we present the outcome of several numerical simulations, in
order to exhibit the effect of the considered type of stochastic noise, in
comparison to Gaussian noise, which has been used in the previous literature.
Our established results lay the ground-work for investigating other problems
associated with the new proposed class of epidemic models, such as asymptotic
behavior analyses, optimal control as well as identification problems, which
primarily rely on the existence and uniqueness of biologically feasible
solutions.


  Cryogenic electron microscopy is widely used in structural biology, but its
resolution is often limited by the dynamics of the macromolecule. Here, we
developed a refinement protocol based on Gaussian mixture models that
integrates particle orientation and conformation estimation, and improves the
alignment for flexible domains of protein structures. We demonstrated this
protocol on multiple datasets, resulting in improved resolution and
resolvability, locally and globally, by visual and quantitative measures.


  A compartmental deterministic model that allows (1) immunity from two stages
of infection and carriage, and (2) disease induced death, is used in studying
the dynamics of meningitis epidemic process in a closed population. It allows
for difference in the transmission rate of infection to a susceptible by a
carrier and an infective. It is generalized to allow a proportion ({\phi}) of
those susceptibles infected to progress directly to infectives in stage I. Both
models are used in this study. The threshold conditions for the spread of
carrier and infectives in stage I are derived for the two models. Sensitivity
analysis is performed on the reproductive number derived from the next
generation matrix. The case-carrier ratio profile for various parameters and
threshold values are shown. So also are the graphs of the total number ever
infected as influenced by {\epsilon} and {\phi}. The infection transmission
rate (\b{eta}), the odds in favor of a carrier, over an infective, in
transmitting an infection to a susceptible ({\epsilon}) and the carrier
conversion rate ({\phi}) to an infective in stage I, are identified as key
parameters that should be subject of attention for any control intervention
strategy. The case-carrier ratio profiles provide evidence of a critical
case-carrier ratio attained before the number of reported cases grows to an
epidemic level. They also provide visual evidence of epidemiological context,
in this case, epidemic incidence (in later part of dry season) and endemic
incidence (during rainy season). Results from total proportion ever infected
suggest that the model, in which {\phi}=0 obtained, can adequately represent,
in essence, the generalized model for this study.


  Population structure has been known to substantially affect evolutionary
dynamics. Networks that promote the spreading of fitter mutants are called
amplifiers of natural selection, and those that suppress the spreading of
fitter mutants are called suppressors. Research in the past two decades has
found various families of amplifiers while suppressors still remain somewhat
elusive. It has also been discovered that most networks are amplifiers under
the birth-death updating combined with uniform initialization, which is a
standard condition assumed widely in the literature. In the present study, we
extend the birth-death processes to temporal (i.e., time-varying) networks. For
the sake of tractability, we restrict ourselves to switching temporal networks,
in which the network structure alternates between two static networks at
constant time intervals. We show that, in a majority of cases, switching
networks are less amplifying than both of the two static networks constituting
the switching networks. Furthermore, most small switching networks are
suppressors, which contrasts to the case of static networks.


  In this paper, the replicator dynamics of the two-locus two-allele system
under weak mutation and weak selection is investigated in a generation-wise
non-overlapping unstructured population of individuals mating at random. Our
main finding is that the dynamics is gradient-like when the point mutations at
the two loci are independent. This is in stark contrast to the case of
one-locus multi-allele where the existence gradient behaviour is contingent on
a specific relationship between the mutation rates. When the mutations are not
independent in the two-locus two-allele system, there is the possibility of
non-convergent outcomes, like asymptotically stable oscillations, through
either the Hopf bifurcation or the Neimark--Sacker bifurcation depending on the
strength of the weak selection. The results can be straightforwardly extended
for multi-locus two-allele systems.


  Representation learning is an important step in the machine learning
pipeline. Given the current biological sequencing data volume, learning an
explicit representation is prohibitive due to the dimensionality of the
resulting feature vectors. Kernel-based methods, e.g., SVM, are a proven
efficient and useful alternative for several machine learning (ML) tasks such
as sequence classification. Three challenges with kernel methods are (i) the
computation time, (ii) the memory usage (storing an $n\times n$ matrix), and
(iii) the usage of kernel matrices limited to kernel-based ML methods
(difficult to generalize on non-kernel classifiers). While (i) can be solved
using approximate methods, challenge (ii) remains for typical kernel methods.
Similarly, although non-kernel-based ML methods can be applied to kernel
matrices by extracting principal components (kernel PCA), it may result in
information loss, while being computationally expensive. In this paper, we
propose a general-purpose representation learning approach that embodies kernel
methods' qualities while avoiding computation, memory, and generalizability
challenges. This involves computing a low-dimensional embedding of each
sequence, using random projections of its $k$-mer frequency vectors,
significantly reducing the computation needed to compute the dot product and
the memory needed to store the resulting representation. Our proposed fast and
alignment-free embedding method can be used as input to any distance (e.g., $k$
nearest neighbors) and non-distance (e.g., decision tree) based ML method for
classification and clustering tasks. Using different forms of biological
sequences as input, we perform a variety of real-world classification tasks,
such as SARS-CoV-2 lineage and gene family classification, outperforming
several state-of-the-art embedding and kernel methods in predictive
performance.


  One of the main reasons we have not been able to cure cancers is that drugs
select for drug-resistant cancer cells. Pest managers face similar challenges
with pesticides selecting for pesticide-resistant organisms. Lessons in pest
management have led to four heuristics that can be translated to controlling
cancers 1. limit use (of chemical controls or modes of action to the lowest
practical levels) 2. diversify use (of modes of action largely through
rotations of chemical controls) 3. partition chemistry (modes of action through
space and time, which in effect is a refuge management strategy) and 4. include
non-chemical methods. These principles are general to all cancers and cancer
drugs, and thus should be employed to improve oncology. We review the parallel
difficulties in controlling the evolution of drug resistance in pests and
cancer cells, and describe the results of single- and multi-drug strategies in
agriculture and oncology. We dissect the methods that pest managers use to
prevent the evolution of pesticide resistance, showing how integrated pest
management inspired the development of adaptive therapy in oncology to
stabilize tumor size, and increase progression-free survival and quality of
life of patients. Finally, we demonstrate these principles in a proposal for
clinical trials in colorectal cancer.


  We revisit the original hawk-dove game with slight modifications to payoff
values while maintaining the fundamental principles of interaction. The
practical robustness of the theoretical tools of game theory is tested on a
simulated population of hawks and doves with varying initial population
distributions and peak growth rates. Additionally, we aim to find conditions in
which the entire community fails or becomes a single-species population. The
results show that the predicted community distribution is established by the
majority of communities but fails to exist in communities with extreme initial
imbalances in species distribution and insufficient growth rates. We also find
that greater growth rates can compensate for more imbalanced initial conditions
and that more balanced initial conditions can compensate for lower growth
rates. Overall, the simple theoretical model is a strong predictor of the
stable behavior of simulated multi-species communities.


  While mathematical models, in particular self-propelled particle (SPP)
models, capture many of the observed properties of large fish schools, they do
not always capture the interactions of smaller shoals. Nor do these models tend
to account for the observation that, when swimming alone or in smaller groups,
many species of fish use intermittent locomotion, often referred to as burst
and coast or burst and glide. Recent empirical studies have suggested that
burst and glide movement is indeed pivotal to the social interactions of
individual fish. In this paper, we propose a model of social burst and glide
motion by combining a well-studied model of neuronal dynamics, the
FitzHugh-Nagumo model, with a model of fish motion. We begin by showing that
the model can capture the motion of a single fish swimming down a channel. By
then extending to a two fish model, where visual stimuli of the position of the
other fish affect the internal burst or glide state of the fish, we find that
our model captures a rich set of swimming dynamics found in many species of
fish. These include: leader-follower behaviour; periodic changes in leadership;
apparently random (i.e. chaotic) leadership change; and pendulum-like
tit-for-tat turn taking. Unlike SPP models, which assume that fish move at a
constant speed, the model produces realistic motion of individual fish.
Moreover, unlike previous studies where a random component is used for
leadership switching to occur, we show that leadership switching, both periodic
and chaotic, can be the result of a deterministic interaction. We give several
empirically testable predictions on how fish interact and discuss our results
in light of recently established correlations between fish locomotion and brain
activity.


  Here we introduce simple structures for the analysis of complex hypergraphs,
hypergraph animals. These structures are designed to describe the local node
neighbourhoods of nodes in hypergraphs. We establish their relationships to
lattice animals and network motifs and develop their combinatorial properties,
for sparse and uncorrelated hypergraphs. Here we can make use of the tight link
of hypergraph animals to partition numbers, which opens up a vast mathematical
framework for the analysis of hypergraph animals. We then study their
abundances in random hypergraphs. Two transferable insights result from this
analysis: (i) it establishes the importance of high-cardinality edges in
ensembles of random hypergraphs that are inspired by the classical
Erdos-Reny\'i random graphs; and (ii) there is a close connection between
degree and hyperedge cardinality in random hypergraphs that shapes the animal
abundances and spectra profoundly. Together these two findings imply that we
need to spend more effort on investigating and developing suitable conditional
ensembles of random hypergraph that can capture real-world structures and their
complex dependency structures.


  The quantitative structure-activity relationship (QSAR) regression model is a
commonly used technique for predicting biological activities of compounds using
their molecular descriptors. Predictions from QSAR models can help, for
example, to optimize molecular structure; prioritize compounds for further
experimental testing; and estimate their toxicity. In addition to the accurate
estimation of the activity, it is highly desirable to obtain some estimate of
the uncertainty associated with the prediction, e.g., calculate a prediction
interval (PI) containing the true molecular activity with a pre-specified
probability, say 70%, 90% or 95%. The challenge is that most machine learning
(ML) algorithms that achieve superior predictive performance require some
add-on methods for estimating uncertainty of their prediction. The development
of these algorithms is an active area of research by statistical and ML
communities but their implementation for QSAR modeling remains limited.
Conformal prediction (CP) is a promising approach. It is agnostic to the
prediction algorithm and can produce valid prediction intervals under some weak
assumptions on the data distribution. We proposed computationally efficient CP
algorithms tailored to the most advanced ML models, including Deep Neural
Networks and Gradient Boosting Machines. The validity and efficiency of
proposed conformal predictors are demonstrated on a diverse collection of QSAR
datasets as well as simulation studies.


  Crocus sativus is a triploid sterile plant with red stigmas belonging to the
family of Iridaceae, and sub-family Crocoideae. Crocin, picrocrocin, and
safranal are three major carotenoid derivatives that are responsible for the
color, taste, and specific aroma of Crocus. Saffron flowers are harvested
manually and used as spice, dye, or medicinal applications. The natural
propagation rate of most geophytes including saffron is relatively low. An in
vitro multiplication technique like micropropagation has been used for the
propagation of saffron. To understand the efficiency of this alternative and
study the molecular basis of apocarotenoid biosynthesis/accumulation, the
RT-PCR method was performed on perianth explants that were cultured on MS
medium to observe the level of expression of zeaxanthin cleavage dioxygenase
(CsZCD) gene during stigma development, and also the impact of soil EC on its
expression. The present study was conducted at Plant molecular and physiology
Lab, Alzahra University, Tehran, Iran during 2011-2013. Stigma-like structures
(SLSs) on calli were collected from immature perianth explants from floral buds
of corms that were collected from Ghaen city, and compared to (Torbat-e
Haidariye, Mardabad, and Shahroud cities) for investigating the impact of
different soil EC on CsZCD expression. The results indicated that the CsZCD
gene was highly expressed in fully developed red SLSs in perianth of cultured
samples of Shahroud with the highest salinity. In this research, a close
relationship between soil EC and second metabolites regulation is studied.
Overall, these results will pave the way for understanding the molecular basis
of apocarotenoid biosynthesis and other aspects of stigma development in C.
sativus.


  We consider the effect of network structure on the evolution of a population.
Models of this kind typically consider a population of fixed size and
distribution. Here we consider eco-evolutionary dynamics where population size
and distribution can change through birth, death and migration, all of which
are separate processes. This allows complex interaction and migration
behaviours that are dependent on competition. For migration, we assume that the
response of individuals to competition is governed by tolerance to their group
members, such that less tolerant individuals are more likely to move away due
to competition. We looked at the success of a mutant in the rare mutation limit
for the complete, cycle and star networks. Unlike models with fixed population
size and distribution, the distribution of the individuals per site is
explicitly modelled by considering the dynamics of the population. This in turn
determines the mutant appearance distribution for each network. Where a mutant
appears impacts its success as it determines the competition it faces. For low
and high migration rates the complete and cycle networks have similar mutant
appearance distributions resulting in similar success levels for an invading
mutant. A higher migration rate in the star network is detrimental for mutant
success because migration results in a crowded central site where a mutant is
more likely to appear.


  Tuberculosis (TB) is among the main public health challenges in Burundi. The
literature lacks mathematical models for key parameter estimates of TB
transmission dynamics in Burundi. In this paper, the
supectible-exposed-infected-recovered (SEIR) model is used to investigate the
transmission dynamics of tuberculosis in Burundi. Using the next generation
method, we calculated the basic reproduction number, R0. The model is
demonstrated to have a disease-free equilibrium (DEF) that is locally and
globally asymptotically stable. When the corresponding reproduction threshold
quantity approaches unity, the model enters an endemic equilibrium (EE). That
means, the disease can be controlled through different interventions in
Burundi. A sensitivity analysis of the model parameters was also investigated.
It shows that the progression rate from latent to becoming infectious had the
highest positive sensitivity, which means that R0 increases and decreases
proportionally with an increase and a decrease of that progression rate.


  Protein structure prediction has reached revolutionary levels of accuracy on
single structures, yet distributional modeling paradigms are needed to capture
the conformational ensembles and flexibility that underlie biological function.
Towards this goal, we develop EigenFold, a diffusion generative modeling
framework for sampling a distribution of structures from a given protein
sequence. We define a diffusion process that models the structure as a system
of harmonic oscillators and which naturally induces a cascading-resolution
generative process along the eigenmodes of the system. On recent CAMEO targets,
EigenFold achieves a median TMScore of 0.84, while providing a more
comprehensive picture of model uncertainty via the ensemble of sampled
structures relative to existing methods. We then assess EigenFold's ability to
model and predict conformational heterogeneity for fold-switching proteins and
ligand-induced conformational change. Code is available at
https://github.com/bjing2016/EigenFold.


  We study the generic dimension of the zero set over $\mathbb{C}^*$,
$\mathbb{R}^*$, and $\mathbb{R}_{>0}$ of parametric polynomial systems
consisting of linear combinations of monomials scaled by free parameters. These
systems generalize sparse systems with fixed monomial support and freely
varying parametric coefficients. As our main result, we establish the
equivalence among three key properties: the existence of nondegenerate zeros,
the zero set having the expected dimension, and the system having zeros for
parameters in a Zariski dense subset of parameter space. Systems of this form
describe the steady states of reaction networks modeled with mass-action
kinetics, and we use our results to answer several fundamental geometric
questions on topics such as generic finiteness, the difference between kinetic
and stoichiometric subspaces, absolute concentration robustness, and
nondegenerate multistationarity.


  In this article, discrete and stochastic changes in (effective) population
size are incorporated into the spectral representation of a biallelic diffusion
process for drift and small mutation rates. A forward algorithm inspired by
Hidden-Markov-Model (HMM) literature is used to compute exact sample allele
frequency spectra for three demographic scenarios: single changes in
(effective) population size, boom-bust dynamics, and stochastic fluctuations in
(effective) population size. An approach for fully agnostic demographic
inference from these sample allele spectra is explored, and sufficient
statistics for step-wise changes in population size are found. Further,
convergence behaviours of the polymorphic sample spectra for population size
changes on different time scales are examined and discussed within the context
of inference of the effective population size. Joint visual assessment of the
sample spectra and the temporal coefficients of the spectral decomposition of
the forward diffusion process is found to be important in determining departure
from equilibrium. Stochastic changes in (effective) population size are shown
to shape sample spectra particularly strongly.


  Biochemical covalent modification networks exhibit a remarkable suite of
steady state and dynamical properties such as multistationarity, oscillations,
ultrasensitivity and absolute concentration robustness. This paper focuses on
conditions required for a network to have a species with absolute concentration
robustness. We find that the robustness in a substrate is endowed by its
interaction with a bifunctional enzyme, which is an enzyme that has different
roles when isolated versus when bound as a substrate-enzyme complex. When
isolated, the bifunctional enzyme promotes production of more molecules of the
robust species while when bound, the same enzyme facilitates degradation of the
robust species. These dual actions produce robustness in the large class of
covalent modification networks. For each network of this type, we find the
network conditions for the presence of robustness, the species that has
robustness, and its robustness value. The unified approach of simultaneously
analyzing a large class of networks for a single property, i.e. absolute
concentration robustness, reveals the underlying mechanism of the action of
bifunctional enzyme while simultaneously providing a precise mathematical
description of bifunctionality.


  In recent years, single cell RNA sequencing has become a widely used
technique to study cellular diversity and function. However, accurately
annotating cell types from single cell data has been a challenging task, as it
requires extensive knowledge of cell biology and gene function. The emergence
of large language models such as ChatGPT and New Bing in 2023 has
revolutionized this process by integrating the scientific literature and
providing accurate annotations of cell types. This breakthrough enables
researchers to conduct literature reviews more efficiently and accurately, and
can potentially uncover new insights into cell type annotation. By using
ChatGPT to annotate single cell data, we can relate rare cell type to their
function and reveal specific differentiation trajectories of cell subtypes that
were previously overlooked. This can have important applications in
understanding cancer progression, mammalian development, and stem cell
differentiation, and can potentially lead to the discovery of key cells that
interrupt the differentiation pathway and solve key problems in the life
sciences. Overall, the future of cell type annotation in single cell data looks
promising and the Large Language model will be an important milestone in the
history of single cell analysis.


  Combining a set of phylogenetic trees into a single phylogenetic network that
explains all of them is a fundamental challenge in evolutionary studies.
Existing methods are computationally expensive and can either handle only small
numbers of phylogenetic trees or are limited to severely restricted classes of
networks. In this paper, we apply the recently-introduced theoretical framework
of cherry picking to design a class of efficient heuristics that are guaranteed
to produce a network containing each of the input trees, for datasets
consisting of binary trees. Some of the heuristics in this framework are based
on the design and training of a machine learning model that captures essential
information on the structure of the input trees and guides the algorithms
towards better solutions. We also propose simple and fast randomised heuristics
that prove to be very effective when run multiple times.
  Unlike the existing exact methods, our heuristics are applicable to datasets
of practical size, and the experimental study we conducted on both simulated
and real data shows that these solutions are qualitatively good, always within
some small constant factor from the optimum. Moreover, our machine-learned
heuristics are one of the first applications of machine learning to
phylogenetics and show its promise.


  The immune response to an acute primary infection is a coupled process of
antigen proliferation, molecular recognition by naive B cells, and their
subsequent proliferation and antibody shedding. This process contains a
fundamental problem: the recognition of an exponentially time-dependent antigen
signal. Here we show that B cells can efficiently recognise new antigens by a
tuned kinetic proofreading mechanism, where the molecular recognition machinery
is adapted to the complexity of the immune repertoire. This process produces
potent, specific and fast recognition of antigens, maintaining a spectrum of
genetically distinct B cell lineages as input for affinity maturation. We show
that the proliferation-recognition dynamics of a primary infection is a
generalised Luria-Delbr\"uck process, akin to the dynamics of the classic
fluctuation experiment. This map establishes a link between signal recognition
dynamics and evolution. We derive the resulting statistics of the activated
immune repertoire: antigen binding affinity, expected size, and frequency of
active B cell clones are related by power laws. Their exponents define the
class of generalised Luria-Delbr\"uck processes; they depend on the antigen and
B cell proliferation rate, the number of proofreading steps, and the lineage
density of the naive repertoire. Empirical data of mouse immune repertoires are
found to be consistent with activation involving at least three proofreading
steps. Our model predicts key clinical characteristics of acute infections,
including the emergence of elite neutralisers and the effects of immune ageing.
More broadly, our results establish acute infections as a new probe into the
global architecture and functional principles of immune repertoires.


  The amount of sequencing data for SARS-CoV-2 is several orders of magnitude
larger than any virus. This will continue to grow geometrically for SARS-CoV-2,
and other viruses, as many countries heavily finance genomic surveillance
efforts. Hence, we need methods for processing large amounts of sequence data
to allow for effective yet timely decision-making. Such data will come from
heterogeneous sources: aligned, unaligned, or even unassembled raw nucleotide
or amino acid sequencing reads pertaining to the whole genome or regions (e.g.,
spike) of interest. In this work, we propose \emph{ViralVectors}, a compact
feature vector generation from virome sequencing data that allows effective
downstream analysis. Such generation is based on \emph{minimizers}, a type of
lightweight "signature" of a sequence, used traditionally in assembly and read
mapping -- to our knowledge, the first use minimizers in this way. We validate
our approach on different types of sequencing data: (a) 2.5M SARS-CoV-2 spike
sequences (to show scalability); (b) 3K Coronaviridae spike sequences (to show
robustness to more genomic variability); and (c) 4K raw WGS reads sets taken
from nasal-swab PCR tests (to show the ability to process unassembled reads).
Our results show that ViralVectors outperforms current benchmarks in most
classification and clustering tasks.


  We study stochastic dynamics of the Prisoner's Dilemma game on random
Erd\"{o}s-R\'{e}nyi and scale-free Barab\'{a}si-Albert networks with a cost of
maintaining a link between interacting players. We show that when the cost
increases, the population of players located Barab\'{a}si-Albert network
undergoes a sharp transition from an ordered state, where almost all players
cooperate, to a state in which both cooperators and defectors coexist. At the
critical cost, the population oscillates in time between these two states.


  Prokaryotic gene expression is dynamic and noisy, and can lead to phenotypic
variations. Previous work has shown that the selection of such variation can be
modeled through basic reaction networks described by O.D.E. displaying a
bistable behavior. While previous mathematical studies have shown that mono- or
bistable behavior depends on the rate of the reactions in the system, no
analytical solution of the curves delineating the actual parameter conditions
that result in mono or bistability has so far been provided. In this work we
provide the first explicit analytical solution for the boundary curve that
separates the parameter space defining domains where double positive and double
negative feedback loops become bistable.


  We study the classic pure selection integrodifferential equation, stemming
from adaptative dynamics, in a measure framework by mean of duality approach.
After providing a well posedness result under fairly general assumptions, we
focus on the asymptotic behavior of various cases, illustrated by some
numerical simulations.


  We study a stochastic epidemic model with multiple patches (locations), where
individuals in each patch are categorized into three compartments, Susceptible,
Infected and Recovered/Removed, and may migrate from one patch to another in
any of the compartments. Each individual is associated with a random
infectivity function which dictates the force of infectivity while the
interactive infection process depends on the age of infection (elapsed time
since infection). We prove a functional law of large number for the epidemic
evolution dynamics including the aggregate infectivity process, the numbers of
susceptible and recovered individuals as well as the number of infected
individuals at each time that have been infected for a certain amount of time.
From the limits, we derive a PDE model for the density of the number of
infected individuals with respect to the infection age, which is a systems of
linear PDE equations with a boundary condition that is determined by a set of
integral equations.


  Causal discovery of genome-scale networks is important for identifying
pathways from genes to observable traits - e.g. differences in cell function,
disease, drug resistance and others. Causal learners based on graphical models
rely on interventional samples to orient edges in the network. However, these
models have not been shown to scale up the size of the genome, which are on the
order of 1e3-1e4 genes. We introduce a new learner, SP-GIES, that jointly
learns from interventional and observational datasets and achieves almost 4x
speedup against an existing learner for 1,000 node networks. SP-GIES achieves
an AUC-PR score of 0.91 on 1,000 node networks, and scales up to 2,000 node
networks - this is 4x larger than existing works. We also show how SP-GIES
improves downstream optimal experimental design strategies for selecting
interventional experiments to perform on the system. This is an important step
forward in realizing causal discovery at scale via autonomous experimental
design.


  Nuclear pore complexes (NPCs) mediate the exchange of materials between the
nucleoplasm and cytoplasm, playing a key role in the separation of nucleic
acids and proteins into their required compartments. The static structure of
the NPC is relatively well defined by recent cryo EM and other studies. The
functional roles of dynamic components in the pore of the NPC,
phenylalanyl-glycyl (FG) repeat rich nucleoporins, is less clear because of our
limited understanding of highly dynamic protein systems. These proteins form a
restrained concentrate which interacts with and concentrates nuclear transport
factors (NTRs) to provide facilitated nucleocytoplasmic transport of cargoes.
Very rapid exchange among FG repeats and NTRs supports extremely fast
facilitated transport, close to the rate of macromolecular diffusion in
cytoplasm, while complexes without specific interactions are entropically
excluded, though details on several aspects of the transport mechanism and FG
repeat behaviors remain to be resolved. However, as discussed here, new
technical approaches combined with more advanced modeling methods will likely
provide an improved dynamic description of NPC transport, potentially at the
atomic level in the near future. Such advances are likely to be of major
benefit in comprehending the roles the malfunctioning NPC plays in cancer,
aging, viral diseases, and neurodegeneration.


  Genome-Wide Association Studies (GWAS) help identify genetic variations in
people with diseases such as Parkinson's disease (PD), which are less common in
those without the disease. Thus, GWAS data can be used to identify genetic
variations associated with the disease. Feature selection and machine learning
approaches can be used to analyze GWAS data and identify potential disease
biomarkers. However, GWAS studies have technical variations that affect the
reproducibility of identified biomarkers, such as differences in genotyping
platforms and selection criteria for individuals to be genotyped. To address
this issue, we collected five GWAS datasets from the database of Genotypes and
Phenotypes (dbGaP) and explored several data integration strategies. We
evaluated the agreement among different strategies in terms of the Single
Nucleotide Polymorphisms (SNPs) that were identified as potential PD
biomarkers. Our results showed a low concordance of biomarkers discovered using
different datasets or integration strategies. However, we identified fifty SNPs
that were identified at least twice, which could potentially serve as novel PD
biomarkers. These SNPs are indirectly linked to PD in the literature but have
not been directly associated with PD before. These findings open up new
potential avenues of investigation.


  We propose a compartmental model for a disease with temporary immunity and
secondary infections. From our assumptions on the parameters involved in the
model, the system naturally evolves in three time scales. We characterize the
equilibria of the system and analyze their stability. We find conditions for
the existence of two endemic equilibria, for some cases in which $\mathcal{R}_0
< 1$. Then, we unravel the interplay of the three time scales, providing
conditions to foresee whether the system evolves in all three scales, or only
in the fast and the intermediate ones. We conclude with numerical simulations
and bifurcation analysis, to complement our analytical results.


  Understanding how proteins structurally interact is crucial to modern
biology, with applications in drug discovery and protein design. Recent machine
learning methods have formulated protein-small molecule docking as a generative
problem with significant performance boosts over both traditional and deep
learning baselines. In this work, we propose a similar approach for rigid
protein-protein docking: DiffDock-PP is a diffusion generative model that
learns to translate and rotate unbound protein structures into their bound
conformations. We achieve state-of-the-art performance on DIPS with a median
C-RMSD of 4.85, outperforming all considered baselines. Additionally,
DiffDock-PP is faster than all search-based methods and generates reliable
confidence estimates for its predictions. Our code is publicly available at
$\texttt{https://github.com/ketatam/DiffDock-PP}$


  In recent years, free energy perturbation (FEP) calculations have garnered
increasing attention as tools to support drug discovery. The lead optimization
mapper (Lomap) was proposed as an algorithm to calculate the relative free
energy between ligands efficiently. However, Lomap requires checking whether
each edge in the FEP graph is removable, which necessitates checking the
constraints for all edges. Consequently, conventional Lomap requires
significant computation time, at least several hours for cases involving
hundreds of compounds, and is impractical for cases with more than tens of
thousands of edges. In this study, we aimed to reduce the computational cost of
Lomap to enable the construction of FEP graphs for hundreds of compounds. We
can reduce the overall number of constraint checks required from an amount
dependent on the number of edges to one dependent on the number of nodes by
using the chunk check process to check the constraints for as many edges as
possible simultaneously. Moreover, the output graph is equivalent to that
obtained using conventional Lomap, enabling direct replacement of the original
Lomap with our method. With our improvement, the execution was tens to hundreds
of times faster than that of the original Lomap.
https://github.com/ohuelab/FastLomap


  Buffaloes are farm animals that contribute to food security by providing high
quality meat and milk. They can better tolerate the adverse effects of global
climate change on their meat and milk production. Despite their advantages,
buffaloes are heavily neglected animals with fewer studies compared to other
farm animals, hence, the real potential of buffaloes has never been realized.
The complete genome sequencing projects of buffaloes are essential to better
understanding the buffalos biology and production since they allow scientists
to identify important genes and understand how the gene networks interact to
determine the critical features of buffaloes. The genome projects are also
valuable for gaining better knowledge of growth, development, maintenance, and
determining factors associated with increased meat and milk production.
Furthermore, having access to a complete genome of high quality and
comprehensive annotations provides a powerful tool in breeding programs. The
current review surveyed the publicly available buffalo genome projects and
studied the impact of incorporating genomic selection into the buffalo breeding
program. Our survey of the publicly available buffalo genome projects showed
the promise of genomic selection in developing water buffalo science and
technology for food security on a global scale.


  Surface plasmon resonance (SPR)-based biosensors are widely used instruments
for characterizing molecular interactions. In theory the SPR signal depends
only on mass changes for interacting molecules of same chemical nature. Whether
conformational changes of interacting molecules also contribute to the SPR
signal is still a subject of lively debates. Works have been published claiming
that conformational changes were detected but all factors contributing to the
SPR signal were not carefully considered, in addition to often using no or
improper controls. In the present work we used a very well-characterized
oligonucleotide, the thrombin-binding DNA aptamer (TBA), which upon binding of
potassium ions folds into a two G-tetrad antiparallel G-quadruplex structure.
All terms contributing to the maximal expected SPR response, Rmax, in
particular the refractive index increment, RII, of both partners and the
fraction of immobilized TBA target available, ca, were experimentally assessed.
The resulting Rmax was then compared to the maximal experimental SPR response
for potassium ions binding to TBA using appropriate controls. Regardless how
the RIIs were measured, by SPR or refractometry, and how much TBA available for
interacting with potassium ions was considered, the theoretical and the
experimental SPR responses never matched, the former being always lower than
the latter. Using a straightforward experimental model system and by thoroughly
taking into account all contributing factors we therefore conclude that
conformational changes can indeed contribute to the measured SPR signal.


  In this proposed work, we identified the significant research issues on lung
cancer risk factors. Capturing and defining symptoms at an early stage is one
of the most difficult phases for patients. Based on the history of patients
records, we reviewed a number of current research studies on lung cancer and
its various stages. We identified that lung cancer is one of the significant
research issues in predicting the early stages of cancer disease. This research
aimed to develop a model that can detect lung cancer with a remarkably high
level of accuracy using the deep learning approach (convolution neural
network). This method considers and resolves significant gaps in previous
studies. We compare the accuracy levels and loss values of our model with
VGG16, InceptionV3, and Resnet50. We found that our model achieved an accuracy
of 94% and a minimum loss of 0.1%. Hence physicians can use our convolution
neural network models for predicting lung cancer risk factors in the real
world. Moreover, this investigation reveals that squamous cell carcinoma,
normal, adenocarcinoma, and large cell carcinoma are the most significant risk
factors. In addition, the remaining attributes are also crucial for achieving
the best performance.


  The telegraph model is the standard model of stochastic gene expression,
which can be solved exactly to obtain the distribution of mature RNA numbers
per cell. A modification of this model also leads to an analytical distribution
of the nascent RNA numbers. These solutions are routinely used for the analysis
of single-cell data, including the inference of transcriptional parameters.
However, these models neglect important mechanistic features of transcription
elongation, such as the stochastic movement of RNA polymerases and their steric
interactions. Here we construct a model of gene expression describing promoter
switching between inactive and active states, binding of RNA polymerases in the
active state, their stochastic movement including steric interactions along the
gene, and their unbinding leading to a mature transcript that subsequently
decays. We derive the steady-state distributions of the nascent and mature RNA
numbers in two important limiting cases: constitutive expression and slow
promoter switching. We show that RNA fluctuations are suppressed by steric
interactions between RNA polymerases, and that this suppression can even lead
to sub-Poissonian fluctuations; these effects are most pronounced for nascent
RNA and less prominent for mature RNA, since the latter is not a direct sensor
of transcription. We find a relationship between the parameters of our
microscopic mechanistic model and those of the standard models that ensures
excellent consistency in their prediction of the first and second RNA number
moments over vast regions of parameter space, encompassing slow, intermediate,
and rapid promoter switching, provided the RNA number distributions are
Poissonian or super-Poissonian. Furthermore, we identify the limitations of
inference from mature RNA data, specifically showing that it cannot
differentiate between highly distinct RNA polymerase traffic patterns on a
gene.


  All possible phenotypes are not equally accessible to evolving populations.
In fact, only phenotypes of large size, i.e. those resulting from many
different genotypes, are found in populations of sequences, presumably because
they are easier to discover and maintain. Genotypes that map to these
phenotypes usually form mostly connected genotype networks that percolate the
space of sequences, thus guaranteeing access to a large set of alternative
phenotypes. Within a given environment, where specific phenotypic traits become
relevant for adaptation, the replicative ability of a phenotype and its overall
fitness (in competition experiments with alternative phenotypes) can be
estimated. Two primary questions arise: how do phenotype size, reproductive
capability and topology of the genotype network affect the fitness of a
phenotype? And, assuming that evolution is only able to access large
phenotypes, what is the range of unattainable fitness values? In order to
address these questions, we quantify the adaptive advantage of phenotypes of
varying size and spectral radius in a two-peak landscape. We derive analytical
relationships between the three variables (size, topology, and replicative
ability) which are then tested through analysis of genotype-phenotype maps and
simulations of population dynamics on such maps. Finally, we analytically show
that the fraction of attainable phenotypes decreases with the length of the
genotype, though its absolute number increases. The fact that most phenotypes
are not visible to evolution very likely forbids the attainment of the highest
peak in the landscape. Nevertheless, our results indicate that the relative
fitness loss due to this limited accessibility is largely inconsequential for
adaptation.


  Networks are a convenient way to represent many interactions among different
entities as they provide an efficient and clear methodology to evaluate and
organize relevant data. While there are many features for characterizing
networks there is a quantity that seems rather elusive: Complexity. The
quantification of the complexity of networks is nowadays a fundamental problem.
Here, we present a novel tool for identifying the complexity of ecological
networks. We compare the behavior of two relevant indices of complexity:
K-complexity and Single value decomposition (SVD) entropy. For that, we use
real data and null models. Both null models consist of randomized networks
built by swapping a controlled number of links of the original ones. We analyze
23 plant-pollinator and 19 host-parasite networks as case studies. Our results
show interesting features in the behavior for the K-complexity and SVD entropy
with clear differences between pollinator-plant and host-parasite networks,
especially when the degree distribution is not preserved. Although SVD entropy
has been widely used to characterize network complexity, our analyses show that
K-complexity is a more reliable tool. Additionally, we show that degree
distribution and density are important drivers of network complexity and should
be accounted for in future studies.


  Integration of heterogeneous and high-dimensional multi-omics data is
becoming increasingly important in understanding genetic data. Each omics
technique only provides a limited view of the underlying biological process and
integrating heterogeneous omics layers simultaneously would lead to a more
comprehensive and detailed understanding of diseases and phenotypes. However,
one obstacle faced when performing multi-omics data integration is the
existence of unpaired multi-omics data due to instrument sensitivity and cost.
Studies may fail if certain aspects of the subjects are missing or incomplete.
In this paper, we propose a deep learning method for multi-omics integration
with incomplete data by Cross-omics Linked unified embedding with Contrastive
Learning and Self Attention (CLCLSA). Utilizing complete multi-omics data as
supervision, the model employs cross-omics autoencoders to learn the feature
representation across different types of biological data. The multi-omics
contrastive learning, which is used to maximize the mutual information between
different types of omics, is employed before latent feature concatenation. In
addition, the feature-level self-attention and omics-level self-attention are
employed to dynamically identify the most informative features for multi-omics
data integration. Extensive experiments were conducted on four public
multi-omics datasets. The experimental results indicated that the proposed
CLCLSA outperformed the state-of-the-art approaches for multi-omics data
classification using incomplete multi-omics data.


  DNA methylation is a significant driver of cell-type heterogeneity and has
been implicated in various regulatory processes ranging from cell
differentiation to imprinting. As the methyl group is embedded in the DNA
molecule, assessing DNA methylation is particularly promising in
liquid-biopsy-based approaches, as cell-free DNA retains information related to
its cell of origin. In this work, I leverage a recently profiled collection of
cell-sorted whole genome bisulfite profiles of 44 healthy cell types. The high
quality and purity of such data provide an ideal basis for discovering and
characterizing discriminative DNA methylation regions that could serve as a
reference for computational deconvolution. First, I characterize differentially
methylated regions between every pair of cell types, obtaining a meaningful
measure of divergence. Pairwise differences were then aggregated to identify a
set of uniquely (de)methylated regions (UMRs) for each cell type. Identified
UMRs are predominantly hypomethylated and their numbers vary significantly
across cell types. They are mostly located in enhancer regions and strongly
support cell-type-specific characteristics. As mapping onto UMRs has proven
unsuitable for deconvolution, I developed a novel approach utilizing the set
cover algorithm to select discriminative regions for this purpose. Based on
these regions, deconvolution was performed in two distinct approaches: a
beta-value-based and a read-level one. Both approaches outperform an existing
deconvolution software modeled on the same data 3-fold in terms of total
deconvolution error. Surprisingly, the beta-based approach slightly
outperformed the read-level one. Overall, I present an adaptable, end-to-end
software solution (methylcover) for obtaining accurate cell type deconvolution,
with possible future applications to non-invasive assays for disease detection
and monitoring.


  Accurately inferring Gene Regulatory Networks (GRNs) is a critical and
challenging task in biology. GRNs model the activatory and inhibitory
interactions between genes and are inherently causal in nature. To accurately
identify GRNs, perturbational data is required. However, most GRN discovery
methods only operate on observational data. Recent advances in neural
network-based causal discovery methods have significantly improved causal
discovery, including handling interventional data, improvements in performance
and scalability. However, applying state-of-the-art (SOTA) causal discovery
methods in biology poses challenges, such as noisy data and a large number of
samples. Thus, adapting the causal discovery methods is necessary to handle
these challenges. In this paper, we introduce DiscoGen, a neural network-based
GRN discovery method that can denoise gene expression measurements and handle
interventional data. We demonstrate that our model outperforms SOTA neural
network-based causal discovery methods.


  Direct reciprocity is a powerful mechanism for evolution of cooperation based
on repeated interactions between the same individuals. But high levels of
cooperation evolve only if the benefit-to-cost ratio exceeds a certain
threshold that depends on memory length. For the best-explored case of
one-round memory, that threshold is two. Here we report that intermediate
mutation rates lead to high levels of cooperation, even if the benefit-to-cost
ratio is only marginally above one, and even if individuals only use a minimum
of past information. This surprising observation is caused by two effects.
First, mutation generates diversity which undermines the evolutionary stability
of defectors. Second, mutation leads to diverse communities of cooperators that
are more resilient than homogeneous ones. This finding is relevant because many
real world opportunities for cooperation have small benefit-to-cost ratios,
which are between one and two, and we describe how direct reciprocity can
attain cooperation in such settings. Our result can be interpreted as showing
that diversity, rather than uniformity, promotes evolution of cooperation.


  Signal transduction and cell function are governed by the spatiotemporal
organization of membrane-associated molecules. Despite significant advances in
visualizing molecular distributions by 3D light microscopy, cell biologists
still have limited quantitative understanding of the processes implicated in
the regulation of molecular signals at the whole cell scale. In particular,
complex and transient cell surface morphologies challenge the complete sampling
of cell geometry, membrane-associated molecular concentration and activity and
the computing of meaningful parameters such as the cofluctuation between
morphology and signals. Here, we introduce u-Unwrap3D, a framework to remap
arbitrarily complex 3D cell surfaces and membrane-associated signals into
equivalent lower dimensional representations. The mappings are bidirectional,
allowing the application of image processing operations in the data
representation best suited for the task and to subsequently present the results
in any of the other representations, including the original 3D cell surface.
Leveraging this surface-guided computing paradigm, we track segmented surface
motifs in 2D to quantify the recruitment of Septin polymers by blebbing events;
we quantify actin enrichment in peripheral ruffles; and we measure the speed of
ruffle movement along topographically complex cell surfaces. Thus, u-Unwrap3D
provides access to spatiotemporal analyses of cell biological parameters on
unconstrained 3D surface geometries and signals.


  Graph neural networks (GNNs) demonstrate great performance in compound
property and activity prediction due to their capability to efficiently learn
complex molecular graph structures. However, two main limitations persist
including compound representation and model interpretability. While atom-level
molecular graph representations are commonly used because of their ability to
capture natural topology, they may not fully express important substructures or
functional groups which significantly influence molecular properties.
Consequently, recent research proposes alternative representations employing
reduction techniques to integrate higher-level information and leverages both
representations for model learning. However, there is still a lack of study
about different molecular graph representations on model learning and
interpretation. Interpretability is also crucial for drug discovery as it can
offer chemical insights and inspiration for optimization. Numerous studies
attempt to include model interpretation to explain the rationale behind
predictions, but most of them focus solely on individual prediction with little
analysis of the interpretation on different molecular graph representations.
This research introduces multiple molecular graph representations that
incorporate higher-level information and investigates their effects on model
learning and interpretation from diverse perspectives. The results indicate
that combining atom graph representation with reduced molecular graph
representation can yield promising model performance. Furthermore, the
interpretation results can provide significant features and potential
substructures consistently aligning with background knowledge. These multiple
molecular graph representations and interpretation analysis can bolster model
comprehension and facilitate relevant applications in drug discovery.


  In stopping the spread of infectious diseases, pathogen genomic data can be
used to reconstruct transmission events and characterize population-level
sources of infection. Most approaches for identifying transmission pairs do not
account for the time that passed since divergence of pathogen variants in
individuals, which is problematic in viruses with high within-host evolutionary
rates. This is prompting us to consider possible transmission pairs in terms of
phylogenetic data and additional estimates of time since infection derived from
clinical biomarkers. We develop Bayesian mixture models with an evolutionary
clock as signal component and additional mixed effects or covariate random
functions describing the mixing weights to classify potential pairs into likely
and unlikely transmission pairs. We demonstrate that although sources cannot be
identified at the individual level with certainty, even with the additional
data on time elapsed, inferences into the population-level sources of
transmission are possible, and more accurate than using only phylogenetic data
without time since infection estimates. We apply the approach to estimate
age-specific sources of HIV infection in Amsterdam MSM transmission networks
between 2010-2021. This study demonstrates that infection time estimates
provide informative data to characterize transmission sources, and shows how
phylogenetic source attribution can then be done with multi-dimensional mixture
models.


  The Cichorium genus offers a unique opportunity to study the sporophytic self
incompatibility (SSI) system, being composed of species characterized by highly
efficient SI (C. intybus) and complete self compatibility (C. endivia). The
chicory genome was used to map 7 previously identified SSI locus-associated
markers. The region containing the S locus was restricted to an 4 M bp window
on chromosome 5. Among the genes predicted in this region, MDIS1 INTERACTING
RECEPTOR LIKE KINASE 2 (MIK2) was promising as a candidate for SSI. Its
ortholog in Arabidopsis is involved in pollen stigma recognition reactions, and
its protein structure is similar to that of S-receptor kinase (SRK), a key
component of the SSI in the Brassica genus. The sequencing of MIK2 in chicory
and endive accessions revealed two contrasting scenarios. In C. endivia, MIK2
was fully conserved even comparing different botanical varieties (smooth and
curly). In C. intybus, 387 SNPs and 3 INDELs were identified when comparing
accessions of different biotypes from the same botanical variety (radicchio).
The SNP distribution throughout the gene was uneven, with hypervariable domains
preferentially localized in the LRR-rich extracellular region, putatively
identified as the receptor domain. The gene was hypothesized to be under
positive selection, as the nonsynonymous mutations were more than double the
synonymous ones (dN / dS = 2.17). An analogous situation was observed analyzing
the first 500 bp of the MIK2 promoter: no SNPs were observed among the endive
samples, whereas 44 SNPs and 6 INDELs were detected among the chicory samples.
Further analyses are needed to confirm the role of MIK2 in SSI and to
demonstrate whether the 23 species-specific nonsynonymous SNPs in the CDS
and/or the species-specific 10 bp INDEL found in a CCAAT box region of the
promoter are responsible for the contrasting sexual behaviors of the two
species.


  Across plant communities worldwide, fire regimes reflect a combination of
climatic factors and plant characteristics. To shed new light on the complex
relationships between plant characteristics and fire regimes, we developed a
new conceptual, mechanistic model that includes plant competition, stochastic
fires, and fire-vegetation feedback. Considering a single standing plant
functional type, we observed that highly flammable and slowly colonizing plants
can persist only when they have a strong fire response, while fast colonizing
and less flammable plants can display a larger range of fire responses. At the
community level, the fire response of the strongest competitor determines the
existence of alternative ecological states, i.e. different plant communities,
under the same environmental conditions. Specifically, when the strongest
competitor had a very strong fire response, such as in Mediterranean forests,
only one ecological state could be achieved. Conversely, when the strongest
competitor was poorly fire-adapted, alternative ecological states emerged, for
example between tropical humid savannas and forests, or between different types
of boreal forests. These findings underline the importance of including the
plant fire response when modeling fire ecosystems, e.g. to predict the
vegetation response to invasive species or to climate change.


  Biochemical reaction networks are an amalgamation of reactions where each
reaction represents the interaction of different species. Generally, these
networks exhibit a multi-scale behavior caused by the high variability in
reaction rates and abundances of species. The so-called jump-diffusion
approximation is a valuable tool in the modeling of such systems. The
approximation is constructed by partitioning the reaction network into a fast
and slow subgroup of fast and slow reactions, respectively. This enables the
modeling of the dynamics using a Langevin equation for the fast group, while a
Markov jump process model is kept for the dynamics of the slow group. Most
often biochemical processes are poorly characterized in terms of parameters and
population states. As a result of this, methods for estimating hidden
quantities are of significant interest. In this paper, we develop a tractable
Bayesian inference algorithm based on Markov chain Monte Carlo. The presented
blocked Gibbs particle smoothing algorithm utilizes a sequential Monte Carlo
method to estimate the latent states and performs distinct Gibbs steps for the
parameters of a biochemical reaction network, by exploiting a jump-diffusion
approximation model. The presented blocked Gibbs sampler is based on the two
distinct steps of state inference and parameter inference. We estimate states
via a continuous-time forward-filtering backward-smoothing procedure in the
state inference step. By utilizing bootstrap particle filtering within a
backward-smoothing procedure, we sample a smoothing trajectory. For estimating
the hidden parameters, we utilize a separate Markov chain Monte Carlo sampler
within the Gibbs sampler that uses the path-wise continuous-time representation
of the reaction counters. Finally, the algorithm is numerically evaluated for a
partially observed multi-scale birth-death process example.


  Missing values are a notable challenge when analysing mass spectrometry-based
proteomics data. While the field is still actively debating on the best
practices, the challenge increased with the emergence of mass
spectrometry-based single-cell proteomics and the dramatic increase in missing
values. A popular approach to deal with missing values is to perform
imputation. Imputation has several drawbacks for which alternatives exist, but
currently imputation is still a practical solution widely adopted in
single-cell proteomics data analysis. This perspective discusses the advantages
and drawbacks of imputation. We also highlight 5 main challenges linked to
missing value management in single-cell proteomics. Future developments should
aim to solve these challenges, whether it is through imputation or data
modelling. The perspective concludes with recommendations for reporting missing
values, for reporting methods that deal with missing values and for proper
encoding of missing values.


  The basic reproduction number, $R_0$, is a well-known quantifier of epidemic
spread. However, a class of existing methods for estimating this quantity from
epidemic incidence data can lead to an over-estimation of this quantity. In
particular, when fitting deterministic models to estimate the rate of spread,
we do not account for the stochastic nature of epidemics and that, given the
same system, some outbreaks may lead to epidemics and some may not. Typically,
an observed epidemic that we wish to control is a major outbreak. This amounts
to implicit selection for major outbreaks which leads to the over-estimation
problem. We show that by conditioning a `deterministic' model on major
outbreaks, we can more reliably estimate the basic reproduction number from an
observed epidemic trajectory.


  Cultures around the world show varying levels of conservatism. While
maintaining traditional ideas prevents wrong ones from being embraced, it also
slows or prevents adaptation to new times. Without exploration there can be no
improvement, but often this effort is wasted as it fails to produce better
results, making it better to exploit the best known option. This tension is
known as the exploration/exploitation issue, and it occurs at the individual
and group levels, whenever decisions are made. As such, it is has been
investigated across many disciplines. In this work, we investigate the balance
between exploration and exploitation in changing environments by thinking of
exploration as mutation in a trait space with a varying fitness function.
Specifically, we study how exploration rates evolves by applying adaptive
dynamics to the replicator-mutator equation, under two types of fitness
functions. For the first, payoffs are accrued from playing a two-player,
two-action symmetric game, we consider representatives of all games in this
class and find exploration rates often evolve downwards, but can also undergo
neutral selection as well. Second, we study time dependent fitness with a
function having a single oscillating peak. By increasing the period, we see a
jump in the optimal exploration rate, which then decreases towards zero. These
results establish several possible evolutionary scenarios for exploration
rates, providing insight into many applications, including why we can see such
diversity in rates of cultural change.


  We study the problem of learning conditional distributions of the form $p(G |
\hat G)$, where $G$ and $\hat G$ are two 3D graphs, using continuous
normalizing flows. We derive a semi-equivariance condition on the flow which
ensures that conditional invariance to rigid motions holds. We demonstrate the
effectiveness of the technique in the molecular setting of receptor-aware
ligand generation.


  Integrating whole-slide images (WSIs) and bulk transcriptomics for predicting
patient survival can improve our understanding of patient prognosis. However,
this multimodal task is particularly challenging due to the different nature of
these data: WSIs represent a very high-dimensional spatial description of a
tumor, while bulk transcriptomics represent a global description of gene
expression levels within that tumor. In this context, our work aims to address
two key challenges: (1) how can we tokenize transcriptomics in a semantically
meaningful and interpretable way?, and (2) how can we capture dense multimodal
interactions between these two modalities? Specifically, we propose to learn
biological pathway tokens from transcriptomics that can encode specific
cellular functions. Together with histology patch tokens that encode the
different morphological patterns in the WSI, we argue that they form
appropriate reasoning units for downstream interpretability analyses. We
propose fusing both modalities using a memory-efficient multimodal Transformer
that can model interactions between pathway and histology patch tokens. Our
proposed model, SURVPATH, achieves state-of-the-art performance when evaluated
against both unimodal and multimodal baselines on five datasets from The Cancer
Genome Atlas. Our interpretability framework identifies key multimodal
prognostic factors, and, as such, can provide valuable insights into the
interaction between genotype and phenotype, enabling a deeper understanding of
the underlying biological mechanisms at play. We make our code public at:
https://github.com/ajv012/SurvPath.


  By formulating data samples' formation as a Markov denoising process,
diffusion models achieve state-of-the-art performances in a collection of
tasks. Recently, many variants of diffusion models have been proposed to enable
controlled sample generation. Most of these existing methods either formulate
the controlling information as an input (i.e.,: conditional representation) for
the noise approximator, or introduce a pre-trained classifier in the test-phase
to guide the Langevin dynamic towards the conditional goal. However, the former
line of methods only work when the controlling information can be formulated as
conditional representations, while the latter requires the pre-trained guidance
classifier to be differentiable. In this paper, we propose a novel framework
named RGDM (Reward-Guided Diffusion Model) that guides the training-phase of
diffusion models via reinforcement learning (RL). The proposed training
framework bridges the objective of weighted log-likelihood and maximum entropy
RL, which enables calculating policy gradients via samples from a pay-off
distribution proportional to exponential scaled rewards, rather than from
policies themselves. Such a framework alleviates the high gradient variances
and enables diffusion models to explore for highly rewarded samples in the
reverse process. Experiments on 3D shape and molecule generation tasks show
significant improvements over existing conditional diffusion models.


  The translation of messenger RNA transcripts to proteins is commonly modeled
as a one-dimensional totally asymmetric exclusion process with extended
particles. Here we focus on the effects of premature termination of translation
through the irreversible detachment of ribosomes. We consider a model where the
detachment is induced by the unsuccessful attempt to move to an occupied site.
The model is exactly solvable in a simplified geometry consisting of the
translation initiation region followed by a single slow site representing a
translation bottleneck. In agreement with recent experimental and computational
studies we find a non-monotonic dependence of the ribosome current on the
initiation rate, but only if the leading particle in a colliding pair detaches.
Simulations show that the effect persists for larger lattices and extended
bottlenecks. In the homogeneous system the ribosome density decays
asymptotically as the inverse square root of the distance to the initiation
site.


  Summary: Shu is a visualization tool that integrates diverse data types into
a metabolic map, with a focus on supporting multiple conditions and visualizing
distributions. The goal is to provide a unified platform for handling the
growing volume of multi-omics data, leveraging the metabolic maps developed by
the metabolic modeling community. Additionally, shu offers a streamlined python
API, based on the Grammar of Graphics, for easy integration with data
pipelines.
  Availability and implementation: Freely available at
https://github.com/biosustain/shu under MIT/Apache 2.0 license. Binaries are
available in the release page of the repository and the web app is deployed at
https://biosustain.github.io/shu.


  Motivation: Microsatellite instability (MSI) is a cancer biomarker associated
with cancer prognosis and chemotherapy sensitivity. Since the discovery of MSI,
polymerase chain reaction (PCR)-based testing has been considered the gold
standard for MSI detection. However, with the decrease in sequencing costs,
software that calculates MSI based on next-generation sequencing (NGS) data has
been widely applied. Results: In this study, we evaluated the performance of
the MSIsensor detection software, focusing on the limitations of the chi-square
test algorithm in determining microsatellite stability under high-depth
sequencing data. We demonstrated that the chi-square test algorithm is
insufficient for accurately as-sessing microsatellite stability in this
context. Furthermore, we explored the application of downsampling techniques to
enhance the accuracy of MSIsensor detection. Our findings provide insight into
the limita-tions of current methods and offer potential improvements for more
reliable MSI detection based on NGS data.


  RNA folding prediction remains challenging, but can be also studied using a
topological mathematical approach. In the present paper, the mathematical
method to compute the topological classification of RNA structures and based on
matrix field theory is shortly reviewed, as well as a computational software,
McGenus, used for topological and folding predictions. Additionally, two types
of analysis are performed: the prediction results from McGenus are compared
with topological information extracted from experimentally-determined RNA
structures, and the topology of RNA structures is investigated for biological
significance, in both evolutionary and functional terms. Lastly, we advocate
for more research efforts to be performed at intersection of
physics-mathematics and biology, and in particular about the possible
contributions that topology can provide to the study of RNA folding and
structure.


  The circadian clock can adapt itself to external cues, but the molecular
mechanisms and regulatory networks governing circadian oscillations' transient
adjustments are still largely unknown. Here we consider the specific case of
circadian oscillations transiently responding to a temperature change. Using a
framework motivated by Floquet theory, we model the mRNA expression level of
the fat body from Drosophila melanogaster following a step change from 25C to
18C. Using the method we infer the adaptation rates of individual genes as they
adapt to the new temperature. To deal with heteroskedastic noise and outliers
present in the expression data we employ quantile regression and wild bootstrap
for significance testing. Model selection with finite-sample corrected Akaike
Information Criterion (AICc) is performed additionally for robust inference. We
identify several genes with fast transition rates as potential sources of
temperature-mediated responses in the circadian system of fruit flies, and the
constructed network suggests that the proteasome may play important roles in
governing these responses.


  In the analysis of complex traits, genetic effects are frequently modelled as
either fixed or random effects. Such assumptions serve as a foundation of
defining heritability and relatedness using genome-wide single nucleotide
polymorphism (SNP) markers. In this work, I propose a genealogical framework
connecting the two assumptions conditional on the ancestral recombination graph
(ARG). It turns out that the reference time point in which the probability is
defined determines whether the effect of a variant to behave as either a fixed
or random effect. This lays a connection between the PC regression and linear
mixed model (LMM) for genetic association study. The framework induces a
genetic relatedness matrix (GRM) in which the elements are a function of the
time to the recent common ancestor. Subsequently, a novel trait variance
decomposition respect to the regions sharing a common genealogy is followed.
The variance decomposition then provides a natural means to define heritability
and a potential method for gene mapping.


  In this paper, Wolbachia infection in a bisexual and diploid population with
a fixed cytoplasmic incompatibility rate $w$ and maternal transmission rate $d$
is studied as an evolution algebra. As the cytoplasmic incompatibility (CI) of
the population causes deaths in the offspring, the evolution algebra of this
model is not baric, and is a dibaric algebra if and only if the cytoplasmic
incompatibility rate $w$ is 1 and $d=1$. The idempotent elements are given in
terms of $d$ and $w$. Moreover, this algebra has no absolute nilpotent elements
when CI expression $w \neq 1$.


  Life is on the razor's edge as resulting from competitive birth and death
random forces. We illustrate this aphorism in the context of three Markov chain
population models where systematic random immigration events promoting growth
are simultaneously balanced with random emigration ones provoking thinning. The
origin of mass removals are either determined by external demands or by aging
leading to different conditions of stability.


  The Lotka-Euler equation is a mathematical expression used to study
population dynamics and growth, particularly in the context of demography and
ecology. The growth rate $\lambda$ is the speed at which an individual produce
their offspring. It is essentially a birth process, and here it is shown that
by reversing the process to a death process, in which individuals die at a rate
$\lambda^{-1}$, the derivation of the Lotka-Euler equation becomes more
intuitive and direct, both in discrete and continuous time.


  Dedicated accelerator hardware has become essential for processing AI-based
workloads, leading to the rise of novel accelerator architectures. Furthermore,
fundamental differences in memory architecture and parallelism have made these
accelerators targets for scientific computing.
  The sequence alignment problem is fundamental in bioinformatics; we have
implemented the $X$-Drop algorithm, a heuristic method for pairwise alignment
that reduces search space, on the Graphcore Intelligence Processor Unit (IPU)
accelerator. The $X$-Drop algorithm has an irregular computational pattern,
which makes it difficult to accelerate due to load balancing.
  Here, we introduce a graph-based partitioning and queue-based batch system to
improve load balancing. Our implementation achieves $10\times$ speedup over a
state-of-the-art GPU implementation and up to $4.65\times$ compared to CPU. In
addition, we introduce a memory-restricted $X$-Drop algorithm that reduces
memory footprint by $55\times$ and efficiently uses the IPU's limited
low-latency SRAM. This optimization further improves the strong scaling
performance by $3.6\times$.


  The issue of reproducibility of computational models and the related FAIR
principles (findable, accessible, interoperable, and reusable) are examined in
a specific test case. I analyze a computational model of the segment polarity
network in Drosophila embryos published in 2000. Despite the high number of
citations to this publication, 23 years later the model is barely accessible,
and consequently not interoperable. Following the text of the original
publication allowed successfully encoding the model for the open source
software COPASI. Subsequently saving the model in the SBML format allowed it to
be reused in other open source software packages. Submission of this SBML
encoding of the model to the BioModels database enables its findability and
accessibility. This demonstrates how the FAIR principles can be successfully
enabled by using open source software, widely adopted standards, and public
repositories, facilitating reproducibility and reuse of computational cell
biology models that will outlive the specific software used.


  Transcription commonly occurs in bursts, with alternating productive (ON) and
quiescent (OFF) periods, governing mRNA production rates. Yet, how
transcription is regulated through bursting dynamics remains unresolved. In
this study, we conduct real-time measurements of endogenous transcriptional
bursting with single-mRNA sensitivity. Leveraging the diverse transcriptional
activities in early fly embryos, we uncover stringent relationships between
bursting parameters. Specifically, we find that the durations of ON and OFF
periods are linked. Regardless of the developmental stage or body-axis
position, gene activity levels predict the average ON and OFF periods of
individual alleles. Lowly transcribing alleles predominantly modulate OFF
durations (burst frequency), while highly transcribing alleles primarily tune
ON durations (burst size). Importantly, these relationships persist even under
perturbation of cis-regulatory elements or trans-factors. This suggests a novel
mechanistic constraint governing bursting dynamics rather than a modular
control of distinct parameters by distinct regulatory processes. Our study
provides a foundation for future investigations into the molecular mechanisms
underpinning spatiotemporal transcriptional control.


  Every interaction of a living organism with its environment involves the
placement of a bet. Armed with partial knowledge about a stochastic world, the
organism must decide its next step or near-term strategy, an act that
implicitly or explicitly involves the assumption of a model of the world.
Better information about environmental statistics can improve the bet quality,
but in practice resources for information gathering are always limited. We
argue that theories of optimal inference dictate that ``complex'' models are
harder to infer with bounded information and lead to larger prediction errors.
Thus, we propose a principle of ``playing it safe'' where, given finite
information gathering capacity, biological systems should be biased towards
simpler models of the world, and thereby to less risky betting strategies. In
the framework of Bayesian inference, we show that there is an optimally safe
adaptation strategy determined by the Bayesian prior. We then demonstrate that,
in the context of stochastic phenotypic switching by bacteria, implementation
of our principle of ``playing it safe'' increases fitness (population growth
rate) of the bacterial collective. We suggest that the principle applies
broadly to problems of adaptation, learning and evolution, and illuminates the
types of environments in which organisms are able to thrive.


  When a biological system robustly corrects component-level errors, the direct
pressure on component performance declines. Components may become less
reliable, maintain more genetic variability, or drift neutrally in design,
creating the basis for new forms of organismal complexity. This article links
the protection-decay dynamic to other aspects of robust and complex systems.
Examples include the hourglass pattern of biological development and Doyle's
hourglass architecture for robustly complex systems in engineering. The deeply
and densely connected wiring architecture in biology's cellular controls and in
machine learning's computational neural networks provide another link. By
unifying these seemingly different aspects into a unified framework, we gain a
new perspective on robust and complex systems.


  The polyvagal theory (PT), offered by Porges (2021), proposes that the
autonomic nervous system (ANS) was repurposed in mammals, via a second vagal
nerve, to suppress defensive strategies and support the expression of
sociality. Three critical assumptions of this theory are that (1) the
transition of the ANS was associated with the evolution of social mammals from
asocial reptiles; (2) the transition enabled mammals, unlike their reptilian
ancestors, to derive a biological benefit from social interactions; and (3) the
transition forces a less parsimonious explanation (convergence) for the
evolution of social behavior in birds and mammals, since birds evolved from a
reptilian lineage. Two recently published reviews, however, provided compelling
evidence that the social asocial dichotomy is overly simplistic, neglects the
diversity of vertebrate social systems, impedes our understanding of the
evolution of social behavior, and perpetuates the erroneous belief that one
group, non-avian reptiles, is incapable of complex social behavior. In the
worst case, if PT depends upon a transition from asocial reptiles to social
mammals, then the ability of PT to explain the evolution of the mammalian ANS
is highly questionable. A great number of social behaviors occur in both
reptiles and mammals. In the best case, PT has misused the terms social and
asocial. Even here, however, the theory would still need to identify a
particular suite of behaviors found in mammals and not reptiles that could be
associated with, or explain, the transition of the ANS, and then replace the
asocial and social labels with more specific descriptors.


  Historically, musings about the structure of ecological communities has
revolved around the structure of pairwise interactions, competition, predation,
mutualism, etc. . . Recently a growing literature acknowledges that the
baseline assumption that the pair of species is not necessarily the
metaphorical molecule of community ecology, and that certain structures
containing three or more species may not be usefully divisible into pairwise
components. Two examples are intransitive competition (species A dominates
species B dominates species C dominates species A), and nonlinear higher-order
effects. While these two processes have been discussed extensively, the
explicit analysis of how the two of them behave when simultaneously part of the
same dynamic system has not yet appeared in the literature. A concrete
situation exists on coffee farms in Puerto Rico in which three ant species, at
least on some farms, form an intransitive competitive triplet, and that triplet
is strongly influenced, nonlinearly, by a fly parasitoid that modifies the
competitive ability of one of the species in the triplet. Using this
arrangement as a template we explore the dynamical consequences with a simple
ODE model. Results are complicated and include include alternative periodic and
chaotic attractors. The qualitative structures of those complications, however,
may be retrieved easily from a reflection on the basic natural history of the
system.


  Tuna and billfish are charismatic pelagic fishes attracting considerable
scientific attention due to their ecophysiological and socioeconomic
importance. However, the knowledge of their basin-wide spawning and larval
habitats, especially in a warming ocean, is limited. We use the largest
available dataset on tuna and billfish larvae in the Pacific Ocean to build a
geostatistical species-distribution model with high explanatory power. The
results reveal the spatial distribution of tuna and billfish larvae through all
seasons across the Pacific. The model also identifies the optimal temperature
ranges for nine major species and assesses the potential impact of ocean
warming on larval distributions. We additionally present evidence that
environmental variables, such as pH, phosphate concentration, and sea-surface
height, exert secondary effects on larval distributions that warrant further
investigation. Our findings make a quantum leap in understanding the
ecophysiology of tuna and billfish, providing valuable information for future
conservation efforts.


  From the smallest biological systems to the largest cosmological structures,
spatial domains undergo expansion and contraction. Within these growing
domains, diffusive transport is a common phenomenon. Mathematical models have
been widely employed to investigate diffusive processes on growing domains.
However, a standard assumption is that the domain growth is spatially uniform.
There are many relevant examples where this is not the case, such as the
colonisation of growing gut tissue by neural crest cells. As such, it is not
straightforward to disentangle the individual roles of heterogeneous growth and
diffusive transport. Here we present exact solutions to models of diffusive
transport on domains undergoing spatially non-uniform growth. The exact
solutions are obtained via a combination of transformation, convolution and
superposition techniques. We verify the accuracy of these solutions via
comparison with simulations of a corresponding lattice-based random walk. We
explore various domain growth functions, including linear growth, exponential
growth and contraction, and oscillatory growth. Provided the domain size
remains positive, we find that the derived solutions are valid. The exact
solutions reveal the relationship between model parameters, such as the
diffusivity and the type and rate of domain growth, and key statistics, such as
the survival and splitting probabilities.