  The lack of recommended drugs or vaccines to deal with the COVID-19 is the
main concern of this pandemic. The approved drugs for similar health problems,
drugs under clinical trials, and molecules from medicinal plants extracts are
investigated randomly to deal with the COVID-19 infection. Molecular docking,
one of the best approach to search therapeutically potent drugs/molecules in
real time with possible hope to apply on COVID-19. In this communication,
molecular docking studies of 18 ligands were carried out with the three
therapeutic target proteins of SARS-CoV-2, i.e., RNA-dependent RNA polymerase
(RdRp), angiotensin-converting enzyme 2 (ACE2) and spike glycoprotein (SGp).
The obtained results revealed that the phytochemicals showed better dock score
in compared to the drugs paracetmol and hydroxychloroquine. Combining the dock
score and medicinal properties, we believe the terpenoids based phytochemicals
limonin and scopadulcic acid B can be further explored for potential use
against COVID-19.


  Computational models for the simulation of the severe acute respiratory
syndrome coronavirus 2 (SARS-CoV-2) epidemic evolution would be extremely
useful to support authorities in designing healthcare policies and lockdown
measures to contain its impact on public health and economy. In Italy, the
devised forecasts have been mostly based on a pure data-driven approach, by
fitting and extrapolating open data on the epidemic evolution collected by the
Italian Civil Protection Center. In this respect, SIR epidemiological models,
which start from the description of the nonlinear interactions between
population compartments, would be a much more desirable approach to understand
and predict the collective emergent response. The present contribution
addresses the fundamental question whether a SIR epidemiological model,
suitably enriched with asymptomatic and dead individual compartments, could be
able to provide reliable predictions on the epidemic evolution. To this aim, a
machine learning approach based on particle swarm optimization (PSO) is
proposed to automatically identify the model parameters based on a training set
of data of progressive increasing size, considering Lombardy in Italy as a case
study. The analysis of the scatter in the forecasts shows that model
predictions are quite sensitive to the size of the dataset used for training,
and that further data are still required to achieve convergent -- and therefore
reliable -- predictions.


  Unraveling patterns of animals' movements is important for understanding the
fundamental basics of biogeography, tracking range shifts resulting from
climate change, predicting and preventing biological invansions. Many
researchers have modeled animals' dispersal studying their behavior under the
assumptions of some movement strategies pre-determined or affected by some
external factor(s) but none of them have compared the efficiency of different
dispersal strategies in providing population survival and fitness. We
hypothesize that 1) successful expansion could result from some evolutionary
stable strategy (ESS) and 2) such strategy could be based particularly on
deferred gain, when animals invest in travel to reach some high-quality habitat
("habitat of dream"). Using simulation model we compare the ecological success
of three strategies: i) "Smart" - choosing the locally optimal cell; ii)
"Random" - random movement between cells without taking into account the
quality of the environment; iii) "Dreamer" - movements that aims to find "a
habitat of dream" with quality much higher than that of the initial and
neighboring cells. The population fitness was measured as survival rate,
dispersal distance, accumulated energy and quality of settled habitat. The most
general conclusion is that while survival and wealth of the population is
affected presumably by overall habitat quality, the dispersal depends mainly on
the behavioral strategy. The "Dreamer" strategy or the strategy of deferred
gain belongs to the Pareto frontier in the Fitness$\times$Dispersal space but
only in optimal and suboptimal habitat and in the relatively mild climate.


  We consider the control of the COVID-19 pandemic through a standard SIR
compartmental model. This control is induced by the aggregation of individuals'
decisions to limit their social interactions: when the epidemic is ongoing, an
individual can diminish his/her contact rate in order to avoid getting
infected, but this effort comes at a social cost. If each individual lowers
his/her contact rate, the epidemic vanishes faster, but the effort cost may be
high. A Mean Field Nash equilibrium at the population level is formed,
resulting in a lower effective transmission rate of the virus. We prove
theoretically that equilibrium exists and compute it numerically. However, this
equilibrium selects a sub-optimal solution in comparison to the societal
optimum (a centralized decision respected fully by all individuals), meaning
that the cost of anarchy is strictly positive. We provide numerical examples
and a sensitivity analysis, as well as an extension to a SEIR compartmental
model to account for the relatively long latent phase of the COVID-19 disease.
In all the scenarii considered, the divergence between the individual and
societal strategies happens both before the peak of the epidemic, due to
individuals' fears, and after, when a significant propagation is still
underway.


  We propose an SEIARD mathematical model to investigate the current outbreak
of coronavirus disease (COVID-19) in Mexico. We conduct a detailed analysis of
this model and demonstrate its application using publicly reported data. We
calculate the basic reproduction number ($R_0$) via the next-generation matrix
method, and we estimate the per day infection, death and recovery rates. We
calibrate the parameters of the SEIARD model to the reported data by minimizing
the sum of squared errors and attempt to forecast the evolution of the outbreak
until June 2020. Our results estimate that the peak of the epidemic in Mexico
will be around May 2, 2020. Our model incorporates the importance of
considering the aysmptomatic infected individuals, because they represent the
majority of the infected population (with symptoms or not) and they could play
a huge role in spreading the virus without any knowledge.


  We propose a methodology for the identification of transcription factors
involved in the deregulation of genes in tumoral cells. This strategy is based
on the inference of a reference gene regulatory network that connects
transcription factors to their downstream targets using gene expression data.
The behavior of genes in tumor samples is then carefully compared to this
network of reference to detect deregulated target genes. A linear model is
finally used to measure the ability of each transcription factor to explain
those deregulations. We assess the performance of our method by numerical
experiments on a breast cancer data set. We show that the information about
deregulation is complementary to the expression data as the combination of the
two improves the supervised classification performance of samples into cancer
subtypes.


  Self-activation coupled to a transport mechanism results in traveling waves
that describe polymerization reactions, forest fires, tumor growth, and even
the spread of epidemics. Diffusion is a simple and commonly used model of
particle transport. Many physical and biological systems are, however, better
described by persistent random walks that switch between multiple states of
ballistic motion. So far, traveling fronts in persistent random walk models
have only been analyzed in special, simplified cases. Here, we formulate the
general model of reaction-transport processes in such systems and show how to
compute the expansion velocity for arbitrary number of states. For the
two-state model, we obtain a closed-form expression for the velocity and report
how it is affected by different transport and replication parameters. We also
show that nonzero death rates result in a discontinuous transition from
quiescence to propagation. We compare our results to a recent observation of a
discontinuous onset of propagation in microtubule asters and comment on the
universal nature of the underlying mechanism.


  We discovered that the time evolution of the inverse fractional daily growth
of new infections, N/dN, in the current outbreak of COVID-19 is accurately
described by a universal function, namely the two-parameter Gumbel cumulative
function, in all countries that we have investigated. While the two Gumbel
parameters, as determined bit fits to the data, vary from country to country
(and even within different regions of the same country), reflecting the
diversity and efficacy of the adopted containment measures, the functional form
of the evolution of N/dN appears to be universal. The result of the fit in a
given region or country appears to be stable against variations of the selected
time interval. This makes it possible to robustly estimate the two parameters
from the data data even over relatively small time periods. In turn, this
allows one to predict with large advance and well-controlled confidence levels,
the time of the peak in the daily new infections, its magnitude and duration
(hence the total infections), as well as the time when the daily new infections
decrease to a pre-set value (e.g. less than about 2 new infections per day per
million people), which can be very useful for planning the reopening of
economic and social activities. We use this formalism to predict and compare
these key features of the evolution of the COVID-19 disease in a number of
countries and provide a quantitative assessment of the degree of success in in
their efforts to countain the outbreak.


  The SARS-CoV-2 virus has resulted in a devastating pandemic of COVID-19.
Exploring compounds that could offer a breakthrough in treatment is the need of
the hour. Re-positioning cheap, freely available and safe drugs is a priority.
The paper proposes evidence for the potential use of diethylcarbamazine (DEC)
in the treatment of COVID-19. DEC has inhibitory effects on arachidonic acid
metabolism to prostaglandins, little known anti-viral effects on animal
retroviruses and demonstrated anti-inflammatory actions in animal models of
lung inflammation indicating the need to explore this hypothesis further. We
believe this is the first time DEC is being proposed to treat COVID-19.


  The standard methods to calculate the Total Fertility Rate require the
reliable age-specific fertility rate including birth data and the related
age-specific women's population data. Historically, the number of births was
often not counted according to the age of the mother, so it is difficult to
estimate the historical total fertility rate with the standard methods. Many
empirical methods have been proposed, but their application is limited to
specific period and place. This paper deduces a new method for calculating the
total fertility rate from the number of birth and the population of women at
childbearing age, can be applied to most of cases. The relative error is
usually less than 5%. It is easier to calculate TFR, and may be applied to
obtain more TFRs for the history.


  To control the current outbreak of the Coronavirus Disease 2019, constant
monitoring of the epidemic is required since, as of today, no vaccines or
antiviral drugs against it are known. We provide daily updated estimates of the
reproduction number over time at
https://stochastik-tu-ilmenau.github.io/COVID-19/. In this document, we
describe the estimator we are using which was developed in (Fraser 2007),
derive its asymptotic properties, and we give details on its implementation.
Furthermore, we validate the estimator on simulated data, demonstrate that
estimates on real data lead to plausible results, and perform a sensitivity
analysis. Finally, we discuss why the estimates obtained need to be interpreted
with care.


  We have proposed an SEIR compartmental mathematical model. The prime
objective of this study is to analyze and forecast the pandemic in Italy for
the upcoming months. The basic reproduction number has been calculated. Based
on the current situation in Italy, in this paper, we will estimate the possible
time for the end of the pandemic in the country. The impact of lockdown and
rapid isolation on the spread of the pandemic are also discussed. Further, we
have studied four of the most pandemic affected regions in Italy. Using the
proposed model, a prediction has been made about the duration of pandemic in
these regions. The variation in the basic reproduction number corresponding to
the sensitive parameters of the model is also examined.


  We create a network model to study the spread of an epidemic through
physically proximate and accidental daily human contacts in a city, and
simulate outcomes for two kinds of agents - poor and non-poor. Under
non-intervention, peak caseload is maximised, but no differences are observed
in infection rates across poor and non-poor. Introducing interventions to
control spread, peak caseloads are reduced, but both cumulative infection rates
and current infection rates are systematically higher for the poor than for
non-poor, across all scenarios. Larger populations, higher fractions of poor,
and longer durations of intervention are found to progressively worsen outcomes
for the poor; and these are of particular concern for economically vulnerable
populations in cities of the developing world. Addressing these challenges
requires a deeper, more rigorous understanding of the relationships between
structural poverty and epidemy, as well as effective utilization of extant
community level infrastructure for primary care in developing cities. Finally,
improving iniquitous outcomes for the poor creates better outcomes for the
whole population, including the non-poor.


  Increasing the infection risk early in an epidemic is individually and
socially optimal under some parameter values. The reason is that the early
patients recover or die before the peak of the epidemic, which flattens the
peak. This improves welfare if the peak exceeds the capacity of the healthcare
system and the social loss rises rapidly enough in the number infected. The
individual incentive to get infected early comes from the greater likelihood of
receiving treatment than at the peak when the disease has overwhelmed
healthcare capacity. Calibration to the Covid-19 pandemic data suggests that
catching the infection at the start was individually optimal and for some loss
functions would have reduced the aggregate loss.


  We construct a zig-zag process targeting a posterior distribution defined on
a hybrid state space consisting of both discrete and continuous variables. The
construction does not require any assumptions on the structure among discrete
variables. We demonstrate our method on two examples in genetics based on the
Kingman coalescent, showing that the zig-zag process can lead to efficiency
gains of up to several orders of magnitude over classical Metropolis-Hastings
algorithms, and that it is well suited to parallel computation. Our
construction resembles existing techniques for Hamiltonian Monte Carlo on a
hybrid state space, which suffers from implementationally and analytically
complex boundary crossings when applied to the coalescent. We demonstrate that
the continuous-time zig-zag process avoids these complications.


  No, they can't. Epidemic spread is characterized by exponentially growing
dynamics, which are intrinsically unpredictable. The time at which the growth
in the number of infected individuals halts and starts decreasing cannot be
calculated with certainty before the turning point is actually attained;
neither can the end of the epidemic after the turning point. An SIR model with
confinement (SCIR) illustrates how lockdown measures inhibit infection spread
only above a threshold that we calculate. The existence of that threshold has
major effects in predictability: A Bayesian fit to the COVID-19 pandemic in
Spain shows that a slow-down in the number of newly infected individuals during
the expansion phase allows to infer neither the precise position of the maximum
nor whether the measures taken will bring the propagation to the inhibition
regime. There is a short horizon for reliable prediction, followed by a
dispersion of the possible trajectories that grows extremely fast. The
impossibility to predict in the mid-term is not due to wrong or incomplete
data, since it persists in error-free, synthetically produced data sets, and
does not necessarily improve by using larger data sets. Our study warns against
precise forecasts of the evolution of epidemics based on mean-field, effective
or phenomenological models, and supports that only probabilities of different
outcomes can be confidently given.


  We present a System Dynamics (SD) model of the Covid-19 pandemic spread in
India. The detailed age-structured compartment-based model endogenously
captures various disease transmission pathways, expanding significantly from
the standard SEIR model. The model is customized for India by using the
appropriate population pyramid, contact rate matrices, external arrivals (as
per actual data), and a few other calibrated fractions based on the reported
cases of Covid-19 in India. Also, we have explicitly modeled, using independent
time-variant levers, the effects of testing, contact tracing, isolating
Covid-positive patients, quarantining, use of mask/better hygiene practices,
social distancing through contact rate reductions at distinct zones of home(H),
work(W), school(S) and other(O) locations. Simulation results show that, even
after an extended lock-down, some non-trivial number of infections (even
asymptomatic) will be left and the pandemic will resurface. Only tools that
work against the pandemic is high rate of testing of those who show Covid-19
like symptoms, isolating them if they are positive and contact tracing all
contacts of positive patients and quarantining them, in combination with use of
face masks and personal hygiene. A wide range of combination of effectiveness
of contact tracing, isolation, quarantining and personal hygiene measures help
minimize the pandemic impact and some imperfections in implementation of one
measure can be compensated by better implementation of other measures.


  We develop here a data-driven approach for disease recognition based on given
symptoms, to be efficient tool for anomaly detection. In a clinical setting and
when presented with a patient with a combination of traits, a doctor may wonder
if a certain combination of symptoms may be especially predictive, such as the
question, "Are fevers more informative in women than men?" The answer to this
question is, yes. We develop here a methodology to enumerate such questions, to
learn what are the stronger warning signs when attempting to diagnose a
disease, called Conditional Predictive Informativity, (CPI), whose ranking we
call CPIR. This simple to use process allows us to identify particularly
informative combinations of symptoms and traits that may help medical field
analysis in general, and possibly to become a new data-driven advised approach
for individual medical diagnosis, as well as for broader public policy
discussion. In particular we have been motivated to develop this tool in the
current enviroment of the pressing world crisis due to the COVID 19 pandemic.
We apply the methods here to data collected from national, provincial, and
municipal health reports, as well as additional information from online, and
then curated to an online publically available Github repository.


  2019-nCoV is a novel coronavirus was isolated and identified in 2019 in
Wuhan, China. On 17th February and according to world health organization, a
number of 71 429 confirmed cases worldwide, among them 2162 new cases recorded
in the last 24 hours. There is no drug or vaccine for human and animal
coronavirus. The inhibition of 3CL hydrolase enzyme provides a promising
therapeutic principle for developing treatments against CoViD-19. The 3CLpro
(Mpro) known for involving in counteracting the host innate immune response.
This work presents the inhibitory effect of some natural compounds against 3CL
hydrolase enzyme, and explain the main interactions in inhibitor-enzyme
complex. Molecular docking study carried out using Autodock Vina. By screening
several molecules, we identified three candidate agents that inhibit the main
protease of coronavirus. Hispidin, lepidine E, and folic acid bound tightly in
the enzyme, strong hydrogen bonds have been formed (1.69-1.80&[Aring]) with the
active site residues. This study provides a possible therapeutic strategy for
CoViD-19.


  We consider the discrete-time migration-recombination equation, a
deterministic, nonlinear dynamical system that describes the evolution of the
genetic type distribution of a population evolving under migration and
recombination in a law of large numbers setting. We relate this dynamics
(forward in time) to a Markov chain, namely a labelled partitioning process,
backward in time. This way, we obtain a stochastic representation of the
solution of the migration-recombination equation. As a consequence, one obtains
an explicit solution of the nonlinear dynamics, simply in terms of powers of
the transition matrix of the Markov chain. Finally, we investigate the limiting
and quasi-limiting behaviour of the Markov chain, which gives immediate access
to the asymptotic behaviour of the dynamical system. We finally sketch the
analogous situation in continuous time.


  In the absence of other tools, monitoring the effects of protective measures,
including social distancing and forecasting the outcome of outbreaks is of
immense interest. Real-time data is noisy and very often hampered by systematic
errors in reporting. Detailed epidemic models may contain a large number of
empirical parameters, which cannot be determined with sufficient accuracy. In
this paper, we show that the cumulative number of deaths can be regarded as a
master variable, and the parameters of the epidemic such as the basic
reproduction number, the size of the susceptible population, and the infection
rate can be determined. In the SIR model, we derive an explicit single variable
differential equation for the evolution of the cumulative number of fatalities.
We show that the epidemic in Spain, Italy, and Hubei Province, China follows
this master equation closely. We discuss the relationship with the logistic
growth model, and we show that it is a good approximation when the basic
reproduction number is less than $2.3$. This condition is valid for the
outbreak in Hubei, but not for the outbreaks in Spain, Italy, and New York. The
difference is in the shorter infectious period in China, probably due to the
separation policy of the infected. For more complex models, with more internal
variables, such as the SEIR model, the equations derived from the SIR model
remain valid approximately, due to the separation of timescales.


  In the midst of pandemic for respiratory illness, the call for
non-pharmaceutical interventions become the highest priority for infectious
disease and public health experts, while the race towards vaccine or medical
intervention are ongoing. Individuals may modify their behavior and take
preventative steps to reduce infection risk in the bid to adhere to the call by
government officials and experts. As a result, the existence of relationship
between the preliminary and the final transmission rates become feeble. This
study evaluates the behavioral changes (mitigation and suppression measures)
proposed by public health experts for COVID-19 which had altered human behavior
and their day to day lives. The dynamics underlying the mitigation and
suppression measures reduces the contacts among citizens and significantly
interfere with their physical and social behavior. The results show all the
measures have a significant impact on the decline of transmission rate.
However, the mitigation measures might prolong the elimination of the
transmission which might lead to a severe economic meltdown, yet, a combination
of the measures show a possibility of rooting out transmission within 30 days
if adhered to in an extreme manner. The result shows a peak period of infection
for Ghana ranges from 64th day to 74th day of infection time period.


  Sweden has adopted far less restrictive social distancing policies than most
countries following the COVID-19 pandemic. This paper uses data on all mobile
phone users, from one major Swedish mobile phone network, to examine the impact
of the Coronavirus outbreak under the Swedish mild recommendations and
restrictions regime on individual mobility and if changes in geographical
mobility vary over different socio-economic strata. Having access to data for
January-March in both 2019 and 2020 enables the estimation of causal effects of
the COVID-19 outbreak by adopting a Difference-in-Differences research design.
The paper reaches four main conclusions: (i) The daytime population in
residential areas increased significantly (64 percent average increase); (ii)
The daytime presence in industrial and commercial areas decreased significantly
(33 percent average decrease); (iii) The distance individuals move from their
homes during a day was substantially reduced (38 percent decrease in the
maximum distance moved and 36 percent increase in share of individuals who move
less than one kilometer from home); (iv) Similar reductions in mobility were
found for residents in areas with different socioeconomic and demographic
characteristics. These results show that mild government policies can compel
people to adopt social distancing behavior.


  Social, supervised, learning from others might amplify individual, possibly
unsupervised, learning by individuals, and might underlie the development and
evolution of culture. We studied a minimal model of the interaction of
individual unsupervised and social supervised learning by interacting agents.
Agents attempted to learn to track a hidden fluctuation "source", which,
linearly mixed with other masking fluctuations, generated observable input
vectors. Learning was driven either solely by direct observation of inputs
(unsupervised, Hebbian) or, in addition, by observation of another agent's
output (supervised, Delta rule). To enhance biological realism, the learning
rules were made slightly connection-inspecific, so that incorrect learning
sometimes occurs. We found that social interaction can foster both correct and
incorrect learning. Useful social learning therefore presumably involves
additional factors some of which we outline.


  CRISPR systems experience off-target effects that interfere with the ability
to accurately perform genetic edits. While empirical models predict off-target
effects in specific platforms, there is a gap for a wide-ranging mechanistic
model of CRISPR systems based on evidence from the essential processes of
target recognition. In this work, we present a first principles model supported
by many experiments performed in vivo and in vitro. First, we observe and model
the conformational changes and discrete stepped DNA unwinding events of
SpCas9-CRISPR target recognition process in a cell-free
transcription-translation system using truncated guide RNAs, confirming
structural and FRET microscopy experiments. Second, we implement an energy
landscape to describe single mismatch effects observed in vivo for spCas9 and
Cas12a CRISPR systems. From the mismatch analysis results we uncover that
CRISPR class 2 systems maintain their kinetic proofreading mechanism by
utilizing intermittent energetic barriers and we show how to leverage the
landscape to improve target specificity.


  We propose a simple SIR model in order to investigate the impact of various
confinement strategies on a most virulent epidemic. Our approach is motivated
by the current COVID-19 pandemic. The main hypothesis is the existence of two
populations of susceptible persons, one which obeys confinement and for which
the infection rate does not exceed 1, and a population which, being non
confined for various imperatives, can be substantially more infective. The
model, initially formulated as a differential system, is discretised following
a specific procedure, the discrete system serving as an integrator for the
differential one. Our model is calibrated so as to correspond to what is
observed in the COVID-19 epidemic.
  Several conclusions can be reached, despite the very simple structure of our
model. First, it is not possible to pinpoint the genesis of the epidemic by
just analysing data from when the epidemic is in full swing. It may well turn
out that the epidemic has reached a sizeable part of the world months before it
became noticeable. Concerning the confinement scenarios, a universal feature of
all our simulations is that relaxing the lockdown constraints leads to a
rekindling of the epidemic. Thus we sought the conditions for the second
epidemic peak to be lower than the first one. This is possible in all the
scenarios considered (abrupt, progressive or stepwise exit) but typically a
progressive exit can start earlier than an abrupt one. However, by the time the
progressive exit is complete, the overall confinement times are not too
different. From our results, the most promising strategy is that of a stepwise
exit. And in fact its implementation could be quite feasible, with the major
part of the population (minus the fragile groups) exiting simultaneously but
obeying rigorous distancing constraints.


  Aim of this manuscript is to show a simple method to infer the time-course of
new COVID-19 infections (the most important information in order to establish
the effect of containment strategies) from available aggregated data, such as
number of deaths and hospitalizations. The method, that was used for HIV-AIDS
and was named `back-calculation', relies on good estimates of the distribution
of the delays between infection and the observed events; assuming that the
epidemic follows a simple SIR model with a known generation interval, we can
then estimate the parameters that define the time-varying contact rate through
maximum likelihood. We show the application of the method to data from Italy
and several of its region; it is found that $R_0$ had decreased consistently
below 1 around March 20, and in the beginning of April it was between 0.5 and
0.8 in the whole Italy and in most regions.


  When effective medical treatment and vaccination are not available,
non-pharmaceutical interventions such as social distancing, home quarantine and
far-reaching shutdown of public life are the only available strategies to
prevent the spread of epidemics. Based on an extended SEIR
(susceptible-exposed-infectious-recovered) model and continuous-time optimal
control theory, we compute the optimal non-pharmaceutical intervention strategy
for the case that a vaccine is never found and complete containment
(eradication of the epidemic) is impossible. In this case, the optimal control
must meet competing requirements: First, the minimization of disease-related
deaths, and, second, the establishment of a sufficient degree of natural
immunity at the end of the measures, in order to exclude a second wave.
Moreover, the socio-economic costs of the intervention shall be kept at a
minimum. The numerically computed optimal control strategy is a
single-intervention scenario that goes beyond heuristically motivated
interventions and simple "flattening of the curve." Careful analysis of the
computed control strategy reveals, however, that the obtained solution is in
fact a tightrope walk close to the stability boundary of the system, where
socio-economic costs and the risk of a new outbreak must be constantly balanced
against one another. The model system is calibrated to reproduce the initial
exponential growth phase of the COVID-19 pandemic in Germany.


  The coronavirus disease 2019 (COVID-19) has changed the world since the World
Health Organization declared its outbreak on 30th January 2020, recognizing the
outbreak as a pandemic on 11th March 2020. As often said by politicians and
scientific advisors, the objective is "to flatten the curve", or "push the peak
down", or similar wording, of the virus spreading. Central to the official
advice are mathematical models and data, which provide estimates on the
evolution of the number of infected, recovered and deaths. The accuracy of the
models is improved day by day by inferring the contact, recovery, and death
rates from data (confirmed cases). A data-driven model trained with {\it both}
data {\it and} first principles is proposed. The model can quickly be
re-trained any time that new data becomes available. The method can be applied
to more detailed epidemic models with virtually no conceptual modification.


  We combine a pedestrian dynamics model with a contact tracing method to
simulate the initial spreading of a highly infectious airborne disease in a
confined environment. We focus on a medium size population (up to 1000 people)
with a small number of infectious people (1 or 2) and the rest of the people is
divided between immune and susceptible. We adopt a space-continuous model that
represents pedestrian dynamics by the forces acting on them, i.e. a microscopic
force-based model. Once discretized, the model results in a high-dimensional
system of second order ordinary differential equations. Before adding the
contact tracing to the pedestrian dynamics model, we calibrate the model
parameters, validate the model against empirical data, and show that pedestrian
self-organization into lanes can be captured. We consider an explicit approach
for contact tracing by introducing a sickness domain around a sick person. A
healthy but susceptible person who remains in the sickness domain for a certain
amount of time may get infected (with a prescribed probability) and become a
so-called secondary contact. As a concrete setting to simulate the onset of
disease spreading, we consider terminals in two US airports: Hobby Airport in
Houston and the Atlanta International Airport. We consider different scenarios
and we quantify the increase in average number of secondary contacts increases
as a given terminal becomes more densely populated, the percentage of immune
people decreases, the number of primary contacts increases, and areas of high
density (such as the boarding buses) are present.


  We propose a strategy based on the site-bond percolation to minimize the
propagation of \textit{Phytophthora} zoospores on plantations, consisting in
introducing physical barriers between neighboring plants. Two clustering
processes are distinguished: i) one of cells with the presence of the pathogen,
detected on soil analysis; and ii) that of diseased plants, revealed from a
visual inspection of the plantation. The former is well described by the
standard site-bond percolation. In the latter, the percolation threshold is
fitted by a Tsallis distribution when no barriers are introduced. We provide,
for both cases, the formulae for the minimal barrier density to prevent the
emergence of the spanning cluster. Though this work is focused on a specific
pathogen, the model presented here can also be applied to prevent the spreading
of other pathogens that disseminate, by other means, from one plant to the
neighboring ones. Finally, the application of this strategy to three types of
commercialy important Mexican chili plants is also shown.


  We present a novel method for automated identification of putative cell types
from single-cell RNA-seq (scRNA-seq) data. By iteratively applying a machine
learning approach to an initial clustering of gene expression profiles of a
given set of cells, we simultaneously identify distinct cell groups and a
weighted list of feature genes for each group. The feature genes, which are
differentially expressed in the particular cell group, jointly discriminate the
given cell group from other cells. Each such group of cells corresponds to a
putative cell type or state, characterised by the feature genes as markers. To
benchmark this approach, we use expert-annotated scRNA-seq datasets from a
range of experiments, as well as comparing to existing cell annotation methods,
which are all based on a pre-existing reference. We show that our method
automatically identifies the 'ground truth' cell assignments with high
accuracy. Moreover, our method, Single Cell Clustering Assessment Framework
(SCCAF) predicts new putative biologically meaningful cell-states in published
data on haematopoiesis and the human cortex. SCCAF is available as an
open-source software package on GitHub (https://github.com/SCCAF/sccaf) and as
a Python package index and has also been implemented as a Galaxy tool in the
Human Cell Atlas.


  The recent spreading of coronavirus made many countries impose restrictions
in order to control its dangerous effect on the citizens. We developed a
theoretical dynamical model based on compartmental SIR system with additional
adjustment taken from Flow network and Markov chain frameworks to illustrate
developments and trends based on publicly available data. Based on this Model,
code in R was written and fed by stamped publicly available data from
authorized governmental websites in Spain and in the Nederlands, to compare
trends. Our results show that the 'peak' of infection is already behind us in
both countries, but also demonstrate that there is a danger of rebound of a
spread. It is obvious that measures imposed are giving the results, but we
should be precarious of near future practices and developments since the
majority of population will still be without immunity.


  In March of 2020, many U.S. state governments encouraged or mandated
restrictions on social interactions to slow the spread of COVID-19, the disease
caused by the novel coronavirus SARS-CoV-2 that has spread to nearly 180
countries. Estimating the effectiveness of these social-distancing strategies
is challenging because surveillance of COVID-19 has been limited, with tests
generally being prioritized for high-risk or hospitalized cases according to
temporally and regionally varying criteria. Here we show that reductions in
mobility across U.S. counties with at least 100 confirmed cases of COVID-19 led
to reductions in fever incidences, as captured by smart thermometers, after a
mean lag of 6.5 days ($90\%$ within 3--10 days) that is consistent with the
incubation period of COVID-19. Furthermore, counties with larger decreases in
mobility subsequently achieved greater reductions in fevers ($p<0.01$), with
the notable exception of New York City and its immediate vicinity. These
results indicate that social distancing has reduced the transmission of
influenza like illnesses, including COVID 19, and support social distancing as
an effective strategy for slowing the spread of COVID-19.


  This work presents a simple and realistic approach to handle the available
data of COVID-19 patients in India and to forecast the scenario. The model
proposed is based on the available facts like the onset of lockdown (as
announced by the Government on 25th day, {\tau}0 and the recovery pattern
dictated by a mean life recovery time of {\tau}1 ( normally said to be around
14 days). The data of infected COVID-19 patients from March 2, to April 16,
2020 has been used to fit the evolution of infected, recovery and death counts.
A slow rising exponential growth, with R0 close to 1/6, is found to represent
the infected counts indicating almost a linear rise. The rest of growth,
saturation and decay of data is comprehensibly modelled by incorporating
lockdown time controlled R0, having a normal error function like behaviour
decaying to zero in some time frame of {\tau}2 . The recovery mean life time
{\tau}1 dictates the peak and decay. The results predicted for coming days are
interesting and optimistic. The introduced time constants based on experimental
data for both the recovery rate as well as for determining the time span of
activity of R0 after the lockdown are subject of debate and provide possibility
to introduce trigger factors to alter these to be more suited to the model. The
model can be extended to other communities with their own R0 and recovery time
parameters.


  While COVID-19 is rapidly propagating around the globe, the need for
providing real-time forecasts of the epidemics pushes fits of dynamical and
statistical models to available data beyond their capabilities. Here we focus
on statistical predictions of COVID-19 infections performed by fitting
asymptotic distributions to actual data. By taking as a case-study the epidemic
evolution of total COVID-19 infections in Chinese provinces and Italian
regions, we find that predictions are characterized by large uncertainties at
the early stages of the epidemic growth. Those uncertainties significantly
reduce after the epidemics peak is reached. Differences in the uncertainty of
the forecasts at a regional level can be used to highlight the delay in the
spread of the virus. Our results warn that long term extrapolation of epidemics
counts must be handled with extreme care as they crucially depend not only on
the quality of data, but also on the stage of the epidemics, due to the
intrinsically non-linear nature of the underlying dynamics. These results
suggest that real-time epidemiological projections should include wide
uncertainty ranges and urge for the needs of compiling high-quality datasets of
infections counts, including asymptomatic patients.


  To model the spread of COVID-19 coronavirus in Russian regions and in Moscow,
a discrete logistic equation describing the increase in the number of cases is
used. To check the adequacy of the mathematical model, the simulation results
were compared with the spread of coronavirus in China, in a number of European
and Asian countries, and the United States. The parameters of the logistics
equation for Russia, Moscow and other large regions were determined in the
interval (01.03 - 08.04). A comparative analysis of growth rates of COVID-19
infected population for different countries and regions is presented. Various
scenarios of the spread of COVID-19 coronavirus in Moscow and in the regions of
Russia are considered. For each scenario, curves for the daily new cases and
graphs for the increase in the total number of cases were obtained, and the
dynamics of infection spread by day was studied. Peak times, epidemic periods,
the number of infected people at the peak and their growth were determined.


  Social distancing remains an important strategy to combat the COVID-19
pandemic in the United States. However, the impacts of specific state-level
policies on mobility and subsequent COVID-19 case trajectories have not been
completely quantified. Using anonymized and aggregated mobility data from
opted-in Google users, we found that state-level emergency declarations
resulted in a 9.9% reduction in time spent away from places of residence.
Implementation of one or more social distancing policies resulted in an
additional 24.5% reduction in mobility the following week, and subsequent
shelter-in-place mandates yielded an additional 29.0% reduction. Decreases in
mobility were associated with substantial reductions in case growth 2 to 4
weeks later. For example, a 10% reduction in mobility was associated with a
17.5% reduction in case growth 2 weeks later. Given the continued reliance on
social distancing policies to limit the spread of COVID-19, these results may
be helpful to public health officials trying to balance infection control with
the economic and social consequences of these policies.


  Messenger RNA (mRNA) vaccines are being used for COVID-19, but still suffer
from the critical issue of mRNA instability and degradation, which is a major
obstacle in the storage, distribution, and efficacy of the vaccine. Previous
work showed that optimizing secondary structure stability lengthens mRNA
half-life, which, together with optimal codons, increases protein expression.
Therefore, a principled mRNA design algorithm must optimize both structural
stability and codon usage to improve mRNA efficiency. However, due to
synonymous codons, the mRNA design space is prohibitively large, e.g., there
are $\sim\!10^{632}$ mRNAs for the SARS-CoV-2 Spike protein, which poses
insurmountable challenges to previous methods. Here we provide a surprisingly
simple solution to this hard problem by reducing it to a classical problem in
computational linguistics, where finding the optimal mRNA is akin to finding
the most likely sentence among similar sounding alternatives. Our algorithm,
named LinearDesign, takes only 11 minutes for the Spike protein, and can
jointly optimize stability and codon usage. Experimentally, without chemical
modification, our designs substantially improve mRNA half-life and protein
expression in vitro, and dramatically increase antibody response by up to
23$\times$ in vivo, compared to the codon-optimized benchmark. Our work enables
the exploration of highly stable and efficient designs that are previously
unreachable and is a timely tool not only for vaccines but also for mRNA
medicine encoding all therapeutic proteins (e.g., monoclonal antibodies and
anti-cancer drugs).


  We report the results of our study of approved drugs as potential treatments
for COVID 19, based on the application of various bioinformatics predictive
methods. The drugs studied include chloroquine, ivermectin, remdesivir,
sofosbuvir, boceprevir, and {\alpha}-difluoromethylornithine (DMFO). Our
results indicate that these small molecules selectively bind to stable,
kinetically active residues and residues adjoining them on the surface of
proteins and inside protein pockets and that some prefer hydrophobic over other
active sites. Our approach is not restricted to viruses and can facilitate
rational drug design, as well as improve our understanding of molecular
interactions, in general.


  We present a new Bayesian inference method for compartmental models that
takes into account the intrinsic stochasticity of the process. We show how to
formulate a SIR-type Markov jump process as the solution of a stochastic
differential equation with respect to a Poisson Random Measure (PRM), and how
to simulate the process trajectory deterministically from a parameter value and
a PRM realisation.
  This forms the basis of our Data Augmented MCMC, which consists in augmenting
parameter space with the unobserved PRM value. The resulting simple
Metropolis-Hastings sampler acts as an efficient simulation-based inference
method, that can easily be transferred from model to model. Compared with a
recent Data Augmentation method based on Gibbs sampling of individual infection
histories, PRM-augmented MCMC scales much better with epidemic size and is far
more flexible.
  PRM-augmented MCMC also yields a posteriori estimates of the PRM, that
represent process stochasticity, and which can be used to validate the model.
If the model is good, the posterior distribution should exhibit no pattern and
be close to the PRM prior distribution. We illustrate this by fitting a
non-seasonal model to some simulated seasonal case count data. Applied to the
Zika epidemic of 2013 in French Polynesia, our approach shows that a simple
SEIR model cannot correctly reproduce both the initial sharp increase in the
number of cases as well as the final proportion of seropositive.
  PRM-augmentation thus provides a coherent story for Stochastic Epidemic Model
inference, where explicitly inferring process stochasticity helps with model
validation.


  The primary cell surface receptor for SARS-CoV-2 is the
angiotensin-converting enzyme 2 (ACE2). Recently it has been noticed that the
viral Spike protein has an RGD motif, suggesting that cell surface integrins
may be co-receptors. We examined the sequences of ACE2 and integrins with the
Eukaryotic Linear Motif resource, ELM, and were presented with candidate short
linear motifs (SLiMs) in their short, unstructured, cytosolic tails with
potential roles in endocytosis, membrane dynamics, autophagy, cytoskeleton and
cell signalling. These SLiM candidates are highly conserved in vertebrates.
They suggest potential interactions with the AP2 mu2 subunit as well as I-BAR,
LC3, PDZ, PTB and SH2 domains found in signalling and regulatory proteins
present in epithelial lung cells. Several motifs overlap in the tail sequences,
suggesting that they may act as molecular switches, often involving tyrosine
phosphorylation status. Candidate LIR motifs are present in the tails of ACE2
and integrin beta3, suggesting that these proteins can directly recruit
autophagy components. We also noticed that the extracellular part of ACE2 has a
conserved MIDAS structural motif, which are commonly used by beta integrins for
ligand binding, potentially supporting the proposal that integrins and ACE2
share common ligands. The findings presented here identify several molecular
links and testable hypotheses that might help uncover the mechanisms of
SARS-CoV-2 attachment, entry and replication, and strengthen the possibility
that it might be possible to develop host-directed therapies to dampen the
efficiency of viral entry and hamper disease progression. The strong sequence
conservation means that these putative SLiMs are good candidates: Nevertheless,
SLiMs must always be validated by experimentation before they can be stated to
be functional.


  Currently the COVID-19 epidemic is developing in the City of Culiac\'n
Sinaloa, Mexico, where up to April 20 of this year there have been 35 deaths
associated with this epidemic. The reproduction number $(R_0)$ of an epidemic
represents the average number of people infected by an infected person during
their period of infection. In this work we use the data published by the
Secretary of Health of the State of Sinaloa on the number of new infected cases
confirmed per day and we estimate that the value of $R_0$ is 1.562 with a 95%
confidence interval given by (1.401,1.742). We also estimate the mortality rate
among the confirmed cases, which turned out to be 16.8%.
  - - - - -
  Actualmente la epidemia COVID-19 se est\'a desarrollando en la Ciudad de
Culiac\'an Sinaloa, M\'exico, donde hasta el 20 de abril del presente a\~no han
ocurrido 35 decesos asociados con esta epidemia. El n\'umero de reproducci\'on
$(R_0)$ de una epidemia representa el n\'umero promedio de personas contagiadas
por una persona infectada durante su periodo de infecci\'on. En este trabajo
usamos los datos publicados por la Secretaria de Salud del Estado de Sinaloa
sobre el n\'umero de nuevos casos infectados confirmados por dia y estimamos
que el valor de $R_0$ es de 1.562 con un intervalo del 95% de confianza dado
por (1.401,1.742). Estimamos tambi\'en la tasa de mortalidad entre los casos
confirmados, la cual result\'o ser de 16.8%.


  In this paper, we propose a new model for the dynamics of COVID-19
infections. Our approach consists of seven phenotypes: the susceptible humans,
exposed humans, infectious humans, the recovered humans, the quarantine
population, the recovered-exposed and deceased population. We proved first
through mathematical approach the positivity, boundness and existence of a
solution to considered model. We also studied the existence of the disease free
equilibrium and corresp onding stability. Hence, the actual reproduction number
was calculated . We analyzed the dependence of basic reproductions R_0 on the
confidence of our model. Our work shows, in particular, that the disease will
decrease out if the number of reproduction was less than one. Moreover, the
impact of the quarantine strategies to reduce the spread of this disease is
discussed. The theoretical results are validated by some numerical simulations
of the system of the epidemic's differential equations.


  We analyze an approach to managing the COVID-19 pandemic without shutting
down the economy while staying within the capacity of the healthcare system. We
base our analysis on a detailed heterogeneous epidemiological model, which
takes into account different population groups and phases of the disease,
including incubation, infection period, hospitalization, and treatment in the
intensive care unit (ICU). We model the healthcare capacity as the total number
of hospital and ICU beds for the whole country. We calibrate the model
parameters to data reported in several recent research papers. For high- and
low-risk population groups, we calculate the number of total and intensive care
hospitalizations, and deaths as functions of time. The main conclusion is that
countries, which enforce reasonable hygienic measures on time can avoid
lockdowns throughout the pandemic provided that the number of spare ICU beds
per million is above the threshold of about 100. In countries where the total
number of ICU beds is below this threshold, a limited period quarantine to
specific high-risk groups of the population suffices. Furthermore, in the case
of an inadequate capacity of the healthcare system, we incorporate a feedback
loop and demonstrate that quantitative impact of the lack of ICU units on the
death curve. In the case of inadequate ICU beds, full- and partial-quarantine
scenarios outcomes are almost identical, making it unnecessary to shut down the
whole economy. We conclude that only a limited-time quarantine of the high-risk
group might be necessary, while the rest of the economy can remain operational.


  The SARS-CoV-2 infectious outbreak has rapidly spread across the globe and
precipitated varying policies to effectuate physical distancing to ameliorate
its impact. In this study, we propose a new hybrid machine learning model,
SIRNet, for forecasting the spread of the COVID-19 pandemic that couples with
the epidemiological models. We use categorized spatiotemporally explicit
cellphone mobility data as surrogate markers for physical distancing, along
with population weighted density and other local data points. We demonstrate at
varying geographical granularity that the spectrum of physical distancing
options currently being discussed among policy leaders have epidemiologically
significant differences in consequences, ranging from viral extinction to near
complete population prevalence. The current mobility inflection points vary
across geographical regions. Experimental results from SIRNet establish
preliminary bounds on such localized mobility that asymptotically induce
containment. The model can support in studying non-pharmacological
interventions and approaches that minimize societal collateral damage and
control mechanisms for an extended period of time.


  Local-ancestry inference (LAI), also referred to as ancestry deconvolution,
provides high-resolution ancestry estimation along the human genome. In both
research and industry, LAI is emerging as a critical step in DNA sequence
analysis with applications extending from polygenic risk scores (used to
predict traits in embryos and disease risk in adults) to genome-wide
association studies, and from pharmacogenomics to inference of human population
history. While many LAI methods have been developed, advances in computing
hardware (GPUs) combined with machine learning techniques, such as neural
networks, are enabling the development of new methods that are fast, robust and
easily shared and stored. In this paper we develop the first neural network
based LAI method, named LAI-Net, providing competitive accuracy with
state-of-the-art methods and robustness to missing or noisy data, while having
a small number of layers.


  We show that precise knowledge of epidemic transmission parameters is not
required to build an informative model of the spread of disease. We propose a
detailed model of the topology of the contact network under various external
control regimes and demonstrate that this is sufficient to capture the salient
dynamical characteristics and to inform decisions. Contact between individuals
in the community is characterised by a contact graph, the structure of that
contact graph is selected to mimic community control measures. Our model of
city-level transmission of an infectious agent (SEIR model) characterises
spread via a (a) scale-free contact network (no control); (b) a random graph
(elimination of mass gatherings); and (c) small world lattice (partial to full
lockdown -- "social" distancing). This model exhibits good qualitative
agreement between simulation and data from the 2020 pandemic spread of
coronavirus. Estimates of the relevant rate parameters of the SEIR model are
obtained and we demonstrate the robustness of our model predictions under
uncertainty of those estimates. The social context and utility of this work is
identified, contributing to a highly effective pandemic response in Western
Australia.


  In the present paper we provide the closed form of the path-like solutions
for the logistic and $\theta$-logistic stochastic differential equations, along
with the exact expressions of both their probability density functions and
their moments. We simulate in addition a few typical sample trajectories, and
we provide a few examples of numerical computation of the said closed formulas
at different noise intensities: this shows in particular that an increasing
randomness - while making the process more unpredictable - asymptotically tends
to suppress in average the logistic growth. These main results are preceded by
a discussion of the noiseless, deterministic versions of these models: a
prologue which turns out to be instrumental - on the basis of a few simplified
but functional hypotheses - to frame the logistic and $\theta$-logistic
equations in a unified context, within which also the Gompertz model emerges
from an anomalous scaling.


  Extreme public health interventions play a critical role in mitigating the
local and global prevalence and pandemic potential of COVID-19. Here, we use
population size for pathogen transmission to measure the intensity of public
health interventions, which is a key characteristic variable for nowcasting and
forecasting of the epidemic. By formulating a hidden Markov dynamic system and
using nonlinear filtering theory, we have developed a stochastic epidemic
dynamic model under public health interventions. The model parameters and
states are estimated in time from internationally available public data by
combining an unscented filter and an interacting multiple model filter.
Moreover, we consider the computability of the population size and provide its
selection criterion. We estimate the mean of the basic reproductive number of
China and the rest of the globe except China (GEC) to be 2.46 (95% CI:
2.41-2.51) and 3.64 (95% CI: (3.55-3.72), respectively. We infer that the
number of latent infections of GEC is about 7.47*10^5 (95% CI:
7.32*10^5-7.62*10^5) as of April 2, 2020. We predict that the peak of
infections in hospitals of GEC may reach 3.00*10^6 on the present trajectory,
i.e., if the population size for pathogen transmission and epidemic parameters
remains unchanged. If the control intensity is strengthened, e.g., 50%
reduction or 75% reduction of the population size for pathogen transmission,
the peak would decline to 1.84*10^6, 1.27*10^6, respectively.


  We show that the COVID-19 pandemic under social distancing exhibits universal
dynamics. The cumulative numbers of both infections and deaths quickly cross
over from exponential growth at early times to a longer period of power law
growth, before eventually slowing. In agreement with a recent statistical
forecasting model by the IHME, we show that this dynamics is well described by
the erf function. Using this functional form, we perform a data collapse across
countries and US states with very different population characteristics and
social distancing policies, confirming the universal behavior of the COVID-19
outbreak. We show that the predictive power of statistical models is limited
until a few days before curves flatten, forecast deaths and infections assuming
current policies continue and compare our predictions to the IHME models. We
present simulations showing this universal dynamics is consistent with disease
transmission on scale-free networks and random networks with non-Markovian
transmission dynamics.


  One of the more widely advocated solutions for slowing down the spread of
COVID-19 has been automated contact tracing. Since proximity data can be
collected by personal mobile devices, the natural proposal has been to use this
for automated contact tracing providing a major gain over a manual
implementation. In this work, we study the characteristics of voluntary and
automated contact tracing and its effectiveness for mapping the spread of a
pandemic due to the spread of SARS-CoV-2. We highlight the infrastructure and
social structures required for automated contact tracing to work. We display
the vulnerabilities of the strategy to inadequate sampling of the population,
which results in the inability to sufficiently determine significant contact
with infected individuals. Of crucial importance will be the participation of a
significant fraction of the population for which we derive a minimum threshold.
We conclude that relying largely on automated contact tracing without
population-wide participation to contain the spread of the SARS-CoV-2 pandemic
can be counterproductive and allow the pandemic to spread unchecked. The
simultaneous implementation of various mitigation methods along with automated
contact tracing is necessary for reaching an optimal solution to contain the
pandemic.


  We propose a compartmental mathematical model for the spread of the COVID-19
disease with special focus on the transmissibility of super-spreaders
individuals. We compute the basic reproduction number threshold, we study the
local stability of the disease free equilibrium in terms of the basic
reproduction number, and we investigate the sensitivity of the model with
respect to the variation of each one of its parameters. Numerical simulations
show the suitability of the proposed COVID-19 model for the outbreak that
occurred in Wuhan, China.


  One of the many unresolved questions that revolves around the Covid-19
pandemic is whether local outbreaks can depend on ambient conditions like
temperature and relative humidity. In this paper, we develop a model that tries
to explain and describe the temperature and relative humidity sensitivity of
respiratory droplets and their possible connection in determining viral
outbreaks. The model has two parts. First, we model the growth rate of the
infected population based on a reaction mechanism - the final equations of
which are similar to the well-known SIR model. The advantage of modeling the
pandemic using the reaction mechanism is that the rate constants have sound
physical interpretation. The infection rate constant is derived using collision
rate theory and shown to be a function of the respiratory droplet lifetime. In
the second part, we have emulated the respiratory droplets responsible for
disease transmission as salt solution droplets and computed their evaporation
time accounting for droplet cooling, heat and mass transfer and finally
crystallization of the salt. The model output favourably compares with the
experimentally obtained evaporation characteristics of levitated droplets of
pure water and salt solution, respectively, ensuring fidelity of the model.
Droplet evaporation/desiccation time is indeed dependent on ambient temperature
and relative humidity, considered at both outdoor and indoor conditions. Since
the droplet evaporation time determines the infection rate constant, ambient
temperature and relative humidity are shown to impact the outbreak growth
rates.


  Horizontal gene transfer is an important factor in bacterial evolution that
can act across species boundaries. Yet, we know little about rate and genomic
targets of cross-lineage gene transfer, and about its effects on the recipient
organism's physiology and fitness. Here, we address these questions in a
parallel evolution experiment with two Bacillus subtilis lineages of 7%
sequence divergence. We observe rapid evolution of hybrid organisms: gene
transfer swaps ~12% of the core genome in just 200 generations, and 60% of core
genes are replaced in at least one population. By genomics, transcriptomics,
fitness assays, and statistical modeling, we show that transfer generates
adaptive evolution and functional alterations in hybrids. Specifically, our
experiments reveal a strong, repeatable fitness increase of evolved populations
in the stationary growth phase. By genomic analysis of the transfer statistics
across replicate populations, we infer that selection on HGT has a broad
genetic basis: 40% of the observed transfers are adaptive. At the level of
functional gene networks, we find signatures of negative and positive
selection, consistent with hybrid incompatibilities and adaptive evolution of
network functions. Our results suggest that gene transfer navigates a complex
cross-lineage fitness landscape, bridging epistatic barriers along multiple
high-fitness paths.


  The public goods game is a famous example illustrating the tragedy of the
commons. In this game cooperating individuals contribute to a pool, which in
turn is distributed to all members of the group, including defectors who reap
the same rewards as cooperators without having made a contribution before. The
question is now, how to incentivize group members to all cooperate as it
maximizes the common good. While costly punishment presents one such method,
the cost of punishment still reduces the common good. Here we show how
group-level selection can be such an incentive, and specifically how even
fractions of group-level selection can overcome the benefits defectors receive.
Further, we show how punishment and group-level selection interact. This work
suggests that a redistribution similar to a basic income that is coupled to the
economic success of the entire group could overcome the tragedy of the commons.


  We propose a model for the transmission of perturbations across the amino
acids of a protein represented as an interaction network. The dynamics consists
of a Susceptible-Infected (SI) model based on the Caputo fractional-order
derivative. We find an upper bound to the analytical solution of this model
which represents the worse-case scenario on the propagation of perturbations
across a protein residue network. This upper bound is expressed in terms of
Mittag-Leffler functions of the adjacency matrix of the network of inter-amino
acids interactions. We then apply this model to the analysis of the propagation
of perturbations produced by inhibitors of the main protease of SARS CoV-2. We
find that the perturbations produced by strong inhibitors of the protease are
propagated far away from the binding site, confirming the long-range nature of
intra-protein communication. On the contrary, the weakest inhibitors only
transmit their perturbations across a close environment around the binding
site. These findings may help to the design of drug candidates against this new
coronavirus.


  Mathematical models are traditionally used to analyze the long-term global
evolution of epidemics, to determine the potential and severity of an outbreak,
and to provide critical information for identifying the type of disease
interventions and intensity. One of the widely used mathematical models of
long-term spreading of epidemics are the so-called deterministic compartmental
models (SIR/SEIR type models). One of the main purposes of applying such models
is to assess how the expensive restriction measures imposed by the authorities
(home and social isolation/quarantine, travel restrictions, etc.) can
effectively reduce the control reproduction number of the disease and its
transmission risk. However the classical SIR/SEIR models have been primarily
studied in what may be called stationary case, where the main parameters, the
Transmission rate Beta (reflecting the virus spread by infected individuals)
and the Removed rate Gamma (reflecting the hospitalization/isolation measures)
remain constant during the whole period of interest. Hence, it is important to
extend the classical SIR/SEIR models by creating new ansatzes for the dynamics
of the transmission rates Beta(t) (which we will call further just Beta) and
removed rates Gamma(t) (which we will call further just Gamma). The main
purpose of the present research is to introduce a spline-based SEIR model with
Time-varying Beta and Gamma parameters, or abbreviated TVBG-SEIR model, which
is used to estimate the practical implications of the public health
interventions and measures. We have designed a Tool based on the TVBG-SEIR
model, which may be used as a Decision Support Tool to assist the health
decision- and policy-makers in creating predictive scenarios.


  Following the emergence of a novel coronavirus (SARS-CoV-2) and its spread
outside of China, Europe has experienced large epidemics. In response, many
European countries have implemented unprecedented non-pharmaceutical
interventions including case isolation, the closure of schools and
universities, banning of mass gatherings and/or public events, and most
recently, wide-scale social distancing including local and national lockdowns.
  In this technical update, we extend a semi-mechanistic Bayesian hierarchical
model that infers the impact of these interventions and estimates the number of
infections over time. Our methods assume that changes in the reproductive
number - a measure of transmission - are an immediate response to these
interventions being implemented rather than broader gradual changes in
behaviour. Our model estimates these changes by calculating backwards from
temporal data on observed to estimate the number of infections and rate of
transmission that occurred several weeks prior, allowing for a probabilistic
time lag between infection and death.
  In this update we extend our original model [Flaxman, Mishra, Gandy et al
2020, Report #13, Imperial College London] to include (a) population saturation
effects, (b) prior uncertainty on the infection fatality ratio, (c) a more
balanced prior on intervention effects and (d) partial pooling of the lockdown
intervention covariate. We also (e) included another 3 countries (Greece, the
Netherlands and Portugal).
  The model code is available at
https://github.com/ImperialCollegeLondon/covid19model/
  We are now reporting the results of our updated model online at
https://mrc-ide.github.io/covid19estimates/
  We estimated parameters jointly for all M=14 countries in a single
hierarchical model. Inference is performed in the probabilistic programming
language Stan using an adaptive Hamiltonian Monte Carlo (HMC) sampler.


  We study the SIR epidemiological model, with a variable contagion rate,
applied to the evolution of COVID19 in Cuba. It is highlighted that an increase
in the predictive character depends on understanding the dynamics for the
temporal evolution of the rate of contagion $\beta^*$. A semi-empirical model
for this dynamics is formulated, where reaching $\beta^*\approx0$ due to
isolation is achieved after the mean duration of the disease $\tau=1/\gamma$,
in which the number of infected in the confined families has decreased. It is
considered that $\beta^*(t)$ should have an abrupt decrease on the day of
initiation of confinement and decrease until canceling at the end of the
interval $\tau$. The analysis describes appropriately the infection curve for
Germany. The model is applied to predict an infection curve for Cuba, which
estimates a maximum number of infected as less than 2000 in the middle of May,
depending on the rigor of the isolation. This is suggested by the ratio between
the daily detected cases and the total. We consider the ratio between the
observed and real infected cases (k) less than unity. The low value of k
decreases the maximum obtained when $\beta^*-\gamma>0$. The observed evolution
is independent of k in the linear region. The value of $\beta^*$ is also
studied by time intervals, adjusting to the data of Cuba, Germany and South
Korea. We compare the extrapolation of the evolution of Cuba with the contagion
rate until 16.04.20 with that obtained by a strict quarantine at the end of
April. This model with variable $\beta^*$ correctly describes the observed
infected evolution curves. We emphasize that the desired maximum of the SIR
infected curve is not the maximum standard with constant $\beta^*$, but one
achieved due to quarantine when $\tilde R_0=\beta^*/\gamma<1$. For the
countries controlling the epidemic the maxima are in the region in which SIR
equations are linear.


  Official counts of COVID-19 deaths have been criticized for potentially
including people who did not die of COVID-19 but merely died with COVID-19. I
address that critique by fitting a generalized additive model to weekly counts
of all deaths registered in England and Wales during the 2010s. The model
produces baseline rates of death registrations expected without the COVID-19
pandemic, and comparing those baselines to recent counts of registered deaths
exposes the emergence of excess deaths late in March 2020. By April's end,
England and Wales registered 45,300 $\pm$ 3200 excess deaths of adults aged
45+. Through 22 May, the last day of available all-deaths data, 56,600 $\pm$
4400 were registered (about 53% of which were of men). Both the ONS's
corresponding count of 43,205 death certificates which mention COVID-19, and
the Department of Health and Social Care's count of 33,671 deaths, are
appreciably less, implying that their counting methods have underestimated, not
overestimated, the pandemic's true death toll. If underreporting rates have
held steady during May, about 59,000 direct and indirect COVID-19 deaths might
have been registered through the end of May but not yet publicly reported in
full.


  Accurate forecasts of COVID-19 is central to resource management and building
strategies to deal with the epidemic. We propose a heterogeneous infection rate
model with human mobility for epidemic modeling, a preliminary version of which
we have successfully used during DARPA Grand Challenge 2014. By linearizing the
model and using weighted least squares, our model is able to quickly adapt to
changing trends and provide extremely accurate predictions of confirmed cases
at the level of countries and states of the United States. We show that during
the earlier part of the epidemic, using travel data increases the predictions.
Training the model to forecast also enables learning characteristics of the
epidemic. In particular, we show that changes in model parameters over time can
help us quantify how well a state or a country has responded to the epidemic.
The variations in parameters also allow us to forecast different scenarios such
as what would happen if we were to disregard social distancing suggestions.


  The emergence of SARS-CoV-2 and the coronavirus infectious disease (COVID-19)
has become a pandemic. Social (physical) distancing is a key non-pharmacologic
control measure to reduce the transmission rate of SARS-COV-2, but high-level
adherence is needed. Using daily travel distance and stay-at-home time derived
from large-scale anonymous mobile phone location data provided by Descartes
Labs and SafeGraph, we quantify the degree to which social distancing mandates
have been followed in the U.S. and its effect on growth of COVID-19 cases. The
correlation between the COVID-19 growth rate and travel distance decay rate and
dwell time at home change rate was -0.586 (95% CI: -0.742 ~ -0.370) and 0.526
(95% CI: 0.293 ~ 0.700), respectively. Increases in state-specific doubling
time of total cases ranged from 1.04 ~ 6.86 days to 3.66 ~ 30.29 days after
social distancing orders were put in place, consistent with mechanistic
epidemic prediction models. Social distancing mandates reduce the spread of
COVID-19 when they are followed.


  The outbreak of COVID-19 since Dec. 2019 has caused severe life and economic
damage worldwide, many countries are trapped by medical resource constraints or
absence of targeted therapeutics, and therefore the implement of systemic
policies to block this pandemic should be prioritized. Based on the
transmission mechanisms and physicochemical properties of betacoronaviruses, we
construct a fine-grained transmission dynamics model (ICRD) to forecast the
crucial information of public concern, therein using dynamical coefficients to
quantify the impact of the implement time and intensity of containment policies
on the spread of epidemic. We find that the comprehensive investigation policy
for susceptible population and the quarantine for suspected cases eminently
contribute to reduce casualties during the phase of the dramatic increase of
diagnosed cases. Statistic evidences strongly suggest that society should take
such forceful public health interventions to cut the infection channels in the
initial stage until the pandemic is interrupted.


  Social distancing, often in the form of lockdown, has been adopted by many
countries as a way to contrast the spreading of COVID-19. We discuss the
temporal aspects of social distancing in contrasting an epidemic diffusion. We
argue that a strategy based uniquely on social distancing requires to maintain
the relative measures for a very long time, while a more articulate strategy,
which also uses early detection and prompt isolation, can be both more
efficient on reducing the epidemic peak and allow to relax the social
distancing measures after a much shorter time. We consider in more detail the
situation in Italy, simulating the effect of different strategies through a
recently introduced SIR-type epidemiological model. The short answer to the
question in the title is: "it depends on what else you do".


  Epidemic parameters are estimated through Bayesian inference using the daily
fatality counts in Mumbai during the period from March 31 to April 14. A
doubling time of 5.5 days (median with 95% CrI of 4.6-6.9 days) is observed. In
the SEIR model this gives the basic reproduction rate R_0 of 3.4 (median with
95% CrI of 2.4-4.8). Using as input the infection fatality rate and the
interval between infection and death, the number of infections in Mumbai is
inferred. It is found that the ratio of the number of test positives to the
total infections is 0.13\% (median), implying that tests are currently finding
1 out of 750 cases of infection. After correcting for different testing rates,
this result is compatible with a measurement of the ratio made recently via
serological testing in the USA. From the estimates of the number of infections
we infer that the first COVID-19 cases were seeded in Mumbai between late
December 2019 and early February 2020. provided the doubling times remained
unchanged since then. We remark on some public health implications if the rate
of growth cannot be controlled in about a week.


  Digital contact tracing approaches based on Bluetooth low energy (BLE) have
the potential to efficiently contain and delay outbreaks of infectious diseases
such as the ongoing SARS-CoV-2 pandemic. In this work we propose a novel
machine learning based approach to reliably detect subjects that have spent
enough time in close proximity to be at risk of being infected. Our study is an
important proof of concept that will aid the battery of epidemiological
policies aiming to slow down the rapid spread of COVID-19.


  Background: Rapid testing for an infection is paramount during a pandemic to
prevent continued viral spread and excess morbidity and mortality. This study
aimed to determine whether alternative testing strategies based on sample
pooling can increase the speed and throughput of screening for SARS-CoV-2.
  Methods: A mathematical modelling approach was chosen to simulate six
different testing strategies based on key input parameters (infection rate,
test characteristics, population size, testing capacity etc.). The situations
in five countries (US, DE, UK, IT and SG) currently experiencing COVID-19
outbreaks were simulated to reflect a broad variety of population sizes and
testing capacities. The primary study outcome measurements that were finalised
prior to any data collection were time and number of tests required; number of
cases identified; and number of false positives.
  Findings: The performance of all tested methods depends on the input
parameters, i.e. the specific circumstances of a screening campaign. To screen
one tenth of each country's population at an infection rate of 1% - e.g. when
prioritising frontline medical staff and public workers -, realistic optimised
testing strategies enable such a campaign to be completed in ca. 29 days in the
US, 71 in the UK, 25 in Singapore, 17 in Italy and 10 in Germany (ca. eight
times faster compared to individual testing). When infection rates are
considerably lower, or when employing an optimal, yet logistically more complex
pooling method, the gains are more pronounced. Pool-based approaches also
reduces the number of false positive diagnoses by 50%.
  Interpretation: The results of this study provide a clear rationale for
adoption of pool-based testing strategies to increase speed and throughput of
testing for SARS-CoV-2. The current individual testing approach unnecessarily
wastes valuable time and resources.


  We are frequently faced with a large collection of antibodies, and want to
select those with highest affinity for their cognate antigen. When developing a
first-line therapeutic for a novel pathogen, for instance, we might look for
such antibodies in patients that have recovered. There exist effective
experimental methods of accomplishing this, such as cell sorting and baiting;
however they are time consuming and expensive. Next generation sequencing of B
cell receptor (BCR) repertoires offers an additional source of sequences that
could be tapped if we had a reliable method of selecting those coding for the
best antibodies. In this paper we introduce a method that uses evolutionary
information from the family of related sequences that share a naive ancestor to
predict the affinity of each resulting antibody for its antigen. When combined
with information on the identity of the antigen, this method should provide a
source of effective new antibodies. We also introduce a method for a related
task: given an antibody of interest and its inferred ancestral lineage, which
branches in the tree are likely to harbor key affinity-increasing mutations?
These methods are implemented as part of continuing development of the partis
BCR inference package, available at https://github.com/psathyrella/partis.


  The hydrogen peroxide is present in the living cell at small concentrations
that increase under the action of the heavy ion beams in the process of
anticancer therapy. The interactions of hydrogen peroxide with DNA, proteins
and other biological molecules are poorly understood. In the present work the
competitive binding of the hydrogen peroxide and water molecules with the DNA
double helix backbone has been studied using the molecular dynamics method. The
simulations have been carried out for the DNA double helix in a water solution
with hydrogen peroxide molecules and Na$^{+}$ counterions. The obtained radial
distribution functions of counterions, H$_2$O$_2$ and H$_2$O molecules with
respect to the oxygen atoms of DNA phosphate groups have been used for the
analysis of the formation of different complexes. The calculated mean residence
times show that a hydrogen peroxide molecule stays at least twice as long near
the phosphate group (up to 7 ps) than a water molecule (about 3 ps). The
hydrogen peroxide molecules form more stable complexes with the phosphate
groups of the DNA backbone than water molecules do.


  We revisit the SICA (Susceptible-Infectious-Chronic-AIDS) mathematical model
for transmission dynamics of the human immunodeficiency virus (HIV) with
varying population size in a homogeneously mixing population. We consider SICA
models given by systems of ordinary differential equations and some
generalizations given by systems with fractional and stochastic differential
operators. Local and global stability results are proved for deterministic,
fractional, and stochastic-type SICA models. Two case studies, in Cape Verde
and Morocco, are investigated.


  Two-dimensional partial covariance mass spectrometry (2D-PC-MS) exploits the
inherent fluctuations of fragment ion abundances across a series of tandem mass
spectra, to identify correlated pairs of fragment ions produced along the same
fragmentation pathway of the same parent (e.g. peptide) ion. Here, we apply
2D-PC-MS to the analysis of intact protein ions in a standard linear ion trap
mass analyzer, using the fact that the fragment-fragment correlation signals
are much more specific to bio-molecular sequence than 1D MS/MS signals at the
same mass accuracy and resolution. We show that from the distribution of
signals on a 2D-PC-MS map it is possible to extract the charge state of both
parent and fragment ions without resolving the isotopic envelope. Furthermore,
the 2D map of fragment-fragment correlations naturally reveals the secondary
decomposition pathways of the fragment ions. We access this spectral
information using an adapted version of the Hough transform. We demonstrate the
successful identification of highly charged, intact protein molecules without
the need for high mass resolution. Using this technique we also perform the in
silico deconvolution of the overlapping fragment ion signals from two
co-isolated and co-fragmented intact protein molecules, demonstrating a viable
new method for the concurrent mass spectrometric identification of a mixture of
intact protein ions from the same fragment ion spectrum.


  Motivation: Detailed mechanistic models of biological processes can pose
significant challenges for analysis and parameter estimations due to the large
number of equations used to track the dynamics of all distinct configurations
in which each involved biochemical species can be found. Model reduction can
help tame such complexity by providing a lower-dimensional model in which each
macro-variable can be directly related to the original variables.
  Results: We present CLUE, an algorithm for exact model reduction of systems
of polynomial differential equations by constrained linear lumping. It computes
the smallest dimensional reduction as a linear mapping of the state space such
that the reduced model preserves the dynamics of user-specified linear
combinations of the original variables. Even though CLUE works with nonlinear
differential equations, it is based on linear algebra tools, which makes it
applicable to high-dimensional models. Using case studies from the literature,
we show how CLUE can substantially lower model dimensionality and help extract
biologically intelligible insights from the reduction.
  Availability: An implementation of the algorithm and relevant resources to
replicate the experiments herein reported are freely available for download at
https://github.com/pogudingleb/CLUE.
  Supplementary information: enclosed.


  With declining sequencing costs a promising and affordable tool is emerging
in cancer diagnostics: genomics. By using association studies, genomic variants
that predispose patients to specific cancers can be identified, while by using
tumor genomics cancer types can be characterized for targeted treatment.
However, a severe disparity is rapidly emerging in this new area of precision
cancer diagnosis and treatment planning, one which separates a few genetically
well-characterized populations (predominantly European) from all other global
populations. Here we discuss the problem of population-specific genetic
associations, which is driving this disparity, and present a novel
solution--coordinate-based local ancestry--for helping to address it. We
demonstrate our boosting-based method on whole genome data from divergent
groups across Africa and in the process observe signals that may stem from the
transcontinental Bantu-expansion.


  The COVID-19 pandemic triggered by SARS-CoV-2 is a worldwide health disaster.
Main protease is an attractive drug target among coronaviruses, due to its
vital role in processing the polyproteins that are translated from the viral
RNA. There is presently no exact drug or treatment for this diseases caused by
SARS-CoV-2. In the present study, we report the potential inhibitory activity
of some FDA approved drugs against SARS-CoV-2 main protease by molecular
docking study to investigate their binding affinity in protease active site.
Docking studies revealed that drug Oseltamivir (anti-H1N1 drug), Rifampin
(anti-TB drug), Maraviroc, Etravirine, Indinavir, Rilpivirine (anti-HIV drugs)
and Atovaquone, Quinidine, Halofantrine, Amodiaquine, Tetracylcine,
Azithromycin, hydroxycholoroquine (anti-malarial drugs) among others binds in
the active site of the protease with similar or higher affinity. However, the
in-silico abilities of the drug molecules tested in this study, further needs
to be validated by carrying out in vitro and in vivo studies. Moreover, this
study spreads the potential use of current drugs to be considered and used to
comprise the fast expanding SARS-CoV-2 infection.


  The discovery of SARS-CoV-2, the responsible virus for the Covid-19 epidemic,
has sparked a global health concern with many countries affected. Developing
models that can interpret the epidemic and give common trend parameters are
useful for prediction purposes by other countries that are at an earlier phase
of the epidemic; it is also useful for future planning against viral
respiratory diseases. One model is developed to interpret the fast-growth phase
of the epidemic and another model for an interpretation of the entire data set.
Both models agree reasonably with the data. It is shown by the first model that
during the fast phase, the number of new infected cases depends on the total
number of cases by a power-law relation with a scaling exponent equal to 0.82.
The second model gives a duplication time in the range 1 to 3 days early in the
start of the epidemic, and another parameter alpha = 0.1-0.5) that deviates the
progress of the epidemic from an exponential growth. Our models may be used for
data interpretation and for guiding predictions regarding this disease, e.g.
the onset of the maximum in the number of new cases.


  We investigated the dependence of Covid-19 attack rates on demographic and
socioeconomic factors for the communities in Suffolk County (Long Island, New
York State), presently the 5th most-affected county in the United States.
Confirming the previous observations that minorities are disproportionately
impacted by the Covid-19 disease, we found that the attack rate is strongly
correlated with the minority population proportion, with an alarmingly high
$\sim4$-fold attack rate increase for Black and Hispanic populations.


  This paper develops an agent-level simulation model, termed ALPS, for
simulating the spread of an infectious disease in a confined community. The
mechanism of transmission is agent-to-agent contact, using parameters reported
for Corona COVID-19 pandemic. The main goal of the ALPS simulation is analyze
effects of preventive measures -- imposition and lifting of lockdown norms --
on the rates of infections, fatalities and recoveries. The model assumptions
and choices represent a balance between competing demands of being realistic
and being efficient for real-time inferences. The model provides quantification
of gains in reducing casualties by imposition and maintenance of restrictive
measures in place.


  In this paper, a mathematical model is proposed to analyze the dynamic
behavior of COVID-19. Based on inter-city networked coupling effects, a
fractional-order SEIHDR system with the real-data from 23 January to 18 March,
2020 of COVID-19 is discussed. Meanwhile, hospitalized individuals and the
mortality rates of three types of individuals (exposed, infected and
hospitalized) are firstly taken into account in the proposed model. And
infectivity of individuals during incubation is also considered in this paper.
By applying least squares method and predictor-correctors scheme, the numerical
solutions of the proposed system in the absence of the inter-city network and
with the inter-city network are stimulated by using the real-data from 23
January to $18-m$ March, 2020 where $m$ is equal to the number of prediction
days. Compared with integer-order system ($\alpha=0$), the fractional-order
model without network is validated to have a better fitting of the data on
Beijing, Shanghai, Wuhan, Huanggang and other cities. In contrast to the case
without network, the results indicate that the inter-city network system may be
not a significant case to virus spreading for China because of the lock down
and quarantine measures, however, it may have an impact on cities that have not
adopted city closure. Meanwhile, the proposed model better fits the data from
24 February to 31, March in Italy, and the peak number of confirmed people is
also predicted by this fraction-order model. Furthermore, the existence and
uniqueness of a bounded solution under the initial condition are considered in
the proposed system. Afterwards, the basic reproduction number $R_0$ is
analyzed and it is found to hold a threshold: the disease-free equilibrium
point is locally asymptotically stable when $R_0\le 1$, which provides a
theoretical basis for whether COVID-19 will become a pandemic in the future.


  Few African countries have reported COVID-19 case numbers above $1\,000$ as
of April 18, 2020, with South Africa reporting $3\,034$ cases being hit hardest
in Sub-Saharan Africa. Several African countries, especially South Africa, have
already taken strong non-pharmaceutical interventions that include physical
distancing, restricted economic, educational and leisure activities and reduced
human mobility options. The required strengths and overall effectiveness of
such interventions, however, are debated because of simultaneous but opposing
interests in most African countries: strongly limited health care capacities
and testing capabilities largely conflict with pressured national economies and
socio-economic hardships on the individual level, limiting compliance to
intervention targets. Here we investigate implications of interventions on the
COVID-19 outbreak dynamics, focusing on South Africa before and after the
national lockdown enacted on March 27, 2020. Our analysis shows that initial
exponential growth of existing case numbers is consistent with doubling times
of about $2.5$ days. After lockdown, the growth remains exponential, now with
doubling times of 18 days, but still in contrast to subexponential growth
reported for Hubei/China after lockdown. Moreover, a scenario analysis of a
computational data-driven agent based mobility model for the Nelson Mandela Bay
Municipality (with $1.14$ million inhabitants) hints that keeping current
levels of intervention measures and compliance until the end of April is of
insufficient length and still too weak, too unspecific or too inconsistently
complied with to not overload local intensive care capacity. Yet, enduring,
slightly stronger, more specific interventions combined with sufficient
compliance may constitute a viable option for interventions for regions in
South Africa and potentially for large parts of the African continent.


  On March 10, 2020, Italy imposed a national lockdown to curtail the spread of
COVID-19. Here we estimate that, fourteen days after the implementation of the
strategy, the net reproduction number has dropped below the epidemic threshold
- estimated range 0.4-0.7. Our findings provide a timeline of the effectiveness
of the implemented lockdown, which is relevant for a large number of countries
that followed Italy in enforcing similar measures.


  In this article, we develop a data assimilation procedure to predict the
evolution of epidemics with data uncertainty, with application to the Covid-19
pandemic. We construct a vademecum of solutions by solving the SIR epidemic
model for a set of data neighboring the estimated real (or official) ones. A
reduced basis is constructed from this vademecum through Proper Orthogonal
Decomposition (POD). The reduced POD base is then applied to assimilate the
pandemic data (infected, recovered, deceased) during the period in which data
are known, by a least squares procedure. The fitted curves are then used to
predict the evolution of the pandemic in the next days. Validation tests for
Andalusia region (Spain), Italy and Spain show accurate predictions for 7 days
that improve as the number of assimilated data increases.


  Researchers have shown that even simple empirical models stemming from
biological growth modeling have the potential to provide useful information on
the development and severity of ongoing epidemics since they can be employed as
tools for carrying out projections on the size of the affected population,
timing of turning points, as well as best- and worst-case scenarios.
Nevertheless, they commonly exhibit considerable sensitivity to some input
parameters' variance which results in large fluctuations in the generated
projections, thus rendering predictions difficult and even risky. In this work
we examine a novel meta-projections-based approach which allows us to evaluate
the model's current trends and assess whether generated projections are at a
transient or stable state. Meta-projections can be extracted from graphs of
successive estimations of model's parameters and resulting projections, over a
sequence of days being gradually added to the employed model. In other words,
projections are carried out on truncated time series of cumulative numbers of
confirmed cases with increased lengths at each successive evaluation. This
allows us to trace the values of model parameters over a certain period of time
and examine their trends which may converge to specific values for
settled-growth cases or exhibit a changing or even an erratic behavior for
cases that undergo epidemiological transitions and/or are inappropriately
described by the current model instance(s). We have computed meta-projections
and compared our findings for countries at different stages of the epidemic
with stable or unstable behaviors and increasing or decreasing numbers of
confirmed cases. Our results indicate that meta-projections can aid researchers
in assessing the appropriateness of their relevant models and in effect
decrease the uncertainty in their estimations of an epidemic's severity and
development.