  Gene genealogies are frequently studied by measuring properties such as their
height ($H$), length ($L$), sum of external branches ($E$), sum of internal
branches ($I$), and mean of their two basal branches ($B$), and the coalescence
times that contribute to the other genealogical features ($T$). These tree
properties and their relationships can provide insight into the effects of
population-genetic processes on genealogies and genetic sequences. Here, under
the coalescent model, we study the 15 correlations among pairs of features of
genealogical trees: $H_n$, $L_n$, $E_n$, $I_n$, $B_n$, and $T_k$ for a sample
of size $n$, with $2 \leq k \leq n$. We report high correlations among $H_n$,
$L_n$, $I_n,$ and $B_n$, with all pairwise correlations of these quantities
having values greater than or equal to $\sqrt{6} [6 \zeta(3) + 6 - \pi^2] / (
\pi \sqrt{18 + 9\pi^2 - \pi^4}) \approx 0.84930$ in the limit as $n \rightarrow
\infty$. Although $E_n$ has an expectation of 2 for all $n$ and $H_n$ has
expectation 2 in the limit as $n \rightarrow \infty$, their limiting
correlation is 0. The results contribute toward understanding features of the
shapes of coalescent trees.


  The 2019 coronavirus disease (COVID-19) became a worldwide pandemic with
currently no effective antiviral drug except treatments for symptomatic
therapy. Flux balance analysis is an efficient method to analyze metabolic
networks. It allows optimizing for a metabolic function and thus e.g.,
predicting the growth rate of a specific cell or the production rate of a
metabolite of interest. Here flux balance analysis was applied on human lung
cells infected with severe acute respiratory syndrome coronavirus 2
(SARS-CoV-2) to reposition metabolic drugs and drug combinations against the
replication of the SARS-CoV-2 virus within the host tissue. Making use of
expression data sets of infected lung tissue, genome-scale COVID-19-specific
metabolic models were reconstructed. Then host-specific essential genes and
gene-pairs were determined through in-silico knockouts that permit reducing the
viral biomass production without affecting the host biomass. Key pathways that
are associated with COVID-19 severity in lung tissue are related to oxidative
stress, as well as ferroptosis, sphingolipid metabolism, cysteine metabolism,
and fat digestion. By in-silico screening of FDA approved drugs on the putative
disease-specific essential genes and gene-pairs, 45 drugs and 99 drug
combinations were predicted as promising candidates for COVID-19 focused drug
repositioning (https://github.com/sysbiolux/DCcov). Among the 45 drug
candidates, six antiviral drugs were found and seven drugs that are being
tested in clinical trials against COVID-19. Other drugs like gemcitabine,
rosuvastatin and acetylcysteine, and drug combinations like
azathioprine-pemetrexed might offer new chances for treating COVID-19.


  There is an urgent need of biosynthetic bone grafts with enhanced osteogenic
capacity. In this study, we describe the design of hierarchical
meso-macroporous 3D-scaffolds based on mesoporous bioactive glasses (MBGs),
enriched with the peptide osteostatin and Zn2+ ions, and their osteogenic
effect on human mesenchymal stem cells (hMSCs) as a preclinical strategy in
bone regeneration. By using additive fabrication techniques, scaffolds
exhibiting hierarchical porosity: mesopores , macropores and big channels, were
prepared. These MBG scaffolds with or without osteostatin were evaluated in
cell cultures of hMSCs. Zinc promoted hMSCs colonization (both the surface and
inside) of MBG scaffolds. Moreover, Zn2+ ions and osteostatin together, but not
independently, in the scaffolds were found to induce the osteoblast
differentiation genes runt related transcription factor-2 (RUNX2) and alkaline
phosphatase (ALP) in hMSCs after 7 d of culture in the absence of an osteogenic
differentiation-promoting medium. These results add credence to the combined
use of zinc and osteostatin as an effective strategy for bone regeneration
applications.


  The estimate of the cumulative incidence of SARS-CoV-2 of 76% in Manaus by
October 2020 by Buss et al. relies on the assumption of an
exponentially-declining probability of sero-reversion over time. We present an
alternative, empirically-supported approach that is based on the observed
dynamics of antibody titers in sero-positive cases. Through this approach we
revise the cumulative incidence estimate to 66% (63.3% - 68.5%) by October
2020. This estimate has implications for the proximity to herd immunity and
future estimates of fitness advantages of virus variants, such as the P.1
variant of concern. This methodology is also relevant for any other sero-survey
analysis that has to adjust for sero-reversion.


  Cultural traits such as words, names, decorative styles, and technical
standards often assume arbitrary values and are thought to evolve neutrally.
But neutral evolution cannot explain why some traits come and go in cycles of
popularity while others become entrenched. Here we study frequency-dependent
selection (FDS)--where a trait's tendency to be copied depends on its current
frequency regardless of the trait value itself. We develop a maximum-likelihood
method to infer the precise form of FDS from time series of trait abundance,
and we apply the method to data on baby names and pet dog breeds over the last
century. We find that the most common names tend to decline by 2%-6% per yr on
average; whereas rare names--1 in 10,000 births--tend to increase by 1%-3% per
yr. This specific form of negative FDS explains patterns of diversity and
replicates across the United States, France, Norway and the Netherlands,
despite cultural, linguistic and demographic variation. We infer a fixed
fitness offset between male and female names that implies different rates of
innovation. We also find a strong selective advantage for biblical names in
every frequency class, which explains their predominance among the most common
names. In purebred dog registrations we infer a form of negative FDS that is
consistent with a preference for novelty, in which each year's newest breeds
outgrow the previous by about 1%/yr, which also recapitulates boom-bust cycles
in dog fanciers. Finally, we define the concept of effective
frequency-dependent selection, which enables a meaningful interpretation of
inferred FDS even for complex mechanisms of evolution. Our analysis generalizes
neutral evolution to incorporate pressures of conformity and anti-conformity as
fundamental forces in social evolution, and our inference procedure provides a
quantitative account of how these forces operate within and across cultures.


  During the COVID-19 pandemic, community containment measures have been
implemented in different parts of the world. However, despite government
recommendations and mandates, a progressive increase in mobility was observed
in several sites. This phenomenon is often referred to as "quarantine fatigue".
Mathematical models have been used for more than a century to study the
dynamics of diseases spread. Since the start of the COVID-19 pandemic, a large
number of models have been developed to investigate different aspects of the
pandemic. The aim of the present work has been to design and explore a
mathematical model of SARS-CoV-2 spread with quarantine fatigue. A
deterministic model represented by a system of differential equations, based on
Susceptible-Infected-Recovered (SIR) dynamics, was developed. The model was
explored by means of computational simulations. The results obtained show that,
in some scenarios, after a first prevalence peak, fatigue can cause a second
peak. Thus, the dynamic behavior of the model suggests that fatigue may have
contributed to the production of the second peak of infected people observed in
different parts of the world between 2020 and early 2021. The explorations
performed also showed that the model is sensitive to infection rate and other
parameters related to propagation dynamics. The developed model is relatively
simple and has important limitations for performing accurate predictions in the
real world. However, our study highlights the importance of considering
quarantine fatigue in the management of the current COVID-19 pandemic, as well
as of possible future pandemics caused by other pathogens.


  Transitive inference (TI) refers to social cognition that facilitates the
discernment of unknown relationships between individuals using known
relationships. It is extensively reported that TI evolves in animals living in
a large group because TI could assess relative rank without deducing all dyadic
relationships, which averts costly fights.The relationships in a large group
become so complex that social cognition may not be developed adequately to
handle such complexity. If all group members apply TI to all possible members
in the group, TI is supposed to require extremely highly developed cognitive
abilities especially in a large group. Instead of developing cognitive
abilities significantly, animals may apply simplified TI we call reference TI
in this study as heuristic approaches. The reference TI allows members to
recognize and remember social interactions only among a set of reference
members rather than all potential members. Our study assumes that information
processes in the reference TI comprises 1) the number of reference members
based on which individuals infer transitively, 2) the number of reference
members shared by the same strategists, and 3) memory capacity. We examined how
information processes evolve in a large group using evolutionary simulations in
the hawk-dove game. Information processes with almost any numbers of reference
members could evolve in a large group as long as the numbers of shared
reference member are high because information from the others' experiences is
shared. TI dominates immediate inference, which assesses relative rank on
direct interactions, because TI could establish social hierarchy more rapidly
applying information from others' experiences.


  At the end of 2020, policy responses to the SARS-CoV-2 outbreak have been
shaken by the emergence of virus variants. The emergence of these more
contagious, more severe, or even vaccine-resistant strains have challenged
worldwide policy interventions. Anticipating the emergence of these mutations
to plan ahead adequate policies, and understanding how human behaviors may
affect the evolution of viruses by coevolution, are key challenges. In this
article, we propose coevolution with genetic algorithms (GAs) as a credible
approach to model this relationship, highlighting its implications, potential
and challenges. We present a dual GA model in which both viruses aiming for
survival and policy measures aiming at minimising infection rates in the
population, competitively evolve. Simulation runs reproduce the emergence of
more contagious variants, and identifies the evolution of policy responses as a
determinant cause of this phenomenon. This coevolution opens new possibilities
to visualise the impact of governments interventions not only on outbreak
dynamics, but also on its evolution, to improve the efficacy of policies.


  Linked-Reads technologies, such as 10x Genomics, combine both the
high-quality and low cost of short-reads sequencing and a long-range
information, through the use of barcodes able to tag reads which originate from
a common long DNA fragment. This technology has been employed in a broad range
of applications including assembly or phasing of genomes, and structural
variant calling. However, to date, no tool or API dedicated to the manipulation
of Linked-Reads data exist. We introduce LRez, a C++ API and toolkit which
allows easy management of Linked-Reads data. LRez includes various
functionalities, for computing number of common barcodes between genomic
regions, extracting barcodes from BAM files, as well as indexing and querying
both BAM and FASTQ files to quickly fetch reads or alignments sharing one or
multiple barcodes. LRez can thus be used in a broad range of applications
requiring barcode processing, in order to improve their performances. LRez is
implemented in C++, supported on Linux platforms, and available under AGPL-3.0
License at https://github.com/morispi/LRez.


  DNA fingerprinting is a cornerstone for human identification in forensics,
where the sequence of highly polymorphic short tandem repeats (STRs) from an
individual is compared against a DNA database. This presents significant
privacy risks to individuals with DNA profiles in the database due to hacking
by malicious attackers who may access the data and misuse it for secondary
purposes. In this paper, we propose a novel cryptographic framework for jointly
encrypting DNA-based fingerprints (STRs) with other biometric data, for
example, facial images, such that the STRs and biometrics information of an
individual are revealed only when a positive match is found, i.e. the STRs act
as decryption keys. Specifically, when a search is performed on the encrypted
database using STR sequences of an individual in the database, a perfect match
generates the facial image and/ or other biometrics of the individual while the
lack of a match returns a null result. By jointly encrypting DNA fingerprints
and other biometrics using the unique STRs generated keys, our approach ensures
perfect privacy of the encrypted information with decryption of only the record
with STRs matching the query. This safeguards the information of other
individuals in the same database. The proposed approach can also be used to
securely authenticate the identity of individuals or biological material in
scenarios beyond forensics including tracking the identity of samples for
clinical genetics and cell therapies.


  Viral fusion proteins are attached to the membrane of enveloped viruses (a
group that includes Coronaviruses, Dengue, HIV and Influenza) and catalyze
fusion between the viral and host membranes, enabling the virus to insert its
genetic material into the host cell. Given the importance of these
biomolecules, this work presents a centralized database containing the most
relevant information on viral fusion proteins, available through a free-to-use
web server accessible through the URL https://viralfp.bio.di.uminho.pt/. This
web application contains several bioinformatic tools, such as Clustal sequence
alignment and Weblogo, including as well a machine learning-based tool capable
of predicting the location of fusion peptides (the component of fusion proteins
that inserts into the host's cell membrane) within the fusion protein sequence.
Given the crucial role of these proteins in viral infection, their importance
as natural targets of our immune system and their potential as therapeutic
targets, this web application aims to foster our ability to fight pathogenic
viruses.


  Several European countries have suspended the inoculation of the AstraZeneca
vaccine out of suspicion of causing deep vein thrombosis. In this letter we
report some Fermi estimates performed using a stochastic model aimed at making
a risk-benefit analysis of the interruption of the delivery of the AstraZeneca
vaccine in France and Italy. Our results clearly show that excess deaths due to
the interruption of the vaccination campaign injections largely overrun those
due to thrombosis even in worst case scenarios of frequency and gravity of the
vaccine side effects.


  At the last step of short read mapping, the candidate locations of the reads
on the reference genome are verified to compute their differences from the
corresponding reference segments using sequence alignment algorithms.
Calculating the similarities and differences between two sequences is still
computationally expensive since approximate string matching techniques
traditionally inherit dynamic programming algorithms with quadratic time and
space complexity. We introduce GateKeeper-GPU, a fast and accurate
pre-alignment filter that efficiently reduces the need for expensive sequence
alignment. GateKeeper-GPU provides two main contributions: first, improving the
filtering accuracy of GateKeeper(state-of-the-art lightweight pre-alignment
filter), second, exploiting the massive parallelism provided by the large
number of GPU threads of modern GPUs to examine numerous sequence pairs rapidly
and concurrently. GateKeeper-GPU accelerates the sequence alignment by up to
2.9x and provides up to 1.4x speedup to the end-to-end execution time of a
comprehensive read mapper (mrFAST). GateKeeper-GPU is available at
https://github.com/BilkentCompGen/GateKeeper-GPU


  Identifying the main environmental drivers of SARS-CoV-2 transmissibility in
the population is crucial for understanding current and potential future
outbursts of COVID-19 and other infectious diseases. To address this problem,
we concentrate on the basic reproduction number $R_0$, which is not sensitive
to testing coverage and represents transmissibility in an absence of social
distancing and in a completely susceptible population. While many variables may
potentially influence $R_0$, a high correlation between these variables may
obscure the result interpretation. Consequently, we combine Principal Component
Analysis with feature selection methods from several regression-based
approaches to identify the main demographic and meteorological drivers behind
$R_0$. We robustly obtain that country's wealth/development (GDP per capita or
Human Development Index) is the most important $R_0$ predictor at the global
level, probably being a good proxy for the overall contact frequency in a
population. This main effect is modulated by built-up area per capita
(crowdedness in indoor space), onset of infection (likely related to increased
awareness of infection risks), net migration, unhealthy living
lifestyle/conditions including pollution, seasonality, and possibly BCG
vaccination prevalence. Also, we argue that several variables that
significantly correlate with transmissibility do not directly influence $R_0$
or affect it differently than suggested by naive analysis.


  Cancers evolve from mutation of a single cell with sequential clonal and
subclonal expansion of somatic mutation acquisition. Inferring clonal and
subclonal structures from bulk or single cell tumor genomic sequencing data has
a huge impact on cancer evolution studies. Clonal state and mutational order
can provide detailed insight into tumor origin and its future development. In
the past decade, a variety of methods have been developed for subclonal
reconstruction using bulk tumor sequencing data. As these methods have been
developed in different programming languages and using different input data
formats, their use and comparison can be problematic. Therefore, we established
a web server for clonal and subclonal structure inference and evolution of
cancer genomic data (COSINE), which included 12 popular subclonal
reconstruction methods. We decomposed each method via a detailed workflow of
single processing steps with a user-friendly interface. To the best of our
knowledge, this is the first web server providing online subclonal inference,
including the most popular subclonal reconstruction methods. COSINE is freely
accessible at www.clab-cosine.net or http://bio.rj.run:48996/cun-web.


  Processing high-throughput DNA sequencing data of individuals or populations
requires stringing together independent software tools with many parameters,
often leading to non-reproducible pipelines and datasets. We developed
grenepipe to streamline this data processing, an all-in-one Snakemake workflow
from raw sequencing data to the end product of a table of individuals'
genotypes or population frequencies. Our pipeline allows users to select among
a range of popular software tools within a single configuration file,
automatically downloads and installs software and dependencies, and runs with
two command calls: to prepare and to run. It is highly optimized for
scalability in cluster environments and parallel computing, splitting data
tasks into manageable genomic sections and automatically consolidating the
outputs. grenepipe is published under the GPL-3 license, and freely available
at https://github.com/moiexpositoalonsolab/grenepipe.


  Recently discovered simple quantitative relations, known as bacterial growth
laws, hint on the existence of simple underlying principles at the heart of
bacterial growth. In this work, we provide a unifying picture on how these
known relations, as well as new relations that we derive, stems from a
universal autocatalytic network common to all bacteria, facilitating balanced
exponential growth of individual cells. We show that the core of the cellular
autocatalytic network is the transcription -- translation machinery -- in
itself an autocatalytic network comprising several coupled autocatalytic
cycles, including the ribosome, RNA polymerase, and tRNA charging cycles. We
derive two types of growth laws per autocatalytic cycle, one relating growth
rate to the relative fraction of the catalyst and its catalysis rate, and the
other relating growth rate to all the time scales in the cycle. The structure
of the autocatalytic network generates numerous regimes in state space,
determined by the limiting components, while the number of growth laws can be
much smaller. We also derive a growth law that accounts for the RNA polymerase
autocatalytic cycle, which we use to explain how growth rate depends on the
inducible expression of the rpoB and rpoC genes, which code for the RpoB and C
protein subunits of RNA polymerase, and how the concentration of rifampicin,
which targets RNA polymerase, affects growth rate without changing the
RNA-to-protein ratio. We derive growth laws for tRNA synthesis and charging,
and predict how growth rate depends on temperature, perturbation to ribosome
assembly, and membrane synthesis.


  Strong positive feedback is considered a necessary condition to observe
abrupt shifts of ecosystems. A few previous studies have shown that demographic
noise -- arising from the probabilistic and discrete nature of birth and death
processes in finite systems -- makes the transitions gradual or continuous. In
this paper, we show that demographic noise may, in fact, promote abrupt
transitions in systems that would otherwise show continuous transitions. We
present our methods and results in a tutorial-like format. We begin with a
simple spatially-explicit individual-based model with local births and deaths
influenced by positive feedback processes. We then derive a stochastic
differential equation that describes how local probabilistic rules scale to
stochastic population dynamics. The infinite-size well-mixed limit of this SDE
for our model is consistent with mean-field models of abrupt regime-shifts.
Finally, we analytically show that as a consequence of demographic noise,
finite-size systems can undergo abrupt shifts even with weak positive
interactions. Numerical simulations of our spatially-explicit model confirm
this prediction. Thus, we predict that small-sized populations and ecosystems
may undergo abrupt collapse even when larger systems - with the same
microscopic interactions - show a smooth response to environmental stress.


  We propose a simple rule of thumb for countries which have embarked on a
vaccination campaign while still facing the need to keep non-pharmaceutical
interventions (NPI) in place because of the ongoing spread of SARS-CoV-2. If
the aim is to keep the death rate from increasing, NPIs can be loosened when it
is possible to vaccinate more than twice the growth rate of new cases. If the
aim is to keep the pressure on hospitals under control, the vaccination rate
has to be about four times higher. These simple rules can be derived from the
observation that the risk of death or a severe course requiring hospitalization
from a COVID-19 infection increases exponentially with age and that the sizes
of age cohorts decrease linearly at the top of the population pyramid.
Protecting the over 60-year-olds, which constitute approximately one-quarter of
the population in Europe (and most OECD countries), reduces the potential loss
of life by 95 percent.


  Simulating the spread of infectious diseases in human communities is critical
for predicting the trajectory of an epidemic and verifying various policies to
control the devastating impacts of the outbreak. Many existing simulators are
based on compartment models that divide people into a few subsets and simulate
the dynamics among those subsets using hypothesized differential equations.
However, these models lack the requisite granularity to study the effect of
intelligent policies that influence every individual in a particular way. In
this work, we introduce a simulator software capable of modeling a population
structure and controlling the disease's propagation at an individualistic
level. In order to estimate the confidence of the conclusions drawn from the
simulator, we employ a comprehensive probabilistic approach where the entire
population is constructed as a hierarchical random variable. This approach
makes the inferred conclusions more robust against sampling artifacts and gives
confidence bounds for decisions based on the simulation results. To showcase
potential applications, the simulator parameters are set based on the formal
statistics of the COVID-19 pandemic, and the outcome of a wide range of control
measures is investigated. Furthermore, the simulator is used as the environment
of a reinforcement learning problem to find the optimal policies to control the
pandemic. The obtained experimental results indicate the simulator's
adaptability and capacity in making sound predictions and a successful policy
derivation example based on real-world data. As an exemplary application, our
results show that the proposed policy discovery method can lead to control
measures that produce significantly fewer infected individuals in the
population and protect the health system against saturation.


  In this paper we show how inter-cellular molecular communication may change
the overall levels of photosynthesis in plants. Individual plant cells respond
to external stimuli, such as illumination levels, to regulate their
photosynthetic output. Here, we present a mathematical model which shows that
by sharing information internally using molecular communication, plants may
increase overall photosynthate production. Numerical results show that higher
mutual information between cells corresponds to an increase in overall
photosynthesis by as much as 25 per cent. This suggests that molecular
communication plays a vital role in maximising the photosynthesis in plants and
therefore suggests new routes to influence plant development in agriculture and
elsewhere.


  Solid estimates describing the clinical course of SARS-CoV-2 infections are
still lacking due to under-ascertainment of asymptomatic and mild-disease
cases. In this work, we quantify age-specific probabilities of transitions
between stages defining the natural history of SARS-CoV-2 infection from 1,965
SARS-CoV-2 positive individuals identified in Italy between March and April
2020 among contacts of confirmed cases. Infected contacts of cases were
confirmed via RT-PCR tests as part of contact tracing activities or
retrospectively via IgG serological tests and followed-up for symptoms and
clinical outcomes. In addition, we provide estimates of time intervals between
key events defining the clinical progression of cases as obtained from a larger
sample, consisting of 95,371 infections ascertained between February and July
2020. We found that being older than 60 years of age was associated with a
39.9% (95%CI: 36.2-43.6%) likelihood of developing respiratory symptoms or
fever >= 37.5 {\deg}C after SARS-CoV-2 infection; the 22.3% (95%CI: 19.3-25.6%)
of the infections in this age group required hospital care and the 1% (95%CI:
0.4-2.1%) were admitted to an intensive care unit (ICU). The corresponding
proportions in individuals younger than 60 years were estimated at 27.9%
(95%CI: 25.4-30.4%), 8.8% (95%CI: 7.3-10.5%) and 0.4% (95%CI: 0.1-0.9%),
respectively. The infection fatality ratio (IFR) ranged from 0.2% (95%CI:
0.0-0.6%) in individuals younger than 60 years to 12.3% (95%CI: 6.9-19.7%) for
those aged 80 years or more; the case fatality ratio (CFR) in these two age
classes was 0.6% (95%CI: 0.1-2%) and 19.2% (95% CI: 10.9-30.1%), respectively.
The median length of stay in hospital was 10 (IQR 3-21) days; the length of
stay in ICU was 11 (IQR 6-19) days. The obtained estimates could be
instrumental to refine mathematical modeling work supporting public health
decisions.


  Prediction of the progression of an infectious disease outbreak is important
for planning and coordinating a response. Differential equations are often used
to model an epidemic outbreak's behaviour but are challenging to parameterise.
Furthermore, these models can suffer from misspecification, which biases
predictions and parameter estimates. Stochastic models can help with
misspecification but are even more expensive to simulate and perform inference
with. Here, we develop an explicitly likelihood-based variation of the
generalised profiling method as a tool for prediction and inference under model
misspecification. Our approach allows us to carry out identifiability analysis
and uncertainty quantification using profile likelihood-based methods without
the need for marginalisation. We provide justification for this approach by
introducing a new interpretation of the model approximation component as a
stochastic constraint. This preserves the rationale for using profiling rather
than integration to remove nuisance parameters while also providing a link back
to stochastic models. We applied an initial version of this method during an
outbreak of measles in Samoa in 2019-2020 and found that it achieved relatively
fast, accurate predictions. Here we present the most recent version of our
method and its application to this measles outbreak, along with additional
validation.


  Pathogens like the SARS-CoV-2 are transmitted not only through violent
expiratory events like coughing, but also through routine activities like
breathing/speaking/singing. We perform direct numerical simulations of the
turbulent transport of potentially infectious aerosols in short conversations.
It is shown that a two-way conversation significantly reduces the aerosol
exposure compared to a relative monologue by one person and relative silence of
the other. This is because the interaction of the jets ejected from the mouth
of each speaker produce a "canceling" effect. Unequal conversation is shown to
significantly increase the risk of infection to the person who talks less.
Interestingly, a small height difference is worse for infection spread, due to
reduced interference between the two speech jets, than two faces at the same
level! For small axial separation, speech jets show large oscillations and
reach the other person intermittently. We suggest a range of lateral
separations between two people to minimize transmission risk. A realistic
estimate of the infection probability is provided by including exposure through
eyes and mouth, in addition to the more common method of using inhaled virions
alone. We expect that our results will provide useful inputs to epidemiological
models and to disease management.


  Sea fishing is a highly mobile activity, favoured by the vastness of the
oceans, the absence of physical boundaries and the abstraction of legislative
boundaries. Understanding and anticipating this mobility is a major challenge
for fisheries management issues, both at the national and international levels.
''FisherMob'' is a free Gama tool designed to study the effect of economic and
biological factors on the dynamics of connected fisheries. It incorporate the
most important processes involved in fisheries dynamics: fish abundance
variability, price of the fishing effort and ex-vessel fish market price that
which depends on the ratio between offer and demand. The tool uses as input a
scheme of a coastal area with delimited fishing sites, fish biological
parameters and fisheries parameters. It runs with a userfriendly graphic
interface and generates output files that can be post-processed easily using
graphic and statistical software.


  This paper is concerned with the design of intermittent non-pharmaceutical
strategies to mitigate the spread of the COVID-19 epidemic exploiting network
epidemiological models. Specifically, by studying a variational equation for
the dynamics of the infected in a network model of the epidemic spread, we
derive, using contractivity arguments, a condition that can be used to
guarantee that, in epidemiological terms, the effective reproduction number is
less than unity. This condition has three advantages: (i) it is easily
computable; (ii) it is directly related to the model parameters; (iii) it can
be used to enforce a scalability condition that prohibits the amplification of
disturbances within the network system. We then include satisfaction of such a
condition as a constraint in a Model Predictive Control problem so as to
mitigate (or suppress) the spread of the epidemic while minimizing the economic
impact of the interventions. A data-driven model of Italy as a network of three
macro-regions (North, Center, and South), whose parameters are identified from
real data, is used to illustrate and evaluate the effectiveness of the proposed
control strategy.


  Biological information processing is generally assumed to be classical.
Measured cellular energy budgets of both prokaryotes and eukaryotes, however,
fall orders of magnitude short of the power required to maintain classical
states of protein conformation and localization at the \AA, fs scales predicted
by single-molecule decoherence calculations and assumed by classical molecular
dynamics models. We suggest that decoherence is limited to the immediate
surroundings of the cell membrane and of intercompartmental boundaries within
the cell, and that bulk cellular biochemistry implements quantum information
processing. Detection of Bell-inequality violations in responses to
perturbation of recently-separated sister cells would provide a sensitive test
of this prediction. If it is correct, modeling both intra- and intercellular
communication requires quantum theory.


  Carlo Cipolla's stupidity quadrant and his five laws of stupidity were
proposed for the first time in 1976. Exposed in a humorous mood by the author,
these concepts nevertheless describe very serious features of the interactions
among human beings. Here, we propose a new interpretation of Cipolla's ideas in
a biophysical framework, using the well-known predator-prey or "Lotka-Volterra"
model. We find that there is indeed a correspondence between Cipolla's
approach, based on economics, and biophysical economics. On the basis of this
examination, we propose a sixth law of stupidity, additional to the five
proposed by Cipolla. The law states that humans are the stupidest species in
the ecosystem.


  The flexibility and stiffness of small DNA play a fundamental role ranging
from several biophysical processes to nano-technological applications. Here, we
estimate the mechanical properties of short double-stranded DNA (dsDNA) having
length ranging from 12 base-pairs (bps) to 56 bps, paranemic crossover (PX)
DNA, and hexagonal DNA nanotubes (DNTs) using two widely used coarse-grain
models $-$ Martini and oxDNA. To calculate the persistence length ($L_p$) and
the stretch modulus ($\gamma$) of the dsDNA, we incorporate the worm-like chain
and elastic rod model, while for DNT, we implement our previously developed
theoretical framework. We compare and contrast all the results with previously
reported all-atom molecular dynamics (MD) simulation and experimental results.
The mechanical properties of dsDNA ($L_p$ $\sim$ 50nm, $\gamma \sim$ 800-1500
pN), PX DNA ($\gamma \sim$ 1600-2000 pN) and DNTs ($L_p \sim 1-10\ \mu$m,
$\gamma \sim$ 6000-8000 pN) estimated using Martini soft elastic network and
oxDNA are in very good agreement with the all-atom MD and experimental values,
while the stiff elastic network Martini reproduces order of magnitude higher
values of $L_p$ and $\gamma$. The high flexibility of small dsDNA is also
depicted in our calculations. However, Martini models proved inadequate to
capture the salt concentration effects on the mechanical properties with
increasing salt molarity. OxDNA captures the salt concentration effect on small
dsDNA mechanics. But it is found to be ineffective to reproduce the
salt-dependent mechanical properties of DNTs. Also, unlike Martini, the time
evolved PX DNA and DNT structures from the oxDNA models are comparable to the
all-atom MD simulated structures. Our findings provide a route to study the
mechanical properties of DNA nanostructures with increased time and length
scales and has a remarkable implication in the context of DNA nanotechnology.


  Molecular structures of RNA molecules reconstructed from X-ray
crystallography frequently contain errors. Motivated by this problem we examine
clustering on a torus since RNA shapes can be described by dihedral angles. A
previously developed clustering method for torus data involves two tuning
parameters and we assess clustering results for different parameter values in
relation to the problem of so-called RNA clashes. This clustering problem is
part of the dynamically evolving field of statistics on manifolds. Statistical
problems on the torus highlight general challenges for statistics on manifolds.
Therefore, the torus PCA and clustering methods we propose make an important
contribution to directional statistics and statistics on manifolds in general.


  Likelihood-based methods are widely considered the best approaches for
reconstructing ancestral states. Although much effort has been made to study
properties of these methods, previous works often assume that both the tree
topology and edge lengths are known. In some scenarios the tree topology might
be reasonably well known for the taxa under study. When sequence length is much
smaller than the number of species, however, edge lengths are not likely to be
accurately estimated. We study the consistency of the maximum likelihood and
empirical Bayes estimators of ancestral state of discrete traits in such
settings under a star tree. We prove that the likelihood-based reconstruction
is consistent under symmetric models but can be inconsistent under
non-symmetric models. We show, however, that a simple consistent estimator for
the ancestral states is available under non-symmetric models. The results
illustrate that likelihood methods can unexpectedly have undesirable properties
as the number of sequences considered get very large. Broader implications of
the results are discussed.


  Objective: Human-curated disease ontologies are widely used for diagnostic
evaluation, treatment and data comparisons over time, and clinical decision
support. The classification principles underlying these ontologies are guided
by the analysis of observable pathological similarities between disorders,
often based on anatomical or histological principles. Although, thanks to
recent advances in molecular biology, disease ontologies are slowly changing to
integrate the etiological and genetic origins of diseases, nosology still
reflects this "reductionist" perspective. Proximity relationships of disease
modules (hereafter DMs) in the human interactome network are now increasingly
used in diagnostics, to identify pathobiologically similar diseases and to
support drug repurposing and discovery. On the other hand, similarity relations
induced from structural proximity of DMs also have several limitations, such as
incomplete knowledge of disease-gene relationships and reliability of clinical
trials to assess their validity. The purpose of the study described in this
paper is to shed more light on disease similarities by analyzing the
relationship between categorical proximity of diseases in human-curated
ontologies and structural proximity of the related DM in the interactome.
Method: We propose a methodology (and related algorithms) to automatically
induce a hierarchical structure from proximity relations between DMs, and to
compare this structure with a human-curated disease taxonomy. Results: We
demonstrate that the proposed methodology allows to systematically analyze
commonalities and differences among structural and categorical similarity of
human diseases, help refine and extend human disease classification systems,
and identify promising network areas where new disease-gene interactions can be
discovered.


  In many biological systems, natural selection acts simultaneously on multiple
levels of organization. This scenario typically presents an evolutionary
conflict between the incentive of individuals to cheat and the collective
incentive to establish cooperation within a group. Generalizing previous work
on multilevel selection in evolutionary game theory, we consider a hyperbolic
PDE model of a group-structured population, in which members within a single
group compete with each other for individual-level replication; while the group
also competes against other groups for group-level replication. We derive a
threshold level of the relative strength of between-group competition such that
defectors take over the population below the threshold while cooperation
persists in the long-time population above the threshold. Under stronger
assumptions on the initial distribution of group compositions, we further prove
that the population converges to a steady state density supporting cooperation
for between-group selection strength above the threshold. We further establish
long-time bounds on the time-average of the collective payoff of the
population, showing that the long-run population cannot outperform the payoff
of a full-cooperator group even in the limit of infinitely-strong between-group
competition. When the group replication rate is maximized by an intermediate
level of within-group cooperation, individual-level selection casts a long
shadow on the dynamics of multilevel selection: no level of between-group
competition can erase the effects of the individual incentive to defect. We
further extend our model to study the case of multiple types of groups, showing
how the games that groups play can coevolve with the level of cooperation.


  Detecting predictive biomarkers from multi-omics data is important for
precision medicine, to improve diagnostics of complex diseases and for better
treatments. This needs substantial experimental efforts that are made difficult
by the heterogeneity of cell lines and huge cost. An effective solution is to
build a computational model over the diverse omics data, including genomic,
molecular, and environmental information. However, choosing informative and
reliable data sources from among the different types of data is a challenging
problem. We propose DIVERSE, a framework of Bayesian importance-weighted tri-
and bi-matrix factorization(DIVERSE3 or DIVERSE2) to predict drug responses
from data of cell lines, drugs, and gene interactions. DIVERSE integrates the
data sources systematically, in a step-wise manner, examining the importance of
each added data set in turn. More specifically, we sequentially integrate five
different data sets, which have not all been combined in earlier bioinformatic
methods for predicting drug responses. Empirical experiments show that DIVERSE
clearly outperformed five other methods including three state-of-the-art
approaches, under cross-validation, particularly in out-of-matrix prediction,
which is closer to the setting of real use cases and more challenging than
simpler in-matrix prediction. Additionally, case studies for discovering new
drugs further confirmed the performance advantage of DIVERSE.


  In Section 1, we revisit the partial differential equation (PDE) for the
probability generating function (PGF) of the time-nonhomogeneous BDI
(birth-and-death-with-immigration) process and derive a closed form solution.
To the best of our knowledge, this is a new mathematical result. We state this
result as Proposition 1. We state as Corollary 1 that the negative binomial
distribution of the time-homogeneous BDI process discussed in Part I extends to
the general time-nonhomogeneous case, provided that the ratio of the
immigration rate to the birth rate is a constant. In section 1.2, we take up
the heuristic approach discussed by Bartlett and Bailey (1964), and carry it
out to completion by arriving at the solution obtained above,.
  In Section 2, we present the results of our extensive simulation experiments
of the time-nonhomogeneous BD process that was analyzed in Part III-A and
confirm our analytic results. In Section 3, we undertake similar simulation
experiments for the BDI process that is analyzed in Section 1.
  As we discuss in Section 4, our stochastic model now seems more promising and
powerful than has been heretofore expected. In Appendix B, a closed form
solution for the M(t)/M(t)/infinity queue is obtained, as a special case of
this BDI process model.


  Clinical trials of a vaccine during an epidemic face particular challenges,
such as the pressure to identify an effective vaccine quickly to control the
epidemic, and the effect that time-space-varying infection incidence has on the
power of a trial. We illustrate how the operating characteristics of different
trial design elements may be evaluated using a network epidemic and trial
simulation model, based on COVID-19 and individually randomised two-arm trials
with a binary outcome. We show that "ring" recruitment strategies, prioritising
participants at high risk of infection, can result in substantial improvement
in terms of power, if sufficiently many contacts of observed cases are at high
risk. In addition, we introduce a novel method to make more efficient use of
the data from the earliest cases of infection observed in the trial, whose
infection may have been too early to be vaccine-preventable. Finally, we
compare several methods of response-adaptive randomisation, discussing their
advantages and disadvantages in this two-arm context and identifying particular
adaptation strategies that preserve power and estimation properties, while
slightly reducing the number of infections, given an effective vaccine.


  We consider the classical problem of optimal portfolio construction with the
constraint that no short position is allowed, or equivalently the valid
equilibria of multispecies Lotka-Volterra equations with self-regulation in the
special case where the interaction matrix is of unit rank, corresponding to
species competing for a common resource. We compute the average number of
solutions and show that its logarithm grows as $N^\alpha$, where $N$ is the
number of assets or species and $\alpha \leq 2/3$ depends on the interaction
matrix distribution. We conjecture that the most likely number of solutions is
much smaller and related to the typical sparsity $m(N)$ of the solutions, which
we compute explicitly. We also find that the solution landscape is similar to
that of spin-glasses, i.e. very different configurations are quasi-degenerate.
Correspondingly, "disorder chaos" is also present in our problem. We discuss
the consequence of such a property for portfolio construction and ecologies,
and question the meaning of rational decisions when there is a very large
number "satisficing" solutions.


  A key question concerning collective decisions is whether a social system can
settle on the best available option when some members learn from others instead
of evaluating the options on their own. This question is challenging to study,
and previous research has reached mixed conclusions, because collective
decision outcomes depend on the insufficiently understood complex system of
cognitive strategies, task properties, and social influence processes. This
study integrates these complex interactions together in one general yet
partially analytically tractable mathematical framework using a dynamical
system model. In particular, it investigates how the interplay of the
proportion of social learners, the relative merit of options, and the type of
conformity response affect collective decision outcomes in a binary choice. The
model predicts that when the proportion of social learners exceeds a critical
threshold, a bi-stable state appears in which the majority can end up favoring
either the higher- or lower-merit option, depending on fluctuations and initial
conditions. Below this threshold, the high-merit option is chosen by the
majority. The critical threshold is determined by the conformity response
function and the relative merits of the two options. The study helps reconcile
disagreements about the effect of social learners on collective performance and
proposes a mathematical framework that can be readily adapted to extensions
investigating a wider variety of dynamics.


  West Nile virus (WNV) is a vector-borne pathogen of global relevance and is
currently the most widely distributed flavivirus of encephalitis worldwide.
This virus infects birds, humans, horses, and other mammals, and its
transmission cycle occurs in urban and rural areas. Climate conditions have
direct and indirect impacts on vector abundance and virus dynamics within the
mosquito. The significance of environmental variables as drivers in WNV
epidemiology is increasing under the current climate change scenario. In this
study, we used a machine learning algorithm to model WNV distributions in South
America. Our model evaluated eight environmental variables (type of biome,
annual temperature, seasonality of temperature, daytime temperature variation,
thermal amplitude, seasonality of precipitation, annual rainfall, and
elevation) for their contribution to the occurrence of WNV since its
introduction in South America (2004). Our results showed that environmental
variables can directly alter the occurrence of WNV, with lower precipitation
and higher temperatures associated with increased virus incidence. High-risk
areas may be modified in the coming years, becoming more evident with high
greenhouse gas emission levels. Countries such as Bolivia and Paraguay will be
greatly affected, drastically changing their current WNV distribution. Several
Brazilian areas will also increase the likelihood of presenting WNV, mainly in
the Northeast and Midwest regions and the Pantanal biome. The Galapagos Islands
will also probably increase their geographic range suitable for WNV occurrence.
It is necessary to develop preventive policies to minimize potential WNV
infection in humans and enhance active epidemiological surveillance in birds,
humans, and other mammals before it becomes a more significant public health
problem in South America.


  A widely used tool for analysing the Covid-19 pandemic is the standard SIR
model. It seems often to be used as a black box, not taking into account that
this model was derived as a special case of the seminal Kermack-McKendrick
theory from 1927. This is our starting point. We explain the setup of the
Kermack-McKendrick theory (passing to a discrete approach) and use medical
information for specializing to a model which we call {\em adapted
K-McK-model}. This includes effects of vaccination, mass testing and mutants.
We demonstrate the use of the model by applying it to the development in
Germany. As a striking application we demonstrate that a comparatively mild
intervention reducing the time until quarantine by one day leads to a drastic
improvement. A similar effect can be obtained by certain mass testings as we
will demonstrate. We discuss possibilities to apply the model both for
predictions and as an analysis tool. We compare the adapted K-McK-model to the
standard SIR model and observe condiderable differences if the contact rates
are not constant. Finally we compare the model reproduction rate with the
empirical reproduction rate determined by the Robert Koch-Institut.


  Prokaryote mobilome genomes rely on host machineries for survival and
replication. Given that mobile genetic elements (MGEs) derive their energy from
host cells, we investigated the diversity of ATP-utilizing proteins in MGE
genomes to determine whether they might be associated with proteins that could
suppress related host proteins that consume host energy. A comprehensive search
of 353 huge phage genomes revealed that up to 9% of the proteins have ATPase
domains. For example, ATPase proteins constitute ~3% of the genomes of Lak
phages with ~550 kbp genomes that occur in the microbiomes of humans and other
animals. Statistical analysis shows the number of ATPase proteins increases
linearly with genome length, consistent with a large sink for host ATP during
replication of megaphages. Using metagenomic data from diverse environments, we
found 505 mobilome proteins with ATPase domains fused to diverse functional
domains. Among these composite ATPase proteins, 61.6% have known functional
domains that could contribute to host energy diversion during the mobilome life
cycle. As many have domains that are known to interact with nucleic acids and
proteins, we infer that numerous ATPase proteins are used during replication
and for protection from host immune systems. We found a set of uncharacterized
ATPase proteins with nuclease and protease activities, displaying unique domain
architectures that are energy intensive based on the presence of multiple
ATPase domains. In many cases, these composite ATPase proteins genomically
co-localize with small proteins in genomic contexts that are reminiscent of
toxin-antitoxin systems. Small proteins that function as inhibitors may be a
common strategy for control of cellular processes, thus could inspire the
development of new nucleic acid and protein manipulation tools, with diverse
biotechnological applications.


  In this white paper we introduce Helix, an AI based solution for missense
pathogenicity prediction. With recent advances in the sequencing of human
genomes, massive amounts of genetic data have become available. This has
shifted the burden of labor for genetic diagnostics and research from the
gathering of data to its interpretation. Helix presents a state of the art
platform for the prediction of pathogenicity in human missense variants. In
addition to offering best-in-class predictive performance, Helix offers a
platform that allows researchers to analyze and interpret variants in depth
that can be accessed at helixlabs.ai.


  Background: Following the outbreak of the coronavirus epidemic in early 2020,
municipalities, regional governments and policymakers worldwide had to plan
their Non-Pharmaceutical Interventions (NPIs) amidst a scenario of great
uncertainty. At this early stage of an epidemic, where no vaccine or medical
treatment is in sight, algorithmic prediction can become a powerful tool to
inform local policymaking. However, when we replicated one prominent
epidemiological model to inform health authorities in a region in the south of
Brazil, we found that this model relied too heavily on manually predetermined
covariates and was too reactive to changes in data trends.
  Methods: Our four proposed variations of the original method allow accessing
data of daily reported infections and take into account the under-reporting of
cases more explicitly. Two of the proposed versions also attempt to model the
delay in test reporting. We simulated weekly forecasting of deaths from the
period from 31/05/2020 until 31/01/2021. That workflow allowed us to run a
lighter version of the model after the first calibration week. Google Mobility
data, weekly updated, were used as covariates to the model at each simulated
run.
  Findings: The changes made the model significantly less reactive and more
rapid in adapting to scenarios after a peak in deaths is observed. Assuming
that reported cases were under-reported greatly benefited the model in its
stability, and modelling retroactively-added data (due to the "hot" nature of
the data used) had a negligible impact on performance.
  Interpretation: Although not as reliable as death statistics, case
statistics, when modelled in conjunction with an "overestimate" parameter,
provide a good alternative for improving the forecasting of models, especially
in long-range predictions and after the peak of an infection wave.


  Local coexistence of species in large ecosystems is traditionally explained
within the broad framework of niche theory. However, its rationale hardly
justifies rich biodiversity observed in nearly homogeneous environments. Here
we consider a consumer-resource model in which a coarse-graining procedure
accounts for a variety of ecological mechanisms and leads to effective spatial
effects which favour species coexistence. Herein, we provide conditions for
several species to live in an environment with very few resources. In fact, the
model displays two different phases depending on whether the number of
surviving species is larger or smaller than the number of resources. We obtain
conditions whereby a species can successfully colonize a pool of coexisting
species. Finally, we analytically compute the distribution of the population
sizes of coexisting species. Numerical simulations as well as empirical
distributions of population sizes support our analytical findings.


  Computational models of biological processes provide one of the most powerful
methods for a detailed analysis of the mechanisms that drive the behavior of
complex systems. Logic-based modeling has enhanced our understanding and
interpretation of those systems. Defining rules that determine how the output
activity of biological entities is regulated by their respective inputs has
proven to be challenging, due to increasingly larger models and the presence of
noise in data, allowing multiple model parameterizations to fit the
experimental observations.
  We present several Boolean function metrics that provide modelers with the
appropriate framework to analyze the impact of a particular model
parameterization. We demonstrate the link between a semantic characterization
of a Boolean function and its consistency with the model's underlying
regulatory structure. We further define the properties that outline such
consistency and show that several of the Boolean functions under study violate
them, questioning their biological plausibility and subsequent use. We also
illustrate that regulatory functions can have major differences with regard to
their asymptotic output behavior, with some of them being biased towards
specific Boolean outcomes when others are dependent on the ratio between
activating and inhibitory regulators.
  Application results show that in a specific signaling cancer network, the
function bias can be used to guide the choice of logical operators for a model
that matches data observations. Moreover, graph analysis indicates that the
standardized Boolean function bias becomes more prominent with increasing
numbers of regulators, confirming the fact that rule specification can
effectively determine regulatory outcome despite the complex dynamics of
biological networks.


  In this study, the mechanical behavior of single tau, dimerized tau, and
tau-microtubule interface subjected to high strain rate is investigated by
molecular dynamics simulation.


  Cooperation in an open dynamic system fundamentally depends upon information
distributed across its components. Yet in an environment with rapidly enlarging
complexity, this information may need to change adaptively to enable not only
cooperation but also the mere survival of an organism. Combining the methods of
evolutionary game theory, agent-based simulation, and statistical physics, we
develop a model of the evolution of cooperation in an ageing population of
artificial decision makers playing spatial tag-mediated prisoner's dilemma
games with their ingroup neighbors and with genetically unrelated immigrant
agents. In our model with six strategies we introduce the concept of
time-varying tags such that the phenotypic features of 'new' agents can change
into 'approved' following variable approval times. Our Monte Carlo simulations
show that ingroup-biased ethnocentric cooperation can dominate only at low
costs and short approval times. In the standard 4-strategy model with fixed
tags, we identified a critical cost $c_{\mathrm{crit}}$ above which cooperation
transitioned abruptly into the phase of pure defection, revealing remarkable
fragility of ingroup-biased generosity. In our generalized 6-strategy model
with time-varying tags, elevated cooperation was observed for a wider region of
the parameter space, peaking at intermediate approval times and cost values
above $c_{\mathrm{crit}}$. Our findings show that in an open system subject to
immigration dynamics, high levels of social cooperation are possible if a
fraction of the population adopts the strategy with an egalitarian generosity
directed towards both native and approved naturalized citizens, regardless of
their actual origin. These findings also suggest that instead of relying upon
arbitrarily fixed approval times, there is an optimal duration of the
naturalization procedure from which the society as a whole can profit most.


  Hyperactivity is one of the hallmakrs of ADHD. Aberrant dopamine signaling is
a major theme in ADHD and dopamine production is directly linked to the
intensity and persistence of hyperactive conduct. The strength and persistence
of hyperactivity responses in Drosophila to startle stimuli were measured in a
study to determine the effects of sugar on dopamine development. A total of
four experimental groups, namely 1%,3%, and 5% glucose, as well as a control
group were taken for the diet of Drosophila, and these four different amounts
of glucose were introduced to the growth medium where Drosophila was cultured.
The movements of Drosophila in the four treatment groups were captured using a
camera. This experiment was carried out five times, each time using a different
batch of Drosophila. Each group's average velocity over time was also reported.
The web adaptation of Drosophila Activity Monitor(DAM) was used to analyze the
captured movies from the camera. Furthermore, when it came to hyperactivity
persistence, all four treatment classes were statistically different (p0.05).
Since the strength and persistence of hyperactive behavior are directly
correlated to dopamine output, this study shows that higher glucose intake is
associated with more hyperactivity, for both the intensity({\Delta}V) and
persistence.


  Single-Cell RNA sequencing (scRNA-seq) measurements have facilitated
genome-scale transcriptomic profiling of individual cells, with the hope of
deconvolving cellular dynamic changes in corresponding cell sub-populations to
better understand molecular mechanisms of different development processes.
Several scRNA-seq analysis methods have been proposed to first identify cell
sub-populations by clustering and then separately perform differential
expression analysis to understand gene expression changes. Their corresponding
statistical models and inference algorithms are often designed disjointly. We
develop a new method -- SimCD -- that explicitly models cell heterogeneity and
dynamic differential changes in one unified hierarchical gamma-negative
binomial (hGNB) model, allowing simultaneous cell clustering and differential
expression analysis for scRNA-seq data. Our method naturally defines cell
heterogeneity by dynamic expression changes, which is expected to help achieve
better performances on the two tasks compared to the existing methods that
perform them separately. In addition, SimCD better models dropout (zero
inflation) in scRNA-seq data by both cell- and gene-level factors and obviates
the need for sophisticated pre-processing steps such as normalization, thanks
to the direct modeling of scRNA-seq count data by the rigorous hGNB model with
an efficient Gibbs sampling inference algorithm. Extensive comparisons with the
state-of-the-art methods on both simulated and real-world scRNA-seq count data
demonstrate the capability of SimCD to discover cell clusters and capture
dynamic expression changes. Furthermore, SimCD helps identify several known
genes affected by food deprivation in hypothalamic neuron cell subtypes as well
as some new potential markers, suggesting the capability of SimCD for
bio-marker discovery.


  In this study we document the unexpected discovery of multiple coronaviruses
and a BSL-3 pathogen in agricultural cotton and rice sequencing datasets. In
particular, we have identified a novel HKU5-related Merbecovirus in a cotton
dataset sequenced by the Huazhong Agricultural University in 2017. We have also
found an infectious clone sequence containing a novel HKU4-related Merbecovirus
related to MERS coronavirus in a rice dataset sequenced by the Huazhong
Agricultural University in early 2020. Another HKU5-related Merbecovirus, as
well as Japanese encephalitis virus, were identified in a cotton dataset
sequenced by the Huazhong Agricultural University in 2018. An HKU3-related
Betacoronavirus was found in a Mus musculus sequencing dataset from the Wuhan
Institute of Virology in 2017. Finally, a SARS-WIV1-like Betacoronavirus was
found in a rice dataset sequenced by the Fujian Agriculture and Forestry
University in 2017. Using the contaminating reads we have extracted from the
above datasets, we were able to assemble complete genomes of two novel
coronaviruses which we disclose herein. In light of our findings, we raise
concerns about biosafety protocol breaches, as indicated by our discovery of
multiple dangerous human pathogens in agricultural sequencing laboratories in
Wuhan and Fouzou City, China.


  Dispersal-induced growth (DIG) occurs when two populations with time-varying
growth rates, each of which, when isolated, would become extinct, are able to
persist and grow exponentially when dispersal among the two populations is
present. This work provides a mathematical exploration of this surprising
phenomenon, in the context of a deterministic model with periodic variation of
growth rates, and characterizes the factors which are important in generating
the DIG effect and the corresponding conditions on the parameters involved.


  The idea of a metapopulation has become canonical in ecology. Its original
mean field form provides the important intuition that migration and extinction
interact to determine the dynamics of a population composed of subpopulations.
From its conception, it has been evident that the very essence of the
metapopulation paradigm centers on the process of local extinction. We note
that there are two qualitatively distinct types of extinction, gradual and
catastrophic, and explore their impact on the dynamics of metapopulation
formation using discrete iterative maps. First, by modifying the classic
logistic map with the addition of the Allee effect, we show that catastrophic
local extinctions in subpopulations are a pre-requisite of metapopulation
formation. When subpopulations experience gradual extinction, increased
migration rates force synchrony and drive the metapopulation below the Allee
point resulting in migration induced destabilization of the system across
parameter space. Second, a sawtooth map (an extension of the Bernoulli bit
shift map) is employed to simultaneously explore the increasing and decreasing
modes of population behavior. We conclude with four generalizations. 1. At low
migration rates, a metapopulation may go extinct faster than completely
unconnected subpopulations. 2. There exists a gradient between stable
metapopulation formation and population synchrony, with critical transitions
from no metapopulation to metapopulation to synchronization, the latter
frequently inducing metapopulation extinction. 3. Synchronization patterns
emerge through time, resulting in synchrony groups and chimeric populations
existing simultaneously. 4. There are two distinct mechanisms of
synchronization: i. extinction and rescue and, ii.) stretch reversals in a
modification of the classic chaotic stretching and folding.


  Microbial communities are widely studied using high-throughput sequencing
techniques, such as 16S rRNA gene sequencing. These techniques have attracted
biologists as they offer powerful tools to explore microbial communities and
investigate their patterns of diversity in biological and biomedical samples at
remarkable resolution. However, the accuracy of these methods can negatively
affected by the presence of contamination. Several studies have recognized that
contamination is a common problem in microbial studies and have offered
promising computational and laboratory-based approaches to assess and remove
contaminants. Here we propose a novel strategy, MI-based (mutual information
based) filtering method, which uses information theoretic functionals and graph
theory to identify and remove contaminants. We applied MI-based filtering
method to a mock community data set and evaluated the amount of information
loss due to filtering taxa. We also compared our method to commonly practice
traditional filtering methods. In a mock community data set, MI-based filtering
approach maintained the true bacteria in the community without significant loss
of information. Our results indicate that MI-based filtering method effectively
identifies and removes contaminants in microbial communities and hence it can
be beneficial as a filtering method to microbiome studies. We believe our
filtering method has two advantages over traditional filtering methods. First,
it does not required an arbitrary choice of threshold and second, it is able to
detect true taxa with low abundance.


  The increasing quantity of multi-omics data, such as methylomic and
transcriptomic profiles, collected on the same specimen, or even on the same
cell, provide a unique opportunity to explore the complex interactions that
define cell phenotype and govern cellular responses to perturbations. We
propose a network approach based on Gaussian Graphical Models (GGMs) that
facilitates the joint analysis of paired omics data. This method, called DRAGON
(Determining Regulatory Associations using Graphical models on multi-Omic
Networks), calibrates its parameters to achieve an optimal trade-off between
the network's complexity and estimation accuracy, while explicitly accounting
for the characteristics of each of the assessed omics "layers." In simulation
studies, we show that DRAGON adapts to edge density and feature size
differences between omics layers, improving model inference and edge recovery
compared to state-of-the-art methods. We further demonstrate in an analysis of
joint transcriptome - methylome data from TCGA breast cancer specimens that
DRAGON can identify key molecular mechanisms such as gene regulation via
promoter methylation. In particular, we identify Transcription Factor AP-2 Beta
(TFAP2B) as a potential multi-omic biomarker for basal-type breast cancer.
DRAGON is available as open-source code in Python through the Network Zoo
package (netZooPy v0.8; netzoo.github.io).


  The anticancer potential of monoamine oxidase (MAO) was observed in
pre-clinical assays conducted with cell cultures and animals. L-Deprenyl (DEP)
causes apoptosis in melanoma, leukemia and mammary cells. High-dose DEP has
shown toxicity in mammary and pituitary cancers, as well as in monoblastic
leukemia, in rats. DEP accounts for immune-stimulant effect capable of
increasing natural killer cell activity, IL-2 generation, as well as of
inhibiting tumor growth. DEP administration in old female rats has increased
IL-2 generation and inverted the age-related depletion of IFN-{\gamma}
generation in the spleen. Co-adjuvant DEP administration helped
preventing/mitigating symptoms associated with peripheral neuropathy in cancer
treatment. It also enhanced the cytotoxic effects of antineoplastic drugs -
such as doxorubicin, cisplatin, among others - in cancer cells while they
protected healthy cells from being damaged. DEP presented effect against
dysfunctions such as debilitating hormone imbalance triggered by pituitary
gland tumor; this gland produces the stimulatory hormone of adrenocorticotropic
hormone which was related to the exacerbation of this disease. Thus, DEP
emerges as an excellent potential drug against several cancer types and it also
presents low toxicity in Parkinson`s disease patients subjected to long
treatment with it.


  This paper investigates the relationship between the spread of the COVID-19
pandemic, the state of community activity, and the financial index performance
across 20 countries. First, we analyze which countries behaved similarly in
2020 with respect to one of three multivariate time series: daily COVID-19
cases, Apple mobility data and national equity index price. Next, we study the
trajectories of all three of these attributes in conjunction to determine which
exhibited greater similarity. Finally, we investigate whether country financial
indices or mobility data responded quicker to surges in COVID-19 cases. Our
results indicate that mobility data and national financial indices exhibited
the most similarity in their trajectories, with financial indices responding
quicker. This suggests that financial market participants may have interpreted
and responded to COVID-19 data more efficiently than governments. Further,
results imply that efforts to study community mobility data as a leading
indicator for financial market performance during the pandemic were misguided.


  Despite major environmental and genetic differences, microbial metabolic
networks are known to generate consistent physiological outcomes across vastly
different organisms. This remarkable robustness suggests that, at least in
bacteria, metabolic activity may be guided by universal principles. The
constrained optimization of evolutionarily-motivated objective functions like
the growth rate has emerged as the key theoretical assumption for the study of
bacterial metabolism. While conceptually and practically useful in many
situations, the idea that certain functions are optimized is hard to validate
in data. Moreover, it is not always clear how optimality can be reconciled with
the high degree of single-cell variability observed in experiments within
microbial populations. To shed light on these issues, we develop an inverse
modeling framework that connects the fitness of a population of cells
(represented by the mean single-cell growth rate) to the underlying metabolic
variability through the Maximum-Entropy inference of the distribution of
metabolic phenotypes from data. While no clear objective function emerges, we
find that, as the medium gets richer, the fitness and inferred variability for
Escherichia coli populations follow and slowly approach the theoretically
optimal bound defined by minimal reduction of variability at given fitness.
These results suggest that bacterial metabolism may be crucially shaped by a
population-level trade-off between growth and heterogeneity.


  Epidemiological models with constant parameters may not capture satisfactory
infection patterns in the presence of pharmaceutical and non-pharmaceutical
mitigation measures during a pandemic, since infectiousness is a function of
time. In this paper, an Epidemiology-Informed Neural Network algorithm is
introduced to learn the time-varying transmission rate for the COVID-19
pandemic in the presence of various mitigation scenarios. There are
asymptomatic infectives, mostly unreported, and the proposed algorithm learns
the proportion of the total infective individuals that are asymptomatic
infectives. Using cumulative and daily reported cases of the symptomatic
infectives, we simulate the impact of non-pharmaceutical mitigation measures
such as early detection of infectives, contact tracing, and social distancing
on the basic reproduction number. We demonstrate the effectiveness of
vaccination on the transmission of COVID-19. The accuracy of the proposed
algorithm is demonstrated using error metrics in the data-driven simulation for
COVID-19 data of Italy, South Korea, the United Kingdom, and the United States.


  With the consolidation of deep learning in drug discovery, several novel
algorithms for learning molecular representations have been proposed. Despite
the interest of the community in developing new methods for learning molecular
embeddings and their theoretical benefits, comparing molecular embeddings with
each other and with traditional representations is not straightforward, which
in turn hinders the process of choosing a suitable representation for QSAR
modeling. A reason behind this issue is the difficulty of conducting a fair and
thorough comparison of the different existing embedding approaches, which
requires numerous experiments on various datasets and training scenarios. To
close this gap, we reviewed the literature on methods for molecular embeddings
and reproduced three unsupervised and two supervised molecular embedding
techniques recently proposed in the literature. We compared these five methods
concerning their performance in QSAR scenarios using different classification
and regression datasets. We also compared these representations to traditional
molecular representations, namely molecular descriptors and fingerprints. As
opposed to the expected outcome, our experimental setup consisting of over
25,000 trained models and statistical tests revealed that the predictive
performance using molecular embeddings did not significantly surpass that of
traditional representations. While supervised embeddings yielded competitive
results compared to those using traditional molecular representations,
unsupervised embeddings tended to perform worse than traditional
representations. Our results highlight the need for conducting a careful
comparison and analysis of the different embedding techniques prior to using
them in drug design tasks, and motivate a discussion about the potential of
molecular embeddings in computer-aided drug design.


  Feed-forward loops (FFLs) are among the most ubiquitously found motifs of
reaction networks in nature. However, little is known about their stochastic
behavior and the variety of network phenotypes they can exhibit. In this study,
we provide full characterizations of the properties of stochastic multimodality
of FFLs, and how switching between different network phenotypes are controlled.
We have computed the exact steady state probability landscapes of all eight
types of coherent and incoherent FFLs using the finite-butter ACME algorithm,
and quantified the exact topological features of their high-dimensional
probability landscapes using persistent homology. Through analysis of the
degree of multimodality for each of a set of 10,812 probability landscapes,
where each landscape resides over 10^5-10^6 microstates, we have constructed
comprehensive phase diagrams of all relevant behavior of FFL multimodality over
broad ranges of input and regulation intensities, as well as different regimes
of promoter binding dynamics. Our results show that with slow binding and
unbinding dynamics of transcription factor to promoter, FFLs exhibit strong
stochastic behavior that is very different from what would be inferred from
deterministic models. In addition, input intensity play major roles in the
phenotypes of FFLs: At weak input intensity, FFL exhibit monomodality, but
strong input intensity may result in up to 6 stable phenotypes. Furthermore, we
found that gene duplication can enlarge stable regions of specific
multimodalities and enrich the phenotypic diversity of FFL networks, providing
means for cells towards better adaptation to changing environment. Our results
are directly applicable to analysis of behavior of FFLs in biological processes
such as stem cell differentiation and for design of synthetic networks when
certain phenotypic behavior is desired.


  To optimize strategies for curbing the transmission of airborne pathogens,
the efficacy of three key controls -- face masks, ventilation, and physical
distancing -- must be well understood. In this study we used the
Quadrature-based model of Respiratory Aerosol and Droplets to quantify the
reduction in exposure to airborne pathogens from various combinations of
controls. For each combination of controls, we simulated thousands of scenarios
that represent the tremendous variability in factors governing airborne
transmission and the efficacy of mitigation strategies. While the efficacy of
any individual control was highly variable among scenarios, combining universal
mask-wearing with distancing of 1~m or more reduced the median exposure by more
than 99\% relative to a close, unmasked conversation, with further reductions
if ventilation is also enhanced. The large reductions in exposure to airborne
pathogens translated to large reductions in the risk of initial infection in a
new host. These findings suggest that layering controls is highly effective for
reducing transmission of airborne pathogens and will be critical for curbing
outbreaks of novel viruses in the future.


  Mid-infrared photothermal microscopy is a new chemical imaging technology in
which a visible beam senses the photothermal effect induced by a pulsed
infrared laser. This technology provides infrared spectroscopic information at
sub-micron spatial resolution and enables infrared spectroscopy and imaging of
living cells and organisms. Yet, current mid-infrared photothermal imaging
sensitivity suffers from a weak dependance of scattering on temperature and the
image quality is vulnerable to the speckles caused by scattering. Here, we
present a novel version of mid-infrared photothermal microscopy in which
thermo-sensitive fluorescent probes are harnessed to sense the mid-infrared
photothermal effect. The fluorescence intensity can be modulated at the level
of 1% per Kelvin, which is 100 times larger than the modulation of scattering
intensity. In addition, fluorescence emission is free of speckles, thus much
improving the image quality. Moreover, fluorophores can target specific
organelles or biomolecules, thus augmenting the specificity of photothermal
imaging. Spectral fidelity is confirmed through fingerprinting a single
bacterium. Finally, the photobleaching issue is successfully addressed through
the development of a wide-field fluorescence-enhanced mid-infrared photothermal
microscope which allows video rate bond-selective imaging of biological
specimens.


  A pandemic caused by a new coronavirus (COVID-19) has spread worldwide,
inducing an epidemic still active in Argentina. In this chapter, we present a
case study using an SEIR (Susceptible-Exposed-Infected-Recovered) diffusion
model of fractional order in time to analyze the evolution of the epidemic in
Buenos Aires and neighboring areas (Regi\'on Metropolitana de Buenos Aires,
(RMBA)) comprising about 15 million inhabitants. In the SEIR model, individuals
are divided into four classes, namely, susceptible (S), exposed (E), infected
(I) and recovered (R). The SEIR model of fractional order allows for the
incorporation of memory, with hereditary properties of the system, being a
generalization of the classic SEIR first-order system, where such effects are
ignored. Furthermore, the fractional model provides one additional parameter to
obtain a better fit of the data. The parameters of the model are calibrated by
using as data the number of casualties officially reported. Since infinite
solutions honour the data, we show a set of cases with different values of the
lockdown parameters, fatality rate, and incubation and infectious periods. The
different reproduction ratios R0 and infection fatality rates (IFR) so obtained
indicate the results may differ from recent reported values, constituting
possible alternative solutions. A comparison with results obtained with the
classic SEIR model is also included. The analysis allows us to study how
isolation and social distancing measures affect the time evolution of the
epidemic.


  Reputation is a powerful mechanism to enforce cooperation among unrelated
individuals through indirect reciprocity, but it suffers from disagreement
originating from private assessment, noise, and incomplete information. In this
work, we investigate stability of cooperation in the donation game by regarding
each player's reputation and behaviour as continuous variables. Through
perturbative calculation, we derive a condition that a social norm should
satisfy to give penalties to its close variants, provided that everyone
initially cooperates with a good reputation, and this result is supported by
numerical simulation. A crucial factor of the condition is whether a
well-reputed player's donation to an ill-reputed co-player is appreciated by
other members of the society, and the condition can be reduced to a threshold
for the benefit-cost ratio of cooperation which depends on the reputational
sensitivity to a donor's behaviour as well as on the behavioural sensitivity to
a recipient's reputation. Our continuum formulation suggests how indirect
reciprocity can work beyond the dichotomy between good and bad even in the
presence of inhomogeneity, noise, and incomplete information.


  Background: Creeping bentgrass (Agrostis soionifera) is a perennial grass of
Gramineae, belonging to cold season turfgrass, but has poor disease resistance.
Up to now, little is known about the induced systemic resistance (ISR)
mechanism, especially the relevant functional proteins, which is important to
disease resistance of turfgrass. Achieving more information of proteins of
infected creeping bentgrass is helpful to understand the ISR mechanism.
Results: With BDO treatment, creeping bentgrass seedlings were grown, and the
ISR response was induced by infecting Rhizoctonia solani. High-quality protein
sequences of creeping bentgrass seedlings were obtained. Some of protein
sequences were functionally annotated according to the database alignment while
a large part of the obtained protein sequences was left non-annotated. To treat
the non-annotated sequences, a prediction model based on convolutional neural
network was established with the dataset from Uniport database in three domains
to acquire good performance, especially the higher false positive control rate.
With established model, the non-annotated protein sequences of creeping
bentgrass were analyzed to annotate proteins relevant to disease-resistance
response and signal transduction. Conclusions: The prediction model based on
convolutional neural network was successfully applied to select good candidates
of the proteins with functions relevant to the ISR mechanism from the protein
sequences which cannot be annotated by database alignment. The waste of
sequence data can be avoided, and research time and labor will be saved in
further research of protein of creeping bentgrass by molecular biology
technology. It also provides reference for other sequence analysis of turfgrass
disease-resistance research.


  The purpose of this paper is to analyze the mechanism for the interplay of
deterministic and stochastic models for contagious diseases. Deterministic
models for contagious diseases are prone to predict global stability. Small
natural birth and death rates in comparison to disease parameters like the
contact rate and the removal rate ensures that the globally stable endemic
equilibrium corresponds to a tiny average proportion of infected individuals.
Asymptotic equilibrium levels corresponding to low numbers of individuals
invalidate the deterministic results.
  Diffusion effects force frequency functions of the stochastic model to
possess similar stability properties as the deterministic model. Particular
simulations of the stochastic model predict, however, oscillatory patterns.
Small and isolated populations show longer periods, more violent oscillations,
and larger probabilities of extinction.
  We prove that evolution maximizes the infectiousness of the disease as
measured by the ability to increase the proportion of infected individuals.
This holds provided the stochastic oscillations are moderate enough to keep the
proportion of susceptible individuals near a deterministic equilibrium.
  We close our paper with a discussion of the herd-immunity concept and stress
its close relation to vaccination-programs.


  Finding synthesis routes for molecules of interest is an essential step in
the discovery of new drugs and materials. To find such routes,
computer-assisted synthesis planning (CASP) methods are employed which rely on
a model of chemical reactivity. In this study, we model single-step
retrosynthesis in a template-based approach using modern Hopfield networks
(MHNs). We adapt MHNs to associate different modalities, reaction templates and
molecules, which allows the model to leverage structural information about
reaction templates. This approach significantly improves the performance of
template relevance prediction, especially for templates with few or zero
training examples. With inference speed several times faster than that of
baseline methods, we improve predictive performance for top-k exact match
accuracy for $\mathrm{k}\geq5$ in the retrosynthesis benchmark USPTO-50k.


  In this paper, we wish to investigate the dynamics of information transfer in
evolutionary dynamics. We use information theoretic tools to track how much
information an evolving population has obtained and managed to retain about
different environments that it is exposed to. By understanding the dynamics of
information gain and loss in a static environment, we predict how that same
evolutionary system would behave when the environment is fluctuating.
Specifically, we anticipate a cross-over between the regime in which
fluctuations improve the ability of the evolutionary system to capture
environmental information and the regime in which the fluctuations inhibit it,
governed by a cross-over in the timescales of information gain and decay.


  Systematic relations between multiple objects that occur in various fields
can be represented as networks. Real-world networks typically exhibit complex
topologies whose structural properties are key factors in characterizing and
further exploring the networks themselves. Uncertainty, modelling procedures
and measurement difficulties raise often insurmountable challenges in fully
characterizing most of the known real-world networks; hence, the necessity to
predict their unknown elements from the limited data currently available in
order to estimate possible future relations and/or to unveil unmeasurable
relations. In this work, we propose a deep learning approach to this problem
based on Graph Convolutional Networks for predicting networks while preserving
their original structural properties. The study reveals that this method can
preserve scale-free and small-world properties of complex networks when
predicting their unknown parts, a feature lacked by the up-to-date conventional
methods. An external validation realized by testing the approach on biological
networks confirms the results, initially obtained on artificial data. Moreover,
this process provides new insights into the retainability of network structure
properties in network prediction. We anticipate that our work could inspire
similar approaches in other research fields as well, where unknown mechanisms
behind complex systems need to be revealed by combining machine-based and
experiment-based methods.


  The dynamics of biofilm lifecycle are deeply influenced by the surrounding
environment and the interactions between sessile and planktonic phenotypes.
Bacterial biofilms typically develop in three distinct stages: attachment of
cells to a surface, growth of cells into colonies, and detachment of cells from
the colony into the surrounding medium. The attachment of planktonic cells
plays a prominent role in the initial phase of biofilm lifecycle as it
initiates the colony formation. During the maturation stage, biofilms harbor
numerous microenvironments which lead to metabolic heterogeneity. Such
microniches provide conditions suitable for the growth of new species, which
are present in the bulk liquid as planktonic cells and can penetrate the porous
biofilm matrix. We present a 1D continuum model on the interaction of sessile
and planktonic phenotypes in biofilm lifestyle which considers both the initial
attachment and colonization phenomena. The model is formulated as a
hyperbolic-elliptic free boundary value problem with vanishing initial value.
Hyperbolic equations reproduce the transport and growth of sessile species,
while elliptic equations model the diffusion and conversion of planktonic cells
and dissolved substrates. The attachment is modelled as a continuous,
deterministic process which depends on the concentrations of the attaching
species. The growth of new species is modelled through a reaction term in the
hyperbolic equations which depends on the concentration of planktonic species
within the biofilm. Existence and uniqueness of solutions are discussed and
proved for the attachment regime. Finally, some numerical examples show that
the proposed model correctly reproduces the growth of new species within the
biofilm and overcomes the ecological restrictions characterizing the
Wanner-Gujer type models.


  We aim at clarifying the controversy about the effects of school openings or
closures on the course of the Covid-19 pandemic. The mathematical analysis of
compartmental models with subpopulations shows that the in-school contact rates
affects the overall course of the pandemic only above a certain threshold that
separates an influence phase from a non-influence one. The threshold, that we
calculate via linear approximation in several cases, seems to appear in all
contexts, including outbreaks or new strains upsurge, lockdowns, and
vaccination campaigns excluding children, albeit with different values. Our
theoretical findings are then confirmed by several data-driven studies that
have previously identified the phase transition in specific cases. Specific
outcomes of this study are: opposite conclusions reached by studies of the same
or similar situations might depend on, possibly small, differences in modeling
or in parameter estimation from the very noisy Covid-19 data, that result in
identifying different phases; it is possible to keep schools open at any stage
of the Covid-19 pandemic, but suitably strict rules must be applied at all
times or else this becomes highly detrimental to virus containment efforts; as
the threshold during vaccination turns out to correspond to the internal
transmission rate that would lead to virus extinction if the school population
was isolated, the needed strict control can be sustained only for very brief
periods; as a result, either schools will have to face a prolonged closure or
children need to be vaccinated as well.


  Spatial constraints such as rigid barriers affect the dynamics of cell
populations, potentially altering the course of natural evolution. In this
paper, we investigate the population genetics of Escherichia coli proliferating
in microchannels with open ends. Our analysis is based on a population model in
which reproducing cells shift entire lanes of cells towards the open ends of
the channel. The model predicts that diversity is lost very rapidly within
lanes, but at a much slower pace among lanes. As a consequence, two mixed,
neutral E. coli strains competing in a microchannel must organize into an
ordered regular stripe pattern in the course of a few generations. These
predictions are in quantitative agreement with our experiments. We also
demonstrate that random mutations appearing in the middle of the channel are
much more likely to reach fixation than those occurring elsewhere. Our results
illustrate fundamental mechanisms of microbial evolution in spatially confined
space.


  Cell-fate transition can be modeled by ordinary differential equations (ODEs)
which describe the behavior of several molecules in interaction, and for which
each stable equilibrium corresponds to a possible phenotype (or 'biological
trait'). In this paper, we focus on simple ODE systems modeling two molecules
which each negatively (or positively) regulate the other. It is well-known that
such models may lead to monostability or multistability, depending on the
selected parameters. However, extensive numerical simulations have led systems
biologists to conjecture that in the vast majority of cases, there cannot be
more than two stable points. Our main result is a proof of this conjecture.
More specifically, we provide a criterion ensuring at most bistability, which
is indeed satisfied by most commonly used functions. This includes Hill
functions, but also a wide family of convex and sigmoid functions. We also
determine which parameters lead to monostability, and which lead to
bistability, by developing a more general framework encompassing all our
results.


  In late-2020, many countries around the world faced another surge in number
of confirmed cases of COVID-19, including United Kingdom, Canada, Brazil,
United States, etc., which resulted in a large nationwide and even worldwide
wave. While there have been indications that precaution fatigue could be a key
factor, no scientific evidence has been provided so far. We used a stochastic
metapopulation model with a hierarchical structure and fitted the model to the
positive cases in the US from the start of outbreak to the end of 2020. We
incorporated non-pharmaceutical interventions (NPIs) into this model by
assuming that the precaution strength grows with positive cases and studied two
types of pandemic fatigue. We found that people in most states and in the whole
US respond to the outbreak in a sublinear manner (with exponent k=0.5), while
only three states (Massachusetts, New York and New Jersey) have linear reaction
(k=1). Case fatigue (decline in people's vigilance to positive cases) is
responsible for 58% of cases, while precaution fatigue (decay of maximal
fraction of vigilant group) accounts for 26% cases. If there were no pandemic
fatigue (no case fatigue and no precaution fatigue), total positive cases would
have reduced by 68% on average. Our study shows that pandemic fatigue is the
major cause of the worsening situation of COVID-19 in United States. Reduced
vigilance is responsible for most positive cases, and higher mortality rate
tends to push local people to react to the outbreak faster and maintain
vigilant for longer time.


  Coronavirus disease 2019 (CoViD-19) is an infectious disease caused by severe
acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Among many symptoms,
cough, fever and tiredness are the most common. People over 60 years old and
with associated comorbidities are most likely to develop a worsening health
condition. This paper proposes a non-integer order model to describe the
dynamics of CoViD-19 in a standard population. The model incorporates the
reinfection rate in the individuals recovered from the disease. Numerical
simulations are performed for different values of the order of the fractional
derivative and of reinfection rate. The results are discussed from a biological
point of view.