  We consider the general problem of describing the dynamics of subnetworks of
larger biochemical reaction networks, e.g. protein interaction networks
involving complex formation and dissociation reactions. We propose the use of
model reduction strategies to understand the 'extrinsic' sources of
stochasticity arising from the rest of the network. Our approaches are based on
subnetwork dynamical equations derived by projection methods and by path
integrals. The results provide a principled derivation of the different
components of the extrinsic noise that is observed experimentally in cellular
biochemical reactions, over and above the intrinsic noise from the
stochasticity of biochemical events in the subnetwork. We explore several
intermediate approximations to assess systematically the relative importance of
different extrinsic noise components, including initial transients, long-time
plateaus, temporal correlations, multiplicative noise terms and nonlinear noise
propagation. The best approximations achieve excellent accuracy in quantitative
tests on a simple protein network and on the epidermal growth factor receptor
signalling network.


  There is increasing evidence that one of the most difficult problems in
trying to control the ongoing COVID-19 epidemic is the presence of a large
cohort of asymptomatic infectives. We develop a SIR-type model taking into
account the presence of asymptomatic, or however undetected, infective, and the
substantially long time these spend being infective and not isolated. We
discuss how a SIR-based prediction of the epidemic course based on early data
but not taking into account the presence of a large set of asymptomatic
infectives would give wrong estimate of very relevant quantities such as the
need of hospital beds, the time to the epidemic peak, and the number of people
which are left untouched by the first wave and thus in danger in case of a
second epidemic wave. In the second part of the note, we apply our model to the
COVID-19 epidemics in Italy. We obtain a good agreement with epidemiological
data; according to the best fit of epidemiological data in terms of this model,
only 10\% of infectives in Italy is symptomatic.


  Using a stochastic Susceptible-Infected-Removed (SIR) meta-population model
of disease transmission, we present analytical calculations and numerical
simulations dissecting the interplay between stochasticity and the division of
a population into mutually independent sub-populations. We show that
subdivision activates two stochastic effects---extinction and
desynchronization---diminishing the overall impact of the outbreak, even when
the total population has already left the stochastic regime and the basic
reproduction number is not altered by the subdivision. Both effects are
quantitatively captured by our theoretical estimates, allowing us to determine
their individual contributions to the observed reduction of the peak of the
epidemic.


  An epidemiological model for COVID-19 was developed and implemented in
MATLAB/GNU Octave for use by public health practitioners, policy makers and the
general public. The model distinguishes four stages in the disease: infected,
sick, seriously sick, and better. The model was preliminarily parameterized
based on observations of the spread of the disease. The model is consistent
with a mortality rate of 1.5 %. Preliminary simulations with the model indicate
that concepts such as "herd immunity" and "flattening the curve" are highly
misleading in the context of this virus. Public policies based on these
concepts are inadequate to protect the population. Only reducing the R0 of the
virus below 1 is an effective strategy for maintaining the death burden of
COVID-19 within the normal range of seasonal flu. As R0 values estimated with
the model range from 2.82 worldwide outside of China and 3.83 in the Western
world in late February - early March 2020, this means social distancing with
effectiveness greater than 65 % (worldwide) or 75 % (Western world) are needed
to combat the virus successfully.


  Here, we focus on the data analysis of the growth of epidemic spread of
Covid-19 in countries where different policies of containment were activated.
It is known that the growth of pandemic spread at its threshold is exponential,
but it is not known how to quantify the success of different containment
policies. We identify that a successful approach gives an arrested phase regime
following the Ostwald growth, where, over the course of time, one phase
transforms into another metastable phase with a similar free energy as observed
in oxygen interstitial diffusion in quantum complex matter and in
crystallization of proteins. We introduce the s factor which provides a
quantitative measure of the efficiency and speed of the adopted containment
policy, which is very helpful not only to monitor the Covid-19 pandemic spread
but also for other countries to choose the best containment policy. The results
show that a policy based on joint confinement, targeted tests, and tracking
positive cases is the most rapid pandemic containment policy; in fact, we found
values of 9, 5, and 31 for the success s factor for China, South Korea, and
Italy, respectively, where the lowest s factor indicates the best containment
policy


  Identifying the nodes that have the potential to influence the state of a
network is a relevant question for many complex systems. In many applications
it is often essential to test the ability of an individual node to control a
specific target subset of the network. In biological networks, this might
provide precious information on how single genes regulate the expression of
specific groups of molecules in the cell. Taking into account these
constraints, we propose an optimized heuristic based on the Kalman rank
condition to quantify the centrality of a node as the number of target nodes it
can control. By introducing a hierarchy among the nodes in the target set, and
performing a step-wise research, we ensure for sparse and directed networks the
identification of a controllable driver-target configuration in a significantly
reduced space and time complexity. We show how the method works for simple
network configurations, then we use it to characterize the inflammatory
pathways in molecular gene networks associated with macrophage dysfunction in
patients with multiple sclerosis. Results indicate that the targeted secreted
molecules can in general be controlled by a large number of driver nodes (51%)
involved in different cell functions, i.e. sensing, signaling and
transcription. However, during the inflammatory response only a moderate
fraction of all the possible driver-target pairs are significantly coactivated,
as measured by gene expression data obtained from human blood samples. Notably,
they differ between multiple sclerosis patients and healthy controls, and we
find that this is related to the presence of dysregulated genes along the
controllable walks. Our method, that we name step-wise target controllability,
represents a practical solution to identify controllable driver-target
configurations in directed complex networks and test their relevance from a
functional perspective.


  The Price equation partitions the change in the expected value of a
population measure. The first component describes the partial change caused by
altered frequencies. The second component describes the partial change caused
by altered measurements. In biology, frequency changes often associate with the
direct effect of natural selection. Measure changes reflect processes during
transmission that alter trait values. More broadly, the two components describe
the direct forces that change population composition and the altered frame of
reference that changes measured values. The classic Price equation is limited
to population statistics that can expressed as the expected value of a measure.
Many statistics cannot be expressed as expected values, such as the harmonic
mean and the family of rescaled diversity measures. We generalize the Price
equation to any population statistic that can be expressed as a function of
frequencies and measurements. We obtain the generalized partition between the
direct forces that cause frequency change and the altered frame of reference
that changes measurements.


  The recent epidemic of Coronavirus (COVID-19) that started in China has
already been "exported" to more than 140 countries in all the continents,
evolving in most of them by local spreading. In this contribution we analyze
the trends of the cases reported in all the Chinese provinces, as well as in
some countries that, until March 15th, 2020, have more than 500 cases reported.
Notably and differently from other epidemics, the provinces did not show an
exponential phase. The data available at the Johns Hopkins University site seem
to fit well an algebraic sub-exponential growing behavior as was pointed out
recently. All the provinces show a clear and consistent pattern of slowing down
with growing exponent going nearly zero, so it can be said that the epidemic
was contained in China. On the other side, the more recent spread in countries
like, Italy, Iran, and Spain show a clear exponential growth, as well as other
European countries. Even more recently, US -which was one of the first
countries to have an individual infected outside China (Jan 21st, 2020)- seems
to follow the same path. We calculate the exponential growth of the most
affected countries, showing the evolution along time after the first local
case. We identify clearly different patterns in the analyzed data and we give
interpretations and possible explanations for them. The analysis and
conclusions of our study can help countries that, after importing some cases,
are not yet in the local spreading phase, or have just started.


  We derive and asymptotically analyze mass-action models for disease spread
that include transient pair formation and dissociation. Populations of unpaired
susceptibles and infecteds are distinguished from the population of three types
of pairs of individuals; both susceptible, one susceptible and one infected,
and both infected. Disease transmission can occur only within a pair consisting
of one susceptible individual and one infected individual. By considering the
fast pair formation and fast pair dissociation limits, we use a perturbation
expansion to formally derive a uniformly valid approximation for the dynamics
of the total infected and susceptible populations. Under different parameter
regimes, we derive uniformly valid effective equations for the total infected
population and compare their results to those of the full mass-action model.
Our results are derived from the fundamental mass-action system without
implicitly imposing transmission mechanisms such as that used in
frequency-dependent models. They provide a new formulation for effective
pairing models and are compared with previous models.


  A deeper understanding of the processes underlying the distribution of
animals in space is crucial for both basic and applied ecology. The Common
buzzard (Buteo buteo) is a highly aggressive, territorial bird of prey that
interacts strongly with its intra- and interspecific competitors. We propose
and use random matrix theory to quantify the strength and range of repulsion as
a function of the buzzard population density, thus providing a novel approach
to model density dependence. As an indicator of territorial behaviour, we
perform a large-scale analysis of the distribution of buzzard nests in an area
of $300$ square kilometres around the Teutoburger Wald, Germany, as gathered
over a period of $20$ years. The nearest and next-to-nearest neighbour spacing
distribution between nests is compared to the two-dimensional Poisson
distribution, originating from uncorrelated random variables, to the complex
eigenvalues of random matrices, which are strongly correlated, and to a
two-dimensional Coulomb gas interpolating between these two. A one-parameter
fit to a time-moving average reveals a significant increase of repulsion
between neighbouring nests, as a function of the observed increase in absolute
population density over the monitored period of time, thereby proving an
unexpected yet simple model for density-dependent spacing of predator
territories. A similar effect is obtained for next-to-nearest neighbours,
albeit with weaker repulsion, indicating a short-range interaction. Our results
show that random matrix theory might be useful in the context of population
ecology.


  In the night of February 20, 2020, the first case of novel coronavirus
disease (COVID-19) was confirmed in the Lombardy Region, Italy. In the week
that followed, Lombardy experienced a very rapid increase in the number of
cases. We analyzed the first 5,830 laboratory-confirmed cases to provide the
first epidemiological characterization of a COVID-19 outbreak in a Western
Country. Epidemiological data were collected through standardized interviews of
confirmed cases and their close contacts. We collected demographic backgrounds,
dates of symptom onset, clinical features, respiratory tract specimen results,
hospitalization, contact tracing. We provide estimates of the reproduction
number and serial interval. The epidemic in Italy started much earlier than
February 20, 2020. At the time of detection of the first COVID-19 case, the
epidemic had already spread in most municipalities of Southern-Lombardy. The
median age for of cases is 69 years (range, 1 month to 101 years). 47% of
positive subjects were hospitalized. Among these, 18% required intensive care.
The mean serial interval is estimated to be 6.6 days (95% CI, 0.7 to 19). We
estimate the basic reproduction number at 3.1 (95% CI, 2.9 to 3.2). We
estimated a decreasing trend in the net reproduction number starting around
February 20, 2020. We did not observe significantly different viral loads in
nasal swabs between symptomatic and asymptomatic. The transmission potential of
COVID-19 is very high and the number of critical cases may become largely
unsustainable for the healthcare system in a very short-time horizon. We
observed a slight decrease of the reproduction number, possibly connected with
an increased population awareness and early effect of interventions. Aggressive
containment strategies are required to control COVID-19 spread and catastrophic
outcomes for the healthcare system.


  Initiation and development of a malignant tumor is a complex phenomenon that
has critical stages determining its long time behavior. This phenomenon is
mathematically described by means of various models: from simple heuristic
models to those employing stochastic processes. In this chapter, we discuss
some aspects of such modeling by analyzing a simple individual-based model, in
which tumor cells are presented as point particles drifting in
$\mathbf{R}_{+}:=[0,+\infty)$ towards the origin with unit speed. At the
origin, each of them splits into two new particles that instantly appear in
$\mathbf{R}_{+}$ at random positions. During their drift the particles are
subject to a random death before splitting. In this model, trait $x\in
\mathbf{R}_{+}$ of a given cell corresponds to time to its division and the
death is caused by therapeutic factors. On its base we demonstrate how to
derive a condition -- involving the therapy related death rate and cell cycle
distribution parameters -- under which the tumor size remains bounded in time,
which practically means combating the disease.


  We present numerical results obtained from the modelling of a stochastic,
highly connected and mobile community. The spread of attributes like health,
disease among the community members is simulated using cellular automata on a
planar 2 dimensional bounded space. With remarkably few assumptions, we are
able to predict the future course of propagation of such disease as a function
of time and the fine-tuning of parameters related to the interaction among the
automata.


  Temporal networks are widely used to represent a vast diversity of systems,
including in particular social interactions, and the spreading processes
unfolding on top of them. The identification of structures playing important
roles in such processes remains largely an open question, despite recent
progresses in the case of static networks. Here, we consider as candidate
structures the recently introduced concept of span-cores: the span-cores
decompose a temporal network into subgraphs of controlled duration and
increasing connectivity, generalizing the core-decomposition of static graphs.
To assess the relevance of such structures, we explore the effectiveness of
strategies aimed either at containing or maximizing the impact of a spread,
based respectively on removing span-cores of high cohesiveness or duration to
decrease the epidemic risk, or on seeding the process from such structures. The
effectiveness of such strategies is assessed in a variety of empirical data
sets and compared to baselines that use only static information on the
centrality of nodes and static concepts of coreness, as well as to a baseline
based on a temporal centrality measure. Our results show that the most stable
and cohesive temporal cores play indeed an important role in epidemic processes
on temporal networks, and that their nodes are likely to represent influential
spreaders.


  In a move described as unprecedented in public health history, starting 24
January 2020, China imposed quarantine and isolation restrictions in Wuhan, a
city of more than 10 million people. This raised the question: is mass
quarantine and isolation effective as a social tool in addition to its
scientific use as a medical tool? In an effort to address this question, using
a epidemiological model driven approach augmented by machine learning, we show
that the quarantine and isolation measures implemented in Wuhan brought down
the effective reproduction number R(t) of the CoVID-19 spread from R(t) > 1 to
R(t) <1 within a month after the imposition of quarantine control measures in
Wuhan, China. This ultimately resulted in a stagnation phase in the infected
case count in Wuhan. Our results indicate that the strict public health
policies implemented in Wuhan may have played a crucial role in halting down
the spread of infection and such measures should potentially be implemented in
other highly affected countries such as South Korea, Italy and Iran to curtail
spread of the disease. Finally, our forecasting results predict a stagnation in
the quarantine control measures implemented in Wuhan towards the end of March
2020; this would lead to a subsequent stagnation in the effective reproduction
number at R(t) <1. We warn that immediate relaxation of the quarantine measures
in Wuhan may lead to a relapse in the infection spread and a subsequent
increase in the effective reproduction number to R(t) >1. Thus, it may be wise
to relax quarantine measures after sufficient time has elapsed, during which
maximum of the quarantined/isolated individuals are recovered.


  The relationship between microscopic observations and macroscopic behavior is
a fundamental open question in biophysical systems. Here, we develop a unified
approach that---in contrast with existing methods---predicts cell type from
macromolecular data even when accounting for the scale of human tissue
diversity and limitations in the available data. We achieve these benefits by
applying a k-nearest-neighbors algorithm after projecting our data onto the
eigenvectors of the correlation matrix inferred from many observations of gene
expression or chromatin conformation. Our approach identifies variations in
epigenotype that impact cell type, thereby supporting the cell type attractor
hypothesis and representing the first step toward model-independent control
strategies in biological systems.


  Attempts to curb the spread of coronavirus by introducing strict quarantine
measures apparently have different effect in different countries: while the
number of new cases has reportedly decreased in China and South Korea, it still
exhibit significant growth in Italy and other countries across Europe. In this
brief note, we endeavour to assess the efficiency of quarantine measures by
means of mathematical modelling. Instead of the classical SIR model, we
introduce a new model of infection progression under the assumption that all
infected individual are isolated after the incubation period in such a way that
they cannot infect other people. Disease progression in this model is
determined by the basic reproduction number $\mathcal{R}_0$ (the number of
newly infected individuals during the incubation period), which is different
compared to that for the standard SIR model. If $\mathcal{R}_0 >1$, then the
number of latently infected individuals exponentially grows. However, if
$\mathcal{R}_0 <1$ (e.g.~due to quarantine measures and contact restrictions
imposed by public authorities), then the number of infected decays
exponentially. We then consider the available data on the disease development
in different countries to show that there are three possible patterns: growth
dynamics, growth-decays dynamics, and patchy dynamics (growth-decay-growth).
Analysis of the data in China and Korea shows that the peak of infection
(maximum of daily cases) is reached about 10 days after the restricting
measures are introduced. During this period of time, the growth rate of the
total number of infected was gradually decreasing. However, the growth rate
remains exponential in Italy. Arguably, it suggests that the introduced
quarantine is not sufficient and stricter measures are needed.


  The goal of this note is to present a simple mathematical model with two
parameters for the number of deaths due to the corona (COVID-19) virus. The
model only requires basic knowledge in differential calculus, and can also be
understood by pupils attending secondary school. The model can easily be
implemented on a computer, and we will illustrate it on the basis of case
studies for different countries.


  We estimate the growth in demand for ICU beds in Chicago during the emerging
COVID-19 epidemic, using state-of-the-art computer simulations calibrated for
the SARS-CoV-2 virus. The questions we address are these:
  (1) Will the ICU capacity in Chicago be exceeded, and if so by how much?
  (2) Can strong mitigation strategies, such as lockdown or shelter in place
order, prevent the overflow of capacity?
  (3) When should such strategies be implemented?
  Our answers are as follows:
  (1) The ICU capacity may be exceeded by a large amount, probably by a factor
of ten.
  (2) Strong mitigation can avert this emergency situation potentially, but
even that will not work if implemented too late.
  (3) If the strong mitigation precedes April 1st, then the growth of COVID-19
can be controlled and the ICU capacity could be adequate. The earlier the
strong mitigation is implemented, the greater the probability that it will be
successful. After around April 1 2020, any strong mitigation will not avert the
emergency situation. In Italy, the lockdown occurred too late and the number of
deaths is still doubling every 2.3 days. It is difficult to be sure about the
precise dates for this window of opportunity, due to the inherent uncertainties
in computer simulation. But there is high confidence in the main conclusion
that it exists and will soon be closed.
  Our conclusion is that, being fully cognizant of the societal trade-offs,
there is a rapidly closing window of opportunity to avert a worst-case scenario
in Chicago, but only with strong mitigation/lockdown implemented in the next
week at the latest. If this window is missed, the epidemic will get worse and
then strong mitigation/lockdown will be required after all, but it will be too
late.


  We characterize the network of COVID-19 spread in India and find that the
transmission rate is 0.43, with daily case growth driven by individuals who
contracted the virus abroad. We explore the question of whether this represents
exponentially decaying dynamics or is simply an artefact of India's testing
strategy. Testing has largely been limited to individuals travelling from
high-risk countries and their immediate contacts, meaning that the network
reflects positive identifications from a biased testing sample. Given generally
low levels of testing and an almost complete absence of testing for community
spread, there is significant risk that we may be missing out on the actual
nature of outbreak. India still has an apparently low current caseload, with
possibly a small window of time to act, and should therefore aggressively and
systematically expand random testing for community spread, including for
asymptomatic cases. This will help understand true transmission characteristics
and plan appropriately for the immediate future.


  Objective: To conduct an early estimation of the Basic Reproduction Number
(BRN) induced by government interference, and to project resulting day to day
number of in-patients, ICU-patients and cumulative number of deaths in a Danish
setting.
  Method: We used the Kermack and McKendrick model with varying basic
reproduction number to estimate number infected and age stratified percentages
to estimate number of in-patients, ICU-patients and cumulative number of
deaths. Changes in basic reproduction number was estimated based on current
in-patient numbers.
  Results: The basic reproductive number in the time period of February 27th to
March 18th was found to be 2.65, however, this number was reduced to 1.99 after
March 18th.
  Keywords: COVID-19, basic reproduction number, Danish population


  Western honeybees (Apis Mellifera) serve extremely important roles in our
ecosystem and economics as they are responsible for pollinating $ 215 billion
dollars annually over the world. Unfortunately, honeybee population and their
colonies have been declined dramatically. The purpose of this article is to
explore how we should model honeybee population with age structure and validate
the model using empirical data so that we can identify different factors that
lead to the survival and healthy of the honeybee colony. Our theoretical study
combined with simulations and data validation suggests that the proper age
structure incorporated in the model and seasonality are important for modeling
honeybee population. Specifically, our work implies that the model assuming
that (1) the adult bees are survived from the {egg population} rather than the
brood population; and (2) seasonality in the queen egg laying rate, give the
better fit than other honeybee models. The related theoretical and numerical
analysis of the most fit model indicate that (a) the survival of honeybee
colonies requires a large queen egg-laying rate and smaller values of the other
life history parameter values in addition to proper initial condition; (b) both
brood and adult bee populations are increasing with respect to the increase in
the {egg-laying rate} and the decreasing in other parameter values; and (c)
seasonality may promote/suppress the survival of the honeybee colony.


  When the Covid-19 pandemic enters dangerous new phase, whether and when to
take aggressive public health interventions to slow down the spread of
COVID-19. To develop the artificial intelligence (AI) inspired methods for
real-time forecasting and evaluating intervention strategies to curb the spread
of Covid-19 in the World. A modified auto-encoder for modeling the transmission
dynamics of the epidemics is developed and applied to the surveillance data of
cumulative and new Covid-19 cases and deaths from WHO, as of March 16, 2020.
The average errors of 5-step forecasting were 2.5%. The total peak number of
cumulative cases and new cases, and the maximum number of cumulative cases in
the world with later intervention (comprehensive public health intervention is
implemented 4 weeks later) could reach 75,249,909, 10,086,085, and 255,392,154,
respectively. The case ending time was January 10, 2021. However, the total
peak number of cumulative cases and new cases and the maximum number of
cumulative cases in the world with one week later intervention were reduced to
951,799, 108,853 and 1,530,276, respectively. Duration time of the Covid-19
spread would be reduced from 356 days to 232 days. The case ending time was
September 8, 2020. We observed that delaying intervention for one month caused
the maximum number of cumulative cases to increase 166.89 times, and the number
of deaths increase from 53,560 to 8,938,725. We will face disastrous
consequences if immediate action to intervene is not taken.


  In late December 2019, a novel strand of Coronavirus (SARS-CoV-2) causing a
severe, potentially fatal respiratory syndrome (COVID-19) was identified in
Wuhan, Hubei Province, China and is causing outbreaks in multiple world
countries, soon becoming a pandemic. Italy has now become the most hit country
outside of Asia: on March 16, 2020, the Italian Civil Protection documented a
total of 27980 confirmed cases and 2158 deaths of people tested positive for
SARS-CoV-2. In the context of an emerging infectious disease outbreak, it is of
paramount importance to predict the trend of the epidemic in order to plan an
effective control strategy and to determine its impact. This paper proposes a
new epidemic model that discriminates between infected individuals depending on
whether they have been diagnosed and on the severity of their symptoms. The
distinction between diagnosed and non-diagnosed is important because
non-diagnosed individuals are more likely to spread the infection than
diagnosed ones, since the latter are typically isolated, and can explain
misperceptions of the case fatality rate and of the seriousness of the epidemic
phenomenon. Being able to predict the amount of patients that will develop
life-threatening symptoms is important since the disease frequently requires
hospitalisation (and even Intensive Care Unit admission) and challenges the
healthcare system capacity. We show how the basic reproduction number can be
redefined in the new framework, thus capturing the potential for epidemic
containment. Simulation results are compared with real data on the COVID-19
epidemic in Italy, to show the validity of the model and compare different
possible predicted scenarios depending on the adopted countermeasures.


  In the advent of the novel coronavirus epidemic since December 2019,
governments and authorities have been struggling to make critical decisions
under high uncertainty at their best efforts. Composite Monte-Carlo (CMC)
simulation is a forecasting method which extrapolates available data which are
broken down from multiple correlated/casual micro-data sources into many
possible future outcomes by drawing random samples from some probability
distributions. For instance, the overall trend and propagation of the infested
cases in China are influenced by the temporal-spatial data of the nearby cities
around the Wuhan city (where the virus is originated from), in terms of the
population density, travel mobility, medical resources such as hospital beds
and the timeliness of quarantine control in each city etc. Hence a CMC is
reliable only up to the closeness of the underlying statistical distribution of
a CMC, that is supposed to represent the behaviour of the future events, and
the correctness of the composite data relationships. In this paper, a case
study of using CMC that is enhanced by deep learning network and fuzzy rule
induction for gaining better stochastic insights about the epidemic development
is experimented. Instead of applying simplistic and uniform assumptions for a
MC which is a common practice, a deep learning-based CMC is used in conjunction
of fuzzy rule induction techniques. As a result, decision makers are benefited
from a better fitted MC outputs complemented by min-max rules that foretell
about the extreme ranges of future possibilities with respect to the epidemic.


  In this note, we describe simple generalizations of the basic SIR model for
epidemic, in case of a multi-region scenario, to be used for predicting the
COVID-19 epidemic spread in Italy.


  In this paper we propose a Susceptible-Infected-Exposed-Recovered-Dead
(SEIRD) differential model for the analysis and forecast of the COVID-19 spread
in some regions of Italy, using the data from the Italian Protezione Civile
from February 24th 2020. In this study investigate an adaptation of the model.
Since several restricting measures have been imposed by the Italian government
at different times, starting from March 8th 2020, we propose a modification of
SEIRD by introducing a time dependent transmitting rate. In the numerical
results we report the maximum infection spread for the three Italian regions
firstly affected by the COVID-19 outbreak(Lombardia, Veneto and Emilia
Romagna). This approach will be successively extended to other Italian regions,
as soon as more data will be available.


  COVID-19 abatement strategies have risks and uncertainties which could lead
to repeating waves of infection. We show -- as proof of concept grounded on
rigorous mathematical evidence -- that periodic, high-frequency alternation of
into, and out-of, lockdown effectively mitigates second-wave effects, while
allowing continued, albeit reduced, economic activity. Periodicity confers (i)
predictability, which is essential for economic sustainability, and (ii)
robustness, since lockdown periods are not activated by uncertain measurements
over short time scales. In turn -- while not eliminating the virus -- this fast
switching policy is sustainable over time, and it mitigates the infection until
a vaccine or treatment becomes available, while alleviating the social costs
associated with long lockdowns. Typically, the policy might be in the form of
1-day of work followed by 6-days of lockdown every week (or perhaps 2 days
working, 5 days off) and it can be modified at a slow-rate based on
measurements filtered over longer time scales. Our results highlight the
potential efficacy of high frequency switching interventions in post lockdown
mitigation. All code is available on Github
(https://github.com/V4p1d/FPSP_Covid19). A software tool has also been
developed so that interested parties can explore the proof-of-concept system.


  In the current COVID19 crisis many national healthcare systems are confronted
with an acute shortage of tests for confirming SARS-CoV-2 infections. For low
overall infection levels in the population, pooling of samples can drastically
amplify the testing efficiency. Here we present a formula to estimate the
optimal pooling size, the efficiency gain (tested persons per test), and the
expected upper bound of missed infections in the pooled testing, all as a
function of the populationwide infection levels and the false negative/positive
rates of the currently used PCR tests. Assuming an infection level of 0.1 % and
a false negative rate of 2 %, the optimal pool size is about 32, the efficiency
gain is about 15 tested persons per test. For an infection level of 1 % the
optimal pool size is 11, the efficiency gain is 5.1 tested persons per test.
For an infection level of 10 % the optimal pool size reduces to about 4, the
efficiency gain is about 1.7 tested persons per test. For infection levels of
30 % and higher there is no more benefit from pooling. To see to what extent
replicates of the pooled tests improve the estimate of the maximal number of
missed infections, we present all results for 1, 3, and 5 replicates.


  Chromatin communities stabilized by protein machinery play essential role in
gene regulation and refine global polymeric folding of the chromatin fiber.
However, treatment of these communities in the framework of the classical
network theory (stochastic block model, SBM) does not take into account
intrinsic linear connectivity of the chromatin loci. Here we propose the
"polymer" block model, paving the way for community detection in polymer
networks. On the basis of this new model we modify the non-backtracking flow
operator and suggest the first protocol for annotation of compartmental domains
in sparse single cell Hi-C matrices. In particular, we prove that our approach
corresponds to the maximum entropy principle. The benchmark analyses
demonstrates that the spectrum of the polymer non-backtracking operator
resolves the true compartmental structure up to the theoretical detectability
threshold, while all commonly used operators fail above it. We test various
operators on real data and conclude that the sizes of the non-backtracking
single cell domains are most close to the sizes of compartments from the
population data. Moreover, the found domains clearly segregate in the gene
density and correlate with the population compartmental mask, corroborating
biological significance of our annotation of the chromatin compartmental
domains in single cells Hi-C matrices.


  We present an empirical algorithm to forecast the evolution of the number of
COVID-19 symptomatic patients in the early stages of the pandemic spread and
after strict social distancing interventions. The algorithm is based on a low
dimensional model for the variation of the exponential growth rate that
decreases after the implementation of strict social distancing measures. From
the observable data given by the number of tested positive, our model estimates
the number of infected hindcast introducing in the model formulation the
incubation time. We also use the model to follow the number of infected
patients who later die using the registered number of deaths and the
distribution time from infection to death. The relationship of the proposed
model with the SIR models is studied. Model parameters fitting is done by
minimizing a quadratic error between the data and the model forecast. An
extended model is also proposed that allows a longer term forecast. An online
implementation of the model is avalaible at www.ctim.es/covid19


  The new coronavirus covid-19 has spread very quickly worldwide, leading the
World Health Organization (WHO) to declare a state of pandemic. One of the
questions many policy makers, and governments are asking themselves is how the
spread is going to evolve in time. In this study, I applied an adapted SIR
model previously used in South Korea to model the MERS outbreak, which is also
caused by a coronavirus, to estimate the evolution of the curve of active cases
in the case of the Portuguese situation, using Italian dara. I then construct
five different scenarios for the evolution of covid-19 in Portugal. In the out
of control scenario, the number of active cases could reach as much as ~40,000
people on the 5th of April (out-of-control scenario). If the self-protective
and control are taken in the same level as what was considered for the South
Korean model, this number could have be reduced to about 800 cases (scenario
1). Considering that this scenario is now unrealistic, three other scenarios
were devised. In all these scenarios, the government measures had a 50\%
effectiveness when compared to the measures in Korea. But in scenario 2 the
transmission rate $\beta$ was effectively reduced to 50\%, In this scenario
active cases could reach circa 7,000 people. In scenario 3, the transmission
rate $\beta$ was reduced to 70\% of its initial value, in which the number of
cases would reach a peak of ~11,000 people. And finally in scenario 4, $\beta$
was reduced to 80\%. In this scenario, the peak would be reached at about
~13,000 cases. This study shows the importance of control and self-protecting
measure to bring down the number of affected people by following the
recommendations of the WHO and health authorities. With the appropriate
measures, this number can be brought down to ~7,000-13,000 people. Hopefully
that will be the case not just in Portugal, but in the rest of the World.


  We analyze the spread of COVID-19 by considering the transmission of the
disease among individuals both within and between regions. A set of regions can
be defined as any partition of a population such that travel/social contact
within each region far exceeds that between them. COVID-19 can be eliminated if
the region-to-region reproductive number---i.e. the average number of other
regions to which a single infected region will transmit the virus---is reduced
to less than one. We find that this region-to-region reproductive number is
proportional to the travel rate between regions and exponential in the length
of the time-delay before region-level control measures are imposed. Thus,
reductions in travel and the speed with which regions take action play decisive
roles in whether COVID-19 is eliminated from a collection of regions. If, on
average, infected regions (including those that become re-infected in the
future) impose social distancing measures shortly after active spreading begins
within them, the number of infected regions, and thus the number of regions in
which such measures are required, will exponentially decrease over time.
Elimination will in this case be a stable fixed point even after the social
distancing measures have been lifted from most of the regions.


  There is a continuing debate on relative benefits of various mitigation and
suppression strategies aimed to control the spread of COVID-19. Here we report
the results of agent-based modelling using a fine-grained computational
simulation of the ongoing COVID-19 pandemic in Australia. This model is
calibrated to match key characteristics of COVID-19 transmission. An important
calibration outcome is the age-dependent fraction of symptomatic cases, with
this fraction for children found to be one-fifth of such fraction for adults.
We apply the model to compare several intervention strategies, including
restrictions on international air travel, case isolation, home quarantine,
social distancing with varying levels of compliance, and school closures.
School closures are not found to bring decisive benefits, unless coupled with
high level of social distancing compliance. We report several trade-offs, and
an important transition across the levels of social distancing compliance, in
the range between 70% and 80% levels, with compliance at the 90% level found to
control the disease within 13--14 weeks, when coupled with effective case
isolation and international travel restrictions.


  In this article, a biallelic reversible mutation model with linear and
quadratic selection is analyzed. The approach reconnects to one proposed by
Kimura ( Possibility of extensive neutral evolution under stabilizing selection
with special reference to nonrandom use of codons (PNAS,1981)), who starts from
a diffusion model and derives its equilibrium distribution up to a constant. We
use a boundary-mutation Moran model, which approximates a general mutation
model for small effective mutation rates, and derive its equilibrium
distribution for polymorphic and monomorphic variants in small to moderately
sized populations. Using this model, we show that biased mutation rates and
linear selection alone can cause patterns of polymorphism rates within and
substitution rates between populations that are usually ascribed to balancing
or overdominant selection. We illustrate this using a data set of short introns
and fourfold degenerate sites from Drosophila simulans and Drosophila
melanogaster.


  Phylogenetic comparative methods correct for shared evolutionary history
among a set of non-independent organisms by modeling sample traits as arising
from a diffusion process along on the branches of a possibly unknown history.
To incorporate such uncertainty, we present a scalable Bayesian inference
framework under a general Gaussian trait evolution model that exploits
Hamiltonian Monte Carlo (HMC). HMC enables efficient sampling of the
constrained model parameters and takes advantage of the tree structure for fast
likelihood and gradient computations, yielding algorithmic complexity linear in
the number of observations. This approach encompasses a wide family of
stochastic processes, including the general Ornstein-Uhlenbeck (OU) process,
with possible missing data and measurement errors. We implement inference tools
for a biologically relevant subset of all these models into the BEAST
phylogenetic software package and develop model comparison through marginal
likelihood estimation. We apply our approach to study the morphological
evolution in the superfamilly of Musteloidea (including weasels and allies) as
well as the heritability of HIV virulence. This second problem furnishes a new
measure of evolutionary heritability that demonstrates its utility through a
targeted simulation study.


  The current outbreak of novel coronavirus disease 2019 (COVID-19) poses an
unprecedented global health and economic threat to interconnected human
societies. Until a vaccine is developed, strategies for controlling the
outbreak rely on aggressive social distancing. These measures largely
disconnect the social network fabric of human societies, especially in urban
areas. Here, we estimate the growth rates and reproductive numbers of COVID-19
in US cities from March 14th through March 19th to reveal a power-law scaling
relationship to city population size. This means that COVID-19 is spreading
faster on average in larger cities with the additional implication that, in an
uncontrolled outbreak, larger fractions of the population are expected to
become infected in more populous urban areas. We discuss the implications of
these observations for controlling the COVID-19 outbreak, emphasizing the need
to implement more aggressive distancing policies in larger cities while also
preserving socioeconomic activity.


  For nearly a century, the initial reproduction number (R0) has been used as a
one number summary to compare outbreaks of infectious disease, yet there is no
`standard' estimator for R0. Difficulties in estimating R0 arise both from how
a disease transmits through a population as well as from differences in
statistical estimation method. We describe eight methods used to estimate R0
and provide a thorough simulation study of how these estimates change in the
presence of different disease parameters. As motivation, we analyze the 2009
outbreak of the H1N1 pandemic influenza in the USA and compare the results from
our eight methods to a previous study. We discuss the most important aspects
from our results which effect the estimation of R0, which include the
population size, time period used, and the initial percent of infectious
individuals. Additionally, we discuss how pre-processing incidence counts may
effect estimates of R0. Finally, we provide guidelines for estimating point
estimates and confidence intervals to create reliable, comparable estimates of
R0.


  Colonies of rod-shaped bacteria constitute a system of colloidal active
matter with nematic properties. As a single initial bacterium multiplies
through repeated divisions, the resulting colony quickly loses long-range
orientational order, but retains locally ordered domains. At the boundaries of
these domains, topological defects emerge, which move around randomly as the
colony grows. In both experiments and simulations, we find that these defects
are created at a rate that corresponds to the exponential growth of the colony,
resulting in a stable defect density. Both this defect density and the colony's
correlation length are regulated by the aspect ratio of the rod-shaped
particles. Moreover, we find that the defect dynamics are well described by a
Gamma distribution, which is due to repeated divisions and subsequent
re-orientations of the bacteria.


  The data analysis on deaths in the Lombardy Region and of both the provinces
of Bergamo and Brescia shows a twofold aspect on the trend of the epidemic: -
all the data show a bias linked to the event of March 10th (day for which the
Lombardy region data is partial) and the subsequent change in the way in which
positive cases and deaths are calculated; - following the containment measures
of the Prime Minister's Decree of March 11th, the spread of the epidemic,
although still exponential in nature, has a reduced multiplication coefficient.
Our analysis concludes that the situation is not yet compatible with a plateau
trend and allows us to predict the trend in the number of deaths in the
Lombardy region. We therefore conclude that the containment measures put in
place by the government on March 11th will allow a reduction in deaths from
around 8000 to just over 6500 for March 27th.


  We propose a simple model for understanding the kinetics of corona virus
transmission. Our model assume spreading of corona virus can happen from one to
another only, if someone without enough protection comes close contact to a
person carrying the corona virus. Therefore this virus spreads on a large scale
within a short time through chains of such events. Using our model we provide
an estimation of the number of people affected by this virus within reasonable
duration of time. We choose values of different parameters of our model by
non-linear least square fit of the real time data and we predict fate of this
corona virus transmission using our model.


  We develop a simple 3-dimensional iterative map model to forecast the global
spread of the coronavirus disease. Our model contains at most two fitting
parameters, which we determine from the data supplied by the world health
organisation for the total number of cases and new cases each day. We find that
our model provides a surprisingly good fit to the currently-available data,
which exhibits a cross-over from exponential to power-law growth, as lock-down
measures begin to take effect. Before these measures, our model predicts
exponential growth from day 30 to 69, starting from the date on which the world
health organisation provided the first `Situation report' (21 January 2020 $-$
day 1). Based on this initial data the disease may be expected to infect
approximately 23% of the global population, i.e. about 1.76 billion people,
taking approximately 83 million lives. Under this scenario, the global number
of new cases is predicted to peak on day 133 (about the middle of May 2020),
with an estimated 60 million new cases per day. If current lock-down measures
can be maintained, our model predicts power law growth from day 69 onward. Such
growth is comparatively slow and would have to continue for several decades
before a sufficient number of people (at least 23% of the global population)
have developed immunity to the disease through being infected. Lock-down
measures appear to be very effective in postponing the unimaginably large peak
in the daily number of new cases that would occur in the absence of any
interventions. However, should these measure be relaxed, the spread of the
disease will most likely revert back to its original exponential growth
pattern. As such, the duration and severity of the lock-down measures should be
carefully timed against their potentially devastating impact on the world
economy.


  There is a rich history of models for the interaction of a biological
contagion like influenza with the spread of related information such as an
influenza vaccination campaign. Recent work on the spread of interacting
contagions on networks has highlighted that these interacting contagions can
have counter-intuitive interplay with network structure. Here we generalize one
of these frameworks to tackle three important features of the spread of
awareness and disease: one, we model the dynamics on highly clustered,
cliquish, networks to mimic the role of workplaces and households; two, the
awareness contagion affects the spread of the biological contagion by reducing
its transmission rate where an aware or vaccinated individual is less likely to
be infected; and three, the biological contagion also affects the spread of the
awareness contagion but by increasing its transmission rate where an infected
individual is more receptive and more likely to share information related to
the disease. Under these conditions, we find that increasing network
clustering, which is known to hinder disease spread, can actually allow them to
sustain larger epidemics of the disease in models with awareness. This
counter-intuitive result goes against the conventional wisdom suggesting that
random networks are justifiable as they provide worst-case scenario forecasts.
To further investigate this result, we provide a closed-form criterion based on
a two-step branching process (i.e., the numbers of expected tertiary
infections) to identify different regions in parameter space where the net
effect of clustering and co-infection varies. Altogether, our results highlight
once again the need to go beyond random networks in disease modeling and
illustrate the type of analysis that is possible even in complex models of
interacting contagions.


  Traditional vaccines against viruses are designed to target their surface
proteins, i.e., antigens, which can trigger the immune system to produce
specific antibodies to capture and neutralize the viruses. However, viruses
often evolve quickly, and their antigens are prone to mutations to avoid
recognition by the antibodies (antigenic drift). This limitation of the
antibody-mediated immunity could be addressed by the T-cell mediated immunity,
which is able to recognize conserved viral HLA peptides presented on
virus-infected cells. Thus, by targeting conserved regions on the genome of a
virus, T-cell epitope-based vaccines are less subjected to mutations and may
work effectively on different strains of the virus. Here we propose a
personalized workflow to identify an optimal set of T-cell epitopes based on
the HLA alleles and the immunopeptidome of an individual person. Specifically,
our workflow trains a machine learning model on the immunopeptidome and then
predicts HLA peptides from conserved regions of a virus that are most likely to
trigger responses from the person T cells. We applied the workflow to identify
T-cell epitopes for the SARS-COV-2 virus, which has caused the recent COVID-19
pandemic in more than 100 countries across the globe.


  The ongoing pandemic of Coronavirus disease (COVID-19) emerged in Wuhan,
China in the end of 2019. It has already affected more than 300,000 people,
with the number of deaths nearing 13000 across the world. As it has been posing
a huge threat to global public health, it is of utmost importance to identify
the rate at which the disease is spreading. In this study, we propose a time
series model to analyze the trend pattern of the incidence of COVID-19
outbreak. We also incorporate information on total or partial lockdown,
wherever available, into the model. The model is concise in structure, and
using appropriate diagnostic measures, we showed that a time-dependent
quadratic trend successfully captures the incidence pattern of the disease. We
also estimate the basic reproduction number across different countries, and
find that it is consistent except for the United States of America. The above
statistical analysis is able to shed light on understanding the trends of the
outbreak, and gives insight on what epidemiological stage a region is in. This
has the potential to help in prompting policies to address COVID-19 pandemic in
different countries.


  The first cases of COVID-19 in France were detected on January 24, 2020. The
number of screening tests carried out and the methodology used to target the
patients tested do not allow for a direct computation of the real number of
cases and the mortality rate.In this report, we develop a
'mechanistic-statistical' approach coupling a SIR ODE model describing the
unobserved epidemiological dynamics, a probabilistic model describing the data
acquisition process and a statistical inference method. The objective of this
model is not to make forecasts but to estimate the real number of people
infected with COVID-19 during the observation window in France and to deduce
the mortality rate associated with the epidemic.Main results. The actual number
of infected cases in France is probably much higher than the observations: we
find here a factor x 15 (95%-CI: 4-33), which leads to a 5.2/1000 mortality
rate (95%-CI: 1.5 / 1000-11.7/ 1000) at the end of the observation period. We
find a R0 of 4.8, a high value which may be linked to the long viral shedding
period of 20 days.


  Epidemic is a rapid and wide spread of infectious disease threatening many
lives and economy damages. It is important to fore-tell the epidemic lifetime
so to decide on timely and remedic actions. These measures include closing
borders, schools, suspending community services and commuters. Resuming such
curfews depends on the momentum of the outbreak and its rate of decay. Being
able to accurately forecast the fate of an epidemic is an extremely important
but difficult task. Due to limited knowledge of the novel disease, the high
uncertainty involved and the complex societal-political factors that influence
the widespread of the new virus, any forecast is anything but reliable. Another
factor is the insufficient amount of available data. Data samples are often
scarce when an epidemic just started. With only few training samples on hand,
finding a forecasting model which offers forecast at the best efforts is a big
challenge in machine learning. In the past, three popular methods have been
proposed, they include 1) augmenting the existing little data, 2) using a panel
selection to pick the best forecasting model from several models, and 3)
fine-tuning the parameters of an individual forecastingmodel for the highest
possible accuracy. In this paper, a methodology that embraces these three
virtues of data mining from a small dataset is proposed...


  During the current COVID-19 pandemic it is imperative to give early warnings
to reduce mortality. However, non-specialist such as authorities and the
general population face several problems to understand the real thread of this
pandemic, and under/overestimation of its risk are a commonplace in the press
and social media. Here we define an index, which we call the COVID-19 Burden
Index, that relates the capacities of the healthcare system of a given country
to treat severe and critical cases. Its value is 0 if there is no extra strain
in the healthcare system, and it reaches 1.0 when the collapse is imminent. As
of 23 March 2020, we show that Chile, the USA, UK, among other countries, must
reduce the rate of infections right now, otherwise, in less than 7 days they
could be in a catastrophic situation such as Italy, Spain and Iran.


  Widespread opinions and discussion exist regarding the efficiency of social
distancing after crucial spread of the SARS-CoV-2 virus during the actual
Covid-19 pandemic. While Germany has released a federal law that prohibits any
type of direct contact for more than 2 people other countries including the US
released curfews.
  People are now wondering whether these measures are helpful to stop or hamper
the Covid-19 pandemic and to limit the spread of the new corona virus.
  A quantitative statement on this question depends on many parameters that are
difficult to grasp mathematically and cannot therefore be made conclusively
(they include consistent adherence to the measures decided, the estimated
number of unreported cases, the possible limitation by test capacities,
possible mutations of the virus, etc ...). However, it turns out that a
reduction in the actual daily new infection rate (actual daily growth rate of
reported cases, in short: infection rate) from the current value of 30-35% in
the US to 10% would be extremely effective in stopping the spread of the virus.
The severe restrictions in Germany which closed any public events, schools and
universities a week ago might already have contributed to a reduction of the
growth rate of reported cases below 30%.


  Cytoskeletons are self organized networks based on polymerized proteins,
actin, tubulin, and driven by motor proteins, such as myosin, kinesin and
dynein. Their positive Darwinian evolution enables them to approach optimized
functionality self organized criticality. Dynein has three distinct titled
subunits, but how these units connect to function as a molecular motor is
mysterious. Dynein binds to tubulin through two coiled coil stalks and a stalk
head. The energy used to alter the head binding and propel cargo along tubulin
is supplied by ATP at a ring 1500 amino acids away.


  Studies on the influence of a modern lifestyle in abetting Coronary Heart
Diseases (CHD) have mostly focused on deterrent health factors, like smoking,
alcohol intake, cheese consumption and average systolic blood pressure, largely
disregarding the impact of a healthy lifestyle in mitigating CHD risk. In this
study, 30+ years' World Health Organization (WHO) data have been analyzed,
using a wide array of advanced Machine Learning techniques, to quantify how
regulated reliance on positive health indicators, e.g. fruits/vegetables,
cereals can offset CHD risk factors over a period of time. Our research ranks
the impact of the negative outliers on CHD and then quantifies the impact of
the positive health factors in mitigating the negative risk-factors. Our
research outcomes, presented through simple mathematical equations, outline the
best CHD prevention strategy using lifestyle control only. We show that a 20%
increase in the intake of fruit/vegetable leads to 3-6% decrease in SBP; or, a
10% increase in cereal intake lowers SBP by 3%; a simultaneous increase of 10%
in fruit-vegetable can further offset the effects of SBP by 6%. Our analysis
establishes gender independence of lifestyle on CHD, refuting long held
assumptions and unqualified beliefs. We show that CHD risk can be lowered with
incremental changes in lifestyle and diet, e.g. fruit-vegetable intake
ameliorating effects of alcohol-smoking-fatty food. Our multivariate data model
also estimates functional relationships amongst lifestyle factors that can
potentially redefine the diagnostics of Framingham score-based CHD-prediction.


  The emerging global infectious COVID-19 coronavirus disease by novel Severe
Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) presents critical threats
to global public health and the economy since it was identified in late
December 2019 in China. The virus has gone through various pathways of
evolution. For understanding the evolution and transmission of SARS-CoV-2,
genotyping of virus isolates is of great importance. We present an accurate
method for effectively genotyping SARS-CoV-2 viruses using complete genomes.
The method employs the multiple sequence alignments of the genome isolates with
the SARS-CoV-2 reference genome. The SNP genotypes are then measured by Jaccard
distances to track the relationship of virus isolates. The genotyping analysis
of SARS-CoV-2 isolates from the globe reveals that specific multiple mutations
are the predominated mutation type during the current epidemic. Our method
serves a promising tool for monitoring and tracking the epidemic of pathogenic
viruses in their gradual and local genetic variations. The genotyping analysis
shows that the genes encoding the S proteins and RNA polymerase, RNA primase,
and nucleoprotein, undergo frequent mutations. These mutations are critical for
vaccine development in disease control.


  Do people discriminate between men and women when they have the option to
punish defectors or reward cooperators? Here we report on four pre-registered
experiments, that shed some light on this question. Study 1 (N=544) shows that
people do not discriminate between genders when they have the option to punish
(reward) defectors (cooperators) in a one-shot prisoner's dilemma with
third-party punishment/reward. Study 2 (N=253) extends Study 1 to a different
method of punishing/rewarding: participants are asked to rate the behaviour of
a defector/cooperator on a scale of 1 to 5 stars. In this case too, we find
that people do not discriminate between genders. Study 3a (N=331) and Study 3b
(N=310) conceptually replicate Study 2 with a slightly different gender
manipulation. These latter studies show that, in situations where they do not
have specific beliefs about the gender of the defector/cooperator's partner,
neither men nor women discriminate between genders.


  In the wake of the COVID-19 virus outbreak, a brief phenomenological
(descriptive, comparative) analysis of the dynamics of the disease spread among
different countries is presented. Results indicate that the infection spread
dynamics is currently the most pronounced in the USA (confirmed cases are
currently doubling every 2.16 days, with a decreasing doubling time tendency),
while other countries with the most confirmed cases show different values, and
tendencies. The reported number of deaths is currently doubled every 2.28 days
in Germany, 2.56 days in France, 2.57 days in Switzerland, 2.59 days in France,
and 2.62 days in USA, while only France and USA are currently exhibiting
further acceleration (diminishing doubling time).


  We propose a novel numerical method able to determine efficiently and
effectively the relationship of complementarity between portions of proteins
surfaces. This innovative and general procedure, based on the representation of
the molecular iso-electron density surface in terms of 2D Zernike polynomials,
allows the rapid and quantitative assessment of the geometrical shape
complementarity between interacting proteins, that was unfeasible with previous
methods. We first tested the method with a large dataset of known protein
complexes obtaining an overall area under the ROC curve of 0.76 in the blind
recognition of binding sites and then applied it to investigate the features of
the interaction between the Spike protein of SARS-Cov-2 and human cellular
receptors. Our results indicate that SARS-CoV-2 uses a dual strategy: its spike
protein could also interact with sialic acid receptors of the cells in the
upper airways, in addition to the known interaction with Angiotensin-converting
enzyme 2.


  In this investigation I used the Logistic Model to fit the COVID-19 pandemic
data for some countries. The data modeled is the death numbers per day in
China, Iran, Italy, South Korea, Spain and United States. Considering the
current growth rate of the pandemic, it was possible to show that the death
toll should be between 3,277-3,327 deaths in China, 2,035-2,107 in Iran,
120-134 in South Korea, 11,227-12,793 in Italy and 6,217-7,405 in Spain. Also,
with this toy model it was possible to show a clear emergence of a new outbreak
within the same country (Iran, China and the United States). The growth rate of
deaths found for South Korea was the lowest among the countries studied
($0.14701\pm0.00923$) and for China ($0.16667\pm0.00284$). Italy
($0.22594\pm0.00599$) and Spain ($0.31213\pm0.02337$) had the highest rates and
in the second wave in Iran ($0.37893\pm0.02712$).


  The 2019 Novel Corona virus infection (COVID 19) is an ongoing public health
emergency of international focus. Significant gaps persist in our knowledge of
COVID 19 epidemiology, transmission dynamics, investigation tools and
management, despite (or possibly because of) the fact that the outbreak is an
unprecedented global threat. On the positive side, enough is currently known
about the epidemic process to permit the construction of mathematical
predictive models. In our work, we adapt a traditional SEIR epidemic model to
the specific dynamic compartments and epidemic parameters of COVID 19, as it
spreads in an age-heterogeneous community. We analyze management strategies of
the epidemic course (as they were implemented through lockdown and reopening
procedures in many of the US states and countries worldwide); however, to more
clearly illustrate ideas, we focus on the example of a small scale college town
community, with the timeline of control measures introduced in the state of New
York. We generate predictions, and assess the efficiency of these control
measures (closures, mobility restrictions, social distancing), in a
sustainability context.


  I estimate the Susceptible-Infected-Recovered (SIR) epidemic model for
Coronavirus Disease 2019 (COVID-19). The transmission rate is heterogeneous
across countries and far exceeds the recovery rate, which enables a fast
spread. In the benchmark model, 28% of the population may be simultaneously
infected at the peak, potentially overwhelming the healthcare system. The peak
reduces to 6.2% under the optimal mitigation policy that controls the timing
and intensity of social distancing. A stylized asset pricing model suggests
that the stock price temporarily decreases by 50% in the benchmark case but
shows a W-shaped, moderate but longer bear market under the optimal policy.


  In this paper we develop a SIR epidemiological model with parameters
calculated according to existing data at the time of writing (24/03/2020); the
data is from Italy, South Korea and Colombia, the model is then used to project
the evolution of the COVID-19 epidemic in Colombia for different scenarios
using the data of population for the country and known initial conditions at
the start of the simulation.


  COVID-19 is a new type of coronavirus disease which is caused by severe acute
respiratory syndrome coronavirus 2 (SARS-CoV-2). It originated in China in the
month of December 2019 and quickly started to spread within the country. On
31st December 2019, it was first reported to country office of World Health
Organization (WHO) in China. Since then, it has spread to most of the countries
around the globe. However, there has been a recent rise in trend in believing
that it would go away during summer days, which has not yet been properly
investigated. In this paper, relationship of daily number of confirmed cases of
COVID-19 with three environmental factors, viz. maximum relative humidity
(RH_max), maximum temperature (T_max) and highest wind speed (WS_max),
considering the incubation period, have been investigated statistically, for
four of the most affected places of China, viz. Beijing, Chongqing, Shanghai,
Wuhan and five of the most affected places of Italy, viz. Bergamo, Cremona,
Lodi, Milano. It has been found that the relationship with maximum relative
humidity and highest wind is mostly negligible, whereas relationship with
maximum temperature is ranging between negligible to moderate.


  We give some numerical observations on the total number of infected by the
SARS-CoV-2 in Italy. The analysis is based on a tanh formula involving two
parameters. A polynomial correlation between the parameters gives an upper
bound for the time of the peak of new infected. A numerical indicator of the
temporal variability of the upper bound is introduced. The result and the
possibility to extend the analysis to other countries are discussed in the
conclusions.


  In this research, we study the propagation patterns of epidemic diseases such
as the COVID-19 coronavirus, from a mathematical modeling perspective. The
study is based on an extensions of the well-known
susceptible-infected-recovered (SIR) family of compartmental models. It is
shown how social measures such as distancing, regional lockdowns, quarantine
and global public health vigilance, influence the model parameters, which can
eventually change the mortality rates and active contaminated cases over time,
in the real world. As with all mathematical models, the predictive ability of
the model is limited by the accuracy of the available data and to the so-called
\textit{level of abstraction} used for modeling the problem. In order to
provide the broader audience of researchers a better understanding of spreading
patterns of epidemic diseases, a short introduction on biological systems
modeling is also presented and the Matlab source codes for the simulations are
provided online.


  We report a statistical analysis of some highly infected countries by the
novel coronavirus (COVID-19). The cumulative infected data were fitted with
various growth models (e.g. Logistic equation, Weibull equation and Hill
equation) and obtained the power index of top ten highly infected countries.
The newly infected data were fitted with Gaussian distribution with the peak at
~40 days for the countries whose infection curves are seem to be saturated. The
similarity in growth kinetics of infected people of different countries
provides first-hand guidelines to take proper precautions to minimize human
damage.


  The rapid increase of transportation infrastructure during the recent decades
has caused a number of effects on bird species, including collision mortality,
habitat loss, fragmentation and noise. This paper investigates the effects of
traffic noise and road-induced fragmentation on breeding bird richness in
Cyprus. Cyprus, situated along one of the main migratory routes for birds, has
a rich and diverse avifauna threatened by an ever-expanding road network and a
road density among the highest in Europe. In this first island-wide study we
used data from 102 breeding birds recorded in 10 km x 10 km grid cells. Within
every cell we calculated road traffic noise and eight road-related properties.
Most of the grid cells are subject to intense fragmentation and traffic noise
with combined impact hotspots located even within protected areas (such as Cape
Greco, and the Troodos Massif). Results from variance partitioning indicated
that road-related properties (total road extent and road length) accounted for
a combined 59% of variation in species richness, followed by
fragmentation-related properties and noise properties. The study posits the
need for further in-depth research on the effects of road networks on birds,
and road construction, particularly in protected areas within Mediterranean
islands.


  The new coronavirus (2019-nCoV or SARS-CoV2), inducing the current pandemic
disease (COVID-19) and causing pneumoniae in humans, is dramatically increasing
in epidemic scale since its first appearance in Wuhan, China, in December 2019.
The first infection from epidemic coronaviruses in 2003 fostered the spread of
an overwhelming amount of related scientific efforts. The manifold aspects that
have been raised, as well as their redundancy offer precious information that
has been underexploited and needs to be critically re-evaluated, appropriately
used and offered to the whole community, from scientists, to medical doctors,
stakeholders and common people. These efforts will favour a holistic view on
the comprehension, prevention and development of strategies (pharmacological,
clinical etc) as well as common intervention against the new coronavirus
spreading. Here we describe a model that emerged from our analysis that was
focused on the Renin Angiotensin System (RAS) and the possible routes linking
it to the viral infection. because the infection is mediated by the viral
receptor on human cell membranes Angiotensin Converting Enzyme (ACE2), which is
a key component in RAS signalling. The model depicts the main pathways
determining the disease and the molecular framework for its establishment, and
can help to shed light on mechanisms involved in the infection. It promptly
gives an answer to some of the controversial, and still open, issues concerning
predisposing conditions and medical treatments that protect from or favour the
severity of the disease (such as the use of ACE inhibitors or ARBs/sartans), or
to the sex related biases in the affected population. The model highlights
novel opportunities for further investigations, diagnosis and appropriate
intervention to understand and fight COVID19.


  We present Coronavirus disease 2019 (COVID-19) statistics in China dataset:
daily statistics of the COVID-19 outbreak in China at the city/county level.
For each city/country, we include the six most important numbers for epidemic
research: daily new infections, accumulated infections, daily new recoveries,
accumulated recoveries, daily new deaths, and accumulated deaths. We cross
validate the dataset and the estimate error rate is about 0.04%. We then give
several examples to show how to trace the spreading in particular cities or
provinces, and also contrast the development of COVID-19 in all cities in China
at the early, middle and late stages. We hope this dataset can help researchers
around the world better understand the spreading dynamics of COVID-19 at a
regional level, to inform intervention and mitigation strategies for
policymakers.


  The third-generation long reads sequencing technologies, such as PacBio and
Nanopore, have great advantages over second-generation Illumina sequencing in
de novo assembly studies. However, due to the inherent low base accuracy,
third-generation sequencing data cannot be used for k-mer counting and
estimating genomic profile based on k-mer frequencies. Thus, in current genome
projects, second-generation data is also necessary for accurately determining
genome size and other genomic characteristics. We show that corrected
third-generation data can be used to count k-mer frequencies and estimate
genome size reliably, in replacement of using second-generation data.
Therefore, future genome projects can depend on only one sequencing technology
to finish both assembly and k-mer analysis, which will largely decrease
sequencing cost in both time and money. Moreover, we present a fast
light-weight tool kmerfreq and use it to perform all the k-mer counting tasks
in this work. We have demonstrated that corrected third-generation sequencing
data can be used to estimate genome size and developed a new open-source C/C++
k-mer counting tool, kmerfreq, which is freely available at
https://github.com/fanagislab/kmerfreq.


  We revisit well-established concepts of epidemiology, the Ising-model, and
percolation theory. Also, we employ a spin $S$ = 1/2 Ising-like model and a
(logistic) Fermi-Dirac-like function to describe the spread of Covid-19. Our
analysis reinforces well-established literature results, namely: \emph{i}) that
the epidemic curves can be described by a Gaussian-type function; \emph{ii})
that the temporal evolution of the accumulative number of infections and
fatalities follow a logistic function, which has some resemblance with a
distorted Fermi-Dirac-like function; \emph{iii}) the key role played by the
quarantine to block the spread of Covid-19 in terms of an \emph{interacting}
parameter, which emulates the contact between infected and non-infected people.
Furthermore, in the frame of elementary percolation theory, we show that:
\emph{i}) the percolation probability can be associated with the probability of
a person being infected with Covid-19; \emph{ii}) the concepts of blocked and
non-blocked connections can be associated, respectively, with a person
respecting or not the social distancing, impacting thus in the probability of
an infected person to infect other people. Increasing the number of infected
people leads to an increase in the number of net connections, giving rise thus
to a higher probability of new infections (percolation). We demonstrate the
importance of social distancing in preventing the spread of Covid-19 in a
pedagogical way. Given the impossibility of making a precise forecast of the
disease spread, we highlight the importance of taking into account additional
factors, such as climate changes and urbanization, in the mathematical
description of epidemics. Yet, we make a connection between the standard
mathematical models employed in epidemics and well-established concepts in
condensed matter Physics, such as the Fermi gas and the Landau Fermi-liquid
picture.


  In this short note we model the region-wise trends of the evolution to
COVID-19 infections using a stochastic SIR model. The SIR dynamics are
expressed using \textit{It\^o-stochastic differential equations}. We first
derive the parameters of the model from the available daily data from European
regions based on a 24-day history of infections, recoveries and deaths. The
derived parameters have been aggregated to project future trends for the Indian
subcontinent, which is currently at an early stage in the infection cycle. The
projections are meant to serve as a guideline for strategizing the
socio-political counter measures to mitigate COVID-19.


  The CoVid-19 is spreading pandemically all over the world. A rapid defeat of
the pandemic requires carrying out on the population a mass screening, able to
separate positive from negative cases. Such a cleaning will free a flow of
productive population. The current rate and cost of testing, performed with the
common PCR (polymerase chain reaction) method and with the available resources,
is forcing a selection of the subjects to be tested. Indeed, each one must be
examined individually at the cost of precious time. Moreover, the exclusion of
potentially positive individuals from screening induces health risks, a broad
slowdown in the effort to curb the viral spread, and the consequent mortality
rates. We present a new procedure, the Purified by Unified Resampling of
Infected Multitudes, in short Purim, able to untangle any massive candidate
sample with inexpensive screening, through the cross-correlated analysis of the
joint speciments. This procedure can reveal and detect most negative patients
and in most cases discover the identity of the few positives already in the
first or few secondary tests. We investigate the the two-dimensional
correlation case in function of the infection probability. The
multi-dimensional topology, the scaled Purim procedure are also considered.
Extensive Purim tests may measure and weight the degree of epidemic: their
outcome may identify focal regions in the early stages. Assuming hundreds or
thousand subjects, the saving both in time and in cost will be remarkable.
Purim may be able to filter scheduled flights, scholar acceptance, popular
international event participants. The optimal extension of Purim outcome is
growing as the inverse of the epidemia expansion. Therefore, the earlier, the
better.


  The number of corona virus (COVID-19) infections grows worldwide. In order to
create short term predictions to prepare for the extent of the global pandemic
we analyze infection data from the top 25 affected countries. It is shown that
all country-specific infection rates follow a power law growth behavior and the
scaling exponents per country are calculated. We find two different growth
patterns: either steady power law growth from the very beginning with moderate
scaling exponents of 3-5 or explosive power law growth with dramatic scaling
exponents of 8-11. In the case of the USA we even find an exponent of 16.59. By
means of data analysis we confirm that instituting strict measures of
lock-downs combined with a strong adherence by the population are effective
means to bring the growth rates down. While many countries have instituted
measures there are only three countries (Denmark, Norway, and South Korea) so
far where such lock-downs led to a significant reduction of the growth rate. In
the case of Denmark we calculate the reduction of the scaling exponents to move
from 6.82 to 1.47.


  In this paper, we have applied the logistic growth regression model and
genetic algorithm to predict the number of coronavirus infected cases that can
be expected in upcoming days in India and also estimated the final size and its
peak time of the coronavirus epidemic in India.


  An urgent problem in controlling COVID-19 spreading is to understand the role
of undocumented infection. We develop a five-state model for COVID-19, taking
into account the unique features of the novel coronavirus, with key parameters
determined by the government reports and mathematical optimization. Tests using
data from China, South Korea, Italy, and Iran indicate that the model is
capable of generating accurate prediction of the daily accumulated number of
confirmed cases and is entirely suitable for real-time prediction. The
drastically disparate testing and diagnostic standards/policies among different
countries lead to large variations in the estimated parameter values such as
the duration of the outbreak, but such uncertainties have little effect on the
occurrence time of the inflection point as predicted by the model, indicating
its reliability and robustness. Model prediction for Italy suggests that
insufficient government action leading to a large fraction of undocumented
infection plays an important role in the abnormally high mortality in that
country. With the data currently available from United Kingdom, our model
predicts catastrophic epidemic scenarios in the country if the government did
not impose strict travel and social distancing restrictions. A key finding is
that, if the percentage of undocumented infection exceeds a threshold, a
non-negligible hidden population can exist even after the the epidemic has been
deemed over, implying the likelihood of future outbreaks should the currently
imposed strict government actions be relaxed. This could make COVID-19 evolving
into a long-term epidemic or a community disease a real possibility, suggesting
the necessity to conduct universal testing and monitoring to identify the
hidden individuals.


  Different ways of calculating mortality ratios during epidemics have yielded
very different results, particularly during the current COVID-19 pandemic. We
formulate both a survival probability model and an associated infection
duration-dependent SIR model to define individual- and population-based
estimates of dynamic mortality ratios. The key parameters that affect the
dynamics of the different mortality estimates are the incubation period and the
time individuals were infected before confirmation of infection. We stress that
none of these ratios are accurately represented by the often misinterpreted
case fatality ratio (CFR), the number of deaths to date divided by the total
number of confirmed infected cases to date. Using data on the recent SARS-CoV-2
outbreaks, we estimate and compare the different dynamic mortality ratios and
highlight their differences. Informed by our modeling, we propose more
systematic methods to determine mortality ratios during epidemic outbreaks and
discuss sensitivity to confounding effects and uncertainties in the data.


  The outbreak of the novel coronavirus, COVID-19, has been declared a pandemic
by the WHO. The structures of social contact critically determine the spread of
the infection and, in the absence of vaccines, the control of these structures
through large-scale social distancing measures appears to be the most effective
means of mitigation. Here we use an age-structured SIR model with social
contact matrices obtained from surveys and Bayesian imputation to study the
progress of the COVID-19 epidemic in India. The basic reproductive ratio R0 and
its time-dependent generalization are computed based on case data, age
distribution and social contact structure. The impact of social distancing
measures - workplace non-attendance, school closure, lockdown - and their
efficacy with durations are then investigated. A three-week lockdown is found
insufficient to prevent a resurgence and, instead, protocols of sustained
lockdown with periodic relaxation are suggested. Forecasts are provided for the
reduction in age-structured morbidity and mortality as a result of these
measures. Our study underlines the importance of age and social contact
structures in assessing the country-specific impact of mitigatory social
distancing.


  From March 23rd, the data for the recovered cases of COVID-19 are missing
from the standard repository maintained by the Johns Hopkins University in
collaboration with the WHO. But since data concerning recovered patients are
extremely important for modelling the COVID-19 outbreak, a method for
estimating the missing data is provided and tested. As a byproduct, it produces
estimates for the probabilities to die $k$ days after confirmation, or to
survive after $d$ days.


  Honeybees play an important role in the production of many agricultural crops
and in sustaining plant diversity in undisturbed ecosystems. The rapid decline
of honeybee populations have sparked great concern worldwide. Previous studies
have shown that the parasitic Varroa mite could be the main reason for colony
losses. In order to understand how mites affect population dynamics of
honeybees and a colony health, we propose a brood-adult bee-mite model in which
the time lag from brood to adult is taken into account. Noting that the
dynamics of a honeybee colony varies with respect to season, we validate the
model and perform parameter estimations under both constant and fluctuating
seasonality scenarios. Our analytical and numerical studies reveal the
following: (a) In the presence of parasite mites, the large time lag from brood
to adult could destabilize population dynamics and drive the colony to
collapse; but the small natural mortality of the adult population can promote a
mite-free colony when time lag is small or at an intermediate level; (b) Small
brood' infestation rates could stabilize all populations at the unique interior
equilibrium under constant seasonality while may drive the mite population to
die out when seasonality is taken into account; (c) High brood' infestation
rates can destabilize the colony dynamics leading to population collapse
depending on initial population size under constant and seasonal conditions;
(d) Results from sensitivity analysis indicate the queen's egg-laying may have
the greatest effect on colony population size. The brood death rate and the
colony size at which brood survivability is the half maximal were also shown to
be highly sensitive with an inverse correlation to the colony population size.
Our results provide insights on the effects of seasonality on the dynamics.


  The world evolution of the Severe acute respiratory syndrome coronavirus 2
(SARS-Cov2 or simply COVID-19) led the World Health Organization to declare it
a pandemic. The disease appeared in China in December 2019, and it has spread
fast around the world, specially in european countries like Italy and Spain.
The first reported case in Brazil was recorded in February 26, and after that
the number of cases growed fast. In order to slow down the initial growth of
the disease through the country, confirmed positive cases were isolated to not
transmit the disease. To better understand the early evolution of COVID-19 in
Brazil, we apply a Susceptible-Infectious-Quarantined-Recovered (SIQR) model to
the analysis of data from the Brazilian Department of Health, obtained from
February 26, 2020 through March 25, 2020. Based on analyical and numerical
results, as well on the data, the basic reproduction number is estimated to
$R_{0}=5.25$. In addition, we estimate that the ratio unidentified infectious
individuals and confirmed cases at the beginning of the epidemic is about $10$,
in agreement with previous studies. We also estimated the epidemic doubling
time to be $2.72$ days.


  The spike glycoprotein of the SARS-CoV-2 virus, which causes COVID-19, has
attracted attention for its vaccine potential and binding capacity to host cell
surface receptors. Much of this research focus has centered on the ectodomain
of the spike protein. The ectodomain is anchored to a transmembrane region,
followed by a cytoplasmic tail. Here we report a distant sequence similarity
between the cysteine-rich cytoplasmic tail of the coronavirus spike protein and
the hepcidin protein that is found in humans and other vertebrates. Hepcidin is
thought to be the key regulator of iron metabolism in humans. An implication of
this preliminary observation is to suggest a potential route of investigation
in the coronavirus research field making use of an already-established
literature on the interplay of local and systemic iron regulation,
cytokine-mediated inflammatory processes, respiratory infections and the
hepcidin protein. The question of possible homology and an evolutionary
connection between the viral spike protein and hepcidin is not assessed in this
report, but some scenarios for its study are discussed.


  The novel coronavirus and its deadly outbreak have posed grand challenges to
human society: as of March 26, 2020, there have been 85,377 confirmed cases and
1,293 reported deaths in the United States; and the World Health Organization
(WHO) characterized coronavirus disease (COVID-19) - which has infected more
than 531,000 people with more than 24,000 deaths in at least 171 countries - a
global pandemic. A growing number of areas reporting local sub-national
community transmission would represent a significant turn for the worse in the
battle against the novel coronavirus, which points to an urgent need for
expanded surveillance so we can better understand the spread of COVID-19 and
thus better respond with actionable strategies for community mitigation. By
advancing capabilities of artificial intelligence (AI) and leveraging the
large-scale and real-time data generated from heterogeneous sources (e.g.,
disease related data from official public health organizations, demographic
data, mobility data, and user geneated data from social media), in this work,
we propose and develop an AI-driven system (named $\alpha$-Satellite}, as an
initial offering, to provide hierarchical community-level risk assessment to
assist with the development of strategies for combating the fast evolving
COVID-19 pandemic. More specifically, given a specific location (either user
input or automatic positioning), the developed system will automatically
provide risk indexes associated with it in a hierarchical manner (e.g., state,
county, city, specific location) to enable individuals to select appropriate
actions for protection while minimizing disruptions to daily life to the extent
possible. The developed system and the generated benchmark datasets have been
made publicly accessible through our website. The system description and
disclaimer are also available in our website.


  Pandemic distribution of COVID-19 in the world has motivated us to discuss
combined effects of network clustering and adaptivity on epidemic spreading. We
address the question concerning the choice of optimal mechanism for most
effective prohibiting disease propagation in a connected network: adaptive
clustering, which mimics self-isolation (SI) in local communities, or sharp
instant clustering, which looks like frontiers closing (FC) between cities and
countries. SI-networks are "adaptively grown" under condition of maximization
of small cliques in the entire network, while FC-networks are "instantly
created". Running the standard SIR model on clustered SI- and FC-networks, we
demonstrate that the adaptive network clustering prohibits the epidemic
spreading better than the instant clustering in the network with similar
parameters. We found that SI model has scale-free property for degree
distribution $P(k)\sim k^{\eta}$ with small critical exponent $-2<\eta<-1$ and
argue that scale-free behavior emerges due to the randomness in the initial
degree distributions and is absent for random regular graphs.


  This paper serves as a framework for designing advanced models for drug
action on metabolism. Drug treatment may affect metabolism by either enhancing
or inhibiting metabolic reactions comprising a metabolic network. We introduce
the concept of \textit{metabolic graphs}, a generalization of hypergraphs
having specialized features common to metabolic networks.
Linear-in-flux-expression (briefly LIFE) is a methodology for analyzing
metabolic networks and simulating virtual patients. We extend LIFE dynamics to
be compatible with metabolic graphs, including the more complex interactions of
enhancer and inhibitor molecules that affect biochemical reactions. We discuss
results considering network structure required for existence and uniqueness of
equilibria on metabolic graphs and show simulations of drug action on
\textit{Mycobacterium tuberculosis} (briefly MTB)


  The recent coronavirus pandemic follows in its early stages an almost
exponential growth, with the number of cases quite well fit in time by
$N(t)\propto e^{\alpha t}$, in many countries. We analyze the rate $\alpha$ for
each country, starting from a threshold of 30 total cases and using the next 12
days, capturing thus the early growth homogeneously. We look for a link between
$\alpha$ and the average temperature $T$ of each country, in the month of the
epidemic growth. We analyze a {\it base} set of 42 countries, which developed
the epidemic earlier, an {\it intermediate} set of 88 countries and an {\it
extended} set of 125 countries, which developed the epidemic more recently.
Applying a linear fit $\alpha(T)$, we find increasing evidence for a decreasing
$\alpha$ as a function of $T$, at $99.66\%$C.L., $99.86\%$C.L. and $99.99995
\%$ C.L. ($p$-value $5 \cdot 10^{-7}$, or 5$\sigma$ detection) in the {\it
base}, {\it intermediate} and {\it extended} dataset, respectively. The
doubling time is expected to increase by $40\%\sim 50\%$, going from $5^\circ$
C to $25^\circ$ C. In the {\it base} set, going beyond a linear model, a peak
at $(7.7\pm 3.6)^\circ C$ seems to be present, but its evidence disappears for
the larger datasets. We also analyzed a possible bias: poor countries, often
located in warm regions, might have less intense testing. By excluding
countries below a given GDP per capita, we find that our conclusions are only
slightly affected and only for the {\it extended} dataset. The significance
remains high, with a $p$-value of $10^{-3}-10^{-4}$ or less. Our findings give
hope that, for northern hemisphere countries, the growth rate should
significantly decrease as a result of both warmer weather and lockdown
policies. In general the propagation should be hopefully stopped by strong
lockdown, testing and tracking policies, before the arrival of the cold season.


  Thermodynamic profiles, based on bioinformatic hydropathicity scales and
membrane-fitted sliding windows, exhibit accurately level sets of hydrophobic
peaks or valleys in the constant fragment Fc of monoclonal antibodies. These
sets provide a roadmap for selected short regions where a few directed amino
acid mutations could produce either accelerated antibody kinetics or increased
stability. These profiles may reveal dramatic differences in antibody activity