  Sample efficiency is a fundamental challenge in de novo molecular design.
Ideally, molecular generative models should learn to satisfy a desired
objective under minimal oracle evaluations (computational prediction or wet-lab
experiment). This problem becomes more apparent when using oracles that can
provide increased predictive accuracy but impose a significant cost.
Consequently, these oracles cannot be directly optimized under a practical
budget. Molecular generative models have shown remarkable sample efficiency
when coupled with reinforcement learning, as demonstrated in the Practical
Molecular Optimization (PMO) benchmark. Here, we propose a novel algorithm
called Augmented Memory that combines data augmentation with experience replay.
We show that scores obtained from oracle calls can be reused to update the
model multiple times. We compare Augmented Memory to previously proposed
algorithms and show significantly enhanced sample efficiency in an exploitation
task and a drug discovery case study requiring both exploration and
exploitation. Our method achieves a new state-of-the-art in the PMO benchmark
which enforces a computational budget, outperforming the previous best
performing method on 19/23 tasks.


  Directed evolution of proteins has been the most effective method for protein
engineering. However, a new paradigm is emerging, fusing the library generation
and screening approaches of traditional directed evolution with computation
through the training of machine learning models on protein sequence fitness
data. This chapter highlights successful applications of machine learning to
protein engineering and directed evolution, organized by the improvements that
have been made with respect to each step of the directed evolution cycle.
Additionally, we provide an outlook for the future based on the current
direction of the field, namely in the development of calibrated models and in
incorporating other modalities, such as protein structure.


  To develop effective strategies for controlling the spread of the virus and
potential future outbreaks, a deep understanding of disease transmission
dynamics is crucial. This study proposes a modification to existing
mathematical models used to describe the transmission dynamics of COVID-19 with
environmental pathogens, incorporating a variable population, and employing
incommensurate fractional order derivatives in ordinary differential equations.
Our analysis accurately computes the basic reproduction number and demonstrates
the global stability of the disease-free equilibrium point. Numerical
simulations fitted to real data from South Africa show the efficacy of our
proposed model, with fractional models enhancing flexibility. We also provide
reliable values for initial conditions, model parameters, and order
derivatives, and examine the sensitivity of model parameters. Our study
provides valuable insights into COVID-19 transmission dynamics and has the
potential to inform the development of effective control measures and
prevention strategies.


  The emergence of an autocatalytic network from an available set of elements
is a fundamental step in early evolutionary processes, such as the origin of
metabolism. Given a set of elements, the reactions between them (chemical or
otherwise), and certain elements catalysing certain reactions, a Reflexively
Autocatalytic F-generated (RAF) set is a subset $R'$ of reactions that is
self-generating from a given food set, and with each reaction in $R'$ being
catalysed from within $R'$. RAF theory has been applied to various phenomena in
theoretical biology, and a key feature of the approach is that it is possible
to efficiently identify and classify RAFs within large systems. This is
possible because RAFs can be described as the (nonempty) subsets of the
reactions that are the fixed points of an (efficiently computable) interior map
that operates on subsets of reactions. Although the main generic results
concerning RAFs can be derived using just this property, we show that for
systems with at least 12 reactions there are generic results concerning RAFs
that cannot be proven using the interior operator property alone.


  We show that the flagellar beat of bull spermatozoa and Chlamydomonas
Reinhardtii can be modelled by a minimal, geometrically nonlinear
reaction-diffusion system. Model solutions are spatio-temporally animated
patterns describing flagellar bending waves, connecting beating patterns of
cilia and flagella with chemical patterns from classical reaction-diffusion
systems. Instead of chemical species freely reacting and diffusing in space,
our system describes the tug-of-war reaction-kinetics of molecular motors that
are anchored in the flagellar structure, but the shear deformation that they
generate can diffuse away via the bending elasticity of the flagellum.
Synchronization of the reaction-kinetics in neighbouring elements occurs via a
sliding-control mechanism. We derive from first principles the
reaction-diffusion basis of animated patterns as a direct consequence of the
high internal energy dissipation by the flagellum relative to the external
dissipation by the surrounding fluid. By fitting nonlinear, large-amplitude
solutions of a specific motor cross-bridge reaction-kinetics, we show that
reaction-diffusion successfully accounts for beating patterns of both bull
sperm and Chlamydomonas (wild-type and mbo2-mutant). Our results suggest that
the flagellar beat occurs far from equilibrium, in the strongly nonlinear
regime, and that a unified mechanism may exist for dynein molecular motor
control that is regulated by axonemal sliding, without requiring
curvature-sensing or the fine-tuning of basal compliance, and only weakly
influenced by hydrodynamic dissipation and the cell body boundary condition.
High internal dissipation allows the emergence of base-to-tip autonomous
travelling waves, independently of the external fluid viscosity. This enables
progressive swimming in low viscosity environments, and may be critical for
external fertilizers and aquatic microorganisms.


  The nanoscale resolution of super-resolution microscopy has now enabled the
use of fluorescent based molecular localization tools to study whole cell
structural biology. Machine learning based analysis of super-resolution data
offers tremendous potential for discovery of new biology, that by definition is
not known and lacks ground truth. Herein, we describe the application of weakly
supervised learning paradigms to super-resolution microscopy and its potential
to enable the accelerated exploration of the molecular architecture of
subcellular macromolecules and organelles.


  In the modeling of parasite transmission dynamics, understanding the
reproductive characteristics of these parasites is crucial.
  This paper presents a mathematical model that explores the reproductive
behavior of dioecious parasites and its impact on transmission dynamics.
  Specifically, the study focuses on the investigation of various reproductive
variables such as the mating probability and the fertilized egg production in
the case of helminth parasites.
  While previous studies have commonly assumed Poisson and negative binomial
distributions to describe the distribution of parasites among hosts, this study
adopts an arbitrary distribution model and examines its consequences on some
reproductive variables. These variables include mean number of fertile females,
mean egg production, mating probability and mean fertilized egg production.
  In addition, the study of these variables takes into account the sex
distribution of the parasites and whether male and female parasites are
considered to be distributed together or separately.
  We show that the models obtained for the case of male and female parasites
distributed separately in the hosts are ecologically unrealistic.
  We present the results obtained for some specific models and we tested the
models obtained in this work using Monte Carlo simulations.


  In ecology, foraging requires animals to expend energy in order to obtain
resources. The cost of foraging can be reduced through kleptoparasitism, the
theft of a resource that another individual has expended effort to acquire.
Thus, kleptoparasitism is one of the most significant feeding techniques in
ecology. In this study, we investigate a two predator one prey paradigm in
which one predator acts as a kleptoparasite and the other as a host. This
research considers the post-kleptoparasitism scenario, which has received
little attention in the literature. Parametric requirements for the existence
as well as local and global stability of biologically viable equilibria have
been proposed. The occurrences of various one parametric bifurcations, such as
saddle-node bifurcation, transcritical bifurcation, and Hopf bifurcation, as
well as two parametric bifurcations, such as Bautin bifurcation, are explored
in depth. Relatively low growth rate of first predator induces a subcritical
Hopf bifurcation although a supercritical Hopf bifurcation occurs at relatively
high growth rate of first predator making coexistence of all three species
possible. Some numerical simulations have been provided for the purpose of
verifying our theoretical conclusions.


  Large-scale manufacturing of induced pluripotent stem cells (iPSCs) is
essential for cell therapies and regenerative medicines. Yet, iPSCs form large
cell aggregates in suspension bioreactors, resulting in insufficient nutrient
supply and extra metabolic waste build-up for the cells located at the core.
Since subtle changes in micro-environment can lead to a heterogeneous cell
population, a novel Biological System-of-Systems (Bio-SoS) framework is
proposed to model cell-to-cell interactions, spatial and metabolic
heterogeneity, and cell response to micro-environmental variation. Building on
stochastic metabolic reaction network, aggregation kinetics, and
reaction-diffusion mechanisms, the Bio-SoS model characterizes causal
interdependencies at individual cell, aggregate, and cell population levels. It
has a modular design that enables data integration and improves predictions for
different monolayer and aggregate culture processes. In addition, a variance
decomposition analysis is derived to quantify the impact of factors (i.e.,
aggregate size) on cell product health and quality heterogeneity.


  Symmetry concepts in parametrized dynamical systems may reduce the number of
external parameters by a suitable normalization prescription. If, under the
action of a symmetry group G, parameter space A becomes a (locally) trivial
principal bundle, A ~ A/G x G, then the normalized dynamics only depends on the
quotient A/G. In this way, the dynamics of fractional variables in homogeneous
epidemic SI(R)S models, with standard incidence, absence of R-susceptibility
and compartment independent birth and death rates, turns out to be isomorphic
to (a marginally extended version of) Hethcote's classic endemic model, first
presented in 1973. The paper studies a 10-parameter master model with constant
and I-linear vaccination rates, vertical transmission and a vaccination rate
for susceptible newborns. As recently shown by the author, all demographic
parameters are redundant. After adjusting time scale, the remaining 5-parameter
model admits a 3-dimensional abelian scaling symmetry. By normalization we end
up with Hethcote's extended 2-parameter model. Thus, in view of symmetry
concepts, reproving theorems on endemic bifurcation and stability in such
models becomes needless.


  The COVID-19 pandemic has initiated a global health emergency, with an
exigent need for effective cure. Progressively, drug repurposing is emerging as
a promising solution as it saves the time, cost and labor. However, the number
of drug candidates that have been identified as being repurposed for the
treatment of COVID-19 are still insufficient, so more effective and thorough
drug repurposing strategies are required. In this study, we joint the molecular
docking with machine learning regression approaches to find some prospective
therapeutic candidates for COVID-19 treatment. We screened the 5903 approved
drugs for their inhibition by targeting the main protease 3CL of SARS-CoV-2,
which is responsible to replicate the virus. Molecular docking is used to
calculate the binding affinities of these drugs to the main protease 3CL. We
employed several machine learning regression approaches for QSAR modeling to
find out some potential drugs with high binding affinity. Out outcomes
demonstrated that the Decision Tree Regression (DTR) model with best scores of
R2 and RMSE, is the most suitable model for drug repurposing. We shortlisted
six favorable drugs and examined their physiochemical and pharmacokinetic
properties of these top-ranked selected drugs and their best binding
interaction for specific target protease 3CLpro. Our study provides an
efficient framework for drug repurposing against COVID-19, and establishes the
potential of combining molecular docking with machine learning regression
approaches to accelerate the identification of potential therapeutic
candidates. Our findings contribute to the larger goal of finding effective
cures for COVID-19, which is an acute global health challenge.


  Inverse protein folding -- the task of predicting a protein sequence from its
backbone atom coordinates -- has surfaced as an important problem in the "top
down", de novo design of proteins. Contemporary approaches have cast this
problem as a conditional generative modelling problem, where a large generative
model over protein sequences is conditioned on the backbone. While these
generative models very rapidly produce promising sequences, independent draws
from generative models may fail to produce sequences that reliably fold to the
correct backbone. Furthermore, it is challenging to adapt pure generative
approaches to other settings, e.g., when constraints exist. In this paper, we
cast the problem of improving generated inverse folds as an optimization
problem that we solve using recent advances in "deep" or "latent space"
Bayesian optimization. Our approach consistently produces protein sequences
with greatly reduced structural error to the target backbone structure as
measured by TM score and RMSD while using fewer computational resources.
Additionally, we demonstrate other advantages of an optimization-based approach
to the problem, such as the ability to handle constraints.


  Recent advancements in conversational large language models (LLMs), such as
ChatGPT, have demonstrated remarkable promise in various domains, including
drug discovery. However, existing works mainly focus on investigating the
capabilities of conversational LLMs on chemical reaction and retrosynthesis.
While drug editing, a critical task in the drug discovery pipeline, remains
largely unexplored. To bridge this gap, we propose ChatDrug, a framework to
facilitate the systematic investigation of drug editing using LLMs. ChatDrug
jointly leverages a prompt module, a retrieval and domain feedback (ReDF)
module, and a conversation module to streamline effective drug editing. We
empirically show that ChatDrug reaches the best performance on 33 out of 39
drug editing tasks, encompassing small molecules, peptides, and proteins. We
further demonstrate, through 10 case studies, that ChatDrug can successfully
identify the key substructures (e.g., the molecule functional groups, peptide
motifs, and protein structures) for manipulation, generating diverse and valid
suggestions for drug editing. Promisingly, we also show that ChatDrug can offer
insightful explanations from a domain-specific perspective, enhancing
interpretability and enabling informed decision-making. This research sheds
light on the potential of ChatGPT and conversational LLMs for drug editing. It
paves the way for a more efficient and collaborative drug discovery pipeline,
contributing to the advancement of pharmaceutical research and development.


  The effects of climate change, urbanisation and agriculture are changing the
way insects occupy habitats. Some species may utilise anthropogenic
microhabitat features for their existence, either because they prefer them to
natural features, or because of no choice. Other species are dependent on
natural microhabitats. Identifying and analysing these insects' use of natural
and anthropogenic microhabitats is important to assess their responses to a
changing environment, for improving pollination and managing invasive pests.
Traditional studies of insect microhabitat use can now be supplemented by
machine learning-based insect image analysis. Typically, research has focused
on automatic insect classification, but valuable data in image backgrounds has
been ignored. In this research, we analysed the image backgrounds available on
the ALA database to determine their microhabitats. We analysed the
microhabitats of three insect species common across Australia: Drone flies,
European honeybees and European wasps. Image backgrounds were classified as
natural or anthropogenic microhabitats using computer vision and machine
learning tools benchmarked against a manual classification algorithm. We found
flies and honeybees in natural microhabitats, confirming their need for natural
havens within cities. Wasps were commonly seen in anthropogenic microhabitats.
Results show these insects are well adapted to survive in cities. Management of
this invasive pest requires a thoughtful reduction of their access to
human-provided resources. The assessment of insect image backgrounds is
instructive to document the use of microhabitats by insects. The method offers
insight that is increasingly vital for biodiversity management as urbanisation
continues to encroach on natural ecosystems and we must consciously provide
resources within built environments to maintain insect biodiversity and manage
invasive pests.


  Molecule pretraining has quickly become the go-to schema to boost the
performance of AI-based drug discovery. Naturally, molecules can be represented
as 2D topological graphs or 3D geometric point clouds. Although most existing
pertaining methods focus on merely the single modality, recent research has
shown that maximizing the mutual information (MI) between such two modalities
enhances the molecule representation ability. Meanwhile, existing molecule
multi-modal pretraining approaches approximate MI based on the representation
space encoded from the topology and geometry, thus resulting in the loss of
critical structural information of molecules. To address this issue, we propose
MoleculeSDE. MoleculeSDE leverages group symmetric (e.g., SE(3)-equivariant and
reflection-antisymmetric) stochastic differential equation models to generate
the 3D geometries from 2D topologies, and vice versa, directly in the input
space. It not only obtains tighter MI bound but also enables prosperous
downstream tasks than the previous work. By comparing with 17 pretraining
baselines, we empirically verify that MoleculeSDE can learn an expressive
representation with state-of-the-art performance on 26 out of 32 downstream
tasks.


  The need for more usable and explainable machine learning models in
healthcare increases the importance of developing and utilizing causal
discovery algorithms, which aim to discover causal relations by analyzing
observational data. Explainable approaches aid clinicians and biologists in
predicting the prognosis of diseases and suggesting proper treatments. However,
very little research has been conducted at the crossroads between causal
discovery, genomics, and breast cancer, and we aim to bridge this gap.
Moreover, evaluation of causal discovery methods on real data is in general
notoriously difficult because ground-truth causal relations are usually
unknown, and accordingly, in this paper, we also propose to address the
evaluation problem with large language models. In particular, we exploit
suitable causal discovery algorithms to investigate how various perturbations
in the genome can affect the survival of patients diagnosed with breast cancer.
We used three main causal discovery algorithms: PC, Greedy Equivalence Search
(GES), and a Generalized Precision Matrix-based one. We experiment with a
subset of The Cancer Genome Atlas, which contains information about mutations,
copy number variations, protein levels, and gene expressions for 705 breast
cancer patients. Our findings reveal important factors related to the vital
status of patients using causal discovery algorithms. However, the reliability
of these results remains a concern in the medical domain. Accordingly, as
another contribution of the work, the results are validated through language
models trained on biomedical literature, such as BlueBERT and other large
language models trained on medical corpora. Our results profess proper
utilization of causal discovery algorithms and language models for revealing
reliable causal relations for clinical applications.


  Recently, molecular relational learning, whose goal is to predict the
interaction behavior between molecular pairs, got a surge of interest in
molecular sciences due to its wide range of applications. In this work, we
propose CMRL that is robust to the distributional shift in molecular relational
learning by detecting the core substructure that is causally related to
chemical reactions. To do so, we first assume a causal relationship based on
the domain knowledge of molecular sciences and construct a structural causal
model (SCM) that reveals the relationship between variables. Based on the SCM,
we introduce a novel conditional intervention framework whose intervention is
conditioned on the paired molecule. With the conditional intervention
framework, our model successfully learns from the causal substructure and
alleviates the confounding effect of shortcut substructures that are spuriously
correlated to chemical reactions. Extensive experiments on various tasks with
real-world and synthetic datasets demonstrate the superiority of CMRL over
state-of-the-art baseline models. Our code is available at
https://github.com/Namkyeong/CMRL.


  This article presents the "PubChemQC B3LYP/6-31G*//PM6" dataset, containing
electronic properties of 85,938,443 molecules. It includes orbitals, orbital
energies, total energies, dipole moments, and other relevant properties. The
dataset encompasses a wide range of molecules, from essential compounds to
biomolecules up to 1000 molecular weight, covering 94.0% of the original
PubChem Compound catalog (as of August 29, 2016). The electronic properties
were calculated using the B3LYP/6-31G* and PM6 methods. The dataset is
available in three formats: (i) GAMESS quantum chemistry program files, (ii)
selected JSON output files, and (iii) a PostgreSQL database, enabling
researchers to query molecular properties. Five sub-datasets offer more
specific data. The first two subsets include molecules with C, H, O, and N,
under 300 and 500 molecular weight respectively. The third and fourth subsets
contain C, H, N, O, P, S, F, and Cl, under 300 and 500 molecular weight
respectively. The fifth subset includes C, H, N, O, P, S, F, Cl, Na, K, Mg, and
Ca, under 500 molecular weight. Coefficients of determination ranged from 0.892
(CHON500) to 0.803 (whole) for the HOMO-LUMO energy gap. These findings
represent extensive investigations and can be utilized for drug discovery,
material science, and other applications. The datasets are available under the
Creative Commons Attribution 4.0 International license at
https://nakatamaho.riken.jp/pubchemqc.riken.jp/b3lyp_pm6_datasets.html.


  Wind turbines (WT) cause bird and bat mortalities which depend on the WT and
landscape features. The effects of WT features and environmental variables at
different spatial scales associated to bat deaths in a mountainous and forested
area in Thrace, NE Greece were investigated. Initially, we sought to quantify
the most lethal WT characteristic between tower height, rotor diameter and
power. The scale of interaction distance between bat deaths and the land cover
characteristics surrounding the WTs was quantified. A statistical model was
trained and validated against bat deaths and WT, land cover and topography
features. Variance partitioning between bat deaths and the explanatory
covariates was conducted. The trained model was used to predict bat deaths
attributed to existing and future wind farm development in the region. Results
indicated that the optimal interaction distance between WT and surrounding land
cover was 5 km, the larger distance than the ones examined. WT power, natural
land cover type and distance from water explained 40 %, 15 % and 11 %
respectively of the total variance in bat deaths by WTs. The model predicted
that operating but not surveyed WTs comprise of 377.8% and licensed but not
operating yet will contribute to 210.2% additional deaths than the ones
recorded. Results indicate that among all WT features and land cover
characteristics, wind turbine power is the most significant factor associated
to bat deaths. Results indicated that WTs located within 5 km buffer comprised
of natural land cover types have substantial higher deaths. More WT power will
result in more deaths. Wind turbines should not be licensed in areas where
natural land cover at a radius of 5km exceeds 50%. These results are discussed
in the climate-land use-biodiversity-energy nexus.


  Compartmental epidemic models with dynamics that evolve over a graph network
have gained considerable importance in recent years but analysis of
  these models is in general difficult due to their complexity. In this paper,
we develop two positive feedback frameworks that are applicable to the
  study of steady-state values in a wide range of compartmental epidemic
models, including both group and networked
  processes.
  In the case of a group (resp. networked) model, we show that the convergence
limit of the susceptible proportion of the population (resp. the
  susceptible proportion in at least one of the subgroups) is upper bounded by
the reciprocal of the basic reproduction number (BRN) of the model. The
  BRN, when it is greater than unity, thus demonstrates the level of
penetration into a subpopulation by the disease. Both non-strict and strict
  bounds on the convergence limits are derived and shown to correspond to
substantially distinct scenarios in the epidemic processes, one in the
  presence of the endemic state and another without. Formulae for calculating
the limits are provided in the latter case. We apply the developed
  framework to examining various group and networked epidemic models commonly
seen in the literature to verify the validity of our conclusions.


  Diffusion-based manifold learning methods have proven useful in
representation learning and dimensionality reduction of modern high
dimensional, high throughput, noisy datasets. Such datasets are especially
present in fields like biology and physics. While it is thought that these
methods preserve underlying manifold structure of data by learning a proxy for
geodesic distances, no specific theoretical links have been established. Here,
we establish such a link via results in Riemannian geometry explicitly
connecting heat diffusion to manifold distances. In this process, we also
formulate a more general heat kernel based manifold embedding method that we
call heat geodesic embeddings. This novel perspective makes clearer the choices
available in manifold learning and denoising. Results show that our method
outperforms existing state of the art in preserving ground truth manifold
distances, and preserving cluster structure in toy datasets. We also showcase
our method on single cell RNA-sequencing datasets with both continuum and
cluster structure, where our method enables interpolation of withheld
timepoints of data. Finally, we show that parameters of our more general method
can be configured to give results similar to PHATE (a state-of-the-art
diffusion based manifold learning method) as well as SNE (an
attraction/repulsion neighborhood based method that forms the basis of t-SNE).


  Our planet is roughly closed to matter, but open to energy input from the
sun. However, to harness this energy, organisms must transform matter from one
chemical (redox) state to another. For example, photosynthetic organisms can
capture light energy by carrying out a pair of electron donor and acceptor
transformations (e.g., water to oxygen, CO$_2$ to organic carbon). Closure of
ecosystems to matter requires that all such transformations are ultimately
balanced, i.e., other organisms must carry out corresponding reverse
transformations, resulting in cycles that are coupled to each other. A
sustainable closed ecosystem thus requires self-organized cycles of matter, in
which every transformation has sufficient thermodynamic favorability to
maintain an adequate number of organisms carrying out that process. Here, we
propose a new conceptual model that explains the self-organization and emergent
features of closed ecosystems. We study this model with varying levels of
metabolic diversity and energy input, finding that several thermodynamic
features converge across ecosystems. Specifically, irrespective of their
species composition, large and metabolically diverse communities self-organize
to extract roughly 10% of the maximum extractable energy, or 100 fold more than
randomized communities. Moreover, distinct communities implement energy
extraction in convergent ways, as indicated by strongly correlated fluxes
through nutrient cycles. As the driving force from light increases, however,
these features -- fluxes and total energy extraction -- become more variable
across communities, indicating that energy limitation imposes tight
thermodynamic constraints on collective metabolism.


  During the last decades macroecology has identified broad-scale patterns of
abundances and diversity of microbial communities and put forward some
potential explanations for them. However, these advances are not paralleled by
a full understanding of the dynamical processes behind them. In particular,
abundance fluctuations of different species are found to be correlated, both
across time and across communities in metagenomic samples. Reproducing such
correlations through appropriate population models remains an open challenge.
The present paper tackles this problem and points to sparse species
interactions as a necessary mechanism to account for them. Specifically, we
discuss several possibilities to include interactions in population models and
recognize Lotka-Volterra constants as a successful ansatz. For this, we design
a Bayesian inference algorithm to extract sets of interaction constants able to
reproduce empirical probability distributions of pairwise correlations for
diverse biomes. Importantly, the inferred models still reproduce well-known
single-species macroecological patterns concerning abundance fluctuations
across both species and communities. Endorsed by the agreement with the
empirically observed phenomenology, our analyses provide insights on the
properties of the networks of microbial interactions, revealing that sparsity
is a crucial feature.


  DNA sequence alignment is important today as it is usually the first step in
finding gene mutation, evolutionary similarities, protein structure, drug
development and cancer treatment. Covid-19 is one recent example. There are
many sequencing algorithms developed over the past decades but the sequence
alignment using expert systems is quite new. To find DNA sequence alignment,
dynamic programming was used initially. Later faster algorithms used small DNA
sequence length of fixed size to find regions of similarity, and then build the
final alignment using these regions. Such systems were not sensitive but were
fast. To improve the sensitivity, we propose a new algorithm which is based on
finding maximal matches between two sequences, find seeds between them, employ
rules to find more seeds of varying length, and then employ a new stitching
algorithm, and weighted seeds to solve the problem


  Tumor recurrence, driven by the evolution of drug resistance is a major
barrier to therapeutic success in cancer. Resistance is often caused by genetic
alterations such as point mutation, which refers to the modification of a
single genomic base pair, or gene amplification, which refers to the
duplication of a region of DNA that contains a gene. Here we investigate the
dependence of tumor recurrence dynamics on these mechanisms of resistance,
using stochastic multi-type branching process models. We derive tumor
extinction probabilities and deterministic estimates for the tumor recurrence
time, defined as the time when an initially drug sensitive tumor surpasses its
original size after developing resistance. For models of amplification-driven
and mutation-driven resistance, we prove law of large numbers results regarding
the convergence of the stochastic recurrence times to their mean. Additionally,
we prove sufficient and necessary conditions for a tumor to escape extinction
under the gene amplification model, discuss behavior under biologically
relevant parameters, and compare the recurrence time and tumor composition in
the mutation and amplification models both analytically and using simulations.
In comparing these mechanisms, we find that the ratio between recurrence times
driven by amplification vs. mutation depends linearly on the number of
amplification events required to acquire the same degree of resistance as a
mutation event, and we find that the relative frequency of amplification and
mutation events plays a key role in determining the mechanism under which
recurrence is more rapid. In the amplification-driven resistance model, we also
observe that increasing drug concentration leads to a stronger initial
reduction in tumor burden, but that the eventual recurrent tumor population is
less heterogeneous, more aggressive, and harbors higher levels of
drug-resistance.


  We make a retrospective review on various control measures taken by 127
countries/territories during the first wave of COVID-19 pandemic until July 7,
2020, and evaluate their impacts on the epidemic dynamics quantitatively. The
SEIR-QD model, as a representative for general compartment models, is used to
fit the epidemic data, enabling the extraction of crucial model parameters and
dynamical features. The mediation effect of the SEIR-QD model is revealed by
using the mediation analysis with structure equation modeling for multiple
mediators operating in parallel. The inherent impacts of these control policies
on the transmission dynamics of COVID-19 epidemics are clarified, and compared
with results derived from both multiple linear regression and
neural-network-based nonlinear regression. Through this data-driven analysis,
the mediation effect of compartment models is confirmed, which provides a
better understanding on the intrinsic correlations among the strength of
control measures and the dynamical features of COVID-19 epidemics.


  One of the most relevant prognostication tests for tumors is cellular
proliferation, which is most commonly measured by the mitotic activity in
routine tumor sections. The goal of this systematic review is to scholarly
analyze the methods and prognostic relevance of histologically measuring
mitotic activity in canine tumors. A total of 137 articles that correlated the
mitotic activity in canine tumors with patient outcome were identified through
a systematic (PubMed and Scopus) and manual (Google Scholar) literature search
and eligibility screening process. These studies determined the mitotic count
(MC, number of mitotic figures per tumor area) in 126 instances, presumably the
mitotic count (method not specified) in 6 instances and the mitotic index (MI,
proportion of mitotic figures per tumor cells) in 5 instances. A particularly
high risk of bias was identified in the available details of the MC methods and
statistical analysis, which often did not quantify the prognostic
discriminative ability of the MC and only reported p-values. A significant
association of the MC with survival was found in 72/109 (66%) studies. However,
survival was evaluated by at least three studies in only 7 tumor types/groups,
of which a prognostic relevance is apparent for mast cell tumors of the skin,
cutaneous melanoma and soft tissue sarcoma of the skin. None of the studies on
the MI found a prognostic relevance. Further studies on the MC and MI with
standardized methods are needed to prove the prognostic benefit of this test
for further tumor types.


  Macrocyclic peptides are an emerging therapeutic modality, yet computational
approaches for accurately sampling their diverse 3D ensembles remain
challenging due to their conformational diversity and geometric constraints.
Here, we introduce RINGER, a diffusion-based transformer model for
sequence-conditioned generation of macrocycle structures based on internal
coordinates. RINGER provides fast backbone sampling while respecting key
structural invariances of cyclic peptides. Through extensive benchmarking and
analysis against gold-standard conformer ensembles of cyclic peptides generated
with metadynamics, we demonstrate how RINGER generates both high-quality and
diverse geometries at a fraction of the computational cost. Our work lays the
foundation for improved sampling of cyclic geometries and the development of
geometric learning methods for peptides.


  The accurate prediction of changes in protein stability under multiple amino
acid substitutions is essential for realising true in-silico protein re-design.
To this purpose, we propose improvements to state-of-the-art Deep learning (DL)
protein stability prediction models, enabling first-of-a-kind predictions for
variable numbers of amino acid substitutions, on structural representations, by
decoupling the atomic and residue scales of protein representations. This was
achieved using E(3)-equivariant graph neural networks (EGNNs) for both atomic
environment (AE) embedding and residue-level scoring tasks. Our AE embedder was
used to featurise a residue-level graph, then trained to score mutant stability
($\Delta\Delta G$). To achieve effective training of this predictive EGNN we
have leveraged the unprecedented scale of a new high-throughput protein
stability experimental data-set, Mega-scale. Finally, we demonstrate the
immediately promising results of this procedure, discuss the current
shortcomings, and highlight potential future strategies.


  Stochasticity plays a central role in nearly every biological process, and
the noise power spectral density (PSD) is a critical tool for understanding
variability and information processing in living systems. In steady-state, many
such processes can be described by stochastic linear time-invariant (LTI)
systems driven by Gaussian white noise, whose PSD is a complex rational
function of the frequency that can be concisely expressed in terms of their
Jacobian, dispersion, and diffusion matrices, fully defining the statistical
properties of the system's dynamics at steady-state. Here, we arrive at compact
element-wise solutions of the rational function coefficients for the auto- and
cross-spectrum that enable the explicit analytical computation of the PSD in
dimensions n=2,3,4. We further present a recursive Leverrier-Faddeev-type
algorithm for the exact computation of the rational function coefficients.
Crucially, both solutions are free of matrix inverses. We illustrate our
element-wise and recursive solutions by considering the stochastic dynamics of
neural systems models, namely Fitzhugh-Nagumo (n=2), Hindmarsh-Rose (n=3),
Wilson-Cowan (n=4), and the Stabilized Supralinear Network (n=22), as well as
an evolutionary game-theoretic model with mutations (n=5, 31).


  A popular approach to protein design is to combine a generative model with a
discriminative model for conditional sampling. The generative model samples
plausible sequences while the discriminative model guides a search for
sequences with high fitness. Given its broad success in conditional sampling,
classifier-guided diffusion modeling is a promising foundation for protein
design, leading many to develop guided diffusion models for structure with
inverse folding to recover sequences. In this work, we propose diffusioN
Optimized Sampling (NOS), a guidance method for discrete diffusion models that
follows gradients in the hidden states of the denoising network. NOS makes it
possible to perform design directly in sequence space, circumventing
significant limitations of structure-based methods, including scarce data and
challenging inverse design. Moreover, we use NOS to generalize LaMBO, a
Bayesian optimization procedure for sequence design that facilitates multiple
objectives and edit-based constraints. The resulting method, LaMBO-2, enables
discrete diffusions and stronger performance with limited edits through a novel
application of saliency maps. We apply LaMBO-2 to a real-world protein design
task, optimizing antibodies for higher expression yield and binding affinity to
several therapeutic targets under locality and developability constraints,
attaining a 99% expression rate and 40% binding rate in exploratory in vitro
experiments.


  Motivated by the question of the impact of selective advantage in populations
with skewed reproduction mechanims, we study a Moran model with selection. We
assume that there are two types of individuals, where the reproductive success
of one type is larger than the other. The higher reproductive success may stem
from either more frequent reproduction, or from larger numbers of offspring,
and is encoded in a measure $\Lambda$ for each of the two types. Our approach
consists of constructing a $\Lambda$-asymmetric Moran model in which
individuals of the two populations compete, rather than considering a Moran
model for each population. Under certain conditions, that we call the "partial
order of adaptation", we can couple these measures. This allows us to construct
the central object of this paper, the $\Lambda-$asymmetric ancestral selection
graph, leading to a pathwise duality of the forward in time
$\Lambda$-asymmetric Moran model with its ancestral process. Interestingly, the
construction also provides a connection to the theory of optimal transport. We
apply the ancestral selection graph in order to obtain scaling limits of the
forward and backward processes, and note that the frequency process converges
to the solution of an SDE with discontinous paths. Finally, we derive a
Griffiths representation for the generator of the SDE and use it to find a
semi-explicit formula for the probability of fixation of the less beneficial of
the two types.


  Muller's ratchet, in its prototype version, models a haploid, asexual
population whose size~$N$ is constant over the generations. Slightly
deleterious mutations are acquired along the lineages at a constant rate, and
individuals carrying less mutations have a selective advantage. The classical
variant considers {\it fitness proportional} selection, but other fitness
schemes are conceivable as well. Inspired by the work of Etheridge et al.
([EPW09]) we propose a parameter scaling which fits well to the
``near-critical'' regime that was in the focus of [EPW09] (and in which the
mutation-selection ratio diverges logarithmically as $N\to \infty$). Using a
Moran model, we investigate the``rule of thumb'' given in [EPW09] for the click
rate of the ``classical ratchet'' by putting it into the context of new results
on the long-time evolution of the size of the best class of the ratchet with
(binary) tournament selection, which (other than that of the classical ratchet)
follows an autonomous dynamics up to the time of its extinction.
  In [GSW23] it was discovered that the tournament ratchet has a hierarchy of
dual processes which can be constructed on top of an Ancestral Selection graph
with a Poisson decoration. For a regime in which the mutation/selection-ratio
remains bounded away from 1, this was used in [GSW23] to reveal the asymptotics
of the click rates as well as that of the type frequency profile between
clicks. We will describe how these ideas can be extended to the near-critical
regime in which the mutation-selection ratio of the tournament ratchet
converges to 1 as $N\to \infty$.


  Gene expression patterns are established by cross-regulating target genes
that interpret morphogen gradients. However, as development progresses,
morphogen activity is reduced, leaving the emergent pattern without stabilizing
positional cues. The pattern then can be deteriorated by the intrinsically
noisy biochemical processes acting at the cellular level. But remarkably, the
established gene expression patterns remain spatially and temporally stable in
many biological systems. Here we combine spatial-stochastic simulations with an
enhanced sampling method and a recently developed stability theory to address
how spatiotemporal integrity of a gene expression pattern is maintained in
developing tissue lacking morphogen gradients. Using a minimal embryo model
consisting of spatially coupled biochemical reactor volumes, we study a stripe
pattern in which weak cross-repression between nearest neighbor domians
alternates with strong repression between next-nearest neighbor domains,
inspired by the gap gene system in the Drosophila embryo. We find that
fine-tuning of the weak repressive interactions to an optimal level can
significantly increase temporal stability of the expression patterns, allowing
for stable patterns over developmentally relevant times in the absence of
morphogen gradients. The numerically determined optimal parameters closely
agree with the predictions of the stability theory. By analizing the dynamics
of factors characterizing pattern integrity, we trace back the stability
enhancement to the emergence of restoring forces, maintaining the pattern in a
meta-stable basin. Altogether our results demonstrate that metastable
attractors can emerge as a property of stochastic gene expression patterns even
without system-wide positional cues, provided that the gene regulatory
interactions shaping the pattern are optimally tuned.


  We propose a novel approach for predicting protein-peptide interactions using
a bi-modal transformer architecture that learns an inter-facial joint
distribution of residual contacts. The current data sets for crystallized
protein-peptide complexes are limited, making it difficult to accurately
predict interactions between proteins and peptides. To address this issue, we
propose augmenting the existing data from PepBDB with pseudo protein-peptide
complexes derived from the PDB. The augmented data set acts as a method to
transfer physics-based contextdependent intra-residue (within a domain)
interactions to the inter-residual (between) domains. We show that the
distributions of inter-facial residue-residue interactions share overlap with
inter residue-residue interactions, enough to increase predictive power of our
bi-modal transformer architecture. In addition, this dataaugmentation allows us
to leverage the vast amount of protein-only data available in the PDB to train
neural networks, in contrast to template-based modeling that acts as a prior


  Logarithmic oscillations superimposed on a power-law trend appear in the
behavior of various complex hierarchical systems. In this paper, we study the
logarithmic oscillations of relaxation curves in p-adic diffusion models that
are used to describe the conformational dynamics of protein. We consider the
case of a purely p-adic diffusion, as well as the case of p-adic diffusion with
a reaction sink. We show that, relaxation curves for large times in these two
cases are described by a power law on which logarithmic oscillations are
superimposed whose period and amplitude are determined by the parameters of the
model. We also provide a physical explanation of the emergence of oscillations
in relaxation curves and discuss the relation of the results to the experiments
on relaxation dynamics of protein.


  Machine learning-based design has gained traction in the sciences, most
notably in the design of small molecules, materials, and proteins, with
societal implications spanning drug development and manufacturing, plastic
degradation, and carbon sequestration. When designing objects to achieve novel
property values with machine learning, one faces a fundamental challenge: how
to push past the frontier of current knowledge, distilled from the training
data into the model, in a manner that rationally controls the risk of failure.
If one trusts learned models too much in extrapolation, one is likely to design
rubbish. In contrast, if one does not extrapolate, one cannot find novelty.
Herein, we ponder how one might strike a useful balance between these two
extremes. We focus in particular on designing proteins with novel property
values, although much of our discussion addresses machine learning-based design
more broadly.


  Digital contact tracing/notification was initially hailed as a promising
strategy to combat SARS-CoV-2, but in most jurisdictions it did not live up to
its promise. To avert a given transmission event, both parties must have
adopted the tech, it must detect the contact, the primary case must be promptly
diagnosed, notifications must be triggered, and the secondary case must change
their behavior to avoid the focal tertiary transmission event. If we
approximate these as independent events, achieving a 26% reduction in R(t)
would require 80% success rates at each of these six points of failure. Here we
review the six failure rates experienced by a variety of digital contact
tracing/notification schemes, including Singapore's TraceTogether, India's
Aarogya Setu, and leading implementations of the Google Apple Exposure
Notification system. This leads to a number of recommendations, e.g. that
tracing/notification apps be multi-functional and integrated with testing,
manual contact tracing, and the gathering of critical scientific data, and that
the narrative be framed in terms of user autonomy rather than user privacy.


  We developed pgMAP, an analysis pipeline to map gRNA sequencing reads from
dual-targeting CRISPR screens. pgMAP output includes a dual gRNA read counts
table and quality control metrics including the proportion of correctly-paired
reads and CRISPR library sequencing coverage across all time points and
samples. pgMAP is implemented using Snakemake and is available open-source
under the MIT license at https://github.com/fredhutch/pgmap_pipeline.


  Antibodies are Y-shaped proteins that neutralize pathogens and constitute the
core of our adaptive immune system. De novo generation of new antibodies that
target specific antigens holds the key to accelerating vaccine discovery.
However, this co-design of the amino acid sequence and the 3D structure
subsumes and accentuates some central challenges from multiple tasks, including
protein folding (sequence to structure), inverse folding (structure to
sequence), and docking (binding). We strive to surmount these challenges with a
new generative model AbODE that extends graph PDEs to accommodate both
contextual information and external interactions. Unlike existing approaches,
AbODE uses a single round of full-shot decoding and elicits continuous
differential attention that encapsulates and evolves with latent interactions
within the antibody as well as those involving the antigen. We unravel
fundamental connections between AbODE and temporal networks as well as
graph-matching networks. The proposed model significantly outperforms existing
methods on standard metrics across benchmarks.


  The investigation of adaptive dynamics, involving many degrees of freedom on
two separated timescales, one for fast changes of state variables and another
for the slow adaptation of parameters controlling the former's dynamics is
crucial for understanding feedback mechanisms underlying evolutionary and
developmental processes. We present a general statistical-physics theory for
such dynamical systems, thus contributing to the development of an exact
analytical treatment of the nonequilibrium phase transitions in large-scale
systems with adaptation. As an illustration, we apply our framework to
biological systems whose fate is determined by the feedback between genotype
and phenotype: fast phenotypic evolution is shaped by stochastic
gene-expression dynamics and is coupled to slow adaptive changes of genotypes
encoding the network of gene regulations. We establish how a robust phase,
where both genotype and phenotype achieve optimal values, can emerge within an
intermediate level of external noise for high reciprocity of network
interactions.


  Many processes in biology and drug discovery involve various 3D interactions
between molecules, such as protein and protein, protein and small molecule,
etc. Given that different molecules are usually represented in different
granularity, existing methods usually encode each type of molecules
independently with different models, leaving it defective to learn the
universal underlying interaction physics. In this paper, we first propose to
universally represent an arbitrary 3D complex as a geometric graph of sets,
shedding light on encoding all types of molecules with one model. We then
propose a Generalist Equivariant Transformer (GET) to effectively capture both
domain-specific hierarchies and domain-agnostic interaction physics. To be
specific, GET consists of a bilevel attention module, a feed-forward module and
a layer normalization module, where each module is E(3) equivariant and
specialized for handling sets of variable sizes. Notably, in contrast to
conventional pooling-based hierarchical models, our GET is able to retain
fine-grained information of all levels. Extensive experiments on the
interactions between proteins, small molecules and RNA/DNAs verify the
effectiveness and generalization capability of our proposed method across
different domains.


  The Kauffman model is the archetypal model of genetic computation. It
highlights the importance of criticality, at which many biological systems seem
poised. In a series of advances, researchers have honed in on how the number of
attractors in the critical regime grows with network size. But a definitive
answer has proved elusive. We prove that, for the critical Kauffman model with
connectivity one, the number of attractors grows at least, and at most, as
$(2/\!\sqrt{e})^N$. This is the first proof that the number of attractors in a
critical Kauffman model grows exponentially.


  This paper presents a portable, privacy-preserving, in-browser platform for
the reproducible assessment of mutational signature detection methods from
sparse sequencing data generated by targeted gene panels. The platform aims to
address the reproducibility challenges in mutational signature research by
adhering to the FAIR principles, making it findable, accessible, interoperable,
and reusable. Our approach focuses on the detection of specific mutational
signatures, such as SBS3, which have been linked to specific mutagenic
processes. The platform relies on publicly available data, simulation,
downsampling techniques, and machine learning algorithms to generate training
data and labels and to train and evaluate models. The key achievement of our
platform is its transparency, reusability, and privacy preservation, enabling
researchers and clinicians to analyze mutational signatures with the guarantee
that no data circulates outside the client machine.


  Deep generative models have recently emerged as a promising de novo drug
design method. In this respect, deep generative conditional variational
autoencoder (CVAE) models are a powerful approach for generating novel
molecules with desired drug-like properties. However, molecular graph-based
models with disentanglement and multivariate explicit latent conditioning have
not been fully elucidated. To address this, we proposed a molecular-graph
$\beta$-CVAE model for de novo drug design. Here, we empirically tuned the
value of disentanglement and assessed its ability to generate molecules with
optimised univariate- or-multivariate properties. In particular, we optimised
the octanol-water partition coefficient (ClogP), molar refractivity (CMR),
quantitative estimate of drug-likeness (QED), and synthetic accessibility score
(SAS). Results suggest that a lower $\beta$ value increases the uniqueness of
generated molecules (exploration). Univariate optimisation results showed our
model generated molecular property averages of ClogP = 41.07% $\pm$ 0.01% and
CMR 66.76% $\pm$ 0.01% by the Ghose filter. Multivariate property optimisation
results showed that our model generated an average of 30.07% $\pm$ 0.01%
molecules for both desired properties. Furthermore, our model improved the QED
and SAS (exploitation) of molecules generated. Together, these results suggest
that the $\beta$-CVAE could balance exploration and exploitation through
disentanglement and is a promising model for de novo drug design, thus
providing a basis for future studies.


  In the study of cancer evolution and therapeutic strategies, scientific
evidence shows that a key dynamics lies in the tumor-environment interaction.
In particular, oxygen concentration plays a central role in the determination
of the phenotypic heterogeneity of cancer cell populations, whose qualitative
and geometric characteristics are predominant factors in the occurrence of
relapses and failure of eradication. We propose a mathematical model able to
describe the eco-evolutionary spatial dynamics of tumour cells in their
adaptation to hypoxic microenvironments. As a main novelty with respect to the
existing literature, we combine a phenotypic indicator reflecting the
experimentally-observed metabolic trade-off between the hypoxia-resistance
ability and the proliferative potential with a 2d geometric domain, without the
constraint of radial symmetry. The model is settled in the mathematical
framework of phenotype-structured population dynamics and it is formulated in
terms of systems of coupled non-linear integro-differential equations. The
computational outcomes demonstrate that hypoxia-induced selection results in a
geometric characterization of phenotypic-defined tumour niches that impact on
tumour aggressiveness and invasive ability. Furthermore, results show how the
knowledge of environmental characteristics provides a predictive advantage on
tumour mass development in terms of size, shape, and composition.


  The field of protein folding research has been greatly advanced by deep
learning methods, with AlphaFold2 (AF2) demonstrating exceptional performance
and atomic-level precision. As co-evolution is integral to protein structure
prediction, AF2's accuracy is significantly influenced by the depth of multiple
sequence alignment (MSA), which requires extensive exploration of a large
protein database for similar sequences. However, not all protein sequences
possess abundant homologous families, and consequently, AF2's performance can
degrade on such queries, at times failing to produce meaningful results. To
address this, we introduce a novel generative language model, MSA-Augmenter,
which leverages protein-specific attention mechanisms and large-scale MSAs to
generate useful, novel protein sequences not currently found in databases.
These sequences supplement shallow MSAs, enhancing the accuracy of structural
property predictions. Our experiments on CASP14 demonstrate that MSA-Augmenter
can generate de novo sequences that retain co-evolutionary information from
inferior MSAs, thereby improving protein structure prediction quality on top of
strong AF2.


  The Ebola Virus Disease (EVD) can persist in some body fluids after clinical
recovery. In Guinea and the Democratic Republic of the Congo, there are
well-documented cases of EVD reemergence that were associated with previous
outbreaks. In many cases, male EVD survivors were associated with the
re-introduction of new outbreaks. This has shown that even after controlling an
EVD outbreak, a close biomedical monitoring of survivors and contacts is
critical to avoiding future outbreaks. Thus, in order to explore the main
features of EVD transmission dynamics in the context of re-emergence, we used a
compartmental model by considering vaccination around EVD contacts. Analytical
and numerical analyses of the model were conducted. The model is mathematically
and epidemiologically well-posed. We computed the reproductive number (R0) and
the disease equilibrium points (disease-free equilibrium and endemic
equilibrium) for the re-emerging outbreak. The stability analysis of the model
around those equilibrium points was performed. The model undergoes a backward
bifurcation at R0 close to 1, regardless R0 < 1, the disease will not be
eradicated. This means that R0 cannot be considered as an intervention control
measure in our model.


  Inferring microbial interaction networks from abundance patterns is an
important approach to advance our understanding of microbial communities in
general and the human microbiome in particular. Here we suggest discriminating
two levels of information contained in microbial abundance data: (1) the
quantitative abundance values and (2) the pattern of presences and absences of
microbial organisms. The latter allows for a binary view on microbiome data and
a novel interpretation of microbial data as attractors, or more precisely as
fixed points, of a Boolean network.
  Starting from these attractors, our aim is to infer an interaction network
between the species present in the microbiome samples. To accomplish this task,
we introduce a novel inference method that combines the previously published
ESABO (Entropy Shifts of Abundance vectors under Boolean Operations) method
with an evolutionary algorithm. The key idea of our approach is that the
inferred network should reproduce the original set of (observed) binary
abundance patterns as attractors.
  We study the accuracy and runtime properties of this evolutionary method, as
well as its behavior under incomplete knowledge of the attractor sets. Based on
this theoretical understanding of the method we then show an application to
empirical data.


  As the importance of high-throughput screening (HTS) continues to grow due to
its value in early stage drug discovery and data generation for training
machine learning models, there is a growing need for robust methods for
pre-screening compounds to identify and prevent false-positive hits. Small,
colloidally aggregating molecules are one of the primary sources of
false-positive hits in high-throughput screens, making them an ideal candidate
to target for removal from libraries using predictive pre-screening tools.
However, a lack of understanding of the causes of molecular aggregation
introduces difficulty in the development of predictive tools for detecting
aggregating molecules. Herein, we present an examination of the molecular
features differentiating datasets of aggregating and non-aggregating molecules,
as well as a machine learning approach to predicting molecular aggregation. Our
method uses explainable graph neural networks and counterfactuals to reliably
predict and explain aggregation, giving additional insights and design rules
for future screening. The integration of this method in HTS approaches will
help combat false positives, providing better lead molecules more rapidly and
thus accelerating drug discovery cycles.


  Geographic barriers prevent migration between populations, thereby
facilitating speciation through allopatry. However, these barriers can exhibit
dynamic behavior in nature, promoting cycles of expansion and contraction of
populations. Such oscillations cause temporal variations in migration that do
not necessarily prevent speciation; on the contrary, they have been suggested
as a driving force for diversification. Here we present a study on a two-island
neutral speciation model in scenarios with intermittent migration driven by
sea-level fluctuations. Mating is constrained to genetically compatible
individuals inhabiting the same island, and offspring inherit nuclear genomes
from both parents with recombination. We observe pulses of speciation that
would not occur in strict isolation or continuous migration. According to the
seabed height, which modulates the duration of the isolation and connection
periods, the maximum richness occurs at different times and in an ephemeral
fashion. The expansion-contraction dynamics can accelerate diversification, but
a long time in isolation can reduce the richness to one species per island,
resembling patterns described by the taxon pulse hypothesis of diversification.
Together with other studies, our results support the relevance of research on
the impact of variable migration on diversification, suggested to be related to
regions of high diversity.


  We propose a novel alternative approach to our previous work (Ben Hammouda et
al., 2023) to improve the efficiency of Monte Carlo (MC) estimators for rare
event probabilities for stochastic reaction networks (SRNs). In the same spirit
of (Ben Hammouda et al., 2023), an efficient path-dependent measure change is
derived based on a connection between determining optimal importance sampling
(IS) parameters within a class of probability measures and a stochastic optimal
control formulation, corresponding to solving a variance minimization problem.
In this work, we propose a novel approach to address the encountered curse of
dimensionality by mapping the problem to a significantly lower-dimensional
space via a Markovian projection (MP) idea. The output of this model reduction
technique is a low-dimensional SRN (potentially even one dimensional) that
preserves the marginal distribution of the original high-dimensional SRN
system. The dynamics of the projected process are obtained by solving a related
optimization problem via a discrete $L^2$ regression. By solving the resulting
projected Hamilton-Jacobi-Bellman (HJB) equations for the reduced-dimensional
SRN, we obtain projected IS parameters, which are then mapped back to the
original full-dimensional SRN system, resulting in an efficient IS-MC estimator
for rare events probabilities of the full-dimensional SRN. Our analysis and
numerical experiments reveal that the proposed MP-HJB-IS approach substantially
reduces the MC estimator variance, resulting in a lower computational
complexity in the rare event regime than standard MC estimators.


  Chlamydiae are bacteria with an interesting unusual developmental cycle. A
single bacterium in its infectious form (elementary body, EB) enters the host
cell, where it converts into its dividing form (reticulate body, RB), and
divides by binary fission. Since only the EB form is infectious, before the
host cell dies, RBs start to convert into EBs. After the host cell dies RBs do
not survive. We model the population growth by a 2-type discrete-time branching
process, where the probability of duplication depends on the state. Maximizing
the EB production leads to a stochastic optimization problem. Simulation study
shows that our novel model is able to reproduce the main features of the
development of the population.


  The dynamic nature of proteins is crucial for determining their biological
functions and properties, and molecular dynamics (MD) simulations stand as a
predominant tool to study such phenomena. By utilizing empirically derived
force fields, MD simulations explore the conformational space through
numerically evolving the system along MD trajectories. However, the high-energy
barrier of the force fields can hamper the exploration of MD, resulting in
inadequately sampled ensemble. In this paper, we propose leveraging score-based
generative models (SGMs) trained on general protein structures to perform
protein conformational sampling to complement traditional MD simulations. We
argue that SGMs can provide a novel framework as an alternative to traditional
enhanced sampling methods by learning multi-level score functions, which
directly sample a diversity-controllable ensemble of conformations. We
demonstrate the effectiveness of our approach on several benchmark systems by
comparing the results with long MD trajectories and state-of-the-art generative
structure prediction models. Our framework provides new insights that SGMs have
the potential to serve as an efficient and simulation-free methods to study
protein dynamics.


  Despite recent advances in the length and the accuracy of long-read data,
building haplotype-resolved genome assemblies from telomere to telomere still
requires considerable computational resources. In this study, we present an
efficient de novo assembly algorithm that combines multiple sequencing
technologies to scale up population-wide telomere-to-telomere assemblies. By
utilizing twenty-two human and two plant genomes, we demonstrate that our
algorithm is around an order of magnitude cheaper than existing methods, while
producing better diploid and haploid assemblies. Notably, our algorithm is the
only feasible solution to the haplotype-resolved assembly of polyploid genomes.


  We consider the Moran model of population genetics with two types, mutation,
and selection, and investigate the line of descent of a randomly-sampled
individual from a contemporary population. We trace this ancestral line back
into the distant past, far beyond the most recent common ancestor of the
population (thus connecting population genetics to phylogeny) and analyse the
mutation process along this line.
  To this end, we use the pruned lookdown ancestral selection graph, which
consists of the set of potential ancestors of the sampled individual at any
given time. A crucial observation is that a mutation on the ancestral line
requires that the ancestral line occupy the top position in the graph just
`before' the event (in forward time). Relative to the neutral case (that is,
without selection), we obtain a general bias towards the beneficial type, an
increase in the beneficial mutation rate, and a decrease in the deleterious
mutation rate. This sheds new light on previous analytical results. We discuss
our findings in the light of a well-known observation at the interface of
phylogeny and population genetics, namely, the difference in the mutation rates
(or, more precisely, mutation fluxes) estimated via phylogenetic methods
relative to those observed in pedigree studies.


  Mathematical models of genetic evolution often come in pairs, connected by a
so-called duality relation. The most seminal example are the Wright-Fisher
diffusion and the Kingman coalescent, where the former describes the stochastic
evolution of neutral allele frequencies in a large population forwards in time,
and the latter describes the genetic ancestry of randomly sampled individuals
from the population backwards in time. As well as providing a richer
description than either model in isolation, duality often yields equations
satisfied by quantities of interest. We employ the so-called Bernoulli factory
- a celebrated tool in simulation-based computing - to derive duality relations
for broad classes of genetics models. As concrete examples, we present
Wright-Fisher diffusions with general drift functions, and Allen-Cahn equations
with general, nonlinear forcing terms. The drift and forcing functions can be
interpreted as the action of frequency-dependent selection. To our knowledge,
this work is the first time a connection has been drawn between Bernoulli
factories and duality in models of population genetics.


  Photosynthesis is a plausible pathway for the sustenance of a substantial
biosphere on an exoplanet. In fact, it is also anticipated to create
distinctive biosignatures detectable by next-generation telescopes. In this
work, we explore the excitation features of photopigments that harvest
electromagnetic radiation by constructing a simple quantum-mechanical model.
Our analysis suggests that the primary Earth-based photopigments for
photosynthesis may not function efficiently at wavelengths $> 1.1$ $\mu$m. In
the context of (hypothetical) extrasolar photopigments, we calculate the
potential number of conjugated $\pi$-electrons ($N_\star$) in the relevant
molecules, which can participate in the absorption of photons. By hypothesizing
that the absorption maxima of photopigments are close to the peak spectral
photon flux of the host star, we utilize the model to estimate $N_\star$. As
per our formalism, $N_\star$ is modulated by the stellar temperature, and is
conceivably higher (lower) for planets orbiting stars cooler (hotter) than the
sun; exoplanets around late-type M-dwarfs might require an $N_\star$ twice that
of the Earth. We conclude the analysis with a brief exposition of how our model
could be empirically tested by future observations.


  Phase-separated biomolecular condensates exhibit a wide range of dynamical
properties, which depend on the sequences of the constituent proteins and RNAs.
However, it is unclear to what extent condensate dynamics can be tuned without
also changing the thermodynamic properties that govern phase separation. Using
coarse-grained simulations of intrinsically disordered proteins, we show that
the dynamics and thermodynamics of homopolymer condensates are strongly
correlated, with increased condensate stability being coincident with low
mobilities and high viscosities. We then apply an "active learning" strategy to
identify heteropolymer sequences that break this correlation. This data-driven
approach and accompanying analysis reveal how heterogeneous amino-acid
compositions and non-uniform sequence patterning map to a range of
independently tunable dynamical and thermodynamic properties of biomolecular
condensates. Our results highlight key molecular determinants governing the
physical properties of biomolecular condensates and establish design rules for
the development of stimuli-responsive biomaterials.


  A novel generalized mean field approximation, called the Small-Coupling
Dynamic Cavity (SCDC) method, for epidemic inference and risk assessment is
presented. The method is developed within a fully Bayesian framework and
accounts for non-causal effects generated by the presence of observations. It
is based on a graphical model representation of the epidemic stochastic process
and utilizes dynamic cavity equations to derive a set of self-consistent
equations for probability marginals defined on the edges of the contact graph.
By performing a small-coupling expansion, a pair of time-dependent cavity
messages is obtained, which capture the probability of individual infection and
the conditioning power of observations. In its efficient formulation, the
computational cost per iteration of the SCDC algorithm is linear in the
duration of the epidemic dynamics and in the number of contacts. The SCDC
method is derived for the Susceptible-Infected (SI) model and straightforwardly
applicable to other Markovian epidemic processes, including recurrent ones. It
exhibits high accuracy in assessing individual risk on par with Belief
Propagation techniques and outperforming heuristic methods based on
individual-based mean-field approximations. Although convergence issues may
arise due to long-range correlations in contact graphs, the estimated marginal
probabilities remain sufficiently accurate for reliable risk estimation. Future
work includes extending the method to non-Markovian recurrent epidemic models
and investigating the role of second-order terms in the small coupling
expansion of the observation-reweighted Dynamic Cavity equations.


  While studying nucleic acids to reveal the weak interactions responsible for
their three-dimensional structure and for their interactions with drugs, we
also contributed to the field of biomolecular mass spectrometry, both in terms
of fundamental understanding and with new methodological developments. A first
goal was to develop mass spectrometry approaches to detect non-covalent
interactions between antitumor drugs and their DNA target. Twenty years ago,
our attention turned towards specific DNA structures such as the G-quadruplex
(a structure formed by guanine-rich strands). Mass spectrometry allows to
discern which molecules interact with one another by measuring the masses of
the complexes, and quantity the affinities by measuring their abundance. The
most important findings came from unexpected masses. For example, we showed the
formation of higher- or lower-order structures by G-quadruplexes used in
traditional biophysical assays. We can also derive complete thermodynamic and
kinetic description of G-quadruplex folding pathways by measuring cation
binding, one at a time. Getting quantitative information required accounting
for nonspecific adduct formation and for the response factors of the different
molecular forms. But with these caveats in mind, the approach is mature enough
for routine biophysical characterization of nucleic acids. A second goal is to
obtain more detailed structural information on each of the complexes separated
by the mass spectrometer. One such approaches is ion mobility spectrometry, and
even today the challenge lies in the structural interpretation of the
measurements. We showed that, although structures such as G-quadruplexes are
well-preserved in the MS conditions, double helices actually get more compact
in the gas phase. These major rearrangements forced us to challenge comfortable
assumptions. Further work is still needed to generalize how to deduce
structures in solution from ion mobility spectrometry data, and in particular
how to account for the electrospray charging mechanisms and for ion internal
energy effects. These studies also called for complementary approaches to ion
mobility. Recently, we applied isotope exchange labeling mass spectrometry to
characterize nucleic acid structures for the first time, and we reported the
first ever circular dichroism ion spectroscopy measurement on mass-selected
trapped ions. Circular dichroism plays a key role in assigning the stacking
topology, and our new method now opens the door to characterizing a wide
variety of chiral molecule by mass spectrometry. In summary, advanced mass
spectrometry approaches to characterize gas-phase structures work well for
G-quadruplexes because they are stiffened by inner cations. The next objective
will be to generalize these methodologies to a wider range of nucleic acid
structures.


  Human telomeric G-quadruplex DNA structures are attractive anticancer drug
targets, but the target's polymorphism complicates the drug design: different
ligands prefer different folds, and very few complexes have been solved at high
resolution. Here we report that Phen-DC3, one of the most prominent
G-quadruplex ligands in terms of high binding affinity and selectivity, causes
dTAGGG(TTAGGG)3 to completely change its fold in KCl solution from a hybrid-1
to an antiparallel chair-type structure, wherein the ligand intercalates
between a two-quartet unit and a pseudo-quartet, thereby ejecting one potassium
ion. This unprecedented high-resolution NMR structure shows for the first time
a true ligand intercalation into an intramolecular G-quadruplex.


  Adopting a $300 \, \mu$s-long molecular dynamics (MD) trajectory of the
reversible folding of villin headpiece (HP35) published by D. E. Shaw Research,
we recently constructed a Markov state model (MSM) of the folding process based
on interresidue contacts [J. Chem. Theory Comput. 2023, ${\bf {19}}$, 3391].
The model reproduces the MD folding times of the system and predicts that both
the native basin and the unfolded region of the free energy landscape are
partitioned into several metastable substates that are structurally well
characterized. Recognizing the need to establish well-defined but nontrivial
benchmark problems, in this Perspective we study to what extent and in what
sense this MSM may be employed as a reference model. To this end, we test the
robustness of the MSM by comparing it to models that use alternative
combinations of features, dimensionality reduction methods and clustering
schemes. The study suggests some main characteristics of the folding of HP35,
which should be reproduced by any other competitive model of the system.
Moreover, the discussion reveals which parts of the MSM workflow matter most
for the considered problem, and illustrates the promises and possible pitfalls
of state-based models for the interpretation of biomolecular simulations.


  The control of opportunistic infections among HIV infected individuals should
be one of the major public health concerns in reducing mortality rate of
individuals living with HIV/AIDS. In this study a deterministic co-infection
mathematical model is employed to provide a quantification of treatment at each
contagious stage against Pneumocystis Pneumonia (PCP) among HIV infected
individuals on ART. The disease-free equilibrium for the HIV/AIDS sub model,
PCP sub model and the co-infection model are shown to be locally asymptotically
stable when their associated disease threshold parameter is less than a unity.
By use of suitable Lyapunov functions, the endemic equilibrium corresponding to
HIV/AIDS and PCP sub models are globally asymptotically stable whenever
$\mathcal{R}_{0H}>1$ and $\mathcal{R}_{0P}>1$ respectively. The sensitivity
analysis results implicate that the effective contact rates are the main
mechanisms fueling the proliferation of the two diseases and on the other hand
treatment efforts play an important role in reducing the incidence. Numerical
simulations show that treatment of PCP at all contagious stages reduces its
burden on HIV/AIDS patients and dual treatment of the co-infected individuals
significantly reduces the burden of the co-infection.


  In recent years, in-silico molecular design has received much attention from
the machine learning community. When designing a new compound for
pharmaceutical applications, there are usually multiple properties of such
molecules that need to be optimised: binding energy to the target,
synthesizability, toxicity, EC50, and so on. While previous approaches have
employed a scalarization scheme to turn the multi-objective problem into a
preference-conditioned single objective, it has been established that this kind
of reduction may produce solutions that tend to slide towards the extreme
points of the objective space when presented with a problem that exhibits a
concave Pareto front. In this work we experiment with an alternative
formulation of goal-conditioned molecular generation to obtain a more
controllable conditional model that can uniformly explore solutions along the
entire Pareto front.


  Directed evolution is a molecular biology technique that is transforming
protein engineering by creating proteins with desirable properties and
functions. However, it is experimentally impossible to perform the deep
mutational scanning of the entire protein library due to the enormous
mutational space, which scales as $20^N$ , where N is the number of amino
acids. This has led to the rapid growth of AI-assisted directed evolution
(AIDE) or AI-assisted protein engineering (AIPE) as an emerging research field.
Aided with advanced natural language processing (NLP) techniques, including
long short-term memory, autoencoder, and transformer, sequence-based embeddings
have been dominant approaches in AIDE and AIPE. Persistent Laplacians, an
emerging technique in topological data analysis (TDA), have made
structure-based embeddings a superb option in AIDE and AIPE. We argue that a
class of persistent topological Laplacians (PTLs), including persistent
Laplacians, persistent path Laplacians, persistent sheaf Laplacians, persistent
hypergraph Laplacians, persistent hyperdigraph Laplacians, and evolutionary de
Rham-Hodge theory, can effectively overcome the limitations of the current TDA
and offer a new generation of more powerful TDA approaches. In the general
framework of topological deep learning, mathematics-assisted directed evolution
(MADE) has a great potential for future protein engineering.


  Pathogens drive changes in host immune systems that in turn exert pressure
for pathogens to evolve. Quantifying and understanding this constant
coevolutionary process has clear practical global health implications. Yet its
relatively easier accessibility compared to macroevolution makes it a
fascinating system to learn about the basic laws of evolution. Focusing on
immune-viral evolution, we present an overview of theoretical and experimental
approaches that have recently started coming together to build the foundations
for a quantitative and predictive co-evolutionary theory.


  Protein-ligand binding affinity (PLBA) prediction is the fundamental task in
drug discovery. Recently, various deep learning-based models predict binding
affinity by incorporating the three-dimensional structure of protein-ligand
complexes as input and achieving astounding progress. However, due to the
scarcity of high-quality training data, the generalization ability of current
models is still limited. In addition, different bioassays use varying affinity
measurement labels (i.e., IC50, Ki, Kd), and different experimental conditions
inevitably introduce systematic noise, which poses a significant challenge to
constructing high-precision affinity prediction models. To address these
issues, we (1) propose Multi-task Bioassay Pre-training (MBP), a pre-training
framework for structure-based PLBA prediction; (2) construct a pre-training
dataset called ChEMBL-Dock with more than 300k experimentally measured affinity
labels and about 2.8M docked three-dimensional structures. By introducing
multi-task pre-training to treat the prediction of different affinity labels as
different tasks and classifying relative rankings between samples from the same
bioassay, MBP learns robust and transferrable structural knowledge from our new
ChEMBL-Dock dataset with varied and noisy labels. Experiments substantiate the
capability of MBP as a general framework that can improve and be tailored to
mainstream structure-based PLBA prediction tasks. To the best of our knowledge,
MBP is the first affinity pre-training model and shows great potential for
future development.


  Given the increasing volume and quality of genomics data, extracting new
insights requires interpretable machine-learning models. This work presents
Genomic Interpreter: a novel architecture for genomic assay prediction. This
model outperforms the state-of-the-art models for genomic assay prediction
tasks. Our model can identify hierarchical dependencies in genomic sites. This
is achieved through the integration of 1D-Swin, a novel Transformer-based block
designed by us for modelling long-range hierarchical data. Evaluated on a
dataset containing 38,171 DNA segments of 17K base pairs, Genomic Interpreter
demonstrates superior performance in chromatin accessibility and gene
expression prediction and unmasks the underlying `syntax' of gene regulation.


  Using anonymized mobility data from Facebook users and publicly available
information on the Italian population, we model the circulation of people in
Italy before and during the early phase of the SARS-CoV-2 pandemic (COVID-19).
We perform a spatial and temporal clustering of the movement network at the
level of fluxes across provinces on a daily basis. The resulting partition in
time successfully identifies the first two lockdowns without any prior
information. Similarly, the spatial clustering returns 11 to 23 clusters
depending on the period ("standard" mobility vs. lockdown) using the greedy
modularity communities clustering method, and 16 to 30 clusters using the
critical variable selection method. Fascinatingly, the spatial clusters
obtained with both methods are strongly reminiscent of the 11 regions into
which emperor Augustus had divided Italy according to Pliny the Elder. This
work introduces and validates a data analysis pipeline that enables us: i) to
assess the reliability of data obtained from a partial and potentially biased
sample of the population in performing estimates of population mobility
nationwide; ii) to identify areas of a Country with well-defined mobility
patterns, and iii) to distinguish different patterns from one another, resolve
them in time and find their optimal spatial extent. The proposed method is
generic and can be applied to other countries, with different geographical
scales, and also to similar networks (e.g. biological networks). The results
can thus represent a relevant step forward in the development of methods and
strategies for the containment of future epidemic phenomena.


  Periodic boundary conditions are commonly applied in molecular dynamics
simulations in the microcanonical (NVE), canonical (NVT) and
isothermal-isobaric (NpT) ensembles. In their simplest application, a
biological system of interest is placed in the middle of a solvation box, which
is chosen 'sufficiently large' to minimize any numerical artefacts associated
with the periodic boundary conditions. This practical approach brings
limitations to the size of biological systems that can be simulated. Here, we
study simulations of effectively infinitely-long nucleic acids, which are
solvated in the directions perpendicular to the polymer chain, while periodic
boundary conditions are also applied along the polymer chain. We study the
effects of these asymmetric periodic boundary conditions (APBC) on the
simulated results, including the mechanical properties of biopolymers and the
properties of the surrounding solvent. To get some further insights into the
advantages of using the APBC, a coarse-grained worm-like chain model is first
studied, illustrating how the persistence length can be extracted from local
properties of the polymer chain, which are less affected by the APBC than some
global averages. This is followed by all-atom molecular dynamics simulations of
DNA in ionic solutions, where we use the APBC to investigate sequence-dependent
properties of DNA molecules and properties of the surrounding solvent.


  Advances in deep learning have greatly improved structure prediction of
molecules. However, many macroscopic observations that are important for
real-world applications are not functions of a single molecular structure, but
rather determined from the equilibrium distribution of structures. Traditional
methods for obtaining these distributions, such as molecular dynamics
simulation, are computationally expensive and often intractable. In this paper,
we introduce a novel deep learning framework, called Distributional Graphormer
(DiG), in an attempt to predict the equilibrium distribution of molecular
systems. Inspired by the annealing process in thermodynamics, DiG employs deep
neural networks to transform a simple distribution towards the equilibrium
distribution, conditioned on a descriptor of a molecular system, such as a
chemical graph or a protein sequence. This framework enables efficient
generation of diverse conformations and provides estimations of state
densities. We demonstrate the performance of DiG on several molecular tasks,
including protein conformation sampling, ligand structure sampling,
catalyst-adsorbate sampling, and property-guided structure generation. DiG
presents a significant advancement in methodology for statistically
understanding molecular systems, opening up new research opportunities in
molecular science.


  Our earlier work in \cite{nguyen2022population} shows that concentrating the
resources on the upstream end tends to maximize the total biomass in a
metapopulation model for a stream species. In this paper, we continue our
research direction by further considering a Lotka-Voletrra competition patch
model for two stream species. We show that the species whose resource
allocations maximize the total biomass has competitive advantage.


  In a previous paper [CSIAM Trans. Appl. Math. 2 (2021), 1-55], the authors
proposed a theoretical framework for the analysis of RNA velocity, which is a
promising concept in scRNA-seq data analysis to reveal the cell
state-transition dynamical processes underlying snapshot data. The current
paper is devoted to the algorithmic study of some key components in RNA
velocity workflow. Four important points are addressed in this paper: (1) We
construct a rational time-scale fixation method which can determine the global
gene-shared latent time for cells. (2) We present an uncertainty quantification
strategy for the inferred parameters obtained through the EM algorithm. (3) We
establish the optimal criterion for the choice of velocity kernel bandwidth
with respect to the sample size in the downstream analysis and discuss its
implications. (4) We propose a temporal distance estimation approach between
two cell clusters along the cellular development path. Some illustrative
numerical tests are also carried out to verify our analysis. These results are
intended to provide tools and insights in further development of RNA velocity
type methods in the future.


  Phylogenetics is now fundamental in life sciences, providing insights into
the earliest branches of life and the origins and spread of epidemics. However,
finding suitable phylogenies from the vast space of possible trees remains
challenging. To address this problem, for the first time, we perform both tree
exploration and inference in a continuous space where the computation of
gradients is possible. This continuous relaxation allows for major leaps across
tree space in both rooted and unrooted trees, and is less susceptible to
convergence to local minima. Our approach outperforms the current best methods
for inference on unrooted trees and, in simulation, accurately infers the tree
and root in ultrametric cases. The approach is effective in cases of empirical
data with negligible amounts of data, which we demonstrate on the phylogeny of
jawed vertebrates. Indeed, only a few genes with an ultrametric signal were
generally sufficient for resolving the major lineages of vertebrates.
Optimisation is possible via automatic differentiation and our method presents
an effective way forwards for exploring the most difficult, data-deficient
phylogenetic questions.


  The mechanisms underlying lifespan evolution in organisms have long been
mysterious. However, recent studies have demonstrated that organisms
evolutionarily gain noncoding RNAs (ncRNAs) that carry endogenous profound
functions in higher organisms, including lifespan. This study unveils ncRNAs as
crucial drivers driving animal lifespan evolution. Species in the animal
kingdom evolutionarily increase their ncRNA length in their genomes, coinciding
with trimming mitochondrial genome length. This leads to lower energy
consumption and ultimately lifespan extension. Notably, during lifespan
extension, species exhibit a gradual acquisition of long-life ncRNA motifs
while concurrently losing short-life motifs. These longevity-associated ncRNA
motifs, such as GGTGCG, are particularly active in key tissues, including the
endometrium, ovary, testis, and cerebral cortex. The activation of ncRNAs in
the ovary and endometrium offers insights into why women generally exhibit
longer lifespans than men. This groundbreaking discovery reveals the pivotal
role of ncRNAs in driving lifespan evolution and provides a fundamental
foundation for the study of longevity and aging.


  In 2006, Warnow, Evans, Ringe, and Nakhleh proposed a stochastic model
(hereafter, the WERN 2006 model) of multi-state linguistic character evolution
that allowed for homoplasy and borrowing. They proved that if there is no
borrowing between languages and homoplastic states are known in advance, then
the phylogenetic tree of a set of languages is statistically identifiable under
this model, and they presented statistically consistent methods for estimating
these phylogenetic trees. However, they left open the question of whether a
phylogenetic network -- which would explicitly model borrowing between
languages that are in contact -- can be estimated under the model of character
evolution. Here, we establish that under some mild additional constraints on
the WERN 2006 model, the phylogenetic network topology is statistically
identifiable, and we present algorithms to infer the phylogenetic network. We
discuss the ramifications for linguistic phylogenetic network estimation in
practice, and suggest directions for future research.


  We study the effect of the mutant's degree on the fixation probability,
extinction and fixation times, in the Moran process on Erd\"{o}s-R\'{e}nyi and
Barab\'{a}si-Albert graphs. We performed stochastic simulations and use
mean-field type approximations to obtain analytical formulas. We showed that
the initial placement of a mutant has a significant impact on the fixation
probability and extinction time, while it has no effect on the fixation time.
In both types of graphs, an increase in the degree of a initial mutant results
in a decreased fixation probability and a shorter time to extinction.


  Ecologists have long investigated how the demographic and movement parameters
of a population determine its spatial spread and the critical habitat size that
can sustain it. Yet, most existing models make oversimplifying assumptions
about individual movement behavior, neglecting how landscape heterogeneity
influences dispersal. We relax this assumption and introduce a
reaction-advection-diffusion model that describes the spatial density
distribution of a population with space-dependent movement bias toward
preferred regions, including avoidance of degraded habitats. In this scenario,
the critical habitat size depends on the spatial location of the habitat edges
with respect to the preferred regions and on the intensity of the movement bias
components. In particular, we identify parameter regions where the critical
habitat size decreases when diffusion increases, a phenomenon known as the
``drift paradox''. We also find that biased movement toward low-quality or
highly populated regions can reduce the population size, therefore creating
ecological traps. Our results emphasize the importance of species-specific
movement behavior and habitat selection as drivers of population dynamics in
fragmented landscapes and, therefore, the need to account for them in the
design of protected areas.


  It was recently shown that a large class of phylogenetic networks, the
`labellable' networks, is in bijection with the set of `expanding' covers of
finite sets. In this paper, we show how several prominent classes of
phylogenetic networks can be characterised purely in terms of properties of
their associated covers. These classes include the tree-based, tree-child,
orchard, tree-sibling, and normal networks.