  Polymer models play the special role of elucidating the elementary features
describing the physics of long molecules and become essential to interpret the
measurements of their magnitudes. In this work the end-to-end distance of an
extensible discrete worm-like chain polymer as a function of the applied force
has been calculated both numerically and analytically, the latter as an
effective approximation. The numerical evaluation uses the Transfer Matrix
formalism to obtain an exact calculation of the partition function, while the
analytic derivations generalize the simple phenomenological formulas largely
used up to now. The obtained formulas are simple enough to be implemented in
the fit analysis of experimental data of semi-flexible extensible polymers,
with the result that the elastic parameters obtained are compatible with
previous measurements, and more, their accuracy strongly improves in a large
range of chain extensibility.


  Motivation: Genome-wide association studies (GWAS) have successfully
identified thousands of genetic risk loci for complex traits and diseases. Most
of these GWAS loci lie in regulatory regions of the genome and the gene through
which each GWAS risk locus exerts its effects is not always clear. Many
computational methods utilizing biological data sources have been proposed to
identify putative casual genes at GWAS loci; however, these methods can be
improved upon. Results: We present the Relations-Maximization Method, a dense
module searching method to identify putative causal genes at GWAS loci through
the generation of candidate sub-networks derived by integrating association
signals from GWAS data into the gene co-regulation network. We employ our
method in a chronic obstructive pulmonary disease GWAS. We perform an
extensive, comparative study of Relations-Maximization Method's performance
against well-established baselines.


  Descriptor generation methods using latent representations of
encoder$-$decoder (ED) models with SMILES as input are useful because of the
continuity of descriptor and restorability to the structure. However, it is not
clear how the structure is recognized in the learning progress of ED models. In
this work, we created ED models of various learning progress and investigated
the relationship between structural information and learning progress. We
showed that compound substructures were learned early in ED models by
monitoring the accuracy of downstream tasks and input$-$output substructure
similarity using substructure$-$based descriptors, which suggests that existing
evaluation methods based on the accuracy of downstream tasks may not be
sensitive enough to evaluate the performance of ED models with SMILES as
descriptor generation methods. On the other hand, we showed that structure
restoration was time$-$consuming, and in particular, insufficient learning led
to the estimation of a larger structure than the actual one. It can be inferred
that determining the endpoint of the structure is a difficult task for the
model. To our knowledge, this is the first study to link the learning progress
of SMILES by ED model to chemical structures for a wide range of chemicals.


  Biological functions of RNAs are determined by their three-dimensional (3D)
structures. Thus, given the limited number of experimentally determined RNA
structures, the prediction of RNA structures will facilitate elucidating RNA
functions and RNA-targeted drug discovery, but remains a challenging task. In
this work, we propose a Graph Neural Network (GNN)-based scoring function
trained only with the atomic types and coordinates on limited solved RNA 3D
structures for distinguishing accurate structural models. The proposed
Physics-aware Multiplex Graph Neural Network (PaxNet) separately models the
local and non-local interactions inspired by molecular mechanics. Furthermore,
PaxNet contains an attention-based fusion module that learns the individual
contribution of each interaction type for the final prediction. We rigorously
evaluate the performance of PaxNet on two benchmarks and compare it with
several state-of-the-art baselines. The results show that PaxNet significantly
outperforms all the baselines overall, and demonstrate the potential of PaxNet
for improving the 3D structure modeling of RNA and other macromolecules. Our
code is available at https://github.com/zetayue/Physics-aware-Multiplex-GNN.


  Geometric deep learning has been revolutionizing the molecular modeling
field. Despite the state-of-the-art neural network models are approaching ab
initio accuracy for molecular property prediction, their applications, such as
drug discovery and molecular dynamics (MD) simulation, have been hindered by
insufficient utilization of geometric information and high computational costs.
Here we propose an equivariant geometry-enhanced graph neural network called
ViSNet, which elegantly extracts geometric features and efficiently models
molecular structures with low computational costs. Our proposed ViSNet
outperforms state-of-the-art approaches on multiple MD benchmarks, including
MD17, revised MD17 and MD22, and achieves excellent chemical property
prediction on QM9 and Molecule3D datasets. Additionally, ViSNet achieved the
top winners of PCQM4Mv2 track in the OGB-LCS@NeurIPS2022 competition.
Furthermore, through a series of simulations and case studies, ViSNet can
efficiently explore the conformational space and provide reasonable
interpretability to map geometric representations to molecular structures.


  A challenge in global change biology is to predict how species will respond
to future environmental change and to manage these responses. To make such
predictions and management actions robust to novel futures, we need to
accurately characterize how organisms experience their environments and the
biological mechanisms by which they respond. All organisms are
thermodynamically connected to their environments through the exchange of heat
and water at fine spatial and temporal scales and this exchange can be captured
with biophysical models. Although mechanistic models based on biophysical
ecology have a long history of development and application, their use in global
change biology remains limited despite their enormous promise and increasingly
accessible software. We contend that greater understanding and training in the
theory and methods of biophysical models is vital to expand their application.
Our review shows how biophysical models can be implemented to understand and
predict climate change impacts on species' behavior, phenology, survival,
distribution, and abundance. We illustrate the types of outputs that can be
generated, and the data inputs required for different implementations. Examples
range from simple calculations of body temperature to more complex analyses of
species' distribution limits based on projected energy and water balances,
accounting for behavior and phenology. We outline challenges that currently
limit the widespread application of biophysical models. We discuss progress and
future developments that could allow these models to be applied to many species
across large spatial extents and timeframes. We highlight how biophysical
models are uniquely suited to solve global change biology problems that involve
predicting and interpreting responses to environmental variability and
extremes, multiple or shifting constraints, and novel abiotic or biotic
environments.


  How did written works evolve, disappear or survive down through the ages? In
this paper, we propose a unified, formal framework for two fundamental
questions in the study of the transmission of texts: how much was lost or
preserved from all works of the past, and why do their genealogies (their
``phylogenetic trees'') present the very peculiar shapes that we observe or,
more precisely, reconstruct? We argue here that these questions share
similarities to those encountered in evolutionary biology, and can be described
in terms of ``genetic'' drift and ``natural'' selection. Through agent-based
models, we show that such properties as have been observed by philologists
since the 1800s can be simulated, and confronted to data gathered for ancient
and medieval texts across Europe, in order to obtain plausible estimations of
the number of works and manuscripts that existed and were lost.


  Structural Maintenance of Chromosomes, SMCs, proteins have long rod like
structures immersed in water. Here we use our hydroanalytic methods based on
amino acid sequences to discuss their dynamics at multiple length scales
identified by evolution. The length scales are 10 to 100 times longer than used
in normal studies of sequence evolution. Their hydropathic profiles exhibit
many features unique to their structure and function.


  Studies of nanoscale biological transport often use a mean-field
approximation that is exact only when the system is at equilibrium and there
are no interactions between particles on different sites in the network. We
explore the limitations of this approximation to describe many-particle
transport in the context of enzyme function and biological transport networks.
Our focus is on three bioenergetic networks: a linear electron transfer chain
(as found in bacterial nanowires), a redox-coupled proton pump (as in complex
IV of respiration), and a near reversible electron bifurcation network (as in
complex III of respiration and other recently discovered structures). Away from
equilibrium and with typical site-site interactions, we find that the
mean-field approximation adequately describes linear transport chains. However,
the mean-field approximation fails catastrophically to describe
energy-transducing systems, as in the redox coupled proton pump and reversible
electron bifurcation reactions. The mean-field approximation fails to capture
the essential correlations that are needed to prevent slippage events and to
produce efficient energy transduction.


  In the context of SARS-CoV-2 pandemic, mathematical modelling has played a
fundamental role for making forecasts, simulating scenarios and evaluating the
impact of preventive political, social and pharmaceutical measures. Optimal
control theory can be a useful tool based on solid mathematical bases to plan
the vaccination campaign in the direction of eradicating the pandemic as fast
as possible. The aim of this work is to explore the optimal prioritisation
order for planning vaccination campaigns able to achieve specific goals, as the
reduction of the amount of infected, deceased and hospitalized in a fixed time
frame, among age classes. For this purpose, we introduce an age stratified
SIR-like epidemic compartmental model settled in an abstract framework for
modelling two-doses vaccination campaigns and conceived with the description of
COVID19 disease. Overall, we formalize an optimal control framework adopting
the model as state problem by acting on the administrations of vaccine-doses.
An extensive campaign of numerical tests, featured in the Italian scenario and
calibrated on available data from Dipartimento di Protezione Civile Italiana,
shows that the presented framework can be a valuable tool to support the
planning of vaccination campaigns minimizing specific goals.


  We study the coevolution of network structure and signaling behavior. We
model agents who can preferentially associate with others in a dynamic network
while they also learn to play a simple sender-receiver game. We have four major
findings. First, signaling interactions in dynamic networks are sufficient to
cause the endogenous formation of distinct signaling groups, even in an
initially homogeneous population. Second, dynamic networks allow the emergence
of novel {\em hybrid} signaling groups that do not converge on a single common
signaling system but are instead composed of different yet complementary
signaling strategies. We show that the presence of these hybrid groups promotes
stable diversity in signaling among other groups in the population. Third, we
find important distinctions in information processing capacity of different
groups: hybrid groups diffuse information more quickly initially but at the
cost of taking longer to reach all group members. Fourth, our findings pertain
to all common interest signaling games, are robust across many parameters, and
mitigate known problems of inefficient communication.


  GNN-based methods have achieved excellent results as a mainstream task in
drug response prediction tasks in recent years. Traditional GNN methods use
only the atoms in a drug molecule as nodes to obtain the representation of the
molecular graph through node information passing, whereas the method using the
transformer can only extract information about the nodes. However, the covalent
bonding and chirality of a drug molecule have a great influence on the
pharmacological properties of the molecule, and these information are implied
in the chemical bonds formed by the edges between the atoms. In addition, CNN
methods for modelling cell lines genomics sequences can only perceive local
rather than global information about the sequence. In order to solve the above
problems, we propose the decoupled dual transformer structure with edge
embedded for drug respond prediction (TransEDRP), which is used for the
representation of cell line genomics and drug respectively. For the drug
branch, we encoded the chemical bond information within the molecule as the
embedding of the edge in the molecular graph, extracted the global structural
and biochemical information of the drug molecule using graph transformer. For
the branch of cell lines genomics, we use the multi-headed attention mechanism
to globally represent the genomics sequence. Finally, the drug and genomics
branches are fused to predict IC50 values through the transformer layer and the
fully connected layer, which two branches are different modalities. Extensive
experiments have shown that our method is better than the current mainstream
approach in all evaluation indicators.


  The COVID-19 pandemic has caused over 6.4 million registered deaths to date
and has had a profound impact on economic activity. Here, we study the
interaction of transmission, mortality, and the economy during the SARS-CoV-2
pandemic from January 2020 to December 2022 across 25 European countries. We
adopt a Bayesian Mixed Effects model with auto-regressive terms. We find that
increases in disease transmission intensity decreases Gross domestic product
(GDP) and increases daily excess deaths, with a longer lasting impact on excess
deaths in comparison to GDP, which recovers more rapidly. Broadly, our results
reinforce the intuitive phenomenon that significant economic activity arises
from diverse person-to-person interactions. We report on the effectiveness of
non-pharmaceutical interventions (NPIs) on transmission intensity, excess
deaths, and changes in GDP, and resulting implications for policy makers. Our
results highlight a complex cost-benefit trade off from individual NPIs. For
example, banning international travel increases GDP and reduces excess deaths.
We consider country random effects and their associations with excess changes
in GDP and excess deaths. For example, more developed countries in Europe
typically had more cautious approaches to the COVID-19 pandemic, prioritising
healthcare, and excess deaths over economic performance. Long term economic
impairments are not fully captured by our model, as well as long term disease
effects (Long Covid). Our results highlight that the impact of disease on a
country is complex and multifaceted, and simple heuristic conclusions to
extract the best outcome from the economy and disease burden are challenging.


  We study SIS epidemic spreading models under population dispersal on
multi-layer networks. We consider a patchy environment in which each patch
comprises individuals belonging to different classes. Individuals disperse to
other patches on a multi-layer network in which each layer corresponds to a
class. The dispersal on each layer is modeled by a Continuous Time Markov Chain
(CTMC). At each time, individuals disperse according to their CTMC and
subsequently interact with the local individuals in the patch according to an
SIS model. We establish the existence of various equilibria under different
parameter regimes and establish their (almost) global asymptotic stability
using Lyapunov techniques. We also derive simple conditions that highlight the
influence of the multi-layer network on the stability of these equilibria. For
this model, we study optimal intervention strategies using a convex
optimization framework. Finally, we numerically illustrate the influence of the
multi-layer network structure and the effectiveness of the optimal intervention
strategies.


  The preferred conditions for evolutionary innovation represent a fundamental
question, but little is known experimentally or theoretically. In this study,
we focused on the potential role of polyploidy in the evolution of novel
traits. We proposed a simple model and demonstrated that the evolutionary rate
of polyploids is similar to more much slower than that of haploids under
neutral selection or during gradual evolution. However, experiments using
polyploid cyanobacteria demonstrated that the probability of achieving
antibiotic resistance increased with the number of chromosomes and implied an
optimal number of chromosomes. Then, we investigated the dynamics of the same
model on a fitness landscape in which cells should jump over a lethal valley to
increase their fitness. The evolutionary rate could be increased in polyploidy,
and the optimal number of chromosomes was identified. Further, we proposed that
the optimization for evolutionary innovation might determine the number of
chromosomes in polyploid bacteria.


  There is a growing interest in modeling and analyzing the spread of diseases
like the SARS-CoV-2 infection using stochastic models. These models are
typically analyzed quantitatively and are not often subject to validation using
formal verification approaches, nor leverage policy syntheses and analysis
techniques developed in formal verification. In this paper, we take a Markovian
stochastic model for the spread of a SARSCoV-2-like infection. A state of this
model represents the number of subjects in different health conditions. The
considered model considers the different parameters that may have an impact on
the spread of the disease and exposes the various decision variables that can
be used to control it. We show that the modeling of the problem within
state-of-the-art model checkers is feasible and it opens several opportunities.
However, there are severe limitations due to i) the espressivity of the
existing stochastic model checkers on one side, and ii) the size of the
resulting Markovian model even for small population sizes.


  Hybridization events complicate the accurate reconstruction of phylogenies,
as they lead to patterns of genetic heritability that are unexpected under
traditional, bifurcating models of species trees. This has led to the
development of methods to infer these varied hybridization events, both methods
that reconstruct networks directly, and summary methods that predict individual
hybridization events. However, a lack of empirical comparisons between methods
- especially pertaining to large networks with varied hybridization scenarios -
hinders their practical use. Here, we provide a comprehensive review of popular
summary methods: TICR, MSCquartets, HyDe, Patterson's D-Statistic (ABBA-BABA),
D3, and Dp. TICR and MSCquartets are based on quartet concordance factors
gathered from gene tree topologies and Patterson's D-Statistic, D3, and Dp use
site pattern frequencies to identify hybridization events. We then use
simulated data to address questions of method accuracy and ideal use scenarios
by testing methods against complex networks which depict gene flow events that
differ in depth (timing), quantity (single vs. multiple, overlapping
hybridizations), and rate of gene flow. We find that deeper or multiple
hybridization events may introduce noise and weaken the signal of
hybridization, leading to higher false negative rates across methods. Despite
some forms of hybridization eluding quartet-based detection methods,
MSCquartets displays high precision in most scenarios. While HyDe results in
high false negative rates when tested on hybridizations involving ghost
lineages, HyDe is the only method to be able to separate hybrid vs parent
signals. Lastly, we test the methods on ultraconserved elements from the bee
subfamily Nomiinae, finding the possibility of hybridization events between
clades which correspond to regions of poor support in the species tree
estimated in the original study.


  The human {\it ether-a-go-go} (hERG) potassium channel (K$_\text{v}11.1$)
plays a critical role in mediating cardiac action potential. The blockade of
this ion channel can potentially lead fatal disorder and/or long QT syndrome.
Many drugs have been withdrawn because of their serious hERG-cardiotoxicity. It
is crucial to assess the hERG blockade activity in the early stage of drug
discovery. We are particularly interested in the hERG-cardiotoxicity of
compounds collected in the DrugBank database considering that many DrugBank
compounds have been approved for therapeutic treatments or have high potential
to become drugs. Machine learning-based in silico tools offer a rapid and
economical platform to virtually screen DrugBank compounds. We design accurate
and robust classifiers for blockers/non-blockers and then build regressors to
quantitatively analyze the binding potency of the DrugBank compounds on the
hERG channel. Molecular sequences are embedded with two natural language
processing (NPL) methods, namely, autoencoder and transformer. Complementary
three-dimensional (3D) molecular structures are embedded with two advanced
mathematical approaches, i.e., topological Laplacians and algebraic graphs.
With our state-of-the-art tools, we reveal that 227 out of the 8641 DrugBank
compounds are potential hERG blockers, suggesting serious drug safety problems.
Our predictions provide guidance for the further experimental interrogation of
DrugBank compounds' hERG-cardiotoxicity .


  Mapper, a topological algorithm, is frequently used as an exploratory tool to
build a graphical representation of data. This representation can help to gain
a better understanding of the intrinsic shape of high-dimensional genomic data
and to retain information that may be lost using standard dimension-reduction
algorithms. We propose a novel workflow to process and analyze RNA-seq data
from tumor and healthy subjects integrating Mapper and differential gene
expression. Precisely, we show that a Gaussian mixture approximation method can
be used to produce graphical structures that successfully separate tumor and
healthy subjects, and produce two subgroups of tumor subjects. A further
analysis using DESeq2, a popular tool for the detection of differentially
expressed genes, shows that these two subgroups of tumor cells bear two
distinct gene regulations, suggesting two discrete paths for forming lung
cancer, which could not be highlighted by other popular clustering methods,
including t-SNE. Although Mapper shows promise in analyzing high-dimensional
data, building tools to statistically analyze Mapper graphical structures is
limited in the existing literature. In this paper, we develop a scoring method
using heat kernel signatures that provides an empirical setting for statistical
inferences such as hypothesis testing, sensitivity analysis, and correlation
analysis.


  Morphogenesis, the establishment and repair of emergent complex anatomy by
groups of cells, is a fascinating and biomedically-relevant problem. One of its
most fascinating aspects is that a developing embryo can reliably recover from
disturbances, such as splitting into twins. While this reliability implies some
type of goal-seeking error minimization over a morphogenic field, there are
many gaps with respect to detailed, constructive models of such a process being
used to implement the collective intelligence of cellular swarms. We describe a
closed-loop negative-feedback system for creating reaction-diffusion (RD)
patterns with high reliability. It uses a cellular automaton to characterize a
morphogen pattern, then compares it to a goal and adjusts accordingly,
providing a framework for modeling anatomical homeostasis and robust generation
of target morphologies. Specifically, we create a RD pattern with N
repetitions, where N is easily changeable. Furthermore, the individual
repetitions of the RD pattern can be easily stretched or shrunk under genetic
control to create, e.g., some morphological features larger than others.
Finally, the cellular automaton uses a computation wave that scans the
morphogen pattern unidirectionally to characterize the features that the
negative feedback then controls. By taking advantage of a prior process
asymmetrically establishing planar polarity (e.g., head vs. tail), our
automaton is greatly simplified. This work contributes to the exciting effort
of understanding design principles of morphological computation, which can be
used to understand evolved developmental mechanisms, manipulate them in
regenerative medicine settings, or embed a degree of synthetic intelligence
into novel bioengineered constructs.


  Genomic signal processing has been used successfully in bioinformatics to
analyze biomolecular sequences and gain varied insights into DNA structure,
gene organization, protein binding, sequence evolution, etc. But challenges
remain in finding the appropriate spectral representation of a biomolecular
sequence, especially when multiple variable-length sequences need to be handled
consistently. In this study, we address this challenge in the context of the
well-studied problem of classifying genomic sequences into different taxonomic
units (strain, phyla, order, etc.). We propose a novel technique that employs
signal processing in tandem with Gaussian mixture models to improve the
spectral representation of a sequence and subsequently the taxonomic
classification accuracies. The sequences are first transformed into spectra,
and projected to a subspace, where sequences belonging to different taxons are
better distinguishable. Our method outperforms a similar state-of-the-art
method on established benchmark datasets by an absolute margin of 6.06%
accuracy.


  Molecular property prediction is essential for drug discovery. In recent
years, deep learning methods have been introduced to this area and achieved
state-of-the-art performances. However, most of existing methods ignore the
intrinsic relations between molecular properties which can be utilized to
improve the performances of corresponding prediction tasks. In this paper, we
propose a new approach, namely Physics properties Enhanced Molecular Property
prediction (PEMP), to utilize relations between molecular properties revealed
by previous physics theory and physical chemistry studies. Specifically, we
enhance the training of the chemical and physiological property predictors with
related physics property prediction tasks. We design two different methods for
PEMP, respectively based on multi-task learning and transfer learning. Both
methods include a model-agnostic molecule representation module and a property
prediction module. In our implementation, we adopt both the state-of-the-art
molecule embedding models under the supervised learning paradigm and the
pretraining paradigm as the molecule representation module of PEMP,
respectively. Experimental results on public benchmark MoleculeNet show that
the proposed methods have the ability to outperform corresponding
state-of-the-art models.


  The physics of how molecules, organelles, and foreign objects move within
living cells has been extensively studied in organisms ranging from bacteria to
human cells. In mammalian cells, in particular, cellular vesicles move across
the cell using motor proteins that carry the vesicle down the cytoskeleton to
their destination. We have recently noted several similarities between the
motion of such vesicles and that in disordered, "glassy", systems, but it
remains unclear whether that is a general observation or something specific to
certain vesicles in one particular cell type. Here we follow the motion of
mitochondria, the organelles responsible for cell energy production, in several
mammalian cell types over timescales ranging from 50 ms up to 70 s. Qualitative
observations show that single mitochondria remain stalled, remaining within a
spatially limited region, for extended periods of time, before moving longer
distances relatively quickly. Analysing this motion quantitatively, we observe
a displacement distribution that is roughly Gaussian for shorter distances
($\lesssim$ 0.05 $\mu$m) but which exhibits exponentially decaying tails at
longer distances (up to 0.40 $\mu$m). We show that this behaviour is
well-described by a model originally developed to describe the motion in glassy
systems. These observations are extended to in total 3 different objects
(mitochondria, lysosomes and nano-sized beads enclosed in vesicles), 3
different mammalian cell types, from 2 different organisms (human and mouse).
We provide further evidence that supports glass-like characteristics of the
motion by showing a difference between the time it takes to move a longer
distance for the first time and subsequent times, as well as a weak ergodicity
breaking of the motion. Overall, we demonstrate the ubiquity of glass-like
motion in mammalian cells, providing a different perspective on intracellular
motion.


  Gradients of probabilistic model likelihoods with respect to their parameters
are essential for modern computational statistics and machine learning. These
calculations are readily available for arbitrary models via automatic
differentiation implemented in general-purpose machine-learning libraries such
as TensorFlow and PyTorch. Although these libraries are highly optimized, it is
not clear if their general-purpose nature will limit their algorithmic
complexity or implementation speed for the phylogenetic case compared to
phylogenetics-specific code. In this paper, we compare six gradient
implementations of the phylogenetic likelihood functions, in isolation and also
as part of a variational inference procedure. We find that although automatic
differentiation can scale approximately linearly in tree size, it is much
slower than the carefully-implemented gradient calculation for tree likelihood
and ratio transformation operations. We conclude that a mixed approach
combining phylogenetic libraries with machine learning libraries will provide
the optimal combination of speed and model flexibility moving forward.


  An approach is introduced for comparing the estimated states of stochastic
compartmental models for an epidemic or biological process with analytically
obtained solutions from the corresponding system of ordinary differential
equations (ODEs). Positive integer valued samples from a stochastic model are
generated numerically at discrete time intervals using either the Reed-Frost
chain Binomial or Gillespie algorithm. The simulated distribution of
realisations is compared with an exact solution obtained analytically from the
ODE model. Using this novel methodology this work demonstrates it is feasible
to check that the realisations from the stochastic compartmental model adhere
to the ODE model they represent. There is no requirement for the model to be in
any particular state or limit. These techniques are developed using the
stochastic compartmental model for a susceptible-infected-recovered (SIR)
epidemic process. The Lotka-Volterra model is then used as an example of the
generality of the principles developed here. This approach presents a way of
testing/benchmarking the numerical solutions of stochastic compartmental
models, e.g. using unit tests, to check that the computer code along with its
corresponding algorithm adheres to the underlying ODE model.


  We introduce a biologically natural, mathematically tractable model of random
phylogenetic network to describe evolution in the presence of hybridization.
One of the features of this model is that the hybridization rate of the
lineages correlates negatively with their phylogenetic distance. We give
formulas / characterizations for quantities of biological interest that make
them straightforward to compute in practice. We show that the appropriately
rescaled network, seen as a metric space, converges to the Brownian continuum
random tree, and that the uniformly rooted network has a local weak limit,
which we describe explicitly.


  Polymers consisting of more than one type of monomer, known as copolymers,
are vital to both living and synthetic systems. Copolymerisation has been
studied theoretically in a number of contexts, often by considering a Markov
process in which monomers are added or removed from the growing tip of a long
copolymer. To date, the analysis of the most general models of this class has
necessitated simulation. We present a general method for analysing such
processes without resorting to simulation. Our method can be applied to models
with an arbitrary network of sub-steps prior to addition or removal of a
monomer, including non-equilibrium kinetic proofreading cycles. Moreover, the
approach allows for a dependency of addition and removal reactions on the
neighbouring site in the copolymer, and thermodynamically self-consistent
models in which all steps are assumed to be microscopically reversible. Using
our approach, thermodynamic quantities such as chemical work; kinetic
quantities such as time taken to grow; and statistical quantities such as the
distribution of monomer types in the growing copolymer can be derived either
analytically or numerically directly from the model definition.


  The field of geometric deep learning has had a profound impact on the
development of innovative and powerful graph neural network architectures.
Disciplines such as computer vision and computational biology have benefited
significantly from such methodological advances, which has led to breakthroughs
in scientific domains such as protein structure prediction and design. In this
work, we introduce GCPNet, a new geometry-complete, SE(3)-equivariant graph
neural network designed for 3D molecular graph representation learning.
Rigorous experiments across four distinct geometric tasks demonstrate that
GCPNet's predictions (1) for protein-ligand binding affinity achieve a
statistically significant correlation of 0.608, more than 5% greater than
current state-of-the-art methods; (2) for protein structure ranking achieve
statistically significant target-local and dataset-global correlations of 0.616
and 0.871, respectively; (3) for Newtownian many-body systems modeling achieve
a task-averaged mean squared error less than 0.01, more than 15% better than
current methods; and (4) for molecular chirality recognition achieve a
state-of-the-art prediction accuracy of 98.7%, better than any other machine
learning method to date. The source code, data, and instructions to train new
models or reproduce our results are freely available at
https://github.com/BioinfoMachineLearning/GCPNet.


  The BIOSCAN project, led by the International Barcode of Life Consortium,
seeks to study changes in biodiversity on a global scale. One component of the
project is focused on studying the species interaction and dynamics of all
insects. In addition to genetically barcoding insects, over 1.5 million images
per year will be collected, each needing taxonomic classification. With the
immense volume of incoming images, relying solely on expert taxonomists to
label the images would be impossible; however, artificial intelligence and
computer vision technology may offer a viable high-throughput solution.
Additional tasks including manually weighing individual insects to determine
biomass, remain tedious and costly. Here again, computer vision may offer an
efficient and compelling alternative. While the use of computer vision methods
is appealing for addressing these problems, significant challenges resulting
from biological factors present themselves. These challenges are formulated in
the context of machine learning in this paper.


  Genomic complexity is a growing field of evolution, with case studies for
comparative evolutionary analyses in model and emerging non-model systems.
Understanding complexity and the functional components of the genome is an
untapped wealth of knowledge ripe for exploration. With the "remarkable lack of
correspondence" between genome size and complexity, there needs to be a way to
quantify complexity across organisms. In this study we use a set of complexity
metrics that allow for evaluation of changes in complexity using TranD. We
ascertain if complexity is increasing or decreasing across transcriptomes and
at what structural level, as complexity is varied. We define three metrics --
TpG, EpT, and EpG in this study to quantify the complexity of the transcriptome
that encapsulate the dynamics of alternative splicing. Here we compare
complexity metrics across 1) whole genome annotations, 2) a filtered subset of
orthologs, and 3) novel genes to elucidate the impacts of ortholog and novel
genes in transcriptome analysis. We also derive a metric from Hong et al.,
2006, Effective Exon Number (EEN), to compare the distribution of exon sizes
within transcripts against random expectations of uniform exon placement. EEN
accounts for differences in exon size, which is important because novel genes
differences in complexity for orthologs and whole transcriptome analyses are
biased towards low complexity genes with few exons and few alternative
transcripts. With our metric analyses, we are able to implement changes in
complexity across diverse lineages with greater precision and accuracy than
previous cross-species comparisons under ortholog conditioning. These analyses
represent a step forward toward whole transcriptome analysis in the emerging
field of non-model evolutionary genomics, with key insights for evolutionary
inference of complexity changes on deep timescales across the tree of life. We
suggest a means to quantify biases generated in ortholog calling and correct
complexity analysis for lineage-specific effects. With these metrics, we
directly assay the quantitative properties of newly formed lineage-specific
genes as they lower complexity in transcriptomes.


  Models that accurately predict properties based on chemical structure are
valuable tools in drug discovery. However, for many properties, public and
private training sets are typically small, and it is difficult for the models
to generalize well outside of the training data. Recently, large language
models have addressed this problem by using self-supervised pretraining on
large unlabeled datasets, followed by fine-tuning on smaller, labeled datasets.
In this paper, we report MolE, a molecular foundation model that adapts the
DeBERTa architecture to be used on molecular graphs together with a two-step
pretraining strategy. The first step of pretraining is a self-supervised
approach focused on learning chemical structures, and the second step is a
massive multi-task approach to learn biological information. We show that
fine-tuning pretrained MolE achieves state-of-the-art results on 9 of the 22
ADMET tasks included in the Therapeutic Data Commons.


  The Zika virus, is an emerging infectious disease causing severe
complications such as microcephaly in infants and Guillain Barre syndrome in
adults. There is no licensed vaccination or approved medicine to treat ZIKV
infection. Therefore, extensive research is being carried out to find compounds
that can be used effectively as therapeutic molecules to treat ZIKV infection.
Oregano, a commonly found herb in the Mediterranean region, has been used
predominantly for culinary purposes. The fact that the members of the Origanum
species are a storehouse of various bioactive compounds gives us a solid reason
to study compounds extracted from it for therapeutic purposes. In this study,
were retrieved 20 Flavonoids found in various Origanum species belonging to the
Mediterranean region from the PubChem database and pharmacological analysis
using SwissADME and Molecular docking using AutoDock Vina 4.0. were carried out
against the NS2B NS3 protease since it serves as an effective drug target owing
to its role in viral replication and immune evasion within the host. The best
hit compounds were subjected to MD simulation at 100 ns using Desmond
Schrodinger to analyze the molecule's stability. We observed Cirsiliol as the
best hit compound against the NS2B NS3 complex with a binding affinity of -8.5
kcal per mol. It also showed good stability during MD simulation. We recommend
the use of Cirsiliol for in vitro and in vivo studies for further investigation
concerning the ZIKA virus.


  Biomolecular binding kinetics including the association (kon) and
dissociation (koff) rates are critical parameters for therapeutic design of
small-molecule drugs, peptides and antibodies. Notably, drug molecule residence
time or dissociation rate has been shown to correlate with their efficacies
better than binding affinities. A wide range of modeling approaches including
quantitative structure-kinetic relationship models, Molecular Dynamics
simulations, enhanced sampling and Machine Learning have been developed to
explore biomolecular binding and dissociation mechanisms and predict binding
kinetic rates. Here, we review recent advances in computational modeling of
biomolecular binding kinetics, with an outlook for future improvements.


  Understanding the spread of diseases through complex networks is of great
interest where realistic, heterogeneous contact patterns play a crucial role in
the spread. Most works have focused on mean-field behavior -- quantifying how
contact patterns affect the emergence and stability of (meta)stable endemic
states in networks. On the other hand, much less is known about longer time
scale dynamics, such as disease extinction, whereby inherent process
stochasticity and contact heterogeneity interact to produce large fluctuations
that result in the spontaneous clearance of infection. Here we show that
heterogeneity in both susceptibility and infectiousness (incoming and outgoing
degree, respectively) has a non-trivial effect on extinction in directed
contact networks, both speeding-up and slowing-down extinction rates depending
on the relative proportion of such edges in a network, and on whether the
heterogeneities in the incoming and outgoing degrees are correlated or
anticorrelated. In particular, we show that weak anticorrelated heterogeneity
can increase the disease stability, whereas strong heterogeneity gives rise to
markedly different results for correlated and anticorrelated heterogeneous
networks. All analytical results are corroborated through various numerical
schemes including network Monte-Carlo simulations.


  High order structures (cavities and cliques) of the gene network of influenza
A virus reveal tight associations among viruses during evolution and are key
signals that indicate viral cross-species infection and cause pandemics. As
indicators for sensing the dynamic changes of viral genes, these higher order
structures have been the focus of attention in the field of virology. However,
the size of the viral gene network is usually huge, and searching these
structures in the networks introduces unacceptable delay. To mitigate this
issue, in this paper, we propose a simple-yet-effective model named HyperSearch
based on deep learning to search cavities in a computable complex network for
influenza virus genetics. Extensive experiments conducted on a public influenza
virus dataset demonstrate the effectiveness of HyperSearch over other advanced
deep-learning methods without any elaborated model crafting. Moreover,
HyperSearch can finish the search works in minutes while 0-1 programming takes
days. Since the proposed method is simple and easy to be transferred to other
complex networks, HyperSearch has the potential to facilitate the monitoring of
dynamic changes in viral genes and help humans keep up with the pace of virus
mutations.


  Nanopore sequencing generates noisy electrical signals that need to be
converted into a standard string of DNA nucleotide bases using a computational
step called basecalling. The accuracy and speed of basecalling have critical
implications for all later steps in genome analysis. Many researchers adopt
complex deep learning-based models to perform basecalling without considering
the compute demands of such models, which leads to slow, inefficient, and
memory-hungry basecallers. Therefore, there is a need to reduce the computation
and memory cost of basecalling while maintaining accuracy. Our goal is to
develop a comprehensive framework for creating deep learning-based basecallers
that provide high efficiency and performance. We introduce RUBICON, a framework
to develop hardware-optimized basecallers. RUBICON consists of two novel
machine-learning techniques that are specifically designed for basecalling.
First, we introduce the first quantization-aware basecalling neural
architecture search (QABAS) framework to specialize the basecalling neural
network architecture for a given hardware acceleration platform while jointly
exploring and finding the best bit-width precision for each neural network
layer. Second, we develop SkipClip, the first technique to remove the skip
connections present in modern basecallers to greatly reduce resource and
storage requirements without any loss in basecalling accuracy. We demonstrate
the benefits of RUBICON by developing RUBICALL, the first hardware-optimized
basecaller that performs fast and accurate basecalling. Compared to the fastest
state-of-the-art basecaller, RUBICALL provides a 3.96x speedup with 2.97%
higher accuracy. We show that RUBICON helps researchers develop
hardware-optimized basecallers that are superior to expert-designed models.


  Protein structure prediction is a critical problem linked to drug design,
mutation detection, and protein synthesis, among other applications. To this
end, evolutionary data has been used to build contact maps which are
traditionally minimized as energy functions via gradient descent based schemes
like the L-BFGS algorithm. In this paper we present what we call the
Alternating Metropolis-Hastings (AMH) algorithm, which (a) significantly
improves the performance of traditional MCMC methods, (b) is inherently
parallelizable allowing significant hardware acceleration using GPU, and (c)
can be integrated with the L-BFGS algorithm to improve its performance. The
algorithm shows an improvement in energy of found structures of 8.17% to 61.04%
(average 38.9%) over traditional MH and 0.53% to 17.75% (average 8.9%) over
traditional MH with intermittent noisy restarts, tested across 9 proteins from
recent CASP competitions. We go on to map the Alternating MH algorithm to a
GPGPU which improves sampling rate by 277x and improves simulation time to a
low energy protein prediction by 7.5x to 26.5x over CPU. We show that our
approach can be incorporated into state-of-the-art protein prediction pipelines
by applying it to both trRosetta2's energy function and the distogram component
of Alphafold1's energy function. Finally, we note that specially designed
probabilistic computers (or p-computers) can provide even better performance
than GPUs for MCMC algorithms like the one discussed here.


  Drug development is a wide scientific field that faces many challenges these
days. Among them are extremely high development costs, long development times,
as well as a low number of new drugs that are approved each year. To solve
these problems, new and innovate technologies are needed that make the drug
discovery process of small-molecules more time and cost-efficient, and which
allow to target previously undruggable target classes such as protein-protein
interactions. Structure-based virtual screenings have become a leading
contender in this context. In this review, we give an introduction to the
foundations of structure-based virtual screenings, and survey their progress in
the past few years. We outline key principles, recent success stories, new
methods, available software, and promising future research directions. Virtual
screenings have an enormous potential for the development of new small-molecule
drugs, and are already starting to transform early-stage drug discovery.


  Effective multicellularity requires both cooperation and competition between
constituent cells. Cooperation involves sacrificing individual fitness in favor
of that of the community, but excessive cooperation makes the community
susceptible to senescence and aging. Competition eliminates unfit senescent
cells via natural selection and thus slows down aging, but excessive
competition makes the community susceptible to cheaters, as exemplified by
cancer and cancer-like phenomena. These observations suggest that an optimal
level of intercellular competition in a multicellular organism maximizes
organismal vitality by delaying the inevitability of aging. We quantify this
idea using a statistical mechanical framework that leads to a generalized
replicator dynamical system for the population of cells that change their
vitality and cooperation due to somatic mutations that make them susceptible to
aging and/or cancer. By accounting for the cost of cooperation and strength of
competition in a minimal setting, we show that our model predicts an optimal
value of competition that maximizes vitality and delays the inevitability of
senescence or cancer.


  The classical infinitesimal model is a simple and robust model for the
inheritance of quantitative traits. In this model, a quantitative trait is
expressed as the sum of a genetic and a non-genetic (environmental) component
and the genetic component of offspring traits within a family follows a normal
distribution around the average of the parents' trait values, and has a
variance that is independent of the trait values of the parents. In previous
work, Barton et al.(2017), we showed that when trait values are determined by
the sum of a large number of Mendelian factors, each of small effect, one can
justify the infinitesimal model as limit of Mendelian inheritance.
  In this paper, we show that the robustness of the infinitesimal model extends
to include dominance. We define the model in terms of classical quantities of
quantitative genetics, before justifying it as a limit of Mendelian inheritance
as the number, M, of underlying loci tends to infinity. As in the additive
case, the multivariate normal distribution of trait values across the pedigree
can be expressed in terms of variance components in an ancestral population and
probabilities of identity by descent determined by the pedigree. In this
setting, it is natural to decompose trait values, not just into the additive
and dominance components, but into a component that is shared by all
individuals within the family and an independent `residual' for each offspring,
which captures the randomness of Mendelian inheritance. We show that, even if
we condition on parental trait values, both the shared component and the
residuals within each family will be asymptotically normally distributed as the
number of loci tends to infinity, with an error of order 1/\sqrt{M}.
  We illustrate our results with some numerical examples.


  Latent variable models such as the Variational Auto-Encoder (VAE) have become
a go-to tool for analyzing biological data, especially in the field of
single-cell genomics. One remaining challenge is the interpretability of latent
variables as biological processes that define a cell's identity. Outside of
biological applications, this problem is commonly referred to as learning
disentangled representations. Although several disentanglement-promoting
variants of the VAE were introduced, and applied to single-cell genomics data,
this task has been shown to be infeasible from independent and identically
distributed measurements, without additional structure. Instead, recent methods
propose to leverage non-stationary data, as well as the sparse mechanism shift
assumption in order to learn disentangled representations with a causal
semantic. Here, we extend the application of these methodological advances to
the analysis of single-cell genomics data with genetic or chemical
perturbations. More precisely, we propose a deep generative model of
single-cell gene expression data for which each perturbation is treated as a
stochastic intervention targeting an unknown, but sparse, subset of latent
variables. We benchmark these methods on simulated single-cell data to evaluate
their performance at latent units recovery, causal target identification and
out-of-domain generalization. Finally, we apply those approaches to two
real-world large-scale gene perturbation data sets and find that models that
exploit the sparse mechanism shift hypothesis surpass contemporary methods on a
transfer learning task. We implement our new model and benchmarks using the
scvi-tools library, and release it as open-source software at
https://github.com/Genentech/sVAE.


  We consider the counting problem of the number of \textit{leaf-labeled
increasing trees}, where internal nodes may have an arbitrary number of
descendants. The set of all such trees is a discrete representation of the
genealogies obtained under certain population-genetical models such as
multiple-merger coalescents. While the combinatorics of the binary trees among
those are well understood, for the number of all trees only an approximate
asymptotic formula is known. In this work, we validate this formula up to
constant terms and compare the asymptotic behavior of the number of all
leaf-labelled increasing trees to that of binary, ternary and quaternary trees.


  Although rare, an increasing number of proteins have been observed to contain
entanglements in their native structures. To gain more insight into the
significance of protein knotting, researchers have been investigating protein
knot formation using both experimental and theoretical methods. Motivated by
the hypothesized folding pathway of $\alpha$-haloacid dehalogenase (DehI)
protein, Flapan, He, and Wong proposed a theory of how protein knots form,
which includes existing folding pathways described by Taylor and B\"olinger et
al. as special cases. In their topological descriptions, two loops in an
unknotted open protein chain containing at most two twists each come close
together, and one end of the protein eventually passes through the two loops.
In this paper, we build on Flapan, He, and Wong's theory where we pay attention
to the crossing signs of the threading process and assume that the unknotted
protein chain may arrange itself into a more complicated configuration before
threading occurs. We then apply tangle calculus, originally developed by Ernst
and Sumners to analyze the action of specific proteins on DNA, to give all
possible knots or knotoids that may be discovered in the future according to
our model and give recipes for engineering specific knots in proteins from
simpler pieces. We show why twists knots are the most likely knots to occur in
proteins. We use chirality to show that the most likely knots to occur in
proteins via Taylor's twisted hairpin model are the knots $+3_1$, $4_1$, and
$-5_2$.


  Single-cell reference atlases are large-scale, cell-level maps that capture
cellular heterogeneity within an organ using single cell genomics. Given their
size and cellular diversity, these atlases serve as high-quality training data
for the transfer of cell type labels to new datasets. Such label transfer,
however, must be robust to domain shifts in gene expression due to measurement
technique, lab specifics and more general batch effects. This requires methods
that provide uncertainty estimates on the cell type predictions to ensure
correct interpretation. Here, for the first time, we introduce uncertainty
quantification methods for cell type classification on single-cell reference
atlases. We benchmark four model classes and show that currently used models
lack calibration, robustness, and actionable uncertainty scores. Furthermore,
we demonstrate how models that quantify uncertainty are better suited to detect
unseen cell types in the setting of atlas-level cell type transfer.


  Collective behaviors exhibited by animal groups, such as fish schools, bird
flocks, or insect swarms are fascinating examples of self-organization in
biology. Concepts and methods from statistical physics have been used to argue
theoretically about the potential consequences of collective effects in such
living systems. In particular, it has been proposed that such collective
systems should operate close to a phase transition, specifically a
(pseudo-)critical point, in order to optimize their capability for collective
computation. In this chapter, we will first review relevant phase transitions
exhibited by animal collectives, pointing out the difficulties of applying
concepts from statistical physics to biological systems. Then we will discuss
the current state of research on the "criticality hypothesis", including
methods for how to measure distance from criticality and specific functional
consequences for animal groups operating near a phase transition. We will
highlight the emerging view that de-emphasizes the optimality of being exactly
at a critical point and instead explores the potential benefits of living
systems being able to tune to an optimal distance from criticality. We will
close by laying out future challenges for studying collective behavior at the
interface of physics and biology.


  The tumor microenvironment (TME) significantly impacts cancer prognosis due
to its immune composition. While therapies for altering the immune composition,
including immunotherapies, have shown exciting results for treating
hematological cancers, they are less effective for immunologically-cold, solid
tumors. Spatial omics technologies capture the spatial organization of the TME
with unprecedented molecular detail, revealing the relationship between immune
cell localization and molecular signals. Here, we formulate T-cell infiltration
prediction as a self-supervised machine learning problem and develop a
counterfactual optimization strategy that leverages large scale spatial omics
profiles of patient tumors to design tumor perturbations predicted to boost
T-cell infiltration. A convolutional neural network predicts T-cell
distribution based on signaling molecules in the TME provided by imaging mass
cytometry. Gradient-based counterfactual generation, then, computes
perturbations predicted to boost T-cell abundance. We apply our framework to
melanoma, colorectal cancer liver metastases, and breast tumor data,
discovering combinatorial perturbations predicted to support T-cell
infiltration across tens to hundreds of patients. This work presents a paradigm
for counterfactual-based prediction and design of cancer therapeutics using
spatial omics data.


  In our previous paper, the logistic curve of the removed number was derived
from SIR and SEIR models in the case of the small basic reproduction number. In
this paper, we derive various logistic curves of the removed, unsusceptible and
infectious numbers respectively from SIS and SIR models in the case of small
and large basic reproduction numbers.


  In this work we study the dynamics in an apocalypse where the individuals can
temporarily become zombies, returning to the living state again or die. We
describe this dynamics using a modified version of the epidemic SIS model. The
zombies can die when in contact with a terminator, which is part of the living
people. To define the possible interactions we use an Erd\"{o}s-R\`{e}nyi
network and we calculate how the absorbing and active phase of the original SIS
model is influenced by the average network degree. We also use the component
and domain distribution of the network to understand how this influence occurs.


  In this work we propose a measure-valued stochastic process representing the
dynamics of a virus population, structured by phenotypic traits and
geographical space, and where viruses are transported between spatial locations
by mechanical vectors. As a first example of the use of this model, we show how
to use this model to infer results on the probability of extinction of the
virus population. Later, by combining various scalings on population sizes,
speed of diffusion of vectors, and other relevant model parameters, we show the
emergence of two systems of integro-differential equations as Macroscopic
descriptions of the system. Under the existence of densities at time zero, we
also show the propagation of this property for later times, and derive the
strong formulation of the limiting systems of IDEs. These strong formulations,
in a sense, correspond to spatial Lotka-Volterra competition models with
mutation and vector-borne dispersal.


  Catch-and-effort data are among the primary sources of information for
assessing the status of terrestrial wildlife and fish. In fishery science,
elaborate stock-assessment models are fitted to such data in order to estimate
fish-population sizes and guide management decisions. Given the importance of
catch-and-effort data, we scoured a comprehensive dataset pertaining to
albacore tuna (Thunnus alalunga) in the north Pacific ocean for novel
ecological information content about this commercially valuable species.
Specifically, we used unsupervised learning based on finite mixture modelling
to reveal that the north Pacific albacore-tuna stock can be divided into four
pseudo-cohorts ranging in age from approximately 3 to 12 years old. We
discovered that smaller size pseudo-cohorts inhabit relatively high --
subtropical to temperate -- latitudes, with hotspots off the coast of Japan.
Larger size pseudo-cohorts inhabit lower -- tropical to subtropical --
latitudes, with hotspots in the western and central north Pacific. These
results offer evidence that albacore tuna prefer different habitats depending
on their size and age, and point to long-term migratory routes for the species
that the current tagging technology is unlikely to capture in full. We discuss
the implications of the results for data-driven modelling of albacore tuna in
the north Pacific, as well as the management of the north Pacific albacore-tuna
fishery.


  Tuberculosis is one of the most lethal contagious diseases caused by
Mycobacterium tuberculosis (MTB), in many cases, the infected did not show any
symptoms, because the bacilli entered the dormant stage in granulomas. The
dormant stage of MTB is also associated with higher resistance to drugs and the
immune system. Among multiple epigenetic regulations critical to MTB stress
responses, DNA methylation is necessary for the survival of MTB in hypoxic
conditions, which is a common stress event during granuloma formation. This
review gathers previous findings and demonstrates a meta-analysis by collecting
hypoxia gene expression data from several articles and perform association
analysis between those genes and methylation site profiles across whole genomes
of representative strains pf lineage 2 and 4. While more data is required for
more conclusive support, our results suggest that methylation sites in the
possible promoter regions may induce differential gene regulation in hypoxia.


  Polyploidization is an important evolutionary process which affects organisms
ranging from plants to fish and fungi. The signal left behind by it is in the
form of a species' ploidy level (number of complete chromosome sets found in a
cell) which is inherently non-treelike. Currently available tools for
reconstructing the evolutionary past of a polyploid dataset generally start
with a multi-labelled tree obtained for a dataset of interest and then derive a
(phylogenetic) network from that tree in some way that reflects that past by
interpreting the networks's vertices of indegree at least two as
polyploidization events. Since obtaining such a tree can be computationally
expensive it is paramount to have alternative approaches available that allow
one to shed light into the reticulate evolutionary past of a polyploid dataset.
SPRINT aims to reconstruct the evolutionary past of a polyploid dataset in
terms of a binary network which realises the dataset's ploidy profile (vector
of ploidy levels of the dataset's taxa) and requires the fewest number of
polyploidization events. It does this by representing the ploidy level of a
species x in terms of the number of directed paths from the root of the network
to the leaf of the network labelled by x. SPRINT is distributed on GitHub:
https://github.com/lmaher1/SPRINT.


  This work presents a population genetic model of evolution, which includes
haploid selection, mutation, recombination, and drift. The mutation-selection
equilibrium can be expressed exactly in closed form for arbitrary fitness
functions without resorting to diffusion approximations. Tractability is
achieved by generating new offspring using n-parent rather than 2-parent
recombination. While this enforces linkage equilibrium among offspring, it
allows analysis of the whole population under linkage disequilibrium. We derive
a general and exact relationship between fitness fluctuations and response to
selection. Our assumptions allow analytical calculation of the stationary
distribution of the model for a variety of non-trivial fitness functions. These
results allow us to speak to genetic architecture, i.e., what stationary
distributions result from different fitness functions. This paper presents
methods for exactly deriving stationary states for finite and infinite
populations. This method can be applied to many fitness functions, and we give
exact calculations for four of these. These results allow us to investigate
metastability, tradeoffs between fitness functions, and even consider
error-correcting codes


  Probabilistic programming frameworks are powerful tools for statistical
modelling and inference. They are not immediately generalisable to phylogenetic
problems due to the particular computational properties of the phylogenetic
tree object. TreeFlow is a software library for probabilistic programming and
automatic differentiation with phylogenetic trees. It implements inference
algorithms for phylogenetic tree times and model parameters given a tree
topology. We demonstrate how TreeFlow can be used to quickly implement and
assess new models. We also show that it provides reasonable performance for
gradient-based inference algorithms compared to specialized computational
libraries for phylogenetics.


  Human mobility contributes to the fast spatio-temporal propagation of
infectious diseases. During an outbreak, monitoring the infection situation on
either side of an international border is very crucial as there is always a
higher risk of disease importation associated with cross-border migration.
Mechanistic models are effective tools to investigate the consequences of
cross-border mobility on disease dynamics and help in designing effective
control strategies. However, in practice, due to the unavailability of
cross-border mobility data, it becomes difficult to propose reliable,
model-based strategies. In this study, we propose a method for estimating
cross-border mobility flux between any pair of regions that share an
international border from the observed difference in the timing of the
infection peak in each region. Assuming the underlying disease dynamics is
governed by a Susceptible-Infected-Recovered (SIR) model, we employ stochastic
simulations to obtain the maximum likelihood cross-border mobility estimate for
any pair of regions where the difference in peak time can be measured. We then
investigate how the estimate of cross-border mobility flux varies depending on
the disease transmission rate, which is a key epidemiological parameter. We
further show that the uncertainty in mobility flux estimates decreases for
higher disease transmission rates and larger observed differences in peak
timing. Finally, as a case study, we apply the method to some selected regions
along the Poland-Germany border which are directly connected through multiple
modes of transportation and quantify the cross-border fluxes from the COVID-19
cases data during the period $20^{\rm th}$ February $2021$ to $20^{\rm th}$
June $2021$.


  A number of studies have shown that cumulative mortality followed a Gompertz
curve in the initial Covid pandemic period, March-April 2020. We show that the
Gompertz curve is incompatible with expected initial logistic growth curves as
predicted by traditional Susceptible-Infected-Recovered (SIR) models, and
propose a new theory which better explains the nature of the mortality
characteristics based on a global biosphere disturbance. Second, we show that
for the Gompertz curve to emerge, the disturbance has to act on everyone
simultaneously, rejecting the possibility of a disease propagation stage.
Third, we connect logistic growth with Gompertzian growth by augmenting the
logistic growth equation with higher order interaction terms, and show that the
SIR model family is compatible with Gompertzian growth only when all nodes in
the transmission network communicate with infinite speed and interaction.
Crucially, this augmentation must be accompanied by a causality-reversal where
the source of growth is not the pool of infected but the pool of susceptible
people. We thus find a novel bridge between logistic and Gompertzian growth,
separate from the existing Richards model (also called $\theta$-logistic
growth).


  The dynamics of a chemical reaction network (CRN) is often modelled under the
assumption of mass action kinetics by a system of ordinary differential
equations (ODEs) with polynomial right-hand sides that describe the time
evolution of concentrations of chemical species involved. Given an arbitrarily
large integer $K \in {\mathbb N}$, we show that there exists a CRN such that
its ODE model has at least $K$ stable limit cycles. Such a CRN can be
constructed with reactions of at most second order provided that the number of
chemical species grows linearly with $K$. Bounds on the minimal number of
chemical species and the minimal number of chemical reactions are presented for
CRNs with $K$ stable limit cycles and at most second order or seventh order
kinetics. We also show that CRNs with only two chemical species can have $K$
stable limit cycles, when the order of chemical reactions grows linearly with
$K$.


  We introduce a mathematical model of savanna vegetation dynamics. The usual
approach of nonequilibrium ecology is extended by including the impact of wet
and dry seasons. We present and rigorously analyze a model describing a mixed
woodland-grassland ecosystem with stochastic environmental noise in the form of
vegetation biomass losses manifesting fires. Both, the probability of ignition
and the strength of these losses depend on the current season (as well as
vegetation growth rates etc.). Formally it requires an introduction and
analysis of a system that is a piecewise deterministic Markov process with
parameters switching between given constant periods of time. We study the long
time behavior of time averages for such processes.


  Probabilistic Boolean Networks have been proposed for estimating the
behaviour of dynamical systems as they combine rule-based modelling with
uncertainty principles. Inferring PBNs directly from gene data is challenging
however, especially when data is costly to collect and/or noisy, e.g., in the
case of gene expression profile data. In this paper, we present a reproducible
method for inferring PBNs directly from real gene expression data measurements
taken when the system was at a steady state. The steady-state dynamics of PBNs
is of special interest in the analysis of biological machinery. The proposed
approach does not rely on reconstructing the state evolution of the network,
which is computationally intractable for larger networks. We demonstrate the
method on samples of real gene expression profiling data from a well-known
study on metastatic melanoma. The pipeline is implemented using Python and we
make it publicly available.


  The education sector has suffered a catastrophic setback due to ongoing
COVID-pandemic, with classrooms being closed indefinitely. The current study
aims to solve the existing dilemma by examining COVID transmission inside a
classroom and providing long-term sustainable solutions. In this work, a
standard 5m x 3m x 5m classroom is considered where 24 students are seated,
accompanied by a teacher. A computational fluid dynamics simulation based on
OpenFOAM is performed using a Eulerian-Lagrangian framework. Based on the
stochastic dose response framework, we have evaluated the infection risk in the
classroom for two distinct cases: (i) certain students are infected (ii) the
teacher is infected. If the teacher is infected, the probability of infection
could reach 100% for certain students. When certain students are infected, the
maximum infection risk for a susceptible person reaches 30%. The commonly used
cloth mask proves to be ineffective in providing protection against infection
transmission reducing the maximum infection probability by approximately 26%
only. Another commonly used solution in the form of shields installed on desks
have also failed to provide adequate protection against infection reducing the
infection risk only by 50%. Furthermore, the shields serves as a source of
fomite mode of infection. Screens suspended from the ceiling, which entrap
droplets, have been proposed as a novel solution that reduces the infection
risk by 90% and 95% compared to the no screen scenario besides being completely
devoid of fomite infection mode. As a result of the screens, the class-time can
be extended by 55 minutes.


  We study a parametric version of the Fermat-Weber problem with respect to an
asymmetric distance function, which occurs naturally in tropical geometry. Our
results yield a method for constructing phylogenetic supertrees.


  Mathematical models have been used to understand the spread patterns of
infectious diseases such as Coronavirus Disease 2019 (COVID-19). The
transmission component of the models can be modelled in an age-dependent manner
via introducing contact matrix for the population, which describes the contact
rates between the age groups. Since social contact patterns vary from country
to country, we can compare and group the countries using the corresponding
contact matrices. In this paper, we present a framework for clustering
countries based on their contact matrices with respect to an underlying
epidemic model. Since the pipeline is generic and modular, we demonstrate its
application in a COVID-19 model from R\"ost et. al. which gives a hint about
which countries can be compared in a pandemic situation, when only
non-pharmaceutical interventions are available.


  Representation learning for proteins has primarily focused on the global
understanding of protein sequences regardless of their length. However, shorter
proteins (known as peptides) take on distinct structures and functions compared
to their longer counterparts. Unfortunately, there are not as many naturally
occurring peptides available to be sequenced and therefore less
peptide-specific data to train with. In this paper, we propose a new peptide
data augmentation scheme, where we train peptide language models on
artificially constructed peptides that are small contiguous subsets of longer,
wild-type proteins; we refer to the training peptides as "chopped proteins". We
evaluate the representation potential of models trained with chopped proteins
versus natural peptides and find that training language models with chopped
proteins results in more generalized embeddings for short protein sequences.
These peptide-specific models also retain information about the original
protein they were derived from better than language models trained on
full-length proteins. We compare masked language model training objectives to
three novel peptide-specific training objectives: next-peptide prediction,
contrastive peptide selection and evolution-weighted MLM. We demonstrate
improved zero-shot learning performance for a deep mutational scan peptides
benchmark.


  Human leukocyte antigen (HLA) genes are associated with a variety of
diseases, however direct typing of HLA is time and cost consuming. Thus various
imputation methods using sequential SNPs data have been proposed based on
statistical or deep learning models, e.g. CNN-based model, named DEEP*HLA.
However, imputation efficiency is not sufficient for in frequent alleles and a
large size of reference panel is required. Here, we developed a
Transformer-based model to impute HLA alleles, named "HLA Reliable IMputatioN
by Transformer (HLARIMNT)" to take advantage of sequential nature of SNPs data.
We validated the performance of HLARIMNT using two different reference panels;
Pan-Asian reference panel (n = 530) and Type 1 Diabetes Genetics Consortium
(T1DGC) reference panel (n = 5,225), as well as the mixture of those two panels
(n = 1,060). HLARIMNT achieved higher accuracy than DEEP*HLA by several
indices, especially for infrequent alleles. We also varied the size of data
used for training, and HLARIMNT imputed more accurately among any size of
training data. These results suggest that Transformer-based model may impute
efficiently not only HLA types but also any other gene types from sequential
SNPs data.


  In the context of reconstructing phylogenetic networks from a collection of
phylogenetic trees, several characterisations and subsequently algorithms have
been established to reconstruct a phylogenetic network that collectively embeds
all trees in the input in some minimum way. For many instances however, the
resulting network also embeds additional phylogenetic trees that are not part
of the input. However, little is known about these inferred trees. In this
paper, we explore the relationships among all phylogenetic trees that are
embedded in a given phylogenetic network. First, we investigate some
combinatorial properties of the collection P of all rooted binary phylogenetic
trees that are embedded in a rooted binary phylogenetic network N. To this end,
we associated a particular graph G, which we call rSPR graph, with the elements
in P and show that, if |P|=2^k, where k is the number of vertices with
in-degree two in N, then G has a Hamiltonian cycle. Second, by exploiting rSPR
graphs and properties of hypercubes, we turn to the well-studied class of
rooted binary level-1 networks and give necessary and sufficient conditions for
when a set of rooted binary phylogenetic trees can be embedded in a level-1
network without inferring any additional trees. Lastly, we show how these
conditions translate into a polynomial-time algorithm to reconstruct such a
network if it exists.


  Rare genetic disorders affect more than 6% of the global population. Reaching
a diagnosis is challenging because rare disorders are very diverse. Many
disorders have recognizable facial features that are hints for clinicians to
diagnose patients. Previous work, such as GestaltMatcher, utilized
representation vectors produced by a DCNN similar to AlexNet to match patients
in high-dimensional feature space to support "unseen" ultra-rare disorders.
However, the architecture and dataset used for transfer learning in
GestaltMatcher have become outdated. Moreover, a way to train the model for
generating better representation vectors for unseen ultra-rare disorders has
not yet been studied. Because of the overall scarcity of patients with
ultra-rare disorders, it is infeasible to directly train a model on them.
Therefore, we first analyzed the influence of replacing GestaltMatcher DCNN
with a state-of-the-art face recognition approach, iResNet with ArcFace.
Additionally, we experimented with different face recognition datasets for
transfer learning. Furthermore, we proposed test-time augmentation, and model
ensembles that mix general face verification models and models specific for
verifying disorders to improve the disorder verification accuracy of unseen
ultra-rare disorders. Our proposed ensemble model achieves state-of-the-art
performance on both seen and unseen disorders.


  The first step in drug discovery is finding drug molecule moieties with
medicinal activity against specific targets. Therefore, it is crucial to
investigate the interaction between drug-target proteins and small chemical
molecules. However, traditional experimental methods for discovering potential
small drug molecules are labor-intensive and time-consuming. There is currently
a lot of interest in building computational models to screen small drug
molecules using drug molecule-related databases. In this paper, we propose a
method for predicting drug-target binding affinity using deep learning models.
This method uses a modified GRU and GNN to extract features from the
drug-target protein sequences and the drug molecule map, respectively, to
obtain their feature vectors. The combined vectors are used as vector
representations of drug-target molecule pairs and then fed into a fully
connected network to predict drug-target binding affinity. This proposed model
demonstrates its accuracy and effectiveness in predicting drug-target binding
affinity on the DAVIS and KIBA datasets.


  Consider a graph whose vertices are populated by identical objects, together
with an algorithm for the time-evolution of the number of objects placed at
each of the vertices. The discrete dynamics of these objects can be observed
and studied using simple and inexpensive laboratory settings. There are many
similarities but also many differences between such population dynamics and the
quantum dynamics of a particle hopping on the same graph. In this work, we show
that a specific decoration of the original graph enables an exact mapping
between the models of population and quantum dynamics. As such, population
dynamics over graphs is yet another classical platform that can simulate
quantum effects. Several examples are used to demonstrate this claim.


  Characterizing the differential excision of mRNA is critical for
understanding the functional complexity of a cell or tissue, from normal
developmental processes to disease pathogenesis. Most transcript reconstruction
methods infer full-length transcripts from high-throughput sequencing data.
However, this is a challenging task due to incomplete annotations and the
differential expression of transcripts across cell-types, tissues, and
experimental conditions. Several recent methods circumvent these difficulties
by considering local splicing events, but these methods lose transcript-level
splicing information and may conflate transcripts. We develop the first
probabilistic model that reconciles the transcript and local splicing
perspectives. First, we formalize the sequence of mRNA excisions (SME)
reconstruction problem, which aims to assemble variable-length sequences of
mRNA excisions from RNA-sequencing data. We then present a novel hierarchical
Bayesian admixture model for the Reconstruction of Excised mRNA (BREM). BREM
interpolates between local splicing events and full-length transcripts and thus
focuses only on SMEs that have high posterior probability. We develop posterior
inference algorithms based on Gibbs sampling and local search of independent
sets and characterize differential SME usage using generalized linear models
based on converged BREM model parameters. We show that BREM achieves higher F1
score for reconstruction tasks and improved accuracy and sensitivity in
differential splicing when compared with four state-of-the-art transcript and
local splicing methods on simulated data. Lastly, we evaluate BREM on both bulk
and scRNA sequencing data based on transcript reconstruction, novelty of
transcripts produced, model sensitivity to hyperparameters, and a functional
analysis of differentially expressed SMEs, demonstrating that BREM captures
relevant biological signal.


  DNA computing, a nontraditional computing mechanism, provides a feasible and
effective method for solving NP-hard problems because of the vast parallelism
and high-density storage of DNA molecules. Although DNA computing has been
exploited to solve various intractable computational problems, such as the
Hamiltonian path problem, SAT problem, and graph coloring problem, there has
been little discussion of designing universal DNA computing-based models, which
can solve a class of problems. In this paper, by leveraging the dynamic and
enzyme-free properties of DNA strand displacement, we propose a universal model
named DCMSubset for solving subset problems in graph theory. The model aims to
find a minimum (or maximum) set satisfying given constraints. For each element
x involved in a given problem, DCMSubset uses an exclusive single-stranded DNA
molecule to model x as well as a specific DNA complex to model the relationship
between x and other elements. Based on the proposed model, we conducted
simulation and biochemical experiments on three kinds of subset problems, a
minimum dominating set, maximum independent set, and minimum vertex cover. We
observed that DCMSubset can also be used to solve the graph coloring problem.
Moreover, we extended DCMSubset to a model for solving the SAT problem. The
results of experiments showed the feasibility and university of the proposed
method. Our results highlighted the potential for DNA strand displacement to
act as a computation tool to solve NP-hard problems.


  DNA methylation is an important epigenetic mark that has been studied
extensively for its regulatory role in biological processes and diseases. WGBS
allows for genome-wide measurements of DNA methylation up to single-base
resolutions, yet poses challenges in identifying significantly different
methylation patterns across distinct biological conditions. We propose a novel
methylome change-point model which describes the joint dynamics of methylation
regimes of a case and a control group and benefits from taking into account the
information of neighbouring methylation sites among all available samples. We
also devise particle filtering and smoothing algorithms to perform efficient
inference of the latent methylation patterns. We illustrate that our approach
can detect and test for very flexible differential methylation signatures with
high power while controlling Type-I error measures.


  Third-generation nanopore sequencers offer a feature called selective
sequencing or 'Read Until' that allows genomic reads to be analyzed in
real-time and abandoned halfway, if not belonging to a genomic region of
'interest'. This selective sequencing opens the door to important applications
such as rapid and low-cost genetic tests. The latency in analyzing should be as
low as possible for selective sequencing to be effective so that unnecessary
reads can be rejected as early as possible. However, existing methods that
employ subsequence Dynamic Time Warping (sDTW) algorithm for this problem are
too computationally intensive that a massive workstation with dozens of CPU
cores still struggles to keep up with the data rate of a mobile phone-sized
MinION sequencer. In this paper, we present Hardware Accelerated Read Until
(HARU), a resource-efficient hardware-software co-design-based method that
exploits a low-cost and portable heterogeneous MPSoC platform with on-chip FPGA
to accelerate the sDTW-based Read Until algorithm. Experimental results show
that HARU on a Xilinx FPGA embedded with a 4-core ARM processor is around 2.5X
faster than a highly optimized multi-threaded software version (around 85X
faster than the existing unoptimized multi-threaded software) running on a
sophisticated server with 36-core Intel Xeon processor for a SARS-CoV-2
dataset. The energy consumption of HARU is two orders of magnitudes lower than
the same application executing on the 36-core server. Source code for HARU sDTW
module is available as open-source at https://github.com/beebdev/HARU and an
example application that utilises HARU is at
https://github.com/beebdev/sigfish-haru.


  In this paper, we tested several sparse optimization algorithms based on the
public dataset of the DREAM5 Gene Regulatory Network Inference Challenge. And
we find that introducing 20% of the regulatory network as a priori known data
can provide a basis for parameter selection of inference algorithms, thus
improving prediction efficiency and accuracy. In addition to testing common
sparse optimization methods, we also developed voting algorithms by bagging
them. Experiments on the DREAM5 dataset show that the sparse optimization-based
inference of the moderation relation works well, achieving better results than
the official DREAM5 results on three datasets. However, the performance of
traditional independent algorithms varies greatly in the face of different
datasets, while our voting algorithm achieves the best results on three of the
four datasets.


  The COVID-19 pandemic has evolved over time through multiple spatial and
temporal dynamics. The varying extent of interactions among different
geographical areas can result to a complex pattern of spreading so that
influences between these areas can be hard to discern. Here, we use
cross-correlation analysis to detect synchronous evolution and potential
inter-influences in the time evolution of new COVID-19 cases at the county
level in the USA. Our analysis identified two main time periods with
distinguishable features in the behavior of correlations. In the first phase,
there were few strong correlations which only emerged between urban areas. In
the second phase of the epidemic, strong correlations became widespread and
there was a clear directionality of influence from urban to rural areas. In
general, the effect of distance between two counties was much weaker than that
of the counties population. Such analysis can provide possible clues on the
evolution of the disease and may identify parts of the country where
intervention may be more efficient in limiting the disease spread.


  New Q-conditional (nonclassical) symmetries and exact solutions of the
hunter-gatherer-farmer population model proposed by Aoki, Shida and Shigesada
(Theor. Popul. Biol. 1996;50:1-17) are constructed. The main method used for
the aforementioned purposes is an extension of the nonclassical method for
system of partial differential equations. An analysis of properties of the
exact solutions obtained and their biological interpretation are carried out.
New results are compared with those derived in recent studies devoted to the
same model.


  Background: There is a lack of population-based studies on the effectiveness
and durability of the SARS-CoV-2-induced immune protection during the Omicron
wave. Methods: This retrospective study included 470 159 cases aged 60 years or
older, who tested positive for SARS-CoV-2 between January 1 and June 30, 2022
in Bavaria, Germany. We examined time to death, measured from the earliest date
of reporting/symptom onset until day 60. Cox models were used to estimate
adjusted hazard ratios (HR) for sex, age, earliest calendar date of documented
infection, and level of vaccine-induced or hybrid immunity (LI). Results: We
observed 3836 deaths (case fatality rate 0.82%). Risk of death was
significantly lower in all LIs than in unvaccinated cases (adjusted HR for a
full primary LI achieved within the last six months: 0.30, 95% confidence
interval (CI) 0.23-0.39; after six months: 0.46, 95% CI 0.35-0.60). Boosting
further decreased the risk of death (within the last three months: HR 0.17, 95%
CI 0.15-0.20; after three months: HR 0.25, 95% CI 0.21-0.29). Conclusion: In
elderly Bavarians, an increasing LI was associated with an increasing
protection against death during the Omicron wave. Protection, however, may
decrease to some extent after several months.


  Long non-coding RNAs (lncRNAs) are important regulators to modulate gene
expression and cell proliferation in the developing human brain. Previous
methods mainly use bulk lncRNA and mRNA expression data to study lncRNA
regulation. However, to analyze lncRNA regulation regarding individual cells,
we focus on single-cell RNA-sequencing (scRNA-seq) data instead of bulk data.
Recent advance in scRNA-seq has provided a way to investigate lncRNA regulation
at single-cell level. We will propose a novel computational method, CSlncR
(cell-specific lncRNA regulation), which combines putative lncRNA-mRNA binding
information with scRNA-seq data including lncRNAs and mRNAs to identify
cell-specific lncRNA-mRNA regulation networks at individual cells. To
understand lncRNA regulation at different development stages, we apply CSlncR
to the scRNA-seq data of human neocortex. Network analysis shows that the
lncRNA regulation is unique in each cell from the different human neocortex
development stages. The comparison results indicate that CSlncR is also an
effective tool for predicting cell-specific lncRNA targets and clustering
single cells, which helps understand cell-cell communication.


  Searching for similar genomic sequences is an essential and fundamental step
in biomedical research and an overwhelming majority of genomic analyses.
State-of-the-art computational methods performing such comparisons fail to cope
with the exponential growth of genomic sequencing data. We introduce the
concept of sparsified genomics where we systematically exclude a large number
of bases from genomic sequences and enable much faster and more
memory-efficient processing of the sparsified, shorter genomic sequences, while
providing similar or even higher accuracy compared to processing non-sparsified
sequences. Sparsified genomics provides significant benefits to many genomic
analyses and has broad applicability. We show that sparsifying genomic
sequences greatly accelerates the state-of-the-art read mapper (minimap2) by
2.57-5.38x, 1.13-2.78x, and 3.52-6.28x using real Illumina, HiFi, and ONT
reads, respectively, while providing up to 2.1x smaller memory footprint, 2x
smaller index size, and more truly detected small and structural variations
compared to minimap2. Sparsifying genomic sequences makes containment search
through very large genomes and very large databases 72.7-75.88x faster and
723.3x more storage-efficient than searching through non-sparsified genomic
sequences (with CMash and KMC3). Sparsifying genomic sequences enables robust
microbiome discovery by providing 54.15-61.88x faster and 720x more
storage-efficient taxonomic profiling of metagenomic samples over the
state-of-art tool (Metalign). We design and open-source a framework called
Genome-on-Diet as an example tool for sparsified genomics, which can be freely
downloaded from https://github.com/CMU-SAFARI/Genome-on-Diet.


  In this work, we discuss the SIR epidemiological model and different
variations of it applied to the propagation of the COVID-19 pandemia; we employ
the data of the state of Guanajuato and of Mexico. We present some
considerations that can improve the predictions made by those models. We
consider a time-dependent infection rate, which we adjust to the data. Starting
from a linear regime where the populations are much smaller that the country or
state population and the population of susceptible (S) can be approximated in
convenient units to S approximately 1, we make fits of the parameters. We also
consider the case when the susceptible starts departing from 1, for this case
we adjust an effective contagion rate. We also explore the ratio of detected
populations and the real ones, obtaining that -for the analyzed case it is of
approximately 10%. We estimate the number of deaths by making a fit versus the
recovered cases, this fit is in first approximation linear, but other powers
can give a good agreement. By predictions to past data, we conclude that
adaptations of the SIR model can be of great use in describing pandemia's
propagation, specially in limited time periods.


  The massive amount of genomic data appearing for SARS-CoV-2 since the
beginning of the COVID-19 pandemic has challenged traditional methods for
studying its dynamics. As a result, new methods such as Pangolin, which can
scale to the millions of samples of SARS-CoV-2 currently available, have
appeared. Such a tool is tailored to take as input assembled, aligned and
curated full-length sequences, such as those found in the GISAID database. As
high-throughput sequencing technologies continue to advance, such assembly,
alignment and curation may become a bottleneck, creating a need for methods
which can process raw sequencing reads directly.
  In this paper, we propose Reads2Vec, an alignment-free embedding approach
that can generate a fixed-length feature vector representation directly from
the raw sequencing reads without requiring assembly. Furthermore, since such an
embedding is a numerical representation, it may be applied to highly optimized
classification and clustering algorithms. Experiments on simulated data show
that our proposed embedding obtains better classification results and better
clustering properties contrary to existing alignment-free baselines. In a study
on real data, we show that alignment-free embeddings have better clustering
properties than the Pangolin tool and that the spike region of the SARS-CoV-2
genome heavily informs the alignment-free clusterings, which is consistent with
current biological knowledge of SARS-CoV-2.


  Predicting the evolution of diseases is challenging, especially when the data
availability is scarce and incomplete. The most popular tools for modelling and
predicting infectious disease epidemics are compartmental models. They stratify
the population into compartments according to health status and model the
dynamics of these compartments using dynamical systems. However, these
predefined systems may not capture the true dynamics of the epidemic due to the
complexity of the disease transmission and human interactions. In order to
overcome this drawback, we propose Sparsity and Delay Embedding based
Forecasting (SPADE4) for predicting epidemics. SPADE4 predicts the future
trajectory of an observable variable without the knowledge of the other
variables or the underlying system. We use random features model with sparse
regression to handle the data scarcity issue and employ Takens' delay embedding
theorem to capture the nature of the underlying system from the observed
variable. We show that our approach outperforms compartmental models when
applied to both simulated and real data.