  Mutual interference and prey refuge are important drivers of predator-prey
dynamics. The "exponent" or degree of mutual interference has been under much
debate in theoretical ecology. In the present work, we investigate the
interplay of the mutual interference exponent, on the behavior of a
predator-prey model with a generalized Holling type functional response. We
investigate stability properties of the system and derive conditions for the
occurrence of saddle-node and Hopf-bifurcations. A sufficient condition for
extinction of the prey species has also been derived for the model. In
addition, we investigate the effect of a prey refuge on the population dynamics
of the model and derive conditions for the prey refuge that would yield
persistence of populations. We provide additional verification our analytical
results via numerical simulations. Our findings are in accordance with
classical experimental results in ecology by Gauss G. F (1934), that show that
extinction of predator and prey populations is possible in a finite time period
- but that bringing in refuge can effectively cause persistence.


  The Protein Data Bank (PDB) today contains more than 153,000 entries with the
3-dimensional structures of biological macromolecules. Using the rich resources
of this repository, it is possible identifying subsets with specific,
interesting properties for different applications. Our research group prepared
an automatically updated list of amyloid- and probably amyloidogenic molecules,
the PDB\_Amyloid collection, which is freely available at the address
\url{http://pitgroup.org/amyloid}. This resource applies exclusively the
geometric properties of the steric structures for identifying amyloids. In the
present contribution, we analyze the starting (i.e., prefix) subsequences of
the characteristic, parallel beta-sheets of the structures in the PDB\_Amyloid
collection, and identify further appearances of these length-5 prefix
subsequences in the whole PDB data set. We have identified this way numerous
proteins, whose normal or irregular functions involve amyloid formation,
structural misfolding, or anti-coagulant properties, simply by containing these
prefixes: including the T-cell receptor (TCR), bound with the major
histocompatibility complexes MHC-1 and MHC-2; the p53 tumor suppressor protein;
a mycobacterial RNA polymerase transcription initialization complex; the human
bridging integrator protein BIN-1; and the tick anti-coagulant peptide TAP.


  In this paper, we propose a dynamical model to describe the transmission of
COVID-19, which is spreading in China and many other countries. To avoid a
larger outbreak in the worldwide, Chinese government carried out a series of
strong strategies to prevent the situation from deteriorating. Home quarantine
is the most important one to prevent the spread of COVID-19. In order to
estimate the effect of population quarantine, we divide the population into
seven categories for simulation. Based on a Least-Squares procedure and
officially published data, the estimation of parameters for the proposed model
is given. Numerical simulations show that the proposed model can describe the
transmission of COVID-19 accurately, the corresponding prediction of the trend
of the disease is given. The home quarantine strategy plays an important role
in controlling the disease spread and speeding up the decline of COVID-19. The
control reproduction number of most provinces in China are analyzed and
discussed adequately. We should pay attention to that, though the epidemic is
in decline in China, the disease still has high risk of human-to-human
transmission continuously. Once the control strategy is removed, COVID-19 may
become a normal epidemic disease just like flu. Further control for the disease
is still necessary, we focus on the relationship between the spread rate of the
virus and the meteorological conditions. A comprehensive meteorological index
is introduced to represent the impact of meteorological factors on both high
and low migration groups. As the progress on the new vaccine, we design detail
vaccination strategies for COVID-19 in different control phases and show the
effectiveness of efficient vaccination. Once the vaccine comes into use, the
numerical simulation provide a promptly prospective research.


  The evolution and emergence of antibiotic resistance is a major public health
concern. The understanding of the within-host microbial dynamics combining
mutational processes, horizontal gene transfer and resource consumption, is one
of the keys to solve this problem. We analyze a generic model to rigorously
describe interactions dynamics of four bacterial strains: one fully sensitive
to the drug, one with mutational resistance only, one with plasmidic resistance
only and one with both resistances. By dening thresholds numbers (i.e. each
strain's eective reproduction and each strain's transition thresholds numbers),
we rst express conditions for the existence of non trivial stationary states.
We nd that these thresholds mainly depend on bacteria quantitative traits such
as nutrient consumption ability, growth conversion factor, death rate, mutation
(forward or reverse) and segregational loss of plasmid probabilities (for
plasmid-bearing strains). Next, with respect to the order in the set of
strain's eective reproduction thresholds numbers, we show that the qualitative
dynamics of the model range from the extinction of all strains, coexistence of
sensitive and mutational resistance strains to the coexistence of all strains
at equilibrium. Finally, we go through some applications of our general
analysis depending on whether bacteria strains interact without or with drug
action (either cytostatic or cytotoxic).


  By-products represent a major disposal problem for the food industry, but
they are also promising sources of bioactive compounds. Olive pomace, one of
the main by-products of olive oil production, is a potential low-cost,
phenol-rich ingredient for the formulation of functional food. In this study,
bakery products enriched with defatted olive pomace powder and their
conventional counterparts were chemical characterized and in vitro digested.
The bioaccessible fractions were supplemented to cultured human intestinal
cells exposed to an inflammatory stimulus, and the anti-inflammatory effect and
metabolome modification were evaluated. Although in all bakery products the
enrichment with olive pomace significantly increased the total phenolic
content, this increase was paralleled by an enhanced anti-inflammatory activity
only in conventionally fermented bread. Therefore, while confirming olive oil
by-products as functional ingredients for bakery food enrichment, our data
highlight that changes in chemical composition cannot predict changes in
functionality. Functionality should be first evaluated in biological in vitro
systems, and then confirmed in human intervention studies.


  Bacteria populate the colon where they replicate and migrate in response to
nutrient availability. Here I model the colon bacterial population as a
sandpile model, the colon-pile. Sand addition mimics bacterial replication and
grains toppling represents bacterial migration coupled to high population
density. The numerical simulations reveal a behaviour similar to
non-conservative sandpile models, approaching a critical state with system wide
avalanches when the death rate becomes negligible. The critical exponents
estimation indicates that the colon-pile belongs to a new universality class.
This work suggest that the colon microbiome is in a self-organised critical
state, where small perturbations can trigger large scale rearrangements,
covering an area comparable to the system size and characterised by a 1/f noise
spectra


  Pepducins are cell-penetrating, membrane-tethered lipopeptides designed to
target the intracellular region of a G protein-coupled receptor (GPCR) in order
to allosterically modulate the receptor's signaling output. In this
proof-of-concept study, we explored the pain-relief potential of a pepducin
series derived from the first intracellular loop of neurotensin receptor type 1
(NTS1), a class A GPCR that mediates many of the effects of the neurotensin
(NT) tridecapeptide, including hypothermia, hypotension and analgesia. We used
BRET-based biosensors to determine the pepducins' ability to engage G protein
signaling pathways associated with NTS1 activation. We observed partial Gq and
G13 activation at a 10 {\mu}M concentration, indicating that these pepducins
may act as allosteric agonists of NTS1. Additionally, we used surface plasmon
resonance (SPR) as a label-free assay to monitor pepducin-induced responses in
CHO-K1 cells stably expressing hNTS1. This whole-cell integrated assay enabled
us to subdivide our pepducin series into three profile response groups. In
order to determine the pepducins' antinociceptive potential, we then screened
the series in an acute pain model (tail-flick test) by measuring tail
withdrawal latencies to a thermal nociceptive stimulus, following intrathecal
pepducin administration (275 nmol/kg). We further evaluated promising pepducins
in a tonic pain model (formalin test), as well as in neuropathic (Chronic
Constriction Injury) and inflammatory (Complete Freund's Adjuvant) chronic pain
models. We report one pepducin, PP-001, that consistently reduced rat
nociceptive behaviors, even in chronic pain paradigm. Altogether, these results
suggest that NTS1-derived pepducins may represent a promising strategy in
pain-relief.


  We consider a stochastic individual-based model for the evolution of a
haploid, asexually reproducing population. The space of possible traits is
given by the vertices of a (possibly directed) finite graph $G=(V,E)$. The
evolution of the population is driven by births, deaths, competition, and
mutations along the edges of $G$. We are interested in the large population
limit under a mutation rate $\mu_K$ given by a negative power of the carrying
capacity $K$ of the system: $\mu_K=K^{-1/\alpha},\alpha>0$. This results in
several mutant traits being present at the same time and competing for invading
the resident population. We describe the time evolution of the orders of
magnitude of each sub-population on the $\log K$ time scale, as $K$ tends to
infinity. Using techniques developed in [Champagnat, M\'el\'eard, Tran, 2019]
we show that these are piecewise affine continuous functions, whose slopes are
given by an algorithm describing the changes in the fitness landscape due to
the succession of new resident or emergent types. This work generalises [Kraut,
Bovier, 2019] to the stochastic setting, and Theorem 3.2 of [Bovier, Coquille,
Smadi, 2018] to any finite mutation graph. We illustrate our theorem by a
series of examples describing surprising phenomena arising from the geometry of
the graph and/or the rate of mutations.


  Variational Autoencoders (VAEs) provide a flexible and scalable framework for
non-linear dimensionality reduction. However, in application domains such as
genomics where data sets are typically tabular and high-dimensional, a
black-box approach to dimensionality reduction does not provide sufficient
insights. Common data analysis workflows additionally use clustering techniques
to identify groups of similar features. This usually leads to a two-stage
process, however, it would be desirable to construct a joint modelling
framework for simultaneous dimensionality reduction and clustering of features.
In this paper, we propose to achieve this through the BasisVAE: a combination
of the VAE and a probabilistic clustering prior, which lets us learn a one-hot
basis function representation as part of the decoder network. Furthermore, for
scenarios where not all features are aligned, we develop an extension to handle
translation-invariant basis functions. We show how a collapsed variational
inference scheme leads to scalable and efficient inference for BasisVAE,
demonstrated on various toy examples as well as on single-cell gene expression
data.


  The recognition of essential proteins not only can help to understand the
mechanism of cell operation, but also help to study the mechanism of biological
evolution. At present, many scholars have been discovering essential proteins
according to the topological structure of protein network and complexes. While
some proteins still can not be recognized. In this paper, we proposed two new
methods complex degree centrality (CDC) and complex in-degree and betweenness
definition (CIBD) which integrate the local character of protein complexes and
topological properties to determine the essentiality of proteins. First, we
give the definitions of complex average centrality (CAC) and complex hybrid
centrality (CHC) which both describe the properties of protein complexes. Then
we propose these new methods CDC and CIBD based on CAC and CHC definitions. In
order to access these two methods, different Protein-Protein Interaction (PPI)
networks of Saccharomyces cerevisiae, DIP, MIPS and YMBD are used as
experimental materials. Experimental results in networks show that the methods
of CDC and CIBD can help to improve the precision of predicting essential
proteins.


  We comparatively analyzed the spatiotemporal fluctuations of the 2019-novel
coronavirus (COVID-19) and SARS outbreaks to understand their epidemiological
characteristics. Methodologically, we introduced TPL (Taylor power law) to
characterize their spatiotemporal heterogeneity/stability and Hubbell (2001)
unified neutral theory of biodiversity (UNTB) [specifically Harris et al.
(2015) HDP-MSN model (hierarchical Dirichlet process multi-site neutral model)]
to approximate the metapopulation of coronavirus infections. First, TPL
analysis suggested that the coronaviruses appear to have a specific
heterogeneity/stability scaling parameter (TPL-b) slightly exceeding 2 for
cumulative infections or exceeding 1 for daily incremental infections,
suggesting their potentially chaotic, unstable outbreaks. Another TPL parameter
(M0) (i.e., infection critical threshold) depends on virus kinds
(COVID-19/SARS), time (disease-stages), space (regions) and public-health
interventions (e.g., quarantines and mobility control). M0 measures the
infection level, at which infections are random (Poisson distribution) and
below which infections follow uniform distribution and may die off if M0
coincides or below the level of Allee effects. It was found that COVID-19
outbreak seems nearly twice more risky than SARS, and the lower infection
threshold may be due to its lower lethality than SARS since lower fatality
rates can facilitate the survival and spread of pathogen. Second, metacommunity
UNTB neutrality testing seems appropriate for approximating metapopulation of
coronavirus infections. Specifically, two parameters {\theta} and M, borrowed
from neutral theory, may be used to assess the relative significance of
infection through local contagion vs. infection through migration, both of
which may depend on time, space, virus kinds, and particularly public-health
interventions.


  In a recent report we have collected some data about the COVID-19 epidemics
in Northern Italy; in this follow-up we analyze how these changed after the
mild restrictive measures taken by the Government two weeks ago and the large
campaign of public awareness developed in the meanwhile.


  In his seminal work in the 1970s Robert May suggested that there was an upper
limit to the number of species that could be sustained in stable equilibrium by
an ecosystem. This deduction was at odds with both intuition and the observed
complexity of many natural ecosystems. The so-called stability-diversity debate
ensued, and the discussion about the factors making an ecosystem stable or
unstable continues to this day. We show in this work that dispersal can be a
destabilising influence. To do this, we combine ideas from Alan Turing's work
on pattern formation with May's random-matrix approach. We demonstrate how a
stable equilibrium in a complex ecosystem with two trophic levels can become
unstable with the introduction of dispersal in space. Conversely, we show that
Turing instabilities can occur more easily in complex ecosystems with many
species than in the case of only a few species. Our work shows that adding more
details to the model of May gives rise to more ways in which an equilibrium can
become unstable. Making May's simple model more realistic is therefore unlikely
to remove the upper bound on complexity.


  Mapping pathways to achieving the sustainable development goals requires
understanding and predicting how social, economic and political factors impact
biodiversity. Trends in demography, economic growth, regional alliances and
consumption behaviours can have profound effects on the environment by driving
resource use and production. While these distant socio-economic drivers impact
species and ecosystems at global scales, for example by driving greenhouse gas
emissions and climate change, the most prevalent human impacts on biodiversity
manifest through habitat loss and land use change decisions at finer scales. We
provide the first integrated ecological-economic analysis pathway capable of
supporting both national policy design challenges and global scale assessment
of biodiversity risks posed by socio-economic drivers such as population
growth, consumption and trade. To achieve this, we provide state-of-the-art
integration of economic, land use, and biodiversity modelling, and illustrate
its application using two case studies. We evaluate the national-level
implications of change in trading conditions under a multi-lateral free trade
agreement for the bird biodiversity of Vietnam. We review the implications for
land-use and biodiversity under coupled socio-economic (Shared Socioeconomic
Pathways) and climate (Resource Concentration Pathways) scenarios for
Australia. Our study provides a roadmap for setting up high dimensional
integrated analyses foe evaluating global priorities for protecting nature and
livelihoods in vulnerable areas with the greatest conflicts for economic,
social and environmental opportunities.


  Networks effectively capture interactions among components of complex
systems, and have thus become a mainstay in many scientific disciplines.
Growing evidence, especially from biology, suggest that networks undergo
changes over time, and in response to external stimuli. In biology and
medicine, these changes have been found to be predictive of complex diseases.
They have also been used to gain insight into mechanisms of disease initiation
and progression. Primarily motivated by biological applications, this article
provides a review of recent statistical machine learning methods for inferring
networks and identifying changes in their structures.


  We analyze ecological systems that are influenced by random environmental
fluctuations. We first provide general conditions which ensure that the species
coexist and the system converges to a unique invariant probability measure
(stationary distribution). Since it is usually impossible to characterize this
invariant probability measure analytically, we develop a powerful method for
numerically approximating invariant probability measures. This allows us to
shed light upon how the various parameters of the ecosystem impact the
stationary distribution.
  We analyze different types of environmental fluctuations. At first we study
ecosystems modeled by stochastic differential equations. In the second setting
we look at piecewise deterministic Markov processes. These are processes where
one follows a system of differential equations for a random time, after which
the environmental state changes, and one follows a different set of
differential equations -- this procedure then gets repeated indefinitely.
Finally, we look at stochastic differential equations with switching, which
take into account both the white noise fluctuations and the random
environmental switches.
  As applications of our theoretical and numerical analysis, we look at
competitive Lotka--Volterra, Beddington-DeAngelis predator-prey, and
rock-paper-scissors dynamics. We highlight new biological insights by analyzing
the stationary distributions of the ecosystems and by seeing how various types
of environmental fluctuations influence the long term fate of populations.


  Coronavirus (COVID-19) outbreak in late 2019 and 2020 comprises a serious and
more likely a pandemic threat worldwide. Given that the disease has not
approved vaccines or drugs up to now, any efforts for drug design and or
clinical trails of old drugs based on their mechanism of action are worthy and
creditable in such circumstances. Experienced docking experiments using the
newly released coordinate structure for COVID-19 protease as a receptor and
thoughtfully selected chemicals among antiviral and antibiotics drugs as
ligands may be leading in this context. We selected nine drugs from HIV-1
protease inhibitors and twenty-one candidates from anti bronchitis drugs based
on their chemical structures and enrolled them in blind and active
site-directed dockings in different modes and in native-like conditions of
interactions. Our findings suggest the binding capacity and the inhibitory
potency of candidates are as follows
Tipranavir>Indinavir>Atazanavir>Darunavir>Ritonavir>Amprenavir for HIV-1
protease inhibitors and Cefditoren>Cefixime>Erythromycin>Clarithromycin for
anti bronchitis medicines. The drugs bioavailability, their hydrophobicity and
the hydrophobic properties of their binding sites and also the rates of their
metabolisms and deactivations in the human body are the next determinants for
their overall effects on viral infections, the net results that should survey
by clinical trials to assess their therapeutic usefulness for coronavirus
infections.


  A library of substituted pyrimidines was synthesized and evaluated for free
radical scavenging, and in vitro cytotoxic activity in 3T3 cells. All compounds
showed good free radical scavenging activity with IC50 values in the range of
42.9 + 0.31 to 438.3 3.3 {\mu}M as compared to the standard butylated
hydroxytoluene having IC50 value of 128.83 2. 1 {\mu}M. The structure
activity-relationship was also established. Selected analogues 1, 2, 3, 5, 6,
7, 8, 9, 10, 12, 13, 15, 19, 20, 21, 24, 25, 26 and 28 were tested for
cytotoxicity in mouse fibroblast 3T3 cell line using MTT assay, and most of the
analogues showed cytotoxicity. This study has identified a number of cytotoxic
novel substituted pyrimidines having free radical scavenging activities that
can be used as inhibitory compounds for those cancer cells whose growth is
mediated by reactive oxygen species.


  Metastatic Melanoma, the fifth most common cancer in the western countries
and the most common malignancy diagnosed in United States present itself as the
most lethal treatment resistant cancer worldwide. In addition to the reactive
oxygen species(ROS), mutations in the genes encoding receptors and non-receptor
tyrosine/serene/threonine protein kinases are known to be involved in its
etiology. Kinases are molecular players of cell survival, growth, and
proliferation and migration that mediate their effects via various signal
transduction pathways. A number of such molecular players have been previously
found to be mutated and hyper phosphorylated in melanoma. Although, several
systemic therapies including cytotoxic chemotherapy, targeted drugs, hormonal
therapy, radiation therapy, bio-chemotherapy, and therapies that inhibit
negative regulation of immune system have been approved from U. S. Food and
Drug Administration (FDA) for metastatic melanoma treatment. However, no
systemic therapy has meaningfully changed its survival end points so far and
surgery still presents primary treatment option for advanced and metastatic
melanomaa due to its highly resistant nature towards systemic drugs, high rate
of severe, life-threatening, or fatal side effects, and un satisfactory overall
response rate. Therefore, there is still a need to develop therapies that
target the unique molecular profile of melanoma tumors.


  This review mainly focuses on the relation between antioxidants with cancer
therapy. Antioxidants have been reported to play an essential role to reduce
free radical species. Free radicals commonly cause oxidative damage which is a
common factor in the aging process, and also the vital factor of formation, and
development of major disease specially cancer. Although, since last many
decades several antioxidants belong to natural and synthetic origin have been
tested in clinical trials against oxidative stress, however these clinical
trials end up with undesirable effects. This review also complied with the most
recent findings of oxidative stress, highlighting of free racial production,
and its related oxidative damage at cellular and molecular level, with the new
and existing natural and synthetic classes of free radical scavenger and their
related clinical trials.


  Delay differential equations are used as a model when the effect of past
states has to be taken into account. In this work we consider delay models of
chemical reaction networks with mass action kinetics. We obtain a sufficient
condition for absolute delay stability of equilibrium concentrations, i.e.,
local asymptotic stability independent of the delay parameters. Several
interesting examples on sequestration networks with delays are presented.


  In this article, we introduce the R package EpiILM, which provides tools for
simulation from, and inference for, discrete-time individual-level models of
infectious disease transmission proposed by Deardon et al. (2010). The
inference is set in a Bayesian framework and is carried out via
Metropolis-Hastings Markov chain Monte Carlo (MCMC). For its fast
implementation, key functions are coded in Fortran. Both spatial and contact
network models are implemented in the package and can be set in either
susceptible-infected (SI) or susceptible-infected-removed (SIR) compartmental
frameworks. The use of the package is demonstrated through examples involving
both simulated and real data.


  Objectives: We aim to assess the impact of temperature and relative humidity
on the transmission of COVID-19 across communities after accounting for
community-level factors such as demographics, socioeconomic status, and human
mobility status. Design: A retrospective cross-sectional regression analysis
via the Fama-MacBeth procedure is adopted. Setting: We use the data for
COVID-19 daily symptom-onset cases for 100 Chinese cities and COVID-19 daily
confirmed cases for 1,005 U.S. counties. Participants: A total of 69,498 cases
in China and 740,843 cases in the U.S. are used for calculating the effective
reproductive numbers. Primary outcome measures: Regression analysis of the
impact of temperature and relative humidity on the effective reproductive
number (R value). Results: Statistically significant negative correlations are
found between temperature/relative humidity and the effective reproductive
number (R value) in both China and the U.S. Conclusions: Higher temperature and
higher relative humidity potentially suppress the transmission of COVID-19.
Specifically, an increase in temperature by 1 degree Celsius is associated with
a reduction in the R value of COVID-19 by 0.026 (95% CI [-0.0395,-0.0125]) in
China and by 0.020 (95% CI [-0.0311, -0.0096]) in the U.S.; an increase in
relative humidity by 1% is associated with a reduction in the R value by 0.0076
(95% CI [-0.0108,-0.0045]) in China and by 0.0080 (95% CI [-0.0150,-0.0010]) in
the U.S. Therefore, the potential impact of temperature/relative humidity on
the effective reproductive number alone is not strong enough to stop the
pandemic.


  Essential protein plays a crucial role in the process of cell life. The
identification of essential proteins can not only promote the development of
drug target technology, but also contribute to the mechanism of biological
evolution. There are plenty of scholars who pay attention to discovering
essential proteins according to the topological structure of protein network
and biological information. The accuracy of protein recognition still demands
to be improved. In this paper, we propose a method which integrate the
clustering coefficient in protein complexes and topological properties to
determine the essentiality of proteins. First, we give the definition of
In-clustering coefficient (IC) to describe the properties of protein complexes.
Then we propose a new method, complex edge and node clustering coefficient
(CENC) to identify essential proteins. Different Protein-Protein Interaction
(PPI) networks of Saccharomyces cerevisiae, MIPS and DIP are used as
experimental materials. Through some experiments of logistic regression model,
the results show that the method of CENC can promote the ability of recognizing
essential proteins, by comparing with the existing methods DC, BC, EC, SC, LAC,
NC and the recent method UC.


  Coronaviruses are a class of virus responsible of the recent outbreak of
Human Severe Acute Respiratory Syndrome. The molecular machinery behind the
viral entry and thus infectivity is based on the formation of the complex of
virus spike protein with the angiotensin-converting enzyme 2 (ACE2). The
detection of putative allosteric sites on the viral spike protein can trace the
path to develop allosteric drugs to weaken the strength of the spike-ACE2
interface and, thus, reduce the viral infectivity. In this work we present
results of the application of the Protein Contact Network (PCN) paradigm to the
complex SARS-CoV spike - ACE2 relative to both 2003 SARS and the recent 2019 -
CoV. Results point to a specific region, present in both structures, that is
predicted to act as allosteric site modulating the binding of the spike protein
with ACE2.


  1. Joint Species Distribution models (JSDMs) explain spatial variation in
community composition by contributions of the environment, biotic associations,
and possibly spatially structured residual covariance. They show great promise
as a general analytical framework for community ecology and macroecology, but
current JSDMs, even when approximated by latent variables, scale poorly on
large datasets, limiting their usefulness for currently emerging big (e.g.,
metabarcoding and metagenomics) community datasets. 2. Here, we present a
novel, more scalable JSDM (sjSDM) that circumvents the need to use latent
variables by using a Monte-Carlo integration of the joint JSDM likelihood and
allows flexible elastic net regularization on all model components. We
implemented sjSDM in PyTorch, a modern machine learning framework that can make
use of CPU and GPU calculations. Using simulated communities with known
species-species associations and different number of species and sites, we
compare sjSDM with state-of-the-art JSDM implementations to determine
computational runtimes and accuracy of the inferred species-species and
species-environmental associations. 3. We find that sjSDM is orders of
magnitude faster than existing JSDM algorithms (even when run on the CPU) and
can be scaled to very large datasets. Despite the dramatically improved speed,
sjSDM produces more accurate estimates of species association structures than
alternative JSDM implementations. We demonstrate the applicability of sjSDM to
big community data using eDNA case study with thousands of fungi operational
taxonomic units (OTU). 4. Our sjSDM approach makes the analysis of JSDMs to
large community datasets with hundreds or thousands of species possible,
substantially extending the applicability of JSDMs in ecology. We provide our
method in an R package to facilitate its applicability for practical data
analysis.


  The cell cycle duration is a variable cellular phenotype that underlies
long-term population growth and age structures. By analyzing the stationary
solutions of a branching process with heritable cell division times, we
demonstrate existence of a phase transition, which can be continuous or
first-order, by which a non-zero fraction of the population becomes localized
at a minimal division time. Just below the transition, we demonstrate
coexistence of localized and delocalized age-structure phases, and power law
decay of correlation functions. Above it, we observe self-synchronization of
cell cycles, collective divisions, and slow 'aging' of population growth rates.


  Current protein forcefields like the ones seen in CHARMM or Xplor-NIH have
many terms that include bonded and non-bonded terms. Yet the forcefields do not
take into account the use of hydrogen bonds which are important for secondary
structure creation and stabilization of proteins. SCOPE is an open-source
program that generates proteins from rotamer space. It then creates a
forcefield that uses only non-bonded and hydrogen bond energy terms to create a
profile for a given protein. The profiles can then be used in an artificial
neural network to create a linear model that is funneled to the true protein
conformation.


  In December 2019, a novel coronavirus was found in a seafood wholesale market
in Wuhan, China. WHO officially named this coronavirus as COVID-19. Since the
first patient was hospitalized on December 12, 2019, China has reported a total
of 78,824 confirmed CONID-19 cases and 2,788 deaths as of February 28, 2020.
Wuhan's cumulative confirmed cases and deaths accounted for 61.1% and 76.5% of
the whole China mainland , making it the priority center for epidemic
prevention and control. Meanwhile, 51 countries and regions outside China have
reported 4,879 confirmed cases and 79 deaths as of February 28, 2020. COVID-19
epidemic does great harm to people's daily life and country's economic
development. This paper adopts three kinds of mathematical models, i.e.,
Logistic model, Bertalanffy model and Gompertz model. The epidemic trends of
SARS were first fitted and analyzed in order to prove the validity of the
existing mathematical models. The results were then used to fit and analyze the
situation of COVID-19. The prediction results of three different mathematical
models are different for different parameters and in different regions. In
general, the fitting effect of Logistic model may be the best among the three
models studied in this paper, while the fitting effect of Gompertz model may be
better than Bertalanffy model. According to the current trend, based on the
three models, the total number of people expected to be infected is 49852-57447
in Wuhan,12972-13405 in non-Hubei areas and 80261-85140 in China respectively.
The total death toll is 2502-5108 in Wuhan, 107-125 in Non-Hubei areas and
3150-6286 in China respetively. COVID-19 will be over p robably in late-April,
2020 in Wuhan and before late-March, 2020 in other areas respectively.


  Broom and Rycht\'{a}\v{r} [Proc. R. Soc. A (2008) 464, 2609--2627] found an
exact solution for the fixation probabilities of the Moran process for a
structured population, in which the interaction structure among individuals is
given by the so-called star graph, i.e. one central vertex and $n$ leaves, the
leaves connecting only to the center. We generalize on their solution by
allowing individuals' fitnesses to depend on the population frequency, and also
by allowing a possible change in the order of reproduction and death draws. In
their cited paper, Broom and Rycht\'{a}\v{r} considered the birth-death (BD)
process, in which at each time step an individual is first drawn for
reproduction and then an individual is selected for death. In the death-birth
(DB) process, the order of the draws is reversed. It may be seen that the order
of the draws makes a big difference in the fixation probabilities. Our solution
method applies to both the BD and the DB cases. As expected, the exact formulae
for the fixation probabilities are complicated. We will also illustrate them
with some examples and provide results on the asymptotic behavior of the
fixation probabilities when the number $n$ of leaves in the graph tends to
infinity.


  The first time a searcher finds a target is called a first passage time
(FPT). In many physical, chemical, and biological processes, the searcher is
"mortal," which means that the searcher might become inactivated (degrade, die,
etc.) before finding the target. In the context of intracellular signaling, an
important recent work discovered that fast inactivation can drastically alter
the conditional FPT distribution of a mortal diffusive searcher, if the
searcher is conditioned to find the target before inactivation. In this paper,
we prove a general theorem which yields an explicit formula for all the moments
of such conditional FPTs in the fast inactivation limit. This formula is quite
universal, as it holds under very general conditions on the diffusive searcher
dynamics, the target, and the spatial domain. These results prove in
significant generality that if inactivation is fast, then the conditional FPT
compared to the FPT without inactivation is (i) much faster, (ii) much less
affected by spatial heterogeneity, and (iii) much less variable. Our results
agree with recent computational and theoretical analysis of a certain discrete
intracellular diffusion model and confirm a conjecture related to the effect of
spatial heterogeneity on intracellular signaling.


  Today, we are all threatened by an unprecedented pandemic: COVID-19. How
different is it from other coronaviruses? Will it be attenuated or become more
virulent? Which animals may be its original host? In this study, we collected
and analyzed nearly thirty thousand publicly available complete genome
sequences for COVID-19 virus from 79 different countries, the previously known
flu-causing coronaviruses (HCov-229E, HCov-OC43, HCov-NL63 and HCov-HKU1) and
the lethal, pathogenic viruses, SARS, MERS, Victoria, Lassa, Yamagata, Ebola,
and Dengue. We found strong similarities between the current circulating
COVID-19 and SARS and MERS, as well as COVID-19 in rhinolophines and pangolins.
On the contrary, COVID-19 shares little similarity with the flu-causing
coronaviruses and the other known viruses. Strikingly, we observed that the
divergence of COVID-19 strains isolated from human hosts has steadily increased
from December 2019 to May 2020, suggesting COVID-19 is actively evolving in
human hosts. In this paper, we first propose a novel MLCS algorithm NP-MLCS1
for the big sequence analysis, which can calculate the common model for
COVID-19 complete genome sequences to provide important information for vaccine
and antibody development. Geographic and time-course analysis of the evolution
trees of the human COVID-19 reveals possible evolutional paths among strains
from 79 countries. This finding has important implications to the management of
COVID-19 and the development of vaccines and medications.


  In this paper, based on the Akaike information criterion, root mean square
error and robustness coefficient, a rational evaluation of various epidemic
models/methods, including seven empirical functions, four statistical inference
methods and five dynamical models, on their forecasting abilities is carried
out. With respect to the outbreak data of COVID-19 epidemics in China, we find
that before the inflection point, all models fail to make a reliable
prediction. The Logistic function consistently underestimates the final
epidemic size, while the Gompertz's function makes an overestimation in all
cases. Towards statistical inference methods, the methods of sequential
Bayesian and time-dependent reproduction number are more accurate at the late
stage of an epidemic. And the transition-like behavior of exponential growth
method from underestimation to overestimation with respect to the inflection
point might be useful for constructing a more reliable forecast. Compared to
ODE-based SIR, SEIR and SEIR-AHQ models, the SEIR-QD and SEIR-PO models
generally show a better performance on studying the COVID-19 epidemics, whose
success we believe could be attributed to a proper trade-off between model
complexity and fitting accuracy. Our findings not only are crucial for the
forecast of COVID-19 epidemics, but also may apply to other infectious
diseases.


  Started in Wuhan, China, the COVID-19 has been spreading all over the world.
We calibrate the logistic growth model, the generalized logistic growth model,
the generalized Richards model and the generalized growth model to the reported
number of infected cases for the whole of China, 29 provinces in China, and 33
countries and regions that have been or are undergoing major outbreaks. We
dissect the development of the epidemics in China and the impact of the drastic
control measures both at the aggregate level and within each province. We
quantitatively document four phases of the outbreak in China with a detailed
analysis on the heterogeneous situations across provinces. The extreme
containment measures implemented by China were very effective with some
instructive variations across provinces. Borrowing from the experience of
China, we made scenario projections on the development of the outbreak in other
countries. We identified that outbreaks in 14 countries (mostly in western
Europe) have ended, while resurgences of cases have been identified in several
among them. The modeling results clearly show longer after-peak trajectories in
western countries, in contrast to most provinces in China where the after-peak
trajectory is characterized by a much faster decay. We identified three groups
of countries in different level of outbreak progress, and provide informative
implications for the current global pandemic.


  The outbreaks of algal blooms occur in both shallow and deep-water bodies. To
compare the characteristics of algal blooms in the shallow and deep water, we
consider the Arabian Gulf and Sea of Oman as a case study. While the Arabian
Gulf is a shallow region dominated by advective features, the Sea of Oman is a
deep channel dominated by numerous eddies. Both these regions is have rarely
been studied due to lack of available data thus preventing the characterization
of phenomenon such as algal blooms in the regions. Nonetheless, a recent unique
and comprehensive dataset of the frequent algal blooms collected over the last
decade in the Arabian Gulf and Oman Sea is utilized in this study to analyze
the spatio-temporal variability of the blooms thereof. These data are also used
to characterize the algal bloom species and analyze the relationship between
the water properties and the algal blooms in the shallow and deep waters. There
is a general decreasing trend of the algal bloom events from 2010 to 2018 in
the Arabian Gulf while in the Sea of Oman there is an increasing trend. We
reveal a clear seasonality with the highest frequency of algal blooms during
winter and spring. We have noticed that algal blooms have the feature of the
annual cycle with initial blooms happening in November-December and
December-January in the Arabian Gulf and Sea of Oman, respectively. The
analysis further demonstrates that the algal blooms grow better at salinity
levels of 39-40 psu/37-37.5 psu, temperature of 23-24 oC, and pH of 8 in the
Arabian Gulf/Oman Sea. Findings of this study provide insight into the
relationship between water properties and algal bloom frequency, and a basis
for future research into the drivers behind these observed spatio-temporal
trends.


  We address the existence of periodic orbits for periodic eco-epidemiological
system with disease in the prey. To do it, we consider three main steps.
Firstly we study a one parameter family of systems and obtain uniform bounds
for the components of any periodic solution of these systems. Next, we make a
suitable change of variables in our family of systems to establish the setting
where we are able to apply Mawhin's continuation Theorem. Finally, we use
Mawhin's continuation Theorem to obtain our result. Later on, we present two
examples that include previous results in the literature and some numerical
simulations to illustrate our results.


  We further clarify the relation between detailed-balanced and
complex-balanced equilibria of reversible chemical reaction networks. Our
results hold for arbitrary kinetics and also for boundary equilibria. Detailed
balance, complex balance, "formal balance", and the new notion of "cycle
balance" are all defined in terms of the underlying graph. This fact allows
elementary graph-theoretic (non-algebraic) proofs of a previous result
(detailed balance = complex balance + formal balance), our main result
(detailed balance = complex balance + cycle balance), and a corresponding
result in the setting of continuous-time Markov chains.


  In this note we analyze the temporal dynamics of the coronavirus disease 2019
outbreak in China, Italy and France in the time window $22/01-09/03/2020$. A
first analysis of simple day-lag maps points to some universality in the
epidemic spreading, suggesting that simple mean-field models can be
meaningfully used to gather a quantitative picture of the epidemic spreading,
and notably the height and time of the peak of confirmed infected individuals.
  The analysis of the same data within a simple susceptible-(confirmed)
infected-recovered-deaths model indicates that the kinetic parameter that
describes the rate of recovery seems to be the same, irrespective of the
country, while the infection and death rates appear to be more variable. The
model places the peak in Italy around March 20$^{\rm th}$ 2020, with a maximum
number of confirmed infected individuals of about 16,000.


  The recent advances in sequencing technologies enables the assembly of
individual genomes to the reference quality. How to integrate multiple genomes
from the same species and to make the integrated representation accessible to
biologists remain an open challenge. Here we propose a graph-based data model
and associated formats to represent multiple genomes while preserving the
coordinate of the linear reference genome. We implemented our ideas in the
minigraph toolkit and demonstrate that we can efficiently construct a pangenome
graph and compactly encode tens of thousands of structural variants missing
from the current reference genome.


  Aquatic invasive species (AIS) cause significant ecological and economic
damages around the world. A major spread mechanism for AIS is traffic of
boaters transporting their watercraft from invaded to uninvaded waterbodies. To
inhibit the spread of AIS, several Canadian provinces and American states set
up watercraft inspection stations at roadsides, where potentially infested
boats are screened for AIS and, if necessary, decontaminated. However, since
budgets for AIS control are limited, watercraft inspection stations can only be
operated at specific locations and daytimes. Though theoretical studies provide
managers with general guidelines for AIS management, more specific results are
needed to determine when and where watercraft inspections would be most
effective. This is the subject of this paper. We show how linear integer
programming techniques can be used to optimize watercraft inspection policies
under budget constraints. We introduce our approach as a general framework and
apply it to the prevention of the spread of zebra and quagga mussels (Dreissena
spp.) to the Canadian province British Columbia. We consider a variety of
scenarios and show how variations in budget constraints, propagule sources, and
model uncertainty affect the optimal policy. Based on these results, we
identify simple, generally applicable principles for optimal AIS management.


  The aim of this study was to assess the effect of primary treatment on the
performance of two pilot-scale high rate algal ponds (HRAPs) treating urban
wastewater, considering their treatment efficiency, biomass productivity,
characteristics and biogas production potential. Results indicated that the
primary treatment did not significantly affect the wastewater treatment
efficiency (NH4+-N removal of 93 and 91% and COD removal of 62 and 65% in HRAP
with and without primary treatment, respectively). The HRAP without primary
treatment had higher biodiversity and productivity (18 vs. 16 g VSS/m2d).
Biomass from both systems presented good settling capacity. Results of
biochemical methane potential test showed that co-digesting microalgae and
primary sludge led to higher methane yields (238 - 258 mL CH4/g VS) compared
with microalgae mono-digestion (189 - 225 mL CH4/g VS). Overall, HRAPs with and
without primary treatment seem to be appropriate alternatives for combining
wastewater treatment and bioenergy recovery.


  The aim of this study was to investigate the cultivation of Nostoc sp.,
Arthrospira platensis and Porphyridium purpureum in industrial wastewater to
produce phycobiliproteins. Initially, light intensity and growth medium
composition were optimized, indicating that light conditions influenced the
phycobiliproteins production more than the medium composition. Conditions were
then selected, according to biomass growth, nutrients removal and
phycobiliproteins production, to cultivate these microalgae in food-industry
wastewater. The three species could efficiently remove up to 98%, 94% and 100%
of COD, inorganic nitrogen and PO43--P, respectively. Phycocyanin,
allophycocyanin and phycoerythrin were successfully extracted from the biomass
reaching concentrations up to 103, 57 and 30 mg/g dry weight, respectively.
Results highlight the potential use of microalgae for industrial wastewater
treatment and related high-value phycobiliproteins recovery.


  The aim of this study was to assess the potential environmental impacts
associated with high rate algal ponds (HRAP) systems for wastewater treatment
and resource recovery in small communities. To this aim, a Life Cycle
Assessment (LCA) and an economic assessment were carried out evaluating two
alternatives: i) a HRAPs system for wastewater treatment where microalgal
biomass is valorised for energy recovery (biogas production); ii) a HRAPs
system for wastewater treatment where microalgal biomass is reused for
nutrients recovery (biofertiliser production). Additionally, both alternatives
were compared to a typical small-sized activated sludge system. The results
showed that HRAPs system coupled with biogas production appeared to be more
environmentally friendly than HRAPs system coupled with biofertiliser
production in the climate change, ozone layer depletion, photochemical oxidant
formation, and fossil depletion impact categories. Different climatic
conditions have strongly influenced the results obtained in the eutrophication
and metal depletion impact categories, with the HRAPs system located where warm
temperatures and high solar radiation are predominant showing lower impact. In
terms of costs, HRAPs systems seemed to be more economically feasible when
combined with biofertiliser production instead of biogas.


  One of the important issues in oncology is finding the genes that
perturbation the cell functionality, and result in cancer propagation. The
genes, namely driver genes, when they mutate in expression, result in cancer
through activation of the mutated proteins. So, many methods have been
introduced to predict this group of genes. These are mostly computational
methods based on the number of mutations of each gene. Recently, some
network-based methods have been proposed to predict Cancer Driver Genes (CDGs).
In this study, we use a network-based approach and relative importance of each
gene in the propagation and absorption of genes anomalies in the network to
recognize CDGs. The experimental results are compared with 19 previous methods
that show our proposed algorithm is better than the others in terms of
accuracy, precision, and the number of recognized CDGs.


  Evolution of distribution of strategies in game theory is an interesting
question that has been studied only for specific cases. Here I develop a
general method to extend analysis of the evolution of continuous strategy
distributions given bi-linear and quadratic payoff functions for any initial
distribution to answer the following question: given the initial distribution
of strategies in a game, how will it evolve over time? I look at several
specific examples, including normal distribution on the entire line, normal
truncated distribution, as well as exponential, uniform and Gamma
distributions. I show that the class of exponential distributions is invariant
with respect to replicator dynamics in games with bi-linear payoff functions. I
show also that the class of normal distributions is invariant with respect to
replicator dynamics in games with quadratic payoff functions. The developed
method can now be applied to a broad class of questions pertaining to evolution
of strategies in games with different payoff functions and different initial
distributions.


  After the sudden outbreak of Coronavirus in Wuhan, continuous and rich data
of the epidemic has been made public as the vital fact for decision support in
control measures and aggressive implementation of containment strategies and
plans. With the further growth and spreading of the virus, future resource
allocation and planning under updated strategies and measures rely on careful
study of the epidemic data and characteristics for accurate prediction and
estimation. By using the SIR model and reported data, key parameters are
obtained from least sum of squared errors for an accurate prediction of
epidemic trend in the last four weeks.


  Melanoma is one of the most dangerous skin malignancies due to its high
metastatic tendency and high mortality. Activation of key signaling pathways
enforcing melanoma progression depends on phosphorylation of tyrosine kinases,
and oxidative stress. We here investigated the effect of the new bis-coumarin
derivative (3,5-DCPBC) on human melanoma cell survival, growth, proliferation,
migration, and intracellular redox state, and deciphered associated signal
pathways. This novel derivative was found to be toxic for melanoma cells, and
non-toxic for their benign counterparts, melanocytes and fibroblasts. 3,5-DCPBC
inhibited cell survival, migration and proliferation of different metastatic,
and non-metastatic melanoma cell lines through the profound suppression of
phosphorylation of the Epidermal Growth Factor receptor, and related downstream
pathways. Suppression of phosphorylation of key downstream transcription
factors and different tyrosine kinases comprise JAK/STAT, SRC kinases, ERK and
MAP kinases (p38alpha), all involved in melanoma progression. Simultaneous and
specific targeting of multiple tyrosine kinases and corresponding key genes in
melanoma cells makes 3,5-DCPBC a highly interesting anti-melanoma, and
anti-metastatic drug candidate which may in the long term hold promise in the
therapy of advanced melanoma.


  COVID-19 is an emergent viral infection which rose in December 2019 in a city
in the Chinese province of Hubei, Wuhan; the viral aetiology of this infection
is now known as COVID-19 virus, which belongs to the Betacoronavirus genus.
This virus produces the syndrome of acute respiratory stress that h as been
witnessed in other coronaviruses, such as that MERS-CoV in Middle East
countries or SARS-CoV which was seen in 2002 and 2003 in China. This virus
mediates its entry through its spike (S) proteins interacting with ACE2
receptors in lung epithelial cells, and may promote an inflammatory response by
means of inflammasome NLRP3 activation and unfolded protein response (these are
possibly consequence of the envelope E protein of COVID-19 virus). Efforts have
been made worldwide to prevent further spread of the disease, but in March 2020
the WHO declared it a pandemic emergency and Mexico started to report its first
cases. In this paper we attempt to summarize the biological features of the
virus and the possible pathophysiological mechanisms of its disease, as well as
a stochastic model characterizing the probability distribution of cases in
Mexico by states and the estimated number of cases in Mexico through a
differential equation model (modified SIR model), thus will we be able to
characterize the disease and its course in Mexico in order to display more
preparedness and promote more logical actions by both the policy makers as well
as the general population.


  Promoting some typical spreading dynamics, for instance, the spreading of
information, commercial message, vaccination guidance, innovation, and
political movement, can bring benefits to all aspects of the socio-economic
systems. In this study, we propose a strategy for promoting the spreading of
the susceptible-infected-recovered model, which is widely applied to describe
these common spreading dynamics in real life. Specifically, we first quantify
the potential influence that the addition of each latent edge (that is, edges
that do not exist before) could cause to the spreading dynamics. Then, we
strategically add the latent edges to the original networks according to the
potential influence of each latent edge. Numerical simulations verify the
effectiveness of our strategy and demonstrate that our strategy outperforms
several static strategies, namely, adding the latent edges between nodes with
the largest degree or eigenvector centrality. This study provides an effective
way of promoting the spreading of the susceptible-infected-recovered model by
modifying the network structure slightly and helps in understanding what a
better network structure for the spreading dynamics is. Besides, the
theoretical framework established in this study provides inspirations for the
further investigations of edge-based promoting strategies for other spreading
models.


  In the mitotic spindle microtubules attach to kinetochores via catch bonds
during metaphase, and microtubule depolymerization forces give rise to
stochastic chromosome oscillations. We investigate the cooperative stochastic
microtubule dynamics in spindle models consisting of ensembles of parallel
microtubules, which attach to a kinetochore via elastic linkers. We include the
dynamic instability of microtubules and forces on microtubules and kinetochores
from elastic linkers. A one-sided model, where an external force acts on the
kinetochore is solved analytically employing a mean-field approach based on
Fokker-Planck equations. The solution establishes a bistable force-velocity
relation of the microtubule ensemble in agreement with stochastic simulations.
We derive constraints on linker stiffness and microtubule number for
bistability. The bistable force-velocity relation of the one-sided spindle
model gives rise to oscillations in the two-sided model, which can explain
stochastic chromosome oscillations in metaphase (directional instability). We
derive constraints on linker stiffness and microtubule number for metaphase
chromosome oscillations. Including poleward microtubule flux into the model we
can provide an explanation for the experimentally observed suppression of
chromosome oscillations in cells with high poleward flux velocities. Chromosome
oscillations persist in the presence of polar ejection forces, however, with a
reduced amplitude and a phase shift between sister kinetochores. Moreover,
polar ejection forces are necessary to align the chromosomes at the spindle
equator and stabilize an alternating oscillation pattern of the two
kinetochores. Finally, we modify the model such that microtubules can only
exert tensile forces on the kinetochore resulting in a tug-of-war between the
two microtubule ensembles. Then, induced microtubule catastrophes after
reaching the...


  Models of adaptive bet-hedging commonly adopt insights from Kelly's famous
work on optimal gambling strategies and the financial value of information. In
particular, such models seek evolutionary solutions that maximize long term
average growth rate of lineages, even in the face of highly stochastic growth
trajectories. Here, we argue for extensive departures from the standard
approach to better account for evolutionary contingencies. Crucially, we
incorporate considerations of volatility minimization, motivated by interim
extinction risk in finite populations, within a finite time horizon approach to
growth maximization. We find that a game-theoretic competitive-optimality
approach best captures these additional constraints, and derive the equilibria
solutions under straightforward fitness payoff functions and extinction risks.
We show that for both maximal growth and minimal time relative payoffs the
log-optimal strategy is a unique pure-strategy symmetric equilibrium, invariant
with evolutionary time horizon and robust to low extinction risks.


  Based on the official data modeling, this paper studies the transmission
process of the Corona Virus Disease 2019 (COVID-19). The error between the
model and the official data curve is within 3%. At the same time, it realized
forward prediction and backward inference of the epidemic situation, and the
relevant analysis help relevant countries to make decisions.


  The present study evaluates the removal capacity of two high rate algae ponds
(HRAPs) to eliminate 12 pharmaceuticals (PhACs) and 26 of their corresponding
main metabolites and transformation products. The efficiency of these ponds,
operating with and without primary treatment, was compared in order to study
their capacity under the best performance conditions (highest solar
irradiance). Concentrations of all the target compounds were determined in both
water and biomass samples. Removal rates ranged from moderate (40-60%) to high
(>60%) for most of them, with the exception of the psychiatric drugs
carbamazepine, the beta-blocking agent metoprolol and its metabolite,
metoprolol acid. O-desmethylvenlafaxine, despite its very low biodegradability
in conventional wastewater treatment plants, was removed to certain extent
(13-39%) Biomass concentrations suggested that bioadsorption/bioaccumulation to
microalgae biomass was decisive regarding the elimination of some
non-biodegradable compounds such as venlafaxine and its main metabolites. HRAP
treatment with and without primary treatment did not yield significant
differences in terms of PhACs removal efficiency. The implementation of HRAPs
as secondary treatment is a viable alternative to CAS in terms of overall
wastewater treatment (including organic micropollutants), with generally higher
removal performances and implying a green, low-cost and more sustainable
technology.


  The 2019-Novel-Coronavirus (COVID-19) has affected 116 countries (By March
12) and out of more than 118,000 confirmed cases. Understanding the
transmission dynamics of the infection in each country which affected on a
daily basis and evaluating the effectiveness of control policies is critical
for our further actions. To date, the statistics of COVID-19 reported cases
show more than 80 percent of infected had a mild case of disease, while around
14 percent of infected experienced a severe one and about 5 percent are
categorized as critical disease victims. Today's report (2020-03-12; daily
updates in the prepared website) shows the confirmed cases of COVID-19 in
China, South Korea, Italy, and Iran are 80932, 7869, 12462 and 10075;
respectively. Calculating the total Case Fatality Rate (CFR) of Italy
(2020-03-04), about 7.9% of confirmed cases passed away. Compared to South
Korea's rate of 0.76% (10 times lower than Italy) and China's 3.8% (50% lower
than Italy), the CFR of Italy is too high. There are some effective policies
that yield significant changes in the trend of cases. The lockdown policy in
China and Italy (the effect observed after 11 days), Shutdown of all
nonessential companies in Hubei (the effect observed after 5 days), combined
policy in South Korea and reducing working hours in Iran.


  The modeling of the spreading of communicable diseases has experienced
significant advances in the last two decades or so. This has been possible due
to the proliferation of data and the development of new methods to gather, mine
and analyze it. A key role has also been played by the latest advances in new
disciplines like network science. Nonetheless, current models still lack a
faithful representation of all possible heterogeneities and features that can
be extracted from data. Here, we bridge a current gap in the mathematical
modeling of infectious diseases and develop a framework that allows to account
simultaneously for both the connectivity of individuals and the age-structure
of the population. We compare different scenarios, namely, i) the homogeneous
mixing setting, ii) one in which only the social mixing is taken into account,
iii) a setting that considers the connectivity of individuals alone, and
finally, iv) a multilayer representation in which both the social mixing and
the number of contacts are included in the model. We analytically show that the
thresholds obtained for these four scenarios are different. In addition, we
conduct extensive numerical simulations and conclude that heterogeneities in
the contact network are important for a proper determination of the epidemic
threshold, whereas the age-structure plays a bigger role beyond the onset of
the outbreak. Altogether, when it comes to evaluate interventions such as
vaccination, both sources of individual heterogeneity are important and should
be concurrently considered. Our results also provide an indication of the
errors incurred in situations in which one cannot access all needed information
in terms of connectivity and age of the population.


  This paper studies if and to which extent COVID-19 epidemics can be
controlled by authorities taking decisions on public health measures on the
basis of daily reports of swab test results, active cases and total cases. A
suitably simplified process model is derived to support the controllability
analysis, highlighting the presence of very significant time delay; the model
is validated with data from several outbreaks. The analysis shows that
suppression strategies can be effective if strong enough and enacted early on.
It also shows how mitigation strategies can fail because of the combination of
delay, unstable dynamics, and uncertainty in the feedback loop; approximate
conditions based on the theory of limitation of linear control are given for
feedback control to be feasible.


  We study the datafor the cumulative as well as daily number of cases in the
Covid-19 outbreak in China. The cumulative data can be fit to an empirical form
obtained from a Susceptible-Infected-Removed (SIR) model studied on an
Euclidean network previously. Plotting the number of cases against the distance
from the epicenter for both China and Italy, we find an approximate power law
variation with an exponent $\sim 1.85$ showing strongly that the spatial
dependence plays a key role, a factor included in the model. We report here
that the SIR model on the Eucledean network can reproduce with a high accuracy
the data for China for given parameter values, and can also predict when the
epidemic, at least locally, can be expected to be over.


  Although displaying genetic correlations, psychiatric disorders are
clinically defined as categorical entities as they each have distinguishing
clinical features and may involve different treatments. Identifying
differential genetic variations between these disorders may reveal how the
disorders differ biologically and help to guide more personalized treatment.
  Here we presented a comprehensive analysis to identify genetic markers
differentially associated with various psychiatric disorders/traits based on
GWAS summary statistics, covering 18 psychiatric traits/disorders and 26
comparisons. We also conducted comprehensive analysis to unravel the genes,
pathways and SNP functional categories involved, and the cell types and tissues
implicated. We also assessed how well one could distinguish between psychiatric
disorders by polygenic risk scores (PRS).
  SNP-based heritabilities (h2SNP) were significantly larger than zero for most
comparisons. Based on current GWAS data, PRS have mostly modest power to
distinguish between psychiatric disorders. For example, we estimated that AUC
for distinguishing schizophrenia from major depressive disorder (MDD), bipolar
disorder (BPD) from MDD and schizophrenia from BPD were 0.694, 0.602 and 0.618
respectively, while the maximum AUC (based on h2SNP) were 0.763, 0.749 and
0.726 respectively. We also uncovered differences in each pair of studied
traits in terms of their differences in genetic correlation with comorbid
traits. For example, clinically-defined MDD appeared to more strongly
genetically correlated with other psychiatric disorders and heart disease, when
compared to non-clinically-defined depression in UK Biobank.
  Our findings highlight genetic differences between psychiatric disorders and
the mechanisms involved. PRS may aid differential diagnosis of selected
psychiatric disorders in the future with larger GWAS samples.


  The Selkov oscillator, a simple description of glycolysis, is a system of two
ordinary differential equations with mass action kinetics. In previous work the
authors established several properties of the solutions of this system. In the
present paper we extend this to prove that this system has solutions which
diverge to infinity in an oscillatory manner at late times. This system was
originally derived from another system with Michaelis-Menten kinetics. It is
shown that the Michaelis-Menten system, like that with mass action, has
solutions which diverge to infinity in a monotone manner. It is also shown to
admit subcritical Hopf bifurcations and thus unstable periodic solutions. We
discuss to what extent the unbounded solutions cast doubt on the biological
relevance of the Selkov oscillator and compare it with other models in the
literature.


  Variant calling, the problem of estimating whether a position in a DNA
sequence differs from a reference sequence, given noisy, redundant, overlapping
short sequences that cover that position, is fundamental to genomics. We
propose a deep averaging network designed specifically for variant calling. Our
model takes into account the independence of each short input read sequence by
transforming individual reads through a series of convolutional layers,
limiting the communication between individual reads to averaging and
concatenating operations. Training and testing on the precisionFDA Truth
Challenge (pFDA), we match state of the art overall 99.89 F1 score. Genome
datasets exhibit extreme skew between easy examples and those on the decision
boundary. We take advantage of this property to converge models at 5x the speed
of standard epoch-based training by skipping easy examples during training. To
facilitate future work, we release our code, trained models and pre-processed
public domain datasets.


  We show that the space of rooted tree-based phylogenetic networks is
connected under rooted nearest-neighbour interchange (rNNI) moves.


  Within a short period of time, COVID-19 grew into a world-wide pandemic.
Transmission by pre-symptomatic and asymptomatic viral carriers rendered
intervention and containment of the disease extremely challenging. Based on
reported infection case studies, we construct an epidemiological model that
focuses on transmission around the symptom onset. The model is calibrated
against incubation period and pairwise transmission statistics during the
initial outbreaks of the pandemic outside Wuhan with minimal non-pharmaceutical
interventions. Mathematical treatment of the model yields explicit expressions
for the size of latent and pre-symptomatic subpopulations during the
exponential growth phase, with the local epidemic growth rate as input. We then
explore reduction of the basic reproduction number R_0 through specific disease
control measures such as contact tracing, testing, social distancing, wearing
masks and sheltering in place. When these measures are implemented in
combination, their effects on R_0 multiply. We also compare our model behaviour
to the first wave of the COVID-19 spreading in various affected regions and
highlight generic and less generic features of the pandemic development.


  At the time of this writing, the novel coronavirus (COVID-19) pandemic
outbreak has already put tremendous strain on many countries' citizens,
resources and economies around the world. Social distancing measures, travel
bans, self-quarantines, and business closures are changing the very fabric of
societies worldwide. With people forced out of public spaces, much conversation
about these phenomena now occurs online, e.g., on social media platforms like
Twitter. In this paper, we describe a multilingual coronavirus (COVID-19)
Twitter dataset that we have been continuously collecting since January 22,
2020. We are making our dataset available to the research community
(https://github.com/echen102/COVID-19-TweetIDs). It is our hope that our
contribution will enable the study of online conversation dynamics in the
context of a planetary-scale epidemic outbreak of unprecedented proportions and
implications. This dataset could also help track scientific coronavirus
misinformation and unverified rumors, or enable the understanding of fear and
panic -- and undoubtedly more. Ultimately, this dataset may contribute towards
enabling informed solutions and prescribing targeted policy interventions to
fight this global crisis.


  Coronaviruses are a famous family of viruses that cause illness in both
humans and animals. The new type of coronavirus COVID-19 was firstly discovered
in Wuhan, China. However, recently, the virus has widely spread in most of the
world and causing a pandemic according to the World Health Organization (WHO).
Further, nowadays, all the world countries are striving to control the
COVID-19. There are many mechanisms to detect coronavirus including clinical
analysis of chest CT scan images and blood test results. The confirmed COVID-19
patient manifests as fever, tiredness, and dry cough. Particularly, several
techniques can be used to detect the initial results of the virus such as
medical detection Kits. However, such devices are incurring huge cost, taking
time to install them and use. Therefore, in this paper, a new framework is
proposed to detect COVID-19 using built-in smartphone sensors. The proposal
provides a low-cost solution, since most of radiologists have already held
smartphones for different daily-purposes. Not only that but also ordinary
people can use the framework on their smartphones for the virus detection
purposes. Nowadays Smartphones are powerful with existing computation-rich
processors, memory space, and large number of sensors including cameras,
microphone, temperature sensor, inertial sensors, proximity, colour-sensor,
humidity-sensor, and wireless chipsets/sensors. The designed Artificial
Intelligence (AI) enabled framework reads the smartphone sensors signal
measurements to predict the grade of severity of the pneumonia as well as
predicting the result of the disease.


  In recent biomedical scientific problems, it is a fundamental issue to
integratively cluster a set of objects from multiple sources of datasets. Such
problems are mostly encountered in genomics, where data is collected from
various sources, and typically represent distinct yet complementary
information. Integrating these data sources for multi-source clustering is
challenging due to their complex dependence structure including directional
dependency. Particularly in genomics studies, it is known that there is certain
directional dependence between DNA expression, DNA methylation, and RNA
expression, widely called The Central Dogma.
  Most of the existing multi-view clustering methods either assume an
independent structure or pair-wise (non-directional) dependency, thereby
ignoring the directional relationship. Motivated by this, we propose a
copula-based multi-view clustering model where a copula enables the model to
accommodate the directional dependence existing in the datasets. We conduct a
simulation experiment where the simulated datasets exhibiting inherent
directional dependence: it turns out that ignoring the directional dependence
negatively affects the clustering performance. As a real application, we
applied our model to the breast cancer tumor samples collected from The Cancer
Genome Altas (TCGA).


  We present a robust data-driven machine learning analysis of the COVID-19
pandemic from its early infection dynamics, specifically infection counts over
time. The goal is to extract actionable public health insights. These insights
include the infectious force, the rate of a mild infection becoming serious,
estimates for asymtomatic infections and predictions of new infections over
time. We focus on USA data starting from the first confirmed infection on
January 20 2020. Our methods reveal significant asymptomatic (hidden)
infection, a lag of about 10 days, and we quantitatively confirm that the
infectious force is strong with about a 0.14% transition from mild to serious
infection. Our methods are efficient, robust and general, being agnostic to the
specific virus and applicable to different populations or cohorts.


  In this paper we introduce a new representation for the multistationarity
region of a reaction network, using polynomial superlevel sets. The advantages
of using this polynomial superlevel set representation over the already
existing representations (cylindrical algebraic decompositions, numeric
sampling, rectangular divisions) is discussed, and algorithms to compute this
new representation are provided. The results are given for the general
mathematical formalism of a parametric system of equations and so may be
applied to other application domains.


  In mid of March 2020, Coronaviruses such as COVID-19 is declared as an
international epidemic. More than 125000 confirmed cases and 4,607 death cases
have been recorded around more than 118 countries. Unfortunately, a coronavirus
vaccine is expected to take at least 18 months if it works at all. Moreover,
COVID -19 epidemics can mutate into a more aggressive form. Day level
information about the COVID -19 spread is crucial to measure the behavior of
this new virus globally. Therefore, this study presents a comparison of day
level forecasting models on COVID-19 affected cases using time series models
and mathematical formulation. The forecasting models and data strongly suggest
that the number of coronavirus cases grows exponentially in countries that do
not mandate quarantines, restrictions on travel and public gatherings, and
closing of schools, universities, and workplaces (Social Distancing).


  A principal component analysis of the TCGA data for 15 cancer localizations
unveils the following qualitative facts about tumors: 1) The state of a tissue
in gene expression space may be described by a few variables. In particular,
there is a single variable describing the progression from a normal tissue to a
tumor. 2) Each cancer localization is characterized by a gene expression
profile, in which genes have specific weights in the definition of the cancer
state. There are no less than 2500 differentially-expressed genes, which lead
to power-like tails in the expression distribution functions. 3) Tumors in
different localizations share hundreds or even thousands of differentially
expressed genes. There are 6 genes common to the 15 studied tumor
localizations. 4) The tumor region is a kind of attractor. Tumors in advanced
stages converge to this region independently of patient age or genetic
variability. 5) There is a landscape of cancer in gene expression space with an
approximate border separating normal tissues from tumors.


  During the COVID-19 outbreak, it is essential to monitor the effectiveness of
measures taken by governments on the course of the epidemic. Here we show that
there is already a sufficient amount of data collected in Italy to predict the
outcome of the process. We show that using the proper metric, the data from
Hubei Province and Italy has striking similarity, which enables us to calculate
the expected number of confirmed cases and the number of deaths by the end of
the process. Our predictions will improve as new data points are generated day
by day, which can help to make further public decisions. The method is based on
the data analysis of logistic growth equations describing the process on the
macroscopic level. At the time of writing of the first version, the number of
fatalities in Italy was expected to be 6000, and the predicted end of the
crisis was April 15, 2020. In this new version, we discuss what changed in the
two weeks which passed since then. The trend changed drastically on March 17,
2020, when the Italian health system reached its capacity limit. Without this
limit, probably 3500 more people would have died. Instead, due to the
limitations, 17.000 people are expected to die now, which is a five-fold
increase. The predicted end of the crisis now shifted to May 8, 2020.


  An explanatory model for the emergence of evolvable units must display
emerging structures that (1) preserve themselves in time (2) self-reproduce and
(3) tolerate a certain amount of variation when reproducing. To tackle this
challenge, here we introduce Combinatory Chemistry, an Algorithmic Artificial
Chemistry based on a minimalistic computational paradigm named Combinatory
Logic. The dynamics of this system comprise very few rules, it is initialised
with an elementary tabula rasa state, and features conservation laws
replicating natural resource constraints. Our experiments show that a single
run of this dynamical system with no external intervention discovers a wide
range of emergent patterns. All these structures rely on acquiring basic
constituents from the environment and decomposing them in a process that is
remarkably similar to biological metabolisms. These patterns include
autopoietic structures that maintain their organisation, recursive ones that
grow in linear chains or binary-branching trees, and most notably, patterns
able to reproduce themselves, duplicating their number at each generation.


  In this paper we will consider a mathematical model that describes, the
tritrophic interaction between plants, herbivores and their natural enemies,
where volatiles organic compounds (VOCs) released by plants play an important
role. We show positivity and boundedness of the system solutions, existence of
positive equilibrium and its local stability, we analyse global stability of
positive equilibrium via the geometrical approach of Li and Muldowney. We pay
attention to parameters in order to discuss different types of bifurcations.
Finally, we present some numerical simulations to justify our analytical
results.


  Protein-ligand interactions are one of the fundamental types of molecular
interactions in living systems. Ligands are small molecules that interact with
protein molecules at specific regions on their surfaces called binding sites.
Tasks such as assessment of protein functional similarity and detection of side
effects of drugs need identification of similar binding sites of disparate
proteins across diverse pathways. Machine learning methods for similarity
assessment require feature descriptors of binding sites. Traditional methods
based on hand engineered motifs and atomic configurations are not scalable
across several thousands of sites. In this regard, deep neural network
algorithms are now deployed which can capture very complex input feature space.
However, one fundamental challenge in applying deep learning to structures of
binding sites is the input representation and the reference frame. We report
here a novel algorithm Site2Vec that derives reference frame invariant vector
embedding of a protein-ligand binding site. The method is based on pairwise
distances between representative points and chemical compositions in terms of
constituent amino acids of a site. The vector embedding serves as a locality
sensitive hash function for proximity queries and determining similar sites.
The method has been the top performer with more than 95% quality scores in
extensive benchmarking studies carried over 10 datasets and against 23 other
site comparison methods. The algorithm serves for high throughput processing
and has been evaluated for stability with respect to reference frame shifts,
coordinate perturbations and residue mutations. We provide Site2Vec as a stand
alone executable and a web service hosted at
\url{http://services.iittp.ac.in/bioinfo/home}.


  An outbreak of respiratory disease caused by a novel coronavirus is ongoing
from December 2019. As of July 22, 2020, it has caused an epidemic outbreak
with more than 15 million confirmed infections and above 6 hundred thousand
reported deaths worldwide. During this period of an epidemic when
human-to-human transmission is established and reported cases of coronavirus
disease 2019 (COVID-19) are rising worldwide, investigation of control
strategies and forecasting are necessary for health care planning. In this
study, we propose and analyze a compartmental epidemic model of COVID-19 to
predict and control the outbreak. The basic reproduction number and control
reproduction number are calculated analytically. A detailed stability analysis
of the model is performed to observe the dynamics of the system. We calibrated
the proposed model to fit daily data from the United Kingdom (UK) where the
situation is still alarming. Our findings suggest that independent
self-sustaining human-to-human spread ($R_0>1$, $R_c>1$) is already present.
Short-term predictions show that the decreasing trend of new COVID-19 cases is
well captured by the model. Further, we found that effective management of
quarantined individuals is more effective than management of isolated
individuals to reduce the disease burden. Thus, if limited resources are
available, then investing on the quarantined individuals will be more fruitful
in terms of reduction of cases.


  Recently, emerging epidemics like COVID-19 and its variants require
predictive mathematical models to implement suitable responses in order to
limit their negative and profound impact on society. The SIR
(Susceptible-InfectedRemoved) system is a straightforward mathematical
formulation to model the dissemination of many infectious diseases. The present
paper reports novel theoretical and analytical results for a perturbed version
of an SIR model with Gamma-distributed delay. Notably, our epidemic model is
represented by It\^o-L'evy stochastic differential equations in order to
simulate sudden and unexpected external phenomena. By using some new and
ameliorated mathematical approaches, we study the long-run characteristics of
the perturbed delayed model. Within this scope, we give sufficient conditions
for two interesting asymptotic proprieties: extinction and persistence of the
epidemic. One of the most interesting results is that the dynamics of the
stochastic model are closely related to the intensities of white noises and
L'evy jumps, which can give us a good insight into the evolution of the
epidemic in some unexpected situations. Our work complements the results of
some previous investigations and provides a new approach to predict and analyze
the dynamic behavior of epidemics with distributed delay. For illustrative
purposes, numerical examples are presented for checking the theoretical study


  A model of reactive social distancing in epidemics is proposed, in which the
infection rate changes with the number infected. The final-size equation for
the total number that the epidemic will infect can be derived analytically, as
can the peak infection proportion. This model could assist planners during for
example the COVID-19 pandemic.


  The goal of this work is to study an infectious disease spreading in a medium
size population occupying a confined environment. For this purpose, we consider
a kinetic theory approach to model crowd dynamics in bounded domains and couple
it to a kinetic equation to model contagion. The interactions of a person with
other pedestrians and the environment are modeled by using tools of game
theory. The pedestrian dynamics model allows to weight between two competing
behaviors: the search for less congested areas and the tendency to follow the
stream unconsciously in a panic situation. Each person in the system has a
contagion level that is affected by their neighborhood. For the numerical
solution of the coupled problem, we propose a numerical algorithm that at every
time step solves one crowd dynamics problem and one contagion problem, i.e.
with no subiterations between the two. We test our coupled model on a problem
involving a small crowd walking through a corridor.


  The new coronavirus known as COVID-19 is spread world-wide since December
2019. Without any vaccination or medicine, the means of controlling it are
limited to quarantine and social distancing. Here we study the spatio-temporal
propagation of the first wave of the COVID-19 virus in China and compare it to
other global locations. We provide a comprehensive picture of the spatial
propagation from Hubei to other provinces in China in terms of distance,
population size, and human mobility and their scaling relations. Since strict
quarantine has been usually applied between cities, more insight about the
temporal evolution of the disease can be obtained by analyzing the epidemic
within cities, especially the time evolution of the infection, death, and
recovery rates which affected by policies. We study and compare the infection
rate in different cities in China and provinces in Italy and find that the
disease spread is characterized by a two-stages process. At early times, at
order of few days, the infection rate is close to a constant probably due to
the lack of means to detect infected individuals before infection symptoms are
observed. Then at later times it decays approximately exponentially due to
quarantines. The time evolution of the death and recovery rates also
distinguish between these two stages and reflect the health system situation
which could be overloaded.


  The fast and untraceable virus mutations take lives of thousands of people
before the immune system can produce the inhibitory antibody. Recent outbreak
of novel coronavirus infected and killed thousands of people in the world.
Rapid methods in finding peptides or antibody sequences that can inhibit the
viral epitopes of COVID-19 will save the life of thousands. In this paper, we
devised a machine learning (ML) model to predict the possible inhibitory
synthetic antibodies for Corona virus. We collected 1933 virus-antibody
sequences and their clinical patient neutralization response and trained an ML
model to predict the antibody response. Using graph featurization with variety
of ML methods, we screened thousands of hypothetical antibody sequences and
found 8 stable antibodies that potentially inhibit COVID-19. We combined
bioinformatics, structural biology, and Molecular Dynamics (MD) simulations to
verify the stability of the candidate antibodies that can inhibit the Corona
virus.


  Motivated with various responses of world governments to COVID-19, here we
develop a toy model of the dependence epidemics spreading on the availability
of tests for disease. Our model, that we call SUDR+K, is based on usual SIR
model, but it splits the total fraction of infected individuals into two
components: those that are undetected and those that are detected through
tests. Moreover, we assume that available tests increase at a constant rate
from the beginning of epidemics but are consumed to detect infected
individuals. Strikingly we find a bi-stable behavior between a phase with a
giant fraction of infected and a phase with a very small fraction. We show that
the separation between these two regimes is governed by a match between the
rate of testing and a rate of infection spread at given time. We also show that
the existence of two phases does not depend on the mathematical choice of the
form of the term describing the rate at which undetected individuals are tested
and detected. Presented research implies that a vigorous early testing
activity, before the epidemics enters into its giant phase, can potentially
keep epidemics under control, and that even a very small change in rate of
testing can increase or decrease the size of the whole epidemics of various
orders of magnitude. For the real application of realistic model to ongoing
epidemics, we would gladly collaborate with field epidemiologists in order to
develop quantitative models of testing process.


  Genome-wide association studies(GWAS) have proven to be highly useful in
revealing the genetic basis of complex diseases. At present, most GWAS are
studies of a particular single disease diagnosis against controls. However, in
practice, an individual is often affected by more than one condition/disorder.
For example, patients with coronary artery disease(CAD) are often comorbid with
diabetes mellitus(DM). Along a similar line, it is often clinically meaningful
to study patients with one disease but without a comorbidity. For example,
obese DM may have different pathophysiology from non-obese DM.
  Here we developed a statistical framework to uncover susceptibility variants
for comorbid disorders (or a disorder without comorbidity), using GWAS summary
statistics only. In essence, we mimicked a case-control GWAS in which the cases
are affected with comorbidities or a disease without a relevant comorbid
condition (in either case, we may consider the cases as those affected by a
specific subtype of disease, as characterized by the presence or absence of
comorbid conditions). We extended our methodology to deal with continuous
traits with clinically meaningful categories (e.g. lipids). In addition, we
illustrated how the analytic framework may be extended to more than two traits.
We verified the feasibility and validity of our method by applying it to
simulated scenarios and four cardiometabolic (CM) traits. We also analyzed the
genes, pathways, cell-types/tissues involved in CM disease subtypes. LD-score
regression analysis revealed some subtypes may indeed be biologically distinct
with low genetic correlations. Further Mendelian randomization analysis found
differential causal effects of different subtypes to relevant complications. We
believe the findings are of both scientific and clinical value, and the
proposed method may open a new avenue to analyzing GWAS data.


  Epigenetic modifications of histones crucially affect the eukaryotic gene
activity, while the epigenetic histone state is largely determined by the
binding of specific factors such as the transcription factors (TFs) to DNA.
Here, the way how the TFs and the histone state are dynamically correlated is
not obvious when the TF synthesis is regulated by the histone state. This type
of feedback regulatory relations are ubiquitous in gene networks to determine
cell fate in differentiation and other cell transformations. To gain insights
into such dynamical feedback regulations, we theoretically analyze a model of
epigenetic gene switching by extending the Doi-Peliti operator formalism of
reaction kinetics to the problem of coupled molecular processes. The spin-1 and
spin-1/2 coherent state representations are introduced to describe stochastic
reactions of histones and binding/unbinding of TF in a unified way, which
provides a concise view of the effects of timescale difference among these
molecular processes; even in the case that binding/unbinding of TF to/from DNA
are adiabatically fast, the slow nonadiabatic histone dynamics give rise to a
distinct circular flow of the probability flux around basins in the landscape
of the gene state distribution, which leads to hysteresis in gene switching. In
contrast to the general belief that the change in the amount of TF precedes the
histone state change, the flux drives histones to be modified prior to the
change in the amount of TF in the self-regulating circuits. The flux-landscape
analyses shed light on the nonlinear nonadiabatic mechanism of epigenetic cell
fate decision making.


  The relative case fatality rates (CFRs) between groups and countries are key
measures of relative risk that guide policy decisions regarding scarce medical
resource allocation during the ongoing COVID-19 pandemic. In the middle of an
active outbreak when surveillance data is the primary source of information,
estimating these quantities involves compensating for competing biases in time
series of deaths, cases, and recoveries. These include time- and severity-
dependent reporting of cases as well as time lags in observed patient outcomes.
In the context of COVID-19 CFR estimation, we survey such biases and their
potential significance. Further, we analyze theoretically the effect of certain
biases, like preferential reporting of fatal cases, on naive estimators of CFR.
We provide a partially corrected estimator of these naive estimates that
accounts for time lag and imperfect reporting of deaths and recoveries. We show
that collection of randomized data by testing the contacts of infectious
individuals regardless of the presence of symptoms would mitigate bias by
limiting the covariance between diagnosis and death. Our analysis is
supplemented by theoretical and numerical results and a simple and fast
open-source codebase at https://github.com/aangelopoulos/cfr-covid-19 .


  The sudden spread of Covid-19 outside China has pushed on March 11 the World
Health Organization to acknowledge the ongoing outbreak as a pandemic. It is
crucial in this phase to understand what should countries which presently lag
behind in the spread of the infection learn from countries where the infection
spread earlier. The choice of this work is to prefer timeliness to
comprehensiveness. By adopting a purely empirical approach, we will limit
ourselves to identifying different phases in the plots of different countries,
based on their different functional behaviour, and to make a comparative
analysis. The comparative analysis of the registered cases curves highlights
remarkable similarities, especially among Western countries, together with some
minor but crucial differences. We highlight how timeliness can largely reduce
the size of the individual national outbreaks, ultimately limiting the final
death toll. Our data suggest that Western governments have not unfortunately
shown the capability to anticipate their decisions, based on the experience of
countries hit earlier by the outbreak.


  We employ methods largely exploited in Physics, in the analysis of the
evolution of dynamical systems, to study the pattern of the Covid-19 infection
in Italy. The techniques we employ are based on the use of logistic function
and of its derivative, namely the Hubbert function. The latter is exploited to
give a prediction on the number of infected per day. We also mention the
possibility of taking advantage from other mathematical tools based e.g. on the
Gompertz equation and make some comparison on the different predictive
capabilities.