  Increasing attention has been drawn to the misuse of statistical methods over
recent years, with particular concern about the prevalence of practices such as
poor experimental design, cherry-picking and inadequate reporting. These
failures are largely unintentional and no more common in ecology than in other
scientific disciplines, with many of them easily remedied given the right
guidance.
  Originating from a discussion at the 2020 International Statistical Ecology
Conference, we show how ecologists can build their research following four
guiding principles for impactful statistical research practices: 1. Define a
focused research question, then plan sampling and analysis to answer it; 2.
Develop a model that accounts for the distribution and dependence of your data;
3. Emphasise effect sizes to replace statistical significance with ecological
relevance; 4. Report your methods and findings in sufficient detail so that
your research is valid and reproducible.
  Listed in approximate order of importance, these principles provide a
framework for experimental design and reporting that guards against unsound
practices. Starting with a well-defined research question allows researchers to
create an efficient study to answer it, and guards against poor research
practices that lead to false positives and poor replicability. Correct and
appropriate statistical models give sound conclusions, good reporting practices
and a focus on ecological relevance make results impactful and replicable.
  Illustrated with an example from a recent study into the impact of
disturbance on upland swamps, this paper explains the rationale for the
selection and use of effective statistical practices and provides practical
guidance for ecologists seeking to improve their use of statistical methods.


  Random Boolean networks have been used widely to explore aspects of gene
regulatory networks. A modified form of the model through which to
systematically explore the effects of increasing the number of gene states has
previously been introduced. In this paper, these discrete dynamical networks
are coevolved within coupled, rugged fitness landscapes to explore their
behaviour. Results suggest the general properties of the Boolean model remain
with higher valued logic regardless of the update scheme or fitness sampling
method. Introducing topological asymmetry in the coevolving networks is seen to
alter behaviour.


  In this paper, we discuss a hierarchical model, based on a Ricker competition
model. The species considered are competing for resources and may be subject to
an Allee effect due to mate limitation, anti-predator vigilance or aggression,
cooperative predation or resource defense, or social thermoregulation. The
species may be classified into a more dominant species and less dominant or
``wimpy" species or just as a predator and prey. The model under consideration
also has components taking into account immigration in both species and a more
structured Holling type II functional response. Local and global stability
analyses are discussed and simulations are provided. We also consider
demographic stochasticity on the species due to environmental fluctuations in
the form of Wiener processes and we show that there are conditions under which
a global solution exists, stationary distributions exists, and strong
persistence in mean of the species is possible. We also use Wasserstein
distance to show empirically that stochasticity can acts as bifurcation
parameter.


  Dengue and Zika incidence data and the latest research have raised questions
about how dengue vaccine strategies might be impacted by the emergence of Zika
virus. Existing antibodies to one virus might temporarily protect or promote
infection by the other through antibody-dependent enhancement (ADE). With this
condition, understanding the dynamics of propagation of these two viruses is of
great importance when implementing vaccines. In this work, we analyze the
effect of vaccination against one strain, in a two-strain model that accounts
for cross-immunity and ADE. Using basic and invasion reproductive numbers, we
examined the dynamics of the model and provide conditions to ensure the
stability of the disease-free equilibrium. We provide conditions on
cross-immunity, ADE and vaccination rate under which the vaccination could
ensure the global stability of the disease-free equilibrium. The results
indicate scenarios in which vaccination against one strain may improve or
worsen the control of the other, as well as contribute to the eradication or
persistence of one or both viruses in the population.


  Molecular representations are of fundamental importance for the modeling and
analysis of molecular systems. Representation models and in general approaches
based on topological data analysis (TDA) have demonstrated great success in
various steps of drug design and materials discovery. Here we develop a
mathematically rigorous computational framework for molecular representation
based on the persistent Dirac operator. The properties of the spectrum of the
discrete weighted and unweighted Dirac matrices are systemically discussed and
used to demonstrate the geometric and topological properties of both
non-homology and homology eigenvectors of real molecular structures. This
allows us to asses the influence of weighting schemes on the information
encoded in the Dirac eigenspectrum. A series of physical persistent attributes,
which characterize the spectrum of the Dirac matrices across a filtration, are
proposed and used as efficient molecular fingerprints. Finally, our persistent
Dirac-based model is used for clustering molecular configurations from nine
types of organic-inorganic halide perovskites. We found that our model can
cluster the structures very well, demonstrating the representation and
featurization power of the current approach.


  The human neurotensin receptor 1 (NTSR1) is a G protein-coupled receptor. The
receptor is activated by a small peptide ligand neurotensin. NTSR1 can be
expressed in HEK cells by stable transfection. Previously we used the
fluorescent protein markers mRuby3 or mNeonGreen fused to NTSR1 for EMCCD-based
structured illumination microscopy (SIM) in living HEK cells. Ligand binding
induced conformational changes in NTSR1 which triggered the intracellular
signaling processes. Recent single-molecule studies revealed a dynamic
monomer/dimer equilibrium of this receptor in artificial lipid bilayers. Here
we report on the oligomerization state of human NTSR1 from living cells by
trapping them into lipid nanodiscs. Briefly, SMALPs (styrene-maleic acid
copolymer lipid nanoparticles) were produced directly from the plasma membranes
of living HEK293T FlpIn cells. SMALPs with a diameter of 15 nm were soluble and
stable. NTSR1 in SMALPs were analyzed by single-molecule intensity measurements
one membrane patch at a time using a custom-built confocal anti-Brownian
electrokinetic trap (ABEL trap) microscope. We found oligomerization changes
before and after stimulation of the receptor with its ligand neurotensin.


  Epidemic propagation on networks represents an important departure from
traditional massaction models. However, the high-dimensionality of the exact
models poses a challenge to both mathematical analysis and parameter inference.
By using mean-field models, such as the pairwise model (PWM), the complexity
becomes tractable. While such models have been used extensively for model
analysis, there is limited work in the context of statistical inference. In
this paper, we explore the extent to which the PWM with the
susceptible-infected-recovered (SIR) epidemic can be used to infer disease- and
network-related parameters. The widely-used MLE approach exhibits several
issues pertaining to parameter unidentifiability and a lack of robustness to
exact knowledge about key quantities such as population size and/or proportion
of under reporting. As an alternative, we considered the recently developed
dynamical survival analysis (DSA). For scenarios in which there is no model
mismatch, such as when data are generated via simulations, both methods perform
well despite strong dependence between parameters. However, for real-world
data, such as foot-and-mouth, H1N1 and COVID19, the DSA method appears more
robust to potential model mismatch and the parameter estimates appear more
epidemiologically plausible. Taken together, however, our findings suggest that
network-based mean-field models can be used to formulate approximate
likelihoods which, coupled with an efficient inference scheme, make it possible
to not only learn about the parameters of the disease dynamics but also that of
the underlying network.


  Recurrent COVID-19 outbreaks have placed immense strain on the hospital
system in Quebec. We develop a Bayesian three-state coupled Markov switching
model to analyze COVID-19 outbreaks across Quebec based on admissions in the 30
largest hospitals. Within each catchment area, we assume the existence of three
states for the disease: absence, a new state meant to account for many zeroes
in some of the smaller areas, endemic and outbreak. Then we assume the disease
switches between the three states in each area through a series of coupled
nonhomogeneous hidden Markov chains. Unlike previous approaches, the transition
probabilities may depend on covariates and the occurrence of outbreaks in
neighboring areas, to account for geographical outbreak spread. Additionally,
to prevent rapid switching between endemic and outbreak periods we introduce
clone states into the model which enforce minimum endemic and outbreak
durations. We make some interesting findings, such as that mobility in retail
and recreation venues had a positive association with the development and
persistence of new COVID-19 outbreaks in Quebec. Based on model comparison our
contributions show promise in improving state estimation retrospectively and in
real-time, especially when there are smaller areas and highly spatially
synchronized outbreaks. Furthermore, our approach offers new and interesting
epidemiological interpretations, such as being able to estimate the effect of
covariates on disease extinction.


  The advances in variational inference are providing promising paths in
Bayesian estimation problems. These advances make variational phylogenetic
inference an alternative approach to Markov Chain Monte Carlo methods for
approximating the phylogenetic posterior. However, one of the main drawbacks of
such approaches is modelling the prior through fixed distributions, which could
bias the posterior approximation if they are distant from the current data
distribution. In this paper, we propose an approach and an implementation
framework to relax the rigidity of the prior densities by learning their
parameters using a gradient-based method and a neural network-based
parameterization. We applied this approach for branch lengths and evolutionary
parameters estimation under several Markov chain substitution models. The
results of performed simulations show that the approach is powerful in
estimating branch lengths and evolutionary model parameters. They also show
that a flexible prior model could provide better results than a predefined
prior model. Finally, the results highlight that using neural networks improves
the initialization of the optimization of the prior density parameters.


  Elucidating exocytosis processes provide insights into cellular
neurotransmission mechanisms, and may have potential in neurodegenerative
diseases research. Amperometry is an established electrochemical method for the
detection of neurotransmitters released from and stored inside cells. An
important aspect of the amperometry method is the sub-millisecond temporal
resolution of the current recordings which leads to several hundreds of
gigabytes of high-quality data. In this study, we present a universal method
for the classification with respect to diverse amperometric datasets using
data-driven approaches in computational science. We demonstrate a very high
prediction accuracy (greater than or equal to 95%). This includes an end-to-end
systematic machine learning workflow for amperometric time series datasets
consisting of pre-processing; feature extraction; model identification;
training and testing; followed by feature importance evaluation - all
implemented. We tested the method on heterogeneous amperometric time series
datasets generated using different experimental approaches, chemical
stimulations, electrode types, and varying recording times. We identified a
certain overarching set of common features across these datasets which enables
accurate predictions. Further, we showed that information relevant for the
classification of amperometric traces are neither in the spiky segments alone,
nor can it be retrieved from just the temporal structure of spikes. In fact,
the transients between spikes and the trace baselines carry essential
information for a successful classification, thereby strongly demonstrating
that an effective feature representation of amperometric time series requires
the full time series. To our knowledge, this is one of the first studies that
propose a scheme for machine learning, and in particular, supervised learning
on full amperometry time series data.


  Amperometry is a commonly used electrochemical method for studying the
process of exocytosis in real-time. Given the high precision of recording that
amperometry procedures offer, the volume of data generated can span over
several hundreds of megabytes to a few gigabytes and therefore necessitates
systematic and reproducible methods for analysis. Though the spike
characteristics of amperometry traces in the time domain hold information about
the dynamics of exocytosis, these biochemical signals are, more often than not,
characterized by time-varying signal properties. Such signals with time-variant
properties may occur at different frequencies and therefore analyzing them in
the frequency domain may provide statistical validation for observations
already established in the time domain. This necessitates the use of
time-variant, frequency-selective signal processing methods as well, which can
adeptly quantify the dominant or mean frequencies in the signal. The Fast
Fourier Transform (FFT) is a well-established computational tool that is
commonly used to find the frequency components of a signal buried in noise. In
this work, we outline a method for spike-based frequency analysis of
amperometry traces using FFT that also provides statistical validation of
observations on spike characteristics in the time domain. We demonstrate the
method by utilizing simulated signals and by subsequently testing it on diverse
amperometry datasets generated from different experiments with various chemical
stimulations. To our knowledge, this is the first fully automated open-source
tool available dedicated to the analysis of spikes extracted from amperometry
signals in the frequency domain.


  We present a differentiable formulation of abstract chemical reaction
networks (CRNs) that can be trained to solve a variety of computational tasks.
Chemical reaction networks are one of the most fundamental computational
substrates used by nature. We study well-mixed single-chamber systems, as well
as systems with multiple chambers separated by membranes, under mass-action
kinetics. We demonstrate that differentiable optimisation, combined with proper
regularisation, can discover non-trivial sparse reaction networks that can
implement various sorts of oscillators and other chemical computing devices.


  Researchers have employed stochastic simulations to determine the validity of
their theoretical findings and to study analytically intractable spreading
dynamics. In both cases, the correctness and efficiency of the simulation
algorithm are of paramount importance. We prove in this article that the Next
Reaction Method and the non-Markovian Gillespie algorithm, two algorithms for
simulating non-Markovian epidemics, are statistically equivalent. We also study
the performance and applicability under various circumstances through
complexity analyses and numerical experiments. In our numerical simulations, we
apply the Next Reaction Method and the Gillespie algorithm to epidemic
simulations on time-varying networks and epidemic simulations with cooperative
infections. Both tasks have only been done using the Gillespie algorithm, while
we show that the Next Reaction Method is a good alternative. We believe this
article may also serve as a guide for choosing simulation algorithms that are
both correct and efficient for researchers from epidemiology and beyond.


  Populations are often subject to catastrophes that cause mass removal of
individuals. Many stochastic growth models have been considered to explain such
dynamics. Among the results reported, it has been considered whether dispersion
strategies, at times of catastrophes, increase the survival probability of the
population. In this paper, we contrast dispersion strategies comparing mean
extinction times of the population when extinction occurs almost surely. In
particular, we consider populations subject to binomial catastrophes, that is,
the population size is reduced according to a binomial law when a catastrophe
occurs. Our findings delineate the optimal strategy (dispersion or
non-dispersion) based on variations in model parameter values.


  In recent years, the sequencing, assembling and annotation of prokaryotic
genomes has become increasingly easy and cheap. Thus it becomes increasingly
feasible and interesting to perform comparative genomics analyses of new
genomes to those of related organisms. Thereby related organisms can be defined
by different criteria, such as taxonomy or phenotype.
  Expectations regarding the contents of genomes are often expressed in
scientific articles describing group of organisms. Evaluating such
expectations, when a new genome becomes available, requires analysing the text
snippets which express such expectations, extracting the logical elements of
the text and enabling a formal expression, more suitable for further automated
analyses.
  Hereby we present a theoretical framework, alongside practical consideration
for expressing expectations about the content of genomes, with the purpose of
enabling such comparative genomics analyses. The components of the framework
include a system for the definition of groups of organisms, supported by a
Prokaryotic Group Types Ontology, a system for the definition of genomic
contents, supported by a Prokaryotic Genomic Contents Definition Ontology.
Finally we discuss how the combination of these two systems may enable an
unambiguous definition of absolute and relative genome content expectation
rules.


  Spatially resolved transcriptomics brings exciting breakthroughs to
single-cell analysis by providing physical locations along with gene
expression. However, as a cost of the extremely high spatial resolution, the
cellular level spatial transcriptomic data suffer significantly from missing
values. While a standard solution is to perform imputation on the missing
values, most existing methods either overlook spatial information or only
incorporate localized spatial context without the ability to capture long-range
spatial information. Using multi-head self-attention mechanisms and positional
encoding, transformer models can readily grasp the relationship between tokens
and encode location information. In this paper, by treating single cells as
spatial tokens, we study how to leverage transformers to facilitate spatial
tanscriptomics imputation. In particular, investigate the following two key
questions: (1) $\textit{how to encode spatial information of cells in
transformers}$, and (2) $\textit{ how to train a transformer for transcriptomic
imputation}$. By answering these two questions, we present a transformer-based
imputation framework, SpaFormer, for cellular-level spatial transcriptomic
data. Extensive experiments demonstrate that SpaFormer outperforms existing
state-of-the-art imputation algorithms on three large-scale datasets.


  Combining a spatiotemporal, multi-agent based model of a foraging ecosystem
with linear, genetically programmed rules for the agents' behaviors results in
implicit, endogenous, objective functions and selection algorithms based on
"natural selection". Use of this implicit optimization of genetic programs for
study of biological systems is tested by application to an artificial foraging
ecosystem, and compared with established biological, ecological, and stochastic
gene diffusion models. Limited program memory and execution time constraints
emulate real-time and concurrent properties of physical and biological systems,
and stress test the optimization algorithms. Relative fitness of the agents'
programs and efficiency of the resultant populations as functions of these
constraints gauge optimization effectiveness and efficiency. Novel solutions
confirm the creativity of the optimization process and provide an unique
opportunity to experimentally test the neutral code bloating hypotheses. Use of
this implicit, endogenous, evolutionary optimization of spatially interacting,
genetically programmed agents is thus shown to be novel, consistent with
biological systems, and effective and efficient in discovering fit and novel
solutions.


  Forecasting disease spread is a critical tool to help public health officials
design and plan public health interventions.However, the expected future state
of an epidemic is not necessarily well defined as disease spread is inherently
stochastic, contact patterns within a population are heterogeneous, and
behaviors change. In this work, we use time-dependent probability generating
functions (PGFs) to capture these characteristics by modeling a stochastic
branching process of the spread of a disease over a network of contacts in
which public health interventions are introduced over time. To achieve this, we
define a general transmissibility equation to account for varying transmission
rates (e.g. masking), recovery rates (e.g. treatment), contact patterns (e.g.
social distancing) and percentage of the population immunized (e.g.
vaccination). The resulting framework allows for a temporal and probabilistic
analysis of an intervention's impact on disease spread, which match
continuous-time stochastic simulations that are much more computationally
expensive.To aid policy making, we then define several metrics over which
temporal and probabilistic intervention forecasts can be compared: Looking at
the expected number of cases and the worst-case scenario over time, as well as
the probability of reaching a critical level of cases and of not seeing any
improvement following an intervention.Given that epidemics do not always follow
their average expected trajectories and that the underlying dynamics can change
over time, our work paves the way for more detailed short-term forecasts of
disease spread and more informed comparison of intervention strategies.


  Autism spectrum disorder (ASD) is a heterogeneous neurodevelopmental disorder
(NDD) that is caused by genetic, epigenetic, and environmental factors. Recent
advances in genomic analysis have uncovered numerous candidate genes with
common and/or rare mutations that increase susceptibility to ASD. In addition,
there is increasing evidence that copy number variations (CNVs), single
nucleotide polymorphisms (SNPs), and unusual de novo variants negatively affect
neurodevelopment pathways in various ways. The overall rate of copy number
variants found in patients with autism is 10%-20%, of which 3%-7% can be
detected cytogenetically. Although the role of submicroscopic CNVs in ASD has
been studied recently, their association with genomic loci and genes has not
been properly studied. In this review, we focus on 47 ASD-associated CNV
regions and their related genes. Here, we identify 1,632 protein-coding genes
and long non-coding RNAs (lncRNAs) within these regions. Among them, 552 are
significantly expressed in the brain. Using a list of ASD-associated genes from
SFARI, we detect 17 regions containing at least one known ASD-associated
protein-coding genes. Of the remaining 30 regions, we identify 24 regions
containing at least one protein-coding genes with brain-enriched expression and
nervous system phenotype in mouse mutant and one lncRNAs with both
brain-enriched expression and upregulation in iPSC to neuron differentiation.
Our analyses highlight the diversity of genetic lesions of CNV regions that
contribute to ASD and provide new genetic evidence that lncRNA genes may
contribute to etiology of ASD. In addition, the discovered CNVs will be a
valuable resource for diagnostic facilities, therapeutic strategies, and
research in terms of variation priority.


  Black-box optimization methods play an important role in many fields of
computational simulation. In particular, such methods are often used in the
design and modelling of biological systems, including proteins and their
complexes with various ligands. This work is mainly focused on the
protein-protein docking that plays a key role in modern drug-design workflows.
We develop a black-box approach for such docking problems using a novel
technique based on the tensor-train decomposition of high-dimensional
interaction functions. Our method shows an advantage in terms of the discovered
global minima and has a high potential for further implementation on a wide
range of devices, including graphical processing units and quantum processing
units.


  A reaction-diffusion Leslie-Gower predator-prey model, incorporating the fear
effect and prey refuge, with Beddington-DeAngelis functional response, is
introduced. A qualitative analysis of the solutions of the model and the
stability analysis of the coexistence equilibrium, are performed. Sufficient
conditions guaranteeing the occurrence of Turing instability have been
determined either in the case of self-diffusion or in the case of
cross-diffusion. Different types of Turing patterns, representing a spatial
redistribution of population in the environment, emerge for different values of
the model parameters.


  Atrazine as a weed control tactic profoundly changed world agriculture.
Long-term use revealed resistant biotypes, R, with a single base pair mutation
of the chloroplast psbA gene. The R phenotype emerged from a sequential cascade
of pleiotropic effects from the plastid to the whole plant. This reorganization
of the R biotype revealed photosynthetic regulation at different levels of
plant organization. The environment affected R plant productivity differently
than in the susceptible, S, biotype. A consistent, differential, pattern of
photosynthesis was observed between R and S over the diurnal light period.
Photosynthetic superiority of a biotype was a function of the time of day,
plant age temperature. Under highly favorable environmental conditions S often
had the advantage over R. Under less favorable, stressful, conditions R can be
at an advantage over S. Pleiotropic reorganization revealed a sun-air-leaf
Shannon communication system, providing insights into the complex interaction
of chloroplast components in photosynthetic regulation. Altered plastid
thylakoid and stomatal function regulate how the R leaf utilizes the sun-air
environment. Movement of sun-air messages demonstrated how sunlight and air are
modified to a usable message for carbon fixation. These insights showed how
agriculture changed weed populations, and how resistant weed populations
changed agriculture. These changes changed herbicide resistance: introduction
of herbicide resistant crops, HRC. The development of HRCs extended the
evolutionary reach of R weeds. R weed biotypes were naturally selected in these
introduced HRCs: an evolutionary spiral of human technology extended the
phenotypic reach of R biotypes.


  Biological systems can be studied at multiple levels of information,
including gene, protein, RNA and different interaction networks levels. The
goal of this work is to explore how the fusion of systems' level information
with temporal dynamics of gene expression can be used in combination with
non-linear approximation power of deep neural networks to predict novel gene
functions in a non-model organism potato \emph{Solanum tuberosum}. We propose
DDeMON (Dynamic Deep learning from temporal Multiplex Ontology-annotated
Networks), an approach for scalable, systems-level inference of function
annotation using time-dependent multiscale biological information. The proposed
method, which is capable of considering billions of potential links between the
genes of interest, was applied on experimental gene expression data and the
background knowledge network to reliably classify genes with unknown function
into five different functional ontology categories, linked to the experimental
data set. Predicted novel functions of genes were validated using extensive
protein domain search approach.


  As the availability of omics data has increased in the last few years, more
multi-omics data have been generated, that is, high-dimensional molecular data
consisting of several types such as genomic, transcriptomic, or proteomic data,
all obtained from the same patients. Such data lend themselves to being used as
covariates in automatic outcome prediction because each omics type may
contribute unique information, possibly improving predictions compared to using
only one omics data type. Frequently, however, in the training data and the
data to which automatic prediction rules should be applied, the test data, the
different omics data types are not available for all patients. We refer to this
type of data as block-wise missing multi-omics data. First, we provide a
literature review on existing prediction methods applicable to such data.
Subsequently, using a collection of 13 publicly available multi-omics data
sets, we compare the predictive performances of several of these approaches for
different block-wise missingness patterns. Finally, we discuss the results of
this empirical comparison study and draw some tentative conclusions.


  Optimal transport (OT) theory focuses, among all maps
$T:\mathbb{R}^d\rightarrow \mathbb{R}^d$ that can morph a probability measure
onto another, on those that are the ``thriftiest'', i.e. such that the averaged
cost $c(x, T(x))$ between $x$ and its image $T(x)$ be as small as possible.
Many computational approaches have been proposed to estimate such Monge maps
when $c$ is the $\ell_2^2$ distance, e.g., using entropic maps [Pooladian'22],
or neural networks [Makkuva'20, Korotin'20]. We propose a new model for
transport maps, built on a family of translation invariant costs $c(x,
y):=h(x-y)$, where $h:=\tfrac{1}{2}\|\cdot\|_2^2+\tau$ and $\tau$ is a
regularizer. We propose a generalization of the entropic map suitable for $h$,
and highlight a surprising link tying it with the Bregman centroids of the
divergence $D_h$ generated by $h$, and the proximal operator of $\tau$. We show
that choosing a sparsity-inducing norm for $\tau$ results in maps that apply
Occam's razor to transport, in the sense that the displacement vectors
$\Delta(x):= T(x)-x$ they induce are sparse, with a sparsity pattern that
varies depending on $x$. We showcase the ability of our method to estimate
meaningful OT maps for high-dimensional single-cell transcription data, in the
$34000$-$d$ space of gene counts for cells, without using dimensionality
reduction, thus retaining the ability to interpret all displacements at the
gene level.


  Leukemia, the cancer of blood cells, originates in the blood-forming cells of
the bone marrow. In Chronic Myeloid Leukemia (CML) conditions, the cells
partially become mature that look like normal white blood cells but do not
resist infection effectively. Early detection of CML is important for effective
treatment, but there is a lack of routine screening tests. Regular check-ups
and monitoring of symptoms are the best way to detect CML in the early stages.
In the study, we developed a multi-layer-perception-based meta-ensemble system
using protein amino acid sequences for early risk prediction of CML. The
deleterious mutation analysis of protein sequences provides 7discriminant
information in amino acid sequences causing CML. The protein sequences are
expressed into molecular descriptors using the values of hydrophobicity and
hydrophilicity of the amino acids. 9 These descriptors are transformed in
various statistical and correlation-based feature spaces. These 10 features
information is given to several diverse types of base learners. The preliminary
predictions of 11 base-learners are employed to develop Multi-Layered
Perceptron (MLP) based meta-ensemble. The 12 proposed learning approach
effectively utilizes the discriminant information to classify CML/non- 13 CML
protein sequences. The proposed prediction system has given improved results
and it can be 14 employed as a potential biomarker for early diagnosis of CML.


  The dynamics of fractional population sizes y_i=Y_i/N in homogeneous
compartment models with time dependent total population N is analyzed. Assuming
constant per capita birth and death rates the vector field Y_i'=V_i(Y)
naturally projects to a vector field F_i(Y) tangent to the leaves of constant
population N. A universal formula for the projected field F_i is given. In this
way, in many SIR-type models with standard incidence all demographic parameters
become redundant for the dynamical system y_i'=F_i(y). They may be put to zero
by shifting remaining parameters appropriately. Normalizing eight examples from
the literature this way, they unexpectedly become isomorphic for corresponding
parameter ranges. Thus, some recently published results turn out to be already
covered by papers 20 years ago.


  Denoising diffusion probabilistic models (DDPMs) have recently taken the
field of generative modeling by storm, pioneering new state-of-the-art results
in disciplines such as computer vision and computational biology for diverse
tasks ranging from text-guided image generation to structure-guided protein
design. Along this latter line of research, methods have recently been proposed
for generating 3D molecules using equivariant graph neural networks (GNNs)
within a DDPM framework. However, such methods are unable to learn important
geometric and physical properties of 3D molecules during molecular graph
generation, as they adopt molecule-agnostic and non-geometric GNNs as their 3D
graph denoising networks, which negatively impacts their ability to effectively
scale to datasets of large 3D molecules. In this work, we address these gaps by
introducing the Geometry-Complete Diffusion Model (GCDM) for 3D molecule
generation, which outperforms existing 3D molecular diffusion models by
significant margins across conditional and unconditional settings for the QM9
dataset as well as for the larger GEOM-Drugs dataset. Importantly, we
demonstrate that the geometry-complete denoising process GCDM learns for 3D
molecule generation allows the model to generate realistic and stable large
molecules at the scale of GEOM-Drugs, whereas previous methods fail to do so
with the features they learn. Additionally, we show that GCDM's geometric
features can effectively be repurposed to directly optimize the geometry and
chemical composition of existing 3D molecules for specific molecular
properties, demonstrating new, real-world versatility of molecular diffusion
models. Our source code, data, and reproducibility instructions are freely
available at https://github.com/BioinfoMachineLearning/Bio-Diffusion.


  Immunotherapies have been proven to have significant therapeutic efficacy in
the treatment of cancer. The last decade has seen adoptive cell therapies, such
as chimeric antigen receptor T-cell (CART-cell) therapy, gain FDA approval
against specific cancers. Additionally, there are numerous clinical trials
ongoing investigating additional designs and targets. Nevertheless, despite the
excitement and promising potential of CART-cell therapy, response rates to
therapy vary greatly between studies, patients, and cancers. There remains an
unmet need to develop computational frameworks that more accurately predict
CART-cell function and clinical efficacy. Here we present a coarse-grained
model simulated with logical rules that demonstrates the evolution of signaling
signatures following the inter-action between CART-cells and tumor cells and
allows for in silico based prediction of CART-cell functionality prior to
experimentation.


  The endoplasmic reticulum (ER) is a dynamic network of interconnected sheets
and tubules that orchestrates the distribution of lipids, ions, and proteins
throughout the cell. The impact of its complex, dynamic morphology on its
function as an intracellular transport hub remains poorly understood. To
elucidate the functional consequences of ER network structure and dynamics, we
quantify how the heterogeneity of the peripheral ER in COS7 cells affects
diffusive protein transport. In vivo imaging of photoactivated ER membrane
proteins demonstrates their non-uniform spreading to adjacent regions, in a
manner consistent with simulations of diffusing particles on extracted network
structures. Using a minimal network model to represent tubule rearrangements,
we demonstrate that ER network dynamics are sufficiently slow to have little
effect on diffusive protein transport. Furthermore, stochastic simulations
reveal a novel consequence of ER network heterogeneity: the existence of 'hot
spots' where sparse diffusive reactants are more likely to find one another.
Intriguingly, ER exit sites are disproportionately found in these highly
accessible regions. Combining in vivo experiments with analytic calculations,
quantitative image analysis, and computational modeling, we demonstrate how
structure guides diffusive protein transport and reactions in the ER.


  In this study, we investigate the dynamics of a spatial and non spatial
prey-predator interaction model that includes the following: (i) fear effect
incorporated in prey birth rate; (ii) group defence of prey against predators;
and (iii) prey refuge. We provide comprehensive mathematical analysis of
extinction and persistence scenarios for both prey and predator species. To
better explore the dynamics of the system, a thorough investigation of
bifurcation analysis has been performed using fear level, prey birth rate, and
prey death rate caused by intra-prey competition as bifurcation parameter. All
potential occurrences of bi-stability dynamics have also been investigated for
some relevant sets of parametric values. Our numerical evaluations show that
high levels of fear can stabilize the prey-predator system by ruling out the
possibility of periodic solutions. Also, our model Hopf bifurcation is
subcritical in contrast to traditional prey-predator models, which ignore the
cost of fear and have supercritical Hopf bifurcations in general. In contrast
to the general trend, predator species go extinct at higher values of prey
birth rates. We have also found that, contrary to the typical tendency for prey
species to go extinct, both prey and predator populations may coexist in the
system as intra-prey competition level grows noticeably. The stability and
Turing instability of associated spatial model have also been investigated
analytically. We also perform the numerical simulation to observe the effect of
different parameters on the density distribution of species. Different types of
spatiotemporal patterns like spot, mixture of spots and stripes have been
observed via variation of time evolution, diffusion coefficient of predator
population, level of fear factor and prey refuge. The fear level parameter (k)
has a great impact on the spatial dynamics of model system.


  Delays and stochasticity have both served as crucially valuable ingredients
in mathematical descriptions of control, physical, and biological systems. In
this work, we investigate how explicitly dynamical stochasticity in delays
modulates the effect of delayed feedback. To do so, we consider a hybrid model
where stochastic delays evolve by a continuous-time Markov chain, and between
switching events, the system of interest evolves via a deterministic delay
equation. Our main contribution is the calculation of an effective delay
equation in the fast switching limit. This effective equation maintains the
influence of all subsystem delays and cannot be replaced with a single
effective delay. To illustrate the relevance of this calculation, we
investigate a simple model of stochastically switching delayed feedback
motivated by gene regulation. We show that sufficiently fast switching between
two oscillatory subsystems can yield stable dynamics.


  Previous studies have inferred robust stability of reaction networks by
utilizing linear programs or iterative algorithms. Such algorithms become
tedious or computationally infeasible for large networks. In addition, they
operate like black-boxes without offering intuition for the structures that are
necessary to maintain stability. In this work, we provide several graphical
criteria for constructing robust stability certificates, checking robust
non-degeneracy, {verifying persistence}, and establishing global stability. By
characterizing a set of stability-preserving graph modifications that includes
the enzymatic modification motif, we show that the stability of arbitrarily
large nonlinear networks can be examined by simple visual inspection. We show
applications of this technique to ubiquitous motifs in systems biology such as
Post-Translational Modification (PTM) cycles, the Ribosome Flow Model (RFM),
T-cell kinetic proofreading and others.


  A very high rate for cyclic electron flow (CEF) around PSI (~180 s-1 or 210
s-1 in minimum medium or in the presence of a carbon source respectively) is
measured in the presence of methyl viologen (MV) in intact cells of
Chlamydomonas reinhardtii under anaerobic conditions. The observation of an
efficient CEF in the presence of methyl viologen is in agreement with the
previous results reports of Asada et al in broken chloroplasts (Plant Cell
Physiol. 31(4) (1990) 557-564). From the analysis of the P700 and PC absorbance
changes, we propose that a confinement between 2 PC molecules, 1 PSI and 1
cytb6f corresponding to a functional supercomplex is responsible for these high
rates of CEF. Supercomplex formation is also observed in the absence of methyl
viologen, but with lower maximal CEF rate (about 100 s-1) suggesting that this
compound facilitates the mediation of electron transfer from PSI acceptors to
the stromal side of cytb6f. Further analysis of CEF in mutants of Chlamydomonas
defective in state transitions shows the requirement of a kinase-driven
transition to state 2 to establish this functional supercomplex configuration.
However, a movement of the LHCII antennae is not involved in this process. We
discuss the possible involvement of auxiliary proteins, among which is a small
cytb6f-associated polypeptide, the PETO protein, which is one of the targets of
the STT7 kinase.


  The critical Kauffman model with connectivity one is the simplest class of
critical Boolean networks. Nevertheless, it exhibits intricate behavior at the
boundary of order and chaos. We introduce a formalism for expressing the
dynamics of multiple loops as a product of the dynamics of individual loops.
Using it, we prove that the number of attractors scales as $2^m$, where $m$ is
the number of nodes in loops - as fast as possible, and much faster than
previously believed.


  Chemical reactions are usually studied under the assumption that both
substrates and catalysts are well mixed (WM) throughout the system. Although
this is often applicable to test-tube experimental conditions, it is not
realistic in cellular environments, where biomolecules can undergo
liquid-liquid phase separation (LLPS) and form condensates, leading to
important functional outcomes, including the modulation of catalytic action.
Similar processes may also play a role in protocellular systems, like primitive
coacervates, or in membrane-assisted prebiotic pathways. Here we explore
whether the de-mixing of catalysts could lead to the formation of
micro-environments that influence the kinetics of a linear (multi-step)
reaction pathway, as compared to a WM system. We implemented a general lattice
model to simulate LLPS of an ensemble of different catalysts and extended it to
include diffusion and a sequence of reactions of small substrates. We carried
out a quantitative analysis of how the phase separation of the catalysts
affects reaction times depending on the affinity between substrates and
catalysts, the length of the reaction pathway, the system size, and the degree
of homogeneity of the condensate. A key aspect underlying the differences
reported between the two scenarios is that the scale invariance observed in the
WM system is broken by condensation processes. The main theoretical
implications of our results for mean-field chemistry are drawn, extending the
mass action kinetics scheme to include substrate initial hitting times to reach
the catalysts condensate. We finally test this approach by considering open
non-linear conditions, where we successfully predict, through microscopic
simulations, that phase separation inhibits chemical oscillatory behaviour,
providing a possible explanation for the marginal role that this complex
dynamic behaviour plays in real metabolisms.


  In biology, stochastic branching processes with a two-stage, hierarchical
structure arise in the study of population dynamics, gene expression, and
phylogenetic inference. These models have been commonly analyzed using
generating functions, the method of characteristics and various perturbative
approximations. Here we describe a general method for analyzing hierarchic
first-order reaction networks using Lie theory. Crucially, we identify the fact
that the Lie group associated to hierarchic reaction networks decomposes as a
wreath product of the groups associated to the subnetworks of the independent
and dependent types. After explaining the general method, we illustrate it on a
model of population dynamics and the so-called two-state or telegraph model of
single-gene transcription. Solutions to such processes provide essential input
to downstream methods designed to attempt to infer parameters of these and
related models.


  Many biological systems approach physical limits to their performance,
motivating the idea that their behavior and underlying mechanisms could be
determined by such optimality. Nevertheless, optimization as a predictive
principle has only been applied in very simplified setups. Here, in contrast,
we explore a mechanistically-detailed class of models for the gap gene network
of the Drosophila embryo, and determine its 50+ parameters by optimizing the
information that gene expression levels convey about nuclear positions, subject
to physical constraints on the number of available molecules. Optimal networks
recapitulate the architecture and spatial gene expression profiles of the real
organism. Our framework makes precise the many tradeoffs involved in maximizing
functional performance, and allows us to explore alternative networks to
address the questions of necessity vs contingency. Multiple solutions to the
optimization problem may be realized in closely related organisms.


  Structure-based drug design, i.e., finding molecules with high affinities to
the target protein pocket, is one of the most critical tasks in drug discovery.
Traditional solutions, like virtual screening, require exhaustively searching
on a large molecular database, which are inefficient and cannot return novel
molecules beyond the database. The pocket-based 3D molecular generation model,
i.e., directly generating a molecule with a 3D structure and binding position
in the pocket, is a new promising way to address this issue. Herein, we propose
VD-Gen, a novel pocket-based 3D molecular generation pipeline. VD-Gen consists
of several carefully designed stages to generate fine-grained 3D molecules with
binding positions in the pocket cavity end-to-end. Rather than directly
generating or sampling atoms with 3D positions in the pocket like in early
attempts, in VD-Gen, we first randomly initialize many virtual particles in the
pocket; then iteratively move these virtual particles, making the distribution
of virtual particles approximate the distribution of molecular atoms. After
virtual particles are stabilized in 3D space, we extract a 3D molecule from
them. Finally, we further refine atoms in the extracted molecule by iterative
movement again, to get a high-quality 3D molecule, and predict a confidence
score for it. Extensive experiment results on pocket-based molecular generation
demonstrate that VD-Gen can generate novel 3D molecules to fill the target
pocket cavity with high binding affinities, significantly outperforming
previous baselines.


  The central problem in cryo-electron microscopy (cryo-EM) is to recover the
3D structure from noisy 2D projection images which requires estimating the
missing projection angles (poses). Recent methods attempted to solve the 3D
reconstruction problem with the autoencoder architecture, which suffers from
the latent vector space sampling problem and frequently produces suboptimal
pose inferences and inferior 3D reconstructions. Here we present an improved
autoencoder architecture called ACE (Asymmetric Complementary autoEncoder),
based on which we designed the ACE-EM method for cryo-EM 3D reconstructions.
Compared to previous methods, ACE-EM reached higher pose space coverage within
the same training time and boosted the reconstruction performance regardless of
the choice of decoders. With this method, the Nyquist resolution (highest
possible resolution) was reached for 3D reconstructions of both simulated and
experimental cryo-EM datasets. Furthermore, ACE-EM is the only amortized
inference method that reached the Nyquist resolution.


  This is the first paper in the sequel studying the Wolbachia-infection in
bisexual populations. This paper considers the behavior of the population as a
discrete dynamical system. The recurrence relation is obtained as a function of
the initial infected male/female frequencies and the cytoplasmic
incompatibility of the population. The experimental data from
Wolbachia-infected terrestrial isopod populations and the model proposed in
Wolbachia-infected mosquitoes from literature is compared with the discrete
dynamical system achieved.


  Studies on somatic mutations in cancer cells DNA, their roles in tumour
growth and progression between successive stages are of high importance for
improving understanding of cancer evolution. Important insights into scenarios
of cancer growth, roles of somatic mutations and types/strengths of
evolutionary forces they introduce are gained by using mathematical and
computer modelling. Previous studies developed mathematical models of cancer
cell evolution with driver and passenger somatic mutations. Driver mutations
were assumed to have a strong advantageous effect on cancer cell population
growth, while passenger mutations were considered as fully neutral or mildly
deleterious.
  In this paper, following several experimental results, we develop models of
cancer evolution with somatic mutations introducing weakly advantageous force
to the evolution of cancer cells. Our models belong to two classes,
deterministic and stochastic. Deterministic models are systems of differential
master equations of balances of average numbers of cells and mutations in
evolution. We, additionally introduce a modification in equations for balances
aimed at accounting for effects on the finite size of the cancer cell
population. The novelty of our deterministic modelling is deriving parameters
of travelling waves of advantageous mutations and quantification of the pattern
of the population size growth. A stochastic model based on the Gillespie
algorithm is used to verify the results of deterministic modelling. We confront
predictions of modelling with some observational data on cancer evolution. We
propose the scenario of evolution driven by a large number of weakly
advantageous mutations in transition from the early stage to the invasive form
of breast cancer. We also fit our models to observed patterns of rates of
growth of tumours.


  Molecular conformation generation (MCG) is a fundamental and important
problem in drug discovery. Many traditional methods have been developed to
solve the MCG problem, such as systematic searching, model-building, random
searching, distance geometry, molecular dynamics, Monte Carlo methods, etc.
However, they have some limitations depending on the molecular structures.
Recently, there are plenty of deep learning based MCG methods, which claim they
largely outperform the traditional methods. However, to our surprise, we design
a simple and cheap algorithm (parameter-free) based on the traditional methods
and find it is comparable to or even outperforms deep learning based MCG
methods in the widely used GEOM-QM9 and GEOM-Drugs benchmarks. In particular,
our design algorithm is simply the clustering of the RDKIT-generated
conformations. We hope our findings can help the community to revise the deep
learning methods for MCG. The code of the proposed algorithm could be found at
https://gist.github.com/ZhouGengmo/5b565f51adafcd911c0bc115b2ef027c.


  Is there a unified model for generating molecules considering different
conditions, such as binding pockets and chemical properties? Although
target-aware generative models have made significant advances in drug design,
they do not consider chemistry conditions and cannot guarantee the desired
chemical properties. Unfortunately, merging the target-aware and chemical-aware
models into a unified model to meet customized requirements may lead to the
problem of negative transfer. Inspired by the success of multi-task learning in
the NLP area, we use prefix embeddings to provide a novel generative model that
considers both the targeted pocket's circumstances and a variety of chemical
properties. All conditional information is represented as learnable features,
which the generative model subsequently employs as a contextual prompt.
Experiments show that our model exhibits good controllability in both single
and multi-conditional molecular generation. The controllability enables us to
outperform previous structure-based drug design methods. More interestingly, we
open up the attention mechanism and reveal coupling relationships between
conditions, providing guidance for multi-conditional molecule generation.


  Molecular docking, given a ligand molecule and a ligand binding site (called
``pocket'') on a protein, predicting the binding mode of the protein-ligand
complex, is a widely used technique in drug design. Many deep learning models
have been developed for molecular docking, while most existing deep learning
models perform docking on the whole protein, rather than on a given pocket as
the traditional molecular docking approaches, which does not match common
needs. What's more, they claim to perform better than traditional molecular
docking, but the approach of comparison is not fair, since traditional methods
are not designed for docking on the whole protein without a given pocket. In
this paper, we design a series of experiments to examine the actual performance
of these deep learning models and traditional methods. For a fair comparison,
we decompose the docking on the whole protein into two steps, pocket searching
and docking on a given pocket, and build pipelines to evaluate traditional
methods and deep learning methods respectively. We find that deep learning
models are actually good at pocket searching, but traditional methods are
better than deep learning models at docking on given pockets. Overall, our work
explicitly reveals some potential problems in current deep learning models for
molecular docking and provides several suggestions for future works.


  The last decade has witnessed a surge of theoretical and computational models
to describe the dynamics of complex gene regulatory networks, and how these
interactions can give rise to multistable and heterogeneous cell populations.
As the use of theoretical modeling to describe genetic and biochemical circuits
becomes more widespread, theoreticians with mathematical and physical
backgrounds routinely apply concepts from statistical physics, non-linear
dynamics, and network theory to biological systems. This review aims at
providing a clear overview of the most important methodologies applied in the
field while highlighting current and future challenges. It also includes
hands-on tutorials to solve and simulate some of the archetypical biological
system models used in the field. Furthermore, we provide concrete examples from
the existing literature for theoreticians that wish to explore this
fast-developing field. Whenever possible, we highlight the similarities and
differences between biochemical and regulatory networks and 'classical' systems
typically studied in non-equilibrium statistical and quantum mechanics.


  Activity cliffs (ACs), which are generally defined as pairs of structurally
similar molecules that are active against the same bio-target but significantly
different in the binding potency, are of great importance to drug discovery. Up
to date, the AC prediction problem, i.e., to predict whether a pair of
molecules exhibit the AC relationship, has not yet been fully explored. In this
paper, we first introduce ACNet, a large-scale dataset for AC prediction. ACNet
curates over 400K Matched Molecular Pairs (MMPs) against 190 targets, including
over 20K MMP-cliffs and 380K non-AC MMPs, and provides five subsets for model
development and evaluation. Then, we propose a baseline framework to benchmark
the predictive performance of molecular representations encoded by deep neural
networks for AC prediction, and 16 models are evaluated in experiments. Our
experimental results show that deep learning models can achieve good
performance when the models are trained on tasks with adequate amount of data,
while the imbalanced, low-data and out-of-distribution features of the ACNet
dataset still make it challenging for deep neural networks to cope with. In
addition, the traditional ECFP method shows a natural advantage on MMP-cliff
prediction, and outperforms other deep learning models on most of the data
subsets. To the best of our knowledge, our work constructs the first
large-scale dataset for AC prediction, which may stimulate the study of AC
prediction models and prompt further breakthroughs in AI-aided drug discovery.
The codes and dataset can be accessed by https://drugai.github.io/ACNet/.


  We introduce a multi-allele Wright-Fisher model with non-recurrent,
reversible mutation and directional selection. In this setting, the allele
frequencies at a single locus track the path of a hybrid jump-diffusion process
with state space given by the vertex and edge set of a graph. Vertices
represent monomorphic population states and edge-positions mark the biallelic
proportions of ancestral and derived alleles during polymorphic segments. We
derive the stationary distribution in mutation-selection-drift equilibrium and
obtain the expected allele frequency spectrum under large population size
scaling. For the extended model with multiple independent loci we derive
rigorous upper bounds for a wide class of associated measures of genetic
variation. Within this framework we present mathematically precise arguments to
conclude that the presence of directional selection reduces the magnitude of
genetic variation, as constrained by the bounds for neutral evolution.


  There are many problems in biochemistry that are difficult to study
experimentally. Simulation methods are appealing due to direct availability of
atomic coordinates as a function of time. However, direct molecular simulations
are challenged by the size of systems and the time scales needed to describe
relevant motions. In theory, enhanced sampling algorithms can help to overcome
some of the limitations of molecular simulations. Here, we discuss a problem in
biochemistry that offers a significant challenge for enhanced sampling methods
and that could, therefore, serve as a benchmark for comparing approaches that
use machine learning to find suitable collective variables. More in particular,
we study the transitions LacI undergoes upon moving between being
non-specifically and specifically bound to DNA. It is found that many degrees
of freedom change during this transition and that the transition does not occur
reversibly in simulations if only a subset of these degrees of freedom are
biased. We also explain why this problem is so important to biologists and the
transformative impact that a simulation of it would have on the understanding
of DNA regulation.


  Ecosystems with a large number of species are often modelled as
Lotka-Volterra dynamical systems built around a large random interaction
matrix. Under some known conditions, a global equilibrium exists and is unique.
In this article, we rigorously study its statistical properties in the large
dimensional regime. Such an equilibrium vector is known to be the solution of a
so-called Linear Complementarity Problem (LCP). We describe its statistical
properties by designing an Approximate Message Passing (AMP) algorithm, a
technique that has recently aroused an intense research effort in the fields of
statistical physics, Machine Learning, or communication theory. Interaction
matrices taken from the Gaussian Orthogonal Ensemble, or following a Wishart
distribution are considered. Beyond these models, the AMP approach developed in
this article has the potential to describe the statistical properties of
equilibria associated to more involved interaction matrix models.


  Discovering novel drug candidate molecules is one of the most fundamental and
critical steps in drug development. Generative deep learning models, which
create synthetic data given a probability distribution, have been developed
with the purpose of picking completely new samples from a partially known
space. Generative models offer high potential for designing de novo molecules;
however, in order for them to be useful in real-life drug development
pipelines, these models should be able to design target-specific molecules,
which is the next step in this field. In this study, we propose DrugGEN, for
the de novo design of drug candidate molecules that interact with selected
target proteins. The proposed system represents compounds and protein
structures as graphs and processes them via serially connected two generative
adversarial networks comprising graph transformers. DrugGEN is trained using a
large dataset of compounds from ChEMBL and target-specific bioactive molecules,
to design effective and specific inhibitory molecules against the AKT1 protein,
which has critical importance for developing treatments against various types
of cancer. On fundamental benchmarks, DrugGEN models have either competitive or
better performance against other methods. To assess the target-specific
generation performance, we conducted further in silico analysis with molecular
docking and deep learning-based bioactivity prediction. Results indicate that
de novo molecules have high potential for interacting with the AKT1 protein
structure in the level of its native ligand. DrugGEN can be used to design
completely novel and effective target-specific drug candidate molecules for any
druggable protein, given target features and a dataset of experimental
bioactivities. Code base, datasets, results and trained models of DrugGEN are
available at https://github.com/HUBioDataLab/DrugGEN


  In some patients of myeloproliferative neoplasm, two genetic mutations can be
found: JAK2 V617F and TET2. When one mutation is present or not, the other
mutation has different effects on regulating gene expressions. Besides, when
both mutations are present, the order of occurrence might make a difference. In
this paper, we build nonlinear ordinary differential equation models and Markov
chain models to explain such phenomena.


  We devise a method for predicting certain receptor-ligand binding behaviors,
based on stochastic dynamical modelling. We consider the dynamics of a receptor
binding to a ligand on the cell membrane, where the receptor and ligand perform
different motions and are thus modeled by stochastic differential equations
with Gaussian noise or non-Gaussian noise. We use neural networks based on
Onsager-Machlup function to compute the probability $P_1$ of the unbounded
receptor diffusing to the cell membrane. Meanwhile, we compute the probability
$P_2$ of extracellular ligand arriving at the cell membrane by solving the
associated Fokker-Planck equation. Then, we could predict the most probable
binding probability by combining $P_1$ and $P_2$. In this way, we conclude with
some indication about where the ligand will most probably encounter the
receptor, contributing to better understanding of cell's response to external
stimuli and communication with other cells.


  Identification of fossil species is crucial to evolutionary studies. Recent
advances from deep learning have shown promising prospects in fossil image
identification. However, the quantity and quality of labeled fossil images are
often limited due to fossil preservation, conditioned sampling, and expensive
and inconsistent label annotation by domain experts, which pose great
challenges to training deep learning based image classification models. To
address these challenges, we follow the idea of the wisdom of crowds and
propose a multiview ensemble framework, which collects Original (O), Gray (G),
and Skeleton (S) views of each fossil image reflecting its different
characteristics to train multiple base models, and then makes the final
decision via soft voting. Experiments on the largest fusulinid dataset with
2400 images show that the proposed OGS consistently outperforms baselines
(using a single model for each view), and obtains superior or comparable
performance compared to OOO (using three base models for three the same
Original views). While considering the identification consistency estimation
with respect to human experts, OGS receives the highest agreement with the
original labels of dataset and with the re-identifications of two human
experts. We conclude that the proposed framework can present state-of-the-art
performance in the fusulinid fossil identification case study. This framework
is designed for general fossil identification and it is expected to see
applications to other fossil datasets in future work. Notably, the result,
which shows more performance gains as train set size decreases or over a
smaller imbalance fossil dataset, suggests the potential application to
identify rare fossil images. The proposed framework also demonstrates its
potential for assessing and resolving inconsistencies in fossil identification.


  Asymmetric damage segregation (ADS) is ubiquitous among unicellular
organisms: After a mother cell divides, its two daughter cells receive
sometimes slightly, sometimes strongly different fractions of damaged proteins
accumulated in the mother cell. Previous studies demonstrated that ADS provides
a selective advantage over symmetrically dividing cells by rejuvenating and
perpetuating the population as a whole. In this work we focus on the
statistical properties of damage in individual lineages and the overall damage
distributions in growing populations for a variety of ADS models with different
rules governing damage accumulation, segregation, and the lifetime dependence
on damage. We show that for a large class of deterministic ADS rules the
trajectories of damage along the lineages are chaotic, and the distributions of
damage in cells born at a given time asymptotically becomes fractal. By
exploiting the analogy of linear ADS models with the Iterated Function Systems
known in chaos theory, we derive the Frobenius-Perron equation for the
stationary damage density distribution and analytically compute the damage
distribution moments and fractal dimensions. We also investigate nonlinear and
stochastic variants of ADS models and show the robustness of the salient
features of the damage distributions.


  In a Bienaym\'{e}-Galton-Watson process for which there is a
positiveprobability for individuals of having no offspring, there is a
subtlebalance and dependence between the sterile nodes (the dead nodes or
leaves)and the prolific ones (the productive nodes) both at and up to the
currentgeneration. We explore the many facets of this problem, especially in
thecontext of an exactly solvable linear-fractional branching mechanism at
allgeneration. Eased asymptotic issues are investigated. Relation of
thisspecial branching process to skip-free to the left and simple random
walks'excursions is then investigated. Mutual statistical information on
theirshapes can be learnt from this association.


  Biological entities are inherently dynamic. As such, various ecological
disciplines use mathematical models to describe temporal evolution. Typically,
growth curves are modelled as sigmoids, with the evolution modelled by ordinary
differential equations. Among the various sigmoid models, the logistic and
Gompertz equations are well established and widely used in fitting growth data
in the fields of biology and ecology. This paper suggests a statistical
interpretation of the logistic equation within the general framework. This
interpretation is based on modelling the population environment relationship,
the mathematical theory of which we discuss in detail. By applying this theory,
we obtain stochastic evolutionary equations, for which the logistic equation is
a limiting case. The prospect of modifying logistic population growth is
discussed. We support our statistical interpretation of population growth
dynamics with test numerical simulations. We show that the Gompertz equation
and other related models can be treated in a similar way.


  We consider a time-continuous Markov branching process of proliferating cells
with a countable collection of types. Among-type transitions are inspired by
the Tug-of-War process introduced in McFarland et al. as a mathematical model
for competition of advantageous driver mutations and deleterious passenger
mutations in cancer cells. We introduce a version of the model in which a
driver mutation pushes the type of the cell $L$-units up, while a passenger
mutation pulls it $1$-unit down. The distribution of time to divisions depends
on the type (fitness) of cell, which is an integer. The extinction probability
given any initial cell type is strictly less than $1$, which allows us to
investigate the transition between types (type transition) in an infinitely
long cell lineage of cells. The analysis leads to the result that under driver
dominance, the type transition process escapes to infinity, while under
passenger dominance, it leads to a limit distribution. Implications in cancer
cell dynamics and population genetics are discussed.


  A physics-informed neural network (PINN) embedded with the
susceptible-infected-removed (SIR) model is devised to understand the temporal
evolution dynamics of infectious diseases. Firstly, the effectiveness of this
approach is demonstrated on synthetic data as generated from the numerical
solution of the susceptible-asymptomatic-infected-recovered-dead (SAIRD) model.
Then, the method is applied to COVID-19 data reported for Germany and shows
that it can accurately identify and predict virus spread trends. The results
indicate that an incomplete physics-informed model can approach more
complicated dynamics efficiently. Thus, the present work demonstrates the high
potential of using machine learning methods, e.g., PINNs, to study and predict
epidemic dynamics in combination with compartmental models.


  Dynamic full-field optical coherence tomography (D-FFOCT) has recently
emerged as an invaluable live label-free and non-invasive imaging modality able
to image subcellular biological structures and their metabolic activity within
complex 3D samples. However, D-FFOCT suffers from fringe artefacts when imaging
nearby reflective surfaces and is highly sensitive to vibrations. Here, we
present interface Self-Referenced (iSR) D-FFOCT, an alternative configuration
to D-FFOCT that takes advantage of the presence of the sample coverslip in
between the sample and the objective by using it as a defocused reference arm,
thus avoiding the aforementioned artefacts. We demonstrate the ability of iSR
D-FFOCT to image 2D fibroblast cell cultures, which are among the flattest
mammalian cells.


  Structural information of phylogenetic tree topologies plays an important
role in phylogenetic inference. However, finding appropriate topological
structures for specific phylogenetic inference tasks often requires significant
design effort and domain expertise. In this paper, we propose a novel
structural representation method for phylogenetic inference based on learnable
topological features. By combining the raw node features that minimize the
Dirichlet energy with modern graph representation learning techniques, our
learnable topological features can provide efficient structural information of
phylogenetic trees that automatically adapts to different downstream tasks
without requiring domain expertise. We demonstrate the effectiveness and
efficiency of our method on a simulated data tree probability estimation task
and a benchmark of challenging real data variational Bayesian phylogenetic
inference problems.


  In May 2022, a cluster of mpox cases were detected in the UK that could not
be traced to recent travel history from an endemic region. Over the coming
months, the outbreak grew, with over 3000 total cases reported in the UK, and
similar outbreaks occurring worldwide. These outbreaks appeared linked to
sexual contact networks between gay, bisexual and other men who have sex with
men. Following the COVID-19 pandemic, local health systems were strained, and
therefore effective surveillance for mpox was essential for managing public
health policy. However, the mpox outbreak in the UK was characterised by
substantial delays in the reporting of the symptom onset date and specimen
collection date for confirmed positive cases. These delays led to substantial
backfilling in the epidemic curve, making it challenging to interpret the
epidemic trajectory in real-time. Many nowcasting models exist to tackle this
challenge in epidemiological data, but these lacked sufficient flexibility. We
have developed a novel nowcasting model using generalised additive models to
correct the mpox epidemic curve in England, and provide real-time
characteristics of the state of the epidemic, including the real-time growth
rate. This model benefited from close collaboration with individuals involved
in collecting and processing the data, enabling temporal changes in the
reporting structure to be built into the model, which improved the robustness
of the nowcasts generated.


  DNA origami is a widely used method to construct nanostructures by
self-assembling designed DNA strands. These structures are often used as
"breadboards" for templated assembly of proteins, gold nanoparticles, aptamers,
and other molecules, with applications ranging from therapeutics and
diagnostics to plasmonics and photonics. Imaging these structures using AFM or
TEM is not capable to capture their full conformation ensemble as they only
show their structure flattened on a surface. However, certain conformations of
the nanostructure can position guest molecules into distances unaccounted for
in their intended design, thus leading to spurious interactions between guest
molecules that are designed to be separated. Here, we use molecular dynamics
simulations to capture conformational ensemble of 2D DNA origami tiles and show
that introducing single-stranded overhangs, which are typically used for
functionalization of the origami with guest molecules, induces a curvature of
the tile structure in the bulk. We show that the shape deformation is of
entropic origin, with implications for design of robust DNA origami breadboards
as well as potential approach to modulate structure shape by introducing
overhangs. We then verify experimentally that the overhangs introduce curvature
into the DNA origami tiles.


  SBcoyote is an open-source cross-platform biochemical reaction viewer and
editor released under the liberal MIT license. It is written in Python and uses
wxPython to implement the GUI and the drawing canvas. It supports the
visualization and editing of compartments, species, and reactions. It includes
many options to stylize each of these components. For instance, species can be
in different colors and shapes. Other core features include the ability to
create alias nodes, alignment of groups of nodes, network zooming, as well as
an interactive bird-eye view of the network to allow easy navigation on large
networks. A unique feature of the tool is the extensive Python plugin API,
where third-party developers can include new functionality. To assist
third-party plugin developers, we provide a variety of sample plugins,
including, random network generation, a simple auto layout tool, export to
Antimony, export SBML, import SBML, etc. Of particular interest are the export
and import SBML plugins since these support the SBML level 3 layout and render
standard, which is exchangeable with other software packages. Plugins are
stored in a GitHub repository, and an included plugin manager can retrieve and
install new plugins from the repository on demand. Plugins have version
metadata associated with them to make it install plugin updates. Availability:
https://github.com/sys-bio/SBcoyote.


  The Unfolded Protein Response is the cell mechanism for maintaining the
balance of properly folded proteins in the endoplasmic reticulum , the
specialized cellular compartment. Although it is largely studied from a
biological point of view, much of the literature lacks a quantitative analysis
of such a central signaling pathway. In this work, we aim to fill this gap by
applying structural controllability analysis of complex networks to several
Unfolded Protein Response networks to identify crucial nodes in the signaling
flow. In particular, we first build different network models of the Unfolded
Protein Response mechanism, relying on data contained in various
protein-protein interaction databases. Then, we identify the driver nodes,
essential for overall network control, i.e., the key proteins on which external
stimulation may be optimally delivered to control network behavior. Our
structural controllability analysis results show that the driver nodes commonly
identified across databases match with known endoplasmic reticulum stress
sensors. This potentially confirms that the theoretically identified drivers
correspond to the biological key proteins associated with fundamental cellular
activities and diseases. In conclusion, we prove that structural
controllability is a reliable quantitative tool to investigate biological
signaling pathways, and it can be potentially applied to networks more complex
and less explored than Unfolded Protein Response.


  An important problem in evolutionary genomics is to investigate whether a
certain trait measured on each sample is associated with the sample
phylogenetic tree. The phylogenetic tree represents the shared evolutionary
history of the samples and it is usually estimated from molecular sequence data
at a locus or from other type of genetic data. We propose a model for trait
evolution inspired by the Chinese Restaurant Process that includes a parameter
that controls the degree of preferential attachment, that is, the tendency of
nodes in the tree to subtend from nodes of the same type. This model with no
preferential attachment is equivalent to a structured coalescent model with
simultaneous migration and coalescence events and serves as a null model. We
derive a test for phylogenetic binary trait association with linear
computational complexity and empirically demonstrate that it is more powerful
than some other methods. We apply our test to study the phylogenetic
association of some traits in swordtail fish, breast cancer, yellow fever virus
and influenza A H1N1 virus. R package implementation of our methods is
available at $\texttt{https://github.com/jyzhang27/CRPTree}$.


  A major challenge in synthetic genetic circuit development is the
inter-dependency between heterologous gene expressions by circuits and host's
growth rate. Increasing heterologous gene expression increases burden to the
host, resulting in host growth reduction; which reduces overall heterologous
protein abundance. Hence, it is difficult to design predictable genetic
circuits. Here, we develop two biophysical models; one for promoter, another
for RBS; to correlate heterologous gene expression and growth reduction. We
model cellular resource allocation in E. coli to describe the burden, as growth
reduction, caused by genetic circuits. To facilitate their uses in genetic
circuit design, inputs to the model are common characteristics of biological
parts [e.g. relative promoter strength (RPU) and relative ribosome binding
sites strength (RRU)]. The models suggest that E. coli's growth rate reduces
linearly with increasing RPU / RRU of the genetic circuits; thus, providing 2
handy models taking parts characteristics as input to estimate growth rate
reduction for fine tuning genetic circuit design in silico prior to
construction. Our promoter model correlates well with experiments using various
genetic circuits, both single and double expression cassettes, up to a relative
promoter unit of 3.7 with a 60% growth rate reduction (average R2 ~ 0.9).


  Phylogenetic networks generalize phylogenetic trees in order to model
reticulation events. Although the comparison of phylogenetic trees is well
studied, and there are multiple ways to do it in an efficient way, the
situation is much different for phylogenetic networks.
  Some classes of phylogenetic networks, mainly tree-child networks, are known
to be classified efficiently by their $\mu$-representation, which essentially
counts, for every node, the number of paths to each leaf. In this paper, we
introduce the extended $\mu$-representation of networks, where the number of
paths to reticulations is also taken into account. This modification allows us
to distinguish orchard networks and to define a sound metric on the space of
such networks that can, moreover, be computed efficiently.
  The class of orchard networks, as well as being one of the classes with
biological significance (one such network can be interpreted as a tree with
extra arcs involving coexisting organisms), is one of the most generic ones (in
mathematical terms) for which such a representation can (conjecturally) exist,
since a slight relaxation of the definition leads to a problem that is Graph
Isomorphism Complete.


  Studying the vocalisations of wild animals can be a challenge due to the
limitations of traditional computational methods, which often are
time-consuming and lack reproducibility. Here, I present pykanto, a new
software package that provides a set of tools to build, manage, and explore
large sound databases. It can automatically find discrete units in animal
vocalisations, perform semi-supervised labelling of individual repertoires with
a new interactive web app, and feed data to deep learning models to study
things like individual signatures and acoustic similarity between individuals
and populations. To demonstrate its capabilities, I put the library to the test
on the vocalisations of male great tits in Wytham Woods, near Oxford, UK. The
results show that the identities of individual birds can be accurately
determined from their songs and that the use of pykanto improves the efficiency
and reproducibility of the process.


  Intrinsically disordered proteins participate in many biological processes by
folding upon binding with other proteins. However, coupled folding and binding
processes are not well understood from an atomistic point of view. One of the
main questions is whether folding occurs prior to or after binding. Here we use
a novel unbiased high-throughput adaptive sampling approach to reconstruct the
binding and folding between the disordered transactivation domain of
\mbox{c-Myb} and the KIX domain of the CREB-binding protein. The reconstructed
long-term dynamical process highlights the binding of a short stretch of amino
acids on \mbox{c-Myb} as a folded $\alpha$-helix. Leucine residues, specially
Leu298 to Leu302, establish initial native contacts that prime the binding and
folding of the rest of the peptide, with a mixture of conformational selection
on the N-terminal region with an induced fit of the C-terminal.


  The study of cultural evolution benefits from detailed analysis of cultural
transmission in specific human domains. Chess provides a platform for
understanding the transmission of knowledge due to its active community of
players, precise behaviors, and long-term records of high-quality data. In this
paper, we perform an analysis of chess in the context of cultural evolution,
describing multiple cultural factors that affect move choice. We then build a
population-level statistical model of move choice in chess, based on the
Dirichlet-multinomial likelihood, to analyze cultural transmission over decades
of recorded games played by leading players. For moves made in specific
positions, we evaluate the relative effects of frequency-dependent bias,
success bias, and prestige bias on the dynamics of move frequencies. We observe
that negative frequency-dependent bias plays a role in the dynamics of certain
moves, and that other moves are compatible with transmission under prestige
bias or success bias. These apparent biases may reflect recent changes, namely
the introduction of computer chess engines and online tournament broadcasts.
Our analysis of chess provides insights into broader questions concerning how
social learning biases affect cultural evolution.


  The dysregulation of transcripts is characterized as one of the main
mechanisms in tumor pathogenesis. The recent discovery developed a new
hypothesis, competitive endogenous RNAs (ceRNAs), which could regulate other
RNA transcripts via competing for their shared miRNAs. The interaction of
elements in ceRNAs network was involved in a large range of biological
reactions and facilitate to cancer progression. In this study, we performed a
comprehensive investigation on the regulatory mechanisms and functional roles
of ceRNAs in prostate cancer (PCa) and constructed a ceRNAs network which could
possess potential value in patient prognosis and be evaluated as therapeutic
targets for PCa.


  Let $G$ be an undirected graph. We say that $G$ contains a ladder of length
$k$ if the $2 \times (k+1)$ grid graph is an induced subgraph of $G$ that is
only connected to the rest of $G$ via its four cornerpoints. We prove that if
all the ladders contained in $G$ are reduced to length 4, the treewidth remains
unchanged (and that this bound is tight). Our result indicates that, when
computing the treewidth of a graph, long ladders can simply be reduced, and
that minimal forbidden minors for bounded treewidth graphs cannot contain long
ladders. Our result also settles an open problem from algorithmic
phylogenetics: the common chain reduction rule, used to simplify the comparison
of two evolutionary trees, is treewidth-preserving in the display graph of the
two trees.


  Learning meaningful protein representation is important for a variety of
biological downstream tasks such as structure-based drug design. Having
witnessed the success of protein sequence pretraining, pretraining for
structural data which is more informative has become a promising research
topic. However, there are three major challenges facing protein structure
pretraining: insufficient sample diversity, physically unrealistic modeling,
and the lack of protein-specific pretext tasks. To try to address these
challenges, we present the 3D Geometric Pretraining. In this paper, we propose
a unified framework for protein pretraining and a 3D geometric-based,
data-efficient, and protein-specific pretext task: RefineDiff (Refine the
Diffused Protein Structure Decoy). After pretraining our geometric-aware model
with this task on limited data(less than 1% of SOTA models), we obtained
informative protein representations that can achieve comparable performance for
various downstream tasks.


  Drug discovery is a complex process with a large molecular space to be
considered. By constraining the search space, the fragment-based drug design is
an approach that can effectively sample the chemical space of interest. Here we
propose a framework of Recurrent Neural Networks (RNN) with an attention model
to sample the chemical space of organophosphorus molecules using the
fragment-based approach. The framework is trained with a ZINC dataset that is
screened for high druglikeness scores. The goal is to predict molecules with
similar biological action modes as organophosphorus pesticides or chemical
warfare agents yet less toxic to humans. The generated molecules contain a
starting fragment of PO2F but have a bulky hydrocarbon side chain limiting its
binding effectiveness to the targeted protein.


  Many reaction-diffusion models produce travelling wave solutions that can be
interpreted as waves of invasion in biological scenarios such as wound healing
or tumour growth. These partial differential equation models have since been
adapted to describe the interactions between cells and extracellular matrix
(ECM), using a variety of different underlying assumptions. In this work, we
derive a system of reaction-diffusion equations, with cross-species
density-dependent diffusion, by coarse-graining an agent-based, volume-filling
model of cell invasion into ECM. We study the resulting travelling wave
solutions both numerically and analytically across various parameter regimes.
Subsequently, we perform a systematic comparison between the behaviours
observed in this model and those predicted by simpler models in the literature
that do not take into account volume-filling effects in the same way. Our study
justifies the use of some of these simpler, more analytically tractable models
in reproducing the qualitative properties of the solutions in some parameter
regimes, but it also reveals some interesting properties arising from the
introduction of cell and ECM volume-filling effects, where standard model
simplifications might not be appropriate.


  The study of common heritability, or co-heritability, among multiple traits
has been widely established in quantitative and molecular genetics. However, in
bacteria, genome-based estimation of heritability has only been considered very
recently and no methods are currently available for considering
co-heritability. Here we introduce such a method and demonstrate its usefulness
by multi-trait analyses of the three major human pathogens \textit{Escherichia
coli}, \textit{Neisseria gonorrhoeae} and \textit{Streprococcus pneumoniae}. We
anticipate that the increased availability of high-throughput genomic and
phenotypic screens of bacterial populations will spawn ample future
opportunities to understand the common molecular basis of different traits in
bacteria.


  In their well-known textbook (Vincent & Brown, 2005), Vincent and Brown
suggested an attractive approach for studying evolutionary dynamics of
populations that are heterogeneous with respect to some strategy that affects
the fitness of individuals in the population. The authors developed a theory,
whose goal was to expand the applicability of mathematical models of population
dynamics by including dynamics of an evolving heritable phenotype trait subject
to natural selection. The authors studied both the case of evolution of
individual traits and of mean traits in the population (or species) and the
dynamics of total population size. The authors consider the developed approach
as (more or less) universally applicable to models with any fitness function
and any initial distribution of strategies, which is symmetric and has small
variance. Here it was shown that the scope of the approach proposed by Vincent
& Brown is unfortunately much more limited. I show that the approach gives
exact results only if the population dynamics linearly depends on the trait;
examples where the approach is incorrect are given.


  An expression quantitative trait locus (eQTL) is a chromosomal region where
genetic variants are associated with the expression levels of certain genes
that can be both nearby or distant. The identifications of eQTLs for different
tissues, cell types, and contexts have led to better understanding of the
dynamic regulations of gene expressions and implications of functional genes
and variants for complex traits and diseases. Although most eQTL studies to
date have been performed on data collected from bulk tissues, recent studies
have demonstrated the importance of cell-type-specific and context-dependent
gene regulations in biological processes and disease mechanisms. In this
review, we discuss statistical methods that have been developed to enable the
detections of cell-type-specific and context-dependent eQTLs from bulk tissues,
purified cell types, and single cells. We also discuss the limitations of the
current methods and future research opportunities.


  In this paper, we use the biological domain knowledge incorporated into
stochastic models for ab initio RNA secondary-structure prediction to improve
the state of the art in joint compression of RNA sequence and structure data
(Liu et al., BMC Bioinformatics, 2008). Moreover, we show that, conversely,
compression ratio can serve as a cheap and robust proxy for comparing the
prediction quality of different stochastic models, which may help guide the
search for better RNA structure prediction models.
  Our results build on expert stochastic context-free grammar models of RNA
secondary structures (Dowell & Eddy, BMC Bioinformatics, 2004; Nebel & Scheid,
Theory in Biosciences, 2011) combined with different (static and adaptive)
models for rule probabilities and arithmetic coding. We provide a prototype
implementation and an extensive empirical evaluation, where we illustrate how
grammar features and probability models affect compression ratios.


  The Escherichia coli chemotaxis signaling pathway has served as a model
system for studying the adaptive sensing of environmental signals by large
protein complexes. The chemoreceptors control the kinase activity of CheA in
response to the extracellular ligand concentration and adapt across a wide
concentration range by undergoing methylation and demethylation. Methylation
shifts the kinase response curve by orders of magnitude in ligand concentration
while incurring a much smaller change in the ligand binding curve. Here, we
show that this asymmetric shift in binding and kinase response is inconsistent
with equilibrium allosteric models regardless of parameter choices. To resolve
this inconsistency, we present a nonequilibrium allosteric model that
explicitly includes the dissipative reaction cycles driven by ATP hydrolysis.
The model successfully explains all existing measurements for both aspartate
and serine receptors. Our results suggest that while ligand binding controls
the equilibrium balance between the ON and OFF states of the kinase, receptor
methylation modulates the kinetic properties (e.g., the phosphorylation rate)
of the ON state. Furthermore, sufficient energy dissipation is necessary for
maintaining and enhancing the sensitivity range and amplitude of the kinase
response. We demonstrate that the nonequilibrium allosteric model is broadly
applicable to other sensor-kinase systems by successfully fitting previously
unexplained data from the DosP bacterial oxygen-sensing system. Overall, this
work provides a new perspective on cooperative sensing by large protein
complexes and opens up new research directions for understanding their
microscopic mechanisms through simultaneous measurements and modeling of ligand
binding and downstream responses.


  Background: Seasonal influenza causes a substantial burden on healthcare
services over the winter period when these systems are already under pressure.
Policies during the COVID-19 pandemic supressed the transmission of season
influenza, making the timing and magnitude of a potential resurgence difficult
to predict.
  Methods: We developed a hierarchical generalised additive model (GAM) for the
short-term forecasting of hospital admissions with a positive test for the
influenza virus sub-regionally across England. The model incorporates a
multi-level structure of spatio-temporal splines, weekly seasonality, and
spatial correlation. Using multiple performance metrics including interval
score, coverage, bias, and median absolute error, the predictive performance is
evaluated for the 2022/23 seasonal wave. Performance is measured against an
autoregressive integrated moving average (ARIMA) time series model.
  Results: The GAM method outperformed the ARIMA model across scoring rules at
both high and low-level geographies, and across the different phases of the
epidemic wave including the turning point. The performance of the GAM with a
14-day forecast horizon was comparable in error to the ARIMA at 7 days. The
performance of the GAM is found to be most sensitive to the flexibility of the
smoothing function that measures the national epidemic trend.
  Interpretation: This study introduces a novel approach to short-term
forecasting of hospital admissions with influenza using hierarchical, spatial,
and temporal components. The model is data-driven and practical to deploy using
information realistically available at time of prediction, addressing key
limitations of epidemic forecasting approaches. This model was used across the
winter for healthcare operational planning by the UK Health Security Agency and
the National Health Service in England.