  Weak magnetic field exposure can affect many biological processes across a
wide range of living organisms. Recently, it has been observed that weak
magnetic fields can modulate reactive oxygen species (ROS) concentration,
affecting regeneration in planaria. These effects show unusual nonlinear
dependence on magnetic field strength, including a sign change. In another
study by the same group, superoxide is identified as the particular ROS being
modulated. We propose a radical pair mechanism based on a flavin-superoxide
radical pair to explain the modulation of superoxide production and its effect
on planarian regeneration. The results of our calculations favor a triplet-born
radical pair. Our yield calculations can reproduce the observed magnetic field
dependence, including the sign change. Moreover, to explain the size of the
effect on ROS concentration, we suggest a simple amplification model inspired
by known biochemical mechanisms and lay out the conditions for such a model to
work. Further, we also make empirical predictions concerning the hypomagnetic
field effects on planarian regeneration.


  Cancer progression and monotonic accumulation models were developed to
discover dependencies in the irreversible acquisition of binary traits from
cross-sectional data. They have been used in computational oncology and
virology to understand dependencies in the accumulation of mutations but also
in widely different problems such as malaria progression.
  These tools have been used to make predictions about future and unobserved
states of the system, identify different routes of feature acquisition in
subsets of the data, and improve patient stratification based on the
evolutionary trajectories.
  But these methods have shortcomings, in particular with respect to key
evolutionary assumptions and interpretations. After an overview of the
available methods, I focus on how and why the methods can fail and why their
inferences might not be about the process we intend. I first examine
interpretations and violations of assumptions when methods are used to infer
within-cell restrictions in the accumulation of events. The discussion uses
fitness landscapes and highlights the difficulties that arise from bulk
sequencing and reciprocal sign epistasis, from the conflation of lines of
descent, path of the maximum, and mutational profiles, and from the ambiguous
use of the concepts of exclusivity and alternative evolutionary paths. Next, I
suggest that users recognize the limitations that arise from their sampling
designs, and I examine how the previous concerns are modified when bulk
sequencing is accounted for. I also highlight how bulk sequencing can create an
opportunity for addressing dependencies that arise from frequency-dependent
selection. This review should help identify major standing issues and research
opportunities, and encourage the use of these methods in other areas with a
possibly better fit between entities and model assumptions.


  Transcriptomic data is a treasure-trove in modern molecular biology, as it
offers a comprehensive viewpoint into the intricate nuances of gene expression
dynamics underlying biological systems. This genetic information must be
utilised to infer biomolecular interaction networks that can provide insights
into the complex regulatory mechanisms underpinning the dynamic cellular
processes. Gene regulatory networks and protein-protein interaction networks
are two major classes of such networks. This chapter thoroughly investigates
the wide range of methodologies used for distilling insightful revelations from
transcriptomic data that include association based methods (based on
correlation among expression vectors), probabilistic models (using Bayesian and
Gaussian models), and interologous methods. We reviewed different approaches
for evaluating the significance of interactions based on the network topology
and biological functions of the interacting molecules, and discuss various
strategies for the identification of functional modules. The chapter concludes
with highlighting network based techniques of prioritising key genes, outlining
the centrality based, diffusion based and subgraph based methods. The chapter
provides a meticulous framework for investigating transcriptomic data to
uncover assembly of complex molecular networks for their adaptable analyses
across a broad spectrum of biological domains.


  Signal propagation in biochemical networks is characterized by the inherent
randomness in gene expression and fluctuations of the environmental components,
commonly known as intrinsic and extrinsic noise, respectively. We present a
theoretical framework for noise propagation in a generic two-step cascade
(S$\rightarrow$X$\rightarrow$Y) in terms of intrinsic and extrinsic noise. We
identify different channels of noise transmission that regulate the individual
as well as the overall noise properties of each component. Our analysis shows
that the intrinsic noise of S alleviates the overall noise and information
transmission capacity along the cascade. On the other hand, the intrinsic noise
of X and Y acts as a bottleneck of information transmission.


  We have shown recently that the notion of poking pairwise interactions along
a chain provides a unifying framework for understanding the formation of both
secondary and the tertiary protein structure based on symmetry and geometry.
$\alpha$-helices and $\beta$-sheets are found to be special geometries that
have systematic poking contacts in a repetitive manner with the contacts being
local along the $\alpha$-helix and non-local along a pair of adjacent strands
within a $\beta$-sheet. Pairwise poking interactions also govern tertiary
structure formation, but they are weaker and there are no special geometrical
constraints as in secondary structure formation. Here we demonstrate that
protein turns, the most prevalent non-repetitive structural element in
proteins, are instances of local (as in $\alpha$-helices) and isolated
(non-repetitive) poking pairwise contacts for which the geometrical constraints
are partially relaxed. This simple and purely geometrical definition of protein
turns (also sometimes known as reverse turns, $\beta$-turns, $\beta$-bends,
hairpin bends, $3_{10}$ bends, kinks, widgets, ...) provides a simple framework
for unifying them. We present the results of a systematic analysis and identify
their structural classes as well as their respective amino acid preferences.


  We present a geometrical analysis of the protrusion statistics of side chains
in more than 4,000 high-resolution protein structures. We find that when the
amino acids are grouped based on chemical similarity, the resulting groups
equally well fit the geometrical attributes in a harmonious manner. Our
analysis underscores the role of minimal frustration in the evolution of
protein sequences that fit within a pre-sculpted menu of putative native state
folds.


  Coastal ecosystems are increasingly experiencing anthropogenic pressures such
as climate heating, CO2 increase, metal and organic pollution, overfishing and
resource extraction. Some resulting stressors are more direct like fisheries,
others more indirect like ocean acidification, yet they jointly affect marine
biota, communities and entire ecosystems. While single-stressor effects have
been widely investigated, the interactive effects of multiple stressors on
ecosystems are less researched. In this study, we review the literature on
multiple stressors and their interactive effects in coastal environments across
organisms. We classify the interactions into three categories: synergistic,
additive, and antagonistic. We found phytoplankton and mollusks to be the most
studied taxonomic groups. The stressor combinations of climate warming, ocean
acidification, eutrophication, and metal pollution are the most critical for
coastal ecosystems as they exacerbate adverse effects on physiological traits
such as growth rate, basal respiration, and size. Phytoplankton appears to be
most sensitive to interactions between metal and nutrient pollution. In
nutrient-enriched environments, the presence of metals considerably affects the
uptake of nutrients, and increases respiration costs and toxin production in
phytoplankton. For mollusks, warming and low pH are the most lethal stressors.
The combined effect of heat stress and ocean acidification leads to decreased
growth rate, shell size, and acid-base regulation capacity in mollusks.
However, for a holistic understanding of how coastal food webs will evolve with
ongoing changes, we suggest more research on ecosystem-level responses. This
can be achieved by combining in-situ observations from controlled environments
(e.g. mesocosm experiments) with modelling approaches.


  Algebraic techniques in phylogenetics have historically been successful at
proving identifiability results and have also led to novel reconstruction
algorithms. In this paper, we study the ideal of phylogenetic invariants of the
Cavender-Farris-Neyman (CFN) model on a phylogenetic network with the goal of
providing a description of the invariants which is useful for network
inference. It was previously shown that to characterize the invariants of any
level-1 network, it suffices to understand all sunlet networks, which are those
consisting of a single cycle with a leaf adjacent to each cycle vertex. We show
that the parameterization of an affine open patch of the CFN sunlet model,
which intersects the probability simplex factors through the space of
skew-symmetric matrices via Pfaffians. We then show that this affine patch is
isomorphic to a determinantal variety and give an explicit Gr{\"o}bner basis
for the associated ideal, which involves only polynomially many coordinates.
Lastly, we show that sunlet networks with at least 6 leaves are identifiable
using only these polynomials and run extensive simulations, which show that
these polynomials can be used to accurately infer the correct network from DNA
sequence data.


  In this note, we use a result by Harris (1957) to show that there always
exists a stationary measure (not necessarily a distribution) on a closed
irreducible component of a stochastic reaction network. This measure might not
be unique. In particular, any weakly reversible stochastic reaction network has
a stationary measure on all closed irreducibe components, irrespective whether
it is compelx balanced or not.


  Over the past 40 years, the discovery and development of therapeutic
antibodies to treat disease has become common practice. However, as therapeutic
antibody constructs are becoming more sophisticated (e.g., multi-specifics),
conventional approaches to optimisation are increasingly inefficient. Machine
learning (ML) promises to open up an in silico route to antibody discovery and
help accelerate the development of drug products using a reduced number of
experiments and hence cost. Over the past few years, we have observed rapid
developments in the field of ML-guided antibody discovery and development
(D&D). However, many of the results are difficult to compare or hard to assess
for utility by other experts in the field due to the high diversity in the
datasets and evaluation techniques and metrics that are across industry and
academia. This limitation of the literature curtails the broad adoption of ML
across the industry and slows down overall progress in the field, highlighting
the need to develop standards and guidelines that may help improve the
reproducibility of ML models across different research groups. To address these
challenges, we set out in this perspective to critically review current
practices, explain common pitfalls, and clearly define a set of method
development and evaluation guidelines that can be applied to different types of
ML-based techniques for therapeutic antibody D&D. Specifically, we address in
an end-to-end analysis, challenges associated with all aspects of the ML
process and recommend a set of best practices for each stage.


  Many protein design applications, such as binder or enzyme design, require
scaffolding a structural motif with high precision. Generative modelling
paradigms based on denoising diffusion processes emerged as a leading candidate
to address this motif scaffolding problem and have shown early experimental
success in some cases. In the diffusion paradigm, motif scaffolding is treated
as a conditional generation task, and several conditional generation protocols
were proposed or imported from the Computer Vision literature. However, most of
these protocols are motivated heuristically, e.g. via analogies to Langevin
dynamics, and lack a unifying framework, obscuring connections between the
different approaches. In this work, we unify conditional training and
conditional sampling procedures under one common framework based on the
mathematically well-understood Doob's h-transform. This new perspective allows
us to draw connections between existing methods and propose a new variation on
existing conditional training protocols. We illustrate the effectiveness of
this new protocol in both, image outpainting and motif scaffolding and find
that it outperforms standard methods.


  In this paper, we consider the inverse problem of determining the
coefficients of interaction terms within some Lotka-Volterra models, with
support from boundary observation of its non-negative solutions. In the
physical background, the solutions to the predator-prey model stand for the
population densities for predator and prey and are non-negative, which is a
critical challenge in our inverse problem study. We mainly focus on the unique
identifiability issue and tackle it with the high-order variation method, a
relatively new technique introduced by the second author and his collaborators.
This method can ensure the positivity of solutions and has broader
applicability in other physical models with non-negativity requirements. Our
study improves this method by choosing a more general solution $(u_0,v_0)$ to
expand around, achieving recovery for all interaction terms. By this means, we
improve on the previous results and apply this to physical models to recover
coefficients concerning compression, prey attack, crowding, carrying capacity,
and many other interaction factors in the system. Finally, we apply our results
to study three specific cases: the hydra-effects model, the Holling-Tanner
model and the classic Lotka-Volterra model.


  Designing reward functions that guide generative molecular design (GMD)
algorithms to desirable areas of chemical space is of critical importance in
AI-driven drug discovery. Traditionally, this has been a manual and error-prone
task; the selection of appropriate computational methods to approximate
biological assays is challenging and the aggregation of computed values into a
single score even more so, leading to potential reliance on trial-and-error
approaches. We propose a novel approach for automated reward configuration that
relies solely on experimental data, mitigating the challenges of manual reward
adjustment on drug discovery projects. Our method achieves this by constructing
a ranking over experimental data based on Pareto dominance over the
multi-objective space, then training a neural network to approximate the reward
function such that rankings determined by the predicted reward correlate with
those determined by the Pareto dominance relation. We validate our method using
two case studies. In the first study we simulate Design-Make-Test-Analyse
(DMTA) cycles by alternating reward function updates and generative runs guided
by that function. We show that the learned function adapts over time to yield
compounds that score highly with respect to evaluation functions taken from the
literature. In the second study we apply our algorithm to historical data from
four real drug discovery projects. We show that our algorithm yields reward
functions that outperform the predictive accuracy of human-defined functions,
achieving an improvement of up to 0.4 in Spearman's correlation against a
ground truth evaluation function that encodes the target drug profile for that
project. Our method provides an efficient data-driven way to configure reward
functions for GMD, and serves as a strong baseline for future research into
transformative approaches for the automation of drug discovery.


  The Poisson-Boltzmann (PB) model governs the electrostatics of solvated
biomolecules, i.e., potential, field, energy, and force. These quantities can
provide useful information about protein properties, functions, and dynamics.
By considering the advantages of current algorithms and computer hardware, we
focus on the parallelization of the treecode-accelerated boundary integral
(TABI) PB solver using the Message Passing Interface (MPI) on CPUs and the
direct-sum boundary integral (DSBI) PB solver using KOKKOS on GPUs. We provide
optimization guidance for users when the DSBI solver on GPU or the TABI solver
with MPI on CPU should be used depending on the size of the problem.
Specifically, when the number of unknowns is smaller than a predetermined
threshold, the GPU-accelerated DSBI solver converges rapidly thus has the
potential to perform PB model-based molecular dynamics or Monte Carlo
simulation. As practical appliations, our parallelized boundary integral PB
solvers are used to solve electrostatics on selected proteins that play
significant roles in the spread, treatment, and prevention of COVID-19 virus
diseases. For each selected protein, the simulation produces the electrostatic
solvation energy as a global measurement and electrostatic surface potential
for local details.


  Phylogeny can be inferred using two sources of data from an organism:
morphological data and molecular data. Historically, phylogenies were usually
inferred using morphological characters, but some morphological features may
not necessarily indicate shared heritage. With the introduction of molecular
phylogenies, the base sequence of genes, or amino acid sequence of proteins can
be compared to find the number of similarities or differences to ascertain
levels of relatedness between species. These two types of phylogenies are to be
taken as a data-driven hypothesis about the evolutionary history of the studied
organisms, and a consensus is drawn from the comparison between the different
phylogenies built from the two sources of data, utilizing different methods.


  Single-cell RNA sequencing (scRNA-seq) has emerged as a transformative
technology, offering unparalleled insights into the intricate landscape of
cellular diversity and gene expression dynamics. The analysis of scRNA-seq data
poses challenges attributed to both sparsity and the extensive number of genes
implicated. An increasing number of computational tools are devised for
analyzing and interpreting scRNA-seq data. We present a multiscale differential
geometry (MDG) strategy to exploit the geometric and biological properties
inherent in scRNA-seq data. We assume that those intrinsic properties of cells
lies on a family of low-dimensional manifolds embedded in the high-dimensional
space of scRNA-seq data. Subsequently, we explore these properties via
multiscale cell-cell interactive manifolds. Our multiscale curvature-based
representation serves as a powerful approach to effectively encapsulate the
complex relationships in the cell-cell network. We showcase the utility of our
novel approach by demonstrating its effectiveness in classifying cell types.
This innovative application of differential geometry in scRNA-seq analysis
opens new avenues for understanding the intricacies of biological networks and
holds great potential for network analysis in other fields.


  Single-cell RNA sequencing (scRNA-seq) technology provides high-throughput
gene expression data to study the cellular heterogeneity and dynamics of
complex organisms. Graph neural networks (GNNs) have been widely used for
automatic cell type classification, which is a fundamental problem to solve in
scRNA-seq analysis. However, existing methods do not sufficiently exploit both
gene-gene and cell-cell relationships, and thus the true potential of GNNs is
not realized. In this work, we propose a bilevel graph representation learning
method, named scBiGNN, to simultaneously mine the relationships at both gene
and cell levels for more accurate single-cell classification. Specifically,
scBiGNN comprises two GNN modules to identify cell types. A gene-level GNN is
established to adaptively learn gene-gene interactions and cell representations
via the self-attention mechanism, and a cell-level GNN builds on the cell-cell
graph that is constructed from the cell representations generated by the
gene-level GNN. To tackle the scalability issue for processing a large number
of cells, scBiGNN adopts an Expectation Maximization (EM) framework in which
the two modules are alternately trained via the E-step and M-step to learn from
each other. Through this interaction, the gene- and cell-level structural
information is integrated to gradually enhance the classification performance
of both GNN modules. Experiments on benchmark datasets demonstrate that our
scBiGNN outperforms a variety of existing methods for cell type classification
from scRNA-seq data.


  Dictionary learning (DL) is commonly used in computational biology to tackle
ubiquitous clustering problems due to its conceptual simplicity and relatively
low computational complexity. However, DL algorithms produce results that lack
interpretability and are not optimized for large-scale graph-structured data.
We propose a novel DL algorithm called online convex network dictionary
learning (online cvxNDL) that can handle extremely large datasets and enables
the interpretation of dictionary elements, which serve as cluster
representatives, through convex combinations of real measurements. Moreover,
the algorithm can be applied to network-structured data via specialized
subnetwork sampling techniques.
  To demonstrate the utility of our approach, we apply cvxNDL on 3D-genome
RNAPII ChIA-Drop data to identify important long-range interaction patterns.
ChIA-Drop probes higher-order interactions, and produces hypergraphs whose
nodes represent genomic fragments. The hyperedges represent observed physical
contacts. Our hypergraph model analysis creates an interpretable dictionary of
long-range interaction patterns that accurately represent global chromatin
physical contact maps. Using dictionary information, one can also associate the
contact maps with RNA transcripts and infer cellular functions.
  Our results offer two key insights. First, we demonstrate that online cvxNDL
retains the accuracy of classical DL methods while simultaneously ensuring
unique interpretability and scalability. Second, we identify distinct
collections of proximal and distal interaction patterns involving chromatin
elements shared by related processes across different chromosomes, as well as
patterns unique to specific chromosomes. To associate the dictionary elements
with biological properties of the corresponding chromatin regions, we employ
Gene Ontology enrichment analysis and perform RNA coexpression studies.


  Deep learning-based drug response prediction (DRP) methods can accelerate the
drug discovery process and reduce R\&D costs. Although the mainstream methods
achieve high accuracy in predicting response regression values, the
regression-aware representations of these methods are fragmented and fail to
capture the continuity of the sample order. This phenomenon leads to models
optimized to sub-optimal solution spaces, reducing generalization ability and
may result in significant wasted costs in the drug discovery phase. In this
paper, we propose \MN, a contrastive learning framework with natural language
supervision for the DRP. The \MN~converts regression labels into text, which is
merged with the captions text of the drug response as a second modality of the
samples compared to the traditional modalities (graph, sequence). In each
batch, two modalities of one sample are considered positive pairs and the other
pairs are considered negative pairs. At the same time, in order to enhance the
continuous representation capability of the numerical text, a common-sense
numerical knowledge graph is introduced. We validated several hundred thousand
samples from the Genomics of Drug Sensitivity in Cancer dataset, observing the
average improvement of the DRP method ranges from 7.8\% to 31.4\% with the
application of our framework. The experiments prove that the \MN~effectively
constrains the samples to a continuous distribution in the representation
space, and achieves impressive prediction performance with only a few epochs of
fine-tuning after pre-training. The code is available at:
\url{https://gitee.com/xiaoyibang/clipdrug.git}.


  Nestedness is a common property of communication, finance, trade, and
ecological networks. In networks with high levels of nestedness, the link
positions of low-degree nodes (those with few links) form nested subsets of the
link positions of high-degree nodes (those with many links), leading to matrix
representations with characteristic upper-triangular or staircase patterns.
Recent theoretical work has connected nestedness to the functionality of
complex systems and has suggested it is a structural by-product of the skewed
degree distributions often seen in empirical data. However, mechanisms for
generating nestedness remain poorly understood, limiting the connections that
can be made between system processes and observed network structures. Here, we
show that a simple probabilistic model based on phenology -- the timing of
co-presences among interaction partners -- can produce nested structures and
correctly predict around two-thirds of interactions in two fish market networks
and around one-third of interactions in 22 plant-pollinator networks. Notably,
the links most readily explained by frequent actor co-presences appear to form
a backbone of nested interactions, with the remaining interactions attributable
to opportunistic interactions or preferences for particular interaction
partners that are not routinely available.


  Neural language models have become powerful tools for learning complex
representations of entities in natural language processing tasks. However,
their interpretability remains a significant challenge, particularly in domains
like computational biology where trust in model predictions is crucial. In this
work, we aim to enhance the interpretability of protein language models,
specifically the state-of-the-art ESM model, by identifying and characterizing
knowledge neurons - components that express understanding of key information.
After fine-tuning the ESM model for the task of enzyme sequence classification,
we compare two knowledge neuron selection methods that preserve a subset of
neurons from the original model. The two methods, activation-based and
integrated gradient-based selection, consistently outperform a random baseline.
In particular, these methods show that there is a high density of knowledge
neurons in the key vector prediction networks of self-attention modules. Given
that key vectors specialize in understanding different features of input
sequences, these knowledge neurons could capture knowledge of different enzyme
sequence motifs. In the future, the types of knowledge captured by each neuron
could be characterized.


  To curb the initial spread of SARS-CoV-2, many countries relied on
nation-wide implementation of non-pharmaceutical intervention measures,
resulting in substantial socio-economic impacts. Potentially, subnational
implementations might have had less of a societal impact, but comparable
epidemiological impact. Here, using the first COVID-19 wave in the Netherlands
as a case in point, we address this issue by developing a high-resolution
analysis framework that uses a demographically-stratified population and a
spatially-explicit, dynamic, individual contact-pattern based epidemiology,
calibrated to hospital admissions data and mobility trends extracted from
mobile phone signals and Google. We demonstrate how a subnational approach
could achieve similar level of epidemiological control in terms of hospital
admissions, while some parts of the country could stay open for a longer
period. Our framework is exportable to other countries and settings, and may be
used to develop policies on subnational approach as a better strategic choice
for controlling future epidemics.


  Social reputations facilitate cooperation: those who help others gain a good
reputation, making them more likely to receive help themselves. But when people
hold private views of one another, this cycle of indirect reciprocity breaks
down, as disagreements lead to the perception of unjustified behavior that
ultimately undermines cooperation. Theoretical studies often assume
population-wide agreement about reputations, invoking rapid gossip as an
endogenous mechanism for reaching consensus. However, the theory of indirect
reciprocity lacks a mechanistic description of how gossip actually generates
consensus. Here we develop a mechanistic model of gossip-based indirect
reciprocity that incorporates two alternative forms of gossip: exchanging
information with randomly selected peers or consulting a single gossip source.
We show that these two forms of gossip are mathematically equivalent under an
appropriate transformation of parameters. We derive an analytical expression
for the minimum amount of gossip required to reach sufficient consensus and
stabilize cooperation. We analyze how the amount of gossip necessary for
cooperation depends on the benefits and costs of cooperation, the assessment
rule (social norm), and errors in reputation assessment, strategy execution,
and gossip transmission. Finally, we show that biased gossip can either
facilitate or hinder cooperation, depending on the direction and magnitude of
the bias. Our results contribute to the growing literature on cooperation
facilitated by communication, and they highlight the need to study strategic
interactions coupled with the spread of social information.


  During the COVID-19 pandemic, a large variance of incidence rates on local
level, e.g. cities and districts, within one country has been observed, while
the same non-pharmaceutical measures have been taken to control the spread of
the virus. This variance in incidence rates triggered the question, if the
spread of incidence rates can be explained only by statistical processes and
the local population statistics or if indeed other factors, e.g. local
information campaigns, have to be considered. Within this paper we study the
expected spread of incidence rates in the German State of Rhineland Palatinate
during the second COVID-19 wave using an agent based simulation and find that
the spread of incidence rates can be solely explained by population statistics
and further statical effects.


  The effect of school closures on the spread of COVID-19 has been discussed
among experts and the general public since those measures have been taken only
a few months after the start of the pandemic in 2020. Within this study, the
JuneGermany framework, is used to quantify the impact of school closures in the
German state Rhineland Palatinate using an agent-based simulation approach. It
was found that the simulations predicts a reduction of the number of
infections, hospitalizations as well as death by a factor of 2.5 compared to
scenarios, where no school closures are enforced, during the second wave
between October 2020 and February 2021


  In the paper we consider a system of differential equations with two delays
describing plankton--fish interaction. We study stability of the equilibrium
point corresponding to the presence of phytoplankton and zooplankton in the
system and the absence of fish. In the case of asymptotic stability, we
indicate conditions for the initial data under which solutions stabilize to
this equilibrium point and establish estimates for solutions characterizing the
stabilization rate at infinity. The results are obtained using
Lyapunov-Krasovskii functional. The obtained theoretical results can be used
for numerical study of behavior of solutions to the considered model.


  Mathematical models of infectious disease transmission typically neglect
within-host dynamics. Yet within-host dynamics - including pathogen
replication, host immune responses, and interactions with microbiota - are
crucial not only for determining the progression of disease at the individual
level, but also for driving within-host evolution and onwards transmission, and
therefore shape dynamics at the population level. Various approaches have been
proposed to model both within- and between-host dynamics, but these typically
require considerable simplifying assumptions to couple processes at contrasting
scales (e.g., the within-host dynamics quickly reach a steady state) or are
computationally intensive. Here we propose a novel, readily adaptable and
broadly applicable method for modelling both within- and between-host processes
which can fully couple dynamics across scales and is both realistic and
computationally efficient. By individually tracking the deterministic
within-host dynamics of infected individuals, and stochastically coupling these
to continuous host state variables at the population-level, we take advantage
of fast numerical methods at both scales while still capturing individual
transient within-host dynamics and stochasticity in transmission between hosts.
Our approach closely agrees with full stochastic individual-based simulations
and is especially useful when the within-host dynamics do not rapidly reach a
steady state or over longer timescales to track pathogen evolution. By applying
our method to different pathogen growth scenarios we show how common
simplifying assumptions fundamentally change epidemiological and evolutionary
dynamics.


  Alzheimer's Disease (AD), a neurodegenerative disorder, is reported as one of
the most severe health and socioeconomic problems in current public health. Tau
proteins are assumed to be a crucial driving factor of AD that detach from
microtubules (MT) and accumulate as neurotoxic aggregates in the brains of AD
patients. Extensive experimental and computational research has observed that
phosphorylation at specific tau residues enhances aggregation, but the exact
mechanisms underlying this phenomenon remain unclear. In this study, we
employed molecular dynamics (MD) simulations on pseudo-phosphorylated tau-MT
complex (residue 199 ~ 312), incorporating structural data from recent
cryo-electron microscopy studies. Simulation results have revealed altered tau
conformations after applying pseudo-phosphorylation. Additionally,
root-mean-square deviation (RMSD) analyses and dimensionality reduction of
dihedral angles revealed key residues responsible for these conformational
shifts


  Many species of microbes cooperate by producing public goods from which they
collectively benefit. However, these populations are under the risk of being
taken over by cheating mutants that do not contribute to the pool of public
goods. Here we present theoretical findings that address how the social
evolution of microbes can be manipulated by external perturbations, to inhibit
or promote the fixation of cheaters. To control social evolution, we determine
the effects of fluid-dynamical properties such as flow rate or boundary
geometry. We also study the social evolutionary consequences of introducing
beneficial or harmful chemicals at steady state and in a time dependent
fashion. We show that by modulating the flow rate and by applying pulsed
chemical signals, we can modulate the spatial structure and dynamics of the
population, in a way that can select for more or less cooperative microbial
populations.


  mRNA translation is a crucial process that leads to protein synthesis in
living cells. Therefore, it is a process that needs to work optimally for a
cell to stay healthy and alive. With advancements in microscopy and novel
experimental techniques, a lot of the intricate details about the translation
mechanism are now known. However, the why and how of this mechanism are still
ill understood, and therefore, is an active area of research. Theoretical
studies of mRNA translation typically view it in terms of the Totally
Asymmetric Simple Exclusion Process or TASEP. Various works have used the TASEP
model in order to study a wide range of phenomena and factors affecting
translation, such as ribosome traffic on an mRNA under noisy (codon-dependent
or otherwise) conditions, ribosome stalling, premature termination, ribosome
reinitiation and dropoff, codon-dependent elongation and competition among mRNA
for ribosomes, among others. In this review, we relay the history and physics
of the translation process in terms of the TASEP framework. In particular, we
discuss the viability and evolution of this model and its limitations while
also formulating the reasons behind its success. Finally, we also identify gaps
in the existing literature and suggest possible extensions and applications
that will lead to a better understanding of ribosome traffic on the mRNA.


  Accurate identification of protein nucleic-acid-binding residues poses a
significant challenge with important implications for various biological
processes and drug design. Many typical computational methods for protein
analysis rely on a single model that could ignore either the semantic context
of the protein or the global 3D geometric information. Consequently, these
approaches may result in incomplete or inaccurate protein analysis. To address
the above issue, in this paper, we present CrossBind, a novel collaborative
cross-modal approach for identifying binding residues by exploiting both
protein geometric structure and its sequence prior knowledge extracted from a
large-scale protein language model. Specifically, our multi-modal approach
leverages a contrastive learning technique and atom-wise attention to capture
the positional relationships between atoms and residues, thereby incorporating
fine-grained local geometric knowledge, for better binding residue prediction.
Extensive experimental results demonstrate that our approach outperforms the
next best state-of-the-art methods, GraphSite and GraphBind, on DNA and RNA
datasets by 10.8/17.3% in terms of the harmonic mean of precision and recall
(F1-Score) and 11.9/24.8% in Matthews correlation coefficient (MCC),
respectively. We release the code at https://github.com/BEAM-Labs/CrossBind.


  Previously, it has been shown that maximum-entropy models of
immune-repertoire sequence can be used to determine a person's vaccination
status. However, this approach has the drawback of requiring a computationally
intensive method to compute each model's partition function ($Z$), the
normalization constant required for calculating the probability that the model
will generate a given sequence. Specifically, the method required generating
approximately $10^{10}$ sequences via Monte-Carlo simulations for each model.
This is impractical for large numbers of models. Here we propose an alternative
method that requires estimating $Z$ this way for only a few models: it then
uses these expensive estimates to estimate $Z$ more efficiently for the
remaining models. We demonstrate that this new method enables the generation of
accurate estimates for 27 models using only three expensive estimates, thereby
reducing the computational cost by an order of magnitude. Importantly, this
gain in efficiency is achieved with only minimal impact on classification
accuracy. Thus, this new method enables larger-scale investigations in
computational immunology and represents a useful contribution to energy-based
modeling more generally.


  The Mountain pine beetle (MPB) is a destructive forest pest that undergoes
boom-bust cycles every 20-40 years. It also has unusual life-history traits,
including a chemosensory affinity for large pine trees, coordinated attacks
facilitated by aggregation and anti-aggregation pheromones, and variable
dispersal that increases in the absence of trees. Using an empirically
calibrated model, we show that MPB's distinctive life-history characteristics
are responsible for its distinctive population dynamics. Specifically, three
life-history traits are necessary for episodic boom-bust dynamics:
density-dependent dispersal, where beetles disperse once most of the large
trees have been killed; an Allee effect, which requires a threshold number of
beetles to overcome tree defenses; and a short generation time, notably in
comparison with their pine tree hosts. In addition to explaining the
qualitative behavior of MPB dynamics, these three ingredients explain residual
tree density, the duration of outbreaks, and the average waiting time between
outbreaks. The peak number of beetles, believed to be the primary factor
driving range expansion, is mostly a consequence of MPB's constitutive high
fecundity. However, two life history traits -- MPB's size-dependent fecundity
and preference for large trees -- are responsible for 25% of the peak number of
beetles. By identifying patterns across the extensive MPB literature, and by
integrating data across multiple sources, we develop a highly accurate
description of MPB's density-dependent interactions with lodgepole pine. We
conclude with a simplified mechanistic model, distilled to two difference
equations, effectively capturing the essence of MPB outbreak dynamics.


  Determining whether a molecule can be synthesized is crucial for many aspects
of chemistry and drug discovery, allowing prioritization of experimental work
and ranking molecules in de novo design tasks. Existing scoring approaches to
assess synthetic feasibility struggle to extrapolate to out-of-distribution
chemical spaces or fail to discriminate based on minor differences such as
chirality that might be obvious to trained chemists. This work aims to address
these limitations by introducing the Focused Synthesizability score (FSscore),
which learns to rank structures based on binary preferences using a graph
attention network. First, a baseline trained on an extensive set of
reactant-product pairs is established that subsequently is fine-tuned with
expert human feedback on a chemical space of interest. Fine-tuning on focused
datasets improves performance on these chemical scopes over the pre-trained
model exhibiting moderate performance and generalizability. This enables
distinguishing hard- from easy-to-synthesize molecules and improving the
synthetic accessibility of generative model outputs. On very complex scopes
with limited labels achieving satisfactory gains remains challenging. The
FSscore showcases how human expert feedback can be utilized to optimize the
assessment of synthetic feasibility for a variety of applications.


  The Price equation provides a formal account of selection building on a
right-total mapping between two classes of individuals, that is usually
interpreted as a parent-offspring relation. This paper presents a new
formulation of the Price equation in terms of fuzzy set-mappings to account for
structures where the targets of selection may vary in the degree to which they
belong to the classes of "parents" or "offspring," and in the degree to which
these two classes of individuals are related. The fuzzy set formulation widens
the scope of the Price equation such that it equally applies to natural
selection, cultural selection, operant selection, or any combination of
different types of selection.


  Learning representations of molecular structures using deep learning is a
fundamental problem in molecular property prediction tasks. Molecules
inherently exist in the real world as three-dimensional structures;
furthermore, they are not static but in continuous motion in the 3D Euclidean
space, forming a potential energy surface. Therefore, it is desirable to
generate multiple conformations in advance and extract molecular
representations using a 4D-QSAR model that incorporates multiple conformations.
However, this approach is impractical for drug and material discovery tasks
because of the computational cost of obtaining multiple conformations. To
address this issue, we propose a pre-training method for molecular GNNs using
an existing dataset of molecular conformations to generate a latent vector
universal to multiple conformations from a 2D molecular graph. Our method,
called Boltzmann GNN, is formulated by maximizing the conditional marginal
likelihood of a conditional generative model for conformations generation. We
show that our model has a better prediction performance for molecular
properties than existing pre-training methods using molecular graphs and
three-dimensional molecular structures.


  Interrogating the evolution of biological changes at early stages of life
requires longitudinal profiling of molecules, such as DNA methylation, which
can be challenging with children. We introduce a probabilistic and longitudinal
machine learning framework based on multi-mean Gaussian processes (GPs),
accounting for individual and gene correlations across time. This method
provides future predictions of DNA methylation status at different individual
ages while accounting for uncertainty. Our model is trained on a birth cohort
of children with methylation profiled at ages 0-4, and we demonstrated that the
status of methylation sites for each child can be accurately predicted at ages
5-7. We show that methylation profiles predicted by multi-mean GPs can be used
to estimate other phenotypes, such as epigenetic age, and enable comparison to
other health measures of interest. This approach encourages epigenetic studies
to move towards longitudinal design for investigating epigenetic changes during
development, ageing and disease progression.


  Survival is a key life history parameter that can inform management decisions
and basic life history research. Because true survival is often confounded with
permanent and temporary emigration from the study area, many studies resort to
estimating apparent survival (i.e., probability of surviving and remaining
inside the study area), which can be much lower than true survival for highly
mobile species. One popular method for estimating true survival is the Barker
joint live-recapture/live-resight (JLRLR) model, which combines capture data
from a study area (hereafter the capture site) with resighting data from a
broader geographic area. This model assumes resightings are spatially
widespread to encompass the species range, an assumption that may not be
realistic for widely dispersing species. Here we use simulation to evaluate
survival bias from a JLRLR model under study design scenarios that differ in
the site selection for resights: global, random, fixed including the capture
site, and fixed excluding the capture site. Simulation results indicate that
fixed designs that included the capture site showed negative survival bias,
whereas fixed designs that excluded the capture site exhibited positive
survival bias. The magnitude of the bias was dependent on movement and
survival, where scenarios with high survival and frequent movement had minimal
bias. We developed a multistate version of the JLRLR and demonstrated
reductions in survival bias compared to the single-state version for most
designs. To minimize bias, our results suggest: 1) using a random resight
design when feasible and global sampling is not possible, 2) using the
multistate JLRLR model when appropriate, 3) including the capture site in the
resight sampling frame when possible, and 4) reporting survival as apparent
survival if fixed sites are used for resight with the single state JLRLR model.


  B cell receptors (BCRs) play a crucial role in recognizing and fighting
foreign antigens. High-throughput sequencing enables in-depth sampling of the
BCRs repertoire after immunization. However, only a minor fraction of BCRs
actively participate in any given infection. To what extent can we accurately
identify antigen-specific sequences directly from BCRs repertoires? We present
a computational method grounded on sequence similarity, aimed at identifying
statistically significant responsive BCRs. This method leverages well-known
characteristics of affinity maturation and expected diversity. We validate its
effectiveness using longitudinally sampled human immune repertoire data
following influenza vaccination and Sars-CoV-2 infections. We show that
different lineages converge to the same responding CDR3, demonstrating
convergent selection within an individual. The outcomes of this method hold
promise for application in vaccine development, personalized medicine, and
antibody-derived therapeutics.


  Reproduction numbers play a fundamental role in population dynamics. For
age-structured models, these quantities are typically defined as spectral
radius of operators acting on infinite dimensional spaces. As a result, their
analytical computation is hardly achievable without additional assumptions on
the model coefficients (e.g., separability of age-specific transmission rates)
and numerical approximations are needed. In this paper we introduce a general
numerical approach, based on pseudospectral collocation of the relevant
operators, for approximating the reproduction numbers of a class of
age-structured models with finite life span. To our knowledge, this is the
first numerical method that allows complete flexibility in the choice of the
``birth'' and ``transition'' processes, which is made possible by working with
an equivalent problem for the integrated state. We discuss applications to
epidemic models with continuous rates, as well as models with piecewise
continuous rates estimated from real data, illustrating how the method can
compute different reproduction numbers-including the basic and the type
reproduction number as special cases-by considering different interpretations
of the age variable (e.g., chronological age, infection age, disease age) and
the transmission terms (e.g., horizontal and vertical transmission).


  Purpose: To assess Generative Pre-trained Transformer version 4's (GPT-4)
ability to classify articles containing functional evidence relevant to
assessments of variant pathogenicity.
  Results: GPT-4 settings and prompts were trained on a set of 45 articles and
genetic variants. A final test set of 72 manually classified articles and
genetic variants were then processed using two prompts. The prompts asked GPT-4
to supply all functional evidence present in an article for a variant or
indicate that no functional evidence is present. For articles with having
functional evidence, a second prompt asked GPT-4 to classify the evidence into
pathogenic, benign, intermediate, and inconclusive categories. The first prompt
identified articles with variant-level functional evidence with 87% sensitivity
and 89% positive predictive value (PPV). Five of 26 articles with no functional
data were indicated as having functional evidence by GPT-4. For variants with
functional assays present as determined by both manual review and GPT-4, the
sensitivity and PPV of GPT-4 prompt concordance was: Pathogenic (92% sensitive
and 73% PPV), Intermediate or Inconclusive (67% sensitive and 93% PPV), Benign
(100% sensitive and 73% PPV).
  Conclusion: The GPT-4 prompts detected the presence or absence of a
functional assay with high sensitivity and PPV, and articles with unambiguous
evidence supporting a benign or pathogenic classification with high sensitivity
and reasonable PPV. Our prompts detected papers with intermediate or
inconclusive evidence with lower sensitivity but high PPV. Our results support
that GPT-4 may be useful in variant classification workflows by enabling
prioritization of articles for review that are likely to have functional
evidence supporting or refuting pathogenicity, but not that GPT-4 is capable of
fully automating the genetics literature review component of variant
classification.


  Reaction-diffusion models are often used to describe biological invasion,
where populations of individuals that undergo random motility and proliferation
lead to moving fronts. Many models of biological invasion are extensions of the
Fisher-KPP model that describes the evolution of a 1D population density as a
result of linear diffusion and logistic growth. In 2020 Fadai introduced a new
model of biological invasion that was formulated as a moving boundary problem
with a nonlinear degenerate diffusive flux. Fadai's model leads to travelling
wave solutions with infinitely steep, well-defined fronts at the moving
boundary, and the model has the mathematical advantage of being analytically
tractable in certain parameter limits. We aim to provide general insight by
first presenting two key extensions by considering: (i) generalised nonlinear
degenerate diffusion with flux; and, (ii) solutions describing both biological
invasion, and biological recession. We establish the existence of travelling
wave solutions for these two extensions, and then consider stability of the
travelling wave solutions by introducing a lateral perturbation of the
travelling wavefront. Full 2D time-dependent level-set numerical solutions
indicate that invasive travelling waves are stable to small lateral
perturbations, whereas receding travelling waves are unstable. These
preliminary numerical observations are corroborated through a linear stability
analysis that gives more formal insight into short time growth/decay of
wavefront perturbation amplitude.


  Despite the great popularity of virtual screening of existing compound
libraries, the search for new potential drug candidates also takes advantage of
generative protocols, where new compound suggestions are enumerated using
various algorithms. To increase the activity potency of generative approaches,
they have recently been coupled with molecular docking, a leading methodology
of structure-based drug design. In this review, we summarize progress since
docking-based generative models emerged. We propose a new taxonomy for these
methods and discuss their importance for the field of computer-aided drug
design. In addition, we discuss the most promising directions for further
development of generative protocols coupled with docking.


  Antifragility characterizes the benefit of a dynamical system derived from
the variability in environmental perturbations. Antifragility carries a precise
definition that quantifies a system's output response to input variability.
Systems may respond poorly to perturbations (fragile) or benefit from
perturbations (antifragile). In this manuscript, we review a range of
applications of antifragility theory in technical systems (e.g., traffic
control, robotics) and natural systems (e.g., cancer therapy, antibiotics).
While there is a broad overlap in methods used to quantify and apply
antifragility across disciplines, there is a need for precisely defining the
scales at which antifragility operates. Thus, we provide a brief general
introduction to the properties of antifragility in applied systems and review
relevant literature for both natural and technical systems' antifragility. We
frame this review within three scales common to technical systems: intrinsic
(input-output nonlinearity), inherited (extrinsic environmental signals), and
interventional (feedback control), with associated counterparts in biological
systems: ecological (homogeneous systems), evolutionary (heterogeneous
systems), and interventional (control). We use the common noun in designing
systems that exhibit antifragile behavior across scales and guide the reader
along the spectrum of
fragility-adaptiveness-resilience-robustness-antifragility, the principles
behind it, and its practical implications.


  The single-cell RNA sequencing (scRNA-seq) technology enables researchers to
study complex biological systems and diseases with high resolution. The central
challenge is synthesizing enough scRNA-seq samples; insufficient samples can
impede downstream analysis and reproducibility. While various methods have been
attempted in past research, the resulting scRNA-seq samples were often of poor
quality or limited in terms of useful specific cell subpopulations. To address
these issues, we propose a novel method called Single-Cell Latent Diffusion
(SCLD) based on the Diffusion Model. This method is capable of synthesizing
large-scale, high-quality scRNA-seq samples, including both 'holistic' or
targeted specific cellular subpopulations within a unified framework. A
pre-guidance mechanism is designed for synthesizing specific cellular
subpopulations, while a post-guidance mechanism aims to enhance the quality of
scRNA-seq samples. The SCLD can synthesize large-scale and high-quality
scRNA-seq samples for various downstream tasks. Our experimental results
demonstrate state-of-the-art performance in cell classification and data
distribution distances when evaluated on two scRNA-seq benchmarks.
Additionally, visualization experiments show the SCLD's capability in
synthesizing specific cellular subpopulations.


  The surge in high-throughput omics data has reshaped the landscape of
biological research, underlining the need for powerful, user-friendly data
analysis and interpretation tools. This paper presents GenoCraft, a web-based
comprehensive software solution designed to handle the entire pipeline of omics
data processing. GenoCraft offers a unified platform featuring advanced
bioinformatics tools, covering all aspects of omics data analysis. It
encompasses a range of functionalities, such as normalization, quality control,
differential analysis, network analysis, pathway analysis, and diverse
visualization techniques. This software makes state-of-the-art omics data
analysis more accessible to a wider range of users. With GenoCraft, researchers
and data scientists have access to an array of cutting-edge bioinformatics
tools under a user-friendly interface, making it a valuable resource for
managing and analyzing large-scale omics data. The API with an interactive web
interface is publicly available at https://genocraft.stanford. edu/. We also
release all the codes in https://github.com/futianfan/GenoCraft.


  Chromatin is a polymer complex of DNA and proteins that regulates gene
expression. The three-dimensional structure and organization of chromatin
controls DNA transcription and replication. High-throughput chromatin
conformation capture techniques generate Hi-C maps that can provide insight
into the 3D structure of chromatin. Hi-C maps can be represented as a symmetric
matrix where each element represents the average contact probability or number
of contacts between two chromatin loci. Previous studies have detected
topologically associating domains (TADs), or self-interacting regions in Hi-C
maps within which the contact probability is greater than that outside the
region. Many algorithms have been developed to identify TADs within Hi-C maps.
However, most TAD identification algorithms are unable to identify nested or
overlapping TADs and for a given Hi-C map there is significant variation in the
location and number of TADs identified by different methods. We develop a novel
method, KerTAD, using a kernel-based technique from computer vision and image
processing that is able to accurately identify nested and overlapping TADs. We
benchmark this method against state-of-the-art TAD identification methods on
both synthetic and experimental data sets. We find that KerTAD consistently has
higher true positive rates (TPR) and lower false discovery rates (FDR) than all
tested methods for both synthetic and manually annotated experimental Hi-C
maps. The TPR for KerTAD is also largely insensitive to increasing noise and
sparsity, in contrast to the other methods. We also find that KerTAD is
consistent in the number and size of TADs identified across replicate
experimental Hi-C maps for several organisms. KerTAD will improve automated TAD
identification and enable researchers to better correlate changes in TADs to
biological phenomena, such as enhancer-promoter interactions and disease
states.


  X-ray free-electron lasers (XFELs) offer unique capabilities for measuring
the structure and dynamics of biomolecules, helping us understand the basic
building blocks of life. Notably, high-repetition-rate XFELs enable single
particle imaging (X-ray SPI) where individual, weakly scattering biomolecules
are imaged under near-physiological conditions with the opportunity to access
fleeting states that cannot be captured in cryogenic or crystallized
conditions. Existing X-ray SPI reconstruction algorithms, which estimate the
unknown orientation of a particle in each captured image as well as its shared
3D structure, are inadequate in handling the massive datasets generated by
these emerging XFELs. Here, we introduce X-RAI, an online reconstruction
framework that estimates the structure of a 3D macromolecule from large X-ray
SPI datasets. X-RAI consists of a convolutional encoder, which amortizes pose
estimation over large datasets, as well as a physics-based decoder, which
employs an implicit neural representation to enable high-quality 3D
reconstruction in an end-to-end, self-supervised manner. We demonstrate that
X-RAI achieves state-of-the-art performance for small-scale datasets in
simulation and challenging experimental settings and demonstrate its
unprecedented ability to process large datasets containing millions of
diffraction images in an online fashion. These abilities signify a paradigm
shift in X-ray SPI towards real-time capture and reconstruction.


  Selection among autocatalytic species fundamentally depends on their growth
law: exponential species, whose number of copies grows exponentially, are
mutually exclusive, while sub-exponential ones, whose number of copies grows
polynomially, can coexist. Here we consider competitions between autocatalytic
species with different growth laws and make the simple yet counterintuitive
observation that sub-exponential species can exclude exponential ones while the
reverse is, in principle, impossible. This observation has implications for
scenarios pertaining to the emergence of natural selection.


  Quarantine measures are one of the first lines of defense against the spread
of infectious diseases. However, maintaining these measures over extended
periods can be challenging due to a phenomenon known as quarantine fatigue. In
this paper, we investigate the impact of quarantine fatigue on the spread of
infectious diseases by using an epidemic model on random networks with cliques.
In our model, susceptible individuals can be quarantined up to $n$ times, after
which they stop complying with quarantine orders due to fatigue. Our results
show that quarantine fatigue may induce a regime in which increasing the
probability of detecting and isolating infected individuals (along with their
close contacts) could subsequently increase the expected number of cases at the
end of an outbreak. Moreover, we observe that quarantine fatigue can trigger an
abrupt phase transition at the critical reproduction number $R_0=1$. Finally,
we explore a scenario where a non-negligible number of individuals are infected
at the beginning of an epidemic, and our results show that, depending on the
value of $n$, an abrupt transition between a controlled epidemic and a large
epidemic event can occur for $R_0<1$.


  PlayMolecule Viewer is a web-based data visualization toolkit designed to
streamline the exploration of data resulting from structural bioinformatics or
computer-aided drug design efforts. By harnessing state-of-the-art web
technologies such as WebAssembly, PlayMolecule Viewer integrates powerful
Python libraries directly within the browser environment, which enhances its
capabilities of managing multiple types of molecular data. With its intuitive
interface, it allows users to easily upload, visualize, select, and manipulate
molecular structures and associated data. The toolkit supports a wide range of
common structural file formats and offers a variety of molecular
representations to cater to different visualization needs. PlayMolecule Viewer
is freely accessible at open.playmolecule.org, ensuring accessibility and
availability to the scientific community and beyond.


  Antimicrobial resistance to drugs (AMR), a global threat to human and animal
health, is often regarded as resulting from a cooperative behaviour, leading to
the coexistence of drug-resistant and drug-sensitive cells in large communities
and static environments. However, microbial populations generally evolve in
volatile environments, yielding sudden changes in the available amount of
nutrients and toxins. These, together with demographic fluctuations (birth and
death events), have been recognised to drastically alter the population size
and strain coexistence, as well as the ability of one strain to fixate and take
over the entire population. Motivated by the need to better understand the
evolution of AMR, we study a population of time-varying size consisting of two
competing strains, one drug-resistant and one drug-sensitive, subject to
demographic and environmental variability. This is modelled by a binary
carrying capacity randomly switching between mild and harsh environmental
conditions, and driving the fluctuating volume (total amount of nutrients and
antimicrobials at fixed concentration), and thus the size of the community
(number of resistant and sensitive cells). We assume that AMR is a shared
public good when the concentration of resistant cells exceeds a fixed
"concentration cooperation threshold", above which the sensitive strain has a
growth advantage, whereas resistant cells dominate below it. Using
computational means, and devising an analytical treatment (built on suitable
quenched and annealed averaging procedures), we fully characterise the
influence of fluctuations on the eco-evolutionary dynamics of AMR, and notably
obtain specific strain fixation and long-lasting coexistence probabilities as a
function of the environmental variation rate and cooperation threshold.


  Efficient allocation is important in nature and human society where
individuals often compete for finite resources. The Minority Game is perhaps
the simplest model that provides deep insights into how human coordinate to
maximize the resource utilization. However, this model assumes the static
strategies that are provided a priori, failing to capture their adaptive
nature. Here, we turn to the paradigm of reinforcement learning, where
individuals' strategies are evolving by evaluating both the past experience and
rewards in the future. Specifically, we adopt the Q-learning algorithm, each
player is endowed with a Q-table that guides their decision-making. We reveal
that the population is able to reach the optimal allocation when individuals
appreciate both the past experience and rewards in the future, and they are
able to balance the exploitation of their Q-tables and the exploration by
randomly acting. The optimal allocation is ruined when individuals tend to use
either exploitation-only or exploration-only, where only partial coordination
and even anti-coordination are observed. Mechanism analysis reveals that a
moderate level of exploration can escape local minimums of metastable periodic
states, and reaches the optimal coordination as the global minimum.
Interestingly, the optimal coordination is underlined by a symmetry-breaking of
action preferences, where nearly half of the population choose one side while
the other half prefer the other side. The emergence of optimal coordination is
robust to the population size and other game parameters. Our work therefore
provides a natural solution to the Minority Game and sheds insights into the
resource allocation problem in general. Besides, our work demonstrates the
potential of the proposed reinforcement learning paradigm in deciphering many
puzzles in the socio-economic context.


  We investigate how a population's natural and vaccine immunity affects the
competitive balance between two strains of an infectious disease with different
epidemiological characteristics. Specifically, we consider the case where one
strain is more transmissible and the other strain is more immune-resistant. The
competition of these strains is modeled by two SIR-type models which
incorporate waning natural immunity and which have distinct mechanisms for
vaccine immunity. Waning immunity is implemented as a gamma-distributed delay,
which is analyzed using the linear chain trick to transform the delay
differential equation systems into a system of ordinary differential equations.
Our analysis shows that vaccination has a significant effect on the competitive
balance between two strains, potentially leading to dramatic flips from one
strain dominating in the population to the other. We also show that which
strain gains an advantage as a population's immunity level increases depends
upon the integration between the mechanisms of natural and vaccine immunity.
The results of this paper are consequently relevant for public policy.


  Genome-Wide Association Studies (GWAS) face unique challenges in the era of
big genomics data, particularly when dealing with ultra-high-dimensional
datasets where the number of genetic features significantly exceeds the
available samples. This paper introduces an extension to the feature selection
methodology proposed by Mirzaei et al. (2020), specifically tailored to tackle
the intricacies associated with ultra-high-dimensional GWAS data. Our extended
approach enhances the original method by introducing a Frobenius norm penalty
into the student network, augmenting its capacity to adapt to scenarios
characterized by a multitude of features and limited samples. Operating
seamlessly in both supervised and unsupervised settings, our method employs two
key neural networks. The first leverages an autoencoder or supervised
autoencoder for dimension reduction, extracting salient features from the
ultra-high-dimensional genomic data. The second network, a regularized
feed-forward model with a single hidden layer, is designed for precise feature
selection. The introduction of the Frobenius norm penalty in the student
network significantly boosts the method's resilience to the challenges posed by
ultra-high-dimensional GWAS datasets. Experimental results showcase the
efficacy of our approach in feature selection for GWAS data. The method not
only handles the inherent complexities of ultra-high-dimensional settings but
also demonstrates superior adaptability to the nuanced structures present in
genomics data. The flexibility and versatility of our proposed methodology are
underscored by its successful performance across a spectrum of experiments.


  Accurate and robust prediction of drug-target interactions (DTIs) plays a
vital role in drug discovery. Despite extensive efforts have been invested in
predicting novel DTIs, existing approaches still suffer from insufficient
labeled data and cold start problems. More importantly, there is currently a
lack of studies focusing on elucidating the mechanism of action (MoA) between
drugs and targets. Distinguishing the activation and inhibition mechanisms is
critical and challenging in drug development. Here, we introduce a unified
framework called DTIAM, which aims to predict interactions, binding affinities,
and activation/inhibition mechanisms between drugs and targets. DTIAM learns
drug and target representations from large amounts of label-free data through
self-supervised pre-training, which accurately extracts the substructure and
contextual information of drugs and targets, and thus benefits the downstream
prediction based on these representations. DTIAM achieves substantial
performance improvement over other state-of-the-art methods in all tasks,
particularly in the cold start scenario. Moreover, independent validation
demonstrates the strong generalization ability of DTIAM. All these results
suggested that DTIAM can provide a practically useful tool for predicting novel
DTIs and further distinguishing the MoA of candidate drugs. DTIAM, for the
first time, provides a unified framework for accurate and robust prediction of
drug-target interactions, binding affinities, and activation/inhibition
mechanisms.


  Dengue incidence forecasting using hybrid models has been surging in the data
rich world. Hybridization of statistical time series forecasting models and
machine learning models are explored for dengue forecasting with different
degrees of success. In this paper, we propose a multivariate expansion of the
hybrid ARIMA-ARNN model. The main motivation is to propose a novel
hybridization and apply it to dengue outbreak prediction. The asymptotic
stationarity of the proposed model has been established. We check the
forecasting capability and robustness of the forecasts through numerical
experiments. State-of-the-art forecasting models for multivariate time series
data are compared with the proposed model using accuracy metrics. Dengue
incidence data from San Juan and Iquitos are utilized along with rainfall as an
exogenous variable. Results indicate that the proposed model improves the
ARIMAX forecasts in some situations and closely follows it otherwise. The
theoretical as well as experimental results reinforce that the proposed model
has the potential to act as a candidate for early warning of dengue outbreaks.
The proposed model can be readily generalized to incorporate more exogenous
variables and also applied to other time series forecasting problems wherever
exogenous variable(s) are available.


  Individuals with suspected rare genetic disorders often undergo multiple
clinical evaluations, imaging studies, laboratory tests and genetic tests, to
find a possible answer over a prolonged period of multiple years. Addressing
this diagnostic odyssey thus have substantial clinical, psychosocial, and
economic benefits. Many rare genetic diseases have distinctive facial features,
which can be used by artificial intelligence algorithms to facilitate clinical
diagnosis, in prioritizing candidate diseases to be further examined by lab
tests or genetic assays, or in helping the phenotype-driven reinterpretation of
genome/exome sequencing data. However, existing methods using frontal facial
photo were built on conventional Convolutional Neural Networks (CNNs), rely
exclusively on facial images, and cannot capture non-facial phenotypic traits
and demographic information essential for guiding accurate diagnoses. Here we
introduce GestaltMML, a multimodal machine learning (MML) approach solely based
on the Transformer architecture. It integrates the facial images, demographic
information (age, sex, ethnicity), and clinical notes of patients to improve
prediction accuracy. Furthermore, we also introduce GestaltGPT, a GPT-based
methodology with few-short learning capacities that exclusively harnesses
textual inputs using a range of large language models (LLMs) including Llama 2,
GPT-J and Falcon. We evaluated these methods on a diverse range of datasets,
including 449 diseases from the GestaltMatcher Database, several in-house
datasets on Beckwith-Wiedemann syndrome, Sotos syndrome, NAA10-related syndrome
(neurodevelopmental syndrome) and others. Our results suggest that
GestaltMML/GestaltGPT effectively incorporate multiple modalities of data,
greatly narrow down candidate genetic diagnosis of rare diseases, and may
facilitate the reinterpretation of genome/exome sequencing data.


  Therapeutic peptides represent a unique class of pharmaceutical agents
crucial for the treatment of human diseases. Recently, deep generative models
have exhibited remarkable potential for generating therapeutic peptides, but
they only utilize sequence or structure information alone, which hinders the
performance in generation. In this study, we propose a Multi-Modal Contrastive
Diffusion model (MMCD), fusing both sequence and structure modalities in a
diffusion framework to co-generate novel peptide sequences and structures.
Specifically, MMCD constructs the sequence-modal and structure-modal diffusion
models, respectively, and devises a multi-modal contrastive learning strategy
with intercontrastive and intra-contrastive in each diffusion timestep, aiming
to capture the consistency between two modalities and boost model performance.
The inter-contrastive aligns sequences and structures of peptides by maximizing
the agreement of their embeddings, while the intra-contrastive differentiates
therapeutic and non-therapeutic peptides by maximizing the disagreement of
their sequence/structure embeddings simultaneously. The extensive experiments
demonstrate that MMCD performs better than other state-of-theart deep
generative methods in generating therapeutic peptides across various metrics,
including antimicrobial/anticancer score, diversity, and peptide-docking.


  The Hill function is relevant for describing enzyme binding and other
processes in gene regulatory networks. Despite its theoretical foundation, it
is often used empirically as a useful fitting function. Theoretical predictions
suggest that the Hill coefficient should be an integer; however, it is often
assigned a fractional value. This study aims to show that the use of fractional
Hill coefficients actually indicates the presence of intermediate processes or,
in some cases, the effect of noise on these systems. The deterministic
approximation of enzymatic processes leads to the derivation of the Hill
function, which can be expanded around the noise magnitude to derive mesoscopic
corrections. This study establishes relationships between the intermediate
processes and the decimal Hill coefficient, both with and without fluctuations.
This outcome contributes to a deeper understanding of the underlying processes
associated with the fractional Hill coefficient, while also enabling the
prediction of an effective value of the Hill coefficient from the underlying
mechanism, allowing us to have a simplified description of complex systems.


  We report results of a restricted mathematical computation that may be
relevant in various contexts, including the evolution of sexual reproduction.
The computation is an iterative probabilistic replication of sets of real
numbers. Given a set of random values sampled from a certain interval, we
replicate it probabilistically by encoding it as its arithmetic mean and then
sampling a random same-size set of real numbers using uniform distribution over
the same interval length around that mean. The encoded set and its code are
considered the parent of the sampled one. We first sample a cohort of several
such sets around an initial code, and probabilistically replicate each set
once, creating a new same-size cohort, which we then subject to the same
process. In a bi-parental variant of the above mono-parental process, we use a
pair of codes rather than a single one for sampling a new set: each parent
contributes half of the elements of the child set. Neither variant depends on
judging properties as harmful, or on selection: all individuals are paired and
reproduced. Comparing the two variants over many iterations, bi-parental
processing shows enhanced fidelity and consistency, where fidelity is defined
as the mean of the distances between a cohort's codes and the initial code, and
consistency is the standard deviation of code values within a cohort. We do not
attempt here to directly mimic natural reproduction, which includes copying and
recombination of DNA codes subject to mutations. Still, the result intuitively
aligns with the idea that sexual reproduction in nature helps preserve, and
perhaps even define, core traits of each species. The cause for the prevalence,
in eukaryotic species, of sexual reproduction over asexual, despite increased
complexity and unpredictability, is subject to debate and diverse theories. The
mathematical phenomenon reported here might be part of the answer.


  Unsupervised learning has become a staple in classical machine learning,
successfully identifying clustering patterns in data across a broad range of
domain applications. Surprisingly, despite its accuracy and elegant simplicity,
unsupervised learning has not been sufficiently exploited in the realm of
phylogenetic tree inference. The main reason for the delay in adoption of
unsupervised learning in phylogenetics is the lack of a meaningful, yet simple,
way of embedding phylogenetic trees into a vector space. Here, we propose the
simple yet powerful split-weight embedding which allows us to fit standard
clustering algorithms to the space of phylogenetic trees. We show that our
split-weight embedded clustering is able to recover meaningful evolutionary
relationships in simulated and real (Adansonia baobabs) data.


  Recent advancements in immune sequencing and experimental techniques are
generating extensive T cell receptor (TCR) repertoire data, enabling the
development of models to predict TCR binding specificity. Despite the
computational challenges due to the vast diversity of TCRs and epitopes,
significant progress has been made. This paper discusses the evolution of the
computational models developed for this task, with a focus on machine learning
efforts, including the early unsupervised clustering approaches, supervised
models, and the more recent applications of Protein Language Models (PLMs). We
critically assess the most prominent models in each category, and discuss
recurrent challenges, such as the lack of generalization to new epitopes,
dataset biases, and biases in the validation design of the models.
  Furthermore, our paper discusses the transformative role of transformer-based
protein models in bioinformatics. These models, pretrained on extensive
collections of unlabeled protein sequences, can convert amino acid sequences
into vectorized embeddings that capture important biological properties. We
discuss recent attempts to leverage PLMs to deliver very competitive
performances in TCR-related tasks. Finally, we address the pressing need for
improved interpretability in these often opaque models, proposing strategies to
amplify their impact in the field.


  Single-cell RNA sequencing (scRNA-seq) enables researchers to analyze gene
expression at single-cell level. One important task in scRNA-seq data analysis
is unsupervised clustering, which helps identify distinct cell types, laying
down the foundation for other downstream analysis tasks. In this paper, we
propose a novel method called Cluster-aware Iterative Contrastive Learning
(CICL in short) for scRNA-seq data clustering, which utilizes an iterative
representation learning and clustering framework to progressively learn the
clustering structure of scRNA-seq data with a cluster-aware contrastive loss.
CICL consists of a Transformer encoder, a clustering head, a projection head
and a contrastive loss module. First, CICL extracts the feature vectors of the
original and augmented data by the Transformer encoder. Then, it computes the
clustering centroids by K-means and employs the student t-distribution to
assign pseudo-labels to all cells in the clustering head. The projection-head
uses a Multi-Layer Perceptron (MLP) to obtain projections of the augmented
data. At last, both pseudo-labels and projections are used in the contrastive
loss to guide the model training. Such a process goes iteratively so that the
clustering result becomes better and better. Extensive experiments on 25 real
world scRNA-seq datasets show that CICL outperforms the SOTA methods.
Concretely, CICL surpasses the existing methods by from 14% to 280%, and from
5% to 133% on average in terms of performance metrics ARI and NMI respectively.


  Recording the undersea three-dimensional bioacoustic sound field in real-time
promises major benefits to marine behavior studies. We describe a novel
hydrophone array -- the hydroambiphone (HAP) -- that adapts ambisonic
spatial-audio theory to sound propagation in ocean waters to realize many of
these benefits through spatial localization and acoustic immersion. Deploying
it to monitor the humpback whales (Megaptera novaeangliae) of southeast Alaska
demonstrates that HAP recording provides a qualitatively-improved experience of
their undersea behaviors; revealing, for example, new aspects of social
coordination during bubble-net feeding. On the practical side, spatialized
hydrophone recording greatly reduces post-field analytical and computational
challenges -- such as the "cocktail party problem" of distinguishing single
sources in a complicated and crowded auditory environment -- that are common to
field recordings. On the scientific side, comparing the HAP's capabilities to
single-hydrophone and nonspatialized recordings yields new insights into the
spatial information that allows animals to thrive in complex acoustic
environments. Spatialized bioacoustics markedly improves access to the
humpbacks' undersea acoustic environment and expands our appreciation of their
rich vocal lives.


  Molecular property prediction (MPP) is a fundamental but challenging task in
the computer-aided drug discovery process. More and more recent works employ
different graph-based models for MPP, which have made considerable progress in
improving prediction performance. However, current models often ignore
relationships between molecules, which could be also helpful for MPP. For this
sake, in this paper we propose a graph structure learning (GSL) based MPP
approach, called GSL-MPP. Specifically, we first apply graph neural network
(GNN) over molecular graphs to extract molecular representations. Then, with
molecular fingerprints, we construct a molecular similarity graph (MSG).
Following that, we conduct graph structure learning on the MSG (i.e.,
molecule-level graph structure learning) to get the final molecular embeddings,
which are the results of fusing both GNN encoded molecular representations and
the relationships among molecules, i.e., combining both intra-molecule and
inter-molecule information. Finally, we use these molecular embeddings to
perform MPP. Extensive experiments on seven various benchmark datasets show
that our method could achieve state-of-the-art performance in most cases,
especially on classification tasks. Further visualization studies also
demonstrate the good molecular representations of our method.


  Recent findings show that single, non-neuronal cells are also able to learn
signalling responses developing cellular memory. In cellular learning nodes of
signalling networks strengthen their interactions e.g. by the conformational
memory of intrinsically disordered proteins, protein translocation, miRNAs,
lncRNAs, chromatin memory and signalling cascades. This can be described by a
generalized, unicellular Hebbian learning process, where those signalling
connections, which participate in learning, become stronger. Here we review
those scenarios, where cellular signalling is not only repeated in a few times
(when learning occurs), but becomes too frequent, too large, or too complex and
overloads the cell. This leads to desensitisation of signalling networks by
decoupling signalling components, receptor internalization, and consequent
downregulation. These molecular processes are examples of anti-Hebbian learning
and forgetting of signalling networks. Stress can be perceived as signalling
overload inducing the desensitisation of signalling pathways. Aging occurs by
the summative effects of cumulative stress downregulating signalling. We
propose that cellular learning desensitisation, stress and aging may be placed
along the same axis of more and more intensive (prolonged or repeated)
signalling. We discuss how cells might discriminate between repeated and
unexpected signals, and highlight the Hebbian and anti-Hebbian mechanisms
behind the fold-change detection in the NF-\k{appa}B signalling pathway. We
list drug design methods using Hebbian learning (such as chemically-induced
proximity) and clinical treatment modalities inducing (cancer, drug allergies)
desensitisation or avoiding drug-induced desensitisation. A better
discrimination between cellular learning, desensitisation and stress may open
novel directions in drug design, e.g., helping to overcome drug-resistance.


  Boolean networks have been widely used in many areas of science and
engineering to represent various dynamical behaviour. In systems biology, they
became useful tools to study the dynamical characteristics of large-scale
biomolecular networks and there have been a number of studies to develop
efficient ways of finding steady states or cycles of Boolean network models. On
the other hand, there has been little attention to analyzing the dynamic
properties of the network structure itself. Here, we present a systematic way
to study such properties by introducing a recursive self-composite of the logic
update rules. Of note, we found that all Boolean update rules actually have
repeated logic structures underneath. This repeated nature of Boolean networks
reveals interesting algebraic properties embedded in the networks. We found
that each converged logic leads to the same states, called kernel states. As a
result, the longest-length period of states cycle turns out to be equal to the
number of converged logics in the logic cycle. Based on this, we propose a
leaping and filling algorithm to avoid any possible large string explosions
during the self-composition procedures. Finally, we demonstrate how the
proposed approach can be used to reveal interesting hidden properties using
Boolean network examples of a simple network with a long feedback structure, a
T-cell receptor network and a cancer network.


  This work focuses on the task of property targeting: that is, generating
molecules conditioned on target chemical properties to expedite candidate
screening for novel drug and materials development. DiGress is a recent
diffusion model for molecular graphs whose distinctive feature is allowing
property targeting through classifier-based (CB) guidance. While CB guidance
may work to generate molecular-like graphs, we hint at the fact that its
assumptions apply poorly to the chemical domain. Based on this insight we
propose a classifier-free DiGress (FreeGress), which works by directly
injecting the conditioning information into the training process. CF guidance
is convenient given its less stringent assumptions and since it does not
require to train an auxiliary property regressor, thus halving the number of
trainable parameters in the model. We empirically show that our model yields up
to 79% improvement in Mean Absolute Error with respect to DiGress on property
targeting tasks on QM9 and ZINC-250k benchmarks. As an additional contribution,
we propose a simple yet powerful approach to improve chemical validity of
generated samples, based on the observation that certain chemical properties
such as molecular weight correlate with the number of atoms in molecules.


  Wright-Fisher diffusions and their dual ancestral graphs occupy a central
role in the study of allele frequency change and genealogical structure, and
they provide expressions, explicit in some special cases but generally
implicit, for the sampling probability, a crucial quantity in inference. Under
a finite-allele mutation model, with possibly parent-dependent mutation, we
consider the asymptotic regime where the selective advantage of one allele
grows to infinity, while the other parameters remain fixed. In this regime, we
show that the Wright-Fisher diffusion can be approximated either by a Gaussian
process or by a process whose components are independent continuous-state
branching processes with immigration, aligning with analogous results for
Wright-Fisher models but employing different methods. While the first process
becomes degenerate at stationarity, the latter does not and provides a simple,
analytic approximation for the leading term of the sampling probability.
Furthermore, using another approach based on a recursion formula, we
characterise all remaining terms to provide a full asymptotic expansion for the
sampling probability. Finally, we study the asymptotic behaviour of the rates
of the block-counting process of the conditional ancestral selection graph and
establish an asymptotic duality relationship between this and the diffusion.


  The Escherichia coli chemoreceptors form an extensive array that achieves
cooperative and adaptive sensing of extracellular signals. The receptors
control the activity of histidine kinase CheA, which drives a non-equilibrium
phosphorylation-dephosphorylation reaction cycle for response regulator CheY.
Recent single-cell FRET measurements revealed that kinase activity of the array
spontaneously switches between active and inactive states, with asymmetric
switching times that signify time-reversal symmetry breaking in the underlying
dynamics. Here, we show that the asymmetric switching dynamics can be explained
by a non-equilibrium lattice model, which considers both the dissipative
reaction cycles of individual core units and the coupling between neighboring
units. The model reveals that large dissipation and near-critical coupling are
required to explain the observed switching dynamics. Microscopically, the
switching time asymmetry originates from irreversible transition paths. The
model shows that strong dissipation enables sensitive and rapid signaling
response by relieving the speed-sensitivity trade-off, which can be tested by
future single-cell experiments. Overall, our model provides a general framework
for studying biological complexes composed of coupled subunits that are
individually driven by dissipative cycles and the rich non-equilibrium physics
within.


  1. Abrupt environmental changes can lead to evolutionary shifts in not only
mean (optimal value), but also variance of descendants in trait evolution.
There are some methods to detect shifts in optimal value but few studies
consider shifts in variance. 2. We use a multi-optima and multi-variance OU
process model to describe the trait evolution process with shifts in both
optimal value and variance and provide analysis of how the covariance between
species changes when shifts in variance occur along the path. 3. We propose a
new method to detect the shifts in both variance and optimal values based on
minimizing the loss function with L1 penalty. We implement our method in a new
R package, ShiVa (Detection of evolutionary shifts in variance). 4. We conduct
simulations to compare our method with the two methods considering only shifts
in optimal values (l1ou; PhylogeneticEM). Our method shows strength in
predictive ability and includes far fewer false positive shifts in optimal
value compared to other methods when shifts in variance actually exist. When
there are only shifts in optimal value, our method performs similarly to other
methods. We applied our method to the cordylid data, ShiVa outperformed l1ou
and phyloEM, exhibiting the highest log-likelihood and lowest BIC.


  Accurately predicting molecular properties is a challenging but essential
task in drug discovery. Recently, many mono-modal deep learning methods have
been successfully applied to molecular property prediction. However, the
inherent limitation of mono-modal learning arises from relying solely on one
modality of molecular representation, which restricts a comprehensive
understanding of drug molecules and hampers their resilience against data
noise. To overcome the limitations, we construct multimodal deep learning
models to cover different molecular representations. We convert drug molecules
into three molecular representations, SMILES-encoded vectors, ECFP
fingerprints, and molecular graphs. To process the modal information,
Transformer-Encoder, bi-directional gated recurrent units (BiGRU), and graph
convolutional network (GCN) are utilized for feature learning respectively,
which can enhance the model capability to acquire complementary and naturally
occurring bioinformatics information. We evaluated our triple-modal model on
six molecule datasets. Different from bi-modal learning models, we adopt five
fusion methods to capture the specific features and leverage the contribution
of each modal information better. Compared with mono-modal models, our
multimodal fused deep learning (MMFDL) models outperform single models in
accuracy, reliability, and resistance capability against noise. Moreover, we
demonstrate its generalization ability in the prediction of binding constants
for protein-ligand complex molecules in the refined set of PDBbind. The
advantage of the multimodal model lies in its ability to process diverse
sources of data using proper models and suitable fusion methods, which would
enhance the noise resistance of the model while obtaining data diversity.


  The problem of evolutionary complexification of life is considered one of the
fundamental aspects in contemporary evolutionary theory. Parasitism is
ubiquitous, inevitable, and arises as soon as the first replicators appear,
even during the prebiotic stages of evolution. Both in theoretical approaches
(computer modeling and analysis) and in real experiments (replication of
biological macromolecules), parasitic processes emerge almost immediately. An
effective way to avoid the elimination of the host-parasite system is through
compartmentalization. In both theory and experiments, the pressure of
parasitism leads to the complexification of the host-parasite system into a
network of cooperative replicators and their parasites. Parasites have the
ability to create niches for new replicators. The co-evolutionary arms race
between defense systems and counter-defense mechanisms among parasites and
hosts can progress for a considerable duration, involving multiple stages, if
not indefinitely.


  Double-strand breaks (DSBs) in DNA are naturally occurring destructive events
in all organisms that may lead to genome instability. Cells employ various
repair methods known as non-homologous end joining (NHEJ), microhomology
mediated end joining (MMEJ), and homology-directed recombination (HDR). These
repair processes may lead to DNA sequence variations (e.g., nucleotide
insertions, deletions, and substitutions) at the location of the break.
Studying DNA DSB repair processes often involves the use of high throughput
sequencing assays to precisely quantify the sequence variations near the break
with software tools. Often methods of assessing and visualizing these data have
not taken into account the full complexity of the sequencing data, such as the
frequency, type, and position of the sequence variations in a single
comprehensive representation. Here we present a method that allows
visualization of the overall variation pattern as well as comparison of these
patterns among experimental setups.