  DNA sequencing is revolutionising the field of medicine. DNA sequencers, the
machines which perform DNA sequencing, have evolved from the size of a fridge
to that of a mobile phone over the last two decades. The cost of sequencing a
human genome also has reduced from billions of dollars to hundreds of dollars.
Despite these improvements, DNA sequencers output hundreds or thousands of
gigabytes of data that must be analysed on computers to discover meaningful
information with biological implications. Unfortunately, the analysis
techniques have not kept the pace with rapidly improving sequencing
technologies. Consequently, even today, the process of DNA analysis is
performed on high-performance computers, just as it was a couple of decades
ago. Such high-performance computers are not portable. Consequently, the full
utility of an ultra-portable sequencer for sequencing in-the-field or at the
point-of-care is limited by the lack of portable lightweight analytic
techniques. This thesis proposes computer architecture-aware optimisation of
DNA analysis software. DNA analysis software is inevitably convoluted due to
the complexity associated with biological data. Modern computer architectures
are also complex. Performing architecture-aware optimisations requires the
synergistic use of knowledge from both domains, (i.e, DNA sequence analysis and
computer architecture). This thesis aims to draw the two domains together. In
this thesis, gold-standard DNA sequence analysis workflows are systematically
examined for algorithmic components that cause performance bottlenecks.
Identified bottlenecks are resolved through architecture-aware optimisations at
different levels, i.e., memory, cache, register and processor. The optimised
software tools are used in complete end-to-end analysis workflows and their
efficacy is demonstrated by running on prototypical embedded systems.


  Amid the ongoing COVID-19 pandemic, public health authorities and the general
population are striving to achieve a balance between safety and normalcy. Ever
changing conditions call for the development of theory and simulation tools to
finely describe multiple strata of society while supporting the evaluation of
"what-if" scenarios. Particularly important is to assess the effectiveness of
potential testing approaches and vaccination strategies. Here, an agent-based
modeling platform is proposed to simulate the spreading of COVID-19 in small
towns and cities, with a single-individual resolution. The platform is
validated on real data from New Rochelle, NY -- one of the first outbreaks
registered in the United States. Supported by expert knowledge and informed by
reported data, the model incorporates detailed elements of the spreading within
a statistically realistic population. Along with pertinent functionality such
as testing, treatment, and vaccination options, the model accounts for the
burden of other illnesses with symptoms similar to COVID-19. Unique to the
model is the possibility to explore different testing approaches -- in
hospitals or drive-through facilities -- and vaccination strategies that could
prioritize vulnerable groups. Decision making by public authorities could
benefit from the model, for its fine-grain resolution, open-source nature, and
wide range of features.


  Graph theoretical analyses of nervous systems usually omit the aspect of
connection polarity, due to data insufficiency. The chemical synapse network of
Caenorhabditis elegans is a well-reconstructed directed network, but the signs
of its connections are yet to be elucidated. Here, we present the gene
expression-based sign prediction of the ionotropic chemical synapse connectome
of C. elegans (3,638 connections and 20,589 synapses total), incorporating
available presynaptic neurotransmitter and postsynaptic receptor gene
expression data for three major neurotransmitter systems. We made predictions
for more than two-thirds of these chemical synapses and observed an
excitatory-inhibitory (E:I) ratio close to 4:1 which was found similar to that
observed in many real-world networks. Our open source tool
(http://EleganSign.linkgroup.hu) is simple but efficient in predicting
polarities by integrating neuronal connectome and gene expression data.


  Liquid-liquid phase separation (LLPS) is important to control a wide range of
reactions from gene expression to protein degradation in a cell-sized space. To
bring a better understanding of the compatibility of such phase-separated
structures with protein synthesis, we study emergent LLPS in a cell-free
transcription-translation (TXTL) reaction. When the TXTL reaction composed of
many proteins is concentrated, the uniformly mixed state becomes unstable, and
membrane-less phases form spontaneously. This LLPS droplet formation is induced
when the TXTL reaction is enclosed in water-in-oil emulsion droplets, in which
water evaporates from the surface. As the emulsion droplets shrink, smaller
LLPS droplets appear inside the emulsion droplets and coalesce into large
phase-separated domains that partition the localization of synthesized reporter
proteins. The presence of PEG in the TXTL reaction is important not only for
versatile cell-free protein synthesis but also for the formation of two large
domains capable of protein partitioning. Our results may shed light on the
dynamic interplay of LLPS formation and cell-free protein synthesis toward the
construction of synthetic organelles.


  Though membrane trafficking of cell junction proteins has been studied
extensively for more than two decades, the accumulated knowledge remains
fragmentary. The goal of this review is to synthesize published studies on the
membrane trafficking of the five major junction transmembrane proteins:
claudins, occludin, and junction adhesion molecules (JAMs) in tight junctions;
cadherins and nectins in adherens junctions; to identify underlying common
mechanisms; to highlight their functional consequences on barrier function; and
to identify knowledge gaps. Clathrin-mediated endocytosis appears to be the
main, but not exclusive, mode of internalization. Caveolin-mediated endocytosis
and macropinocytosis are employed less frequently. PDZ-domain binding is the
predominant mode of interaction between junction protein cytoplasmic tails and
scaffold proteins. It is shared by claudins, the largest family of junction
integral proteins, by junction adhesion molecules A, B, and C, and by the three
nectins. All eight proteins are destined to either recycling via Rab4/Rab11
GTPases or to degradation. The sorting mechanisms that underlie the specificity
of their endocytic pathways and determine their fates are not fully known. New
data is presented to introduce an emerging role of junction-associated scaffold
proteins in claudin membrane trafficking.


  Background: Advanced biological techniques have helped produce more
insightful findings on the genetic etiology of infertility that may lead to
better management of the condition. This review provides an update on genes
predisposing to syndromic and nonsyndromic infertility.
  Main body: The review identified 65 genes linked with infertility and
infertility-related disorders. These genes regulate fertility. However,
mutational loss of the functions of the genes predisposes to infertility.
Twenty-three (23) genes representing 35% were linked with syndromic
infertility, while 42 genes (65%) cause nonsyndromic infertility. Of the 42
nonsyndromic genes, 26 predispose to spermatogenic failure and sperm
morphological abnormalities, 11 cause ovarian failures, and 5 cause sex
reversal and puberty delay. Overall, 31 genes (48%) predispose to male
infertility, 15 genes (23%) cause female infertility, and 19 genes (29%)
predispose to both. The common feature of male infertility was spermatogenic
failure and sperm morphology abnormalities, while ovarian failure has been the
most frequently reported among infertile females. The mechanisms leading to
these pathologies are gene-specific, which, if targeted in the affected, may
lead to improved treatment.
  Conclusions: Mutational loss of the functions of some genes involved in the
development and maintenance of fertility may predispose to syndromic or
nonsyndromic infertility via gene-specific mechanisms. A treatment procedure
that targets the affected gene(s) in individuals expressing infertility may
lead to improved treatment.


  The spatial organization of the DNA in the cell nucleus plays an important
role for gene regulation, DNA replication, and genomic integrity. Through the
development of chromosome conformation capture experiments (such as 3C, 4C,
Hi-C) it is now possible to obtain the contact frequencies of the DNA at the
whole-genome level. In this paper, we study the problem of reconstructing the
3D organization of the genome from such whole-genome contact frequencies. A
standard approach is to transform the contact frequencies into noisy distance
measurements and then apply semidefinite programming (SDP) formulations to
obtain the 3D configuration. However, neglected in such reconstructions is the
fact that most eukaryotes including humans are diploid and therefore contain
two copies of each genomic locus. We prove that the 3D organization of the DNA
is not identifiable from distance measurements derived from contact frequencies
in diploid organisms. In fact, there are infinitely many solutions even in the
noise-free setting. We then discuss various additional biologically relevant
and experimentally measurable constraints (including distances between
neighboring genomic loci and higher-order interactions) and prove
identifiability under these conditions. Furthermore, we provide SDP
formulations for computing the 3D embedding of the DNA with these additional
constraints and show that we can recover the true 3D embedding with high
accuracy from both noiseless and noisy measurements. Finally, we apply our
algorithm to real pairwise and higher-order contact frequency data and show
that we can recover known genome organization patterns.


  Working in the multi-type Galton-Watson branching-process framework we
analyse the spread of a pandemic via a general multi-type random contact graph.
Our model consists of several communities, and takes, as input, parameters that
outline the contacts between individuals in distinct communities. Given these
parameters, we determine whether there will be an outbreak and if yes, we
calculate the size of the giant connected component of the graph, thereby,
determining the fraction of the population of each type that would be infected
before it ends. We show that the pandemic spread has a natural evolution
direction given by the Perron-Frobenius eigenvector of a matrix whose entries
encode the average number of individuals of one type expected to be infected by
an individual of another type. The corresponding eigenvalue is the basic
reproduction number of the pandemic. We perform numerical simulations that
compare homogeneous and heterogeneous spread graphs and quantify the difference
between them. We elaborate on the difference between herd immunity and the end
of the pandemic and the effect of countermeasures on the fraction of infected
population.


  We present a generic epidemic model with stochastic parameters, in which the
dynamics self-organize to a critical state with suppressed exponential growth.
More precisely, the dynamics evolve into a quasi-steady-state, where the
effective reproduction rate fluctuates close to the critical value one, as
observed for different epidemics. The main assumptions underlying the model are
that the rate at which each individual becomes infected changes stochastically
in time with a heavy-tailed steady state. The critical regime is characterized
by an extremely long duration of the epidemic. Its stability is analyzed both
numerically and analytically.


  We present the GeneScore, a concept of feature reduction for Machine Learning
analysis of biomedical data. Using expert knowledge, the GeneScore integrates
different molecular data types into a single score. We show that the GeneScore
is superior to a binary matrix in the classification of cancer entities from
SNV, Indel, CNV, gene fusion and gene expression data. The GeneScore is a
straightforward way to facilitate state-of-the-art analysis, while making use
of the available scientific knowledge on the nature of molecular data features
used.


  It has been recently discovered that the measles virus can wipe out the
adaptive immune system, destroying B lymphocytes and reducing the diversity of
non-specific B cells of the infected host. In particular, this implies that
previously acquired immunization from vaccination or direct exposition to other
pathogens could be erased in a phenomenon named "immune amnesia", whose effects
can become particularly worrisome given the actual rise of anti-vaccination
movements. Here we present the first attempt to incorporate immune amnesia into
standard models of epidemic spreading. In particular, we analyze diverse
variants of a model that describes the spreading of two concurrent pathogens
causing measles and another generic disease: the SIR-IA model. Analytical and
computational studies confirm that immune amnesia can indeed have important
consequences for epidemic spreading, significantly altering the vaccination
coverage required to reach herd-immunity for concurring infectious diseases.
More specifically, we uncover the existence of novel propagating and endemic
phases which are induced by immune amnesia, that appear both in fully-connected
and more structured networks, such as random networks and power-law
degree-distributed ones. In particular, the transitions from a quiescent state
into these novel phases can become rather abrupt in some cases that we
specifically analyze. Furthermore, we discuss the meaning and consequences of
our results and their relation with, e.g., immunization strategies, together
with the possibility that explosive types of transitions may emerge, making
immune-amnesia effects particularly dramatic. This work opens the door to
further developments and analyses of immune amnesia effects, contributing, more
generally, to the theory of interacting epidemics on complex networks.


  In this paper, we investigate compact ultrametric measure spaces which form a
subset $\mathcal{U}^w$ of the collection of all metric measure spaces
$\mathcal{M}^w$. Similar as for the ultrametric Gromov-Hausdorff distance on
the collection of ultrametric spaces $\mathcal{U}$, we define ultrametric
versions of two metrics on $\mathcal{U}^w$, namely of Sturm's distance of order
$p$ and of the Gromov-Wasserstein distance of order $p$. We study the basic
topological and geometric properties of these distances as well as their
relation and derive for $p=\infty$ a polynomial time algorithm for their
calculation. Further, several lower bounds for both distances are derived and
some of our results are generalized to the case of finite ultra-dissimilarity
spaces.


  Decoding the genome confers the capability to predict characteristics of the
organism(phenotype) from DNA (genotype). We describe the present status and
future prospects of genomic prediction of complex traits in humans. Some highly
heritable complex phenotypes such as height and other quantitative traits can
already be predicted with reasonable accuracy from DNA alone. For many
diseases, including important common conditions such as coronary artery
disease, breast cancer, type I and II diabetes, individuals with outlier
polygenic scores (e.g., top few percent) have been shown to have 5 or even 10
times higher risk than average. Several psychiatric conditions such as
schizophrenia and autism also fall into this category. We discuss related
topics such as the genetic architecture of complex traits, sibling validation
of polygenic scores, and applications to adult health, in vitro fertilization
(embryo selection), and genetic engineering.


  A great issue of discussion of an infectious disease is its basic
reproduction number R0, which provides an estimation of the contagiousness of
the disease. When R0 > 1, the disease spread will potentially lead to an
outbreak, such that of the ongoing COVID-19 pandemics. During the evolution of
an outbreak, various non-pharmaceutical interventions are employed, the impact
of which is frequently assessed by the reduction that they introduce to the
effective reproduction number Rt; reduction below 1 is an indication of
eventual dying out of the disease spread. Motivated by the fact that R0
essentially expresses the stability of the disease-free equilibrium, in this
work, Rt was examined in the view of timescale analysis. It was shown that
during the evolution of the COVID-19 outbreak in Italy, when various
interventions were in place, Rt had a clear relation with the explosive
timescale characterizing the dynamics of the outbreak. In particular, it is
shown that the existence of an explosive timescale during the progression of
the epidemics implies Rt > 1, while its absence implies Rt < 1. In addition, as
this timescale converges/diverges with the immediately slowest one, Rt
approaches to/withdraws from its threshold value 1. These results suggest that
timescale analysis can be utilized for the assessment of the impact of various
interventions, since it reflects the insight provided by the effective
reproduction number, without being hindered by the selection of the population
model, nor the parameter estimation process followed for model calibration.


  We use categorical methods to define a new flavor of Petri nets where
transitions can only fire a limited number of times, specified by a quantity
that we call mana. We do so with chemistry in mind, looking at ways of
modelling the behavior of chemical reactions that depend on enzymes to work. We
prove that such nets can be either obtained as a result of a comonadic
construction, or by enriching them with extra information encoded into a
functor. We then use a well-established categorical result to prove that the
two constructions are equivalent, and generalize them to the case where the
firing of some transitions can "regenerate" the mana of others. This allows us
to represent the action of catalysts and also of biochemical processes where
the byproducts of some chemical reaction are exactly the enzymes that another
reaction needs to work.


  We may attempt to encapsulate what we know about a physical system by a model
structure, $S$. This collection of related models is defined by parametric
relationships between system features; say observables (outputs), unobservable
variables (states), and applied inputs. Each parameter vector in some parameter
space is associated with a completely specified model in $S$. Before choosing a
model in $S$ to predict system behaviour, we must estimate its parameters from
system observations. Inconveniently, multiple models (associated with distinct
parameter estimates) may approximate data equally well. Yet, if these equally
valid alternatives produce dissimilar predictions of unobserved quantities,
then we cannot confidently make predictions. Thus, our study may not yield any
useful result. We may anticipate the non-uniqueness of parameter estimates
ahead of data collection by testing $S$ for structural global identifiability
(SGI). Here we will provide an overview of the importance of SGI, some
essential theory and distinctions, and demonstrate these in testing some
examples.


  The possibility of targeting the causal genes along with the mechanisms of
pathogenically complex diseases has led to numerous studies on the genetic
etiology of some diseases. In particular, studies have added more genes to the
list of type 1 diabetes mellitus (T1DM) suspect genes, necessitating an update
for the interest of all stakeholders. Therefore this review articulates T1DM
suspect genes and their pathophysiology. Notable electronic databases,
including Medline, Scopus, PubMed, and Google-Scholar were searched for
relevant information. The search identified over 73 genes suspected in the
pathogenesis of T1DM, with human leukocyte antigen, insulin gene, and cytotoxic
T lymphocyte-associated antigen 4 accounting for most of the cases. Mutations
in these genes, along with environmental factors, may produce a defective
immune response in the pancreas, resulting in \b{eta}-cell autoimmunity,
insulin deficiency, and hyperglycemia. The mechanisms leading to these cellular
reactions are gene-specific and, if targeted in diabetic individuals, may lead
to improved treatment. Medical practitioners are advised to formulate treatment
procedures that target these genes in patients with T1DM.


  Single gene mutations have been implicated in the pathogenesis of a form of
diabetes mellitus (DM) known as the maturity-onset diabetes of the young
(MODY). However, there are diverse opinions on the suspect genes and
pathophysiology, necessitating the need to review and communicate the genes to
raise public awareness. We used the Google search engine to retrieve relevant
information from reputable sources such as PubMed and Google Scholar. We
identified 14 classified MODY genes as well as three new and unclassified genes
linked with MODY. These genes are fundamentally embedded in the beta cells, the
most common of which are HNF1A, HNF4A, HNF1B, and GCK genes. Mutations in these
genes cause beta-cell dysfunction, resulting in decreased insulin production
and hyperglycemia. MODY genes have distinct mechanisms of action and phenotypic
presentations compared with type 1 and type 2 DM and other forms of DM.
Healthcare professionals are therefore advised to formulate drugs and treatment
based on the causal genes rather than the current generalized treatment for all
types of DM. This will increase the effectiveness of diabetes drugs and
treatment and reduce the burden of the disease.


  Motility is a ubiquitous feature of microbial life on Earth, and is widely
regarded as a promising biosignature candidate. In the search for motile
organisms, it is therefore valuable to have rough estimates for the number of
such microbes that one may expect to find in a given area or volume. In this
work, we explore this question by employing a simple theoretical model that
takes into account the amount of free energy available in a given environment
and the energetic cost of motility. We present heuristic upper bounds for the
average biomass density and the number density of motile lifeforms for the
Martian subsurface and the ocean of Enceladus by presuming that the motile
microbes in question derive their energy from methanogenesis. We consequently
demonstrate that the resultant densities of motile organisms might be
potentially comparable to, or much lower than, the total microbial densities
documented in various extreme environments on Earth.


  The human proteome is enriched in proteins that do not fold into a stable 3D
structure. These intrinsically disordered proteins (IDPs) spontaneously
fluctuate between a large number of configurations in their native form.
Remarkably, the disorder does not lead to dysfunction as with denatured folded
proteins. In fact, unlike denatured proteins, recent evidences strongly suggest
that multiple biological functions stem from such structural plasticity. Here,
focusing on the nanoscopic length-scale, we review the latest advances in IDP
research and discuss some of the future directions in this highly promising
field.


  2-colored best match graphs (2-BMGs) form a subclass of sink-free
bi-transitive graphs that appears in phylogenetic combinatorics. There, 2-BMGs
describe evolutionarily most closely related genes between a pair of species.
They are explained by a unique least resolved tree (LRT). Introducing the
concept of support vertices we derive an $O(|V|+|E|\log^2|V|)$-time algorithm
to recognize 2-BMGs and to construct its LRT. The approach can be extended to
also recognize binary-explainable 2-BMGs with the same complexity. An empirical
comparison emphasizes the efficiency of the new algorithm.


  Self-organized spatial patterns of vegetation are frequent in drylands and,
because pattern shape correlates with water availability, they have been
suggested as important indicators of ecosystem health. However, the mechanisms
underlying pattern emergence remain unclear. Some theories hypothesize that
patterns could result from a water-mediated scale-dependent feedback (SDF)
whereby interactions favoring plant growth dominate at short distances and
growth-inhibitory interactions dominate in the long range. However, we know
little about how the presence of a focal plant affects the fitness of its
neighbors as a function of the inter-individual distance, which is expected to
be highly ecosystem-dependent. This lack of empirical knowledge and system
dependency challenge the relevance of SDF as a unifying theory for vegetation
pattern formation. Assuming that plant interactions are always inhibitory and
only their intensity is scale-dependent, alternative theories also recover the
typical vegetation patterns observed in nature. Importantly, although these
alternative hypotheses lead to visually indistinguishable patterns, they
predict contrasting desertification dynamics, which questions the potential use
of vegetation patterns as ecosystem-state indicators. To help resolve this
issue, we first review existing theories for vegetation self-organization and
their conflicting predictions about desertification dynamics. Second, we
discuss potential empirical tests via manipulative experiments to identify
pattern-forming mechanisms and link them to specific desertification dynamics.
A comprehensive view of models, the mechanisms they intend to capture, and
experiments to test them in the field will help to better understand both how
patterns emerge and improve predictions on the fate of the ecosystems where
they form.


  The Tilman model of consumers competing for resources has some aspects that
appear counterintuitive. Within the standard Tilman model for a single consumer
and a single resource, when a very efficient consumer rapidly eats all
available food, the resource density becomes zero, but if there is no food, how
can the consumer survive? The paradox can be lifted by realising that on the
short term indeed rapid consumption may lead to starvation and a decline in the
consumer population, but in the long term a finite resource and consumer
density remain. A single consumer living on two essential nutrients leaves the
density of the non-limiting nutrient above its critical level, so what is done
with the extra food? We explain that this does not imply extra food is
consumed, just for fun, but that in fact the food is left untouched because
there is no use for it. For a model with two consumers competing for a single
nutrient one of the consumers will eventually disappear, even if the food
economics of the surviving one is worse and there is still much food wasted.
Both are inherent aspects of the model, and the paradox can be avoided if the
difference between the consumers is small, in which case it will take very long
to reach equilibrium. We argue that in general an extended stability analysis,
in which not only the asymptotically stable state is considered, but also the
unstable steady states and all time scales involved in the transient dynamics,
can help in avoiding apparent contradictions in ecological models, or accepting
them.


  Engineering simple, artificial models of living cells allows synthetic
biologists to study cellular functions under well-controlled conditions.
Reconstituting multicellular behaviors with synthetic cell-mimics is still a
challenge because it requires efficient communication between individual
compartments in large populations. This protocol presents a microfluidic method
to produce large quantities of cell-mimics with highly porous, stable and
chemically modifiable polymer membranes that can be programmed on demand with
nucleus-like DNA-hydrogel compartments for gene expression. We describe
expression of genes encoded in the hydrogel compartment and communication
between neighboring cell-mimics through diffusive protein signals.


  Background: Many studies have been conducted on the genetic and epigenetic
etiology of gestational diabetes mellitus (GDM) in the last two decades because
of the diseases increasing prevalence and role in the global diabetes mellitus
(DM) explosion. An update on the genetic and epigenetic etiology of GDM then
becomes imperative to better understand and stem the rising incidence of the
disease. This review, therefore, articulated GDM candidate genes and their
pathophysiology for the awareness of stakeholders. Main body (genetic and
epigenetic etiology, GDM): The search discovered 83 GDM candidate genes, of
which TCF7L2, MTNR1B, CDKAL1, IRS1, and KCNQ1 are the most prevalent. Certain
polymorphisms of these genes can modulate beta-cell dysfunction, adiposity,
obesity, and insulin resistance through several mechanisms. Environmental
triggers such as diets, pollutants, and microbes may also cause epigenetic
changes in these genes, resulting in a loss of insulin-boosting and glucose
metabolism functions. Early detection and adequate management may resolve the
condition after delivery; otherwise, it will progress to maternal type 2
diabetes mellitus (T2DM) and fetal configuration to future obesity and DM. This
shows that GDM is a strong risk factor for T2DM and, in rare cases, type 1
diabetes mellitus (T1DM) and maturity-onset diabetes of the young (MODY). This
further shows that GDM significantly contributes to the rising incidence and
burden of DM worldwide and its prevention may reverse the trend. Conclusion:
Mutations and epigenetic changes in certain genes are strong risk factors for
GDM. For affected individuals with such etiologies, medical practitioners
should formulate drugs and treatment procedures that target these genes and
their pathophysiology.


  Since the end of 2019, COVID-19 has significantly affected the lives of
people around the world. Towards the end of 2020, several COVID-19 vaccine
candidates with relatively high efficacy have been reported in the final phase
of clinical trials. Vaccines have been considered as critical tools for opening
up social and economic activities, thereby lessening the impact of this disease
on the society. This paper presents a simulation of COVID-19 spread using
modified Susceptible-Infected-Removed (SIR) model under vaccine intervention in
several localities of Malaysia, i.e. those cities or states with high
relatively COVID-19 cases such as Kuala Lumpur, Penang, Sabah, and Sarawak. The
results show that at different vaccine efficacy levels (0.75, 0.85, and 0.95),
the curves of active infection vary slightly, indicating that vaccines with
efficacy above 0.75 would produce the herd immunity required to level the
curves. In addition, disparity is significant between implementing or not
implementing a vaccination program. Simulation results also show that lowering
the reproduction number, R0 is necessary to keep the infection curve flat
despite vaccination. This is due to the assumption that vaccination is mostly
carried out gradually at the assumed fixed rate. The statement is based on our
simulation results with two values of R0: 1.1 and 1.2, indicative of reduction
of R0 by social distancing. The lower R0 shows a smaller peak amplitude about
half the value simulated with R0=1.2. In conclusion, the simulation model
suggests a two-pronged strategy to combat the COVID-19 pandemic in Malaysia:
vaccination and compliance with standard operating procedure issued by the
World Health Organization (e.g. social distancing).


  Identification of genes that initiate cell anomalies and cause cancer in
humans is among the important fields in the oncology researches. The mutation
and development of anomalies in these genes are then transferred to other genes
in the cell and therefore disrupt the normal functionality of the cell. These
genes are known as cancer driver genes (CDGs). Various methods have been
proposed for predicting CDGs, most of which based on genomic data and based on
computational methods. Therefore, some researchers have developed novel
bioinformatics approaches. In this study, we propose an algorithm, which is
able to calculate the effectiveness and strength of each gene and rank them by
using the gene regulatory networks and the stochastic analysis of regulatory
linking structures between genes. To do so, firstly we constructed the
regulatory network using gene expression data and the list of regulatory
interactions. Then, using biological and topological features of the network,
we weighted the regulatory interactions. After that, the obtained regulatory
interactions weight was used in interaction structure analysis process.
Interaction analysis was achieved using two separate Markov chains on the
bipartite graph obtained from the main graph of the gene network. To do so, the
stochastic approach for link-structure analysis has been implemented. The
proposed algorithm categorizes higher-ranked genes as driver genes. The
efficiency of the proposed algorithm, regarding the F-measure value and number
of identified driver genes, was compared with 23 other computational and
network-based methods.


  Mesoscale molecular assemblies on the cell surface, such as cilia and
filopodia, integrate information, control transport and amplify signals.
Synthetic devices mimicking these structures could sensitively monitor these
cellular functions and direct new ones. The challenges in creating such
devices, however are that they must be integrated with cells in a precise
kinetically controlled process and a device's structure and its precisely
structured cell interface must then be maintained during active cellular
function. Here we report the ability to integrate synthetic micro-scale
filaments, DNA nanotubes, into a cell's architecture by anchoring them by their
ends to specific receptors on the surfaces of mammalian cells. These filaments
can act as shear stress meters: how anchored nanotubes bend at the cell surface
quantitatively indicates the magnitude of shear stresses between 0-2 dyn per
cm2, a regime important for cell signaling. Nanotubes can also grow while
anchored to cells, thus acting as dynamic components of cells. This approach to
cell surface engineering, in which synthetic biomolecular assemblies are
organized within existing cellular architecture, could make it possible to
build new types of sensors, machines and scaffolds that can interface with,
control and measure properties of cells.


  In this work, we address a multicoupled dynamics on complex networks with
tunable structural segregation. Specifically, we work on a networked epidemic
spreading under a vaccination campaign with agents in favor and against the
vaccine. Our results show that such coupled dynamics exhibits a myriad of
phenomena such as nonequilibrium transitions accompanied by bistability.
Besides we observe the emergence of an intermediate optimal segregation level
where the community structure enhances negative opinions over vaccination but
counterintuitively hinders - rather than favoring - the global disease
spreading. Thus, our results hint vaccination campaigns should avoid policies
that end up segregating excessively anti-vaccine groups so that they
effectively work as echo chambers in which individuals look to confirmation
without jeopardising the safety of the whole population.


  Seasonal influenza presents an ongoing challenge to public health. The rapid
evolution of the flu virus necessitates annual vaccination campaigns, but the
decision to get vaccinated or not in a given year is largely voluntary, at
least in the United States, and many people decide against it. In early
attempts to model these yearly flu vaccine decisions, it was often assumed that
individuals behave rationally, and do so with perfect information --
assumptions that allowed the techniques of classical economics and game theory
to be applied. However, the usual assumptions are contradicted by the emerging
empirical evidence about human decision-making behavior in this context. We
develop a simple model of coupled disease spread and vaccination dynamics that
instead incorporates experimental observations from social psychology to model
annual vaccine decision-making more realistically. We investigate
population-level effects of these new decision-making assumptions, with the
goal of understanding whether the population can self-organize into a state of
herd immunity, and if so, under what conditions. Our model agrees with
established results while also revealing more subtle population-level behavior,
including biennial oscillations about the herd immunity threshold.


  Manipulation of temperature can be used to actuate DNA origami nano-hinges
containing gold nanoparticles. We develop a physical model of this system that
uses partition function analysis of the interaction between the nano-hinge and
nanoparticle to predict the probability that the nano-hinge is open at a given
temperature. The model agrees well with experimental data and predicts
experimental conditions that allow the actuation temperature of the nano-hinge
to be tuned over a range of temperatures from $30$${}^{\circ}\mathrm{C}$ to
$45$${}^{\circ}\mathrm{C}$. Additionally, the model reveals surprising physical
constraints on the system. This combination of physical insight and predictive
potential is likely to inform future designs that integrate nanoparticles into
dynamic DNA origami structures. Furthermore, our modeling approach could be
expanded to consider the incorporation, stability, and actuation of other types
of functional elements or actuation mechanisms integrated into nucleic acid
devices.


  Living systems operate in a critical dynamical regime -- between order and
chaos -- where they are both resilient to perturbation, and flexible enough to
evolve. To characterize such critical dynamics, the established 'structural
theory' of criticality uses automata network connectivity and node bias (to be
on or off) as tuning parameters. This parsimony in the number of parameters
needed sometimes leads to uncertain predictions about the dynamical regime of
both random and systems biology models of biochemical regulation. We derive a
more accurate theory of criticality by accounting for canalization, the
existence of redundancy that buffers automata response to inputs -- a known
mechanism that buffers the expression of traits, keeping them close to optimal
states despite genetic and environmental perturbations. The new 'canalization
theory' of criticality is based on a measure of effective connectivity. It
contributes to resolving the problem of finding precise ways to design or
control network models of biochemical regulation for desired dynamical
behavior. Our analyses reveal that effective connectivity significantly
improves the prediction of critical behavior in random automata network
ensembles. We also show that the average effective connectivity of a large
battery of systems biology models is much lower than the connectivity of their
original interaction structure. This suggests that canalization has been
selected to dynamically reduce and homogenize the seemingly heterogeneous
connectivity of biochemical networks.


  A common problem in bioinformatics is related to identifying gene regulatory
regions marked by relatively high frequencies of motifs, or deoxyribonucleic
acid sequences that often code for transcription and enhancer proteins.
Predicting alignment scores between subsequence k-mers and a given motif
enables the identification of candidate regulatory regions in a gene, which
correspond to the transcription of these proteins. We propose a one-dimensional
(1-D) Convolution Neural Network trained on k-mer formatted sequences
interspaced with the given motif pattern to predict pairwise alignment scores
between the consensus motif and subsequence k-mers. Our model consists of
fifteen layers with three rounds of a one-dimensional convolution layer, a
batch normalization layer, a dense layer, and a 1-D maximum pooling layer. We
train the model using mean squared error loss on four different data sets each
with a different motif pattern randomly inserted in DNA sequences: the first
three data sets have zero, one, and two mutations applied on each inserted
motif, and the fourth data set represents the inserted motif as a
position-specific probability matrix. We use a novel proposed metric in order
to evaluate the model's performance, $S_{\alpha}$, which is based on the
Jaccard Index. We use 10-fold cross validation to evaluate out model. Using
$S_{\alpha}$, we measure the accuracy of the model by identifying the 15
highest-scoring 15-mer indices of the predicted scores that agree with that of
the actual scores within a selected $\alpha$ region. For the best performing
data set, our results indicate on average 99.3% of the top 15 motifs were
identified correctly within a one base pair stride ($\alpha = 1$) in the out of
sample data. To the best of our knowledge, this is a novel approach that
illustrates how data formatted in an intelligent way can be extrapolated using
machine learning.


  We would like to reply to Neuray et al. who report a series of five new
patients from four unrelated families with bi-allelic mutations of GRN. Their
work nicely completes the few existing reports of similar cases, and refers to
our recent publication describing six homozygous GRN pathogenic variant
carriers with divergent phenotypes and ages at onset (Huin et al., 2020). In
summary, the Letter from Neuray et al., reports valuable findings that lead to
better define CLN11 due to bi-allelic GRN pathogenic variants. Despite the
small sample number that does not allow statistical analysis, the authors
underlined the occurrence of cognitive deterioration and epilepsy. Further
study of the CLN11 families with functional brain imaging and
neuropsychological examinations may be highly informative for the understanding
and the clinical characterization of this rare disease.


  In an earlier paper we proposed a recursive model for epidemics; in the
present paper we generalize this model to include the asymptomatic or
unrecorded symptomatic people, which we call {\em dark people} (dark sector).
We call this the SEPAR$_d$-model. A delay differential equation version of the
model is added; it allows a better comparison to other models. We carry this
out by a comparison with the classical SIR model and indicate why we believe
that the SEPAR$_d$ model may work better for Covid-19 than other approaches. In
the second part of the paper we explain how to deal with the data provided by
the JHU, in particular we explain how to derive central model parameters from
the data. Other parameters, like the size of the dark sector, are less
accessible and have to be estimated more roughly, at best by results of
representative serological studies which are accessible, however, only for a
few countries. We start our country studies with Switzerland where such data
are available. Then we apply the model to a collection of other countries,
three European ones (Germany, France, Sweden), the three most stricken
countries from three other continents (USA, Brazil, India). Finally we show
that even the aggregated world data can be well represented by our approach. At
the end of the paper we discuss the use of the model. Perhaps the most striking
application is that it allows a quantitative analysis of the influence of the
time until people are sent to quarantine or hospital. This suggests that
imposing means to shorten this time is a powerful tool to flatten the curves


  The existence of doublets is a key confounder in single-cell RNA sequencing
(scRNA-seq) data analysis. Computational methods have been developed for
detecting doublets from scRNA-seq data. We developed an R package
DoubletCollection to integrate the installation and execution of eight
doublet-detection methods. DoubletCollection also provides a unified interface
to perform and visualize downstream analysis after doublet detection. Here, we
present a protocol of using DoubletCollection to benchmark doublet-detection
methods. This protocol can automatically accommodate new doublet-detection
methods in the fast-growing scRNA-seq field.


  Chemistry42 is a software platform for de novo small molecule design that
integrates Artificial Intelligence (AI) techniques with computational and
medicinal chemistry methods. Chemistry42 is unique in its ability to generate
novel molecular structures with predefined properties validated through in
vitro and in vivo studies. Chemistry42 is a core component of Insilico Medicine
Pharma.ai drug discovery suite that also includes target discovery and
multi-omics data analysis (PandaOmics) and clinical trial outcomes predictions
(InClinico).


  The Covid-19 pandemic has two key properties: (i) asymptomatic cases (both
detected and undetected) that can result in new infections, and (ii)
time-varying characteristics due to new variants, Non-Pharmaceutical
Interventions etc. We develop a model called SUTRA (Susceptible, Undetected
though infected, Tested positive, and Removed Analysis) that takes into account
both of these two key properties. While applying the model to a region, two
parameters of the model can be learnt from the number of daily new cases found
in the region. Using the learnt values of the parameters the model can predict
the number of daily new cases so long as the learnt parameters do not change
substantially. Whenever any of the two parameters changes due to the key
property (ii) above, the SUTRA model can detect that the values of one or both
of the parameters have changed. Further, the model has the capability to
relearn the changed parameter values, and then use these to carry out the
prediction of the trajectory of the pandemic for the region of concern. The
SUTRA approach can be applied at various levels of granularity, from an entire
country to a district, more specifically, to any large enough region for which
the data of daily new cases are available.
  We have applied the SUTRA model to thirty-two countries, covering more than
half of the world's population. Our conclusions are: (i) The model is able to
capture the past trajectories very well. Moreover, the parameter values, which
we can estimate robustly, help quantify the impact of changes in the pandemic
characteristics. (ii) Unless the pandemic characteristics change significantly,
the model has good predictive capability. (iii) Natural immunity provides
significantly better protection against infection than the currently available
vaccines.


  Recently, [8] has proposed that heterogeneity of infectiousness (and
susceptibility) across individuals in infectious diseases, plays a major role
in affecting the Herd Immunity Threshold (HIT). Such heterogeneity has been
observed in COVID-19 and is recognized as overdispersion (or
"super-spreading"). The model of [8] suggests that super-spreaders contribute
significantly to the effective reproduction factor, R, and that they are likely
to get infected and immune early in the process. Consequently, under R_0 = 3
(attributed to COVID-19), the Herd Immunity Threshold (HIT) is as low as 5%, in
contrast to 67% according to the traditional models [1, 2, 4, 10]. This work
follows up on [8] and proposes that heterogeneity of infectiousness
(susceptibility) has two "faces" whose mix affects dramatically the HIT: (1)
Personal-Trait-, and (2) Event-Based- Infectiousness (Susceptibility). The
former is a personal trait of specific individuals (super-spreaders) and is
nullified once those individuals are immune (as in [8]). The latter is
event-based (e.g cultural super-spreading events) and remains effective
throughout the process, even after the super-spreaders immune. We extend [8]'s
model to account for these two factors, analyze it and conclude that the HIT is
very sensitive to the mix between (1) and (2), and under R_0 = 3 it can vary
between 5% and 67%. Preliminary data from COVID-19 suggests that herd immunity
is not reached at 5%. We address operational aspects and analyze the effects of
lockdown strategies on the spread of a disease. We find that herd immunity (and
HIT) is very sensitive to the lockdown type. While some lockdowns affect
positively the disease blocking and increase herd immunity, others have adverse
effects and reduce the herd immunity.


  Lck (lymphocyte-specific protein tyrosine kinase) is an enzyme which plays a
number of important roles in the function of immune cells. It belongs to the
Src family of kinases which are known to undergo autophosphorylation. It turns
out that this leads to a remarkable variety of dynamical behaviour which can
occur during their activation. We prove that in the presence of
autophosphorylation one phenomenon, bistability, already occurs in a
mathematical model for a protein with a single phosphorylation site. We further
show that a certain model of Lck exhibits oscillations. Finally we discuss the
relations of these results to models in the literature which involve Lck and
describe specific biological processes, such as the early stages of T cell
activation and the stimulation of T cell responses resulting from the
suppression of PD-1 signalling which is important in immune checkpoint therapy
for cancer.


  We investigated progestin and corticosteroid activation of the progesterone
receptor (PR) from elephant shark (Callorhinchus milii), a cartilaginous fish
belonging to the oldest group of jawed vertebrates. Comparison with human PR
experiments provides insights into the evolution of steroid activation of human
PR. At 1 nM steroid, elephant shark PR is activated by progesterone,
17-hydroxy-progesterone, 20beta-hydroxy-progesterone, 11-deoxycorticosterone
(21-hydroxyprogesterone) and 11-deoxycortisol. At 1 nM steroid, human PR is
activated only by progesterone and11-deoxycorticosterone indicating increased
specificity for progestins and corticosteroids during the evolution of human
PR. RU486, an important clinical antagonist of human PR, did not inhibit
progesterone activation of elephant shark PR. Cys-528 in elephant shark PR
corresponds to Gly-722 in human PR, which is essential for RU486 inhibition of
human PR. Confirming the importance of this site on elephant shark PR, RU486
inhibited progesterone activation of the Cys528Gly mutant PR. There also was a
decline in activation of elephant shark Cys528Gly PR by 11-deoxycortisol,
17-hydroxy-progesterone and 20beta-hydroxy-progesterone and an increase in
activation of human Gly722Cys PR by 11-deoxycortisol and decreased activation
by corticosterone. One or more of these changes may have selected for the
mutation corresponding to human glycine-722 PR that first evolved in platypus
PR, a basal mammal.


  Recently, the selection-recombination equation with a single selected site
and an arbitrary number of neutral sites was solved by means of the ancestral
selection-recombination graph. Here, we introduce a more accessible approach,
namely the ancestral initiation graph. The construction is based on a
discretisation of the selection-recombination equation. We apply our method to
systematically explain a long-standing observation concerning the dynamics of
linkage disequilibrium between two neutral loci hitchhiking along with a
selected one. In particular, this clarifies the nontrivial dependence on the
position of the selected site.


  Structural flexibility and/or dynamic interactions with other molecules is a
critical aspect of protein function. CryoEM provides direct visualization of
individual macromolecules sampling different conformational and compositional
states. While numerous methods are available for computational classification
of discrete states, characterization of continuous conformational changes or
large numbers of discrete state without human supervision remains challenging.
Here we present e2gmm, a machine learning algorithm to determine a
conformational landscape for proteins or complexes using a 3-D Gaussian mixture
model mapped onto 2-D particle images in known orientations. Using a deep
neural network architecture, e2gmm can automatically resolve the structural
heterogeneity within the protein complex and map particles onto a small latent
space describing conformational and compositional changes. This system presents
a more intuitive and flexible representation than other manifold methods
currently in use. We demonstrate this method on both simulated data as well as
three biological systems, to explore compositional and conformational changes
at a range of scales. The software is distributed as part of EMAN2.


  Automated classification of animal vocalisations is a potentially powerful
wildlife monitoring tool. Training robust classifiers requires sizable
annotated datasets, which are not easily recorded in the wild. To circumvent
this problem, we recorded four primate species under semi-natural conditions in
a wildlife sanctuary in Cameroon with the objective to train a classifier
capable of detecting species in the wild. Here, we introduce the collected
dataset, describe our approach and initial results of classifier development.
To increase the efficiency of the annotation process, we condensed the
recordings with an energy/change based automatic vocalisation detection.
Segmenting the annotated chunks into training, validation and test sets,
initial results reveal up to 82% unweighted average recall (UAR) test set
performance in four-class primate species classification.


  Predation driven Allee effects play an important role in the dynamics of
small population, however, such predation-driven Allee effects cannot occur for
the model with type I functional response. It generally occurs when a
generalist predator targets some specific prey. However, apart from the lethal
effects of predation, there are some non-lethal effects in the presence of
predator. Due to the fear of predation, positive density dependence growth may
be observed at low population density, because of reduced foraging activities.
Moreover, this non-lethal effect can be carried over generations. In the
present manuscript, we investigate the role of predation fear and its
carry-over effects in prey-predator model. First, we study the single species
model in global perspective. We have shown that depending on the birth rate,
our single species model describes three types of growth dynamics, namely,
strong Allee dynamics, weak Allee dynamics and logistic dynamics. Then we
consider the explicit dynamics of predator, with type I functional response.
Basic dynamical properties, as well as global stability of each equilibria have
been discussed. From our analysis, we can observe that both the fear and its
carry-over effects have significant role in the stability of the coexistence
equilibrium, even if for the model with type I functional response. The
phenomenon paradox of enrichment can be observed in our model, which cannot be
observed in the classical prey-predator model with type I functional response.
However, we can see that such phenomenon can be ruled out by choosing suitable
non-lethal effect parameters. Therefore, our study shows how non-lethal effects
change the dynamics of a competition model, and has important biological
insights, specially for the understanding of the dynamics of small populations.


  Motivation: Hi-C matrices are cornerstones for qualitative and quantitative
studies of genome folding, from its territorial organization to compartments
and topological domains. The high dynamic range of genomic distances probed in
Hi-C assays reflects in an inherent stochastic background of the interactions
matrices, which inevitably convolve the features of interest with largely
aspecific ones. Results: Here we introduce a discuss essHi-C, a method to
isolate the specific, or essential component of Hi-C matrices from the
aspecific portion of the spectrum that is compatible with random matrices.
Systematic comparisons show that essHi-C improves the clarity of the
interaction patterns, enhances the robustness against sequencing depth, allows
the unsupervised clustering of experiments in different cell lines and recovers
the cell-cycle phasing of single-cells based on Hi-C data. Thus, essHi-C
provides means for isolating significant biological and physical features from
Hi-C matrices.


  Polyploidization, whereby an organism inherits multiple copies of the genome
of their parents, is an important evolutionary event that has been observed in
plants and animals. One way to study such events is in terms of the ploidy
number of the species that make up a dataset of interest. It is therefore
natural to ask: How much information about the evolutionary past of the set of
species that form a dataset can be gleaned from the ploidy numbers of the
species? To help answer this question, we introduce and study the novel concept
of a ploidy profile which allows us to formalize it in terms of a multiplicity
vector indexed by the species the dataset is comprised of. Using the framework
of a phylogenetic network, we present a closed formula for computing the hybrid
number (i.e. the minimal number of polyploidization events required to explain
a ploidy profile) of a large class of ploidy profiles. This formula relies on
the construction of a certain phylogenetic network from the simplification
sequence of a ploidy profile and the hybrid number of the ploidy profile with
which this construction is initialized. Both of them can be computed easily in
case the ploidy numbers that make up the ploidy profile are not too large. To
help illustrate the applicability of our approach, we apply it to a simplified
version of a publicly available Viola dataset.


  This paper is concerned with a reaction--diffusion system modeling the
fixation and the invasion in a population of a gene drive (an allele biasing
inheritance, increasing its own transmission to offspring). In our model, the
gene drive has a negative effect on the fitness of individuals carrying it, and
is therefore susceptible of decreasing the total carrying capacity of the
population locally in space. This tends to generate an opposing demographic
advection that the gene drive has to overcome in order to invade. While
previous reaction--diffusion models neglected this aspect, here we focus on it
and try to predict the sign of the traveling wave speed. It turns out to be an
analytical challenge, only partial results being within reach, and we complete
our theoretical analysis by numerical simulations. Our results indicate that
taking into account the interplay between population dynamics and population
genetics might actually be crucial, as it can effectively reverse the direction
of the invasion and lead to failure. Our findings can be extended to other
bistable systems, such as the spread of cytoplasmic incompatibilities caused by
Wolbachia.


  Cancer cells obtain mutations which rely on the production of diffusible
growth factors to confer a fitness benefit. These mutations can be considered
cooperative, and studied as public goods games within the framework of
evolutionary game theory. The population structure, benefit function and update
rule all influence the evolutionary success of cooperators. We model the
evolution of cooperation in epithelial cells using the Voronoi tessellation
model. Unlike traditional evolutionary graph theory, this allows us to
implement global updating, for which birth and death events are spatially
decoupled. We compare, for a sigmoid benefit function, the conditions for
cooperation to be favoured and/or beneficial for well mixed and structured
populations. We find that when population structure is combined with global
updating, cooperation is more successful than if there were local updating or
the population were well-mixed. Interestingly, the qualitative behaviour for
the well-mixed population and the Voronoi tessellation model is remarkably
similar, but the latter case requires significantly lower incentives to ensure
cooperation.


  In Part 1, we introduced a stochastic model of an infectious disease, based
on the BDI (birth and death with immigration) process. We showed that random
processes defined by this model can capture the essence of the stochastic,
often erratic, behavior of the infection process. The most significant finding
was that it is negative binomial distributed with small r, hence it is
geometrically distributed with an exceedingly long tail. This leads to a much
larger disparity in the epidemic patterns than has been known to the modeling
community. In Part 2 we conduct simulation experiments by implementing an
event-driven simulator. Several independent runs are presented to verify the
findings reported in Part 1. The enormous variations among the sample paths
obtained from several consecutive and independent runs with statistically
identical conditions confirm our analysis. Note that the epidemic pattern that
we observe is merely one sample path taken out of infinitely many paths. Thus,
one such path cannot represent the ensemble of the underlying random process.
By plotting the simulation runs in the semi-log scale, we can see that the
probabilistic chances in the early phase of the infection process determine the
behavior of the process. Once the process has reached a considerable number,
the weak law of large numbers sets in, and the process behaves less erratically
than in the early phase. One important implication of our findings is that it
would be a futile effort to attempt to identify all plausible causes of the
epidemic patterns. Mere luck may play a more significant role than most people
believe. It would be worth remarking that the BDI process might be applicable
to explain the disparity in wealth among individuals of similar earning power,
expenditure, and investment portfolio.


  Infectious diseases are a threat for human health with tremendous impact on
our society at large. The recent COVID-19 pandemic, caused by the SARS-CoV-2,
is the latest example of a highly infectious disease ravaging the world, since
late 2019. It is therefore imperative to develop efficient mathematical models,
able to substantially curb the damages of a pandemic by unveiling disease
spreading dynamics and symmetries. This will help inform (non)-pharmaceutical
prevention strategies. For the reasons above we wrote this report that goes at
the heart of mathematical modelling of infectious disease diffusion by
simultaneously investigating the underlying microscopic dynamics in terms of
percolation models, effective description via compartmental models and the
employment of temporal symmetries naturally encoded in the mathematical
language of critical phenomena. Our report reviews these approaches and
determines their common denominators, relevant for theoretical epidemiology and
its link to important concepts in theoretical physics. We show that the
different frameworks exhibit common features such as criticality and
self-similarity under time rescaling. These features are naturally encoded
within the unifying field theoretical approach. The latter leads to an
efficient description of the time evolution of the disease via a framework in
which (near) time-dilation invariance is explicitly realised. As important test
of the relevance of symmetries we show how to mathematically account for
observed phenomena such as multi-wave dynamics. The models presented here are
of immediate relevance for different realms of scientific enquiry from medical
applications to the understanding of human behaviour. Our review offers novel
perspectives on how to model, capture, organise and understand epidemiological
data and disease dynamics for modelling real-world phenomena.


  Now that SARS-CoV-2 (COVID-19) vaccines are developed, it is very important
to plan its distribution strategy. In this paper, we formulated a
multi-objective linear programming model to optimize vaccine distribution and
applied it to the agent-based version of our age-stratified and
quarantine-modified SEIR with non-linear incidence rates (ASQ-SEIR-NLIR)
compartmental model. Simulations were performed using COVID-19 data from Quezon
City and results were analyzed under various scenarios: (1) no vaccination, (2)
base vaccination (prioritizing essential workers and vulnerable population),
(3) prioritizing mobile workforce, (4) prioritizing elderly, and (5)
prioritizing mobile workforce and elderly; in terms of (a) reducing infection
rates and (b) reducing mortality incidence. After 10 simulations on
distributing 500,000 vaccine courses, results show that prioritizing mobile
workforce minimizes further infections by 24.14%, which is better than other
scenarios. On the other hand, prioritizing the elderly yields the highest
protection (439%) for the Quezon City population compared to other scenarios.
This could be due to younger people, when contracted the disease, has higher
chances of recovery than the elderly. Thus, this leads to reduction of
mortality cases.


  Parrondo's paradox indicates a paradoxical situation in which a winning
expectation may occur in sequences of losing games. There are many versions of
the original Parrondo's games in the literature, but the games are played by
two players in all of them. We introduce a new extended version of games played
by three players and a three-sided biased dice instead of two players and a
biased coin in this work. In the first step, we find the part of the parameters
space where the games are played fairly. After adding noise to fair
probabilities, we combine two games randomly, periodically, and nonlinearly and
obtain the conditions under which the paradox can occur. This generalized model
can be applied in all science and engineering fields. It can also be used for
genetic switches. Genetic switches are often made by two reactive elements, but
the existence of more elements can lead to more existing decisions for cells.
Each genetic switch can be considered a game in which the reactive elements
compete to increase their molecular concentrations. We present three genetic
networks based on a new generalized Parrondo's games model, consisting of two
noisy genetic switches. The combination of them can increase network robustness
to noise. Each switch can also be used as an initial pattern to construct a
synthetic switch to change undesirable cells' fate.


  The COVID-19 pandemic has had worldwide devastating effects on human lives,
highlighting the need for tools to predict its development. Dynamics of such
public-health threats can often be efficiently analysed through simple models
that help to make quantitative timely policy decisions. We benchmark a minimal
version of a Susceptible-Infected-Removed model for infectious diseases (SIR)
coupled with a simple least-squares Statistical Heuristic Regression (SHR)
based on a lognormal distribution. We derived the three free parameters for
both models in several cases and tested them against the amount of data needed
to bring accuracy in predictions. The SHR model is approximately +/- 2%
accurate about 20 days past the second inflexion point in the daily curve of
cases, while the SIR model reaches a similar accuracy a fortnight before. All
the analyzed cases assert the utility of SHR and SIR approximants as a useful
tool to forecast the evolution of the disease. Finally, we have studied
simulated stochastic individual-based SIR dynamics, which yields a detailed
spatial and temporal view of the disease that cannot be given by SIR or SHR
methods.


  In 2018, West Nile Virus (WNV) was detected for the first time in Germany.
Since the first detection, 36 human cases and 175 cases in horses and birds are
detected. The transmission cycle of West Nile Virus includes birds and
mosquitoes and -- as dead-end hosts -- people and horses. Spatial dissemination
of the disease is caused by the movements of birds and mosquitoes. While the
activity and movement of mosquitoes are depending mainly on temperature, in the
birds there is a complex movement pattern caused by local birds and long range
dispersal birds. To this end, we have developed a metapopulation network model
framework to delineate the potential spatial distribution and spread of WNV
across Germany as well as to evaluate the risk throughout our proposed network
model. Our model facilitates the interconnection amongst the vector, local
birds and long range dispersal birds contact networks. We have assumed
different distance dispersal kernels models for the vector and avian
populations with the intention to include short and long range dispersal. The
model includes spatial variation of mosquito abundance and the movements to
resemble the reality.


  Crimean-Congo haemorrhagic fever (CCHF) is a tick-borne zoonotic disease
caused by the Crimean-Congo hemorrhagic fever virus (CCHFV). Ticks belonging to
the genus \textit{Hyalomma} are the main vectors and reservoir for the virus.
It is maintained in nature in an endemic vertebrate-tick-vertebrate cycle.
CCHFV is prevalent in wide geographical areas including Asia, Africa,
South-Eastern Europe and the Middle East. Over the last decade, several
outbreaks of CCHFV have been observed in Europe, mainly in Mediterranean
countries. Due to the high case/fatality ratio of CCHFV in human sometimes, it
is of great importance for public health. Climate change and the invasion of
CCHFV vectors in Central Europe suggest that the establishment of the
transmission in Central Europe may be possible in future. We developed a
compartment-based nonlinear Ordinary Differential Equation (ODE) system to
model the disease transmission cycle including blood sucking ticks, livestock
and human. Sensitivity analysis of the basic reproduction number $R_0$ shows
that decreasing in the tick survival time is an efficient method to eradicate
the disease. The model supports us in understanding the influence of different
model parameters on the spread of CCHFV. Tick to tick transmission through
co-feeding and the CCHFV circulation through trasstadial and transovarial
stages are important factors to sustain the disease cycle. The proposed model
dynamics are calibrated through an empirical multi-country analysis and
multidimensional scaling reveals the disease-parameter sets of different
countries burdened with CCHF are different. This necessary information may help
us to select most efficient control strategies.


  Quantifying and comparing patterns of dynamical ecological systems require
averaging over measurable quantities. For example, to infer variation in
movement and behavior, metrics such as step length and velocity are averaged
over large ensembles. Yet, in nonergodic systems such averaging is
inconsistent; thus, identifying ergodicity breaking is essential in ecology.
Using rich high-resolution movement datasets ($>\! 7 \times 10^7$
localizations) from 70 individuals and continuous-time random walk modeling, we
find subdiffusive behavior and ergodicity breaking in the localized movement of
three species of avian predators. Small-scale, within-patch movement was found
to be qualitatively different, not inferrable and separated from large-scale
inter-patch movement. Local search is characterized by long
power-law-distributed waiting times with diverging mean, giving rise to
ergodicity breaking in the form of considerable variability uniquely observed
at this scale. This implies that wild animal movement is scale specific with no
typical waiting time at the local scale.


  A two-type two-sex branching process is introduced with the aim of describing
the interaction of predator and prey populations with sexual reproduction and
promiscuous mating. In each generation and in each species the total number of
individuals which mate and produce offspring is controlled by a binomial
distribution with size given by this number of individuals and probability of
success depending on the density of preys per predator. The resulting model
enables us to depict the typical cyclic behaviour of predator-prey systems
under some mild assumptions on the shape of the function that characterises the
probability of survival of the previous binomial distribution. We present some
basic results about fixation and extinction of both species as well as
conditions for the coexistence of both of them. We also analyse the suitability
of the process to model real ecosystems comparing our model with a real
dataset.


  We seek to understand the interplay between amino acid sequence and local
structure in proteins. Are some amino acids unique in their ability to fit
harmoniously into certain local structures? What is the role of sequence in
sculpting the putative native state folds from myriad possible conformations?
In order to address these questions, we represent the local structure of each
C-alpha atom of a protein by just two angles, theta and mu, and we analyze a
set of more than 4000 protein structures from the PDB. We use a hierarchical
clustering scheme to divide the 20 amino acids into six distinct groups based
on their similarity to each other in fitting local structural space. We present
the results of a detailed analysis of patterns of amino acid specificity in
adopting local structural conformations and show that the sequence-structure
correlation is not very strong compared to a random assignment of sequence to
structure. Yet, our analysis may be useful to determine an effective scoring
rubric for quantifying the match of an amino acid to its putative local
structure.


  Network epidemiology has become a vital tool in understanding the effects of
high-degree vertices, geographic and demographic communities, and other
inhomogeneities in social structure on the spread of disease. However, many
networks derived from modern datasets are quite dense, such as mobility
networks where each location has links to a large number of potential
destinations. One way to reduce the computational effort of simulating
epidemics on these networks is sparsification, where we select a representative
subset of edges based on some measure of their importance. Recently an approach
was proposed using an algorithm based on the effective resistance of the edges.
We explore how effective resistance is correlated with the probability that an
edge transmits disease in the SI model. We find that in some cases these two
notions of edge importance are well correlated, making effective resistance a
computationally efficient proxy for the importance of an edge to epidemic
spread. In other cases, the correlation is weaker, and we discuss situations in
which effective resistance is not a good proxy for epidemic importance.


  Aerosolized droplets play a central role in the transmission of various
infectious diseases, including Legionnaire's disease, gastroenteritis-causing
norovirus, and most recently COVID-19. Respiratory droplets are known to be the
most prominent source of transmission for COVID-19, however, alternative routes
may exist given the discovery of small numbers of viable viruses in urine and
stool samples. Flushing biomatter can lead to the aerosolization of
microorganisms, thus, there is a likelihood that bioaerosols generated in
public restrooms may pose a concern for the transmission of COVID-19,
especially since these areas are relatively confined, experience heavy foot
traffic, and may suffer from inadequate ventilation. To quantify the extent of
aerosolization, we measure the size and number of droplets generated by
flushing toilets and urinals in a public restroom. The results indicate that
the particular designs tested in the study generate a large number of droplets
in the size range 0.3$\mu m$ to 3$\mu m$, which can reach heights of at least
1.52m. Covering the toilet reduced aerosol levels but did not eliminate them
completely, suggesting that aerosolized droplets escaped through small gaps
between the cover and the seat. In addition to consistent increases in aerosol
levels immediately after flushing, there was a notable rise in ambient aerosol
levels due to the accumulation of droplets from multiple flushes conducted
during the tests. This highlights the need for incorporating adequate
ventilation in the design and operation of public spaces, which can help
prevent aerosol accumulation in high occupancy areas and mitigate the risk of
airborne disease transmission.


  Diffusion comes anytime and everywhere. If there is a gradient or a potential
difference of a quantity a diffusion process happens and this ends if an
equilibrium is reached only. The concentration of a species maybe such
quantity, or the voltage. An electric currant will be driven by a voltage
difference for example. In this COVID-19 pandemic one observes both regions
with low incidence and other ones with high incidence. The local different
people density could be a reason for that. In populous areas like big cities or
congested urban areas higher COVID-19 incidences could be observed than in
rural regions. The aim of this paper consists in the application of a diffusion
concept to describe one possible issue of the the COVID-19 propagation. This
will be discussed for the German situation based on the quite different
incidence data for the different federal states of Germany.


  On the basis of a compartment model, the epidemic curve is investigated when
the net rate $\lambda$ of change of the number of infected individuals $I$ is
given by an ellipse in the $\lambda$-$I$ plane which is supported in
$[I_{\ell}, I_h]$. With $a \equiv (I_h - I_{\ell})/(I_h + I_{\ell})$, it is
shown that (1) when $a < 1$ or $I_{\ell} >0$, oscillation of the infection
curve is self-organized and the period of the oscillation is in proportion to
the ratio of the difference $ (I_h - I_{\ell})$ and the geometric mean
$\sqrt{I_h I_{\ell}}$ of $I_h$ and $I_{\ell}$, (2) when $a = 1$, the infection
curve shows a critical behavior where it decays obeying a power law function
with exponent $-2$ in the long time limit after a peak, and (3) when $a > 1$,
the infection curve decays exponentially in the long time limit after a peak.
The present result indicates that the pandemic can be controlled by a measure
which makes $I_{\ell} < 0$.


  Agent-based models of disease transmission involve stochastic rules that
specify how a number of individuals would infect one another, recover or be
removed from the population. Common yet stringent assumptions stipulate
interchangeability of agents and that all pairwise contact are equally likely.
Under these assumptions, the population can be summarized by counting the
number of susceptible and infected individuals, which greatly facilitates
statistical inference. We consider the task of inference without such
simplifying assumptions, in which case, the population cannot be summarized by
low-dimensional counts. We design improved particle filters, where each
particle corresponds to a specific configuration of the population of agents,
that take either the next or all future observations into account when
proposing population configurations. Using simulated data sets, we illustrate
that orders of magnitude improvements are possible over bootstrap particle
filters. We also provide theoretical support for the approximations employed to
make the algorithms practical.


  Motivation: Peptides have attracted the attention in this century due to
their remarkable therapeutic properties. Computational tools are being
developed to take advantage of existing information, encapsulating knowledge
and making it available in a simple way for general public use. However, these
are property-specific redundant data systems, and usually do not display the
data in a clear way. In some cases, information download is not even possible.
This data needs to be available in a simple form for drug design and other
biotechnological applications.
  Results: We developed Peptipedia, a user-friendly database and web
application to search, characterise and analyse peptide sequences. Our tool
integrates the information from thirty previously reported databases, making it
the largest repository of peptides with recorded activities so far. Besides, we
implemented a variety of services to increase our tool's usability. The
significant differences of our tools with other existing alternatives becomes a
substantial contribution to develop biotechnological and bioengineering
applications for peptides.
  Availability: Peptipedia is available for non-commercial use as an
open-access software, licensed under the GNU General Public License, version
GPL 3.0. The web platform is publicly available at pesb2.cl/peptipedia. Both
the source code and sample datasets are available in the GitHub repository
https://github.com/CristoferQ/PeptideDatabase.
  Contact: david.medina@cebib.cl, ana.sanchez@ing.uchile.cl


  It is well recognized that population heterogeneity plays an important role
in the spread of epidemics. While individual variations in social activity are
often assumed to be persistent, i.e. constant in time, here we discuss the
consequences of dynamic heterogeneity. By integrating the stochastic dynamics
of social activity into traditional epidemiological models we demonstrate the
emergence of a new long timescale governing the epidemic in broad agreement
with empirical data. Our model captures multiple features of real-life
epidemics such as COVID-19, including prolonged plateaus and multiple waves,
which are transiently suppressed due to the dynamic nature of social activity.
The existence of the long timescale due to the interplay between epidemic and
social dynamics provides a unifying picture of how a fast-paced epidemic
typically will transition to the endemic state.


  The complexities involved in modelling the transmission dynamics of COVID-19
has been a roadblock in achieving predictability in the spread and containment
of the disease. In addition to understanding the modes of transmission, the
effectiveness of the mitigation methods also needs to be built into any
effective model for making such predictions. We show that such complexities can
be circumvented by appealing to scaling principles which lead to the emergence
of universality in the transmission dynamics of the disease. The ensuing data
collapse renders the transmission dynamics largely independent of geopolitical
variations, the effectiveness of various mitigation strategies, population
demographics, etc. We propose a simple two-parameter model -- the Blue Sky
model -- and show that one class of transmission dynamics can be explained by a
solution that lives at the edge of a blue sky bifurcation. In addition, the
data collapse leads to an enhanced degree of predictability in the disease
spread for several geographical scales which can also be realized in a
model-independent manner as we show using a deep neural network. The
methodology adopted in this work can potentially be applied to the transmission
of other infectious diseases and new universality classes may be found. The
predictability in transmission dynamics and the simplicity of our methodology
can help in building policies for exit strategies and mitigation methods during
a pandemic.


  It is conceivable that an RNA virus could use a polysome, that is, a string
of ribosomes covering the RNA strand, to protect the genetic material from
degradation inside a host cell. This paper discusses how such a virus might
operate, and how its presence might be detected by ribosome profiling. There
are two possible forms for such a polysomally protected virus, depending upon
whether just the forward strand or both the forward and complementary strands
can be encased by ribosomes (these will be termed type 1 and type 2,
respectively). It is argued that in the type 2 case the viral RNA would evolve
an ambigrammatic property, whereby the viral genes are free of stop codons in a
reverse reading frame (with forward and reverse codons aligned). Recent
observations of ribosome profiles of ambigrammatic narnavirus sequences are
consistent with our predictions for the type 2 case.


  The success of modern civilization is built upon widespread cooperation in
human society, deciphering the mechanisms behind has being a major goal for
centuries. A crucial fact is, however, largely missing in most prior studies
that games in the real world are typically played simultaneously and
interactively rather than separately as assumed. Here we introduce the idea of
interacting games that different games coevolve and influence each other's
decision-making. We show that as the game-game interaction becomes important,
the cooperation phase transition dramatically improves, a fairly high level of
cooperation is reached for all involved games when interaction goes to be
strong. A mean-field theory indicates that a new mechanism -- \emph{the
dynamical reciprocity}, as a counterpart to the well-known network reciprocity,
is at work to foster cooperation, which is confirmed by the detailed analysis.
This revealed reciprocity is robust against variations in the game type, the
population structure, and the updating rules etc, and more games generally
yield a higher level of cooperation. Our findings point out the great potential
towards high cooperation for many issues are interwoven with each other in the
real world, and also the possibility of sustaining decent cooperation even in
extremely adverse circumstances.


  We study the evolution of two mutually interacting games with both pairwise
games as well as the public goods game on different topologies. On 2d square
lattices, we reveal that the game-game interaction can promote the cooperation
prevalence in all cases, and the cooperation-defection phase transitions even
become absent and fairly high cooperation is expected when the interaction goes
to be very strong. A mean-field theory is developed that points out new
dynamical routes arising therein. Detailed analysis shows indeed that there are
rich categories of interactions in either individual or bulk scenario:
invasion, neutral, and catalyzed types; their combination puts cooperators at a
persistent advantage position, which boosts the cooperation. The robustness of
the revealed reciprocity is strengthened by the studies of model variants,
including asymmetrical or time-varying interactions, games of different types,
games with time-scale separation, different updating rules etc. The structural
complexities of the underlying population, such as Newman--Watts small world
networks, Erd\H{o}s--R\'enyi random networks, and Barab\'asi--Albert networks,
also do not alter the working of the dynamical reciprocity. In particular, as
the number of games engaged increases, the cooperation level continuously
improves in general. Our work thus uncovers a new class of cooperation
mechanism and indicates the great potential for human cooperation where
concurrent issues are so often seen in the real world.


  There is a long-time quest for understanding physical mechanisms of weak
magnetic field interaction with biological matter. Two factors impeded the
development of such mechanisms: first, a high (room) temperature of a cellular
environment, where a weak, static magnetic field induces a (classically) zero
equilibrium response. Second, the friction in the cellular environment is
large, preventing a weak field to alter non-equilibrium processes such as a
free diffusion of charges. Here we study a class of non-equilibrium steady
states of a cellular ion in a confining potential, where the response to a
(weak, homogeneous, static) magnetic field survives strong friction and thermal
fluctuations. The magnetic field induces a rotational motion of the ion that
proceeds with the cyclotron frequency. Such non-equilibrium states are
generated by a white noise acting on the ion additionally to the non-local
(memory-containing) friction and noise generated by an equilibrium thermal
bath. The intensity of this white noise can be weak, i.e. much smaller than the
thermal noise intensity.


  Accurate and robust prediction of patient's response to drug treatments is
critical for developing precision medicine. However, it is often difficult to
obtain a sufficient amount of coherent drug response data from patients
directly for training a generalized machine learning model. Although the
utilization of rich cell line data provides an alternative solution, it is
challenging to transfer the knowledge obtained from cell lines to patients due
to various confounding factors. Few existing transfer learning methods can
reliably disentangle common intrinsic biological signals from confounding
factors in the cell line and patient data. In this paper, we develop a Coherent
Deconfounding Autoencoder (CODE-AE) that can extract both common biological
signals shared by incoherent samples and private representations unique to each
data set, transfer knowledge learned from cell line data to tissue data, and
separate confounding factors from them. Extensive studies on multiple data sets
demonstrate that CODE-AE significantly improves the accuracy and robustness
over state-of-the-art methods in both predicting patient drug response and
de-confounding biological signals. Thus, CODE-AE provides a useful framework to
take advantage of in vitro omics data for developing generalized patient
predictive models. The source code is available at
https://github.com/XieResearchGroup/CODE-AE.


  Glucose biosensors play an important role in the diagnosis and continued
monitoring of the disease, diabetes mellitus. This report proposes the
development of a novel enzymatic electrochemical glucose biosensor based on
TiO$_2$ nanotubes modified by AgO and Prussian blue (PB) nanoparticles (NPs),
which has an additional advantage of possessing antimicrobial properties for
implantable biosensor applications. In this study, we developed two high
performance glucose biosensors based on the immobilization of glucose oxidase
(GOx) onto Prussian blue (PB) modified TiO$_2$ nanotube arrays functionalized
by Au and AgO NPs. AgO-deposited TiO$_2$ nanotubes were synthesized through an
electrochemical anodization process followed by Ag electroplating process in
the same electrolyte. Deposition of PB particles was performed from an acidic
ferricyanide solution. The surface morphology and elemental composition of the
two fabricated biosensors were investigated by scanning electron microscopy
(SEM) and energy-dispersive X-ray spectroscopy (EDS) which indicate the
successful deposition of Au and AgO nanoparticles as well as PB nanocrystals.
Cyclic voltammetry and chronoamperometry were used to investigate the
performance of the modified electrochemical biosensors. The results show that
the developed electrochemical biosensors display excellent properties in terms
of electron transmission, low detection limit as well as high stability for the
determination of glucose. Under the optimized conditions, the amperometric
response shows a linear dependence on the glucose concentration to a detection
limit down to 4.91 $\mu$M with sensitivity of 185.1 mA M$^{-1}$ cm$^{-2]$ in Au
modified biosensor and detection limit of 58.7 $\mu$M with 29.1 mA M$^{-1}$
cm$^{-2]$ sensitivity in AgO modified biosensor.


  For classic Lotka-Volterra systems governing many interacting species, we
establish an exclusion principle that rules out the existence of linearly
asymptotically stable steady states in subcommunitites of communities that
admit a stable state which is internally D-stable. This type of stability is
known to be ensured, e.g., by diagonal dominance or Volterra-Lyapunov stability
conditions. By consequence, the number of stable steady states of this type is
bounded by Sperner's lemma on anti-chains in a poset. The number of stable
steady states can nevertheless be very large if there are many groups of
species that strongly inhibit outsiders but have weak interactions among
themselves. By examples we also show that in general it is possible for a
stable community to contain a stable subcommunity consisting of a single
species. Thus a recent empirical finding to the contrary, in a study of random
competitive systems by Lischke and L\"offler (Theo. Pop. Biol. 115 (2017)
24--34), does not hold without qualification.


  Our work is concerned with the generation and targeted design of RNA, a type
of genetic macromolecule that can adopt complex structures which influence
their cellular activities and functions. The design of large scale and complex
biological structures spurs dedicated graph-based deep generative modeling
techniques, which represents a key but underappreciated aspect of computational
drug discovery. In this work, we investigate the principles behind representing
and generating different RNA structural modalities, and propose a flexible
framework to jointly embed and generate these molecular structures along with
their sequence in a meaningful latent space. Equipped with a deep understanding
of RNA molecular structures, our most sophisticated encoding and decoding
methods operate on the molecular graph as well as the junction tree hierarchy,
integrating strong inductive bias about RNA structural regularity and folding
mechanism such that high structural validity, stability and diversity of
generated RNAs are achieved. Also, we seek to adequately organize the latent
space of RNA molecular embeddings with regard to the interaction with proteins,
and targeted optimization is used to navigate in this latent space to search
for desired novel RNA molecules.


  The deadly coronavirus disease 2019 (COVID-19) pandemic caused by severe
acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has gone out of control
globally. Despite much effort by scientists, medical experts, and society in
general, the slow progress on drug discovery and antibody therapeutic
development, the unknown possible side effects of the existing vaccines, and
the high transmission rate of the SARS-CoV-2, remind us of the sad reality that
our current understanding of the transmission, infectivity, and evolution of
SARS-CoV-2 is unfortunately very limited. The major limitation is the lack of
mechanistic understanding of viral-host cell interactions, the viral
regulation, protein-protein interactions, including antibody-antigen binding,
protein-drug binding, host immune response, etc. This limitation will likely
haunt the scientific community for a long time and have a devastating
consequence in combating COVID-19 and other pathogens. Notably, compared to the
long-cycle, highly cost, and safety-demanding molecular-level experiments, the
theoretical and computational studies are economical, speedy, and easy to
perform. There exists a tsunami of the literature on molecular modeling,
simulation, and prediction of SARS-CoV-2 that has become impossible to fully be
covered in a review. To provide the reader a quick update about the status of
molecular modeling, simulation, and prediction of SARS-CoV-2, we present a
comprehensive and systematic methodology-centered narrative in the nick of
time. Aspects such as molecular modeling, Monte Carlo (MC) methods, structural
bioinformatics, machine learning, deep learning, and mathematical approaches
are included in this review. This review will be beneficial to researchers who
are looking for ways to contribute to SARS-CoV-2 studies and those who are
assessing the current status in the field.


  The paper was suggested by a brief note of the second author about the
application of the Hubbert curve to predict decay of resource exploitation. A
further suggestion came from the interpretation of the Hubbert curve in terms
of a specific Lotka Volterra (LV) equation. The link with population dynamics
was obvious as logistic function and LV equation were proposed within the
demography science field. Mathematical population dynamics has a history of
about two centuries. The first principle and model of population dynamics can
be regarded the exponential law of Malthus. In the XIX century, the Malthusian
demographic model was first refined to include mortality rate by Gompertz. In
the early XIX century the model was further refined by Verhulst by introducing
the standard logistic function. The previous models only concern the population
of a single species. In the early XX century, the American demographer Lotka
and the Italian mathematician Volterra proposed a pair of state equations which
describe the population dynamics of two competing species, the predator and the
prey. This paper is concerned with the single and two-species fundamental
equations: the logistic and LV equation. The paper starts with the generalized
logistic equation whose free response is derived together with equilibrium
points and stability properties. The parameter estimation of the logistic
function is applied to the raw data of the US crude oil production. The paper
proceeds with the Lotka Volterra equation of the competition between two
species, with the goal of applying it to resource exploitation. At the end, a
limiting version of the LV equation is studied since it describes a competition
model between the production rate of exploited resources and the relevant
capital stock employed in the exploitation.


  This paper introduces the Gene Mover's Distance, a measure of similarity
between a pair of cells based on their gene expression profiles obtained via
single-cell RNA sequencing. The underlying idea of the proposed distance is to
interpret the gene expression array of a single cell as a discrete probability
measure. The distance between two cells is hence computed by solving an Optimal
Transport problem between the two corresponding discrete measures. In the
Optimal Transport model, we use two types of cost function for measuring the
distance between a pair of genes. The first cost function exploits a gene
embedding, called gene2vec, which is used to map each gene to a high
dimensional vector: the cost of moving a unit of mass of gene expression from a
gene to another is set to the Euclidean distance between the corresponding
embedded vectors. The second cost function is based on a Pearson distance among
pairs of genes. In both cost functions, the more two genes are correlated, the
lower is their distance. We exploit the Gene Mover's Distance to solve two
classification problems: the classification of cells according to their
condition and according to their type. To assess the impact of our new metric,
we compare the performances of a $k$-Nearest Neighbor classifier using
different distances. The computational results show that the Gene Mover's
Distance is competitive with the state-of-the-art distances used in the
literature.


  This paper introduces new methods to study the changing dynamics of COVID-19
cases and deaths among the 50 worst-affected countries throughout 2020. First,
we analyse the trajectories and turning points of rolling mortality rates to
understand at which times the disease was most lethal. We demonstrate five
characteristic classes of mortality rate trajectories and determine structural
similarity in mortality trends over time. Next, we introduce a class of
\emph{virulence matrices} to study the evolution of COVID-19 cases and deaths
on a global scale. Finally, we introduce \emph{three-way inconsistency
analysis} to determine anomalous countries with respect to three attributes:
countries' COVID-19 cases, deaths and human development indices. We demonstrate
the most anomalous countries across these three measures are Pakistan, the
United States and the United Arab Emirates.