  The large-scale evolution of the SARS-CoV-2 virus has been marked by rapid
turnover of genetic clades. New variants show intrinsic changes, notably
increased transmissibility, as well as antigenic changes that reduce the
cross-immunity induced by previous infections or vaccinations. How this
functional variation shapes the global evolutionary dynamics has remained
unclear. Here we show that selection induced by vaccination impacts on the
recent antigenic evolution of SARS-CoV-2; other relevant forces include
intrinsic selection and antigenic selection induced by previous infections. We
obtain these results from a fitness model with intrinsic and antigenic fitness
components. To infer model parameters, we combine time-resolved sequence data,
epidemiological records, and cross-neutralisation assays. This model accurately
captures the large-scale evolutionary dynamics of SARS-CoV-2 in multiple
geographical regions. In particular, it quantifies how recent vaccinations and
infections affect the speed of frequency shifts between viral variants. Our
results show that timely neutralisation data can be harvested to identify
hotspots of antigenic selection and to predict the impact of vaccination on
viral evolution.


  Diagnostic testing followed by isolation of identified cases with subsequent
tracing and quarantine of close contacts - often referred to as
test-trace-isolate-and-quarantine (TTIQ) strategy - is one of the cornerstone
measures of infectious disease control. The COVID-19 pandemic has highlighted
that an appropriate response to outbreaks requires us to be aware about the
effectiveness of such containment strategies. This can be evaluated using
mathematical models. We present a delay differential equation model of TTIQ
interventions for infectious disease control. Our model incorporates a detailed
mechanistic description of the state-dependent dynamics induced by limited TTIQ
capacities. In addition, we account for transmission during the early phase of
SARS-CoV-2 infection, including presymptomatic transmission, which may be
particularly adverse to a TTIQ based control. Numerical experiments, inspired
by the early spread of COVID-19 in Germany, reveal the effectiveness of TTIQ in
a scenario where immunity within the population is low and pharmaceutical
interventions are absent - representative of a typical situation during the
(re-)emergence of infectious diseases for which therapeutic drugs or vaccines
are not yet available. Stability and sensitivity analyses emphasize factors,
partially related to the specific disease, which impede or enhance the success
of TTIQ. Studying the diminishing effectiveness of TTIQ along simulations of an
epidemic wave we highlight consequences for intervention strategies.


  Three-dimensional volumetric imaging of cells allows for in situ
visualization, thus preserving contextual insights into cellular processes.
Despite recent advances in machine learning methods, morphological analysis of
sub-nuclear structures have proven challenging due to both the shallow contrast
profile and the technical limitation in feature detection. Here, we present a
convolutional neural network, supervised deep learning-based approach which can
identify sub-nuclear structures with 90% accuracy. We develop and apply this
model to C. elegans gonads imaged using focused ion beam milling combined with
scanning electron microscopy resulting in the accurate identification and
segmentation of all sub-nuclear structures including entire chromosomes. We
discuss in depth the architecture, parameterization, and optimization of the
deep learning model, as well as provide evaluation metrics to assess the
quality of the network prediction. Lastly, we highlight specific aspects of the
model that can be optimized for its broad application to other volumetric
imaging data as well as in situ cryo-electron tomography.


  The surprisingly mercurial Covid-19 pandemic has highlighted the need to not
only accelerate research on infectious disease, but to also study them using
novel techniques and perspectives. A major contributor to the difficulty of
containing the current pandemic is due to the highly asymptomatic nature of the
disease. In this investigation, we develop a modeling framework to study the
spatio-temporal evolution of diseases with high rates of asymptomatic
transmission, and we apply this framework to a hypothetical country with
mathematically tractable geography; namely, square counties uniformly organized
into a rectangle. We first derive a model for the temporal dynamics of
susceptible, infected, and recovered populations, which is applied at the
county level. Next we use likelihood-based parameter estimation to derive
temporally varying disease transmission parameters on the state-wide level.
While these two methods give us some spatial structure and show the effects of
behavioral and policy changes, they miss the evolution of hot zones that have
caused significant difficulties in resource allocation during the current
pandemic. It is evident that the distribution of cases will not be stagnantly
based on the population density, as with many other diseases, but will
continuously evolve. We model this as a diffusive process where the diffusivity
is spatially varying based on the population distribution, and temporally
varying based on the current number of simulated asymptomatic cases. With this
final addition coupled to the SIR model with temporally varying transmission
parameters, we capture the evolution of "hot zones" in our hypothetical setup.


  Mechanistic simulators are an indispensable tool for epidemiology to explore
the behavior of complex, dynamic infections under varying conditions and
navigate uncertain environments. Agent-based models (ABMs) are an increasingly
popular simulation paradigm that can represent the heterogeneity of contact
interactions with granular detail and agency of individual behavior. However,
conventional ABM frameworks are not differentiable and present challenges in
scalability; due to which it is non-trivial to connect them to auxiliary data
sources. In this paper, we introduce GradABM: a scalable, differentiable design
for agent-based modeling that is amenable to gradient-based learning with
automatic differentiation. GradABM can quickly simulate million-size
populations in few seconds on commodity hardware, integrate with deep neural
networks and ingest heterogeneous data sources. This provides an array of
practical benefits for calibration, forecasting, and evaluating policy
interventions. We demonstrate the efficacy of GradABM via extensive experiments
with real COVID-19 and influenza datasets.


  The identification of the transmission parameters of a virus is fundamental
to identify the optimal public health strategy. These parameters can present
significant changes over time caused by genetic mutations or viral
recombination, making their continuous monitoring fundamental. Here we present
a method, suitable for this task, which uses as unique information the daily
number of reported cases. The method is based on a time since infection model
where transmission parameters are obtained by means of an efficient
maximization procedure of the likelihood. Applying the method to SARS-CoV-2
data in Italy we find an average generation time $z = 3.2 \pm 0.8$ days, during
the temporal window when the majority of infections can be attributed to the
Omicron variants. At the same time we find a significantly larger value $z =
6.2 \pm 1.1$ days, in the temporal window when spreading was dominated by the
Delta variant. We are also able to show that the presence of the Omicron
variant, characterized by a shorter $z$, was already detectable in the first
weeks of December 2021, in full agreement with results provided by sequences of
SARS-CoV-2 genomes reported in national databases. Our results therefore
indicate that the novel approach can indicate the existence of virus variants
resulting particularly useful in situations when information about genomic
sequencing is not yet available.


  Profile hidden Markov models (pHMMs) are widely employed in various
bioinformatics applications to identify similarities between biological
sequences, such as DNA or protein sequences. In pHMMs, sequences are
represented as graph structures. These probabilities are subsequently used to
compute the similarity score between a sequence and a pHMM graph. The
Baum-Welch algorithm, a prevalent and highly accurate method, utilizes these
probabilities to optimize and compute similarity scores. However, the
Baum-Welch algorithm is computationally intensive, and existing solutions offer
either software-only or hardware-only approaches with fixed pHMM designs. We
identify an urgent need for a flexible, high-performance, and energy-efficient
HW/SW co-design to address the major inefficiencies in the Baum-Welch algorithm
for pHMMs.
  We introduce ApHMM, the first flexible acceleration framework designed to
significantly reduce both computational and energy overheads associated with
the Baum-Welch algorithm for pHMMs. ApHMM tackles the major inefficiencies in
the Baum-Welch algorithm by 1) designing flexible hardware to accommodate
various pHMM designs, 2) exploiting predictable data dependency patterns
through on-chip memory with memoization techniques, 3) rapidly filtering out
negligible computations using a hardware-based filter, and 4) minimizing
redundant computations.
  ApHMM achieves substantial speedups of 15.55x - 260.03x, 1.83x - 5.34x, and
27.97x when compared to CPU, GPU, and FPGA implementations of the Baum-Welch
algorithm, respectively. ApHMM outperforms state-of-the-art CPU implementations
in three key bioinformatics applications: 1) error correction, 2) protein
family search, and 3) multiple sequence alignment, by 1.29x - 59.94x, 1.03x -
1.75x, and 1.03x - 1.95x, respectively, while improving their energy efficiency
by 64.24x - 115.46x, 1.75x, 1.96x.


  A pandemic, the worldwide spread of a disease, can threaten human beings from
the social as well as biological perspectives and paralyze existing living
habits. To stave off the more devastating disaster and return to a normal life,
people make tremendous efforts at multiscale levels from individual to
worldwide: paying attention to hand hygiene, developing social policies such as
wearing masks, social distancing, quarantine, and inventing vaccines and
remedy. Regarding the current severe pandemic, namely the coronavirus disease
2019, we explore the spreading-suppression effect when adopting the
aforementioned efforts. Especially the quarantine and vaccination are
considered since they are representative primary treatments for block spreading
and prevention at the government level. We establish a compartment model
consisting of susceptible (S), vaccination (V), exposed (E), infected (I),
quarantined (Q), and recovered (R) compartments, called SVEIQR model. We look
into the infected cases in Seoul and consider three kinds of vaccines, Pfizer,
Moderna, and AstraZeneca. The values of the relevant parameters are obtained
from empirical data from Seoul and clinical data for vaccines and estimated by
Bayesian inference. After confirming that our SVEIQR model is plausible, we
test the various scenarios by adjusting the associated parameters with the
quarantine and vaccination policies around the current values. The quantitative
result obtained from our model could suggest a guideline for policy making on
effective vaccination and social policies.


  Intracellular transport processes are essential to the healthy development of
many organisms as well as more generally to healthy cellular function. The
complex dynamics and interactions between protein molecules and filaments on
different time and spatial scales generate many opportunities for mathematical
modeling and analysis that can provide insights into protein sorting and
overall cellular organization. Systems of advection-reaction-diffusion partial
differential equations and stochastic state-switching particle models are
especially useful in understanding spatiotemporal protein dynamics inside
cells.


  Encouraged by decision makers' appetite for future information on topics
ranging from elections to pandemics, and enabled by the explosion of data and
computational methods, model based forecasts have garnered increasing influence
on a breadth of decisions in modern society. Using several classic examples
from fisheries management, I demonstrate that selecting the model or models
that produce the most accurate and precise forecast (measured by statistical
scores) can sometimes lead to worse outcomes (measured by real-world
objectives). This can create a forecast trap, in which the outcomes such as
fish biomass or economic yield decline while the manager becomes increasingly
convinced that these actions are consistent with the best models and data
available. The forecast trap is not unique to this example, but a fundamental
consequence of non-uniqueness of models. Existing practices promoting a broader
set of models are the best way to avoid the trap.


  Carbon fixation is the process by which CO2 is converted from a gas into
biomass. The Calvin Benson Bassham (CBB) cycle is the dominant carbon fixation
pathway on earth, driving >99.5% of the ~120 billion tons of carbon that are
"fixed" as sugar, by plants, algae and cyanobacteria. The carboxylase enzyme in
the CBB, ribulose-1,5-bisphosphate carboxylase/oxygenase (rubisco), fixes one
CO2 molecule per turn of the cycle. Despite being critical to the assimilation
of carbon, rubisco's kinetic rate is not very fast and it is a bottleneck in
flux through the pathway. This presents a paradox - why hasn't rubisco evolved
to be a better catalyst? Many hypothesize that the catalytic mechanism of
rubisco is subject to one or more trade-offs, and that rubisco variants have
been optimized for their native physiological environment. Here we review the
evolution and biochemistry of rubisco through the lens of structure and
mechanism in order to understand what trade-offs limit its improvement. We also
review the many attempts to improve rubisco itself and, thereby, promote plant
growth.


  While there is a general focus on predictions of values, mathematically more
appropriate is prediction of probability distributions: with additional
possibilities like prediction of uncertainty, higher moments and quantiles. For
the purpose of the computer-aided drug design field, this article applies
Hierarchical Correlation Reconstruction approach, previously applied in the
analysis of demographic, financial and astronomical data. Instead of a single
linear regression to predict values, it uses multiple linear regressions to
independently predict multiple moments, finally combining them into predicted
probability distribution, here of several ADMET properties based on
substructural fingerprint developed by Klekota\&Roth. Discussed application
example is inexpensive selection of a percentage of molecules with properties
nearly certain to be in a predicted or chosen range during virtual screening.
Such an approach can facilitate the interpretation of the results as the
predictions characterized by high rate of uncertainty are automatically
detected. In addition, for each of the investigated predictive problems, we
detected crucial structural features, which should be carefully considered when
optimizing compounds towards particular property. The whole methodology
developed in the study constitutes therefore a great support for medicinal
chemists, as it enable fast rejection of compounds with the lowest potential of
desired physicochemical/ADMET characteristic and guides the compound
optimization process.


  The paper deals with the setting where two viruses (say virus 1 and virus 2)
coexist in a population, and they are not necessarily mutually exclusive, in
the sense that infection due to one virus does not preclude the possibility of
simultaneous infection due to the other. We develop a coupled bi-virus
susceptible-infected-susceptible (SIS) model from a 4n-state Markov chain
model, where n is the number of agents (i.e., individuals or subpopulation) in
the population. We identify a sufficient condition for both viruses to
eventually die out, and a sufficient condition for the existence, uniqueness
and asymptotic stability of the endemic equilibrium of each virus. We establish
a sufficient condition and multiple necessary conditions for local exponential
convergence to the boundary equilibrium (i.e., one virus persists, the other
one dies out) of each virus. Under mild assumptions on the healing rate, we
show that there cannot exist a coexisting equilibrium where for each node there
is a nonzero fraction infected only by virus 1; a nonzero fraction infected
only by virus 2; but no fraction that is infected by both viruses 1 and 2.
Likewise, assuming that healing rates are strictly positive, a coexisting
equilibrium where for each node there is a nonzero fraction infected by both
viruses 1 and 2, but no fraction is infected only by virus 1 (resp. virus 2)
does not exist. Further, we provide a necessary condition for the existence of
certain other kinds of coexisting equilibria. We show that, unlike the
competitive bivirus model, the coupled bivirus model is not monotone. Finally,
we illustrate our theoretical findings using an extensive set of in-depth
simulations.


  Gonadotrophin-releasing hormone receptor (GnRH1R) is a promising therapeutic
target for the treatment of uterine diseases. To date, several GnRH1R
antagonists are available in clinical investigation without satisfying multiple
property constraints. To fill this gap, we aim to develop a deep learning-based
framework to facilitate the effective and efficient discovery of a new orally
active small-molecule drug targeting GnRH1R with desirable properties. In the
present work, a ligand-and-structure combined model, namely LS-MolGen, was
firstly proposed for molecular generation by fully utilizing the information on
the known active compounds and the structure of the target protein, which was
demonstrated by its superior performance than ligand- or structure-based
methods separately. Then, a in silico screening including activity prediction,
ADMET evaluation, molecular docking and FEP calculation was conducted, where
~30,000 generated novel molecules were narrowed down to 8 for experimental
synthesis and validation. In vitro and in vivo experiments showed that three of
them exhibited potent inhibition activities (compound 5 IC50 = 0.856 nM,
compound 6 IC50 = 0.901 nM, compound 7 IC50 = 2.54 nM) against GnRH1R, and
compound 5 performed well in fundamental PK properties, such as half-life, oral
bioavailability, and PPB, etc. We believed that the proposed
ligand-and-structure combined molecular generative model and the whole
computer-aided workflow can potentially be extended to similar tasks for de
novo drug design or lead optimization.


  A common sample descriptor in human genomics studies is that of 'genetic
ancestry group', with terms such as 'European genetic ancestry' or 'East Asian
genetic ancestry' frequently used in publications to describe the genetics of
groups of individuals based on the analysis of their genotypes. In this
Perspective, I argue that these terms are imprecise and potentially misleading
and that, for most applications, simple statements of genetic similarity
represent a more accurate description.


  Calcium oscillation is an important calcium homeostasis, imbalance of which
is the key mechanism of initiation and progression of many major diseases. The
formation and maintenance of calcium homeostasis are closely related to the
spatial distribution of calcium channels. In the current paper, a theoretical
framework is established by abstracting the spatial distribution of the calcium
channels as a nonlinear biological complex network with calcium channels as
nodes and Ca$^{2+}$ as edges. A dynamical model for a RyR is adopted to
investigate the effect of spatial distribution on calcium oscillation. The
mean-field model can be well reproduced from the complete graph and dense
Erd\"os-R\'enyi network. The synchronization of RyRs is found important to
generate a global calcium oscillation. The clique graph with a cluster
structure can not produce a global oscillation due to the failure of
synchronization between clusters. A more realistic geometric network is
constructed in a two-dimensional plane based on the experimental information
about the RyR arrangement of clusters and the frequency distribution of cluster
sizes. Different from the clique graph, the global oscillation can be generated
with reasonable parameters on the geometric network. The simulation also
suggests that existence of small clusters and rogue RyR's plays an important
role in the maintenance of global calcium oscillation through keeping
synchronization between large clusters. Such results support the heterogeneous
distribution of RyR's with different-size clusters, which is helpful to
understand recent observations with super resolution nanoscale imaging
techniques. The current theoretical framework can also be extent to investigate
other phenomena in calcium signal transduction.


  This paper proposes a data-driven approximate Bayesian computation framework
for parameter estimation and uncertainty quantification of epidemic models,
which incorporates two novelties: (i) the identification of the initial
conditions by using plausible dynamic states that are compatible with
observational data; (ii) learning of an informative prior distribution for the
model parameters via the cross-entropy method. The new methodology's
effectiveness is illustrated with the aid of actual data from the COVID-19
epidemic in Rio de Janeiro city in Brazil, employing an ordinary differential
equation-based model with a generalized SEIR mechanistic structure that
includes time-dependent transmission rate, asymptomatics, and hospitalizations.
A minimization problem with two cost terms (number of hospitalizations and
deaths) is formulated, and twelve parameters are identified. The calibrated
model provides a consistent description of the available data, able to
extrapolate forecasts over a few weeks, making the proposed methodology very
appealing for real-time epidemic modeling.


  A major goal in genomics is to properly capture the complex dynamical
behaviors of gene regulatory networks (GRNs). This includes inferring the
complex interactions between genes, which can be used for a wide range of
genomics analyses, including diagnosis or prognosis of diseases and finding
effective treatments for chronic diseases such as cancer. Boolean networks have
emerged as a successful class of models for capturing the behavior of GRNs. In
most practical settings, inference of GRNs should be achieved through limited
and temporally sparse genomics data. A large number of genes in GRNs leads to a
large possible topology candidate space, which often cannot be exhaustively
searched due to the limitation in computational resources. This paper develops
a scalable and efficient topology inference for GRNs using Bayesian
optimization and kernel-based methods. Rather than an exhaustive search over
possible topologies, the proposed method constructs a Gaussian Process (GP)
with a topology-inspired kernel function to account for correlation in the
likelihood function. Then, using the posterior distribution of the GP model,
the Bayesian optimization efficiently searches for the topology with the
highest likelihood value by optimally balancing between exploration and
exploitation. The performance of the proposed method is demonstrated through
comprehensive numerical experiments using a well-known mammalian cell-cycle
network.


  Polygnicity refers to the phenomenon that multiple genetic variants have a
non-zero effect on a complex trait. It is defined as the proportion of genetic
variants that have a nonzero effect on the trait. Evaluation of polygenicity
can provide valuable insights into the genetic architecture of the trait.
Several recent works have attempted to estimate polygenicity at the SNP level.
However, evaluating polygenicity at the gene level can be biologically more
meaningful. We propose the notion of gene-level polygenicity, defined as the
proportion of genes having a non-zero effect on the trait under the framework
of transcriptome-wide association study. We introduce a Bayesian approach
polygene to estimate this quantity for a trait. The method is based on spike
and slab prior and simultaneously provides an optimal subset of non-null genes.
Our simulation study shows that polygene efficiently estimates gene-level
polygenicity. The method produces downward bias for small choices of trait
heritability due to a non-null gene, which diminishes rapidly with an increase
in the GWAS sample size. While identifying the optimal subset of non-null
genes, polygene offers a high level of specificity and an overall good level of
sensitivity -- the sensitivity increases as the sample size of the reference
panel expression and GWAS data increase. We applied the method to seven
phenotypes in the UK Biobank, integrating expression data. We find height to be
most polygenic and asthma to be the least polygenic. Our analysis suggests that
both HDL and triglycerides are more polygenic than LDL.


  Gene set collections are a common ground to study the enrichment of genes for
specific phenotypic traits. Gene set enrichment analysis aims to identify genes
that are over-represented in gene sets collections and might be associated with
a specific phenotypic trait. However, as this involves a massive number of
hypothesis testing, it is often questionable whether a pre-processing step to
reduce gene sets collections' sizes is helpful. Moreover, the often highly
overlapping gene sets and the consequent low interpretability of gene sets'
collections demand for a reduction of the included gene sets. Inspired by this
bioinformatics context, we propose a method to rank sets within a family of
sets based on the distribution of the singletons and their size. We obtain
sets' importance scores by computing Shapley values without incurring into the
usual exponential number of evaluations of the value function. Moreover, we
address the challenge of including a redundancy awareness in the rankings
obtained where, in our case, sets are redundant if they show prominent
intersections. We finally evaluate our approach for gene sets collections; the
rankings obtained show low redundancy and high coverage of the genes. The
unsupervised nature of the proposed ranking does not allow for an evident
increase in the number of significant gene sets for specific phenotypic traits
when reducing the size of the collections. However, we believe that the
rankings proposed are of use in bioinformatics to increase interpretability of
the gene sets collections and a step forward to include redundancy into Shapley
values computations.


  In this paper, we investigate the existence of stability-changing
bifurcations in epidemiological models used to study the spread of zombiism
through a human population. These bifurcations show that although linear
instability of disease-free equilibria may exist in a model, perturbations of
model parameters may result in stability. Thus, we show that humans can survive
a zombie outbreak.


  There has been a growing acknowledgement of the involvement of the gut
microbiome - the collection of microbes that reside in our gut - in regulating
our mood and behaviour. This phenomenon is referred to as the
microbiome-gut-brain axis. While our techniques to measure the presence and
abundance of these microbes have been steadily improving, the analysis of
microbiome data is non-trivial.
  Here, we present a perspective on the concepts and foundations of data
analysis and interpretation of microbiome experiments with a focus on the
microbiome-gut-brain axis domain. We give an overview of foundational
considerations prior to commencing analysis alongside the core microbiome
analysis approaches of alpha diversity, beta diversity, differential feature
abundance and functional inference. We emphasize the compositional data
analysis (CoDA) paradigm.
  Further, this perspective features an extensive and heavily annotated
microbiome analysis in R in the supplementary materials, as a resource for new
and experienced bioinformaticians alike.


  The approval success rate of drug candidates is very low with the majority of
failure due to safety and efficacy. Increasingly available high dimensional
information on targets, drug molecules and indications provides an opportunity
for ML methods to integrate multiple data modalities and better predict
clinically promising drug targets. Notably, drug targets with human genetics
evidence are shown to have better odds to succeed. However, a recent tensor
factorization-based approach found that additional information on targets and
indications might not necessarily improve the predictive accuracy. Here we
revisit this approach by integrating different types of human genetics evidence
collated from publicly available sources to support each target-indication
pair. We use Bayesian tensor factorization to show that models incorporating
all available human genetics evidence (rare disease, gene burden, common
disease) modestly improves the clinical outcome prediction over models using
single line of genetics evidence. We provide additional insight into the
relative predictive power of different types of human genetics evidence for
predicting the success of clinical outcomes.


  Ancestral sequence reconstruction is a key task in computational biology. It
consists in inferring a molecular sequence at an ancestral species of a known
phylogeny, given descendant sequences at the tip of the tree. In addition to
its many biological applications, it has played a key role in elucidating the
statistical performance of phylogeny estimation methods. Here we establish a
formal connection to another important bioinformatics problem, multiple
sequence alignment, where one attempts to best align a collection of molecular
sequences under some mismatch penalty score by inserting gaps. Our result is
counter-intuitive: we show that perfect pairwise sequence alignment with high
probability is possible in principle at arbitrary large evolutionary distances
- provided the phylogeny is known and dense enough. We use techniques from
ancestral sequence reconstruction in the taxon-rich setting together with the
probabilistic analysis of sequence evolution models involving insertions and
deletions.


  The problem of how to genetically modify cells in order to maximize a certain
cellular phenotype has taken center stage in drug development over the last few
years (with, for example, genetically edited CAR-T, CAR-NK, and CAR-NKT cells
entering cancer clinical trials). Exhausting the search space for all possible
genetic edits (perturbations) or combinations thereof is infeasible due to cost
and experimental limitations. This work provides a theoretically sound
framework for iteratively exploring the space of perturbations in pooled
batches in order to maximize a target phenotype under an experimental budget.
Inspired by this application domain, we study the problem of batch query bandit
optimization and introduce the Optimistic Arm Elimination ($\mathrm{OAE}$)
principle designed to find an almost optimal arm under different functional
relationships between the queries (arms) and the outputs (rewards). We analyze
the convergence properties of $\mathrm{OAE}$ by relating it to the Eluder
dimension of the algorithm's function class and validate that $\mathrm{OAE}$
outperforms other strategies in finding optimal actions in experiments on
simulated problems, public datasets well-studied in bandit contexts, and in
genetic perturbation datasets when the regression model is a deep neural
network. OAE also outperforms the benchmark algorithms in 3 of 4 datasets in
the GeneDisco experimental planning challenge.


  Chronic infections of the human immunodeficiency virus (HIV) create a very
complex co-evolutionary process, where the virus tries to escape the
continuously adapting host immune system. Quantitative details of this process
are largely unknown and could help in disease treatment and vaccine
development. Here we study a longitudinal dataset of ten HIV-infected people,
where both the B-cell receptors and the virus are deeply sequenced. We focus on
simple measures of turnover, which quantify how much the composition of the
viral strains and the immune repertoire change between time points. At the
single-patient level, the viral-host turnover rates do not show any
statistically significant correlation, however they correlate if the
information is aggregated across patients. In particular, we identify an
anti-correlation: large changes in the viral pool composition come with small
changes in the B-cell receptor repertoire. This result seems to contradict the
naive expectation that when the virus mutates quickly, the immune repertoire
needs to change to keep up. However, we show that the observed anti-correlation
naturally emerges and can be understood in terms of simple population-genetics
models.


  Reconstructing the ancestral state of a group of species helps answer many
important questions in evolutionary biology. Therefore, it is crucial to
understand when we can estimate the ancestral state accurately. Previous works
provide a necessary and sufficient condition, called the big bang condition,
for the existence of an accurate reconstruction method under discrete trait
evolution models and the Brownian motion model. In this paper, we extend this
result to a wide range of continuous trait evolution models. In particular, we
consider a general setting where continuous traits evolve along the tree
according to stochastic processes that satisfy some regularity conditions. We
verify these conditions for popular continuous trait evolution models including
Ornstein-Uhlenbeck, reflected Brownian Motion, and Cox-Ingersoll-Ross.


  With the characterisations of potentially habitable planetary atmospheres on
the horizon, the search for biosignatures is set to become a major area of
research in the coming decades. To understand the atmospheric characteristics
that might indicate alien life we must understand the abiotic characteristics
of a planet and how life interacts with its environment. In the field of
biogeochemistry, sophisticated models of life-environment coupled systems
demonstrate that many assumptions specific to Earth-based life, e.g. specific
ATP maintenance costs, are unnecessary to accurately model a biosphere. We
explore a simple model of a single-species microbial biosphere that produces
CH4 as a byproduct of the microbes' energy extraction - known as a type I
biosignature. We demonstrate that although significantly changing the
biological parameters has a large impact on the biosphere's total population,
such changes have only a minimal impact on the strength of the resulting
biosignature, while the biosphere is limited by H2 availability. We extend the
model to include more accurate microbial energy harvesting and show that
adjusting microbe parameters can lead to a regime change where the biosphere
becomes limited by energy availability and no longer fully exploits the
available H2, impacting the strength of the resulting biosignature. We
demonstrate that, for a nutrient limited biosphere, identifying the limiting
nutrient, understanding the abiotic processes that control its abundance, and
determining the biosphere's ability to exploit it, are more fundamental for
making type I biosignature predictions than the details of the population
dynamics of the biosphere.


  We propose a new mathematical model to study the coinfection dynamics of
COVID-19 and bacterial pneumonia. Our model includes two infection ways for
pneumonia, corresponding to community-acquired and hospital-acquired
infections. We show that the existence and local stability of equilibria depend
on three different parameters, which are interpreted as the basic reproduction
numbers of COVID-19, bacterial pneumonia, and bacterial population. Numerical
simulations are performed to complement our theoretical analysis, and we show
that both diseases can persist if the basic reproduction number of COVID-19 is
greater than one.


  Boolean networks (BNs) are discrete dynamical systems with applications to
the modeling of cellular behaviors. In this paper, we demonstrate how the
software BoNesis can be employed to exhaustively identify combinations of
perturbations which enforce properties on their fixed points and attractors. We
consider marker properties, which specify that some components are fixed to a
specific value. We study 4 variants of the marker reprogramming problem: the
reprogramming of fixed points, of minimal trap spaces, and of fixed points and
minimal trap spaces reachable from a given initial configuration with the most
permissive update mode. The perturbations consist of fixing a set of components
to a fixed value. They can destroy and create new attractors. In each case, we
give an upper bound on their theoretical computational complexity, and give an
implementation of the resolution using the BoNesis Python framework. Finally,
we lift the reprogramming problems to ensembles of BNs, as supported by
BoNesis, bringing insight on possible and universal reprogramming strategies.
This paper can be executed and modified interactively.


  This contribution focuses on the fascinating RNA molecule, its
sequence-dependent folding driven by base-pairing interactions, the interplay
between these interactions and natural evolution, and its multiple regulatory
roles. The four of us have dug into these topics using the tools and the spirit
of the statistical physics of disordered systems, and in particular the concept
of a disordered (energy/fitness) landscape. After an introduction to RNA
molecules and the perspectives they open not only in evolutionary and synthetic
biology but also in medicine, we will introduce the important notions of energy
and fitness landscapes for these molecules. In Section III we will review some
models and algorithms for RNA sequence-to-secondary-structure mapping. Section
IV discusses how the secondary-structure energy landscape can be derived from
unzipping data. Section V deals with the inference of RNA structure from
evolutionary sequence data sampled in different organisms. This will shift the
focus from the `sequence-to-structure' mapping described in Section III to a
`sequence-to-function' landscape that can be inferred from laboratory
evolutionary data on DNA aptamers. Finally, in Section VI, we shall discuss the
rich theoretical picture linking networks of interacting RNA molecules to the
organization of robust, systemic regulatory programs. Along this path, we will
therefore explore phenomena across multiple scales in space, number of
molecules and time, showing how the biological complexity of the RNA world can
be captured by the unifying concepts of statistical physics.


  We analyzed the first 255 PCR-confirmed cases of monkeypox occurred in Italy
in 2022. Preliminary estimates indicate: mean incubation period of 9.1 days
(95%CI of the mean: 6.5-10.9); mean generation time of 12.5 days (95%CI of the
mean: 7.5-17.3); reproduction number in the MSM community of 2.43 (95%CI:
1.82-3.26).


  A dramatic increase in the number of outbreaks of Dengue has recently been
reported, and climate change is likely to extend the geographical spread of the
disease. In this context, this paper shows how a neural network approach can
incorporate Dengue and COVID-19 data as well as external factors (such as
social behaviour or climate variables), to develop predictive models that could
improve our knowledge and provide useful tools for health policy makers.
Through the use of neural networks with different social and natural
parameters, in this paper we define a Correlation Model through which we show
that the number of cases of COVID-19 and Dengue have very similar trends. We
then illustrate the relevance of our model by extending it to a Long short-term
memory model (LSTM) that incorporates both diseases, and using this to estimate
Dengue infections via COVID-19 data in countries that lack sufficient Dengue
data.


  The reduction of the computational effort is desirable for the simulation of
marine ecosystem models. Using a marine ecosystem model, the assessment and the
validation of annual periodic solutions (i.e., steady annual cycles) against
observational data are crucial to identify biogeochemical processes, which, for
example, influence the global carbon cycle. For marine ecosystem models, the
transport matrix method (TMM) already lowers the runtime of the simulation
significantly and enables the application of larger time steps
straightforwardly. However, the selection of an appropriate time step is a
challenging compromise between accuracy and shortening the runtime. Using an
automatic time step adjustment during the computation of a steady annual cycle
with the TMM, we present in this paper different algorithms applying either an
adaptive step size control or decreasing time steps in order to use the time
step always as large as possible without any manual selection. For these
methods and a variety of marine ecosystem models of different complexity, the
accuracy of the computed steady annual cycle achieved the same accuracy as
solutions obtained with a fixed time step. Depending on the complexity of the
marine ecosystem model, the application of the methods shortened the runtime
significantly. Due to the certain overhead of the adaptive method, the
computational effort may be higher in special cases using the adaptive step
size control. The presented methods represent computational efficient methods
for the simulation of marine ecosystem models using the TMM but without any
manual selection of the time step.


  AI-based protein structure prediction pipelines, such as AlphaFold2, have
achieved near-experimental accuracy. These advanced pipelines mainly rely on
Multiple Sequence Alignments (MSAs) as inputs to learn the co-evolution
information from the homologous sequences. Nonetheless, searching MSAs from
protein databases is time-consuming, usually taking dozens of minutes.
Consequently, we attempt to explore the limits of fast protein structure
prediction by using only primary sequences of proteins. HelixFold-Single is
proposed to combine a large-scale protein language model with the superior
geometric learning capability of AlphaFold2. Our proposed method,
HelixFold-Single, first pre-trains a large-scale protein language model (PLM)
with thousands of millions of primary sequences utilizing the self-supervised
learning paradigm, which will be used as an alternative to MSAs for learning
the co-evolution information. Then, by combining the pre-trained PLM and the
essential components of AlphaFold2, we obtain an end-to-end differentiable
model to predict the 3D coordinates of atoms from only the primary sequence.
HelixFold-Single is validated in datasets CASP14 and CAMEO, achieving
competitive accuracy with the MSA-based methods on the targets with large
homologous families. Furthermore, HelixFold-Single consumes much less time than
the mainstream pipelines for protein structure prediction, demonstrating its
potential in tasks requiring many predictions. The code of HelixFold-Single is
available at
https://github.com/PaddlePaddle/PaddleHelix/tree/dev/apps/protein_folding/helixfold-single,
and we also provide stable web services on
https://paddlehelix.baidu.com/app/drug/protein-single/forecast.


  People tend to have their social interactions with members of their own
community. Such group-structured interactions can have a profound impact on the
behaviors that evolve. Group structure affects the way people cooperate, and
how they reciprocate each other's cooperative actions. Past work has shown that
population structure and reciprocity can both promote the evolution of
cooperation. Yet the impact of these mechanisms has been typically studied in
isolation. In this work, we study how the two mechanisms interact. Using a
game-theoretic model, we explore how people engage in reciprocal cooperation in
group-structured populations, compared to well-mixed populations of equal size.
To derive analytical results, we focus on two scenarios. In the first scenario,
we assume a complete separation of time scales. Mutations are rare compared to
between-group comparisons, which themselves are rare compared to within-group
comparisons. In the second scenario, there is a partial separation of time
scales, where mutations and between-group comparisons occur at a comparable
rate. In both scenarios, we find that the effect of population structure
depends on the benefit of cooperation. When this benefit is small,
group-structured populations are more cooperative. But when the benefit is
large, well-mixed populations result in more cooperation. Overall, our results
reveal how group structure can sometimes enhance and sometimes suppress the
evolution of cooperation.


  Single-cell RNA-seq data allow the quantification of cell type differences
across a growing set of biological contexts. However, pinpointing a small
subset of genomic features explaining this variability can be ill-defined and
computationally intractable. Here we introduce MarkerMap, a generative model
for selecting minimal gene sets which are maximally informative of cell type
origin and enable whole transcriptome reconstruction. MarkerMap provides a
scalable framework for both supervised marker selection, aimed at identifying
specific cell type populations, and unsupervised marker selection, aimed at
gene expression imputation and reconstruction. We benchmark MarkerMap's
competitive performance against previously published approaches on real single
cell gene expression data sets. MarkerMap is available as a pip installable
package, as a community resource aimed at developing explainable machine
learning techniques for enhancing interpretability in single-cell studies.


  Previous work indicates that tropical forest can exist as an alternative
stable state to savanna. Therefore, perturbation by climate change or human
impact may lead to crossing of a tipping point beyond which there is rapid
forest dieback that is not easily reversed. A hypothesised mechanism for such
bistability is a feedback between fire and vegetation, where fire spreads as a
contagion process on grass patches. Theoretical models have largely implemented
this mechanism implicitly, by assuming a threshold dependence of fire spread on
flammable vegetation. Here, we show how the nonlinear dynamics and bistability
emerge spontaneously, without assuming equations or thresholds for fire spread.
We find that the forest geometry causes the nonlinearity that induces
bistability. We demonstrate this in three steps. First, we model forest and
fire as interacting contagion processes on grass patches, %using a cellular
automaton showing that spatial structure emerges due to two counteracting
effects on the forest perimeter: forest expansion by dispersal, and forest
erosion by fires originating in adjacent grassland. Then, we derive a
landscape-scale balance equation in which these two effects link forest
geometry and dynamics: forest expands proportionally to its perimeter, while it
shrinks proportionally to its perimeter weighted by adjacent grassland area.
Finally, we show that these perimeter quantities introduce nonlinearity in our
balance equation and lead to bistability. Relying on the link between structure
and dynamics, we propose a forest resilience indicator that could be used for
targeted conservation or restoration.


  Preeclampsia is a leading cause of maternal and fetal morbidity and
mortality. Currently, the only definitive treatment of preeclampsia is delivery
of the placenta, which is central to the pathogenesis of the disease.
Transcriptional profiling of human placenta from pregnancies complicated by
preeclampsia has been extensively performed to identify differentially
expressed genes (DEGs). The decisions to investigate DEGs experimentally are
biased by many factors, causing many DEGs to remain uninvestigated. A set of
DEGs which are associated with a disease experimentally, but which have no
known association to the disease in the literature are known as the ignorome.
Preeclampsia has an extensive body of scientific literature, a large pool of
DEG data, and only one definitive treatment. Tools facilitating knowledge-based
analyses, which are capable of combining disparate data from many sources in
order to suggest underlying mechanisms of action, may be a valuable resource to
support discovery and improve our understanding of this disease. In this work
we demonstrate how a biomedical knowledge graph (KG) can be used to identify
novel preeclampsia molecular mechanisms. Existing open source biomedical
resources and publicly available high-throughput transcriptional profiling data
were used to identify and annotate the function of currently uninvestigated
preeclampsia-associated DEGs. Experimentally investigated genes associated with
preeclampsia were identified from PubMed abstracts using text-mining
methodologies. The relative complement of the text-mined- and
meta-analysis-derived lists were identified as the uninvestigated
preeclampsia-associated DEGs (n=445), i.e., the preeclampsia ignorome. Using
the KG to investigate relevant DEGs revealed 53 novel clinically relevant and
biologically actionable mechanistic associations.


  Recent work on dissimilarity-based hierarchical clustering has led to the
introduction of global objective functions for this classical problem. Several
standard approaches, such as average linkage, as well as some new heuristics
have been shown to provide approximation guarantees. Here we introduce a broad
new class of objective functions which satisfy desirable properties studied in
prior work. Many common agglomerative and divisive clustering methods are shown
to be greedy algorithms for these objectives, which are inspired by related
concepts in phylogenetics.


  Biological genotypes do not code directly for phenotypes; developmental
physiology is the control layer that separates genomes from capacities
ascertained by selection. A key aspect is competency, as cells are not a
passive material but descendants of unicellular organisms with complex
context-sensitive capabilities. We used an evolutionary simulation in the
context of minimal artificial embryogeny to probe the effects of different
degrees of cellular competency on evolutionary dynamics. Virtual embryos
consisted of a single axis of positional information values provided by cells'
genomes, operated upon by an evolutionary cycle in which embryos' fitness was
proportional to monotonicity of the axial gradient. Evolutionary dynamics were
evaluated in two modes: hardwired "mosaic" development (genotype directly
encodes phenotype), and a more realistic mode in which cells interact prior to
evaluation by the fitness function ("regulative" development). Even minimal
competency with respect to improving their position in the embryo results in
better performance of the evolutionary search. Crucially, we observed that as
competency of cells masks the raw fitness of the genomes, the phenotypic
fitness gains are then mostly due to improvements of cells' developmental
problem-solving capacities, not the structural genome. This suggests the
existence of a powerful ratchet mechanism: evolution progressively becomes
locked in to improvements in the intelligence of its agential substrate, with
reduced pressure on the structural genome. A feedback loop in which evolution
increasingly puts more effort into the developmental software than perfecting
the hardware explains the very puzzling divergence of genome from anatomy in
species like planaria, identifies a possible drive for scaling intelligence
over time, and suggests strategies for engineering novel systems in silico and
in bioengineering.


  The Kimura equation is a degenerated partial differential equation of
drift-diffusion type used in population genetics. Its solution is required to
satisfy not only the equation but a series of conservation laws formulated as
integral constraints. In this work, we consider a population of two types
evolving without mutation or selection, the so-called neutral evolution. We
obtain explicit solutions in terms of Gegenbauer polynomials. To satisfy the
integral constraints it is necessary to prove new relations satisfied by the
Gegenbauer polynomials. The long-term in time asymptotic is also studied.


  The literature posits that an introduced predator population, is able to
drive it's target pest population extinct, if supplemented with high quality
additional food of quantity $\xi > \xi_{critical}$, \cite{SP11, SPV18, SPD17,
SPM13}. We show this approach leads to infinite time blow-up of the predator
population. We propose an alternate model in which the additional food induces
predator competition. Analysis herein indicates that there are threshold values
$c^{*}_{1} < c^{*}_{2} < c^{*}_{3}$ of the competition parameter $c$, s.t. when
$c < c^{*}_{1}$, the pest free state is globally stable, when $c^{*}_{2} < c <
c^{*}_{3}$, bi-stability is possible, and when $c^{*}_{3} < c$, up to three
interior equilibriums could exist. As $c$ and $\xi$-$c$ are varied, standard
co-dimension one and co-dimension two bifurcations are observed. The recent
dynamical systems literature involving predator competition, report several
non-standard bifurcations such as the saddle-node-transcritical bifurcation
(SNTC) occurring in co-dimension two \cite{KSV10, BS07}, and cusp-transcritical
bifurcation (CPTC) in co-dimension three, \cite{D20, BS07}. We show that in our
model structural symmetries can be exploited to construct a SNTC in
co-dimension two, and a CPTC also in co-dimension two. We further use these
symmetries to construct a novel pitchfork-transcritical bifurcation (PTC) in
co-dimension two, thus explicitly characterizing a new organizing center of the
model. Dynamics such as homoclinic orbits, concurrently occurring limit cycles,
and competition driven Turing patterns are also observed. Our findings indicate
that increasing additional food in predator-pest models, can hinder
bio-control, contrary to some of the literature. However, additional food that
also induces predator competition, leads to novel bio-control scenarios, and
complements the work in \cite{H21, B98, K04, D20, BS07, VH19}.


  Brucellosis is a zoonotic bacterial infectious disease that affects livestock
and wildlife. The Greater Yellowstone Ecosystem is the last area in the United
States where cattle are regularly infected with brucellosis. Even though
livestock are vaccinated, interactions with reservoir species still result in
spill-over cases to cattle. The National Academy of Sciences, Engineering, and
Medicine has indicated that modeling efforts should focus on the transmission
between elk and cattle and should be conducted to better understand the effects
of land-use change and landscape configuration on disease risk. This chapter
determines how the landscape's configuration, using the shape and amount of
habitat overlap between elk and cattle, contributes to cross-species
brucellosis transmission, and how land-use change translates to disease
prevalence. A mathematical-epidemiological model is combined with landscape
ecology metrics to estimate transmission rates between the species and model
disease spread. The results of this study can ultimately help stakeholders
develop policy for controlling brucellosis transmission between livestock and
elk in the Greater Yellowstone Area. In turn, this could lead to less disease
prevalence, reduce associated costs, and assist in population management.


  The Dynamical Survival Analysis (DSA) is a framework for modeling epidemics
based on mean field dynamics applied to individual (agent) level history of
infection and recovery. Recently, DSA has been shown to be an effective tool in
analyzing complex non-Markovian epidemic processes that are otherwise difficult
to handle using standard methods. One of the advantages of DSA is its
representation of typical epidemic data in a simple although not explicit form
that involves solutions of certain differential equations. In this work we
describe how a complex non-Markovian DSA model may be applied to a specific
data set with the help of appropriate numerical and statistical schemes. The
ideas are illustrated with a data example of the COVID-19 epidemic in Ohio.


  Diseased conditions are a consequence of some abnormality that are associated
with clinical conditions in numerous cells and tissues affecting various
organs. The common role of EBV (Epstein-Barr virus) in causing infectious
mononucleosis (IM) affecting B-cells and epithelial cells and the development
of EBV-associated cancers has been an area of active research. Investigating
such significant interactions may help discover new therapeutic targets for
certain EBV-associated lymphoproliferative (Burkitt's Lymphoma and Hodgkin's
Lymphoma) and non-lymphoproliferative diseases (Gastric cancer and
Nasopharyngeal cancer). Based on the DisGeNET (v7.0) data set, we constructed a
disease-gene network bipartite graph to identify genes that are involved in
various carcinomas namely, gastric cancer (GC), nasopharyngeal cancer (NPC),
Hodgkin's lymphoma (HL) and Burkitt's lymphoma (BL). Using the community
detection algorithm (Louvain method), we identified communities followed by
functional enrichment using over-representation analysis methodology. In this
study, we identified the modular communities to explore the relation of this
common causative pathogen (EBV) with different carcinomas such as GC, NPC, HL
and BL. We could identify the top 10 genes as CASP10, BRAF, NFKBIA, IFNA2,
GSTP1, CSF3, GATA3, UBR5, AXIN2 and POLE based on their degree of distribution.
Further over-representation analysis showed that the ABL1 gene was
significantly over-represented in 3 out of 9 critical biological processes. As
a result, we can infer that the EBV pathogen is selective in targeting critical
pathways to bring about cellular growth arrest/apoptosis and interfering with
vital biological processes, including the TP53 network of genes that leads to
further proliferation of damage to vital cellular activities.


  Plasticity-led evolution is a form of evolution where a change in the
environment induces novel traits via phenotypic plasticity, after which the
novel traits are genetically accommodated over generations under the novel
environment. This mode of evolution is expected to resolve the problem of
gradualism (i.e., evolution by the slow accumulation of mutations that induce
phenotypic variation) implied by the Modern Evolutionary Synthesis, in the face
of a large environmental change. While experimental works are essential for
validating that plasticity-led evolution indeed happened, we need computational
models to gain insight into its underlying mechanisms and make qualitative
predictions. Such computational models should include the developmental process
and gene-environment interactions in addition to genetics and natural
selection. We point out that gene regulatory network models can incorporate all
the above notions. In this review, we highlight results from computational
modelling of gene regulatory networks that consolidate the criteria of
plasticity-led evolution. Since gene regulatory networks are mathematically
equivalent to artificial recurrent neural networks, we also discuss their
analogies and discrepancies, which may help further understand the mechanisms
underlying plasticity-led evolution.


  We interpret the Moran model of natural selection and drift as an algorithm
for learning features of a simplified fitness landscape, specifically genotype
superiority. This algorithm's efficiency in extracting these characteristics is
evaluated by comparing it to a novel Bayesian learning algorithm developed
using information-theoretic tools. This algorithm makes use of a communication
channel analogy between an environment and an evolving population. We use the
associated channel-rate to determine an informative population-sampling
procedure. We find that the algorithm can identify genotype superiority faster
than the Moran model but at the cost of larger fluctuations in uncertainty.


  We present a novel method to estimate the dominant eigenvalue and eigenvector
pair of any non-negative real matrix via graph infection. The key idea in our
technique lies in approximating the solution to the first-order matrix ordinary
differential equation (ODE) with the Euler method. Graphs, which can be
weighted, directed, and with loops, are first converted to its adjacency matrix
A. Then by a naive infection model for graphs, we establish the corresponding
first-order matrix ODE, through which A's dominant eigenvalue is revealed by
the fastest growing term. When there are multiple dominant eigenvalues of the
same magnitude, the classical power iteration method can fail. In contrast, our
method can converge to the dominant eigenvalue even when same-magnitude
counterparts exist, be it complex or opposite in sign. We conduct several
experiments comparing the convergence between our method and power iteration.
Our results show clear advantages over power iteration for tree graphs,
bipartite graphs, directed graphs with periods, and Markov chains with
spider-traps. To our knowledge, this is the first work that estimates dominant
eigenvalue and eigenvector pair from the perspective of a dynamical system and
matrix ODE. We believe our method can be adopted as an alternative to power
iteration, especially for graphs.


  This paper models brucellosis transmission between elk, cattle, and bison, of
high conservation value, in the Greater Yellowstone Ecosystem (GYE). It aims to
show how landscape changes in the GYE concomitantly impact brucellosis
prevalence in the three species. The approach allows us to see how landscape
changes in one location influence disease prevalence in populations elsewhere.
The model uses the fact that the populations are configured in such a way that
elk are an intermediary for disease transmission between cattle and bison.
Using landscape ecology metrics applied to the habitat overlaps between elk and
cattle and between elk and bison, the landscape parameters are varied to
determine how disease propagates throughout the ecosystem as land-use change
occurs. Result aim to provide insights into how land management used for the
control of disease spread between cattle and elk, may impact (and be impacted
by) disease prevalence in bison.


  A variety of approaches using compartmental models have been used to study
the COVID-19 pandemic and the usage of machine learning methods with these
models has had particularly notable success. We present here an approach toward
analyzing accessible data on Covid-19's U.S. development using a variation of
the "Physics Informed Neural Networks" (PINN) which is capable of using the
knowledge of the model to aid learning. We illustrate the challenges of using
the standard PINN approach, then how with appropriate and novel modifications
to the loss function the network can perform well even in our case of
incomplete information. Aspects of identifiability of the model parameters are
also assessed, as well as methods of denoising available data using a wavelet
transform. Finally, we discuss the capability of the neural network methodology
to work with models of varying parameter values, as well as a concrete
application in estimating how effectively cases are being tested for in a
population, providing a ranking of U.S. states by means of their respective
testing.


  Electroporation is a very useful tool for drug delivery into various diseased
tissues of the human body. This technique helps to improve the clinical
treatment by transferring drugs into the targeted cells rapidly. In
electroporation, drug particles enter easily into the intracellular compartment
through the temporarily permeabilized cell membrane due to the applied electric
field. In this work, a mathematical model of drug delivery focusing on
reversible tissue electroporation is presented. In addition, the thermal
effects on the tissue, which is an outcome of Joule heating, are also
considered. This model introduces a time-dependent mass transfer coefficient,
which is significant to drug transport. Multiple pulses with low voltage are
applied to reach sufficient drugs into the targeted cells. According to the
physical circumstances, a set of differential equations are considered and
solved. The changes in drug concentration with different parameters (e.g.,
diffusion coefficient, drug permeability, pulse length, and pulse number) are
analyzed. The model optimizes the electroporation parameters to uptake drugs
sufficiently into the cells with no thermal damage. This model can be used in
clinical experiments to predict drug uptake into the infected cells by
controlling the model parameters according to the nature of infections.


  Drug delivery through electroporation could be highly beneficial for the
treatment of different types of diseased tissues within the human body. In this
work, a mathematical model of reversible tissue electroporation is presented
for injecting drug into the diseased cells. The model emphasizes the tissue
boundary where the drug is injected as a point source. Drug loss from the
tissue boundaries through extracellular space is studied. Multiple pulses are
applied to deliver a sufficient amount of drug into the targeted cells. The set
of differential equations that model the physical circumstances are solved
numerically. This model obtains a mass transfer coefficient in terms of pore
fraction coefficient and drug permeability. It controls the drug transport from
extracellular to intracellular space. The drug penetration throughout the
tissue is captured for the application of different pulses. The boundary
effects on drug concentration are highlighted in this study. The advocated
model is able to perform homogeneous drug transport into the cells so that the
affected tissue is treated completely. This model can be applied to optimize
clinical experiments by avoiding the lengthy and costly in vivo and in vitro
experiments.


  We present a scenario for the origin of biological coding. In this context,
coding is a semiotic relationship between chemical information stored in one
location that links to chemical information stored in a separate location.
Coding originated by the cooperative interaction of two, originally separate
collectively autocatalytic sets, one for nucleic acids and one for peptides.
When these two sets interacted, a series of RNA-folding-directed processes led
to their joint cooperativity. The amino acyl adenylate, today amino acid-AMP,
was the first covalent association made by these two collectively autocatalytic
sets and solidified their interdependence. This molecule is a palimpsest of
this era, and is a relic of the original semiotic, and thus coding,
relationship between RNA and proteins. More defined coding was driven by
selection pressure to eliminate waste in the collective autocatalytic sets.
Eventually a 1:1 relationship between single amino acids and short RNA pieces
(e.g., three nucleotides) was established, leading to what is today known as
the genetic code. Transfer RNA aminoacylating enzymes, or aaRSs, arose
concomitantly with the advent of specific coding. The two classes of aaRS
enzymes are remnants of the duality of complementary information in two nucleic
acid strands, as originally postulated by Rodin and Ohno. Every stage in the
evolution of coding was driven by the downward selection on the components of a
system to satisfy the Kantian whole. Coding was ultimately forced because there
were at least two chemically distinct classes of polymers needed for open-ended
evolution; systems with only one polymer cannot exhibit this characteristic.
Coding is thus synonymous with life as we know it, and can be thought of as a
phase transition in the history of the universe.


  In the analysis of complex ecosystems it is common to use random interaction
coefficients, often assumed to be such that all species are statistically
equivalent. In this work we relax this assumption by choosing interactions
according to the cascade model, which we incorporate into a generalised
Lotka-Volterra dynamical system. These interactions impose a hierarchy in the
community. Species benefit more, on average, from interactions with species
further below them in the hierarchy than from interactions with those above.
Using dynamic mean-field theory, we demonstrate that a strong hierarchical
structure is stabilising, but that it reduces the number of species in the
surviving community, as well as their abundances. Additionally, we show that
increased heterogeneity in the variances of the interaction coefficients across
positions in the hierarchy is destabilising. We also comment on the structure
of the surviving community and demonstrate that the abundance and probability
of survival of a species is dependent on its position in the hierarchy.


  Rapid discovery of new reactions and molecules in recent years has been
facilitated by the advancements in high throughput screening, accessibility to
a much more complex chemical design space, and the development of accurate
molecular modeling frameworks. A holistic study of the growing chemistry
literature is, therefore, required that focuses on understanding the recent
trends and extrapolating them into possible future trajectories. To this end,
several network theory-based studies have been reported that use a directed
graph representation of chemical reactions. Here, we perform a study based on
representing chemical reactions as hypergraphs where the hyperedges represent
chemical reactions and nodes represent the participating molecules. We use a
standard reactions dataset to construct a hypernetwork and report its
statistics such as degree distributions, average path length, assortativity or
degree correlations, PageRank centrality, and graph-based clusters (or
communities). We also compute each statistic for an equivalent directed graph
representation of reactions to draw parallels and highlight differences between
the two. To demonstrate the AI applicability of hypergraph reaction
representation, we generate dense hypergraph embeddings and use them in the
reaction classification problem. We conclude that the hypernetwork
representation is flexible, preserves reaction context, and uncovers hidden
insights that are otherwise not apparent in a traditional directed graph
representation of chemical reactions.


  A nonlinear partial differential equation (PDE) based compartmental model of
COVID-19 provides a continuous trace of infection over space and time. Finer
resolutions in the spatial discretization, the inclusion of additional model
compartments and model stratifications based on clinically relevant categories
contribute to an increase in the number of unknowns to the order of millions.
We adopt a parallel scalable solver allowing faster solutions for these high
fidelity models. The solver combines domain decomposition and algebraic
multigrid preconditioners at multiple levels to achieve the desired strong and
weak scalability. As a numerical illustration of this general methodology, a
five-compartment susceptible-exposed-infected-recovered-deceased (SEIRD) model
of COVID-19 is used to demonstrate the scalability and effectiveness of the
proposed solver for a large geographical domain (Southern Ontario). It is
possible to predict the infections up to three months for a system size of 92
million (using 1780 processes) within 7 hours saving months of computational
effort needed for the conventional solvers.


  Stochastic models of sequential mutation acquisition are widely used to
quantify cancer and bacterial evolution. Across manifold scenarios, recurrent
research questions are: how many cells are there with $n$ alterations, and how
long will it take for these cells to appear. For exponentially growing
populations, these questions have been tackled only in special cases so far.
Here, within a multitype branching process framework, we consider a general
mutational path where mutations may be advantageous, neutral or deleterious. In
the biologically relevant limiting regimes of large times and small mutation
rates, we derive probability distributions for the number, and arrival time, of
cells with n mutations. Surprisingly, the two quantities respectively follow
Mittag-Leffler and logistic distributions regardless of $n$ or the mutations'
selective effects. Our results provide a rapid method to assess how altering
the fundamental division, death, and mutation rates impacts the arrival time,
and number, of mutant cells. We highlight consequences for mutation rate
inference in fluctuation assays.


  Vaccination is essential for the management of infectious diseases, many of
which continue to pose devastating public health and economic challenges across
the world. However, many vaccines are imperfect having only a partial
protective effect in decreasing disease transmission and/or favouring recovery
of infected individuals, and possibly exhibiting trade-off between these two
properties. Furthermore, population turnover, that is the rate at which
individuals enter and exit the population, is another key factor determining
the epidemiological dynamics. While these factors have yet been studied
separately, we investigate the interplay between the efficiency and property of
an imperfect vaccine and population turnover. We build a mathematical model
with frequency incidence rate, a recovered compartment, and an heterogeneous
host population with respect to vaccination. We first compute the basic
reproduction number $\mathcal{R}_0$ and study the global stability of the
equilibrium points. Using a sensitivity analysis, we then assess the most
influential parameters determining the total number of infected and
$\mathcal{R}_0$ over time. We derive analytically and numerically conditions
for the vaccination coverage and efficiency to achieve disease eradication
($\mathcal{R}_0 < 1$) assuming different intensity of the population turnover
(weak and strong), vaccine properties (transmission and/or recovery) and
trade-off between the latter. We show that the minimum vaccination coverage
increases with lower population turnover, decreases with higher vaccine
efficiency (transmission or recovery), and is increased/decreased by up to 15\%
depending on the trade-off between the vaccine properties. We conclude that the
coverage target for vaccination campaigns should be evaluated based on the
interplay between these factors.


  Predicting drug-target interaction is key for drug discovery. Recent deep
learning-based methods show promising performance but two challenges remain:
(i) how to explicitly model and learn local interactions between drugs and
targets for better prediction and interpretation; (ii) how to generalize
prediction performance on novel drug-target pairs from different distribution.
In this work, we propose DrugBAN, a deep bilinear attention network (BAN)
framework with domain adaptation to explicitly learn pair-wise local
interactions between drugs and targets, and adapt on out-of-distribution data.
DrugBAN works on drug molecular graphs and target protein sequences to perform
prediction, with conditional domain adversarial learning to align learned
interaction representations across different distributions for better
generalization on novel drug-target pairs. Experiments on three benchmark
datasets under both in-domain and cross-domain settings show that DrugBAN
achieves the best overall performance against five state-of-the-art baselines.
Moreover, visualizing the learned bilinear attention map provides interpretable
insights from prediction results.


  COVID-19 pandemic has had a disruptive and irreversible impact globally, yet
traditional epidemiological modeling approaches such as the
susceptible-infected-recovered (SIR) model have exhibited limited effectiveness
in forecasting of the up-to-date pandemic situation. In this work,
susceptible-vaccinated-exposed-infected-dead-recovered (SVEIDR) model and its
variants -- aged and vaccination-structured SVEIDR models -- are introduced to
encode the effect of social contact for different age groups and vaccination
status. Then, we implement the physics-informed neural network (PiNN) on both
simulated and real-world data. The PiNN model enables robust analysis of the
dynamic spread, prediction, and parameter optimization of the COVID-19
compartmental models. The models exhibit relative root mean square error
(RRMSE) of $<4\%$ for all components and provide incubation, death, and
recovery rates of $\gamma= 0.0130$, $\lambda=0.0001$, and $\rho=0.0037$,
respectively, for the first 310 days of the epidemic in the US with RRMSE of
$<0.35\%$ for all components. To further improve the model performance,
temporally varying parameters can be included, such as vaccination,
transmission, and incubation rates. Our implementation highlights PiNN as a
reliable candidate approach for forecasting real-world data and can be applied
to other compartmental model variants of interest.


  The Lotka-Volterra predator-prey model still represents the paradigm for the
description of the competition in population dynamics. Despite its extreme
simplicity, it does not admit an analytical solution, and for this reason,
numerical integration methods are usually adopted to apply it to various fields
of science. The aim of the present work is to investigate the existence of new
predator-prey models sharing the broad features of the standard Lotka-Volterra
model and, at the same time, offer the advantage of possessing exact analytical
solutions. To this purpose, a general Hamiltonian formalism, which is suitable
for treating a large class of predator-prey models in population dynamics
within the same framework, has been developed as a first step. The only
existing model having the property of admitting a simple exact analytical
solution, is identified within the above class of models. The solution of this
special predator-prey model is obtained explicitly, in terms of known
elementary functions, and its main properties are studied. Finally, the
generalization of this model, based on the concept of power-law competition, as
well as its extension to the case of $N$-component competition systems, are
considered.


  The advent of Zero-Determinant (ZD) strategies has reshaped the study of
reciprocity and cooperation in the iterated Prisoner's Dilemma games. The
ramification of ZD strategies has been demonstrated through their ability to
unilaterally enforce a linear relationship between their own average payoff and
that of their co-player. Common practice conveniently represents this
relationship by a straight line in the parametric plot of pairwise payoffs. Yet
little attention has been paid to studying the actual geometry of the strategy
space of all admissible ZD strategies. Here, our work offers intuitive
geometric relationships between different classes of ZD strategies as well as
nontrivial geometric interpretations of their specific parameterizations.
Adaptive dynamics of ZD strategies further reveals the unforeseen connection
between general ZD strategies and the so-called equalizers that can set any
co-player's payoff to a fixed value. We show that the class of equalizers
forming a hyperplane is the critical equilibrium manifold, only part of which
is stable. The same hyperplane is also a separatrix of the
cooperation-enhancing region where the optimum response is to increase
cooperation for each of the four payoff outcomes. Our results shed light on the
simple but elegant geometry of ZD strategies that is previously overlooked.


  A new deep neural network based on the WaveNet architecture (WNN) is
presented, which is designed to grasp specific patterns in the NMR spectra.
When trained at a fixed non-uniform sampling (NUS) schedule, the WNN benefits
from pattern recognition of the corresponding point spread function (PSF)
pattern produced by each spectral peak resulting in the highest quality and
robust reconstruction of the NUS spectra as demonstrated in simulations and
exemplified in this work on 2D 1H-15N correlation spectra of three
representative globular proteins with different sizes: Ubiquitin (8.6 kDa),
Azurin (14 kDa), and Malt1 (44 kDa). The pattern recognition by WNN is also
demonstrated for successful virtual homo-decoupling in a 2D methyl 1H-13 HMQC
spectrum of MALT1. We demonstrate using WNN that prior knowledge about the NUS
schedule, which so far was not fully exploited, can be used for designing new
powerful NMR processing techniques that surpass the existing algorithmic
methods.


  This paper studied the relationship between the decomposition rate of fungi
and temperature, humidity, fungus elongation, moisture tolerance and fungus
density in a given volume in the presence of a variety of fungi, and
established a series of models to describe the decomposition of fungi in
different states. Since the volume of soil was given in this case, the latter
two characteristics could be attributed to the influence of the number of
fungal population on the decomposition rate. Based on the Logistic model, the
relationship between the number of population and time was established, and
finally the number of fungi in the steady state was obtained The interaction
between different species of fungi was analyzed by Lotka-Volterra model, and
the decomposition rate of various fungal combinations in different environments
was obtained. After studying the one and two cases, we can extrapher from one
to the other, and the community consisting of n fungal populations will be
similar to the community consisting of n+1 fungal populations. After the study,
we substituted the collected data into the model and found that the fungal
community composed of two kinds of fungi had a lower decomposition rate of
ground decomposition or wooden fiber than that of a single kind of fungus for
the same kind of substance. We found that the fungus in warm and humid
environment of decomposition rate is highest, the change of the atmospheric
cause some fungal population growth rate decreases, there are also some will
increase, which is associated with the nature of fungi.We analyzed the
influence of environmental factors, namely temperature and humidity, on the
model.


  Over the past several decades, algebraic geometry has provided innovative
approaches to biological experimental design that resolved theoretical
questions and improved computational efficiency. However, guaranteeing
uniqueness and perfect recovery of models are still open problems. In this work
we study the problem of uniqueness of wiring diagrams. We use as a modeling
framework polynomial dynamical systems and utilize the correspondence between
simplicial complexes and square-free monomial ideals from Stanley-Reisner
theory to develop theory and construct an algorithm for identifying input data
sets $V\subset \mathbb F_p^n$ that are guaranteed to correspond to a unique
minimal wiring diagram regardless of the experimental output. We apply the
results on a tumor-suppression network mediated by epidermal derived growth
factor receptor and demonstrate how careful experimental design decisions can
lead to a unique minimal wiring diagram identification. One of the insights of
the theoretical work is the connection between the uniqueness of a wiring
diagram for a given $V\subset \mathbb F_p^n$ and the uniqueness of the reduced
Gr\"obner basis of the polynomial ideal $I(V)\subset \mathbb F_p[x_1,\ldots,
x_n]$. We discuss existing results and introduce a new necessary condition on
the points in $V$ for uniqueness of the reduced Gr\"obner basis of $I(V)$.
These results also point to the importance of the relative proximity of the
experimental input points on the number of minimal wiring diagrams, which we
then study computationally. We find that there is a concrete heuristic way to
generate data that tends to result in fewer minimal wiring diagrams.


  The Levinthal paradox exposes many critical questions on the protein folding
problem, among which we could point out why proteins can reach their native
state in a biologically reasonable time. A proper answer to this question is of
foremost importance for evolutive biology since it enables us to understand
life as we know it. Preliminary results, based on the upper bound protein
marginal-stability limit, together with transition state theory arguments, lead
us to show that two-state proteins must reach their native state in, at most,
seconds rather than ($\sim 10^{27}$) years -- as indicated by a naive solution
of the Levinthal paradox. This outcome -- added to the amide hydrogen-exchange
protection factors analysis -- makes it possible for us to suggest how a
protein point mutations and/or post-translational modifications impact its
folding time scales but not its upper bound limit that obeys the physics ruling
the process. Noteworthy for almost 50 years, the protein folding problem --mas
the Levinthal paradox -- has been a topic of passionate debate because
Anfinsen's challenge -- how a sequence encodes its folding -- remains unsolved
despite the smashing success of accurately predicting the protein
tridimensional structures by state-of-the-art numerical-methods. Aimed to
unlock this long-standing challenge, we propose a new perspective of protein
folding, specifically, as a problem that should be devised as an 'analytic
whole' -- a Leibniz & Kant's notion. This viewpoint might help us decode
Anfinsen's challenge and, thus, open new avenues for future research in the
protein folding field.


  Isoforms are mRNAs produced from the same gene site in the phenomenon called
Alternative Splicing. Studies have shown that more than 95% of human multi-exon
genes have undergone alternative splicing. Although there are few changes in
mRNA sequence, They may have a systematic effect on cell function and
regulation. It is widely reported that isoforms of a gene have distinct or even
contrasting functions. Most studies have shown that alternative splicing plays
a significant role in human health and disease. Despite the wide range of gene
function studies, there is little information about isoforms' functionalities.
Recently, some computational methods based on Multiple Instance Learning have
been proposed to predict isoform function using gene function and gene
expression profile. However, their performance is not desirable due to the lack
of labeled training data. In addition, probabilistic models such as Conditional
Random Field (CRF) have been used to model the relation between isoforms. This
project uses all the data and valuable information such as isoform sequences,
expression profiles, and gene ontology graphs and proposes a comprehensive
model based on Deep Neural Networks. The UniProt Gene Ontology (GO) database is
used as a standard reference for gene functions. The NCBI RefSeq database is
used for extracting gene and isoform sequences, and the NCBI SRA database is
used for expression profile data. Metrics such as Receiver Operating
Characteristic Area Under the Curve (ROC AUC) and Precision-Recall Under the
Curve (PR AUC) are used to measure the prediction accuracy.


  Consider a migration process based on a closed network of N stations with K_N
customers. Each station is a ./M/\infty queue with service (migration) rate mu.
Upon departure, a customer is routed at random to another station. In addition
to migration, these customers are subject to an SIS (Susceptible, Infected,
Susceptible) dynamics: customers are either I for infected, or S for
susceptible. They can swap their state either from I to S or from S to I only
in stations. At any station, each S customer becomes I with rate alpha Y if
there are Y infected customers in the station, and each I customer recovers and
becomes S with rate beta. We let N tend to infinity and assume that lim_{N\to
infty} K_N/N= eta>0. The main problem is about the set of parameters for which
there exists a stationary regime where the epidemic survives in the
thermodynamic limit. We establish several structural properties of the system,
which allow us to give the phase transition diagram of this thermodynamic limit
w.r.t. eta. The analysis of the SIS model reduces to that of a wave-type PDE
for which we found no explicit solution. This SIS model is one among several
companion stochastic processes with migration and contagion. Two of them are
discussed as they provide some bounds and approximations to SIS. These two
variants are the DOCS (Departure On Change of State) and the AIR (Averaged
Infection Rate), which both admit closed-form solutions. The AIR system is a
mean-field model where the infection mechanism is based on the empirical
average of the number of infected customers in all stations. The latter admits
a product-form solution. DOCS features accelerated migration in that each
change of SIS state implies an immediate departure. It leads to another
wave-type PDE that admits a closed-form solution.


  Avian flocks display a wide variety of flight behaviors, including steady
directed translation of center of mass, rapid change of overall morphology,
re-shuffling of positions of individuals within a persistent form, etc. These
behaviors may be viewed as flock-scale strategies, emerging from interactions
between individuals, accomplishing some collective adaptive purpose such as
finding a roost, or mitigating the danger from predator attacks. While we do
not conceive the flock as a single cognitive agent, the moment-to-moment
decisions of individuals, influenced by their neighbors, appear as if to
realize collective strategies that are cognizant of purpose. In this paper, we
identify the actions of the flock as allocation of energetic resources, and
thereby associate a cognitive cost to behavior. Our notion of cognitive cost
reflects the burden arising from rapid re-allocation of resource. Using a
recently developed natural geometric approach to kinetic energy allocation, we
map the flock behavior to a temporal signature on the standard (probability)
simplex. Given the signature of a flocking event, we calculate the cognitive
cost as a solution to an optimal control problem based on a game-theoretic
model. Alternatively, one can associate to a signature an entropic cost. These
two cost measures, when applied to data on starling flocks, show a consistent
spread in value across events, and we suggest the possibility that higher cost
may arise from predator attacks.


  We develop numerical and analytical approaches to calculate mutual
information between complete paths of two molecular components embedded into a
larger reaction network. In particular, we focus on a continuous-time Markov
chain formalism, frequently used to describe intracellular processes involving
lowly abundant molecular species. Previously, we have shown how the path mutual
information can be calculated for such systems when two molecular components
interact directly with one another with no intermediate molecular components
being present. In this work, we generalize this approach to biochemical
networks involving an arbitrary number of molecular components. We present an
efficient Monte Carlo method as well as an analytical approximation to
calculate the path mutual information and show how it can be decomposed into a
pair of transfer entropies that capture the causal flow of information between
two network components. We apply our methodology to study information transfer
in a simple three-node feedforward network, as well as a more complex positive
feedback system that switches stochastically between two metastable modes.


  For the reaction networks with zero-one stoichiometric coefficients (or
simply zero-one networks), we prove that if a network admits a Hopf
bifurcation, then the rank of the stoichiometric matrix is at least four. As a
corollary, we show that if a zero-one network admits a Hopf bifurcation, then
it contains at least four species and five reactions. As applications, we show
that there exist rank-four subnetworks, which have the capacity for Hopf
bifurcations/oscillations, in two biologically significant networks: the MAPK
cascades and the ERK network. We provide a computational tool for computing all
four-species, five-reaction, zero-one networks that have the capacity for Hopf
bifurcations.


  RNAs self-interact through hydrogen-bond base-pairing between nucleotides and
fold into specific, stable structures that substantially govern their
biochemical behaviour. Experimental characterization of these structures
remains difficult, hence the desire to predict them computationally from
sequence information. However, correctly predicting even the base pairs
involved in the folded structure of an RNA, known as secondary structure, from
its sequence using minimum free energy models is understood to be NP-hard.
Classical approaches rely on heuristics or avoid considering pseudoknots in
order to render this problem more tractable, with the cost of inexactness or
excluding an entire class of important RNA structures. Given their prospective
and demonstrable advantages in certain domains, including combinatorial
optimization, quantum computing approaches by contrast have the potential to
compute the full RNA folding problem while remaining more feasible and exact.
Herein, we present a physically-motivated QUBO model of the RNA folding problem
amenable to both quantum annealers and circuit-model quantum computers and
compare the performance of this formulation versus current RNA folding QUBOs
after tuning the parameters of all against known RNA structures using an
approach we call "variational hybrid quantum annealing".


  Branching process inspired models are widely used to estimate the effective
reproduction number -- a useful summary statistic describing an infectious
disease outbreak -- using counts of new cases. Case data is a real-time
indicator of changes in the reproduction number, but is challenging to work
with because cases fluctuate due to factors unrelated to the number of new
infections. We develop a new model that incorporates the number of diagnostic
tests as a surveillance model covariate. Using simulated data and data from the
SARS-CoV-2 pandemic in California, we demonstrate that incorporating tests
leads to improved performance over the state-of-the-art.


  Land conversion and the resulting contact between domesticated and wild
species has arguably been the single largest contributor to the emergence of
novel epizootic and zoonotic diseases in the past century. An unintended
consequence of these interactions is zoonotic or epizootic disease spillovers
from wild species to humans and their domesticates. Disease spillovers are edge
effects of land conversion and are sensitive to the size and shape of converted
areas. We combine spatial metrics from landscape ecology with theoretical
epidemiological models to understand how the size and shape of land conversion
affect epizootic and zoonotic disease transmission of single and two species
populations. We show that the less compact the converted area, and the the
greater the depth of the contact zone, the more rapidly will an introduced
disease spread through the domesticated population.


  Electroporation method is a useful tool for delivering drugs into various
diseased tissues in the human body. As a result of an applied electric field,
drug particles enter the intracellular compartment through the temporarily
permeabilized cell membrane. Consequently, electroporation method allows better
penetration of the drug into the diseased tissue and improves treatment
clinically. In this study, a more generalized model of drug transport in a
single-cell is proposed. The model is able to capture non-homogeneous drug
transport in the cell due to non-uniform cell membrane permeabilization.
Several numerical experiments are conducted to understand the effects of
electric field and drug permeability on drug uptake into the cell. Through
investigation, the appropriate electric field and drug permeability are
identified that lead to sufficient drug uptake into the cell. This model can be
used by experimentalists to get information prior to conduct any experiment,
and it may help reduce the number of actual experiments that might be conducted
otherwise.


  The GC content of prokaryotic genomes is species-specific and takes values
from 16 to 77 percent. There are currently no accepted explanations for this
diversity of selection for GC content. We analyzed the known correlations
between GC content, genome size and amino acid cost in thousands of
prokaryotes, together with new or recently compiled data on cell shape and
volume, duplication times, motility, nutrient assimilation, sporulation,
defense mechanisms and Gram staining. GC content integrates well with these
traits into r/K selection theory when phenotypic plasticity is considered. High
GC content prokaryotes are r-strategists with cheaper descendants thanks to a
lower average amino acid metabolic cost and a smaller cell volume, colonize
unstable environments thanks to flagella and a bacillus form and are
generalists in terms of resource opportunism and the ability to defend
themselves from various hazards. Low GC content prokaryotes are K-strategists
specialized for stable environments that maintain homeostasis via a high-cost
outer cell membrane and endospore formation as a response to nutrient
deprivation and attain a higher nutrient-to-biomass yield. The lower proteome
cost of high GC content prokaryotes is driven by the association between
GC-rich codons and cheaper amino acids in the genetic code, while the
correlation between GC content and genome size may be partly due to a shift in
the functional repertoire of genomes driven by r/K selection. In all, we show
that molecular diversity in the GC content of prokaryotes and the corresponding
species diversity may be a consequence of ecological r/K selection