  Most organisms grow in space, whether they are viruses spreading within a
host tissue or invasive species colonizing a new continent. Evolution typically
selects for higher expansion rates during spatial growth, but it has been
suggested that slower expanders can take over under certain conditions. Here,
we report an experimental observation of such population dynamics. We
demonstrate that the slower mutants win not only when the two types are
intermixed at the front but also when they are spatially segregated into
sectors. The latter was thought to be impossible because previous studies
focused exclusively on the global competitions mediated by expansion velocities
but overlooked the local competitions at sector boundaries. We developed a
theory of sector geometry that accounts for both local and global competitions
and describes all possible sector shapes. In particular, the theory predicted
that a slower, but more competitive, mutant forms a dented V-shaped sector as
it takes over the expansion front. Such sectors were indeed observed
experimentally and their shapes matched up quantitatively with the theory. In
simulations, we further explored several mechanism that could provide slow
expanders with a local competitive advantage and showed that they are all
well-described by our theory. Taken together, our results shed light on
previously unexplored outcomes of spatial competition and establish a universal
framework to understand evolutionary and ecological dynamics in expanding
populations.


  While atom tracking with isotope-labeled compounds is an essential and
sophisticated wet-lab tool in order to, e.g., illuminate reaction mechanisms,
there exists only a limited amount of formal methods to approach the problem.
Specifically when large (bio-)chemical networks are considered where reactions
are stereo-specific, rigorous techniques are inevitable. We present an approach
using the right Cayley graph of a monoid in order to track atoms concurrently
through sequences of reactions and predict their potential location in product
molecules. This can not only be used to systematically build hypothesis or
reject reaction mechanisms (we will use the ANRORC mechanism "Addition of the
Nucleophile, Ring Opening, and Ring Closure" as an example), but also to infer
naturally occurring subsystems of (bio-)chemical systems. Our results include
the analysis of the carbon traces within the TCA cycle and infer subsystems
based on projections of the right Cayley graph onto a set of relevant atoms.


  Probabilistic biological network growth models have been utilized for many
tasks including but not limited to capturing mechanism and dynamics of
biological growth activities, null model representation, capturing anomalies,
etc. Well-known examples of these probabilistic models are Kronecker model,
preferential attachment model, and duplication-based model. However, we should
frequently keep developing new models to better fit and explain the observed
network features while new networks are being observed. Additionally, it is
difficult to develop a growth model each time we study a new network. In this
paper, we propose BioCode, a framework to automatically discover novel
biological growth models matching user-specified graph attributes in directed
and undirected biological graphs. BioCode designs a basic set of instructions
which are common enough to model a number of well-known biological graph growth
models. We combine such instruction-wise representation with a genetic
algorithm based optimization procedure to encode models for various biological
networks. We mainly evaluate the performance of BioCode in discovering models
for biological collaboration networks, gene regulatory networks, metabolic
networks, and protein interaction networks which features such as
assortativity, clustering coefficient, degree distribution closely match with
the true ones in the corresponding real biological networks. As shown by the
tests on the simulated graphs, the variance of the distributions of biological
networks generated by BioCode is similar to the known models' variance for
these biological network types.


  Cooperation and competition between pathogens can alter the amount of
individuals affected by a co-infection. Nonetheless, the evolution of the
pathogens' behavior has been overlooked. Here, we consider a co-evolutionary
model where the simultaneous spreading is described by a two-pathogen
susceptible-infected-recovered model in an either synergistic or competitive
manner. At the end of each epidemic season, the pathogens species reproduce
according to their fitness that, in turn, depends on the payoff accumulated
during the spreading season in a hawk-and-dove game. This co-evolutionary model
displays a rich set of features. Specifically, the evolution of the pathogens'
strategy induces abrupt transitions in the epidemic prevalence. Furthermore, we
observe that the long-term dynamics results in a single, surviving pathogen
species, and that the cooperative behavior of pathogens can emerge even under
unfavorable conditions.


  Evolutionary models used for describing molecular sequence variation suppose
that at a non-recombining genomic segment, sequences share ancestry that can be
represented as a genealogy--a rooted, binary, timed tree, with tips
corresponding to individual sequences. Under the infinitely-many-sites mutation
model, mutations are randomly superimposed along the branches of the genealogy,
so that every mutation occurs at a chromosomal site that has not previously
mutated; if a mutation occurs at an interior branch, then all individuals
descending from that branch carry the mutation. The implication is that
observed patterns of molecular variation from this model impose combinatorial
constraints on the hidden state space of genealogies. In particular, observed
molecular variation can be represented in the form of a perfect phylogeny, a
tree structure that fully encodes the mutational differences among sequences.
For a sample of n sequences, a perfect phylogeny might not possess n distinct
leaves, and hence might be compatible with many possible binary tree structures
that could describe the evolutionary relationships among the n sequences. Here,
we investigate enumerative properties of the set of binary ranked and unranked
tree shapes that are compatible with a perfect phylogeny, and hence, the binary
ranked and unranked tree shapes conditioned on an observed pattern of mutations
under the infinitely-many-sites mutation model. We provide a recursive
enumeration of these shapes. We consider both perfect phylogenies that can be
represented as binary and those that are multifurcating. The results have
implications for computational aspects of the statistical inference of
evolutionary parameters that underlie sets of molecular sequences.


  We analyze evolutionary dynamics in a confluent, branching cellular
population, such as in a growing duct, vasculature, or in a branching microbial
colony. We focus on the coarse-grained features of the evolution and build a
statistical model that captures the essential features of the dynamics. Using
simulations and analytic approaches, we show that the survival probability of
strains within the growing population is sensitive to the branching geometry:
Branch bifurcations enhance survival probability due to an overall population
growth (i.e., "inflation"), while branch termination and the small effective
population size at the growing branch tips increase the probability of strain
extinction. We show that the evolutionary dynamics may be captured on a wide
range of branch geometries parameterized just by the branch diameter $N_0$ and
branching rate $b$. We find that the survival probability of neutral cell
strains is largest at an "optimal" branching rate, which balances the effects
of inflation and branch termination. We find that increasing the selective
advantage $s$ of the cell strain mitigates the inflationary effect by
decreasing the average time at which the mutant cell fate is determined. For
sufficiently large selective advantages, the survival probability of the
advantageous mutant decreases monotonically with the branching rate.


  According to the National Cancer Institute, there were 9.5 million
cancer-related deaths in 2018. A challenge in improving treatment is resistance
in genetically unstable cells. The purpose of this study is to evaluate
unsupervised machine learning on classifying treatment-resistant phenotypes in
heterogeneous tumors through analysis of single cell RNA sequencing(scRNAseq)
data with a pipeline and evaluation metrics. scRNAseq quantifies mRNA in cells
and characterizes cell phenotypes. One scRNAseq dataset was analyzed
(tumor/non-tumor cells of different molecular subtypes and patient
identifications). The pipeline consisted of data filtering, dimensionality
reduction with Principal Component Analysis, projection with Uniform Manifold
Approximation and Projection, clustering with nine approaches (Ward, BIRCH,
Gaussian Mixture Model, DBSCAN, Spectral, Affinity Propagation, Agglomerative
Clustering, Mean Shift, and K-Means), and evaluation. Seven models divided
tumor versus non-tumor cells and molecular subtype while six models classified
different patient identification (13 of which were presented in the dataset);
K-Means, Ward, and BIRCH often ranked highest with ~80% accuracy on the tumor
versus non-tumor task and ~60% for molecular subtype and patient ID. An
optimized classification pipeline using K-Means, Ward, and BIRCH models was
evaluated to be most effective for further analysis. In clinical research where
there is currently no standard protocol for scRNAseq analysis, clusters
generated from this pipeline can be used to understand cancer cell behavior and
malignant growth, directly affecting the success of treatment.


  Plankton are effective indicators of environmental change and ecosystem
health in freshwater habitats, but collection of plankton data using manual
microscopic methods is extremely labor-intensive and expensive. Automated
plankton imaging offers a promising way forward to monitor plankton communities
with high frequency and accuracy in real-time. Yet, manual annotation of
millions of images proposes a serious challenge to taxonomists. Deep learning
classifiers have been successfully applied in various fields and provided
encouraging results when used to categorize marine plankton images. Here, we
present a set of deep learning models developed for the identification of lake
plankton, and study several strategies to obtain optimal performances,which
lead to operational prescriptions for users. To this aim, we annotated into 35
classes over 17900 images of zooplankton and large phytoplankton colonies,
detected in Lake Greifensee (Switzerland) with the Dual Scripps Plankton
Camera. Our best models were based on transfer learning and ensembling, which
classified plankton images with 98% accuracy and 93% F1 score. When tested on
freely available plankton datasets produced by other automated imaging tools
(ZooScan, FlowCytobot and ISIIS), our models performed better than previously
used models. Our annotated data, code and classification models are freely
available online.


  For millennia, cyanobacteria have evolved strategies for acclimatization to
dynamic environments- this requires a molecular response system. Phosphate is
essential for cellular integrity and metabolism, yet inorganic phosphate is
often the limiting nutrient for most aquatic environments. This study
characterizes two periplasmic phosphate-binding proteins, namely SphX and SphZ,
which are essential for the functional response of the SphS-SphR signal
transduction system in Synechocystis sp. PCC 6803 and incorporates both into a
model. SphX offers a competitive role so as to ensure inorganic phosphate is
not depleted from the environment at a faster rate than is necessary for
metabolism and, therefore, hinders luxury uptake; whereas, SphZ, encoded by
sll0540, is the auxiliary sensor necessary for the SphS response under
phosphate limitation in tandem with a protein complex association with the
phosphate-specific transport 2 system and the negative regulator, SphU.
Cross-regulation between the pho regulon and the photosynthetic apparatus is
presented as well as a parameterized analysis for low-temperature fluorescence
emission spectroscopy.


  The Reverse Transcription Polymerase Chain Reaction (RTPCR)} test is the
silver bullet diagnostic test to discern COVID infection. Rapid antigen
detection is a screening test to identify COVID positive patients in little as
15 minutes, but has a lower sensitivity than the PCR tests. Besides having
multiple standardized test kits, many people are getting infected and either
recovering or dying even before the test due to the shortage and cost of kits,
lack of indispensable specialists and labs, time-consuming result compared to
bulk population especially in developing and underdeveloped countries.
Intrigued by the parametric deviations in immunological and hematological
profile of a COVID patient, this research work leveraged the concept of
COVID-19 detection by proposing a risk-free and highly accurate Stacked
Ensemble Machine Learning model to identify a COVID patient from communally
available-widespread-cheap routine blood tests which gives a promising
accuracy, precision, recall and F1-score of 100%. Analysis from R-curve also
shows the preciseness of the risk-free model to be implemented. The proposed
method has the potential for large scale ubiquitous low-cost screening
application. This can add an extra layer of protection in keeping the number of
infected cases to a minimum and control the pandemic by identifying
asymptomatic or pre-symptomatic people early.


  We present a novel communication-free algorithm for individual-based
probabilistic neutral biodiversity simulations. The algorithm transforms a
neutral Moran ecosystem model into an embarrassingly parallel problem by
trading off inter-process communication at the cost of some redundant
computation.
  Specifically, by careful design of the random number generator that drives
the simulation, we arrange for evolutionary parent-child interactions to be
modelled without requiring knowledge of the interaction, its participants, or
which processor is performing the computation. Critically, this means that
every individual can be simulated entirely independently. The simulation is
thus fully reproducible irrespective of the number of processors it is
distributed over. With our novel algorithm, a simulation can be (1) split up
into independent batch jobs and (2) simulated across any number of
heterogeneous machines - all without affecting the simulation result.
  We use the Rust programming language to build the extensible and statically
checked simulation package $\texttt{necsim-rust}$. We evaluate our
parallelisation approach by comparing three traditional simulation algorithms
against a CPU and GPU implementation of our Independent algorithm. These
experiments show that as long as some local state is maintained to cull
redundant individuals, our Independent algorithm is as efficient as existing
sequential solutions. The GPU implementation further outperforms all algorithms
on the CPU by a factor ranging from $\sim 2$ to $\sim 80$, depending on the
model parameterisation and the analysis that is performed. Amongst the parallel
algorithms we have investigated, our Independent algorithm provides the only
non-approximate parallelisation strategy that can scale to large simulation
domains.


  Artificial synthesis of DNA molecules is an essential part of the study of
biological mechanisms. The design of a synthetic DNA molecule usually involves
many objectives. One of the important objectives is to eliminate short sequence
patterns that correspond to binding sites of restriction enzymes or
transcription factors. While many design tools address this problem, no
adequate formal solution exists for the pattern elimination problem. In this
work, we present a formal description of the elimination problem and suggest
efficient algorithms that eliminate unwanted patterns and allow optimization of
other objectives with minimal interference to the desired DNA functionality.
Our approach is flexible, efficient, and straightforward, and therefore can be
easily incorporated in existing DNA design tools, making them considerably more
powerful.


  Antibiotic resistance is a topical problem for both humans and animals and
has been the subject of special monitoring for two decades. Several recent
studies, including ours, have shown that this phenomenon is accentuated by the
transfer of certain genetic elements and cross resistance acquisition, the
latter or both resulting from the misuse of these antimicrobial drugs. A more
careful use of antimicrobials, the search for new antibacterial compounds
(including probiotics and phages) are the most recommended alternatives to
overcome this situation. However, tests for modulations of antimicrobial
activity can also play a major role. The main goal of synergy studies is to
assess whether substances with antibacterial properties can improve the
effectiveness of existing antimicrobials or give them a second life against
resistant germs. Moreover, recent studies have demonstrated the ability of
silver nanoparticles and extracts of certain plants to boost the effectiveness
of certain antibiotics. such as ampicillin, benzylpenicillin, cefazolin,
ciprofloxacin, nitrofurantoin, and kanamycin. Yet, from these studies we found
that there was a serious problem with the interpretation of the results when
using the disk method with determination of the increase in fold area. Indeed,
the use of this method firstly requires the determination of the diameter of
inhibition of the antibiotic alone; follow by the determination of the
combination of antibiotic + modulating substance (MS) or extract.....


  Inevitably, almost all cancer patients develop resistance to targeted
therapy. Intratumor heterogeneity (ITH) is a major cause of drug resistance.
Mathematical models that explain experiments quantitatively is useful in
understanding the origin of ITH, which then could be used to explore scenarios
for efficacious therapy. Here, we develop a mathematical model to investigate
ITH in breast cancer by exploiting the observation that HER2+ and HER2- cells
could divide symmetrically or asymmetrically. Our predictions for the evolution
of cell fractions are in quantitative agreement with single-cell experiments.
Remarkably, the colony size of HER2+ cells emerging from a single HER2- cell
(or vice versa), which occurs in about four cell doublings, agrees perfectly
with experimental results, without tweaking any parameter in the model. The
theory quantitatively explains experimental data on the responses of breast
cancer tumor under different treatment protocols. We then used the model to
predict that, not only the order of two drugs, but also the treatment period
for each drug and the tumor cell plasticity could be manipulated to improve the
treatment efficacy. Mathematical models, when integrated with data on patients,
make possible exploration of a broad range of parameters readily, which might
provide insights in devising effective therapies.


  Biological evolution of a population is governed by the fitness landscape,
which is a map from genotype to fitness. However, a fitness landscape depends
on the organisms environment, and evolution in changing environments is still
poorly understood. We study a particular model of antibiotic resistance
evolution in bacteria where the antibiotic concentration is an environmental
parameter and the fitness landscapes incorporate tradeoffs between adaptation
to low and high antibiotic concentration. With evolutionary dynamics that
follow fitness gradients, the evolution of the system under slowly changing
antibiotic concentration resembles the athermal dynamics of disordered physical
systems under external drives. Exploiting this resemblance, we show that our
model can be described as a system with interacting hysteretic elements. As in
the case of the driven disordered systems, adaptive evolution under antibiotic
concentration cycling is found to exhibit hysteresis loops and memory
formation. We derive a number of analytical results for quasistatic
concentration changes. We also perform numerical simulations to study how these
effects are modified under driving protocols in which the concentration is
changed in discrete steps. Our approach provides a general framework for
studying motifs of evolutionary dynamics in biological systems in a changing
environment.


  In recent decades freshwater lakes have seen an increase in human presence. A
common byproduct of this human presence is eutrophication, which readily
results in harmful cyanobacteria blooms. In this work we propose a model that
couples the socio-economic and ecological dynamics related to cyanobacteria
systems. The socio-economic dynamics considers the choices a human population
makes regarding whether or not to mitigate their pollution levels. These
choices are based on various costs related to social ostracism, social norms,
environmental concern and financial burden. The coupled model exhibits bistable
dynamics, with one stable state corresponding to high mitigation efforts and
low CB abundance, and the other to low mitigation efforts and high CB
abundance. Furthermore, we consider social interactions among a network of
lakes and present dynamic outcomes pertaining to various associated costs and
social situations. In each case we show the potential for regime shifts between
levels of cooperation and CB abundance. Social ostracism and pressure are shown
to be driving factors in causing such regime shifts.


  With the increasing affordability and availability of patient data, hospitals
tend to outsource their data to cloud service providers (CSPs) for the purpose
of storage and analytics. However, the concern of data privacy significantly
limits the data owners' choice. In this work, we propose the first solution, to
the best of our knowledge, that allows a CSP to perform efficient
identification of target patients (e.g., pre-processing for a genome-wide
association study - GWAS) over multi-tenant encrypted phenotype data (owned by
multiple hospitals or data owners). We first propose an encryption mechanism
for phenotype data, where each data owner is allowed to encrypt its data with a
unique secret key. Moreover, the ciphertext supports privacy-preserving search
and, consequently, enables the selection of the target group of patients (e.g.,
case and control groups). In addition, we provide a per-query based
authorization mechanism for a client to access and operate on the data stored
at the CSP. Based on the identified patients, the proposed scheme can either
(i) directly conduct GWAS (i.e., computation of statistics about genomic
variants) at the CSP or (ii) provide the identified groups to the client to
directly query the corresponding data owners and conduct GWAS using existing
distributed solutions. We implement the proposed scheme and run experiments
over a real-life genomic dataset to show its effectiveness. The result shows
that the proposed solution is capable to efficiently identify the case/control
groups in a privacy-preserving way.


  We consider the hypothesis that Evolution promotes population-wide genome
patterns that, under randomization, ensures the External Validity of
adaptations across population members. An adaptation is Externally Valid (EV)
if its effect holds under a wide range of population genetic variations. A
prediction following the hypothesis is that pairwise base substitutions in
segregating regions must be 'random' as in Erdos-Renyi-Gilbert random graphs,
but with edge probabilities derived from Experimental-Design concepts. We
demonstrate these probabilities, and consequent mutation rates, in the
full-genomes of 2504 humans, 1135 flowering plants, 1170 flies, 453 domestic
sheep and 1223 brown rats.


  The MinION is a recent-to-market handheld nanopore sequencer. It can be used
to determine the whole genome of a target virus in a biological sample. Its
Read Until feature allows us to skip sequencing a majority of non-target reads
(DNA/RNA fragments), which constitutes more than 99% of all reads in a typical
sample. However, it does not have any on-board computing, which significantly
limits its portability.
  We analyze the performance of a Read Until metagenomic pipeline for detecting
target viruses and identifying strain-specific mutations. We find new sources
of performance bottlenecks (basecaller in classification of a read) that are
not addressed by past genomics accelerators.
  We present SquiggleFilter, a novel hardware accelerated dynamic time warping
(DTW) based filter that directly analyzes MinION's raw squiggles and filters
everything except target viral reads, thereby avoiding the expensive
basecalling step. We show that our 14.3W 13.25mm2 accelerator has 274X greater
throughput and 3481X lower latency than existing GPU-based solutions while
consuming half the power, enabling Read Until for the next generation of
nanopore sequencers.


  The investigation of topological properties of ovarian cancer network (OCN)
and the roles of hubs involved in it by digging the network at various levels
of organization are important to understand how OCN is organized to understand
disease states. The OCN constructed from the experimentally verified genes
exhibits fractal nature in the topological properties of the network and deeply
rooted communities. Also, the network properties at all levels of organization
obey one parameter scaling law which lacks centrality lethality rule. We then
showed that $\langle k\rangle$ can be taken as a scaling parameter, where,
power law exponent can be estimated from the ratio of network diameters,
$\lambda=1+\rho\left[\frac{d_c}{d_k}\right]$. The betweenness centrality shows
two distinct behaviors one shown by high degree hubs and the other by
segregated low degree nodes. The $C_B$ distribution follows power law behavior
with the exponent connected to exponents of distributions of high and low
degree nodes by,
$\theta\sim\left[1+\frac{1}{2}\left(\frac{1}{\epsilon_s}+\frac{1}{\epsilon_h}\right)(\lambda-1)\right]$.
Absence of rich-club formation leads to the missing of a number of attractors
in the network causing formation of weakly tied diverse functional modules to
keep optimal network efficiency. The hubs knockout experiment shows that
provincial hubs take major responsibility to keep network integrity and
organization. The identified key regulators are found to be provincial and
connector hubs. Further, two key regulators, EPCAM and CD44 are found to be
maximally over expressed at various cancer stages (II-IV). They are also
positively correlated with immune infiltrates (CD4+ T cells). Finally, few
potential drugs are identified related to the key regulators.


  Genetic studies of human traits have revolutionized our understanding of the
variation between individuals, and opened the door for numerous breakthroughs
in biology, medicine and other scientific fields. And yet, the ultimate promise
of this area of research is still not fully realized. In this review, we
highlight the major open problems that need to be solved to improve our
understanding of the genetic variation underlying human traits, and by
discussing these challenges provide a primer to the field. Our focus is on
concrete analytical problems, both conceptual and technical in nature. We cover
general issues in genetic studies such as population structure, epistasis and
gene-environment interactions, data-related issues such as ethnic diversity and
rare genetic variants, and specific challenges related to heritability
estimates, genetic association studies and polygenic risk scores. We emphasize
the interconnectedness of these open problems and suggest promising avenues to
address them.


  We analyze an optimal control version of a simple SIR epidemiology model that
includes a partially specified vaccination policy and takes into account
fatigue from protracted application of social distancing measures. The model
assumes demographic (age related) categories, and is otherwise homogeneous.
Maximum capacity for vaccination is exogenous and the authorities select its
allocation to the categories. They can also adopt measures to uniformly
diminish the contact rate between infected and susceptible individuals at an
economic cost. Total or partial immunity is modeled, while the contact rate is
allowed to exhibit seasonality. We apply optimal control methods to minimize a
cost index while guaranteeing that health system capacity is not exceeded. A
reasonable parameter selection leads to policies specifying that vaccination
priority should be given to the category with a higher demand for the limited
health care resources regardless of cost variations among categories.
Vaccination priority reverses if some demographic category has significantly
higher mobility while for some parameter values vaccination alternates among
categories. Optimal social distancing policies exhibit seasonality, and hence
reducing susceptibles below the average spontaneous disease extinction level
does not necessarily lead to a repeal of social distancing.


  Innovation in synthetic biology often still depends on large-scale
experimental trial-and-error, domain expertise, and ingenuity. The application
of rational design engineering methods promise to make this more efficient,
faster, cheaper and safer. But this requires mathematical models of cellular
systems. And for these models we then have to determine if they can meet our
intended target behaviour. Here we develop two complementary approaches that
allow us to determine whether a given molecular circuit, represented by a
mathematical model, is capable of fulfilling our design objectives. We discuss
algebraic methods that are capable of identifying general principles
guaranteeing desired behaviour; and we provide an overview over Bayesian design
approaches that allow us to choose from a set of models, that model which has
the highest probability of fulfilling our design objectives. We discuss their
uses in the context of biochemical adaptation, and then consider how robustness
can and should affect our design approach.


  Protein is linked to almost every life process. Therefore, analyzing the
biological structure and property of protein sequences is critical to the
exploration of life, as well as disease detection and drug discovery.
Traditional protein analysis methods tend to be labor-intensive and
time-consuming. The emergence of deep learning models makes modeling data
patterns in large quantities of data possible. Interdisciplinary researchers
have begun to leverage deep learning methods to model large biological
datasets, e.g. using long short-term memory and convolutional neural network
for protein sequence classification. After millions of years of evolution,
evolutionary information is encoded in protein sequences. Inspired by the
similarity between natural language and protein sequences, we use large-scale
language models to model evolutionary-scale protein sequences, encoding protein
biology information in representation. Significant improvements are observed in
both token-level and sequence-level tasks, demonstrating that our large-scale
model can accurately capture evolution information from pretraining on
evolutionary-scale individual sequences. Our code and model are available at
https://github.com/THUDM/ProteinLM.


  An ecological system with multiple stable equilibria is prone to undergo
catastrophic change or regime shift from one steady-state to another. It should
be noted that, if one of the steady states is an extinction state, the
catastrophic change may lead to extinction. A suitable manual measure may
control the prevention of catastrophic changes of different species from one
equilibrium to another. We consider two stochastic models with linear and
nonlinear harvesting terms. We inspect either density regulation or Allee type
density regulated models [Saha et al., Ecological Modelling, 2013], which have
substantial applications in the herring fish population's viability study. Both
the deterministic models we consider here contain bi-stability under certain
restrictions, and in that case, one of the stable states is the extinction
state. We assume that the dynamical system under consideration is closed, i.e.,
immigration and emigration are absent. The demographic noise is introduced in
the system by substituting an ordinary differential equation with a stochastic
differential equation model, where the birth and death rates of the
deterministic process are used to obtain the instantaneous mean and variance in
the stochastic differential equation. Our study reveals that, the catastrophic
changes can be avoided manually by a suitable choice of handling time that will
eventually help to prevent the sudden extinction of the harvested population.
The entire study is illustrated through the herring population size data
obtained from the Global Population Dynamics Database (GPDD) and simulation
experiment.


  All human diseases involve proteins, yet our current tools to characterize
and quantify them are limited. To better elucidate proteins across space, time,
and molecular composition, we provide provocative projections for technologies
to meet the challenges that protein biology presents. With a broad perspective,
we discuss grand opportunities to transition the science of proteomics into a
more propulsive enterprise. Extrapolating recent trends, we offer potential
futures for a next generation of disruptive approaches to define, quantify and
visualize the multiple dimensions of the proteome, thereby transforming our
understanding and interactions with human disease in the coming decade.


  India's mass vaccination efforts have been slow due to high levels of vaccine
hesitancy. This study uses data from an online discrete choice experiment with
1371 respondents to rigorously examine the factors shaping vaccine preference
in the country. We find that vaccine efficacy, presence of side effects,
protection duration, distance to vaccination centre and vaccination rates
within social networks play a critical role in determining vaccine demand. We
apply a non-parametric model to uncover the heterogeneity in the effects of
these factors. We derive two novel insights from this analysis. First, even
though, on average, domestically developed vaccines are preferred, around 30
percent of the sample favours foreign-developed vaccines. Second, the vaccine
preference of around 15 percent of the sample is highly sensitive to the
presence of side effects and vaccination uptake among their peer group. These
results provide insights into the ongoing policy debate around vaccine adoption
in India.


  The origin of new genes is among the most fundamental processes underlying
genetic innovation. The substrate of new genetic material available defines the
outcomes of evolutionary processes in nature. Historically, the field of
genetic novelty has commonly invoked new mutations at the DNA level to explain
the ways that new genes might originate. In this work, we explore a
fundamentally different source of epistatic interactions that can create new
gene sequences in hybrids. We observe "bursts" of new gene creation in F1
hybrids of D. yakuba and D. santomea, a species complex known to hybridize in
nature. The number of new genes is higher in the gonads than soma. We observe
asymmetry in new gene creation based on the direction of the cross. Greater
numbers of new transcripts form in the testes of F1 male offspring in D.
santomea female x D. yakuba male crosses and greater numbers of new transcripts
appear in ovaries of F1 female offspring of D. yakuba female x D. santomea
male. These loci represent wholly new transcripts expressed in hybrids, but not
in either parental reference strain of the cross. We further observe allelic
activation, where transcripts silenced in one lineage are activated by the
transcriptional machinery of the other genome. These results point to a
fundamentally new model of new gene creation that does not rely on the
formation of new mutations in the DNA. These results suggest that bursts of
genetic novelty can appear in response to hybridization or introgression in a
single generation. Ultimately these processes are expected to contribute to the
substrate of genetic novelty available in nature, with broad impacts on our
understanding new gene formation and on hybrid phenotypes in nature.


  SARS-CoV-2, like any other virus, continues to mutate as it spreads,
according to an evolutionary process. Unlike any other virus, the number of
currently available sequences of SARS-CoV-2 in public databases such as GISAID
is already several million. This amount of data has the potential to uncover
the evolutionary dynamics of a virus like never before. However, a million is
already several orders of magnitude beyond what can be processed by the
traditional methods designed to reconstruct a virus's evolutionary history,
such as those that build a phylogenetic tree. Hence, new and scalable methods
will need to be devised in order to make use of the ever increasing number of
viral sequences being collected.
  Since identifying variants is an important part of understanding the
evolution of a virus, in this paper, we propose an approach based on clustering
sequences to identify the current major SARS-CoV-2 variants. Using a $k$-mer
based feature vector generation and efficient feature selection methods, our
approach is effective in identifying variants, as well as being efficient and
scalable to millions of sequences. Such a clustering method allows us to show
the relative proportion of each variant over time, giving the rate of spread of
each variant in different locations -- something which is important for vaccine
development and distribution. We also compute the importance of each amino acid
position of the spike protein in identifying a given variant in terms of
information gain. Positions of high variant-specific importance tend to agree
with those reported by the USA's Centers for Disease Control and Prevention
(CDC), further demonstrating our approach.


  Metagenomic datasets from pangolin tissue specimens have previously yielded
SARS-related coronaviruses which show high homology in their receptor binding
domain to SARS-CoV-2, suggesting a potential zoonotic source for this feature
of the human virus, possibly via recombination (Liu et al. 2019, Lam et al.
2020, Xiao et al. 2020, Liu et al. 2020). Here we re-examine these published
datasets. We report that only a few pangolin samples were found to contain
coronavirus reads, and even then in low abundance, while other non-pangolin
hosted viruses were present in higher abundance. We also discovered extensive
contamination with human, rodent, and other mammalian gene sequences, which was
a surprising finding. Furthermore, we uncovered a number of pangolin CoV
sequences embedded in standard laboratory cloning vectors, which suggests the
pangolin specimens could have been contaminated with sequences derived from
synthetic biology experiments. Finally, we discover a third pangolin dataset
(He et al. 2022) with low levels of SARSr-CoV sequences and unambiguous
extensive contamination of several pangolin samples. For these reasons, we find
it unlikely that the pangolins in question had a coronavirus infection while
alive, and all current versions of the cited papers claiming a zoonotic
infection of pangolins with a SARS-r CoV require substantial corrections and
should be retracted until such corrections are made.


  Sterile insect technique has been successfully applied in the control of
agricultural pests, however, it has a limited ability to control mosquitoes. A
promising alternative approach is Trojan Y Chromosome strategy, which works by
manipulating the sex ratio of a population through the introduction of
feminized Y Y supermales that guarantee male offspring. A combined Trojan Y
chromosome strategy and sterile insect technique (TYC-SIT) strategy is modeled
with ordinary differential equations that allow the kinetics of the female
population decline of mosquitoes to be evaluated under identical modeling
conditions. The dynamical analysis leads to results on both local and global
stabilities of this combined model. Optimal control analysis is also
implemented to investigate the optimal mechanisms for extinction of mosquitoes.
In particular, the numerical results affirm that the combined TYC-SIT enables
near elimination of mosquitoes. The conclusion has great significance for pest
controls.


  Single species population models and discrete stochastic gene frequency
models are two standards of mathematical biology important for the evolution of
populations. An agent based model is presented which reproduces these models
and then explores where these models agree and disagree under relaxed
specifications. For the population models, the requirement of homogeneous
mixing prevents prediction of extinctions due to local resource depletion.
These models also suggest equilibrium based on attainment of constant
population levels though underlying population characteristics may be nowhere
close to equilibrium. The discrete stochastic gene frequency models assume well
mixed populations at constant levels. The models' predictions for non-constant
populations in strongly oscillating and chaotic regimes are surprisingly good,
only diverging from the ABM at the most chaotic levels.


  The perceived risk and reward for a given situation can vary depending on
resource availability, accumulated wealth, and other extrinsic factors such as
individual backgrounds. Based on this general aspect of everyday life, here we
use evolutionary game theory to model a scenario with randomly perturbed
payoffs in a prisoner's dilemma game. The perception diversity is modeled by
adding a zero-average random noise in the payoff entries and a Monte-Carlo
simulation is used to obtain the population dynamics. This payoff heterogeneity
can promote and maintain cooperation in a competitive scenario where only
defectors would survive otherwise. In this work, we give a step further
understanding the role of heterogeneity by investigating the effects of
quenched disorder in the critical properties of random games. We observe that
payoff fluctuations induce a very slow dynamic, making the cooperation decay
behave as power laws with varying exponents, instead of the usual exponential
decay after the critical point, showing the emergence of a Griffiths phase. We
also find a symmetric Griffiths phase near the defector's extinction point when
fluctuations are present, indicating that Griffiths phases may be frequent in
evolutionary game dynamics and play a role in the coexistence of different
strategies.


  We consider a spatial model of cancer in which cells are points on the
$d$-dimensional torus $\mathcal{T}=[0,L]^d$, and each cell with $k-1$ mutations
acquires a $k$th mutation at rate $\mu_k$. We will assume that the mutation
rates $\mu_k$ are increasing, and we find the asymptotic waiting time for the
first cell to acquire $k$ mutations as the torus volume tends to infinity. This
paper generalizes results on waiting for $k\geq 3$ mutations by Foo, Leder, and
Schweinsberg, who considered the case in which all of the mutation rates
$\mu_k$ were the same. In addition, we find the limiting distribution of the
spatial distances between mutations for certain values of the mutation rates.


  We study the multistationarity for the reaction networks with one-dimensional
stoichiometric subspaces, and we focus on the networks admitting finitely many
positive steady states. We prove that if a network admits multistationarity,
then network has an embedded one-species network with arrow diagram (->,<-) and
another with arrow diagram (<-,->). The inverse is also true if there exist two
reactions in the network such that the subnetwork consisting of the two
reactions admits at least one and finitely many positive steady states. We also
prove that if a network admits at least three positive steady states, then it
contains at least three bi-arrow diagrams. More than that, we completely
characterize the bi-reaction networks that admit at least three positive steady
states.


  The global burden of neurological diseases, the second leading cause of death
after heart dis-eases constitutes one of the major challenges of modern
medicine. Ayurveda, the traditional Indian medicinal systemenrooted in the
Vedic literature and considered as a schema for the holistic management of
health, characterizes various neurological diseases disorders (NDDs) and
prescribes several herbs, formulations, and bio-cleansing regimes for their
care and cure. In this work, we examined neuro-phytoregulatory potential of
34,472 phytochemicals among 3,038 herbs (including their varieties) mentioned
in Ayurveda using network pharmacology approach and found that 45% of these
Ayurvedic phytochemicals (APCs) have regulatory associations with 1,643
approved protein targets. Metabolite interconversion enzymes and protein
modifying enzymes were found to be the major target classes of APCs against
NDDs. The study further suggests that the actions of Ayurvedic herbs in
managing NDDs were majorly via regulating signalling processes, like, G-protein
signaling, acetylcholine signaling, chemokine signaling pathway and GnRH
signaling. A high confidence network specific to 219 pharmaceutically relevant
neuro-phytoregulators (NPRs) from 1,197 Ayurvedic herbs against 102 approved
protein-targets involved in NDDs was developed and analyzed for gaining
mechanistic insights. The key protein targets of NPRs to elicit their
neuro-regulatory effect were highlighted as CYP and TRPA, while estradiol and
melatonin were identified as the NPRs with high multi-targeting ability. 32
herbs enriched in NPRs were identified that include some of the well-known
Ayurvedic neurological recommendations, like, Papaver somniferum, Glycyrrhiza
glabra, Citrus aurantium, Cannabis sativa etc. Herbs enriched in NPRs may be
used as a chemical source library for drug-discovery against NDDs from systems
medicine perspectives.


  Protein-ligand complex structures have been utilised to design benchmark
machine learning methods that perform important tasks related to drug design
such as receptor binding site detection, small molecule docking and binding
affinity prediction. However, these methods are usually trained on only ligand
bound (or holo) conformations of the protein and therefore are not guaranteed
to perform well when the protein structure is in its native unbound
conformation (or apo), which is usually the conformation available for a newly
identified receptor. A primary reason for this is that the local structure of
the binding site usually changes upon ligand binding. To facilitate solutions
for this problem, we propose a dataset called APObind that aims to provide apo
conformations of proteins present in the PDBbind dataset, a popular dataset
used in drug design. Furthermore, we explore the performance of methods
specific to three use cases on this dataset, through which, the importance of
validating them on the APObind dataset is demonstrated.


  Cancer cell population dynamics often exhibit remarkably replicable,
universal laws despite their underlying heterogeneity. Mechanistic explanations
of universal cell population growth remain partly unresolved to this day,
whereby population feedback between the microscopic and mesoscopic
configurations can lead to macroscopic growth laws. We here present a
unification under density-dependent birth events via contact inhibition. We
consider five classical tumor growth laws: exponential, generalized logistic,
Gompertz, radial growth, and fractal growth, which can be seen as
manifestations of a single microscopic model. Our theory is substantiated by
agent based simulations and can explain growth curve differences in
experimental data from in vitro cancer cell population growth. Thus, our
framework offers a possible explanation for the large number of mean-field laws
that can adequately capture seemingly unrelated cancer or microbial growth
dynamics.


  The emergence of SARS-CoV-2 variants of concern endangers the long-term
control of COVID-19, especially in countries with limited genomic surveillance.
In this work, we explored genomic drivers of contagion in Chile. We sequenced
3443 SARS-CoV-2 genomes collected between January and July 2021, where the
Gamma (P.1), Lambda (C.37), Alpha (B.1.1.7), B.1.1.348, and B.1.1 lineages were
predominant. Using a Bayesian model tailored for limited genomic surveillance,
we found that Lambda and Gamma variants' reproduction numbers were about 5% and
16% larger than Alpha's, respectively. We observed an overabundance of
mutations in the Spike gene, strongly correlated with the variant's
transmissibility. Furthermore, the variants' mutational signatures featured a
breakpoint concurrent with the beginning of vaccination (mostly CoronaVac, an
inactivated virus vaccine), indicating an additional putative selective
pressure. Thus, our work provides a reliable method for quantifying novel
variants' transmissibility under subsampling (as newly-reported Delta,
B.1.617.2) and highlights the importance of continuous genomic surveillance.


  The COVID-19 pandemic has posed significant challenges in modeling its
complex epidemic transmissions, infection and contagion, which are very
different from known epidemics. The challenges in quantifying COVID-19
complexities include effectively modeling its process and data uncertainties.
The uncertainties are embedded in implicit and high-proportional undocumented
infections, asymptomatic contagion, social reinforcement of infections, and
various quality issues in the reported data. These uncertainties become even
more apparent in the first two months of the COVID-19 pandemic, when the
relevant knowledge, case reporting and testing were all limited. Here we
introduce a novel hybrid approach Susceptible-Undocumented infected-Documented
infected-Recovered (SUDR) model. First, SUDR (1) characterizes and
distinguishes Undocumented (U) and Documented (D) infections commonly seen
during COVID-19 incubation periods and asymptomatic infections. Second, SUDR
characterizes the probabilistic density of infections by capturing exogenous
processes. Lastly, SUDR approximates the density likelihood of COVID-19
prevalence over time by incorporating Bayesian inference into SUDR. Different
from existing COVID-19 models, SUDR characterizes the undocumented infections
during unknown transmission processes. To capture the uncertainties of temporal
transmission and social reinforcement during COVID-19 contagion, the
transmission rate is modeled by a time-varying density function of undocumented
infectious cases. By sampling from the mean-field posterior distribution with
reasonable priors, SUDR handles the randomness, noise and sparsity of COVID-19
observations widely seen in the public COVID-19 case data. The results
demonstrate a deeper quantitative understanding of the above uncertainties, in
comparison with classic SIR, time-dependent SIR, and probabilistic SIR models.


  First proposed as an empirical rule over half a century ago, the Richards
growth equation has been frequently invoked in population modeling and pandemic
forecasting. Central to this model is the advent of a fractional exponent
$\gamma$, typically fitted to the data. While various motivations for this
non-analytical form have been proposed, it is still considered foremost an
empirical fitting procedure. Here, we find that Richards-like growth laws
emerge naturally from generic analytical growth rules in a distributed
population, upon inclusion of {\bf (i)} migration (spatial diffusion) amongst
different locales, and {\bf (ii)} stochasticity in the growth rate, also known
as "seascape noise." The latter leads to a wide (power-law) distribution in
local population number that, while smoothened through the former, can still
result in a fractional growth law for the overall population. This
justification of the Richards growth law thus provides a testable connection to
the distribution of constituents of the population.


  Direct reciprocity is a well-known mechanism that could explain how
cooperation emerges and prevails in an evolving population. Numerous prior
researches have studied the emergence of cooperation in multiplayer games.
However, most of them use numerical or experimental methods, not theoretical
analysis. This lack of theoretical works on the evolution of cooperation is due
to the high complexity of calculating payoffs. In this paper, we propose a new
method, namely, the state-clustering method to calculate the long-term payoffs
in repeated games. Using this method, in an $n$-player repeated game, the
computing complexity is reduced from $O(2^n)$ to $O(n^2)$, which makes it
possible to compute a large-scale repeated game's payoff. We explore the
evolution of cooperation in both infinitely and finitely repeated public goods
games as an example to show the effectiveness of our method. In both cases, we
find that when the synergy factor is sufficiently large, the increasing number
of participants in a game is detrimental to the evolution of cooperation. Our
work provides a theoretical approach to study the evolution of cooperation in
repeated multiplayer games.


  The unexpected Ebola virus outbreak in West Africa in 2014 involving the
Zaire ebolavirus made clear that other regions outside Central Africa, its
previously documented niche, were at risk of future epidemics. The complex
transmission cycle and a lack of epidemiological data make mapping areas at
risk of the disease challenging. We used a Geographic Information System-based
multicriteria evaluation (GIS-MCE), a knowledge-based approach, to identify
areas suitable for Ebola virus spillover to humans in regions of Guinea, Congo
and Gabon where Ebola viruses already emerged. We identified environmental,
climatic and anthropogenic risk factors and potential hosts from a literature
review. Geographical data layers, representing risk factors, were combined to
produce suitability maps of Ebola virus spillover at the landscape scale. Our
maps show high spatial and temporal variability in the suitability for Ebola
virus spillover at a fine regional scale. Reported spillover events fell in
areas of intermediate to high suitability in our maps, and a sensitivity
analysis showed that the maps produced were robust. There are still important
gaps in our knowledge about what factors are associated with the risk of Ebola
virus spillover. As more information becomes available, maps produced using the
GIS-MCE approach can be easily updated to improve surveillance and the
prevention of future outbreaks.


  Wolbachia is a natural bacterium that can infect mosquitoes and reduce their
ability to transmit mosquito-borne diseases, such as dengue fever, Zika, and
chikungunya. Field trials and modeling studies have shown that the fraction of
infection among the mosquitoes must exceed a threshold level for the infection
to persist. To capture this threshold, it is critical to consider the spatial
heterogeneity in the distributions of the infected and uninfected mosquitoes,
which is created by the local release of the infected mosquitoes. We develop
and analyze partial differential equation (PDE) models to study the invasion
dynamics of Wolbachia infection among mosquitoes in the field. Our
reaction-diffusion-type models account for both the complex vertical
transmission and the spatial mosquito dispersion. We characterize the threshold
for a successful invasion, which is a bubble-shaped profile, called the
"critical bubble". The critical bubble is optimal in its release size compared
to other spatial profiles in a one-dimensional landscape. The fraction of
infection near the release center is higher than the threshold level for the
corresponding homogeneously mixing ODE models. We show that the proposed
spatial models give rise to the traveling waves of Wolbachia infection when
above the threshold. We quantify how the threshold condition and traveling-wave
velocity depend on the diffusion coefficients and other model parameters.
Numerical studies for different scenarios are presented to inform the design of
release strategies.


  Contagious respiratory diseases, such as COVID-19, depend on sufficiently
prolonged exposures for the successful transmission of the underlying pathogen.
It is important for organizations to evaluate the efficacy of interventions
aiming at mitigating viral transmission among their personnel. We have
developed a operational risk assessment simulation framework that couples a
spatial agent-based model of movement with a SIR epidemiological model to
assess the relative risks of different intervention strategies. By applying our
model on MIT's STATA building, we assess the impacts of three possible
dimensions of intervention: one-way vs unrestricted movement, population size
allowed onsite, and frequency of leaving designated work location for breaks.
We find that there is no significant impact made by one-way movement
restrictions over unrestricted movement. Instead, we find that a combination of
lowering the number of individuals admitted below the current recommendations
and advising individuals to reduce the frequency at which they leave their
workstations lowers the likelihood of highly connected individuals within the
contact networks that emerge, which in turn lowers the overall risk of
infection. We discover three classes of possible interventions based on their
epidemiological effects. By assuming a direct relationship between data on
secondary attack rates and transmissibility in the SIR model, we compare
relative infection risk of four respiratory diseases, MERS, SARS, COVID-19, and
Measles, within the simulated area, and recommend appropriate intervention
guidelines.


  A mathematical model of Zika virus transmission incorporating human movement
between rural areas and nearby forests is presented to investigate the role of
human movement in the spread of Zika virus infections in human and mosquito
populations. Proportions of both susceptible and infected humans living in
rural areas are assumed to move to nearby forest areas. Direct, indirect and
vertical transmission routes are incorporated for all populations. Mathematical
analysis of the proposed model has been presented. The analysis starts with
normalizing the proposed model. Positivity and boundedness of solutions to the
normalized model have been then addressed. The basic reproduction number has
been calculated using the next generation matrix method and its relation to the
three routes of disease transmission has been presented. The sensitivity
analysis of the basic reproduction number to all model parameters has been
investigated. The analysis also includes existence and stability of disease
free and endemic equilibrium points. Bifurcation analysis has been also carried
out. Finally, numerical solutions to the normalized model have been obtained to
confirm the theoretical results and to demonstrate the impact of human movement
in the disease transmission in human and mosquito populations.


  The ongoing COVID-19 pandemic highlights the essential role of mathematical
models in understanding the spread of the virus along with a quantifiable and
science-based prediction of the impact of various mitigation measures. Numerous
types of models have been employed with various levels of success. This leads
to the question of what kind of a mathematical model is most appropriate for a
given situation. We consider two widely used types of models: equation-based
models (such as standard compartmental epidemiological models) and agent-based
models. We assess their performance by modeling the spread of COVID-19 on the
Hawaiian island of Oahu under different scenarios. We show that when it comes
to information crucial to decision making, both models produce very similar
results. At the same time, the two types of models exhibit very different
characteristics when considering their computational and conceptual complexity.
Consequently, we conclude that choosing the model should be mostly guided by
available computational and human resources.


  Deep generative chemistry models emerge as powerful tools to expedite drug
discovery. However, the immense size and complexity of the structural space of
all possible drug-like molecules pose significant obstacles, which could be
overcome with hybrid architectures combining quantum computers with deep
classical networks. As the first step toward this goal, we built a compact
discrete variational autoencoder (DVAE) with a Restricted Boltzmann Machine
(RBM) of reduced size in its latent layer. The size of the proposed model was
small enough to fit on a state-of-the-art D-Wave quantum annealer and allowed
training on a subset of the ChEMBL dataset of biologically active compounds.
Finally, we generated 2331 novel chemical structures with medicinal chemistry
and synthetic accessibility properties in the ranges typical for molecules from
ChEMBL. The presented results demonstrate the feasibility of using already
existing or soon-to-be-available quantum computing devices as testbeds for
future drug discovery applications.


  Machine learning methods can detect complex relationships between variables,
but usually do not exploit domain knowledge. This is a limitation because in
many scientific disciplines, such as systems biology, domain knowledge is
available in the form of graphs or networks, and its use can improve model
performance. We need network-based algorithms that are versatile and applicable
in many research areas. In this work, we demonstrate subnetwork detection based
on multi-modal node features using a novel Greedy Decision Forest with inherent
interpretability. The latter will be a crucial factor to retain experts and
gain their trust in such algorithms. To demonstrate a concrete application
example, we focus on bioinformatics, systems biology and particularly
biomedicine, but the presented methodology is applicable in many other domains
as well. Systems biology is a good example of a field in which statistical
data-driven machine learning enables the analysis of large amounts of
multi-modal biomedical data. This is important to reach the future goal of
precision medicine, where the complexity of patients is modeled on a system
level to best tailor medical decisions, health practices and therapies to the
individual patient. Our proposed approach can help to uncover disease-causing
network modules from multi-omics data to better understand complex diseases
such as cancer.


  Adenocarcinoma and squamous cell carcinoma constitute approximately 40% and
30% of all lung cancer subtypes, respectively, and display broad heterogeneity
in terms of clinical and molecular responses to therapy. Molecular subtyping
has enabled precision medicine to overcome these challenges and provide
significant biological insights to predict prognosis and improve clinical
decision making. Over the past decade, conventional ML algorithms and DL-based
CNNs have been espoused for the classification of cancer subtypes from gene
expression datasets. However, these methods are potentially biased toward
identification of cancer biomarkers. Recently proposed transformer-based
architectures that leverage the self-attention mechanism encode high throughput
gene expressions and learn representations that are computationally complex and
parametrically expensive. However, compared to the datasets for natural
language processing applications, gene expression consists of several hundreds
of thousands of genes from a limited number of observations, making it
difficult to efficiently train transformers for bioinformatics applications.
Hence, we propose an end-to-end deep learning approach, Gene Transformer, which
addresses the complexity of high-dimensional gene expression with a multi-head
self-attention module by identifying relevant biomarkers across multiple cancer
subtypes without requiring feature selection as a prerequisite for the current
classification algorithms. The proposed architecture achieved an overall
improved performance for all evaluation metrics and had fewer misclassification
errors than the commonly used traditional classification algorithms. The
classification results show that Gene Transformer can be an efficient approach
for classifying cancer subtypes, indicating that any improvement in deep
learning models in computational biology can also be reflected well in this
domain.


  In this study, we develop and analyze a nested multi-scale model for COVID
-19 disease that integrates within-host scale and between-host scale
sub-models. First, the well-posedness of the multi-scale model is discussed,
followed by the stability analysis of the equilibrium points. The disease-free
equilibrium point is shown to be globally asymptotically stable for $R_0 < 1$.
When $R_0$ exceeds unity, a unique infected equilibrium exists, and the system
is found to undergo a forward (trans-critical) bifurcation at $R_0=1$. Two
parameter heat plots are also done to find the parameter combinations for which
the equilibrium points are stable. The parameters $\beta, \pi$ and $\Lambda$
are found to be most sensitive to $R_0$. The influence of within-host sub-model
parameter on the between-host sub-model variables is numerically illustrated.
The spread of infection in a community is shown to be influenced by within-host
level sub-model parameters, such as the production of viral particles by
infected cells $(\alpha)$, the clearance rate of infected cells by the immune
system $(x)$, and the clearance rate of viral particles by the immune system
$(y)$. The comparative effectiveness of the three health interventions
(antiviral drugs, immunomodulators, and generalized social distancing) for
COVID-19 infection was examined using the effective reproductive number $R_E$
as an indicator of the effectiveness of the interventions. The results suggest
that a combined strategy of antiviral drugs, immunomodulators and generalized
social distancing would be the best strategy to implement to contain the spread
of infection in the community.


  This paper is an exploratory study of two epidemiological questions on a
worldwide basis. How fast is the disease spreading? Are the restrictions
(especially mobility restrictions) for people bring the expected effect? To
answer the first question, we propose a tool for estimating the reproduction
number of epidemic (the number of secondary infections $R_t$) based on the SEIR
model and compare it with an non-model $R_t$ estimation. To measure the $R_t$
of COVID-19 for different countries, real-time data on coronavirus daily cases
of infections, recoveries, deaths are retrieved from the Center for Systems
Science and Engineering (CSSE) at Johns Hopkins University. To assess the
effectiveness of mobility restrictions for the COVID-19 pandemic in 2020, the
correlations between the $R_t$ and people's mobility (based on the Apple
mobility index) are presented. The correlations were considered for 12
countries and for most of them, the correlations are negative. This shows a
delay in the implementation of mobility restrictions - the countries imposed
them in response to growth of new COVID-19 cases, rather than preventively.


  Trans-splicing of leader sequences onto the 59ends of mRNAs is a widespread
phenomenon in protozoa, nematodes and some chordates. Using parallel sequencing
we have developed a method to simultaneously map 59splice sites and analyze the
corresponding gene expression profile, that we term spliced leader trapping
(SLT). The method can be applied to any organism with a sequenced genome and
trans-splicing of a conserved leader sequence. We analyzed the expression
profiles and splicing patterns of bloodstream and insect forms of the parasite
Trypanosoma brucei. We detected the 59splice sites of 85% of the annotated
protein-coding genes and, contrary to previous reports, found up to 40% of
transcripts to be differentially expressed. Furthermore, we discovered more
than 2500 alternative splicing events, many of which appear to be
stage-regulated. Based on our findings we hypothesize that alternatively
spliced transcripts present a new means of regulating gene expression and could
potentially contribute to protein diversity in the parasite. The entire dataset
can be accessed online at TriTrypDB or through: http://splicer.unibe.ch/.


  Conformational sampling of biomolecules using molecular dynamics simulations
often produces large amount of high dimensional data that makes it difficult to
interpret using conventional analysis techniques. Dimensionality reduction
methods are thus required to extract useful and relevant information. Here we
devise a machine learning method, Gaussian mixture variational autoencoder
(GMVAE) that can simultaneously perform dimensionality reduction and clustering
of biomolecular conformations in an unsupervised way. We show that GMVAE can
learn a reduced representation of the free energy landscape of protein folding
with highly separated clusters that correspond to the metastable states during
folding. Since GMVAE uses a mixture of Gaussians as the prior, it can directly
acknowledge the multi-basin nature of protein folding free-energy landscape. To
make the model end-to-end differentialble, we use a Gumbel-softmax
distribution. We test the model on three long-timescale protein folding
trajectories and show that GMVAE embedding resembles the folding funnel with
folded states down the funnel and unfolded states outer in the funnel path.
Additionally, we show that the latent space of GMVAE can be used for kinetic
analysis and Markov state models built on this embedding produce folding and
unfolding timescales that are in close agreement with other rigorous dynamical
embeddings such as time independent component analysis (TICA).


  It is generally accepted that "diversity" is associated with success in
evolutionary algorithms. However, diversity is a broad concept that can be
measured and defined in a multitude of ways. To date, most evolutionary
computation research has measured diversity using the richness and/or evenness
of a particular genotypic or phenotypic property. While these metrics are
informative, we hypothesize that other diversity metrics are more strongly
predictive of success. Phylogenetic diversity metrics are a class of metrics
popularly used in biology, which take into account the evolutionary history of
a population. Here, we investigate the extent to which 1) these metrics provide
different information than those traditionally used in evolutionary
computation, and 2) these metrics better predict the long-term success of a run
of evolutionary computation. We find that, in most cases, phylogenetic metrics
behave meaningfully differently from other diversity metrics. Moreover, our
results suggest that phylogenetic diversity is indeed a better predictor of
success.


  Evolutionary game theory is a powerful mathematical framework to study how
intelligent individuals adjust their strategies in collective interactions. It
has been widely believed that it is impossible to unilaterally control players'
payoffs in games, since payoffs are jointly determined by all players. Until
recently, a class of so-called zero-determinant strategies are revealed, which
enables a player to make a unilateral payoff control over her partners in
two-action repeated games with a constant continuation probability. The
existing methods, however, lead to the curse of dimensionality when the
complexity of games increases. In this paper, we propose a new mathematical
framework to study ruling strategies (with which a player unilaterally makes a
linear relation rule on players' payoffs) in repeated games with an arbitrary
number of actions or players, and arbitrary continuation probability. We
establish an existence theorem of ruling strategies and develop an algorithm to
find them. In particular, we prove that strict Markov ruling strategy exists
only if either the repeated game proceeds for an infinite number of rounds, or
every round is repeated with the same probability. The proposed mathematical
framework also enables the search of collaborative ruling strategies for an
alliance to control outsiders. Our method provides novel theoretical insights
into payoff control in complex repeated games, which overcomes the curse of
dimensionality.


  The quarantine of identified close contacts has been vital to reducing
transmission rates and averting secondary infection risk before symptom onset
and by asymptomatic cases. The effectiveness of this contact tracing strategy
to mitigate transmission is sensitive to the adherence to quarantines, which
may be lower for longer quarantine periods or in vaccinated populations (where
perceptions of risk are reduced). This study develops a simulation model to
evaluate contact tracing strategies based on the sequential testing of
identified contacts after exposure as an alternative to quarantines, in which
contacts are isolated only after confirmation by a positive test. The analysis
considers different number and types of tests (PCR and lateral flow antigen
tests (LFA)) to identify the cost-effective testing policies that minimize the
expected infecting days post-exposure considering different levels of testing
capacity. This analysis suggests that even a limited number of tests can be
effective at reducing secondary infection risk: two LFA tests (with optimal
timing) avert infectiousness at a level that is comparable to 14-day quarantine
with 80-90% adherence, or equivalently, 7-9 day quarantine with full adherence
(depending on the sensitivity of the LFA test). Adding a third test (PCR or
LFA) reaches the efficiency of a 14-day quarantine with 90-100% adherence.
These results are robust to the exposure dates of the contact, test sensitivity
of LFA and alternative models of viral load evolution, which suggests that
simple testing rules can be effective for improving contact tracing in settings
where strict quarantine adherence is difficult to implement.


  The Mediterranean Sea is one of the most heavily invaded marine regions. This
work focuses on the dynamics and potential policy options for ballast
water-mediated nonindigenous species to the Mediterranean. Specifically, we (1)
estimated port risks in years 2012, 2015, and 2018, (2) identified hub ports
that connect many clusters, and (3) evaluated four regulatory scenarios. The
risk results show that Gibraltar, Suez, and Istanbul remained high-risk ports
from 2012-2018, and they served as hub ports that connected several spread
clusters. With policy scenario analysis, we found that regulating the high-risk
hub ports can disproportionately reduce the overall risk to the Mediterranean:
the average risk to all ports was reduced by 5-10% by regulating one high-risk
hub port, while the average risk to all ports was only reduced by 0.2% by
regulating one average-risk Mediterranean port. We also found that only
regulating high-risk ports cannot reduce their risks effectively.


  Genomic surveillance of SARS-CoV-2 has been instrumental in tracking the
spread and evolution of the virus during the pandemic. The availability of
SARS-CoV-2 molecular sequences isolated from infected individuals, coupled with
phylodynamic methods, have provided insights into the origin of the virus, its
evolutionary rate, the timing of introductions, the patterns of transmission,
and the rise of novel variants that have spread through populations. Despite
enormous global efforts of governments, laboratories, and researchers to
collect and sequence molecular data, many challenges remain in analyzing and
interpreting the data collected. Here, we describe the models and methods
currently used to monitor the spread of SARS-CoV-2, discuss long-standing and
new statistical challenges, and propose a method for tracking the rise of novel
variants during the epidemic.


  Cardiovascular disease, especially heart failure is one of the major health
hazard issues of our time and is a leading cause of death worldwide.
Advancement in data mining techniques using machine learning (ML) models is
paving promising prediction approaches. Data mining is the process of
converting massive volumes of raw data created by the healthcare institutions
into meaningful information that can aid in making predictions and crucial
decisions. Collecting various follow-up data from patients who have had heart
failures, analyzing those data, and utilizing several ML models to predict the
survival possibility of cardiovascular patients is the key aim of this study.
Due to the imbalance of the classes in the dataset, Synthetic Minority
Oversampling Technique (SMOTE) has been implemented. Two unsupervised models
(K-Means and Fuzzy C-Means clustering) and three supervised classifiers (Random
Forest, XGBoost and Decision Tree) have been used in our study. After thorough
investigation, our results demonstrate a superior performance of the supervised
ML algorithms over unsupervised models. Moreover, we designed and propose a
supervised stacked ensemble learning model that can achieve an accuracy,
precision, recall and F1 score of 99.98%. Our study shows that only certain
attributes collected from the patients are imperative to successfully predict
the surviving possibility post heart failure, using supervised ML algorithms.


  Ovarian Cancer (OC) is type of female reproductive malignancy which can be
found among young girls and mostly the women in their fertile or reproductive.
There are few number of cysts are dangerous and may it cause cancer. So, it is
very important to predict and it can be from different types of screening are
used for this detection using Transvaginal Ultrasonography (TVUS) screening. In
this research, we employed an actual datasets called PLCO with TVUS screening
and three machine learning (ML) techniques, respectively Random Forest KNN, and
XGBoost within three target variables. We obtained a best performance from this
algorithms as far as accuracy, recall, f1 score and precision with the
approximations of 99.50%, 99.50%, 99.49% and 99.50% individually. The AUC score
of 99.87%, 98.97% and 99.88% are observed in these Random Forest, KNN and XGB
algorithms .This approach helps assist physicians and suspects in identifying
ovarian risks early on, reducing ovarian malignancy-related complications and
deaths.


  A nonlinear cross-diffusion epidemic with a time-dependent
Susceptible-Infected-Recovered-Died system is proposed in this paper. This
system is derived from kinetic theory model by multiscale approach, which leads
to an equivalent system coupled the microscopic and macroscopic equations.
Subsequently, numerical investigations to design asymptotic preserving scheme
property is developed and validated by various numerical tests. Finally, the
numerical computational results of the proposed system are discussed in two
dimensional space using the finite volume method.


  Estimating animal abundance and density are fundamental goals of many
wildlife monitoring programs. Camera trapping has become an increasingly
popular tool to achieve these monitoring goals due to recent advances in
modeling approaches and the capacity to simultaneously collect data on multiple
species. However, estimating the density of unmarked populations continues to
be problematic due to the difficulty in implementing complex modeling
approaches, low precision of estimates, and absence of rigor in testing of
model assumptions and their influence on results. Here, we describe a novel
approach that uses still image camera traps to estimate animal density without
the need for individual identification, based on the Time spent In Front of the
Camera (TIFC). Using results from a large-scale multi-species monitoring
program with nearly 3,000 cameras deployed over six years in Alberta, Canada,
we provide a reproducible methodology to estimate parameters and we test key
assumptions of the TIFC model. We compare moose (Alces alces) density estimates
from aerial surveys and TIFC, including incorporating correction factors for
known TIFC assumption violations. The resulting corrected TIFC density
estimates are comparable to aerial density estimates. We discuss the
limitations of the TIFC method and areas needing further investigation,
including the need for long-term monitoring of assumption violations and the
number of cameras necessary to provide precise estimates. Despite the
challenges of assumption violations and high measurement error, cameras and the
TIFC method can provide useful alternative or complementary animal density
estimates for multi-species monitoring when compared to traditional monitoring
methods.


  In this paper, a nonuniform size modified Poisson-Boltzmann ion channel
(nuSMPBIC) model is presented as a nonlinear system of an electrostatic
potential and multiple ionic concentrations. It mixes nonlinear algebraic
equations with a Poisson boundary value problem involving Dirichlet-Neumann
mixed boundary value conditions and a membrane surface charge density to
reflect the effects of ion sizes and membrane charges on electrostatics and
ionic concentrations. To overcome the difficulties of strong singularities and
exponential nonlinearities, it is split into three submodels with a solution of
Model 1 collecting all the singular points and Models 2 and 3 much easier to
solve numerically than the original nuSMPBIC model. A damped two-block
iterative method is then presented to solve Model 3, along with a novel
modified Newton iterative scheme for solving each related nonlinear algebraic
system. To this end, an effective nuSMPBIC finite element solver is derived and
then implemented as a program package that works for an ion channel protein
with a three-dimensional molecular structure and a mixture solution of multiple
ionic species. Numerical results for a voltage-dependent anion channel (VDAC)
in a mixture of four ionic species demonstrate a fast convergence rate of the
damped two-block iterative method, the high performance of the software
package, and the importance of considering nonuniform ion sizes. Moreover, the
nuSMPBIC model is validated by the anion selectivity property of VDAC.


  Population inhomogeneity, in the variation of the individual social contact
networks and the individual infectious-recovery rates, renders the dynamics of
infectious disease spreading uncertain. As a consequence the overlaying
economical production network with its proper collaboration components is to
extent impacted unpredictably. Our model proposes a `vector logistic' dynamical
approach to SIS dynamics in a social contact network interacting with its
economic capacity network. The probabilistic interpretation of the graph state
in the vector logistic description provides a method to assess the effect of
mean and variance of the infected on the production capacity and allows the
strategic planning of social connectivity regulation. The impact of the
epidemic mean effects and fluctuations on the production capacity is assessed
according `cumulative',`majority' and `fragility' proxy measures.


  A probabilistic approach to the epidemic evolution on realistic
social-contact networks allows for characteristic differences among subjects,
including the individual number and structure of social contacts, and the
heterogeneity of the infection and recovery rates according to age or medical
preconditions. Within our probabilistic Susceptible-Infectious-Removed (SIR)
model on social-contact networks, we evaluate the `infection load' or
`activation margin' of various control scenarios; by confinement, by
vaccination, and by their combination. We compare the epidemic burden for
subpopulations which apply competing or co-operative control strategies. The
simulation experiments are conducted on randomised social-contact graphs that
are designed to exhibit realistic person-person contact characteristics and
which follow near `homogeneous' or `block-localised' subpopulation spreading.
The scalarization method is used for the multi-objective optimization problem
in which both the infection load is minimized and the extent to which each
subpopulation's control strategy preference ranking is adhered to is maximized.
We obtain the compounded payoff matrices for two subpopulations which impose
contrasting control strategies, each according to their proper ranked control
strategy preferences. The Nash equilibria, according to each subpopulation's
compounded objective, and according to their proper ranking intensity, are
discussed. Finally, the interaction effects of the control strategies are
discussed and related to the type of spreading of the two subpopulations.


  Here we describe the coupled standardisation of two complementary
fluorescence imaging techniques and apply it to liquid-liquid phase separated
condensates formed from an EGFP fluorescent reporter of Flowering Control Locus
A (FCA), a protein that associates with chromosomal DNA in plants during
epigenetic regulation of the flowering process. First, we use home-built
single-molecule Slimfield microscopy to establish a fluorescent protein
standard. This sample comprises live yeast cells expressing Mig1 protein, a
metabolic regulator which localises to the nucleus under conditions of high
glucose, fused to the same type of EGFP label as for the FCA fusion construct.
Then we employ commercial confocal AiryScan microscopy to study the same
standard. Finally, we demonstrate how to quantify FCA-EGFP nuclear condensates
in intact root tips at rapid timescales and apply this calibration. This method
is a valuable approach to obtaining single-molecule precise stoichiometry and
copy number estimates of protein condensates that are integrated into the
chromosome architecture of plants, using confocal instrumentation that lacks de
facto single-molecule detection sensitivity.


  Far from being a passive information store, the genome is a mechanically
dynamic and diverse system in which torsion and tension fluctuate and combine
to determine structure and help regulate gene expression. Much of this
mechanical perturbation is due to molecular machines such as topoisomerases
which must stretch and twist DNA as part of various functions including DNA
repair and replication. While the broad-scale mechanical response of nucleic
acids to tension and torsion is well characterized, detail at the single base
pair level is beyond the limits of even super-resolution imaging. Here, we
present a straightforward, flexible, and extensible umbrella-sampling protocol
to twist and stretch nucleic acids in silico using the popular biomolecular
simulation package Amber -- though the principles we describe are applicable
also to other packages such as GROMACS. We discuss how to set up the simulation
system, decide forcefields and solvation models, and equilibrate. We then
introduce the torsionally-constrained stretching protocol, and finally we
present some analysis techniques we have used to characterize structural motif
formation. Rather than define forces or fictional pseudoatoms, we instead
define a fixed translation of specified atoms between each umbrella sampling
step, which allows comparison with experiment without needing to estimate
applied forces by simply using the fractional end-to-end displacement as a
comparison metric. We hope that this easy to implement solution will be
valuable for interrogating optical and magnetic tweezers data on nucleic acids
at base pair resolution.


  More than 198 million cases of severe acute respiratory syndrome coronavirus
2 (SARS-CoV-2) has been reported that result in no fewer than 4.2 million
deaths globally. The rapid spread of the disease coupled with the lack of
specific registered drugs for its treatment pose a great challenge that
necessitate the development of therapeutic agents from a variety of sources. In
this study, we employed an in-silico method to screen natural compounds with a
view to identify inhibitors of the human transmembrane protease serine type 2
(TMPRSS2). The activity of this enzyme is essential for viral access into the
host cells via angiotensin-converting enzyme 2 (ACE-2). Inhibiting the activity
of this enzyme is therefore highly crucial for preventing viral fusion with
ACE-2 thus shielding SARS-CoV-2 infectivity. 3D model of TMPRSS2 was
constructed using I-TASSER, refined by GalaxyRefine, validated by Ramachandran
plot server and overall model quality was checked by ProSA. 95 natural
compounds from microalgae were virtually screened against the modeled protein
that led to the identification 17 best leads capable of binding to TMPRSS2 with
a good binding score comparable, greater or a bit lower than that of the
standard inhibitor (camostat). Physicochemical properties, ADME (absorption,
distribution, metabolism, excretion) and toxicity analysis revealed top 4
compounds including the reference drug with good pharmacokinetic and
pharmacodynamic profiles. These compounds bind to the same pocket of the
protein with a binding energy of -7.8 kcal/mol, -7.6 kcal/mol, -7.4 kcal/mol
and -7.4 kcal/mol each for camostat, apigenin, catechin and epicatechin
respectively. This study shed light on the potential of microalgal compounds
against SARS-CoV-2. In vivo and invitro studies are required to developed
SARS-CoV-2 drugs based on the structures of the compounds identified in this
study.


  Population-wide vaccination is critical for containing the SARS-CoV-2
(Covid-19) pandemic when combined with restrictive and prevention measures. In
this study, we introduce SAIVR, a mathematical model able to forecast the
Covid-19 epidemic evolution during the vaccination campaign. SAIVR extends the
widely used Susceptible-Infectious-Removed (SIR) model by considering the
Asymptomatic (A) and Vaccinated (V) compartments. The model contains several
parameters and initial conditions that are estimated by employing a
semi-supervised machine learning procedure. After training an unsupervised
neural network to solve the SAIVR differential equations, a supervised
framework then estimates the optimal conditions and parameters that best fit
recent infectious curves of 27 countries. Instructed by these results, we
performed an extensive study on the temporal evolution of the pandemic under
varying values of roll-out daily rates, vaccine efficacy, and a broad range of
societal vaccine hesitancy/denial levels. The concept of herd immunity is
questioned by studying future scenarios which involve different vaccination
efforts and more infectious Covid-19 variants.


  Determinants of COVID-19 clinical severity are commonly assessed by
transverse or longitudinal studies of the fatality counts. However, the
fatality counts depend both on disease clinical severity and transmissibility,
as more infected also lead to more deaths. Moreover, fatality counts (and
related measures such as Case Fatality Rate) are dynamic quantities, as they
appear with a delay to infections, while different geographic regions generally
belong to different points on the epidemics curve. Instead, we use
epidemiological modeling to propose a disease severity measure, which accounts
for the underlying disease dynamics. The measure corresponds to the ratio of
population averaged mortality and recovery rates (m/r). It is independent of
the disease transmission dynamics (i.e., the basic reproduction number) and has
a direct mechanistic interpretation. We use this measure to assess demographic,
medical, meteorological and environmental factors associated with the disease
severity. For this, we employ an ecological regression study design and analyze
different US states during the first disease outbreak. Principal Component
Analysis, followed by univariate and multivariate analyses based on machine
learning techniques, is used for selecting important predictors. Without using
prior knowledge from clinical studies, we recover significant predictors known
to influence disease severity, in particular age, chronic diseases, and racial
factors. Additionally, we identify long-term pollution exposure and population
density as not widely recognized (though for the pollution previously
hypothesized) predictors of the disease severity. Overall, the proposed measure
is useful for inferring severity determinants of COVID-19 and other infectious
diseases, and the obtained results may aid a better understanding of COVID-19
risks.


  The novel coronavirus SARS-CoV-2 has resulted in a global pandemic with
worldwide 6-digital infection rates and thousands death tolls daily. Enormeous
effords are undertaken to achieve high coverage of immunization in order to
reach herd immunity to stop spreading of SARS-CoV-2 infection. Several
SARS-CoV-2 vaccines, based either on mRNA, viral vectors, or inactivated
SARS-CoV-2 virus have been approved and are being applied worldwide. However,
recently increased numbers of normally very rare types of thromboses associated
with thrombocytopenia have been reported in particular in the context of the
adenoviral vector vaccine ChAdOx1 nCoV-19 from Astra Zeneca. While statistical
prevalence of these side effects seem to correlate with this particular vaccine
type, i.e. adenonoviral vector based vaccines, the exact molecular mechanisms
are still not clear. The present review summarizes current data and hypotheses
for molecular and cellular mechanisms into one integrated hypothesis indicating
that coagulopathies, including thromboses, thrombocytopenia and other related
side effects are correlated to an interplay of the two components in the
vaccine, i.e. the spike antigen and the adenoviral vector, with the innate and
immune system which under certain circumstances can imitate the picture of a
limited COVID-19 pathological picture.


  The soybean aphid, Aphis glycines (Hemiptera: Aphididae), is an invasive pest
that can cause severe yield loss to soybeans in the northcentral United States.
A tactic to counter this pest is the use of aphid-resistant soybean varieties.
However, the occurrence of virulent biotypes can alter plant physiology and
impair the use of this management strategy. Soybean aphids can alter soybean
physiology primarily by two mechanisms, feeding facilitation and the obviation
of resistance, favoring subsequent colonization by additional conspecifics. We
developed a non-local, differential equation population model, to explore the
dynamics of these biological mechanisms on soybean plants co-infested with
virulent and avirulent aphids. We then use demographic parameters from
laboratory experiments to perform numerical simulations via the model. These
simulations successfully mimic various aphid dynamics observed in the field.
Our model showed an increase in colonization of virulent aphids increases the
likelihood that aphid-resistance is suppressed, subsequently increasing the
survival of avirulent aphids, producing an indirect, positive interaction
between the biotypes. These results suggest the potential for a "within plant"
refuge that could contribute to the sustainable use of aphid resistant
soybeans.


  Natural and synthetic phenazines are widely used in biomedical sciences. In
dehydrogenase histochemistry, phenazine methosulfate (PMS) is applied as a
redox reagent for coupling reduced coenzymes to the reduction of tetrazolium
salts into colored formazans. PMS is also currently used for cytotoxicity and
viability assays of cell cultures using sulfonated tetrazoliums. Under UV (340
nm) excitation, aqueous solutions of the cationic PMS show green fluorescence
({\lambda}em: 526 nm), whereas the reduced hydrophobic derivative
(methyl-phenazine, MPH) shows blue fluorescence ({\lambda}em: 465 nm). Under UV
(365 nm) excitation, cultured cells (LM2, IGROV-1, BGC-1, and 3T3-L1
adipocytes) treated with PMS (5 ug/mL, 30 min) showed cytoplasmic granules with
bright blue fluorescence, which correspond to lipid droplets labeled by the
lipophilic methyl-phenazine. After formaldehyde fixation blue-fluorescing
droplets could be stained with oil red O. Interestingly, PMS-treated 3T3-L1
adipocytes observed under UV excitation 24 h after labeling showed large lipid
droplets with a weak green emission within a diffuse pale blue-fluorescing
cytoplasm, whereas a strong green emission was observed in small lipid
droplets. This fluorescence change from blue to green indicates that
reoxidation of methyl-phenazine to PMS can occur. Regarding cell uptake and
labeling mechanisms, QSAR models predict that the hydrophilic PMS is not
significantly membrane-permeant, so most PMS reduction is expected to be
extracellular and associated with a plasma membrane NAD(P)H reductase. Once
formed, the lipophilic and blue-fluorescing methyl-phenazine enters live cells
and mainly accumulates in lipid droplets. Overall, the results reported here
indicate that PMS is an excellent fluorescent probe to investigate labeling and
redox dynamics of lipid droplets in cultured cells.


  For several decades optical tweezers have proven to be an invaluable tool in
the study and analysis of a myriad biological responses and applications.
However, as every tool, it can have undesirable or damaging effects upon the
very sample it is helping to study. In this review the main negative effects of
optical tweezers upon biostructures and living systems will be presented. Three
are the main areas on which the review will focus: linear optical excitation
within the tweezers, non-linear photonic effects, and thermal load upon the
sampled volume. Additional information is provided on negative mechanical
effects of optical traps on biological structures. Strategies to avoid or, in
the least, minimize these negative effects will be introduced. Finally, all
these effects, undesirable for the most, can have positive applications under
the right conditions. Some hints in this direction will also be discussed.


  In this paper we introduce an agent-based model together with a particle
filter approach for studying the spread of COVID-19. Investigations are
performed on the metropolis of Tokyo, but other cities, regions or countries
could have been equally chosen. A novel method for evaluating the effective
reproduction number is one of the main outcome of our approach. Other unknown
parameters and unknown populations are also evaluated. Uncertain quantities, as
for example the ratio of symptomatic / asymptomatic agents, are tested and
discussed, and the stability of our computations is examined. Detailed
explanations are provided for the model and for the assimilation process.