  In this paper, we propose and analyze a reaction-diffusion
susceptible-infected-susceptible (SIS) epidemic patch model. The individuals
are assumed to reside in different patches, where they are able to move inside
and among the patches. The movement of individuals inside the patches is
descried by diffusion terms, and the movement pattern among patches is modeled
by an essentially nonnegative matrix. We define a basic reproduction number
$\mathcal{R}_0$ for the model and show that it is a threshold value for disease
extinction versus persistence. The monotone dependence of $\mathcal{R}_0$ on
the movement rates of infected individuals is proved when the dispersal pattern
is symmetric or non-symmetric. Numerical simulations are performed to
illustrate the impact of the movement of individuals inside and among patches
on the transmission of the disease.


  The expression of one gene might be regulated by its corresponding protein,
which is called autoregulation. Although gene regulation is a central topic in
biology, autoregulation is much less studied. In general, it is extremely
difficult to determine the existence of autoregulation with direct biochemical
approaches. Nevertheless, some papers have observed that certain types of
autoregulations are linked to noise levels in gene expression. We generalize
these results by two propositions on discrete-state continuous-time Markov
chains. These two propositions form a simple but robust method to infer the
existence of autoregulation in certain scenarios from gene expression data.
This method only depends on the Fano factor, namely the ratio of variance and
mean of the gene expression level. Compared to other methods for inferring
autoregulation, our method only requires non-interventional one-time data, and
does not need to estimate parameters. Besides, our method has few restrictions
on the model. We apply this method to four groups of experimental data and find
some genes that might have autoregulation. Some inferred autoregulations have
been verified by experiments or other theoretical works.


  We study generalised rock-paper-scissors models with an arbitrary odd number
N \geq 5 of species, among which n are weak, with 2 \leq n \leq (N-1)/2.
Because of the species' weakness, the probability of individuals conquering
territory in the cyclic spatial game is low. Running stochastic simulations, we
study the role of unevenness in the rock-paper-scissors game in spatial
patterns and population dynamics, considering diverse models where the weak
species are in different positions in the cyclic game order. Studying systems
with five and seven species, we discover that the individuals' spatial
organisation arising from the pattern formation process determines the
stability of the cyclic game with multiple weak species. Our outcomes show that
the presence of species unbalances the spatial distribution of organisms of the
same species bringing consequences on territorial dominance, with the
predominant species being determined by the position in the cyclic game order.
Our simulations elucidate that, in general, the further apart the regions
inhabited by different weak species are, the less the coexistence between the
species is jeopardised. We show that if multiple weak species occupy adjacent
spatial domains, the unevenness in the cyclic game is reinforced, maximising
the chances of biodiversity loss. Our discoveries may also be helpful to
biologists in comprehending systems where weak species unbalance biodiversity
stability.


  Understanding the timing of the peak of a disease outbreak forms an important
part of epidemic forecasting. In many cases, such information is essential for
planning increased hospital bed demand and for designing of public health
interventions. The time taken for an outbreak to become large is inherently
stochastic, and therefore uncertain, but after a sufficient number of
infections has been reached the subsequent dynamics can be modelled accurately
using ordinary differential equations. Here, we present analytical and
numerical methods for approximating the time at which a stochastic model of a
disease outbreak reaches a large number of cases and for quantifying the
uncertainty arising from demographic stochasticity around that time. We then
project this uncertainty forwards in time using an ordinary differential
equation model in order to obtain a distribution for the peak timing of the
epidemic that agrees closely with large simulations but that, for error
tolerances relevant to most realistic applications, requires a fraction of the
computational cost of full Monte Carlo approaches.


  In therapeutic design, balancing various physiochemical properties is crucial
for molecule development, similar to how Multiparameter Optimization (MPO)
evaluates multiple variables to meet a primary goal. While many molecular
features can now be predicted using \textit{in silico} methods, aiding early
drug development, the vast data generated from high throughput virtual
screening challenges the practicality of traditional MPO approaches. Addressing
this, we introduce a scaffold focused graph-based Markov chain Monte Carlo
framework (ScaMARS) built to generate molecules with optimal properties. This
innovative framework is capable of self-training and handling a wider array of
properties, sampling different chemical spaces according to the starting
scaffold. The benchmark analysis on several properties shows that ScaMARS has a
diversity score of 84.6\% and has a much higher success rate of 99.5\% compared
to conditional models. The integration of new features into MPO significantly
enhances its adaptability and effectiveness in therapeutic design, facilitating
the discovery of candidates that efficiently optimize multiple properties.


  In growing populations, the fate of mutations depends on their competitive
ability against the ancestor and their ability to colonize new territory. Here
we present a theory that integrates both aspects of mutant fitness by coupling
the classic description of one-dimensional competition (Fisher equation) to the
minimal model of front shape (KPZ equation). We solved these equations and
found three regimes, which are controlled solely by the expansion rates, solely
by the competitive abilities, or by both. Collectively, our results provide a
simple framework to study spatial competition.


  Despite advances in clinical genetic testing, including the introduction of
exome sequencing (ES), more than 50% of individuals with a suspected Mendelian
condition lack a precise molecular diagnosis. Clinical evaluation is
increasingly undertaken by specialists outside of clinical genetics, often
occurring in a tiered fashion and typically ending after ES. The current
diagnostic rate reflects multiple factors, including technical limitations,
incomplete understanding of variant pathogenicity, missing genotype-phenotype
associations, complex gene-environment interactions, and reporting differences
between clinical labs. Maintaining a clear understanding of the rapidly
evolving landscape of diagnostic tests beyond ES, and their limitations,
presents a challenge for non-genetics professionals. Newer tests, such as
short-read genome or RNA sequencing, can be challenging to order and emerging
technologies, such as optical genome mapping and long-read DNA or RNA
sequencing, are not available clinically. Furthermore, there is no clear
guidance on the next best steps after inconclusive evaluation. Here, we review
why a clinical genetic evaluation may be negative, discuss questions to be
asked in this setting, and provide a framework for further investigation,
including the advantages and disadvantages of new approaches that are nascent
in the clinical sphere. We present a guide for the next best steps after
inconclusive molecular testing based upon phenotype and prior evaluation,
including when to consider referral to a consortium such as GREGoR, which is
focused on elucidating the underlying cause of rare unsolved genetic disorders.


  Understanding the interplay between recombination and resampling is a
significant challenge in mathematical population genetics and of great
practical relevance. Asymptotic results about the distribution of samples when
recombination is strong compared to resampling are often based on the
approximate solution of certain recursions, which is technically hard and
offers little conceptual insight. This work generalises an elegant
probabilistic argument, based on the coupling of ancestral processes but so far
only available in the case of two sites, to the multilocus setting. This offers
an alternative route to, and slightly generalises, a classical result of
Bhaskar and Song.


  Biological and social scientists have long been interested in understanding
how to reconcile individual and collective interests in iterated Prisoner's
Dilemma. Many effective strategies have been proposed, and they are often
categorized into one of two classes, `partners' and `rivals.' More recently,
another class, `friendly rivals,' has been identified in longer-memory strategy
spaces. Friendly rivals qualify as both partners and rivals: They fully
cooperate with themselves, like partners, but never allow their co-players to
earn higher payoffs, like rivals. Although they have appealing theoretical
properties, it is unclear whether they would emerge in evolving population
because most previous works focus on memory-one strategy space, where no
friendly rival strategy exists. To investigate this issue, we have conducted
large-scale evolutionary simulations in well-mixed and group-structured
populations and compared the evolutionary dynamics between memory-one and
memory-three strategy spaces. In a well-mixed population, the memory length
does not make a major difference, and the key factors are the population size
and the benefit of cooperation. Friendly rivals play a minor role because being
a partner or a rival is often good enough in a given environment. It is in a
group-structured population that memory length makes a stark difference: When
memory-three strategies are available, friendly rivals become dominant, and the
cooperation level nearly reaches a maximum, even when the benefit of
cooperation is so low that cooperation would not be achieved in a well-mixed
population. This result highlights the important interaction between group
structure and memory lengths that drive the evolution of cooperation.


  The goal of Protein Structure Prediction (PSP) problem is to predict a
protein's 3D structure (confirmation) from its amino acid sequence. The problem
has been a 'holy grail' of science since the Noble prize-winning work of
Anfinsen demonstrated that protein conformation was determined by sequence. A
recent and important step towards this goal was the development of AlphaFold2,
currently the best PSP method. AlphaFold2 is probably the highest profile
application of AI to science. Both AlphaFold2 and RoseTTAFold (another
impressive PSP method) have been published and placed in the public domain
(code & models). Stacking is a form of ensemble machine learning ML in which
multiple baseline models are first learnt, then a meta-model is learnt using
the outputs of the baseline level model to form a model that outperforms the
base models. Stacking has been successful in many applications. We developed
the ARStack PSP method by stacking AlphaFold2 and RoseTTAFold. ARStack
significantly outperforms AlphaFold2. We rigorously demonstrate this using two
sets of non-homologous proteins, and a test set of protein structures published
after that of AlphaFold2 and RoseTTAFold. As more high quality prediction
methods are published it is likely that ensemble methods will increasingly
outperform any single method.


  The adaptation rate in theoretical models of biological evolution increases
with the mutation rate but only to a point when mutations into lethal states
cause extinction. One would expect that removing such states should be
beneficial for evolution. We show here that, counter-intuitively, lethal
mutations speed up adaptation on rugged fitness landscapes with many fitness
maxima and minima, if strong competition for resources exist. We consider a
modified stochastic version of the quasispecies model with two types of
genotypes, viable and lethal, and show that increasing the rate of lethal
mutations decreases the time to evolve the best-fit genotype. This can be
explained by an increased rate of crossing fitness valleys, facilitated by
reduced selection against less-fit variants.


  Large-scale events like the UEFA Euro~2020 football (soccer) championship
offer a unique opportunity to quantify the impact of gatherings on the spread
of COVID-19, as the number and dates of matches played by participating
countries resembles a randomized study. Using Bayesian modeling and the gender
imbalance in COVID-19 data, we attribute 840,000 (95% CI: [0.39M, 1.26M])
COVID-19 cases across 12 countries to the championship. The impact depends
non-linearly on the initial incidence, the reproduction number $R$, and the
number of matches played. The strongest effects are seen in Scotland and
England, where as much as 10,000 primary cases per million inhabitants occur
from championship-related gatherings. The average match-induced increase in $R$
was 0.46 [0.18, 0.75] on match days, but important matches caused an increase
as large as +3. Altogether, our results provide quantitative insights that help
judge and mitigate the impact of large-scale events on pandemic spread.


  Since the foundation of population genetics, it is believed that directional
selection should reduce genetic diversity. We present an exactly solvable
population model which contradicts this intuition. The population is modelled
as a cloud of particles evolving in a 1-dimensional fitness space (fitness
wave). We show the existence of a phase transition which separates the
parameter space into a weak and a strong selection regimes. We find that
genetic diversity is highly non-monotone in the selection strength and, in
contrast with the common intuition, our model predicts that genetic diversity
is typically higher in the strong selection regime. This apparent paradox is
resolved by observing that a higher selection strength increases the absolute
fitness of the wave, but typically generate lower relative fitness between
individuals within the wave. These findings entail that inferring the magnitude
of natural selection from genetic data may raise some serious conceptual
issues.
  Along the way, we uncover a new phase transition in front propagation.
Namely, we show that the transition from weak to strong selection can be
reformulated in terms of a transition from fully-pulled to semi-pulled waves.
This transition is the pulled analog to the semi-pushed to fully-pushed regimes
observed in noisy-FKPP travelling waves in the presence of Allee effect.


  Cell membranes phase separate into ordered ${\rm L_o}$ and disordered ${\rm
L_d}$ domains depending on their compositions. This membrane
compartmentalization is heterogeneous and regulates the localization of
specific proteins related to cell signaling and trafficking. However, it is
unclear how the heterogeneity of the membranes affects the diffusion and
localization of proteins in ${\rm L_o}$ and ${\rm L_d}$ domains. Here, using
Langevin dynamics simulations coupled with the phase-field (LDPF) method, we
investigate several tens of milliseconds-scale diffusion and localization of
proteins in heterogeneous biological membrane models showing phase separation
into ${\rm L_o}$ and ${\rm L_d}$ domains. The diffusivity of proteins exhibits
temporal fluctuations depending on the field composition. Increases in
molecular concentrations and domain preference of the molecule induce
subdiffusive behavior due to molecular collisions by crowding and confinement
effects, respectively. Moreover, we quantitatively demonstrate that the protein
partitioning into the ${\rm L_o}$ domain is determined by the difference in
molecular diffusivity between domains, molecular preference of domain, and
molecular concentration. These results pave the way for understanding how
biological reactions caused by molecular partitioning may be controlled in
heterogeneous media. Moreover, the methodology proposed here is applicable not
only to biological membrane systems but also to the study of diffusion and
localization phenomena of molecules in various heterogeneous systems.


  Genomic surveillance of infectious diseases allows monitoring circulating and
emerging variants and quantifying their epidemic potential. However, due to the
high costs associated with genomic sequencing, only a limited number of samples
can be analysed. Thus, it is critical to understand how sampling impacts the
information generated. Here, we combine a compartmental model for the spread of
COVID-19 (distinguishing several SARS-CoV-2 variants) with different sampling
strategies to assess their impact on genomic surveillance. In particular, we
compare adaptive sampling, i.e., dynamically reallocating resources between
screening at points of entry and inside communities, and constant sampling,
i.e., assigning fixed resources to the two locations. We show that adaptive
sampling uncovers new variants up to five weeks earlier than constant sampling,
significantly reducing detection delays and estimation errors. This advantage
is most prominent at low sequencing rates. Although increasing the sequencing
rate has a similar effect, the marginal benefits of doing so may not always
justify the associated costs. Consequently, it is convenient for countries with
comparatively few resources to operate at lower sequencing rates, thereby
profiting the most from adaptive sampling. Finally, our methodology can be
readily adapted to study undersampling in other dynamical systems.


  During cellular translation, incorporation errors occur. It is the addition
of amino acid residues not corresponding to the mRNA code. With an increase in
the number of residues in the synthesized molecule, the probability of failure
in at least one link increases rapidly, which leads to improper folding and
loss of functionality of the entire molecule. A simple mathematical model is
presented, which shows that the amplification factor equals approximately the
number of links in the synthesised sequence. We assume that the enzymatic
processes of recognition of amino acids and their addition to the synthesized
molecule include the formation of intermediate pairs of radicals with
spin-correlated electrons. If a weak external magnetic field slightly changes
the rate of quantum singlet-triplet conversion, then a significant change in
the probability of occurrence of incorrect sequences of ribosomal translation
occurs.


  Discrete time, spatially extended models play an important role in ecology,
modelling population dynamics of species ranging from micro-organisms to birds.
An important question is how 'bottom up', individual-based models can be
approximated by 'top down' models of dynamics. Here, we study a class of
spatially explicit individual-based models with contest competition: where
species compete for space in local cells and then disperse to nearby cells. We
start by describing simulations of the model, which exhibit large-scale
discrete oscillations and characterise these oscillations by measuring spatial
correlations. We then develop two new approximate descriptions of the resulting
spatial population dynamics. The first is based on local interactions of the
individuals and allows us to give a difference equation approximation of the
system over small dispersal distances. The second approximates the long-range
interactions of the individual-based model. These approximations capture
demographic stochasticity from the individual-based model and show that
dispersal stabilizes population dynamics. We calculate extinction probability
for the individual-based model and show convergence between the local
approximation and the non-spatial global approximation of the individual-based
model as dispersal distance and population size simultaneously tend to
infinity. Our results provide new approximate analytical descriptions of a
complex bottom-up model and deepen understanding of spatial population
dynamics.


  Wolbachia infection in Anopheles albimanus mosquitoes can render mosquitoes
less capable of spreading malaria. We develop and analyze an ordinary
differential equation model to evaluate the effectiveness of Wolbachia-based
vector control strategies among wild Anopheles mosquitoes in Haiti. The model
tracks the mosquito life stages, including egg, larva, and adult (male and
female). It also accounts for critical biological effects, such as the maternal
transmission of Wolbachia through infected females and cytoplasmic
incompatibility, which effectively sterilizes uninfected females when they mate
with infected males. We derived and interpreted dimensionless numbers,
including the basic reproductive number and next-generation numbers. The
proposed system presents backward bifurcation, which indicates a threshold
infection that needs to be exceeded to establish a stable Wolbachia infection.
The sensitivity analysis ranks the relative importance of the epidemiological
parameters at the baseline. We simulate different intervention scenarios,
including pre-release mitigation using larviciding and thermal fogging before
the release, multiple releases of infected populations, and different release
timing. Our simulations show that the most efficient approach to establishing
Wolbachia is to release all the infected mosquitoes immediately after the
pre-release mitigation process. Also, the model predicts that it is more
efficient to release during the dry season than the wet season.


  Mapping spatial distributions of transcriptomic cell types is essential to
understanding the brain, with its exceptional cellular heterogeneity and the
functional significance of its spatial organization. Spatial transcriptomics
techniques are hoped to accomplish these measurements, but each method uses
different experimental and computational protocols, with different trade-offs
and optimizations. In 2017, the SpaceTx Consortium was formed to compare these
methods and determine their suitability for large-scale spatial transcriptomic
atlases. SpaceTx work included progress in tissue processing, taxonomy
development, gene selection, image processing and data standardization, cell
segmentation, cell type assignments, and visualization. Although the landscape
of experimental methods has changed dramatically since the beginning of
SpaceTx, the need for quantitative and detailed benchmarking of spatial
transcriptomics methods in the brain is still unmet. Here, we summarize the
work of SpaceTx and highlight outstanding challenges as spatial transcriptomics
grows into a mature field. We also discuss how our progress provides a roadmap
for benchmarking spatial transcriptomics methods in the future. Data and
analyses from this consortium, along with code and methods are publicly
available at https://spacetx.github.io/.


  One possible explanation for magnetosensing in biology, such as avian
magnetoreception, is based on the spin dynamics of certain chemical reactions
that involve radical pairs. Radical pairs have been suggested to also play a
role in anesthesia, hyperactivity, neurogenesis, circadian clock rhythm,
microtubule assembly, etc. It thus seems critical to probe the credibility of
such models. One way to do so is through isotope effects with different nuclear
spins. Here we briefly review the papers involving spin-related isotope effects
in biology. We suggest studying isotope effects can be an interesting avenue
for quantum biology.


  Various non pharmaceutical interventions have been settled to minimise the
burden of the COVID-19 outbreak. We build a framework to analyse the dynamics
of non pharmaceutical interventions, to distinguish between mitigations
measures leading to objective scientific improvements and mitigations based on
both political and scientific considerations. We analyse two possible
strategies within this framework. Namely, we consider mitigations driven by the
limited resources of the health system and mitigations where a constant set of
measures is applied at different moments. We describe the optimal interventions
for these scenarios. Our approach is mathematical and involves sir differential
systems, it is qualitative and geometrical rather than computational. Along
with the analysis of these scenarios, we collect several results that may be
useful on their own, in particular on the ground when the variables are not
known in real time.


  The use of biodegradation as a method for cleaning up soil that has been
contaminated by spilt petroleum can be an effective strategy. So, this study
investigated the existence of the wild microorganism in soil contaminated with
oil and study their ability to degrade petroleum in vitro. Nineteen samples
were collected from various locations near Taq Taq (TTOPCO) natural seeps in
the Kurdistan Region of Iraq. Morphological, cultural, biochemical tests and
molecular identification were used to identify the microbial communities, in
addition, spore texture and the colour of the fungal isolates were investigated
on the fungal isolates. Out of the19 samples, 17 indigenous bacterial strains
and 5 fungal strains were successfully isolated. From the absorption
spectrophotometry, Bacillus anthracis, Bacillus cereus, Achromobacter sp. and
Pseudomonas aeruginosa for the bacterial isolates grew well on a minimal salt
medium supplemented with 1% crude oil. Results showed that these isolates
mentioned above had a strong ability to degrade crude oil by reducing the
colour of 2,6-dichlorophenol indophenol (DCPIP) from deep blue to colourless.
However, for the fractions of hydrocarbon, the bacterial isolates failed and
did not affect the colour of any of the fractions. The results for fungi showed
that Aspergillus lentulus and Rhizopus arrhizus had a strong ability to degrade
both crude oil and fraction F1 by reducing the colour of DCPIP. Each fungal
isolates also had a great tolerance to different concentrations of crude oil
when grown on solid MSM. This study showed these microorganisms have a strong
ability to degrade crude oil and can be used to clean up soil and the
environment.


  With the increasing affordability and accessibility of genome sequencing
data, de novo genome assembly is an important first step to a wide variety of
downstream studies and analyses. Therefore, bioinformatics tools that enable
the generation of high-quality genome assemblies in a computationally efficient
manner are essential. Recent developments in long-read sequencing technologies
have greatly benefited genome assembly work, including scaffolding, by
providing long-range evidence that can aid in resolving the challenging
repetitive regions of complex genomes. ntLink is a flexible and
resource-efficient genome scaffolding tool that utilizes long-read sequencing
data to improve upon draft genome assemblies built from any sequencing
technologies, including the same long reads. Instead of using read alignments
to identify candidate joins, ntLink utilizes minimizer-based mappings to infer
how input sequences should be ordered and oriented into scaffolds. Recent
improvements to ntLink have added important features such as overlap detection,
gap-filling and in-code scaffolding iterations. Here, we present three basic
protocols demonstrating how to use each of these new features to yield highly
contiguous genome assemblies, while still maintaining ntLink's proven
computational efficiency. Further, as we illustrate in the alternate protocols,
the lightweight minimizer-based mappings that enable ntLink scaffolding can
also be utilized for other downstream applications, such as misassembly
detection. With its modularity and multiple modes of execution, ntLink has
broad benefit to the genomics community, from genome scaffolding and beyond.
ntLink is an open-source project and is freely available from
https://github.com/bcgsc/ntLink.


  Recent outbreaks of monkeypox and Ebola, and worrying waves of COVID-19,
influenza and respiratory syncytial virus, have all led to a sharp increase in
the use of epidemiological models to estimate key epidemiological parameters.
The feasibility of this estimation task is known as the practical
identifiability (PI) problem. Here, we investigate the PI of eight commonly
reported statistics of the classic Susceptible-Infectious-Recovered model using
a new measure that shows how much a researcher can expect to learn in a
model-based Bayesian analysis of prevalence data. Our findings show that the
basic reproductive number and final outbreak size are often poorly identified,
with learning exceeding that of individual model parameters only in the early
stages of an outbreak. The peak intensity, peak timing, and initial growth rate
are better identified, being in expectation over 20 times more probable having
seen the data by the time the underlying outbreak peaks. We then test PI for a
variety of true parameter combinations, and find that PI is especially
problematic in slow-growing or less-severe outbreaks. These results add to the
growing body of literature questioning the reliability of inferences from
epidemiological models when limited data are available.


  The main objective of this work is to test whether some stochastic models
typically used in financial markets could be applied to the COVID-19 pandemic.
To this end we have implemented the ARIMAX and Cox-Ingersoll-Ross (CIR) models
originally designed for interest rate pricing but transformed by us into a
forecasting tool. For the latter, which we denoted CIR*, both the
Euler-Maruyama method and the Milstein method were used. Forecasts obtained
with the maximum likelihood method have been validated with 95\% confidence
intervals and with statistical measures of goodness of fit, such as the root
mean square error (RMSE). We demonstrate that the accuracy of the obtained
results is consistent with the observations and sufficiently accurate to the
point that the proposed CIR* framework could be considered a valid alternative
to the classical ARIMAX for modelling pandemics.


  Despite widespread and striking examples of physiological oscillations, their
functional role is often unclear. Even glycolysis, the paradigm example of
oscillatory biochemistry, has seen questions about its oscillatory function.
Here, we take a systems approach to summarize evidence that oscillations play
critical physiological roles. Oscillatory behavior enables systems to avoid
desensitization, to avoid chronically high and therefore toxic levels of
chemicals, and to become more resistant to noise. Oscillation also enables
complex physiological systems to reconcile incompatible conditions such as
oxidation and reduction, by cycling between them, and to synchronize the
oscillations of many small units into one large effect. In pancreatic beta
cells, glycolytic oscillations are in synchrony with calcium and mitochondrial
oscillations to drive pulsatile insulin release, which is pivotal for the liver
to regulate blood glucose dynamics. In addition, oscillation can keep
biological time, essential for embryonic development in promoting cell
diversity and pattern formation. The functional importance of oscillatory
processes requires a rethinking of the traditional doctrine of homeostasis,
holding that physiological quantities are maintained at constant equilibrium
values, a view that has largely failed us in the clinic. A more dynamic
approach will enable us to view health and disease through a new light and
initiate a paradigm shift in treating diseases, including depression and
cancer. This modern synthesis also takes a deeper look into the mechanisms that
create, sustain and abolish oscillatory processes, which requires the language
of nonlinear dynamics, well beyond the linearization techniques of equilibrium
control theory.


  Nanopore sequencers generate electrical raw signals in real-time while
sequencing long genomic strands. These raw signals can be analyzed as they are
generated, providing an opportunity for real-time genome analysis. An important
feature of nanopore sequencing, Read Until, can eject strands from sequencers
without fully sequencing them, which provides opportunities to computationally
reduce the sequencing time and cost. However, existing works utilizing Read
Until either 1) require powerful computational resources that may not be
available for portable sequencers or 2) lack scalability for large genomes,
rendering them inaccurate or ineffective.
  We propose RawHash, the first mechanism that can accurately and efficiently
perform real-time analysis of nanopore raw signals for large genomes using a
hash-based similarity search. To enable this, RawHash ensures the signals
corresponding to the same DNA content lead to the same hash value, regardless
of the slight variations in these signals. RawHash achieves an accurate
hash-based similarity search via an effective quantization of the raw signals
such that signals corresponding to the same DNA content have the same quantized
value and, subsequently, the same hash value.
  We evaluate RawHash on three applications: 1) read mapping, 2) relative
abundance estimation, and 3) contamination analysis. Our evaluations show that
RawHash is the only tool that can provide high accuracy and high throughput for
analyzing large genomes in real-time. When compared to the state-of-the-art
techniques, UNCALLED and Sigmap, RawHash provides 1) 25.8x and 3.4x better
average throughput and 2) significantly better accuracy for large genomes,
respectively. Source code is available at
https://github.com/CMU-SAFARI/RawHash.


  Non-standard Lagrangians do not display any discernible energy-like terms,
yet they give the same equations of motion as standard Lagrangians, which have
easily identifiable energy-like terms. A new method to derive non-standard
Lagrangians for second-order nonlinear differential equations with damping is
developed and the limitations of this method are explored. It is shown that the
limitations do not exist only for those nonlinear dynamical systems that can be
converted into linear ones. The obtained results are applied to selected
population dynamics models for which non-standard Lagrangians and their
corresponding null Lagrangians and gauge functions are derived, and their roles
in the population dynamics are discussed.


  Facing an alarming continuing decline of wild sea cucumber resources,
management strategies were developed over the past three decades to sustainably
promote development, maintenance, or regeneration of wild sea cucumber
fisheries. In New Caledonia (South Pacific), dedicated management efforts via
restocking and sea ranching programs were implemented to cope with the
overharvesting of the sandfish Holothuria scabra and the recent loss of known
populations. In order to investigate genetic implications of a major H. scabra
restocking program, we assessed the genetic diversity and structure of wild
stocks and hatchery-produced sandfish and compared the genetic outcomes of
consecutive spawning and juvenile production events. For this, 1358 sandfish
collected at four sites along the northwestern coasts of New Caledonia, as well
as during five different restocking events in the Tiabet Bay, were genotyped
using nine polymorphic microsatellite markers. We found that wild H. scabra
populations from the northwestern coast of New Caledonia likely belonged to one
panmictic population with high level of gene flow observed along the study
scale. Further, this panmictic population displayed an effective size of
breeders large enough to ensure the feasibility of appropriate breeding
programs for restocking. In contrast, hatchery-produced samples did suffer from
an important reduction in the effective population size: the effective
population size were so small that genetic drift was detectable over one
generation, with the presence of inbred individuals, as well as more related
dyads than in wild populations. All these results suggest that dedicated
efforts in hatcheries are further needed to maintain genetic diversity of
hatchery-produced individuals in order to unbalance any negative impact during
this artificial selection.


  The processes of ecological interactions, dispersal and mutations shape the
dynamics of biological communities, and analogous eco-evolutionary processes
acting upon economic entities have been proposed to explain economic change.
This hypothesis is compelling because it explains economic change through
endogenous mechanisms, but it has not been quantitatively tested at the global
economy level. Here, we use an inverse modelling technique and 59 years of
economic data covering 77 countries to test whether the collective dynamics of
national economic activities can be characterised by eco-evolutionary
processes. We estimate the statistical support of dynamic community models in
which the dynamics of economic activities are coupled with positive and
negative interactions between the activities, the spatial dispersal of the
activities, and their transformations into other economic activities. We find
strong support for the models capturing positive interactions between economic
activities and spatial dispersal of the activities across countries. These
results suggest that processes akin to those occurring in ecosystems play a
significant role in the dynamics of economic systems. The strength-of-evidence
obtained for each model varies across countries and may be caused by
differences in the distance between countries, specific institutional contexts,
and historical contingencies. Overall, our study provides a new quantitative,
biologically inspired framework to study the forces shaping economic change.


  In this work, we propose a general inversion framework to non-uniquely invert
a very large class of ordinary differential equations (ODEs) into chemical
reaction networks. A thorough treatment of the relevant chemical reaction
network theory from the literature is given. Various simulation results are
provided to augment the selection procedure for the inverse framework, where a
previously known kineticization strategy is shown to be deterministically
excellent but undesirable in chemical simulations. The utility of the framework
is verified by simulating reaction network forms of meaningful ODE systems, and
their time series are analyzed. In particular, we provide simulations of
deterministic chaotic attractors whose newly discovered reaction networks are
non-equivalent with any existing chemical interpretations within the
literature, as well as presenting exemplary figures which may form a roadmap to
the successful biochemical implementation of the integration of ODE systems.


  According to the evolutionary death-birth protocol, a player is chosen
randomly to die and neighbors compete for the available position proportional
to their fitness. Hence, the status of the focal player is completely ignored
and has no impact on the strategy update. In this work, we revisit and
generalize this rule by introducing a weight factor to compare the payoff
values of the focal and invading neighbors. By means of evolutionary graph
theory, we analyze the model on joint transitive graphs to explore the possible
consequences of the presence of a weight factor. We find that focal weight
always hinders cooperation under weak selection strength. Surprisingly, the
results show a non-trivial tipping point of the weight factor where the
threshold of cooperation success shifts from positive to negative infinity.
Once focal weight exceeds this tipping point, cooperation becomes unreachable.
Our theoretical predictions are confirmed by Monte Carlo simulations on a
square lattice of different sizes. We also verify the robustness of the
conclusions to arbitrary two-player prisoner's dilemmas, to dispersal graphs
with arbitrary edge weights, and to interaction and dispersal graphs
overlapping arbitrarily.


  For reaction networks arising in systems biology, the capacity for two or
more steady states, that is, multistationarity, is an important property that
underlies biochemical switches. Another property receiving much attention
recently is absolute concentration robustness (ACR), which means that some
species concentration is the same at all positive steady states. In this work,
we investigate the prevalence of each property while paying close attention to
when the properties occur together. Specifically, we consider a stochastic
block framework for generating random networks, and prove edge-probability
thresholds at which - with high probability - multistationarity appears and ACR
becomes rare. We also show that the small window in which both properties occur
only appears in networks with many species. Taken together, our results confirm
that, in random reversible networks, ACR and multistationarity together, or
even ACR on its own, is highly atypical. Our proofs rely on two prior results,
one pertaining to the prevalence of networks with deficiency zero, and the
other "lifting" multistationarity from small networks to larger ones.


  How the architecture of gene regulatory networks ultimately shapes gene
expression patterns is an open question, which has been approached from a
multitude of angles. The dominant strategy has been to identify non-random
features in these networks and then argue for the function of these features
using mechanistic modelling. Here we establish the foundation of an alternative
approach by studying the correlation of eigenvectors with synthetic gene
expression data simulated with a basic and popular model of gene expression
dynamics -- attractors of Boolean threshold dynamics in signed directed graphs.
  Eigenvectors of the graph Laplacian are known to explain collective dynamical
states (stationary patterns) in Turing dynamics on graphs. In this study, we
show that eigenvectors can also predict collective states (attractors) for a
markedly different type of dynamics, Boolean threshold dynamics, and category
of graphs, signed directed graphs. However, the overall predictive power
depends on details of the network architecture, in a predictable fashion.
  Our results are a set of statistical observations, providing the first
systematic step towards a further theoretical understanding of the role of
eigenvectors in dynamics on graphs.


  Lycopene, rich in red, yellow, or orange-colored fruits and vegetables, is
the most potent antioxidant among the other carotenoids available in human
blood plasma. It is evident that regular lycopene intake can prevent chronic
diseases like cardiovascular diseases, type-2 diabetes, hypertension, kidney
diseases and cancer. However, thermal processing, light, oxygen, and enzymes in
gastrointestinal tract (GIT) compromise the bioaccessibility and
bioavailability of lycopene ingested through diet. Nanoencapsulation provides a
potential platform to prevent lycopene from light, air oxygen, thermal
processing and enzymatic activity of the human digestive system.
Physicochemical properties evidenced to be the potential indicator for
determining the bioaccessibility of encapsulated bioactive compounds like
lycopene. By manipulating the size or hydrodynamic diameter, zeta potential
value or stability, polydispersity index or homogeneity and functional activity
or retention of antioxidant properties observed to be the most prominent
physicochemical properties to evaluate beneficial effect of implementation of
nanotechnology on bioaccessibility study. Moreover, the molecular mechanism of
the bioavailability of nanoparticles is not yet to be understood due to lack of
comprehensive design to identify nanoparticles' behaviors if ingested through
oral route as functional food ingredients. This review paper aims to study and
leverage existing techniques about how nanotechnology can be used and verified
to identify the bioaccessibility of lycopene before using it as a functional
food ingredient for therapeutic treatments.


  The major histocompatibility complex (MHC) molecules, which bind peptides for
presentation on the cell surface, play an important role in cell-mediated
immunity. In light of developing databases and technologies over the years,
significant progress has been made in research on peptide binding affinity
calculation. Several in techniques have been developed to predict peptide
binding to MHC class I. Most of the research on MHC Class I due to its nature
brings better performance and more. Considering the use of different methods
and different technologies, and the approach of similar methods on different
proteins, a classification was created according to the binding affinity of
protein peptides.


  We present a proof based on SEIR-models that shows it is impossible to
identify the level of under reporting based on traditional observables of the
disease dynamics alone. This means that the true attack rate must be determined
through other means.


  While artificial intelligence has made remarkable strides in revealing the
relationship between biological macromolecules' primary sequence and tertiary
structure, designing RNA sequences based on specified tertiary structures
remains challenging. Though existing approaches in protein design have
thoroughly explored structure-to-sequence dependencies in proteins, RNA design
still confronts difficulties due to structural complexity and data scarcity.
Adding to the problem, direct transplantation of protein design methodologies
into RNA design fails to achieve satisfactory outcomes although sharing similar
structural components. In this study, we aim to systematically construct a
data-driven RNA design pipeline. We crafted a large, well-curated benchmark
dataset and designed a comprehensive structural modeling approach to represent
the complex RNA tertiary structure. More importantly, we proposed a
hierarchical data-efficient representation learning framework that learns
structural representations through contrastive learning at both cluster-level
and sample-level to fully leverage the limited data. By constraining data
representations within a limited hyperspherical space, the intrinsic
relationships between data points could be explicitly imposed. Moreover, we
incorporated extracted secondary structures with base pairs as prior knowledge
to facilitate the RNA design process. Extensive experiments demonstrate the
effectiveness of our proposed method, providing a reliable baseline for future
RNA design tasks. The source code and benchmark dataset will be released
publicly.


  Protein-ligand binding prediction is a fundamental problem in AI-driven drug
discovery. Prior work focused on supervised learning methods using a large set
of binding affinity data for small molecules, but it is hard to apply the same
strategy to other drug classes like antibodies as labelled data is limited. In
this paper, we explore unsupervised approaches and reformulate binding energy
prediction as a generative modeling task. Specifically, we train an
energy-based model on a set of unlabelled protein-ligand complexes using SE(3)
denoising score matching and interpret its log-likelihood as binding affinity.
Our key contribution is a new equivariant rotation prediction network called
Neural Euler's Rotation Equations (NERE) for SE(3) score matching. It predicts
a rotation by modeling the force and torque between protein and ligand atoms,
where the force is defined as the gradient of an energy function with respect
to atom coordinates. We evaluate NERE on protein-ligand and antibody-antigen
binding affinity prediction benchmarks. Our model outperforms all unsupervised
baselines (physics-based and statistical potentials) and matches supervised
learning methods in the antibody case.


  Semi-directed phylogenetic networks have recently emerged as a class of
phylogenetic networks sitting between rooted (directed) and unrooted
(undirected) phylogenetic networks as they contain both directed as well as
undirected edges. While the spaces of rooted phylogenetic networks and unrooted
phylogenetic networks have been analyzed in recent years and various
rearrangement moves to traverse these spaces have been introduced, the results
do not immediately carry over to semi-directed phylogenetic networks. Here, we
propose a simple rearrangement move for semi-directed phylogenetic networks,
called cut edge transfer (CET), and show that the space of semi-directed
level-$1$ networks with precisely $k$ reticulations is connected under CET.
This level-$1$ space is currently the predominantly used search space for most
algorithms that reconstruct semi-directed phylogenetic networks. Hence, every
semi-directed level-$1$ network with a fixed number of reticulations and leaf
set can be reached from any other such network by a sequence of CETs. By
introducing two additional moves, CET$^+$ and CET$^-$, that allow for the
addition or deletion of reticulations, we then establish connectedness for the
space of all semi-directed phylogenetic networks on a fixed leaf set. As a
byproduct of our results for semi-directed phylogenetic networks, we also show
that the space of rooted level-$1$ networks with a fixed number of
reticulations and leaf set is connected under CET, when translated into the
rooted setting.


  The use of biodegradable polymers simplifies the development of therapeutic
devices with regards to transient implants and three-dimensional platform
suitable for tissue engineering. Further advances have also occurred in the
controlled released mechanism of bioactive compounds encapsulated in
biodegradable polymers. This application requires the understanding of the
physicochemical properties of the polymeric materials and their inherent impact
on the delivery of encapsulated bioactive. Hence, the objective of this study
was to evaluate the effect of surfactant and sonication time on the
bio-accessibility of lycopene encapsulated polymeric nanoparticles. The
emulsion evaporation method was used to encapsulate lycopene in poly-lactic
co-glycolic acid (PLGA) with surfactant concentration, sonication time and
polymer concentration as independent variables. Physicochemical and
morphological characteristics were measured with a zetasizer and SEM, while the
encapsulation efficiency and controlled release kinetics with
spectrophotometric, and the dialysis method was used to estimate
bioaccessibility. The results have shown sonication time to have significantly
(p < 0.05) influenced the encapsulation efficiency. Hence, the sonication time
of 4 min yield an encapsulation efficiency of 78% and increased to 97% with
increase sonication time (6 min). Increased sonication time had a decreasing
effect on the hydrodynamic diameter and stability of the encapsulated
nanoparticles. The slow release of lycopene was observed during the first 12
days, followed by a burst release of about 44% on the 13th day in-vitro. The
study will have significant impact on the manufacturing of functional food with
encapsulated ingredients and provide an understanding of their inherent control
release mechanism in the GIT.


  There is significant interest in using existing repositories of biological
entities, relationships, and models to automate biological model assembly and
extension. Current methods aggregate human-curated biological information into
executable, simulatable models, but these models do not resemble human curated
models and do not recapitulate experimental results. Here, we outline the
process of automated model assembly and extension, while demonstrating it on
both synthetic models and human-curated models of biological signaling
networks. We begin with an iterative, greedy, and combinatoric approach to
automated assembly and demonstrate the key difficulties inherent to contextless
assembly. We publicly release the software used in this paper to enable further
exploration of this problem.


  Nanoencapsulation has become a widespread technique to improve the
bioavailability of bioactive compounds like lycopene. Consumption of lycopene
rich foods is effective in preventing cancer, diabetes, and cardiovascular
diseases due to its strong, oxygen-quenching ability. But the functional
activity of lycopene is compromised by light, oxygen, and heat. A biodegradable
polymer such as PLA or PLGA is the most effective carrier for encapsulating
lycopene due to its excellent biodegradability, biocompatibility, and
nontoxigenic effect on the human metabolic system. Hence, the primary objective
of this study was to evaluate the effect of pasteurization on bioaccessibility
and control release kinetics of encapsulated lycopene nanoparticles in invitro
human GIT. In the first objective, sonication time, surfactant, and polymer
concentration were considered as the three factors to synthesize polymeric
lycopene nanoparticles (LNP) whereas in the second objective, type of
pasteurization, encapsulation, and juice concentration were three factors. In
Objectives 4 and 5, the type of encapsulation and pasteurization, and digestion
time were evaluated as the three factors used to evaluate the in-vitro
bioaccessibility of encapsulated lycopene NP. The study evidenced that
encapsulation improved the lycopene's bioaccessibility by 70% and more than 60%
for conventional pasteurized (CP) and microwave pasteurized (MP) nanoemulsions,
respectively without compromising the physicochemical properties of both PLA
and PLGA-lycopene nanoparticles. The invitro bioaccessibility study also showed
that CP reduced the functional activity of PLALNP by 20% whereas MP had no
significant effect on the bioaccessibility of PLA LNP. It is unlikely that the
PLGA LNP was more sensitive against MP nanoemulsion than the CP. In conclusion,
it was shown that the PLA LNP treated with MP provided the highest
bioaccessibility.


  A two-dimensional homomorphic logistic map that preserves features of
Lotka-Volterra Equations was proposed. In order to examine the Lotka-Volterra
chaos, in addition to ordinary iteration plots of population, Lyapunov
exponents either calculated directly from eigenvalues of Jacobian of the $2$D
logistic mapping, or from time-series algorithms of both Rosenstein and Eckmann
et al. were calculated, among which discrepancies were compared. Bifurcation
diagrams may be divided into five categories depending on different topological
shapes, among which flip bifurcation and Neimark-Sacker bifurcation were
observed, the latter showing closed orbits around limit circles in the phase
portrait and phase space diagram. Our model restored the $1$D logistic map of
the prey at the absence of the predator, as well as the normal competing
behavior between two species when the initial population of the two is equal.
In spite of the possibility for two species going into chaos simultaneously, it
is also possible that with the same inter-species parameters as normal but with
predator population $10$ times more than that of the prey, under certain growth
rate the latter becomes chaotic, and former dramatically reduces to zero,
referring to total annihilation of the predator species. Interpreting humans as
the predator and natural resources as the prey in the ecological system, the
aforementioned conclusion may imply that not only excessive consumption of the
natural resources, but its chaotic state triggered by overpopulation of humans
may also backfire in a manner of total extinction on human species.
Fortunately, a little chance may exist for survival of human race, as isolated
fixed points in bifurcation diagram of the predator reveals.


  We study communities emerging from generalised random Lotka--Volterra
dynamics with a large number of species with interactions determined by the
degree of niche overlap. Each species is endowed with a number of traits, and
competition between pairs of species increases with their similarity in trait
space. This leads to a model with random Hopfield-like interactions. We use
tools from the theory of disordered systems, notably dynamic mean field theory,
to characterise the statistics of the resulting communities at stable fixed
points and determine analytically when stability breaks down. Two distinct
types of transition are identified in this way, both marked by diverging
abundances, but differing in the behaviour of the integrated response function.
At fixed points only a fraction of the initial pool of species survives. We
numerically study the eigenvalue spectra of the interaction matrix between
extant species. We find evidence that the two types of dynamical transition
are, respectively, associated with the bulk spectrum or an outlier eigenvalue
crossing into the right half of the complex plane.


  In medical genetics, each genetic variant is evaluated as an independent
entity regarding its clinical importance. However, in most complex diseases,
variant combinations in specific gene networks, rather than the presence of a
particular single variant, predominates. In the case of complex diseases,
disease status can be evaluated by considering the success level of a team of
specific variants. We propose a high dimensional modelling based method to
analyse all the variants in a gene network together. To evaluate our method, we
selected two gene networks, mTOR and TGF-Beta. For each pathway, we generated
400 control and 400 patient group samples. mTOR and TGF-? pathways contain 31
and 93 genes of varying sizes, respectively. We produced Chaos Game
Representation images for each gene sequence to obtain 2-D binary patterns.
These patterns were arranged in succession, and a 3-D tensor structure was
achieved for each gene network. Features for each data sample were acquired by
exploiting Enhanced Multivariance Products Representation to 3-D data. Features
were split as training and testing vectors. Training vectors were employed to
train a Support Vector Machines classification model. We achieved more than 96%
and 99% classification accuracies for mTOR and TGF-Beta networks, respectively,
using a limited amount of training samples.


  The 3-dimensional (3D) structure of the genome is of significant importance
for many cellular processes. In this paper, we study the problem of
reconstructing the 3D structure of chromosomes from Hi-C data of diploid
organisms, which poses additional challenges compared to the better-studied
haploid setting. With the help of techniques from algebraic geometry, we prove
that a small amount of phased data is sufficient to ensure finite
identifiability, both for noiseless and noisy data. In the light of these
results, we propose a new 3D reconstruction method based on semidefinite
programming, paired with numerical algebraic geometry and local optimization.
The performance of this method is tested on several simulated datasets under
different noise levels and with different amounts of phased data. We also apply
it to a real dataset from mouse X chromosomes, and we are then able to recover
previously known structural features.


  Livestock mobility, particularly that of small and large ruminants, is one of
the main pillars of production and trade in West Africa: livestock is moved
around in search of better grazing or sold in markets for domestic consumption
and for festival-related activities. These movements cover several thousand
kilometers and have the capability of connecting the whole West African region
thus facilitating the diffusion of many animal and zoonotic diseases. Several
factors shape mobility patterns even in normal years and surveillance systems
need to account for such changes. In this paper, we present a procedure based
on temporal network theory to identify possible sentinel locations using two
indicators: vulnerability (i.e. the probability of being reached by the
disease) and time of infection (i.e. the time of first arrival of the disease).
Using these indicators in our structural analysis of the changing network
enabled us to identify a set of nodes that could be used in an early warning
system. As a case study we simulated the introduction of F.A.S.T. (Foot and
Mouth Similar Transboundary) diseases in Senegal and used data taken from 2020
Sanitary certificates (LPS, laissez-passer sanitaire) issued by the Senegalese
Veterinary Services to reconstruct the national mobility network. Our analysis
showed that a static approach can significantly overestimate the speed and the
extent of disease propagation, whereas temporal analysis revealed that the
reachability and vulnerability of the different administrative departments
(used as nodes of the mobility network) change over the course of the year. For
this reason, several sets of sentinel nodes were identified in different
periods of the year, underlining the role of temporality in shaping patterns of
disease diffusion.


  Genealogical networks (i.e. family trees) are of growing interest, with the
largest known data sets now including well over one billion individuals.
Interest in family history also supports an 8.5 billion dollar industry whose
size is projected to double within 7 years (FutureWise report HC1137). Yet
little mathematical attention has been paid to the complex network properties
of genealogical networks, especially at large scales.
  The structure of genealogical networks is of particular interest due to the
practice of forming unions, e.g. marriages, that are typically well outside
one's immediate family. In most other networks, including other social
networks, no equivalent restriction exists on the distance at which
relationships form. To study the effect this has on genealogical networks we
use persistent homology to identify and compare the structure of 101
genealogical and 31 other social networks. Specifically, we introduce the
notion of a network's persistence curve, which encodes the network's set of
persistence intervals. We find that the persistence curves of genealogical
networks have a distinct structure when compared to other social networks. This
difference in structure also extends to subnetworks of genealogical and social
networks suggesting that, even with incomplete data, persistent homology can be
used to meaningfully analyze genealogical networks. Here we also describe how
concepts from genealogical networks, such as common ancestor cycles, are
represented using persistent homology. We expect that persistent homology tools
will become increasingly important in genealogical exploration as popular
interest in ancestry research continues to expand.


  Models of strategy evolution on static networks help us understand how
population structure can promote the spread of traits like cooperation. One key
mechanism is the formation of altruistic spatial clusters, where neighbors of a
cooperative individual are likely to reciprocate, which protects prosocial
traits from exploitation. But most real-world interactions are ephemeral and
subject to exogenous restructuring, so that social networks change over time.
Strategic behavior on dynamic networks is difficult to study, and much less is
known about the resulting evolutionary dynamics. Here, we provide an analytical
treatment of cooperation on dynamic networks, allowing for arbitrary spatial
and temporal heterogeneity. We show that transitions among a large class of
network structures can favor the spread of cooperation, even if each individual
social network would inhibit cooperation when static. Furthermore, we show that
spatial heterogeneity tends to inhibit cooperation, whereas temporal
heterogeneity tends to promote it. Dynamic networks can have profound effects
on the evolution of prosocial traits, even when individuals have no agency over
network structures.


  Current protein language models (PLMs) learn protein representations mainly
based on their sequences, thereby well capturing co-evolutionary information,
but they are unable to explicitly acquire protein functions, which is the end
goal of protein representation learning. Fortunately, for many proteins, their
textual property descriptions are available, where their various functions are
also described. Motivated by this fact, we first build the ProtDescribe dataset
to augment protein sequences with text descriptions of their functions and
other important properties. Based on this dataset, we propose the ProtST
framework to enhance Protein Sequence pre-training and understanding by
biomedical Texts. During pre-training, we design three types of tasks, i.e.,
unimodal mask prediction, multimodal representation alignment and multimodal
mask prediction, to enhance a PLM with protein property information with
different granularities and, at the same time, preserve the PLM's original
representation power. On downstream tasks, ProtST enables both supervised
learning and zero-shot prediction. We verify the superiority of ProtST-induced
PLMs over previous ones on diverse representation learning benchmarks. Under
the zero-shot setting, we show the effectiveness of ProtST on zero-shot protein
classification, and ProtST also enables functional protein retrieval from a
large-scale database without any function annotation.


  The reaction center consists of atoms in the product whose local properties
are not identical to the corresponding atoms in the reactants. Prior studies on
reaction center identification are mainly on semi-templated retrosynthesis
methods. Moreover, they are limited to single reaction center identification.
However, many reaction centers are comprised of multiple bonds or atoms in
reality. We refer to it as the multiple reaction center. This paper presents
RCsearcher, a unified framework for single and multiple reaction center
identification that combines the advantages of the graph neural network and
deep reinforcement learning. The critical insight in this framework is that the
single or multiple reaction center must be a node-induced subgraph of the
molecular product graph. At each step, it considers choosing one node in the
molecular product graph and adding it to the explored node-induced subgraph as
an action. Comprehensive experiments demonstrate that RCsearcher consistently
outperforms other baselines and can extrapolate the reaction center patterns
that have not appeared in the training set. Ablation experiments verify the
effectiveness of individual components, including the beam search and one-hop
constraint of action space.


  Antibodies are vital proteins offering robust protection for the human body
from pathogens. The development of general protein and antibody-specific
pre-trained language models both facilitate antibody prediction tasks. However,
there have been limited studies that comprehensively explore the representation
capability of distinct pre-trained language models on different antibody tasks.
To investigate the problem, we aim to answer several key questions in this
paper, such as how pre-trained language models perform in antibody tasks with
different specificity and how introducing specific biological mechanisms to the
pre-training process can benefit the model. Additionally, we evaluate if the
learned antibody pre-trained representations can be applied to real-world
antibody problems, like drug discovery and immune process understanding.
Previously, no benchmark available largely hindered the study to answer these
questions. To aid in our investigation, we provide an AnTibody Understanding
Evaluation (ATUE) benchmark. We comprehensively evaluate the performance of
protein pre-trained language models by empirical study along with conclusions
and new insights. Our ATUE and code are released at
https://github.com/dqwang122/EATLM.


  Motivation: As viruses that mainly infect bacteria, phages are key players
across a wide range of ecosystems. Analyzing phage proteins is indispensable
for understanding phages' functions and roles in microbiomes. High-throughput
sequencing enables us to obtain phages in different microbiomes with low cost.
However, compared to the fast accumulation of newly identified phages, phage
protein classification remains difficult. In particular, a fundamental need is
to annotate virion proteins, the structural proteins such as major tail,
baseplate etc. Although there are experimental methods for virion protein
identification, they are too expensive or time-consuming, leaving a large
number of proteins unclassified. Thus, there is a great demand to develop a
computational method for fast and accurate phage virion protein classification.
Results: In this work, we adapted the state-of-the-art image classification
model, Vision Transformer, to conduct virion protein classification. By
encoding protein sequences into unique images using chaos gaming
representation, we can leverage Vision Transformer to learn both local and
global features from sequence ``images''. Our method, PhaVIP, has two main
functions: classifying PVP and non-PVP sequences and annotating the types of
PVP, such as capsid and tail. We tested PhaVIP on several datasets with
increasing difficulty and benchmarked it against alternative tools. The
experimental results show that PhaVIP has superior performance. After
validating the performance of PhaVIP, we investigated two applications that can
use the output of PhaVIP: phage taxonomy classification and phage host
prediction. The results show the benefit of using classified proteins rather
than all proteins.


  Our study is based on an epidemiological compartmental model, the SIRS model.
In the SIRS model, each individual is in one of the states susceptible (S),
infected(I) or recovered (R), depending on its state of health. In compartment
R, an individual is assumed to stay immune within a finite time interval only
and then transfers back to the S compartment. We extend the model and allow for
a feedback control of the infection rate by mitigation measures which are
related to the number of infections. A finite response time of the feedback
mechanism is supposed that changes the low-dimensional SIRS model into an
infinite-dimensional set of integro-differential (delay-differential)
equations. It turns out that the retarded feedback renders the originally
stable endemic equilibrium of SIRS (stable focus) into an unstable focus if the
delay exceeds a certain critical value. Nonlinear solutions show persistent
regular oscillations of the number of infected and susceptible individuals. In
the last part we include noise effects from the environment and allow for a
fluctuating infection rate. This results in multiplicative noise terms and our
model turns into a set of stochastic nonlinear integro-differential equations.
Numerical solutions reveal an irregular behavior of repeated disease outbreaks
in the form of infection waves with a variety of frequencies and amplitudes.


  Proteins power a vast array of functional processes in living cells. The
capability to create new proteins with designed structures and functions would
thus enable the engineering of cellular behavior and development of
protein-based therapeutics and materials. Structure-based protein design aims
to find structures that are designable (can be realized by a protein sequence),
novel (have dissimilar geometry from natural proteins), and diverse (span a
wide range of geometries). While advances in protein structure prediction have
made it possible to predict structures of novel protein sequences, the
combinatorially large space of sequences and structures limits the practicality
of search-based methods. Generative models provide a compelling alternative, by
implicitly learning the low-dimensional structure of complex data
distributions. Here, we leverage recent advances in denoising diffusion
probabilistic models and equivariant neural networks to develop Genie, a
generative model of protein structures that performs discrete-time diffusion
using a cloud of oriented reference frames in 3D space. Through in silico
evaluations, we demonstrate that Genie generates protein backbones that are
more designable, novel, and diverse than existing models. This indicates that
Genie is capturing key aspects of the distribution of protein structure space
and facilitates protein design with high success rates. Code for generating new
proteins and training new versions of Genie is available at
https://github.com/aqlaboratory/genie.


  Bacteriophages are central to microbial ecosystems for balancing bacterial
populations and promoting evolution by applying strong selection pressure. Here
we review some of the known aspects that modulate phage-bacteria interaction in
a way that naturally promotes their coexistence. We focus on the modulations
that arise from structural, physical, or physiological constraints. We argue
they should play roles in many phage-bacteria systems providing sustainable
diversity.


  The soil seed bank is a major component of plant communities. However,
long-term analyses of the dynamics of the seed bank and the ensuing vegetation
are rare. Here, we studied the dynamics in plant communities with high
dominance of annuals in Mediterranean, semiarid, and arid ecosystems for nine
consecutive years. For annuals, we hypothesized that the density of the seed
bank would be more stable than the density of the standing herbaceous
vegetation. Moreover, we predicted that differences in temporal variability
between the seed bank and the vegetation would increase with aridity, where
year-to-year rainfall variability is higher. We found that the temporal
variability at the population level (assessed as the standard deviation of the
loge-transformed density) of the nine dominant annuals in each site did not
differ between the seed bank and the ensuing vegetation in any of the sites.
For the total density of annuals, patterns depended on aridity. In the
Mediterranean site, the temporal variability was similar in the seed bank and
the vegetation (0.40 vs. 0.40). Still, in the semiarid and arid sites,
variability in the seed bank was lower than in the vegetation (0.49 vs. 1.01
and 0.63 vs. 1.38, respectively). This difference between the population-level
patterns and the total density of annuals can be related to the lower
population synchrony in their seed bank. In contrast, for the herbaceous
perennials (all species combined), the seed bank variability was higher than in
the vegetation. Overall, our results highlight the role of the seed bank in
buffering the annual vegetation density with increasing climatic uncertainty
typical in aridity gradients. This role is crucial under the increasing
uncertainty imposed by climatic change in the region.


  The COVID-19 pandemic has given rise to numerous articles from different
scientific fields (epidemiology, virology, immunology, airflow physics...)
without any effort to link these different insights. In this review, we aim to
establish relationships between epidemiological data and the characteristics of
the virus strain responsible for the epidemic wave concerned. We have carried
out this study on the Wuhan, Alpha, Delta and Omicron strains allowing us to
illustrate the evolution of the relationships we have highlighted according to
these different viral strains. We addressed the following questions: 1) How can
the mean infectious dose (one quantum, by definition in epidemiology) be
measured and expressed as an amount of viral RNA molecules (in genome units,
GU) or as a number of replicative viral particles (in plaque-forming units,
PFU)? 2) How many infectious quanta are exhaled by an infected person per unit
of time? 3) How many infectious quanta are exhaled, on average, integrated over
the whole contagious period? 4) How do these quantities relate to the epidemic
reproduction rate R as measured in epidemiology, and to the viral load, as
measured by molecular biological methods? 5) How has the infectious dose
evolved with the different strains of SARS-CoV-2? We make use of
state-of-the-art modelling, reviewed and explained in the appendix of the
article (Supplemental Information, SI), to answer these questions using data
from the literature in both epidemiology and virology. We have considered the
modification of these relationships according to the vaccination status of the
population. We hope that this work will allow a better integration of data from
different fields (virology, epidemiology, and immunology) to anticipate the
evolution of the epidemic in the case of COVID-19, but also in respiratory
pathologies transmissible in an airborne manner.


  Cross-links are considered to be a key component of collagen fibrils as they
can change the fibrillar behavior in various ways. Advanced-Glycation
Endproducts (AGEs), one particular type of cross-links, have been shown to
accumulate and impair the mechanical properties of collageneous tissues,
whereas enzymatic cross-links (ECLs) are known for stabilizing the structure of
the fibril. However, the reasons for whether a given type of cross-link
improves or impairs the material properties remain unknown. Here, we use
coarse-grained steered molecular models to evaluate the effect of AGEs and ECLs
cross-links content on the deformation and failure properties of collagen
fibrils. Our simulations show that the collagen fibrils stiffen at high strain
levels when the AGEs content exceeds a critical value. In addition, the
strength of the fibril increases with AGEs accumulation. By analyzing the
forces within the different types of cross-links (AGEs and ECLs) as well as
their failure, we demonstrate that a change of deformation mechanism is at the
origin of these observations. A high AGEs content reinforces force transfer
through AGEs cross-links rather than through friction between sliding
tropocollagen molecules. We show that this failure mechanism, which is
associated with lower energy dissipation, results in more abrupt failure of the
collagen fibril. Our results provide a direct and causal link between increased
AGEs content, inhibited intra-fibrillar sliding, increased stiffness, and
abrupt fibril fracture. Therefore, they explain the mechanical origin of bone
brittleness as commonly observed in elderly and diabetic populations. Our
findings contribute to a better understanding of the mechanisms underlying
impaired tissue behaviour due to elevated AGEs content and could enable
targeted measures regarding the reduction of specific collagen cross-linking
levels.


  Many important problems in Bioinformatics (e.g., assembly or multi-assembly)
admit multiple solutions, while the final objective is to report only one. A
common approach to deal with this uncertainty is finding safe partial solutions
(e.g., contigs) which are common to all solutions. Previous research on safety
has focused on polynomially-time solvable problems, whereas many successful and
natural models are NP-hard to solve, leaving a lack of "safety tools" for such
problems. We propose the first method for computing all safe solutions for an
NP-hard problem, minimum flow decomposition. We obtain our results by
developing a "safety test" for paths based on a general Integer Linear
Programming (ILP) formulation. Moreover, we provide implementations with
practical optimizations aimed to reduce the total ILP time, the most efficient
of these being based on a recursive group-testing procedure.
  Results: Experimental results on the transcriptome datasets of Shao and
Kingsford (TCBB, 2017) show that all safe paths for minimum flow decompositions
correctly recover up to 90% of the full RNA transcripts, which is at least 25%
more than previously known safe paths, such as (Caceres et al. TCBB, 2021),
(Zheng et al., RECOMB 2021), (Khan et al., RECOMB 2022, ESA 2022). Moreover,
despite the NP-hardness of the problem, we can report all safe paths for 99.8%
of the over 27,000 non-trivial graphs of this dataset in only 1.5 hours. Our
results suggest that, on perfect data, there is less ambiguity than thought in
the notoriously hard RNA assembly problem.
  Availability: https://github.com/algbio/mfd-safety


  Predator-prey models have been shown to exhibit resonance-like behaviour, in
which random fluctuations in the number of organisms (demographic noise) are
amplified when their frequency is close to the natural oscillatory frequency of
the system. This behaviour has been traditionally studied in models with
exponentially distributed replication and death times. Here we consider a
biologically more realistic model, in which organisms replicate
quasi-synchronously such that the distribution of replication times has a
narrow maximum at some $T>0$ corresponding to the mean doubling time. We show
that when the frequency of replication $f=1/T$ is tuned to the natural
oscillatory frequency of the predator-prey model, the system exhibits
oscillations that are much stronger than in the model with Poissonian
(non-synchronous) replication and death. The effect can be explained by
resonant amplification of coloured noise generated by quasi-synchronous
replication events. To show this, we consider a single-species model with
quasi-synchronous replication. We calculate the spectrum and the amplitude of
demographic noise in this model, and use these results to obtain these
quantities for the two-species model.


  A core task in precision agriculture is the identification of climatic and
ecological conditions that are advantageous for a given crop. The most succinct
approach is geolocation, which is concerned with locating the native region of
a given sample based on its genetic makeup. Here, we investigate genomic
geolocation of Populus trichocarpa, or poplar, which has been identified by the
US Department of Energy as a fast-rotation biofuel crop to be harvested
nationwide. In particular, we approach geolocation from a reference-free
perspective, circumventing the need for compute-intensive processes such as
variant calling and alignment. Our model, MashNet, predicts latitude and
longitude for poplar trees from randomly-sampled, unaligned sequence fragments.
We show that our model performs comparably to Locator, a state-of-the-art
method based on aligned whole-genome sequence data. MashNet achieves an error
of 34.0 km^2 compared to Locator's 22.1 km^2. MashNet allows growers to quickly
and efficiently identify natural varieties that will be most productive in
their growth environment based on genotype. This paper explores geolocation for
precision agriculture while providing a framework and data source for further
development by the machine learning community.


  Pairs of similar compounds that only differ by a small structural
modification but exhibit a large difference in their binding affinity for a
given target are known as activity cliffs (ACs). It has been hypothesised that
quantitative structure-activity relationship (QSAR) models struggle to predict
ACs and that ACs thus form a major source of prediction error. However, a study
to explore the AC-prediction power of modern QSAR methods and its relationship
to general QSAR-prediction performance is lacking. We systematically construct
nine distinct QSAR models by combining three molecular representation methods
(extended-connectivity fingerprints, physicochemical-descriptor vectors and
graph isomorphism networks) with three regression techniques (random forests,
k-nearest neighbours and multilayer perceptrons); we then use each resulting
model to classify pairs of similar compounds as ACs or non-ACs and to predict
the activities of individual molecules in three case studies: dopamine receptor
D2, factor Xa, and SARS-CoV-2 main protease. We observe low AC-sensitivity
amongst the tested models when the activities of both compounds are unknown,
but a substantial increase in AC-sensitivity when the actual activity of one of
the compounds is given. Graph isomorphism features are found to be competitive
with or superior to classical molecular representations for AC-classification
and can thus be employed as baseline AC-prediction models or simple
compound-optimisation tools. For general QSAR-prediction, however,
extended-connectivity fingerprints still consistently deliver the best
performance. Our results provide strong support for the hypothesis that indeed
QSAR methods frequently fail to predict ACs. We propose twin-network training
for deep learning models as a potential future pathway to increase
AC-sensitivity and thus overall QSAR performance.


  Subtyping neuropsychiatric disorders like schizophrenia is essential for
improving the diagnosis and treatment of complex diseases. Subtyping
schizophrenia is challenging because it is polygenic and genetically
heterogeneous, rendering the standard symptom-based diagnosis often unreliable
and unrepeatable. We developed a novel network-based machine-learning approach,
netMoST, to subtyping psychiatric disorders. NetMoST identifies polygenic risk
SNP-allele modules from genome-wide genotyping data as polygenic haplotype
biomarkers (PHBs) for disease subtyping. We applied netMoST to subtype a cohort
of schizophrenia subjects into three distinct biotypes with differentiable
genetic, neuroimaging and functional characteristics. The PHBs of the first
biotype (36.9% of all patients) were related to neurodevelopment and cognition,
the PHBs of the second biotype (28.4%) were enriched for neuroimmune functions,
and the PHBs of the third biotype (34.7%) were associated with the transport of
calcium ions and neurotransmitters. Neuroimaging patterns provided additional
support to the new biotypes, with unique regional homogeneity (ReHo) patterns
observed in the brains of each biotype compared with healthy controls. Our
findings demonstrated netMoST's capability for uncovering novel biotypes of
complex diseases such as schizophrenia. The results also showed the power of
exploring polygenic allelic patterns that transcend the conventional GWAS
approaches.


  Antibody design is an essential yet challenging task in various domains like
therapeutics and biology. There are two major defects in current learning-based
methods: 1) tackling only a certain subtask of the whole antibody design
pipeline, making them suboptimal or resource-intensive. 2) omitting either the
framework regions or side chains, thus incapable of capturing the full-atom
geometry. To address these pitfalls, we propose dynamic Multi-channel
Equivariant grAph Network (dyMEAN), an end-to-end full-atom model for
E(3)-equivariant antibody design given the epitope and the incomplete sequence
of the antibody. Specifically, we first explore structural initialization as a
knowledgeable guess of the antibody structure and then propose shadow paratope
to bridge the epitope-antibody connections. Both 1D sequences and 3D structures
are updated via an adaptive multi-channel equivariant encoder that is able to
process protein residues of variable sizes when considering full atoms.
Finally, the updated antibody is docked to the epitope via the alignment of the
shadow paratope. Experiments on epitope-binding CDR-H3 design, complex
structure prediction, and affinity optimization demonstrate the superiority of
our end-to-end framework and full-atom modeling.


  Related groups of microbes are widely distributed across Earth's habitats,
implying numerous dispersal and adaptation events over evolutionary time.
However, to date, relatively little is known about the characteristics and
mechanisms of these habitat transitions, particularly for populations that
reside in animal microbiomes. Here, we review the existing literature
concerning habitat transitions among a variety of bacterial and archaeal
lineages, considering the frequency of migration events, potential
environmental barriers, and mechanisms of adaptation to new physicochemical
conditions, including the modification of protein inventories and other genomic
characteristics. Cells dependent on microbial hosts, particularly bacteria from
the Candidate Phyla Radiation (CPR), have undergone repeated habitat
transitions from environmental sources into animal microbiomes. We compare
their trajectories to those of both free-living cells - including the
Melainabacteria, Elusimicrobia, and methanogenic archaea - as well as cellular
endosymbionts and bacteriophages, which have made similar transitions. We
conclude by highlighting major related topics that may be worthy of future
study.


  Designing molecular structures with desired chemical properties is an
essential task in drug discovery and material design. However, finding
molecules with the optimized desired properties is still a challenging task due
to combinatorial explosion of candidate space of molecules. Here we propose a
novel \emph{decomposition-and-reassembling} based approach, which does not
include any optimization in hidden space and our generation process is highly
interpretable. Our method is a two-step procedure: In the first decomposition
step, we apply frequent subgraph mining to a molecular database to collect
smaller size of subgraphs as building blocks of molecules. In the second
reassembling step, we search desirable building blocks guided via reinforcement
learning and combine them to generate new molecules. Our experiments show that
not only can our method find better molecules in terms of two standard
criteria, the penalized $\log P$ and drug-likeness, but also generate drug
molecules with showing the valid intermediate molecules.


  The inception of spatial transcriptomics has allowed improved comprehension
of tissue architectures and the disentanglement of complex underlying
biological, physiological, and pathological processes through their positional
contexts. Recently, these contexts, and by extension the field, have seen much
promise and elucidation with the application of graph learning approaches. In
particular, neural operators have risen in regards to learning the mapping
between infinite-dimensional function spaces. With basic to deep neural network
architectures being data-driven, i.e. dependent on quality data for prediction,
neural operators provide robustness by offering generalization among different
resolutions despite low quality data. Graph neural operators are a variant that
utilize graph networks to learn this mapping between function spaces. The aim
of this research is to identify robust machine learning architectures that
integrate spatial information to predict tissue types. Under this notion, we
propose a study incorporating various graph neural network approaches to
validate the efficacy of applying neural operators towards prediction of brain
regions in mouse brain tissue samples as a proof of concept towards our
purpose. We were able to achieve an F1 score of nearly 72% for the graph neural
operator approach which outperformed all baseline and other graph network
approaches.


  Most packages for the analysis of fMRI-based functional connectivity (FC) and
genomic data are used with a programming language interface, lacking an
easy-to-navigate GUI frontend. This exacerbates two problems found in these
types of data: demographic confounds and quality control in the face of high
dimensionality of features. The reason is that it is too slow and cumbersome to
use a programming interface to create all the necessary visualizations required
to identify all correlations, confounding effects, or quality control problems
in a dataset. To remedy this situation, we have developed ImageNomer, a data
visualization and analysis tool that allows inspection of both subject-level
and cohort-level demographic, genomic, and imaging features. The software is
Python-based, runs in a self-contained Docker image, and contains a
browser-based GUI frontend. We demonstrate the usefulness of ImageNomer by
identifying an unexpected race confound when predicting achievement scores in
the Philadelphia Neurodevelopmental Cohort (PNC) dataset. In the past, many
studies have attempted to use FC to identify achievement-related features in
fMRI. Using ImageNomer, we find a clear potential for confounding effects of
race. Using correlation analysis in the ImageNomer software, we show that FCs
correlated with Wide Range Achievement Test (WRAT) score are in fact more
highly correlated with race. Investigating further, we find that whereas both
FC and SNP (genomic) features can account for 10-15\% of WRAT score variation,
this predictive ability disappears when controlling for race. In this work, we
demonstrate the advantage of our ImageNomer GUI tool in data exploration and
confound detection. Additionally, this work identifies race as a strong
confound in FC data and casts doubt on the possibility of finding unbiased
achievement-related features in fMRI and SNP data of healthy adolescents.


  Exact solutions of the SEIR epidemic model are derived, and various
properties of solutions are obtained directly from the exact solution. In this
paper Abel differential equations play an important role in establishing the
exact solution of SEIR differential system, in particular the number of
infected individuals can be represented in a simple form by using a positive
solution of an Abel differential equation. It is shown that the parametric form
of the exact solution satisfies some linear differential system including a
positive solution of an Abel differential equation.


  Molecular dynamic simulations are important in computational physics,
chemistry, material, and biology. Machine learning-based methods have shown
strong abilities in predicting molecular energy and properties and are much
faster than DFT calculations. Molecular energy is at least related to atoms,
bonds, bond angles, torsion angles, and nonbonding atom pairs. Previous
Transformer models only use atoms as inputs which lack explicit modeling of the
aforementioned factors. To alleviate this limitation, we propose Moleformer, a
novel Transformer architecture that takes nodes (atoms) and edges (bonds and
nonbonding atom pairs) as inputs and models the interactions among them using
rotational and translational invariant geometry-aware spatial encoding.
Proposed spatial encoding calculates relative position information including
distances and angles among nodes and edges. We benchmark Moleformer on OC20 and
QM9 datasets, and our model achieves state-of-the-art on the initial state to
relaxed energy prediction of OC20 and is very competitive in QM9 on predicting
quantum chemical properties compared to other Transformer and Graph Neural
Network (GNN) methods which proves the effectiveness of the proposed
geometry-aware spatial encoding in Moleformer.


  Glutamate-gated kainate receptors (KARs) are ubiquitous in the central
nervous system of vertebrates, mediate synaptic transmission on post-synapse,
and modulate transmitter release on pre-synapse. In the brain, the trafficking,
gating kinetics, and pharmacology of KARs are tightly regulated by Neuropilin
and tolloid-like proteins (Netos). Here we report cryo-EM structures of
homo-tetrameric GluK2 in complex with Neto2 at inhibited and desensitized
states, illustrating variable stoichiometry of GluK2-Neto2 complexes, with one
or two Neto2 subunits associate with the GluK2. We find that Neto2 accesses
only two broad faces of KARs, intermolecularly crosslinking the lower-lobe of
ATDA/C, upper-lobe of LBDB/D, and lower-lobe of LBDA/C, illustrating how Neto2
regulates receptor-gating kinetics. The transmembrane helix of Neto2 is
positioned proximal to the selectivity filter and competes with the amphiphilic
H1-helix after M4 for interacting with an ICD formed by the M1-M2 linkers of
the receptor, revealing how rectification is regulated by Neto2.


  Frequent emergence of communicable diseases has been a major concern
worldwide. Lack of sufficient resources to mitigate the disease-burden makes
the situation even more challenging for lower-income countries. Hence, strategy
development towards disease eradication and optimal management of the social
and economic burden has garnered a lot of attention in recent years. In this
context, we quantify the optimal fraction of resources that can be allocated to
two major intervention measures, namely reduction of disease transmission and
improvement of healthcare infrastructure. Our results demonstrate that the
effectiveness of each of the interventions has a significant impact on the
optimal resource allocation in both long-term disease dynamics and outbreak
scenarios. Often allocating resources to both strategies is optimal. For
long-term dynamics, a non-monotonic behavior of optimal resource allocation
with intervention effectiveness is observed which is different from the more
intuitive strategy recommended in case of outbreaks. Further, our result
indicates that the relationship between investment into interventions and the
corresponding outcomes has a decisive role in determining optimal strategies.
Intervention programs with decreasing returns promote the necessity for
resource sharing. Our study provides a fundamental insight into determining the
best response strategy in case of controlling epidemics under
resource-constrained situations.


  We investigate sub-leading orders of the classic SEIR-model using contact
matrices from modeling of the Omicron and Delta variants of COVID-19 in
Denmark. The goal of this is to illustrate when the growth rate, and by
extension the infectiousness, can be accurately measured in a new outbreak,
e.g. after introduction of a new variant of a virus. We find that as long as
susceptible depletion is a minor effect, the transients are gone within around
4 generations.


  This paper will argue that one of the biggest challenges for livestock
genomics is to make whole-genome sequencing and functional genomics applicable
to breeding practice. It discusses potential explanations for why it is so
difficult to consistently improve the accuracy of genomic prediction by means
of whole-genome sequence data, and three potential attacks on the problem.


  Maintaining stability in an uncertain environment is essential for proper
functioning of living systems. Robust perfect adaptation (RPA) is a property of
a system that generates an output at a fixed level even after fluctuations in
input stimulus without fine-tuning parameters, and it is important to
understand how this feature is implemented through biochemical networks. The
existing literature has mainly focused on RPA of the concentration of a chosen
chemical species, and no generic analysis has been made on RPA of reaction
fluxes, that play an equally important role. Here, we identify structural
conditions on reaction networks under which all the reaction fluxes exhibit RPA
against the perturbation of the parameters inside a subnetwork. Based on this
understanding, we give a recipe for obtaining a simpler reaction network, from
which we can fully recover the steady-state reaction fluxes of the original
system. This helps us identify key parameters that determine the fluxes and
study the properties of complex reaction networks using a smaller one without
losing any information about steady-state reaction fluxes.