  Introduction: Crohn's disease and ulcerative colitis, both under the umbrella
of inflammatory bowel diseases (IBD), involve many distinct molecular
processes. The difference in their molecular processes is studied by using the
different genes involved in each disease, and it is explored further for drug
targeting and drug repurposing. Methods: The initial set of genes was obtained
by mining published literature and several curated databases. The identified
genes were then subject to Systems and Network analysis to reveal their
molecular processes and shed some light on their pathogenesis. Such
methodologies have identified newer targets and drugs that can be repurposed.
Results: We use a Systems and Network Medicine approach to understand the
mechanism of actions of genes involved in IBD. From an initial set of genes
mined from literature and curated databases, we used the Multi-Steiner Tree
algorithm implemented within the CoVex systems medicine platform to expand each
disease module by incorporating candidate genes with significant connections to
the disease-related seed genes. Such expanded disease modules will identify a
larger set of potential targets and drugs. We used the Closeness Centrality
algorithm implemented within CoVex to search for newer targets and repurposable
drugs. Through a network medicine approach, we provide a mechanistic view of
the diseases and point to newer drugs and targets. Conclusion: We demonstrate
that the Systems and Network Medicine approach is a powerful way to understand
diseases and understand their mechanisms of action.


  The infections and fatalities due to SARS-CoV-2 virus for cases specific to
India have been studied using a deterministic
susceptible-exposed-infected-recovered-dead (SEIRD) compartmental model. One of
the most significant epidemiological parameter, namely the effective
reproduction number of the infection is extracted from the daily growth rate
data of reported infections and it is included in the model with a time
variation. We evaluate the effect of control interventions implemented till now
and estimate the case numbers for infections and deaths averted by these
restrictive measures. We further provide a forecast on the extent of the future
Covid-19 transmission in India and predict the probable numbers of infections
and fatalities under various potential scenarios.


  Summary: Polymerase chain reaction-based assays are the current gold standard
for detecting and diagnosing SARS-CoV-2. However, as SARS-CoV-2 mutates, we
need to constantly assess whether existing PCR-based assays will continue to
detect all known viral strains. To enable the continuous monitoring of
SARS-CoV-2 assays, we have developed a web-based assay validation algorithm
that checks existing PCR-based assays against the ever-expanding genome
databases for SARS-CoV-2 using both thermodynamic and edit-distance metrics.
The assay screening results are displayed as a heatmap, showing the number of
mismatches between each detection and each SARS-CoV-2 genome sequence. Using a
mismatch threshold to define detection failure, assay performance is summarized
with the true positive rate (recall) to simplify assay comparisons.
Availability: https://covid19.edgebioinformatics.org/#/assayValidation.
Contact: Jason Gans (jgans@lanl.gov) and Patrick Chain (pchain@lanl.gov)


  Generative adversarial networks (GANs) are pairs of artificial neural
networks that are trained one against each other. The outputs from a generator
are mixed with the real-world inputs to the discriminator and both networks are
trained until an equilibrium is reached, where the discriminator cannot
distinguish generated inputs from real ones. Since their introduction, GANs
have allowed for the generation of impressive imitations of real-life films,
images and texts, whose fakeness is barely noticeable to humans. Despite their
impressive performance, training GANs remains to this day more of an art than a
reliable procedure, in a large part due to training process stability.
Generators are susceptible to mode dropping and convergence to random patterns,
which have to be mitigated by computationally expensive multiple restarts.
Curiously, GANs bear an uncanny similarity to a co-evolution of a pathogen and
its host's immune system in biology. In a biological context, the majority of
potential pathogens indeed never make it and are kept at bay by the hots'
immune system. Yet some are efficient enough to present a risk of a serious
condition and recurrent infections. Here, we explore that similarity to propose
a more robust algorithm for GANs training. We empirically show the increased
stability and a better ability to generate high-quality images while using less
computational power.


  Human to human transmissible infectious diseases spread in a population using
human interactions as its transmission vector. The early stages of such an
outbreak can be modeled by a graph whose edges encode these interactions
between individuals, the vertices. This article attempts to account for the
case when each individual entails in different kinds of interactions which have
therefore different probabilities of transmitting the disease. The majority of
these results can be also stated in the language of percolation theory.
  The main contributions of the article are: (1) Extend to this setting some
results which were previously known in the case when each individual has only
one kind of interactions. (2) Find an explicit formula for the basic
reproduction number $R_0$ which depends only on the probabilities of
transmitting the disease along the different edges and the first two moments of
the degree distributions of the associated graphs. (3) Motivated by the recent
Covid-19 pandemic, we use the framework developed to compute the $R_0$ of a
model disease spreading in populations whose trees and degree distributions are
adjusted to several different countries. In this setting, we shall also compute
the probability that the outbreak will not lead to an epidemic. In all cases we
find such probability to be very low if no interventions are put in place.


  In this paper, it is shown that the Moyal distribution is an excelent tool to
study the SARS-Cov-II (Covid-19) epidemiological associated curves and its
propagation. The Moyal parameters give all the information to describe the form
and the impact of the illness outbreak in the different affected countries and
its global impact. We checked that the Moyal distribution can accurately fit
the daily report of {\it{new confirmed cases of infected people}} (NCC) per
country, in that places where the contagion is reaching their final phase,
describing the beginning, the most intense phase and the descend of the
contagion, simultaneously . In order to achieve the purpose of this work, it is
important to work with a complete and well compilated set of the data to be
used to fit the curves. Data from European countries like France, Spain, Italy
Belgium, Sweden, United Kingdom, Denmark and others like USA and China, have
been used. Also, the correlation between the parameters of the Moyal
distribution fitting and the general public health measures imposed in each
country, have been discussed. A relation between those policies and the
features of the Moyal distribution, in terms of their parameters and critical
points, is shown; from that, it can be seen that the knowledge of the time
evolution of the epidemiological curve, their critical points, superposition
properties and rates of the rising and the ending, could help to find a way to
estimate the efficiency of social distancing measures, imposed in each country,
and anticipate the different phases of the pandemic.


  Quarantine measure is a commonly used non-pharmaceutical intervention during
the outbreak of infectious diseases. A key problem for implementing quarantine
measure is to determine the duration of quarantine. In this paper, a policy
with optimal quarantine duration is developed. The policy suggests different
quarantine durations for every individual with different characteristic. The
policy is optimal in the sense that it minimizes the average quarantine
duration of uninfected people with the constraint that the probability of
symptom presentation for infected people attains the given value closing to 1.
The optimal solution for the quarantine duration is obtained and estimated by
some statistic methods with application to analyzing COVID-19 data.


  COVID-19 has forced quarantine measures in several countries across the
world. These measures have proven to be effective in significantly reducing the
prevalence of the virus. To date, no effective treatment or vaccine is
available. In the effort of preserving both public health as well as the
economical and social textures, France and Italy governments have partially
released lockdown measures. Here we extrapolate the long-term behavior of the
epidemics in both countries using a Susceptible-Exposed-Infected-Recovered
(SEIR) model where parameters are stochastically perturbed to handle the
uncertainty in the estimates of COVID-19 prevalence. Our results suggest that
uncertainties in both parameters and initial conditions rapidly propagate in
the model and can result in different outcomes of the epidemics leading or not
to a second wave of infections. Using actual knowledge, asymptotic estimates of
COVID-19 prevalence can fluctuate of order of ten millions units in both
countries.


  Species introductions promote secondary contacts between taxa with long
histories of allopatric divergence. Anthropogenic contact zones thus offer
valuable contrasts to speciation studies in natural systems where past spatial
isolations may have been brief or intermittent. Investigations of anthropogenic
hybridization are rare for marine animals, which have high fecundity and high
dispersal ability, characteristics that contrast to most terrestrial animals.
Genomic studies indicate that gene flow can still occur after millions of years
of divergence, as illustrated by invasive mussels and tunicates. In this
context, we highlight three issues: 1) the effects of high propagule pressure
and demographic asymmetries on introgression directionality, 2) the role of
hybridization in preventing introduced species spread, and 3) the importance of
postzygotic barriers in maintaining reproductive isolation. Anthropogenic
contact zones offer evolutionary biologists unprecedented large scale
hybridization experiments. In addition to breaking the highly effective
reproductive isolating barrier of spatial segregation, they allow researchers
to explore unusual demographic contexts with strong asymmetries. The outcomes
are diverse from introgression swamping to strong barriers to gene flow, and
lead to local containment or widespread invasion. These outcomes should not be
neglected in management policies of marine invasive species.


  Spatially-resolved omics-data enable researchers to precisely distinguish
cell types in tissue and explore their spatial interactions, enabling deep
understanding of tissue functionality. To understand what causes or
deteriorates a disease and identify related biomarkers, clinical researchers
regularly perform large-scale cohort studies, requiring the comparison of such
data at cellular level. In such studies, with little a-priori knowledge of what
to expect in the data, explorative data analysis is a necessity. Here, we
present an interactive visual analysis workflow for the comparison of cohorts
of spatially-resolved omics-data. Our workflow allows the comparative analysis
of two cohorts based on multiple levels-of-detail, from simple abundance of
contained cell types over complex co-localization patterns to individual
comparison of complete tissue images. As a result, the workflow enables the
identification of cohort-differentiating features, as well as outlier samples
at any stage of the workflow. During the development of the workflow, we
continuously consulted with domain experts. To show the effectiveness of the
workflow we conducted multiple case studies with domain experts from different
application areas and with different data modalities.


  In times of outbreaks, an essential requirement for better monitoring is the
evaluation of the number of undiagnosed infected individuals. An accurate
estimate of this fraction is crucial for the assessment of the situation and
the establishment of protective measures. In most current studies using
epidemics models, the total number of infected is either approximated by the
number of diagnosed individuals or is dependent on the model parameters and
assumptions, which are often debated. We here study the relationship between
the fraction of diagnosed infected out of all infected, and the fraction of
infected with known contaminator out of all diagnosed infected. We show that
those two are approximately the same in exponential models and across most
models currently used in the study of epidemics, independently of the model
parameters. As an application, we compute an estimate of the effective number
of infected by the SARS-CoV-2 virus in various countries.


  Measures to reduce the impact of the COVID-19 pandemic require a mix of
logistic, political and social capacity. Depending on the country, different
approaches to increase hospitalization capacity or to properly apply lock-downs
are observed. In order to better understand the impact of these measures we
have developed a compartmental model which, on the one hand allows to calibrate
the reduction of movement of people within and among different areas, and on
the other hand it incorporates a hospitalization dynamics that differentiates
the available kinds of treatment that infected people can receive. By bounding
the hospitalization capacity, we are able to study in detail the interplay
between mobility and hospitalization capacity.


  The Bacille Calmette-Gu\'erin (BCG) tuberculosis vaccine has immunity
benefits against respiratory infections. Accordingly, it has been hypothesized
that it may have a protective effect against COVID-19. Recent research found
that countries with universal Bacillus Calmette-Gu\'erin (BCG) childhood
vaccination policies tend to be less affected by the COVID-19 pandemic.
However, such ecological studies are biased by numerous confounders. Instead,
this paper takes advantage of a rare nationwide natural experiment that took
place in Sweden in 1975, where discontinuation of newborns BCG vaccination led
to a dramatic fall of the BCG coverage rate from 92% to 2% , thus allowing us
to estimate the BCG's effect without all the biases associated with
cross-country comparisons. Numbers of COVID-19 cases and hospitalizations were
recorded for birth cohorts born just before and just after that change,
representing 1,026,304 and 1,018,544 individuals, respectively. We used
regression discontinuity to assess the effect of BCG vaccination on Covid-19
related outcomes. This method used on such a large population allows for a high
precision that would be hard to achieve using a classical randomized controlled
trial. The odds ratio for Covid-19 cases and Covid-19 related hospitalizations
were 0.9997 (CI95: [0.8002-1.1992]) and 1.1931 (CI95: [0.7558-1.6304]),
respectively. We can thus reject with 95\% confidence that universal BCG
vaccination reduces the number of cases by more than 20% and the number of
hospitalizations by more than 24%. While the effect of a recent vaccination
must be evaluated, we provide strong evidence that receiving the BCG vaccine at
birth does not have a protective effect against COVID-19.


  During the COVID-19 pandemic, the scientific community developed predictive
models to evaluate potential governmental interventions. However, the analysis
of the effects these interventions had is less advanced. Here, we propose a
data-driven framework to assess these effects retrospectively. We use a
regularized regression to find a parsimonious model that fits the data with the
least changes in the Rt parameter. Then, we postulate each jump in Rt as the
effect of an intervention. Following the do-operator prescriptions, we simulate
the counterfactual case by forcing Rt to stay at the pre-jump value. We then
attribute a value to the intervention from the difference between true
evolution and simulated counterfactual. We show that the recommendation to use
facemasks for all activities would reduce the number of cases by 170000 (95% CI
160000 to 180000) in Connecticut, Massachusetts, and New York State. The
framework presented here might be used in any case where cause and effects are
sparse in time.


  The outbreak of Coronavirus Disease 2019 (COVID-19) is an ongoing pandemic
affecting over 200 countries and regions. Inference about the transmission
dynamics of COVID-19 can provide important insights into the speed of disease
spread and the effects of mitigation policies. We develop a novel Bayesian
approach to such inference based on a probabilistic compartmental model using
data of daily confirmed COVID-19 cases. In particular, we consider a
probabilistic extension of the classical susceptible-infectious-recovered
model, which takes into account undocumented infections and allows the
epidemiological parameters to vary over time. We estimate the disease
transmission rate via a Gaussian process prior, which captures nonlinear
changes over time without the need of specific parametric assumptions. We
utilize a parallel-tempering Markov chain Monte Carlo algorithm to efficiently
sample from the highly correlated posterior space. Predictions for future
observations are done by sampling from their posterior predictive
distributions. Performance of the proposed approach is assessed using simulated
datasets. Finally, our approach is applied to COVID-19 data from four states of
the United States: Washington, New York, California, and Illinois. An R package
BaySIR is made available at https://github.com/tianjianzhou/BaySIR for the
public to conduct independent analysis or reproduce the results in this paper.


  This article describes a simple Susceptible Infected Recovered (SIR) model
fitting with COVID-19 data for the month of march 2020 in New York (NY) state.
The model is a classical SIR, but is non-autonomous; the rate of susceptible
people becoming infected is adjusted over time in order to fit the available
data. The death rate is also secondarily adjusted. Our fitting is made under
the assumption that due to limiting number of tests, a large part of the
infected population has not been tested positive. In the last part, we extend
the model to take into account the daily fluxes between New Jersey (NJ) and NY
states and fit the data for both states. Our simple model fits the available
data, and illustrates typical dynamics of the disease: exponential increase,
apex and decrease. The model highlights a decrease in the transmission rate
over the period which gives a quantitative illustration about how lockdown
policies reduce the spread of the pandemic. The coupled model with NY and NJ
states shows a wave in NJ following the NY wave, illustrating the mechanism of
spread from one attractive hot spot to its neighbor. }


  Microtubule (MT) radial arrays or asters establish the internal topology of a
cell by interacting with organelles and molecular motors. We proceed to
understand the general pattern forming potential of aster-motor systems using a
computational model of multiple MT asters interacting with motors in a cellular
confinement. In this model dynein motors are attached to the cell cortex and
plus-ended motors resembling kinesin-5 diffuse in the cell interior. The
introduction of 'noise' in the form of MT length fluctuations spontaneously
results in the emergence of coordinated, achiral vortex-like rotation of
asters. The coherence and persistence of rotation requires a threshold density
of both cortical dyneins and coupling kinesins, while the onset of rotation if
diffusion-limited with relation to cortical dynein mobility. The coordinated
rotational motion arises due to the resolution of the 'tug-of-war' of the
rotational component due to cortical motors by 'noise' in the form of MT
dynamic instability, and such transient symmetry breaking is amplified by local
coupling by kinesin-5 complexes. The lack of widespread aster rotation across
cell types suggests biophysical mechanisms that suppress such intrinsic
dynamics may have evolved. This model is analogous to more general models of
locally coupled self-propelled particles (SPP) that spontaneously undergo
collective transport in presence of 'noise' that have been invoked to explain
swarming in birds and fish. However, the aster-motor system is distinct from
SPP models with regard to particle density and 'noise' dependence, providing a
set of experimentally testable predictions for a novel sub-cellular pattern
forming system.


  Until a vaccine or therapy is found against the SARS-CoV-2 coronavirus,
reaching herd immunity appears to be the only mid-term option. However, if the
number of infected individuals decreases and eventually fades only beyond this
threshold, a significant proportion of susceptible may still be infected until
the epidemic is over. A containment strategy is likely the best policy in the
worst case where no vaccine or therapy is found. In order to keep the number of
newly infected persons to a minimum, a possible strategy is to apply strict
containment measures, so that the number of susceptible individuals remains
close to herd immunity. Such an action is unrealistic since containment can
only last for a finite amount of time and is never total. In this article,
using a classical SIR model, we determine the (partial or total) containment
strategy on a given finite time interval that maximizes the number of
susceptible individuals over an infinite horizon, or equivalently that
minimizes the total infection burden during the curse of the epidemic. The
existence and uniqueness of the optimal strategy is proved and the latter is
fully characterized. If applicable in practice, such a strategy would lead
theoretically to an increase by 30% of the proportion of susceptible on an
infinite horizon, for a containment level corresponding to the sanitary
measures put in place in France from March to May 2020. We also analyze the
minimum intervention time to reach a fixed distance from herd immunity, and
show the relationship with the previous problem. Simulations are provided that
illustrate and validate the theoretical results.


  The case fatality ratio (CFR) for COVID-19 is difficult to estimate. One
difficulty is due to ignoring or overestimating time delay between reporting
and death. We claim that all of these cause large errors and artificial time
dependence of the CFR. We find that for each country, there is a unique value
of the time lag between reported cases and deaths versus time, that yields the
optimal correlation between them is a specific sense. We find that the
resulting corrected CFR (deaths shifted back by this time lag, divided by
cases) is actually constant over many months, for many countries, but also for
the entire world. This optimal time lag and constant CFR for each country can
be found through a simple data driven algorithm. The traditional CFR (ignoring
time lag) is spuriously time-dependent and its evolution is hard to quantify.
Our corrected CFR is constant over time, therefore an important index of the
pandemic in each country, and can be inferred from data earlier on,
facilitating improved early estimates of COVID-19 mortality.


  We analyze the late-time relaxation dynamics for a general contagion model.
In this model, nodes are either active or failed. Active nodes can fail either
"spontaneously" at any time or "externally" if their neighborhoods are
sufficiently damaged. Failed nodes may always recover spontaneously. At late
times, the breaking of time-translation-invariance is a necessary condition for
physical aging. We observe that time-translational invariance is lost for
initial conditions that lie between the basins of attraction of the model's two
stable stationary states. Based on corresponding mean-field predictions, we
characterize the observed model behavior in terms of a phase diagram spanned by
the fractions of spontaneously and externally failed nodes. For the square
lattice, the phases in which the dynamics approaches one of the two stable
stationary states are not linearly separable due to spatial correlation
effects. Our results provide new insights into aging and relaxation phenomena
that are observable in a model of social contagion processes.


  We study equilibrium shapes and shape transformations of a confined
semiflexible chain inside a soft lipid tubule using simulations and continuum
theories. The deformed tubular shapes and chain conformations depend on the
relative magnitude of their bending moduli. We characterise the collapsed
macromolecular shapes by computing statistical quantities that probe the
polymer properties at small length scales and report a prolate to toroidal coil
transition for stiff chains. Deformed tubular shapes, calculated using elastic
theories, agree with simulations. In conjunction with scattering studies, our
work may provide a mechanistic understanding of gene encapsulation in soft
structures.


  A compartmental epidemic model is proposed to predict the Covid-19 virus
spread. It considers: both detected and undetected infected populations,
medical quarantine and social sequestration, release from sequestration, plus
possible reinfection. The coefficients in the model are evaluated by fitting to
empirical data for eight US states: Arizona, California, Florida, Illinois,
Louisiana, New Jersey, New York State, and Texas. Together these states make up
43% of the US population; some of these states appear to have handled their
initial outbreaks well, while others appear to be emerging hotspots. The
evolution of Covid-19 is fairly similar among the states: variations in contact
and recovery rates remain below 5%; however, not surprisingly, variations are
larger in death rate, reinfection rate, stay-at-home effect, and release rate
from sequestration. The results reveal that outbreaks may have been well
underway in several states before first detected and that California might have
seen more than one influx of the pandemic. Our projections based on the current
situation indicate that Covid-19 will become endemic, spreading for more than
two years. Should states fully relax stay-at-home orders, most states may
experience a secondary peak in 2021. If lockdowns had been kept in place, the
number of Covid-19 deaths so far could have been significantly lower in most
states that opened up. Additionally, our model predicts that decreasing contact
rate by 10%, or increasing testing by approximately 15%, or doubling lockdown
compliance (from the current $\sim$ 15% to $\sim$ 30%) will eradicate
infections in the state of Texas within a year. Extending our fits for all of
the US states, we predict about 11 million total infections (including
undetected), 8 million cumulative confirmed cases, and 630,000 cumulative
deaths by November 1, 2020.


  Given the growth of urbanization and emerging pandemic threats, more
sophisticated models are required to understand disease propagation and
investigate the impacts of intervention strategies across various city types.
We introduce a fully mechanistic, activity-based and highly spatio-temporally
resolved epidemiological model which leverages on person-trajectories obtained
from integrated mobility demand and supply models in full-scale cities.
Simulating COVID-19 evolution in two full-scale cities with representative
synthetic populations and mobility patterns, we analyze activity-based contact
networks. We observe that transit contacts are scale-free in both cities, work
contacts are Weibull distributed, and shopping or leisure contacts are
exponentially distributed. We also investigate the impact of the transit
network, finding that its removal dampens disease propagation, while work is
also critical to post-peak disease spreading. Our framework, validated against
existing case and mortality data, demonstrates the potential for tracking and
tracing, along with detailed socio-demographic and mobility analyses of
epidemic control strategies.


  Protein chains of the (FG)$_n$ ($n \simeq$ 300) type cap the cytoplasmatic
side of the nucleopore complex, which connects the nucleus to the remainder of
an eukaryotic cell. We study the properties of three fundamental polymer models
that represent these filaments using Monte Carlo computer simulations. Random
walks and the worm like chain model cannot account for the unusual size
selectivity of the pore, while a two-dimensional arrangement of intrinsically
disordered block copolymers with a high content of $\alpha$-helices is in
agreement with the biochemical findings. We predict a linear increase of the
free energy barrier of protein transport through the pore with increasing
protein diameter, which can be probed experimentally using atomic force
microscopy or optical tweezers.


  Antibodies with high titer and affinity to small molecule are critical in the
field for the development of vaccines against drugs of abuse, antidotes to
toxins and immunoassays for compounds. However, little is known regarding how
properties of small molecule influence and which chemical descriptor could
indicate the degree of the antibody response. Based on our previous study, we
designed and synthesized two groups of small molecules, called haptens, with
varied hydrophobicities to investigate the relationship between properties of
small molecules and antibody response in term of titer and affinity. We found
that the magnitude of the antibody response is positively correlated with the
degree of molecular hydrophobicity and related chemical descriptors. This study
provides insight into the immunological characteristics of small molecules
themselves and useful clues to produce high quality antibodies against small
molecules.


  Far-UV circular dichroism (CD) spectroscopy provides a rapid, sensitive,
nondestructive tool to analyze protein conformation by monitoring secondary
structure composition. Originally intended for educational purposes, a
spreadsheet-based program that implements a rudimentary routine for fitting and
simulating far-UV protein CD spectra became quite used in research papers too,
as it allowed very quick deconvolution of spectra into secondary structure
compositions and easy simulation of spectra expected for defined secondary
structure contents. To make such software more readily available, I present
here an online version that runs directly on web browsers allowing even faster
analyses on any modern device and without the need of any spreadsheet or third
programs. The new version further extends the original capabilities to fit and
simulate alpha, beta and random coil contents now including also beta turns; it
allows to quickly select the effective spectral window used for fitting and
enables interactive exploration of the effects of changes in secondary
structure composition on the resulting spectra. The web app allows ubiquitous
implementation in biophysics courses for example right on student smartphones,
and seamless, rapid tests in research settings before moving to more advanced
analysis programs (a few proposed here too). The web app is freely available
without registration at
http://lucianoabriata.altervista.org/jsinscience/cd/cd3.html


  We investigate a modified spatial stochastic Lotka-Volterra formulation of
the rock-paper-scissors model using off-lattice stochastic simulations. In this
model one of the species moves preferentially in a specific direction -- the
level of preference being controlled by a noise strength parameter $\eta \in
[0, 1]$ ($\eta = 0$ and $\eta = 1$ corresponding to total preference and no
preference, respectively) -- while the other two species have no referred
direction of motion. We study the behaviour of the system starting from random
initial conditions, showing that the species with asymmetric mobility has
always an advantage over its predator. We also determine the optimal value of
the noise strength parameter which gives the maximum advantage to that species.
Finally, we find that the critical number of individuals, below which the
probability of extinction becomes significant, decreases as the noise level
increases, thus showing that the addition of a preferred mobility direction
studied in the present paper does not favour coexistence.


  In this paper we propose an epidemiological model for the spread of COVID-19.
The dynamics of the spread is based on four fundamental categories of people in
a population: Tested and infected, Non-Tested but infected, Tested but not
infected, and non-Tested and not infected. The model is based on two levels of
dynamics of spread in the population: at local level and at the global level.
The local level growth is described with data and parameters which include
testing statistics for COVID-19, preventive measures such as nationwide
lockdown, and the migration of people across neighboring locations. In the
context of India, the local locations are considered as districts and migration
or traffic flow across districts are defined by normalized edge weight of the
metapopulation network of districts which are infected with COVID-19. Based on
this local growth, state level predictions for number of people tested with
COVID-19 positive are made. Further, considering the local locations as states,
prediction is made for the country level. The values of the model parameters
are determined using grid search and minimizing an error function while
training the model with real data. The predictions are made based on the
present statistics of testing, and certain linear and log-linear growth of
testing at state and country level. Finally, it is shown that the spread can be
contained if number of testing can be increased linearly or log-linearly by
certain factors along with the preventive measures in near future. This is also
necessary to prevent the sharp growth in the count of infected and to get rid
of the second wave of pandemic.


  The investigation focuses on an on-off protocol for controlling the COVID-19
widespread. The protocol establishes a working period of 4 days for all the
citizens, followed by 8 days of lock-down. We further propose splitting people
into smaller groups that undergo the on-off protocol, but shifted in time. This
procedure is expected to regularize the overall economic activity. Our results
show that either the on-of protocol and the splitting into groups reduces the
amount of infected people. However, the latter seems to be better for economic
reasons. Our simulations further show that the start-up time is a key issue for
the success of the implementation.


  A large percentage of the deaths from COVID-19 occur among residents of
long-term care facilities. There are two possible reasons for this phenomenon.
First, the structural features of such settings may lead to death.
Alternatively, it is possible that individuals in these facilities are in
poorer health than those living elsewhere, and that these individuals would
have died even if they had not been in these facilities. Our findings show
that, controlling for the population density and the percentage of older adults
in the population, there is a significant positive association between the
number of long-term care beds per capita and COVID-19 mortality rates. This
finding provides support for the claim that long-term care living arrangements
(of older people) are a significant risk factor for dying from COVID-19.


  We consider two-player iterated survival games in which players may switch
from a more cooperative behavior to a less cooperative one at some step of the
game. Payoffs are survival probabilities and lone individuals have to finish
the game on their own. We explore the potential of these games to support
cooperation, focusing on the case in which each single step is a Prisoner's
Dilemma. We find that incentives for or against cooperation depend on the
number of defections at the end of the game, as opposed to the number of steps
in the game. Broadly, cooperation is supported when the survival prospects of
lone individuals are relatively bleak. Specifically, we find three critical
values or cutoffs for the loner survival probability which, in concert with
other survival parameters, determine the incentives for or against cooperation.
One cutoff determines the existence of an optimal number of defections against
a fully cooperative partner, one determines whether additional defections
eventually become disfavored as the number of defections by the partner
increases, and one determines whether additional cooperations eventually become
favored as the number of defections by the partner increases. We obtain
expressions for these switch-points and for optimal numbers of defections
against partners with various strategies. These typically involve small numbers
of defections even in very long games. We show that potentially long stretches
of equilibria may exist, in which there is no incentive to defect more or
cooperate more. We describe how individuals find equilibria in best-response
walks among strategies, and establish that evolutionary stability requires
there be just one such equilibrium. Otherwise, equilibria are not protected
against invasion by strategies with fewer defections.


  The history of the Earth has been marked by major ecological transitions,
driven by metabolic innovation, that radically reshaped the composition of the
oceans and atmosphere. The nature and magnitude of the earliest transitions,
hundreds of million years before photosynthesis evolved, remain poorly
understood. Using a novel ecosystem-planetary model, we find that
pre-photosynthetic methane-cycling microbial ecosystems are much less
productive than previously thought. In spite of their low productivity, the
evolution of methanogenic metabolisms strongly modifies the atmospheric
composition, leading to a warmer but less resilient climate. As the abiotic
carbon cycle responds, further metabolic evolution (anaerobic methanotrophy)
may feed back to the atmosphere and destabilize the climate, triggering a
transient global glaciation. Although early metabolic evolution may cause
strong climatic instability, a low CO:CH4 atmospheric ratio emerges as a robust
signature of simple methane-cycling ecosystems on a globally reduced planet
such as the late Hadean/early Archean Earth.


  SARS-CoV-2, the causative agent of the disease known as Covid-19, has so far
reported around 3,435,000 cases of human infections, including more than
239,000 deaths in 187 countries, with no effective treatment currently
available. For this reason, it is necessary to explore new approaches for the
development of a drug capable of inhibiting the entry of the virus into the
host cell. Therefore, this work includes the exploration of potential
inhibitory compounds for the Spike protein of SARS-CoV-2 (PDB ID: 6VSB), which
were obtained from The Patogen Box. Later, they were filtered through virtual
screening and molecular docking techniques, thus obtaining a top of 1000
compounds, which were used against a binding site located in the Receptor
Binding Domain (RBD) and a cryptic site located in the N-Terminal Domain (NTD),
resulting in good pharmaceutical targets for the blocking the infection. From
the top 1000, the best compound (TCMDC-124223) was selected for the binding
site. It interacts with specific residues that intervene in the recognition and
subsequent entry into the host cell, resulting in a more favorable binding free
energy in comparison to the control compounds (Hesperidine and Emodine). In the
same way, the compound TCMDC-133766 was selected for the cryptic site. These
identified compounds are potential inhibitors that can be used for the
development of new drugs that allow effective treatment for the disease.


  Quantifying the influence of microscopic details on the dynamics of
development of the overall structure of a filamentous network is important in a
number of biologically relevant contexts, but it is not obvious what order
parameters can be used to adequately describe this complex process. In this
paper, we investigated the role of multivalent actin-binding proteins (ABPs) in
reorganizing actin filaments into higher-order complex networks via a computer
model of semiflexible filaments. We characterize the importance of local
connectivity among actin filaments as well as the global features of actomyosin
networks. We first map the networks into local graph representations and then,
using principles from network-theory order parameters, combine properties from
these representations to gain insight on the heterogeneous morphologies of
actomyosin networks at a global level. We find that ABPs with a valency greater
than two promote filament bundles and large filament clusters to a much greater
extent than bivalent multilinkers. We also show that active myosin-like motor
proteins promote the formation of dendritic branches from a stalk of actin
bundles. Our work motivates future studies to embrace network theory as a tool
to characterize complex morphologies of actomyosin detected by experiments,
leading to a quantitative understanding of the role of ABPs in manipulating the
self-assembly of actin filaments into unique architectures that underlie the
structural scaffold of a cell relating to its mobility and shape.


  Global coronavirus disease pandemic (COVID-19) caused by newly identified
SARS- CoV-2 coronavirus continues to claim the lives of thousands of people
worldwide. The unavailability of specific medications to treat COVID-19 has led
to drug repositioning efforts using various approaches, including computational
analyses. Such analyses mostly rely on molecular docking and require the 3D
structure of the target protein to be available. In this study, we utilized a
set of machine learning algorithms and trained them on a dataset of
RNA-dependent RNA polymerase (RdRp) inhibitors to run inference analyses on
antiviral and anti-inflammatory drugs solely based on the ligand information.
We also performed virtual screening analysis of the drug candidates predicted
by machine learning models and docked them against the active site of SARS-
CoV-2 RdRp, a key component of the virus replication machinery. Based on the
ligand information of RdRp inhibitors, the machine learning models were able to
identify candidates such as remdesivir and baloxavir marboxil, molecules with
documented activity against RdRp of the novel coronavirus. Among the other
identified drug candidates were beclabuvir, a non-nucleoside inhibitor of the
hepatitis C virus (HCV) RdRp enzyme, and HCV protease inhibitors paritaprevir
and faldaprevir. Further analysis of these candidates using molecular docking
against the SARS-CoV-2 RdRp revealed low binding energies against the enzyme
active site. Our approach also identified anti-inflammatory drugs lupeol,
lifitegrast, antrafenine, betulinic acid, and ursolic acid to have potential
activity against SARS-CoV-2 RdRp. We propose that the results of this study are
considered for further validation as potential therapeutic options against
COVID-19.


  We sample aggravated cases following age-structured probabilities from
confirmed cases and use ICU occupation data to find a subnotification factor. A
logistic fit is then employed to project the progression of the COVID-19
epidemic with plateau scenarios taken from locations that have reached this
stage. Finally, the logistic curve found is corrected by the subnotification
factor and sampled to project the future demand for ICU beds.


  India is one of the worst affected countries by the Covid-19 pandemic at
present. We studied publicly available data of the Covid-19 patients in India
and analyzed possible impacts of quarantine and social distancing within the
stochastic framework of the SEQIR model to illustrate the controlling strategy
of the pandemic. Our simulation results clearly show that proper quarantine and
social distancing should be maintained from an early time just at the start of
the pandemic and should be continued till its end to effectively control the
pandemic. This calls for a more socially disciplined lifestyle in this
perspective in future. The demographic stochasticity, which is quite visible in
the system dynamics, has a critical role in regulating and controlling the
pandemic.


  While measurement advances now allow extensive surveys of gene activity
(large numbers of genes across many samples), interpretation of these data is
often confounded by noise -- expression counts can differ strongly across
samples due to variation of both biological and experimental origin.
Complimentary to perturbation approaches, we extract functionally related
groups of genes by analyzing the standing variation within a sampled
population. To distinguish biologically meaningful patterns from
uninterpretable noise, we focus on correlated variation and develop a novel
density-based clustering approach that takes advantage of a percolation
transition generically arising in random, uncorrelated data. We apply our
approach to two contrasting RNA sequencing data sets that sample individual
variation -- across single cells of fission yeast and whole animals of C.
elegans worms -- and demonstrate robust applicability and versatility in
revealing correlated gene clusters of diverse biological origin, including cell
cycle phase, development/reproduction, tissue-specific functions, and feeding
history. Our technique exploits generic features of noisy high-dimensional data
and is applicable, beyond gene expression, to feature-rich data that sample
population-level variability in the presence of noise.


  A major ambition of systems science is to uncover the building blocks of any
biological network to decipher how cellular function emerges from their
interactions. Here, we introduce a graph representation of the information flow
in these networks as a set of input trees, one for each node, which contains
all pathways along which information can be transmitted in the network. In this
representation, we find remarkable symmetries in the input trees that
deconstruct the network into functional building blocks called fibers. Nodes in
a fiber have isomorphic input trees and thus process equivalent dynamics and
synchronize their activity. Each fiber can then be collapsed into a single
representative base node through an information-preserving transformation
called 'symmetry fibration', introduced by Grothendieck in the context of
algebraic geometry. We exemplify the symmetry fibrations in gene regulatory
networks and then show that they universally apply across species and domains
from biology to social and infrastructure networks. The building blocks are
classified into topological classes of input trees characterized by integer
branching ratios and fractal golden ratios of Fibonacci sequences representing
cycles of information. Thus, symmetry fibrations describe how complex networks
are built from the bottom up to process information through the synchronization
of their constitutive building blocks.


  The new corona virus disease -- COVID-2019 -- is rapidly spreading through
the world. The availability of unbiased timely statistics of trends in disease
events are a key to effective responses. But due to reporting delays, the most
recently reported numbers are frequently underestimating of the total number of
infections, hospitalizations and deaths creating an illusion of a downward
trend. Here we describe a statistical methodology for predicting true daily
quantities and their uncertainty, estimated using historical reporting delays.
The methodology takes into account the observed distribution pattern of the
lag. It is derived from the removal method, a well-established estimation
framework in the field of ecology.


  In epidemiological modelings, the spectral radius of the next generation
matrix evaluated at the trivial equilibrium was considered as the basic
reproduction number. Also, the global stability of the trivial equilibrium
point was determined by the left eigenvector associated to that next generation
matrix. More recently, the fraction of susceptible individuals was also
obtained from the next generation matrix. By revisiting drug sensitive and
resistant tuberculosis model, the gross reproduction number and the fraction of
susceptible individuals are calculated. Hence, the next generation matrices
shed light to the evolution of the dynamics: the beginning of the epidemics via
the reproduction number and the approaching to the epidemics level via the
asymptotic fraction of susceptible individuals.


  Decision-making about pandemic mitigation often relies upon simulation
modelling. Models of disease transmission through networks of contacts--between
individuals or between population centres--are increasingly used for these
purposes. Real-world contact networks are rich in structural features that
influence infection transmission, such as tightly-knit local communities that
are weakly connected to one another. In this paper, we propose a new flow-based
edge-betweenness centrality method for detecting bottleneck edges that connect
nodes in contact networks. In particular, we utilize convex optimization
formulations based on the idea of diffusion with p-norm network flow. Using
simulation models of COVID-19 transmission through real network data at both
individual and county levels, we demonstrate that targeting bottleneck edges
identified by the proposed method reduces the number of infected cases by up to
10% more than state-of-the-art edge-betweenness methods. Furthermore, the
proposed method is orders of magnitude faster than existing methods.


  Stochastic reaction networks are a fundamental model to describe interactions
between species where random fluctuations are relevant. The master equation
provides the evolution of the probability distribution across the discrete
state space consisting of vectors of population counts for each species.
However, since its exact solution is often elusive, several analytical
approximations have been proposed. The deterministic rate equation (DRE) gives
a macroscopic approximation as a compact system of differential equations that
estimate the average populations for each species, but it may be inaccurate in
the case of nonlinear interaction dynamics. Here we propose finite state
expansion (FSE), an analytical method mediating between the microscopic and the
macroscopic interpretations of a stochastic reaction network by coupling the
master equation dynamics of a chosen subset of the discrete state space with
the mean population dynamics of the DRE. An algorithm translates a network into
an expanded one where each discrete state is represented as a further distinct
species. This translation exactly preserves the stochastic dynamics, but the
DRE of the expanded network can be interpreted as a correction to the original
one. The effectiveness of FSE is demonstrated in models that challenge
state-of-the-art techniques due to intrinsic noise, multi-scale populations,
and multi-stability.


  Typing methods are widely used in the surveillance of infectious diseases,
outbreaks investigation and studies of the natural history of an infection. And
their use is becoming standard, in particular with the introduction of High
Throughput Sequencing (HTS). On the other hand, the data being generated is
massive and many algorithms have been proposed for phylogenetic analysis of
typing data, addressing both correctness and scalability issues. Most of the
distance-based algorithms for inferring phylogenetic trees follow the
closest-pair joining scheme. This is one of the approaches used in hierarchical
clustering. And although phylogenetic inference algorithms may seem rather
different, the main difference among them resides on how one defines cluster
proximity and on which optimization criterion is used. Both cluster proximity
and optimization criteria rely often on a model of evolution. In this work we
review, and we provide an unified view of these algorithms. This is an
important step not only to better understand such algorithms, but also to
identify possible computational bottlenecks and improvements, important to deal
with large data sets.


  The rapid spread of the Coronavirus SARS-2 is a major challenge that led
almost all governments worldwide to take drastic measures to respond to the
tragedy. Chief among those measures is the massive lockdown of entire countries
and cities, which beyond its global economic impact has created some deep
social and psychological tensions within populations. While the adopted
mitigation measures (including the lockdown) have generally proven useful,
policymakers are now facing a critical question: how and when to lift the
mitigation measures? A carefully-planned exit strategy is indeed necessary to
recover from the pandemic without risking a new outbreak. Classically, exit
strategies rely on mathematical modeling to predict the effect of public health
interventions. Such models are unfortunately known to be sensitive to some key
parameters, which are usually set based on rules-of-thumb.In this paper, we
propose to augment epidemiological forecasting with actual data-driven models
that will learn to fine-tune predictions for different contexts (e.g., per
country). We have therefore built a pandemic simulation and forecasting toolkit
that combines a deep learning estimation of the epidemiological parameters of
the disease in order to predict the cases and deaths, and a genetic algorithm
component searching for optimal trade-offs/policies between constraints and
objectives set by decision-makers. Replaying pandemic evolution in various
countries, we experimentally show that our approach yields predictions with
much lower error rates than pure epidemiological models in 75% of the cases and
achieves a 95% R2 score when the learning is transferred and tested on unseen
countries. When used for forecasting, this approach provides actionable
insights into the impact of individual measures and strategies.


  The Susceptible-Infected-Recovered (SIR) epidemic model is extensively used
for the study of the spread of infectious diseases. Even that the exact
solution of the model can be obtained in an exact parametric form, in order to
perform the comparison with the epidemiological data a simple but highly
accurate representation of the time evolution of the SIR compartments would be
very useful. In the present paper we obtain a series representation of the
solution of the SIR model by using the Laplace-Adomian Decomposition Method to
solve the basic evolution equation of the model. The solutions are expressed in
the form of infinite series. The series representations of the time evolution
of the SIR compartments are compared with the exact numerical solutions of the
model. We find that there is a good agreement between the Laplace-Adomian
semianalytical solutions containing only three terms, and the numerical
results.


  Aims: Nonpharmaceutical interventions against the spread of SARS-CoV-2 in
Germany included the cancellation of mass events (from March 8), closures of
schools and child day care facilities (from March 16) as well as a "lockdown"
(from March 23). This study attempts to assess the effectiveness of these
interventions in terms of revealing their impact on infections over time.
Methods: Dates of infections were estimated from official German case data by
incorporating the incubation period and an empirical reporting delay.
Exponential growth models for infections and reproduction numbers were
estimated and investigated with respect to change points in the time series.
Results: A significant decline of daily and cumulative infections as well as
reproduction numbers is found at March 8 (CI [7, 9]), March 10 (CI [9, 11] and
March 3 (CI [2, 4]), respectively. Further declines and stabilizations are
found in the end of March. There is also a change point in new infections at
April 19 (CI [18, 20]), but daily infections still show a negative growth. From
March 19 (CI [18, 20]), the reproduction numbers fluctuate on a level below
one. Conclusions: The decline of infections in early March 2020 can be
attributed to relatively small interventions and voluntary behavioural changes.
Additional effects of later interventions cannot be detected clearly.
Liberalizations of measures did not induce a re-increase of infections. Thus,
the effectiveness of most German interventions remains questionable. Moreover,
assessing of interventions is impeded by the estimation of true infection dates
and the influence of test volume.


  In constant parameter compartmental models an early onset of herd immunity is
at odds with estimates of R values from early stage growth. This paper utilizes
a result from the theory of interest rate modeling, namely a bond pricing
formula of Vasicek, and an approach inspired by a foundational result in
statistics, de Finetti's Theorem, to show how the modeling discrepancy can be
explained. Moreover the difference between predictions of classic constant
parameter epidemiological models and those with variation and stochastic
evolution can be reduced to simple "convexity" formulas. A novel feature of
this approach is that we do not attempt to locate a true model but only a model
that is equivalent after permutations. Convexity adjustments can also be used
for cross sectional comparisons and we derive easy to use rules of thumb for
estimating threshold infection level in one region given knowledge of threshold
infection in another.


  Breaking a complex bio-social phenomenon (epidemic) into its components,
considering the processes that determine its dynamics, formalizing the accepted
hypotheses in mathematical equations, selecting appropriate experimental and
statistical material, and constructing a mathematical model - those are typical
tasks of scientific research. A specific data processing method (balanced
identification) and appropriate information technology made it possible to
consider a number of models, determine the general biological laws of the
virus-human interaction (common to all populations), and the country specific
social features of epidemic management in the countries (or cities) under
consideration. As the initial data, only new cases are used. Data from
different countries is taken from official sources and processed in a uniform
way. The obtained estimates of the number of undetected infected are lower
estimates.


  The study of transcription remains one of the centerpieces of modern biology
with implications in settings from development to metabolism to evolution to
disease. Precision measurements using a host of different techniques including
fluorescence and sequencing readouts have raised the bar for what it means to
quantitatively understand transcriptional regulation. In particular our
understanding of the simplest genetic circuit is sufficiently refined both
experimentally and theoretically that it has become possible to carefully
discriminate between different conceptual pictures of how this regulatory
system works. This regulatory motif, originally posited by Jacob and Monod in
the 1960s, consists of a single transcriptional repressor binding to a promoter
site and inhibiting transcription. In this paper, we show how seven distinct
models of this so-called simple-repression motif, based both on equilibrium and
kinetic thinking, can be used to derive the predicted levels of gene expression
and shed light on the often surprising past success of the equilbrium models.
These different models are then invoked to confront a variety of different data
on mean, variance and full gene expression distributions, illustrating the
extent to which such models can and cannot be distinguished, and suggesting a
two-state model with a distribution of burst sizes as the most potent of the
seven for describing the simple-repression motif.


  COVID-19 is an infectious respiratory disease that the World Health
Organization has declared a pandemic. Although a global phenomenon, there is
the need to react locally with swift and informed actions, as some of the
essential plans are highly dependent on factors such as culture, geography,
laws, and customs. This paper presents our approach to mapping the geospatial
spread of COVID-19 in Mexico at the state and municipal level, the highest
allowed possible resolution in the publicly available dataset provided by the
Health Ministry. To visualize the magnitude of the infection, we offer a map
for the confirmed positive, pending, and deceases cases, while to support
mobility, we provide a geospatial visualization of $R_t$, the basic
reproduction number. This document describes the structure of the dataset, the
software tools employed, and a description of the functionality of the maps,
which definition we make publicly available.


  Clinical predictions using clinical data by computational methods are common
in bioinformatics. However, clinical predictions using information from
genomics datasets as well is not a frequently observed phenomenon in research.
Precision medicine research requires information from all available datasets to
provide intelligent clinical solutions. In this paper, we have attempted to
create a prediction model which uses information from both clinical and
genomics datasets. We have demonstrated multiclass disease predictions based on
combined clinical and genomics datasets using machine learning methods. We have
created an integrated dataset, using a clinical (ClinVar) and a genomics (gene
expression) dataset, and trained it using instance-based learner to predict
clinical diseases. We have used an innovative but simple way for multiclass
classification, where the number of output classes is as high as 75. We have
used Principal Component Analysis for feature selection. The classifier
predicted diseases with 73\% accuracy on the integrated dataset. The results
were consistent and competent when compared with other classification models.
The results show that genomics information can be reliably included in datasets
for clinical predictions and it can prove to be valuable in clinical
diagnostics and precision medicine.


  High impact epidemics constitute one of the largest threats humanity is
facing in the 21st century. Testing, contact tracing and quarantining are
critical in slowing down epidemic dynamics, but may prove insufficient for
highly contagious diseases. In the absence of pharmaceutical interventions,
physical distancing measures remain as the last resort to avoid a widespread
outbreak. Here we show that such combined countermeasures drastically change
the rules of the epidemic transition if testing capacities are limited: Instead
of continuous the response to countermeasures becomes discontinuous and rather
than following the conventional exponential growth, the outbreak accelerates
and scales super-exponentially during an intermediate period. As a consequence,
containment measures either suffice to stop the outbreak at low total case
numbers or fail catastrophically if marginally too weak, thus implying large
uncertainties in reliably estimating overall epidemic dynamics, both during
initial phases and during second wave scenarios.


  Genomics has become an essential technology for surveilling emerging
infectious disease outbreaks. A wide range of technologies and strategies for
pathogen genome enrichment and sequencing are being used by laboratories
worldwide, together with different, and sometimes ad hoc, analytical procedures
for generating genome sequences. As a result, public repositories now contain
non-standard entries of varying quality. A standardized analytical process for
consensus genome sequence determination, particularly for outbreaks such as the
ongoing COVID-19 pandemic, is critical to provide a solid genomic basis for
epidemiological analyses and well-informed decision making. To address this
need, we have developed a bioinformatic workflow to standardize the analysis of
SARS-CoV-2 sequencing data generated with either the Illumina or Oxford
Nanopore platforms. Using an intuitive web-based interface, this workflow
automates SARS-CoV-2 reference-based genome assembly, variant calling, lineage
determination, and provides the ability to submit the consensus sequence and
necessary metadata to GenBank or GISAID. Given a raw Illumina or Oxford
Nanopore FASTQ read file, this web-based platform enables non-bioinformatics
experts to automatically produce a SARS-CoV-2 genome that is ready for
submission to GISAID or GenBank.
  Availability:https://edge-covid19.edgebioinformatics.org;https://github.com/LANL-Bioinformatics/EDGE/tree/SARS-CoV2


  Inherent stochasticity in gene expression leads to distributions of mRNA copy
numbers in a population of identical cells. These distributions are determined
primarily by the multitude of states of a gene promoter, each driving
transcription at a different rate. In an era where single-cell mRNA copy number
data are more and more available, there is an increasing need for fast
computations of mRNA distributions. In this paper, we present a method for
computing separate distributions for each species of mRNA molecules, i. e.
mRNAs that have been either partially or fully processed post-transcription.
The method involves the integration over all possible realizations of promoter
states, which we cast into a set of linear ordinary differential equations of
dimension $M\times n_j$, where $M$ is the number of available promoter states
and $n_j$ is the mRNA copy number of species $j$ up to which one wishes to
compute the probability distribution. This approach is superior to solving the
Master equation (ME) directly in two ways: a) the number of coupled
differential equations in the ME approach is
$M\times\Lambda_1\times\Lambda_2\times ...\times\Lambda_L$, where $\Lambda_j$
is the cutoff for the probability of the $j^{\text{th}}$ species of mRNA; and
b) the ME must be solved up to the cutoffs $\Lambda_j$, which are {\it ad hoc}
and must be selected {\it a priori}. In our approach, the equation for the
probability to observe $n$ mRNAs of any species depends only on the the
probability of observing $n-1$ mRNAs of that species, thus yielding a correct
probability distribution up to an arbitrary $n$. To demonstrate the validity of
our derivations, we compare our results with Gillespie simulations for ten
randomly selected system parameters.


  Under high complexity - given by pervasive interdependence between
constituent elements of a decision in an NK landscape - our algorithm obtains
fitness superior to that reported in extant research. We distribute the
decision elements comprising a decision into clusters. When a change in value
of a decision element is considered, a forward move is made if the aggregate
fitness of the cluster members residing alongside the decision element is
higher. The decision configuration with the highest fitness in the path is
selected. Increasing the number of clusters obtains even higher fitness.
Further, implementing moves comprising of up to two changes in a cluster also
obtains higher fitness. Our algorithm obtains superior outcomes by enabling
more extensive search, allowing inspection of more distant configurations. We
name this algorithm the muddling through algorithm, in memory of Charles
Lindblom who spotted the efficacy of the process long before sophisticated
computer simulations came into being.


  In this paper, we propose a continuous-time stochastic intensity model,
namely, two-phase dynamic contagion process(2P-DCP), for modelling the epidemic
contagion of COVID-19 and investigating the lockdown effect based on the
dynamic contagion model introduced by Dassios and Zhao (2011). It allows
randomness to the infectivity of individuals rather than a constant
reproduction number as assumed by standard models. Key epidemiological
quantities, such as the distribution of final epidemic size and expected
epidemic duration, are derived and estimated based on real data for various
regions and countries. The associated time lag of the effect of intervention in
each country or region is estimated. Our results are consistent with the
incubation time of COVID-19 found by recent medical study. We demonstrate that
our model could potentially be a valuable tool in the modeling of COVID-19.
More importantly, the proposed model of 2P-DCP could also be used as an
important tool in epidemiological modelling as this type of contagion models
with very simple structures is adequate to describe the evolution of regional
epidemic and worldwide pandemic.


  Background. The paper concerns the SARS-CoV2 (COVID-19) pandemic that,
starting from the end of February 2020, began spreading along the Italian
peninsula, by first attacking small communities in north regions, and then
extending to the center and south of Italy, including the two main islands.
  Objective. The creation of a forecast model that manages to alert the
decision-making bodies and, in particular, the healthcare system, to hinder the
emergence of any other pandemic outbreaks, or the arrival of subsequent
pandemic waves.
  Methods. A new mathematical model to describe the pandemic is given. The
model includes the class of undiagnosed infected people, and has a multi-region
extension, to cope with the in-time and in-space heterogeneity of the epidemic.
  Results. We obtain a robust and reliable tool for the forecast of the total
and active cases, which can be also used to simulate different scenarios.
  Conclusions. We are able to address a number of issues, such as assessing the
adoption of the lockdown in Italy, started from 11 March 2020, and how to
employ a rapid screening test campaign for containing the epidemic.


  We propose a dynamical model for describing the spread of epidemics. This
model is an extension of the SIQR (susceptible-infected-quarantined-recovered)
and SIRP (susceptible-infected-recovered-pathogen) models used earlier to
describe various scenarios of epidemic spreading. As compared to the basic SIR
model, our model takes into account two possible routes of contagion
transmission: direct from the infected compartment to the susceptible
compartment and indirect via some intermediate medium or fomites. Transmission
rates are estimated in terms of average distances between the individuals in
selected social environments and characteristic time spans for which the
individuals stay in each of these environments. We also introduce a collective
economic resource associated with the average amount of money or income per
individual to describe the socioeconomic interplay between the spreading
process and the resource available to infected individuals. The
epidemic-resource coupling is supposed to be of activation type, with the
recovery rate governed by the Arrhenius-like law. Our model brings an advantage
of building various control strategies to mitigate the effect of epidemic and
can be applied, in particular, to modeling the spread of COVID-19.


  While many epidemiological models have being proposed to understand and
handle COVID-19, too little has been invested to understand how the virus
replicates in the human body and potential antiviral can be used to control the
replication cycle. In this work, using a control theoretical approach,
validated mathematical models of SARS-CoV-2 in humans are properly
characterized. A complete analysis of the main dynamic characteristic is
developed based on the reproduction number. The equilibrium regions of the
system are fully characterized, and the stability of such a regions, formally
established. Mathematical analysis highlights critical conditions to decrease
monotonically SARS-CoV-2 in the host, such conditions are relevant to tailor
future antiviral treatments. Simulation results show the potential benefits of
the aforementioned system characterization.


  We quantified the probability of developing symptoms (respiratory or fever
\geq 37.5 {\deg}C) and critical disease (requiring intensive care or resulting
in death) of SARS-CoV-2 positive subjects. 5,484 contacts of SARS-CoV-2 index
cases detected in Lombardy, Italy were analyzed, and positive subjects were
ascertained via nasal swabs and serological assays. 73.9% of all infected
individuals aged less than 60 years did not develop symptoms (95% confidence
interval: 71.8-75.9%). The risk of symptoms increased with age. 6.6% of
infected subjects older than 60 years had critical disease, with males at
significantly higher risk.


  Testing, tracking and tracing abilities have been identified as pivotal in
helping countries to safely reopen activities after the first wave of the
COVID-19 virus. Contact tracing apps give the unprecedented possibility to
reconstruct graphs of daily contacts, so the question is who should be tested?
As human contact networks are known to exhibit community structure, in this
paper we show that the Kemeny constant of a graph can be used to identify and
analyze bridges between communities in a graph. Our "Kemeny indicator" is the
change in Kemeny constant when a node or edge is removed from the graph. We
show that testing individuals who are associated with large values of the
Kemeny indicator can help in efficiently intercepting new virus outbreaks, when
they are still in their early stage. Extensive simulations provide promising
results in early identification and in blocking possible "super-spreaders"
links that transmit disease between different communities.


  In this paper, we aim to synthesize cell microscopy images under different
molecular interventions, motivated by practical applications to drug
development. Building on the recent success of graph neural networks for
learning molecular embeddings and flow-based models for image generation, we
propose Mol2Image: a flow-based generative model for molecule to cell image
synthesis. To generate cell features at different resolutions and scale to
high-resolution images, we develop a novel multi-scale flow architecture based
on a Haar wavelet image pyramid. To maximize the mutual information between the
generated images and the molecular interventions, we devise a training strategy
based on contrastive learning. To evaluate our model, we propose a new set of
metrics for biological image generation that are robust, interpretable, and
relevant to practitioners. We show quantitatively that our method learns a
meaningful embedding of the molecular intervention, which is translated into an
image representation reflecting the biological effects of the intervention.


  According to the standard protocol of spatial public goods game, a cooperator
player invests not only into his own game but also into the games organized by
neighboring partners. In this work, we relax this assumption by allowing
cooperators to decide which neighboring group to prefer instead of supporting
them uniformly. In particular, we assume that they select their most successful
neighbor and focus external investments exclusively into the related group. We
show that this very simple alteration of the dynamical rule results in a
surprisingly positive evolutionary outcome -- cooperators prevail even in harsh
environment represented by small values of the synergy factor in the game. The
microscopic mechanism behind the reported success of the cooperator strategy
can be explained by a blocking mechanism which affects the propagations of
competing strategies in a biased way. Our results, which remain intact by using
different interaction topologies, reveal that it could be beneficial to
concentrate individual efforts to reach a higher global wellbeing.


  Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a novel virus
of beta coronavirus genus originated in Wuhan, China in December 2019, resulted
in the pandemic spread of Corona Virus Disease 2019 (COVID-19) worldwide. This
genus also contains SARS-CoV (originated in China 2002-2003) and MERS-CoV
(found in Saudi Arabia 2012). The nucleotide sequences of SARS-CoV-2 is closer
to SARS-CoV (with approximately 80% identity) than to MERS-CoV. Despite these
similarities, SARS-CoV-2 has two main features over the other coronaviruses.
The first is the high contagious rate and the second is the immune response
evasion. The higher transmission ability makes this virus quickly spread
worldwide with a high mortality rate and more economic losses. This review will
provide an overview of the current knowledge on the role of some viral and host
cell factors in higher transmission and contagious spread of SARS-CoV-2.


  Toric differential inclusions occur as key dynamical systems in the context
of the Global Attractor Conjecture. We introduce the notions of minimal
invariant regions and minimal globally attracting regions for toric
differential inclusions. We describe a procedure for constructing explicitly
the minimal invariant and minimal globally attracting regions for
two-dimensional toric differential inclusions. In particular, we obtain
invariant regions and globally attracting regions for two-dimensional weakly
reversible or endotactic dynamical systems (even if they have time-dependent
parameters).


  We discuss several issues of statistical design, data collection, analysis,
communication, and decision making that have arisen in recent and ongoing
coronavirus studies, focusing on tools for assessment and propagation of
uncertainty. This paper does not purport to be a comprehensive survey of the
research literature; rather, we use examples to illustrate statistical points
that we think are important.


  On 19th March, the World Health Organisation declared a pandemic. Through
this global spread, many nations have witnessed exponential growth of confirmed
cases brought under control by severe mass quarantine or lockdown measures.
However, some have, through a different timeline of actions, prevented this
exponential growth. Currently as some continue to tackle growth, others attempt
to safely lift restrictions whilst avoiding a resurgence. This study seeks to
quantify the impact of government actions in mitigating viral transmission of
SARS-CoV-2 by a novel soft computing approach that makes concurrent use of a
neural network model, to predict the daily slope increase of cumulative
infected, and an optimiser, with a parametrisation of the government
restriction time series, to understanding the best set of mitigating actions.
Data for two territories, Italy and Taiwan, have been gathered to model
government restrictions in traveling, testing and enforcement of social
distance measures as well as people connectivity and adherence to government
actions. It is found that a larger and earlier testing campaign with tighter
entry restrictions benefit both regions, resulting in significantly less
confirmed cases. Interestingly, this scenario couples with an earlier but
milder implementation of nationwide restrictions for Italy, thus supporting
Taiwan's lack of nationwide lockdown. The results, found with a purely
data-driven approach, are in line with the main findings of mathematical
epidemiological models, proving that the proposed approach has value and that
the data alone contains valuable knowledge to inform decision makers.


  Public cloud computing environments, such as Amazon AWS, Microsoft Azure, and
the Google Cloud Platform, have achieved remarkable improvements in
computational performance in recent years, and are also expected to be able to
perform massively parallel computing. As the cloud enables users to use
thousands of CPU cores and GPU accelerators casually, and various software
types can be used very easily by cloud images, the cloud is beginning to be
used in the field of bioinformatics. In this study, we ported the original
protein-protein interaction prediction (protein-protein docking) software,
MEGADOCK, into Microsoft Azure as an example of an HPC cloud environment. A
cloud parallel computing environment with up to 1,600 CPU cores and 960 GPUs
was constructed using four CPU instance types and two GPU instance types, and
the parallel computing performance was evaluated. Our MEGADOCK on Azure system
showed a strong scaling value of 0.93 for the CPU instance when H16 instance
with 100 instances were used compared to 50, and a strong scaling value of 0.89
for the GPU instance when NC24 instance with 20 were used compared to 5.
Moreover, the results of the usage fee and total computation time supported
that using a GPU instance reduced the computation time of MEGADOCK and the
cloud usage fee required for the computation. The developed environment
deployed on the cloud is highly portable, making it suitable for applications
in which an on-demand and large-scale HPC environment is desirable.


  For a pair consisting of a gene tree and a species tree, the ancestral
configurations at an internal node of the species tree are the distinct sets of
gene lineages that can be present at that node. Ancestral configurations appear
in computations of gene tree probabilities under evolutionary models
conditional on fixed species trees, and the enumeration of root ancestral
configurations -- ancestral configurations at the root of the species tree --
assists in describing the complexity of these computations. In the case that
the gene tree matches the species tree in topology, we study the distribution
of the number of root ancestral configurations of a random labeled tree
topology under each of two models.


  A new coronavirus identified as SARS-CoV-2 virus has brought the world to a
state of crisis, causing a major pandemic, claiming more than 433,000 lives and
instigating major financial damage to the global economy. Despite current
efforts, developing safe and effective treatments remains a major challenge.
Moreover, new strains of the virus are likely to emerge in the future. To
prevent future pandemics, several drugs with various mechanisms of action are
required. Drug discovery efforts against the virus fall into two main
categories: (a) monoclonal antibodies targeting the spike protein of the virus
and blocking it from entry; (b) small molecule inhibitors targeting key
proteins of the virus, interfering with replication and translation of the
virus. In this study, we are presenting a computational investigation of a
potential drug candidate that targets SARS-CoV-2 protease, a viral protein
critical for replication and translation of the virus.


  The first mitigation response to the Covid-19 pandemic was to limit
person-to-person interaction as much as possible. This was implemented by the
temporary closing of many workplaces and people were required to follow social
distancing. Networks are a great way to represent interactions among people and
the temporary severing of these interactions. Here, we present a network model
of human-human interactions that could be mediators of disease spread. The
nodes of this network are individuals and different types of edges denote
family cliques, workplace interactions, interactions arising from essential
needs, and social interactions. Each individual can be in one of four states:
susceptible, infected, immune, and dead. The network and the disease parameters
are informed by the existing literature on Covid-19. Using this model, we
simulate the spread of an infectious disease in the presence of various
mitigation scenarios. For example, lockdown is implemented by deleting edges
that denote non-essential interactions. We validate the simulation results with
the real data by matching the basic and effective reproduction numbers during
different phases of the spread. We also simulate different possibilities of the
slow lifting of the lockdown by varying the transmission rate as facilities are
slowly opened but people follow prevention measures like wearing masks etc. We
make predictions on the probability and intensity of a second wave of infection
in each of these scenarios.


  We consider a model due to Piero Poletti and collaborators that adds
spontaneous human behavioral change to the standard SIR epidemic model. In its
simplest form, the Poletti model adds one differential equation, motivated by
evolutionary game theory, to the SIR model. The new equation describes the
evolution of a variable $x$ that represents the fraction of the population
using normal behavior. The remaining fraction $1-x$ uses altered behavior such
as staying home, social isolation, mask wearing, etc. Normal behavior offers a
higher payoff when the number of infectives is low; altered behavior offers a
higher payoff when the number is high. We show that the entry-exit function of
geometric singular perturbation theory can be used to analyze the model in the
limit in which behavior changes on a much faster time scale than that of the
epidemic. In particular, behavior does not change as soon as a different
behavior has a higher payoff; current behavior is sticky. The delay until
behavior changes in predicted by the entry-exit function.


  Predicting the structure of multi-protein complexes is a grand challenge in
biochemistry, with major implications for basic science and drug discovery.
Computational structure prediction methods generally leverage pre-defined
structural features to distinguish accurate structural models from less
accurate ones. This raises the question of whether it is possible to learn
characteristics of accurate models directly from atomic coordinates of protein
complexes, with no prior assumptions. Here we introduce a machine learning
method that learns directly from the 3D positions of all atoms to identify
accurate models of protein complexes, without using any pre-computed
physics-inspired or statistical terms. Our neural network architecture combines
multiple ingredients that together enable end-to-end learning from molecular
structures containing tens of thousands of atoms: a point-based representation
of atoms, equivariance with respect to rotation and translation, local
convolutions, and hierarchical subsampling operations. When used in combination
with previously developed scoring functions, our network substantially improves
the identification of accurate structural models among a large set of possible
models. Our network can also be used to predict the accuracy of a given
structural model in absolute terms. The architecture we present is readily
applicable to other tasks involving learning on 3D structures of large atomic
systems.


  The synthesis of ATP, life's 'universal energy currency', is the most
prevalent chemical reaction in biological systems, and is responsible for
fueling nearly all cellular processes, from nerve impulse propagation to DNA
synthesis. ATP synthases, the family of enzymes that carry out this endless
task, are nearly as ubiquitous as the energy-laden molecule they are
responsible for making. The F-type ATP synthase (F-ATPase) is found in every
domain of life, and is believed to predate the divergence of these lineages
over 1.5 billion years ago. These enzymes have therefore facilitated the
survival of organisms in a wide range of habitats, ranging from the deep-sea
thermal vents to the human intestine. In this review, we present an overview of
the current knowledge of the structure and function of F-type ATPases,
highlighting several adaptations that have been characterized across taxa. We
emphasize the importance of studying these features within the context of the
enzyme's particular lipid environment: Just as the interactions between an
organism and its physical environment shape its evolutionary trajectory,
ATPases are impacted by the membranes within which they reside. We argue that a
comprehensive understanding of the structure, function, and evolution of
membrane proteins -- including ATP synthase -- requires such an integrative
approach.


  This technical report addresses a pressing issue in the trajectory of the
coronavirus outbreak; namely, the rate at which effective immunity is lost
following the first wave of the pandemic. This is a crucial epidemiological
parameter that speaks to both the consequences of relaxing lockdown and the
propensity for a second wave of infections. Using a dynamic causal model of
reported cases and deaths from multiple countries, we evaluated the evidence
models of progressively longer periods of immunity. The results speak to an
effective population immunity of about three months that, under the model,
defers any second wave for approximately six months in most countries. This may
have implications for the window of opportunity for tracking and tracing, as
well as for developing vaccination programmes, and other therapeutic
interventions.


  This paper is concerned with nonlinear modeling and analysis of the COVID-19
pandemic currently ravaging the planet. There are two objectives: to arrive at
an appropriate model that captures the collected data faithfully, and to use
that as a basis to explore the nonlinear behavior. We use a nonlinear SEIR
(Susceptible, Exposed, Infectious & Removed) transmission model with added
behavioral and government policy dynamics. We develop a genetic algorithm
technique to identify key model parameters employing COVID19 data from South
Korea. Stability, bifurcations and dynamic behavior are analyzed. Parametric
analysis reveals conditions for sustained epidemic equilibria to occur. This
work points to the value of nonlinear dynamic analysis in pandemic modeling and
demonstrates the dramatic influence of social and government behavior on
disease dynamics.


  Deep learning methods have permeated into the research area of computer-aided
drug design. The deep learning generative model and classical algorithm can be
simultaneously used for three-dimensional (3D) drug design in the 3D pocket of
the receptor. Here, three aspects of MolAICal are illustrated for drug design:
in the first part, the MolAICal uses the genetic algorithm, Vinardo score and
deep learning generative model trained by generative adversarial net (GAN) for
drug design. In the second part, the deep learning generative model is trained
by drug-like molecules from the drug database such as ZINC database. The
MolAICal invokes the deep learning generative model and molecular docking for
drug virtual screening automatically. In the third part, the useful drug tools
are added for calculating the relative properties such as Pan-assay
interference compounds (PAINS), Lipinski's rule of five, synthetic
accessibility (SA), and so on. Besides, the structural similarity search and
quantitative structure-activity relationship (QSAR), etc are also embedded for
the calculations of drug properties in the MolAICal. MolAICal will constantly
optimize and develop the current and new modules for drug design. The MolAICal
can help the scientists, pharmacists and biologists to design the rational 3D
drugs in the receptor pocket through the deep learning model and classical
programming. MolAICal is free of charge for any academic and educational
purposes, and it can be downloaded from the website https://molaical.github.io.


  Evolutionary games provide the theoretical backbone for many aspects of our
social life: from cooperation to crime, from climate inaction to imperfect
vaccination and epidemic spreading, from antibiotics overuse to biodiversity
preservation. An important, and so far overlooked, aspect of reality is the
diverse strategic identities of individuals. While applying the same strategy
to all interaction partners may be an acceptable assumption for simpler forms
of life, this fails to account} for the behavior of more complex living beings.
For instance, we humans act differently around different people. Here we show
that allowing individuals to adopt different strategies with different partners
yields a very rich evolutionary dynamics, including time-dependent coexistence
of cooperation and defection, system-wide shifts in the dominant strategy, and
maturation in individual choices. Our results are robust to variations in
network type and size, and strategy updating rules. Accounting for diverse
strategic identities thus has far-reaching implications in the mathematical
modeling of social games.


  An analytic solution is obtained to the SEIR Epidemic Model. The solution is
created by constructing a single second-order nonlinear differential equation
in $\ln S$ and analytically continuing its divergent power series solution such
that it matches the correct long-time exponential damping of the epidemic
model. This is achieved through an asymptotic approximant (Barlow et. al, 2017,
Q. Jl Mech. Appl. Math, 70 (1), 21-48) in the form of a modified symmetric
Pad\'e approximant that incorporates this damping. The utility of the
analytical form is demonstrated through its application to the COVID-19
pandemic.


  We show that the dynamics of the number of deaths due to Covid in different
countries is to a large extent universal once the origin of time is chosen to
be the start of the lockdown, and the number of death is rescaled by the total
number of deaths after the lockdown, itself a proxy of the number of infections
at the start of the lockdown. Such a curve collapse is much less convincing
when normalizing by the total population. Sweden, with its no-lockdown,
light-touch approach, is the only outlier that deviates considerably from the
average behavior. We argue that these model-free findings provide strong
support for the effectiveness of the lockdowns in mitigating the lethality of
the virus.


  SKCM is the most dangerous one of skin cancer, its high degree of malignant,
is the leading cause of skin cancer. And the level of radiation treatment and
chemical treatment is minimal, so the mortality is high. Because of its complex
molecular and cellular heterogeneity, the existing prediction model of skin
cancer risk is not ideal. In this study, we developed an immune-related lncRNAs
model to predict the prognosis of patients with SKCM. Screening for
SKCM-related differential expression of lncRNA from TCGA. Identified
immune-related lncRNAs and lncRNA-related mRNA based on the co-expression
method. Through univariate and multivariate analysis, an immune-related lncRNA
model is established to analyze the prognosis of SKCM patients. A 4-lncRNA skin
cancer prediction model was constructed, including MIR155HG, AL137003.2,
AC011374.2, and AC009495.2. According to the model, SKCM samples were divided
into a high-risk group and low-risk group, and predict the survival of the two
groups in 30 years. The area under the ROC curve is 0.749, which shows that the
model has excellent performance. We constructed a 4-lncRNA model to predict the
prognosis of patients with SKCM, indicating that these lncRNAs may play a
unique role in the carcinogenesis of SKCM.


  Spatial models for occupancy data are used to estimate and map the true
presence of a species, which may depend on biotic and abiotic factors as well
as spatial autocorrelation. Traditionally researchers have accounted for
spatial autocorrelation in occupancy data by using a correlated normally
distributed site-level random effect, which might be incapable of identifying
nontraditional spatial dependence such as discontinuities and abrupt
transitions. Machine learning approaches have the potential to identify and
model nontraditional spatial dependence, but these approaches do not account
for observer errors such as false absences. By combining the flexibility of
Bayesian hierarchal modeling and machine learning approaches, we present a
general framework to model occupancy data that accounts for both traditional
and nontraditional spatial dependence as well as false absences. We demonstrate
our framework using six synthetic occupancy data sets and two real data sets.
Our results demonstrate how to identify and model both traditional and
nontraditional spatial dependence in occupancy data which enables a broader
class of spatial occupancy models that can be used to improve predictive
accuracy and model adequacy.


  The Trojan Y Chromosome strategy (TYC) is a genetic biocontrol strategy
designed to alter the sex ratio of a target invasive population by reducing the
number of females over time. Recently an alternative strategy is introduced,
that mimics the TYC strategy by harvesting females whilst stocking males . We
consider the FHMS strategy, with a weak Allee effect, and show that the
extinction boundary need note be hyperbolic. To the best of our knowledge, this
is the first example of a non-hyperbolic extinction boundary in mating models,
structured by sex. Next, we consider the spatially explicit model and show that
the weak Allee effect is both sufficient and necessary for Turing patterns to
occur. We discuss the applicability of our results to large scale biocontrol,
as well as compare and contrast our results to the case with a strong Allee
effect.


  The growing literature on the propagation of COVID-19 relies on various
dynamic SIR-type models (Susceptible-Infected-Recovered) which yield
model-dependent results. For transparency and ease of comparing the results, we
introduce a common representation of the SIR-type stochastic epidemiological
models. This representation is a discrete time transition model, which allows
us to classify the epidemiological models with respect to the number of states
(compartments) and their interpretation. Additionally, the transition model
eliminates several limitations of the deterministic continuous time
epidemiological models which are pointed out in the paper. We also show that
all SIR-type models have a nonlinear (pseudo) state space representation and
are easily estimable from an extended Kalman filter.


  Recent literature has suggested that climate conditions have considerably
significant influences on the transmission of coronavirus COVID-19. However,
there is a lack of comprehensive study that investigates the relationships
between multiple weather factors and the development of COVID-19 pandemic while
excluding the impact of social factors. In this paper, we study the
relationships between six main weather factors and the infection statistics of
COVID-19 on 250 cities in Mainland China. Our correlation analysis using
weather and infection statistics indicates that all the studied weather factors
are correlated with the spread of COVID-19, where precipitation shows the
strongest correlation. We also build a weather-aware predictive model that
forecasts the number of infected cases should there be a second wave of the
outbreak in Mainland China. Our predicted results show that cities located in
different geographical areas are likely to be challenged with the second wave
of COVID-19 at very different time periods and the severity of the outbreak
varies to a large degree, in correspondence with the varying weather
conditions.