  Background. Pneumococcal vaccination (PCV13) for children (as well as older
adults) in Russia was introduced in 2014, with no prior PCV7 use. While
pneumonia hospitalization rates, both in children and adults did not decrease
in the following years, respiratory mortality rates did decline with time.
Moreover, there is a strong association between antibiotic, including multidrug
resistance and PCV13 vaccine serotypes for S. pneumoniae in children in Russia,
and presence of S. pneumoniae among sepsis cases in Russia during the recent
years has been very low. Annual variability in influenza circulation may affect
some of the changes in respiratory mortality rates. Methods. We applied the
inference method from our recent study of influenza-associated mortality in
Russia to relate monthly respiratory mortality rates between 09/2010 and
08/2019 to monthly indices of influenza circulation and baseline rates of
mortality not associated with influenza, allowing for a potential trend in the
baseline rates starting 2015. Results. Baseline rates of respiratory mortality
not associated with influenza decreased with time starting from the 2014/15
season (with each season running from September to August), with the 2018/19
baseline rates of respiratory mortality per 100,000 individuals being lower by
12.41 (95% CI (10.6,14.2)) compared to the 2010-2014 period. Conclusions. While
improvement in care might have contributed to the reduction in the rates of
respiratory mortality not associated with influenza in Russia, the above
temporally consistent reduction is compatible with the gradual replacement of
vaccine serotypes in the transmission and carriage of S. pneumoniae following
the introduction of PCV13. Further work is needed to better understand the
impact of PCV13 on the epidemiology of respiratory infections and related
mortality in Russia.


  Adoptive T cell based immunotherapy is gaining significant traction in cancer
treatment. Despite its limited success, so far, in treating solid cancers, it
is increasingly successful, demonstrating to have a broader therapeutic
potential. In this paper we develop a mathematical model to study the efficacy
of engineered T cell receptor (TCR) T cell therapy targeting the E7 antigen in
cervical cancer cell lines. We consider a dynamical system that follows the
population of cancer cells, TCR T cells, and IL-2. We demonstrate that there
exists a TCR T cell dosage window for a successful cancer elimination that can
be expressed in terms of the initial tumor size. We obtain the TCR T cell dose
for two cervical cancer cell lines: 4050 and CaSki. Finally, a combination
therapy of TCR T cell and IL-2 treatment is studied. We show that certain
treatment protocols can improve therapy responses in the 4050 cell line, but
not in the CaSki cell line.


  Most of the Middle Palaeolithic evidence of Central Italy still lacks a
reliable chrono-cultural framework mainly due to research history. In this
context Grotta dei Santi, a wide cave located on Monte Argentario, on the
southern coast of Tuscany, is particularly relevant as it contains a very well
preserved sequence including several Mousterian layers.


  Analysis of the structure of biological networks often uses statistical tests
to establish the over-representation of motifs, which are thought to be
important building blocks of such networks, related to their biological
functions. However, there is disagreement as to the statistical significance of
these motifs, and there are potential problems with standard methods for
estimating this significance. Exponential random graph models (ERGMs) are a
class of statistical model that can overcome some of the shortcomings of
commonly used methods for testing the statistical significance of motifs. ERGMs
were first introduced into the bioinformatics literature over ten years ago but
have had limited application to biological networks, possibly due to the
practical difficulty of estimating model parameters. Advances in estimation
algorithms now afford analysis of much larger networks in practical time. We
illustrate the application of ERGM to both an undirected protein-protein
interaction (PPI) network and directed gene regulatory networks. ERGM models
indicate over-representation of triangles in the PPI network, and confirm
results from previous research as to over-representation of transitive
triangles (feed-forward loop) in an E. coli and a yeast regulatory network. We
also confirm, using ERGMs, previous research showing that under-representation
of the cyclic triangle (feedback loop) can be explained as a consequence of
other topological features.


  We study the interaction between epidemic spreading and a vaccination
process. We assume that, similar to the disease spreading, also the vaccination
process occurs through direct contact, i.e., it follows the standard
susceptible-infected-susceptible (SIS) dynamics. The two competing processes
are asymmetrically coupled as vaccinated nodes can directly become infected at
a reduced rate with respect to susceptible ones. We study analytically the
model in the framework of mean-field theory finding a rich phase-diagram. When
vaccination provides little protection toward infection, two continuous
transitions separate a disease-free immunized state from vaccinated-free
epidemic state, with an intermediate mixed state where susceptible, infected
and vaccinated individuals coexist. As vaccine efficiency increases, a
tricritical point leads to a bistable regime and discontinuous phase
transitions emerge. Numerical simulations for homogeneous random networks agree
very well with analytical predictions.


  Having control over species abundances and community resilience is of great
interest for experimental, agricultural, industrial and conservation purposes.
Here, we theoretically explore the possibility of manipulating ecological
communities by modifying pairwise interactions. Specifically, we establish
which interaction values should be modified, and by how much, in order to alter
the composition or resilience of a community towards a favorable direction.
While doing so, we also take into account the experimental difficulties in
making such modifications by including in our optimization process, a cost
parameter, which penalizes large modifications. In addition to prescribing what
changes should be made to interspecies interactions given some modification
cost, our approach also serves to establish the limits of community control,
i.e. how well can one approach an ecological goal at best, even when not
constrained by cost.


  Heterogeneity is a hallmark of all cancers. Tumor heterogeneity is found at
different levels -- interpatient, intrapatient, and intratumor heterogeneity.
All of them pose challenges for clinical treatments. The latter two scenarios
can also increase the risk of developing drug resistance. Although the
existence of tumor heterogeneity has been known for two centuries, a clear
understanding of its origin is still elusive, especially at the level of
intratumor heterogeneity (ITH). The coexistence of different subpopulations
within a single tumor has been shown to play crucial roles during all stages of
carcinogenesis. Here, using concepts from evolutionary game theory and public
goods game, often invoked in the context of the tragedy of commons, we explore
how the interactions among subclone populations influence the establishment of
ITH. By using an evolutionary model, which unifies several experimental results
in distinct cancer types, we develop quantitative theoretical models for
explaining data from {\it in vitro} experiments involving pancreatic cancer as
well as {\it vivo} data in glioblastoma multiforme. Such physical and
mathematical models complement experimental studies, and could optimistically
also provide new ideas for the design of efficacious therapies for cancer
patients.


  MicroRNAs are small non-coding nucleotide sequences that regulate target
protein expression at post-transcriptional levels. Biogenesis of microRNA is a
highly regulated multi-step pathway. Regulation of miRNA biogenesis can be
caused directly by the components of the biogenesis pathway or indirectly by
other regulators. In this study, we have built a detailed mathematical model of
microRNA biogenesis to investigate the regulatory role of biogenesis pathway
components. We extended a previous model to incorporate Microprocessor
regulation of DGCR8 synthesis, exportin-mediated transport to the cytoplasm,
and positive auto-regulation catalysed by mature miRNA translocation into the
nucleus. Our simulation results lead to three hypotheses (i) Biogenesis is
robust to Dicer protein levels at higher Exportin protein levels; (ii) Higher
miRNA transcript formation may lead to lower RISC levels: an optimal level of
both precursor miRNA and Dicer is required for optimal miRNA formation at lower
levels of Exportin protein; and (iii) The positive auto-regulation by mature
miRNA translocation into the nucleus can decrease the net functional
cytoplasmic miRNA. Wherever possible, we compare these results to experimental
observations not used in the model construction or calibration.


  Molecular processes of neuronal learning have been well-described. However,
learning mechanisms of non-neuronal cells have not been fully understood at the
molecular level. Here, we discuss molecular mechanisms of cellular learning,
including conformational memory of intrinsically disordered proteins and
prions, signaling cascades, protein translocation, RNAs (microRNA and lncRNA),
and chromatin memory. We hypothesize that these processes constitute the
learning of signaling networks and correspond to a generalized Hebbian learning
process of single, non-neuronal cells, and discuss how cellular learning may
open novel directions in drug design and inspire new artificial intelligence
methods.


  Many real networks feature the property of nestedness, i.e. the neighbours of
nodes with a few connections are hierarchically nested within the neighbours of
nodes with more connections. Despite the abstract simplicity of this notion,
different mathematical definitions of nestedness have been proposed, sometimes
giving contrasting results. Moreover, there is an ongoing debate on the
statistical significance of nestedness, since even random networks where the
number of connections (degree) of each node is fixed to its empirical value are
typically as nested as real-world ones. Here we propose a clarification that
exploits the recent finding that random networks where the degrees are enforced
as hard constraints (microcanonical ensembles) are thermodynamically different
from random networks where the degrees are enforced as soft constraints
(canonical ensembles). We show that if the real network is perfectly nested,
then the two ensembles are trivially equivalent and the observed nestedness,
independently of its definition, is indeed an unavoidable consequence of the
empirical degrees. On the other hand, if the real network is not perfectly
nested, then the two ensembles are not equivalent and alternative definitions
of nestedness can be even positively correlated in the canonical ensemble and
negatively correlated in the microcanonical one. This result disentangles
distinct notions of nestedness captured by different metrics and highlights the
importance of making a principled choice between hard and soft constraints in
null models of ecological networks.


  We introduce and study a learning theory which is roughly automatic, that is,
it does not require but a minimum of initial programming, and is based on the
potential computational phenomenon of self-reference, (i.e. the potential
ability of an algorithm to have its program as an input).
  The conclusions agree with scientific findings in both biology and
neuroscience and provide a plethora of explanations both (in conjunction with
Darwinism) about evolution, as well as for the functionality and learning
capabilities of human brain, (most importantly), as we perceive them in
ourselves.


  Many proteins remain functionally unannotated. Sequence alignment (SA)
uncovers missing annotations by transferring functional knowledge between
species' sequence-conserved regions. Because SA is imperfect, network alignment
(NA) complements SA by transferring functional knowledge between conserved
biological network, rather than just sequence, regions of different species.
Existing NA assumes that it is topological similarity (isomorphic-like
matching) between network regions that corresponds to the regions' functional
relatedness. However, we recently found that functionally unrelated proteins
are almost as topologically similar as functionally related proteins. So, we
redefined NA as a data-driven framework, TARA, which learns from network and
protein functional data what kind of topological relatedness (rather than
similarity) between proteins corresponds to the proteins' functional
relatedness. TARA used topological information (within each network) but not
sequence information (between proteins across networks). Yet, its alignments
yielded higher protein functional prediction accuracy than alignments of
existing NA methods, even those that used both topological and sequence
information. Here, we propose TARA++ that is also data-driven, like TARA and
unlike other existing methods, but that uses across-network sequence
information on top of within-network topological information, unlike TARA. To
deal with the within-and-across-network analysis, we adapt social network
embedding to the problem of biological NA. TARA++ outperforms protein
functional prediction accuracy of existing methods.


  In exponentially proliferating populations of microbes, the population
typically doubles at a rate less than the average doubling time of a
single-cell due to variability at the single-cell level. It is known that the
distribution of generation times obtained from a single lineage is, in general,
insufficient to determine a population's growth rate. Is there an explicit
relationship between observables obtained from a single lineage and the
population growth rate? We show that a population's growth rate can be
represented in terms of averages over isolated lineages. This lineage
representation is related to a large deviation principle that is a generic
feature of exponentially proliferating populations. Due to the large deviation
structure of growing populations, the number of lineages needed to obtain an
accurate estimate of the growth rate depends exponentially on the duration of
the lineages, leading to a non-monotonic convergence of the estimate, which we
verify in both synthetic and experimental data sets.


  OBJECTIVE: To better understand clinician needs and preferences for the
display of pharmacogenomic (PGx) information in clinical decision support (CDS)
tools.
  MATERIALS AND METHODS: We developed a semi-structured interview to collect
feedback and preferences in six key areas of PGx CDS design, from clinicians
who had prior experience with live PGx CDS tools. Eight clinicians from
Northwestern Medicine's (NM) General Internal Medicine clinic participated in
the study.
  RESULTS: Clinicians expressed preference for interruptive pop-up alerts
during order entry, brief descriptions of relevant drug-gene interactions, and
a clear and specific recommended alternative course of action when a medication
is contraindicated. They did not wish to see detailed genetic data, preferring
phenotypic information predicted from the genotype. Nor did they wish to be
interrupted when genetic test results do not indicate a change in treatment
plan. Clinicians reported little familiarity with Clinical Pharmacogenetic
Implementation Consortium prescribing recommendations but reported trusting
recommendations of their professional societies and resources like UpToDate.
Analysis of unstructured comments concurred with structured results, indicating
a general uncertainty among participants around how to interpret and apply PGx
information in practice.
  DISCUSSION: Results point to several underlying principles that can inform
future PGx CDS alert designs: Be Specific and Actionable; Be Brief; Display
Phenotypes not Genotypes; Rely on Sources Clinicians Already Trust; and, Be
Adaptable to Learning Effects.
  CONCLUSION: This study is part of a broader socio-technical design approach
to PGx CDS design underway at NM and provides a baseline for future PGx CDS
development. Designs based on these results have the potential to improve
clinician education and adherence levels, and to improve patient outcomes.


  OBJECTIVE: To design and evaluate new pharmacogenomic (PGx) clinical decision
support (CDS) alerts, built to adhere to PGx CDS design principles developed
through socio-technical approaches.
  MATERIALS AND METHODS: Based on previously identified design principles, we
created 11 new PGx CDS alert designs and developed an interactive web
application containing realistic clinical scenarios and user workflows that
mimicked a real-world EHR system. We recruited General Internal Medicine and
Cardiology clinicians from Northwestern Medicine and recorded their
interactions with the original and new designs. We measured clinician response,
satisfaction, speed, and confidence through questionnaires and analysis of the
recordings.
  RESULTS: The study included 12 clinicians. Participants were significantly
more satisfied (p=0.0000001), faster (p=0.009), and more confident (p<.05) with
the new designs than the original ones. The study lacked statistical power to
determine whether prescribing accuracy was improved, but participants were no
less accurate, and clinical actions were more concordant with alert
interactions (p=0.004) with the new designs. We found a significant learning
curve associated with the original designs, which was eliminated with the new
designs.
  DISCUSSION: This study successfully demonstrates that socio-technical and
user-centered design techniques can improve PGx CDS alert designs. Best
practices for PGx CDS design are limited in the literature, with few
effectiveness studies available. These results can help guide future PGx CDS
implementations to be more clinician friendly and less time-consuming.
  CONCLUSION: The results of this study support the PGx CDS design principles
we proposed in previous work. As a next step, the new designs should be
implemented in a live setting for further validation.


  In this work, two mathematical models for malaria under resistance are
presented. More precisely, the first model shows the interaction between humans
and mosquitoes inside a patch under infection of malaria when the human
population is resistant to antimalarial drug and mosquitoes population is
resistant to insecticides. For the second model, human-mosquitoes population
movements in two patches is analyzed under the same malaria transmission
dynamic established in one patch. For a single patch, existence and stability
conditions for the equilibrium solutions in terms of the local basic
reproductive number are developed. These results reveal the existence of a
forward bifurcation and the global stability of disease-free equilibrium. In
the case of two patches, a theoretical and numerical framework on sensitivity
analysis of parameters is presented. After that, the use of antimalarial drugs
and insecticides are incorporated as control strategies and an optimal control
problem is formulated. Numerical experiments are carried out in both models to
show the feasibility of our theoretical results.


  We investigate the measurement of DNA supercoiling density ($\sigma$) along
chromosomes using interaction frequencies between DNA and DNA-anchored clusters
of proteins. Specifically, we show how the physics of DNA supercoiling leads,
in bacteria, to the quantitative modeling of binding properties of ParB
proteins around their centromere-like site, {\it parS}. Using this framework,
we provide an upper bound for $\sigma$ in the {\it Escherichia coli}
chromosome, consistent with plasmid values, and offer a proof of concept for a
high accuracy measurement. To reach these conclusions, we revisit the problem
of the formation of ParB clusters. We predict, in particular, that they result
from a non-equilibrium, stationary balance between an influx of produced
proteins and an outflux of excess proteins, i.e., they behave like liquid-like
protein condensates with unconventional ``leaky'' boundaries.


  The focus of this viewpoint is to identify, in the era of atomistic
resolution cryo-electron microscopy data, the areas in which computational
modelling and molecular simulations will bring valuable contributions to
structural biologists and to give an overview of some of the existing efforts
in this direction.


  Mikhail Shapiro's lab investigated a promising new method for non-invasive
control of calcium currents in individual cells in the nervous system by the
selective heating of nanoparticles and show that simple physical laws, properly
applied, explain what is happening, and so can be a foundation for constructing
improved methods and techniques. They use custom built instrumentation and
detailed quantitative measurement to show that special surface properties are
not needed to explain their results. Their work is taken as an example of the
need for quantitative measurement in biophysics and biology in general.
Classical examples are cited and the argument is made that the success of
structural and molecular biology has hidden the need for quantitative
measurements and controls. Semiconductor and computational electronics is cited
as a science even more successful than structural and molecular biology that
depends on quantitative measurement, controls, and analysis.


  Reconstructing pathogen dynamics from genetic data as they become available
during an outbreak or epidemic represents an important statistical scenario in
which observations arrive sequentially in time and one is interested in
performing inference in an 'online' fashion. Widely-used Bayesian phylogenetic
inference packages are not set up for this purpose, generally requiring one to
recompute trees and evolutionary model parameters de novo when new data arrive.
To accommodate increasing data flow in a Bayesian phylogenetic framework, we
introduce a methodology to efficiently update the posterior distribution with
newly available genetic data. Our procedure is implemented in the BEAST 1.10
software package, and relies on a distance-based measure to insert new taxa
into the current estimate of the phylogeny and imputes plausible values for new
model parameters to accommodate growing dimensionality. This augmentation
creates informed starting values and re-uses optimally tuned transition kernels
for posterior exploration of growing data sets, reducing the time necessary to
converge to target posterior distributions. We apply our framework to data from
the recent West African Ebola virus epidemic and demonstrate a considerable
reduction in time required to obtain posterior estimates at different time
points of the outbreak. Beyond epidemic monitoring, this framework easily finds
other applications within the phylogenetics community, where changes in the
data -- in terms of alignment changes, sequence addition or removal -- present
common scenarios that can benefit from online inference.


  Understanding how genotypes map onto phenotypes, fitness, and eventually
organisms is arguably the next major missing piece in a fully predictive theory
of evolution. We refer to this generally as the problem of the
genotype-phenotype map. Though we are still far from achieving a complete
picture of these relationships, our current understanding of simpler questions,
such as the structure induced in the space of genotypes by sequences mapped to
molecular structures, has revealed important facts that deeply affect the
dynamical description of evolutionary processes. Empirical evidence supporting
the fundamental relevance of features such as phenotypic bias is mounting as
well, while the synthesis of conceptual and experimental progress leads to
questioning current assumptions on the nature of evolutionary dynamics-cancer
progression models or synthetic biology approaches being notable examples. This
work delves into a critical and constructive attitude in our current knowledge
of how genotypes map onto molecular phenotypes and organismal functions, and
discusses theoretical and empirical avenues to broaden and improve this
comprehension. As a final goal, this community should aim at deriving an
updated picture of evolutionary processes soundly relying on the structural
properties of genotype spaces, as revealed by modern techniques of molecular
and functional analysis.


  In this paper, we propose a novel dynamical system with time delay to
describe the outbreak of 2019-nCoV in China. One typical feature of this
epidemic is that it can spread in latent period, which is therefore described
by the time delay process in the differential equations. The accumulated
numbers of classified populations are employed as variables, which is
consistent with the official data and facilitates the parameter identification.
The numerical methods for the prediction of outbreak of 2019-nCoV and parameter
identification are provided, and the numerical results show that the novel
dynamic system can well predict the outbreak trend so far. Based on the
numerical simulations, we suggest that the transmission of individuals should
be greatly controlled with high isolation rate by the government.


  In order to predict the development trend of the 2019 coronavirus
(2019-nCov), we established an prediction model to predict the number of
diagnoses case in China except Hubei Province. From January 25 to January 29,
2020, we optimized 6 prediction models, 5 of them based on the number of
medical observations to predicts the peak time of confirmed diagnosis will
appear on the period of morning of January 29 from 24:00 to February 2 before 5
o'clock 24:00. Then we tracked the data from 24 o'clock on January 29 to 24
o'clock on January 31, and found that the predicted value of the data on the
3rd has a small deviation from the actual value, and the actual value has
always remained within the range predicted by the comprehensive prediction
model 6. Therefore we discloses this finding and will continue to track whether
this pattern can be maintained for longer. We believe that the changes medical
observation case number may help to judge the trend of the epidemic situation
in advance.


  Diffusion-limited association reactions are ubiquitous in nature. They are
particularly important for biological reactions, where the reaction rates are
often determined by the diffusive transport of the molecules on two-dimensional
surfaces, such as the cell membrane. The peculiarities of diffusion on
two-dimensional surfaces may lead to nontrivial reaction kinetics, such as
concentration dependent rate of association between two molecules. However,
traditionally, the kinetics of biomolecular association reactions has been
modeled using the law of mass action, which assumes that the rate of reaction
is a concentration independent constant. In this paper, using multiscale
molecular simulation, we investigate the concentration dependence of
diffusion-limited association reactions on 2D surfaces. In particular, we
quantify the influence of short-ranged pair interactions on the concentration
dependence of the reaction rates and codify it in an empirical law. Using this
law in a chemical kinetic model, we find that the the steady state behaviors of
simple chemical systems are drastically modified by the presence of
concentration dependent rates. In particular, we find that it leads to
suppression of intrinsic noise in dimerization reaction and destabilizes robust
oscillation in Lotka-Volterra predator-prey systems. In fact, we see a
transition from robust to fine-tuned behavior in the latter. In addition, we
show that concentration dependent reaction rates arise naturally in stochastic
predator-prey systems due to intrinsic noise. We comment on the consequences of
these results and discuss their implications in the modeling of complex
chemical and biological systems. In particular, we comment on the range of
validity of the law of mass action, which is a staple in all theoretical
modeling of these systems.


  Phylogenetic methods have broad potential in linguistics beyond tree
inference. Here, we show how a phylogenetic approach opens the possibility of
gaining historical insights from entirely new kinds of linguistic data--in this
instance, statistical phonotactics. We extract phonotactic data from 111
Pama-Nyungan vocabularies and apply tests for phylogenetic signal, quantifying
the degree to which the data reflect phylogenetic history. We test three
datasets: (1) binary variables recording the presence or absence of biphones
(two-segment sequences) in a lexicon (2) frequencies of transitions between
segments, and (3) frequencies of transitions between natural sound classes.
Australian languages have been characterized as having a high degree of
phonotactic homogeneity. Nevertheless, we detect phylogenetic signal in all
datasets. Phylogenetic signal is greater in finer-grained frequency data than
in binary data, and greatest in natural-class-based data. These results
demonstrate the viability of employing a new source of readily extractable data
in historical and comparative linguistics.


  Nestedness is a property of interaction networks widely observed in natural
mutualistic communities. Despite a widespread interest on this pattern, no
general consensus exists on how to measure it. Instead, several metrics aiming
at quantifying nestedness, based on different but not necessarily independent
properties of the networks, coexist in the literature blurring the comparison
between ecosystems. In this work, we present a detailed critical study of the
behavior of six popular nestedness metrics and the variants of two of them. In
order to evaluate their performance, we compare the obtained values of the
nestedness of a large set of real networks among them and against a maximum
entropy and maximum likelihood null model. Our results point out, first, that
the metrics do not rank the degree of nestedness of networks universally.
Furthermore, several metrics show significant undesired dependencies on the
network properties considered. The study of these dependencies allows us to
understand some of the systematic shifts between the real values of nestedness
and the average over the null model. This paper intends to provide readers with
a critical guide on how to measure nestedness patterns, by explaining the
functioning of six standard metrics and two of its variants, and then
disclosing its qualities and flaws. By doing so, we also aim to extend the
application of the recently proposed null models based on maximum entropy to
the still largely unexplored area of ecological networks. Finally, to
complement the guide, we provide a fully-documented repository named nullnest
which gathers the codes to produce the null model and calculate the nestedness
index -- both the real value and the null expectation -- using the studied
metrics. The repository contains, moreover, the main results of the null model
applied to a large dataset of more than 200 bipartite networks.


  NeuroEvolution is one of the most competitive evolutionary learning
frameworks for designing novel neural networks for use in specific tasks, such
as logic circuit design and digital gaming. However, the application of
benchmark methods such as the NeuroEvolution of Augmenting Topologies (NEAT)
remains a challenge, in terms of their computational cost and search time
inefficiency. This paper advances a method which incorporates a type of
topological edge coding, named Reverse Encoding Tree (RET), for evolving
scalable neural networks efficiently. Using RET, two types of approaches --
NEAT with Binary search encoding (Bi-NEAT) and NEAT with Golden-Section search
encoding (GS-NEAT) -- have been designed to solve problems in benchmark
continuous learning environments such as logic gates, Cartpole, and Lunar
Lander, and tested against classical NEAT and FS-NEAT as baselines.
Additionally, we conduct a robustness test to evaluate the resilience of the
proposed NEAT algorithms. The results show that the two proposed strategies
deliver improved performance, characterized by (1) a higher accumulated reward
within a finite number of time steps; (2) using fewer episodes to solve
problems in targeted environments, and (3) maintaining adaptive robustness
under noisy perturbations, which outperform the baselines in all tested cases.
Our analysis also demonstrates that RET expends potential future research
directions in dynamic environments. Code is available from
https://github.com/HaolingZHANG/ReverseEncodingTree.


  We present the second generation of a rule-based language called Biochemical
Space Language (BCSL) that combines the advantages of different approaches and
thus makes an effort to overcome several problems with existing solutions. The
key aspect of the language is the level of abstraction it uses, which allows
scalable and compact hierarchical specification of biochemical entities. This
abstraction enables unique analysis techniques to reason about properties of
models written in the language on the semantic and syntactic level.


  Compartmentalization of biochemical processes underlies all biological
systems, from the organelle to the tissue scale. Theoretical models to study
the interplay between noisy reaction dynamics and compartmentalization are
sparse, and typically very challenging to analyze computationally. Recent
studies have made progress towards addressing this problem in the context of
concrete biological systems but general approaches remain lacking. In this work
we propose a mathematical framework based on counting processes that allows us
to study compartment populations with arbitrary interactions and internal
biochemistry. We provide an efficient description of the population dynamics in
terms of differential equations which capture moments of the population and
their variability. We demonstrate the relevance of our approach using several
case studies inspired by biological systems at different scales.


  Disease and pest control largely rely on pesticides use and progress still
remains to be made towards more sustainable practices. Pest Monitoring Networks
(PMNs) can provide useful information for improving crop protection by
restricting pesticide use to the situations that best require it. However, the
efficacy of a PMN to control pests may depend on its spatial density and
space/time sampling balance. Furthermore the best trade-off between the
monitoring effort and the impact of the PMN information may be pest dependent.
We developed a generic simulation model that links PMN information to treatment
decisions and pest dynamics. We derived the number of treatments, the epidemic
extension and the global gross margin for different families of pests. For
soil-borne pathogens and weeds, we found that increasing the spatial density of
a PMN significantly decreased the number of treatments (up to 67\%), with an
only marginal increase in infection. Considering past observations had a
second-order effect (up to a 13\% decrease). For the spatial scale of our
study, the PMN information had practically no influence in the case of insects.
The next step is to go beyond PMN analysis to design and chose among
sustainable management strategies at the landscape scale.


  Today's society faces widening disagreement and conflicts among constituents
with incompatible views. Escalated views and opinions are seen not only in
radical ideology or extremism but also in many other scenes of our everyday
life. Here we show that widening disagreement among groups may be linked to the
advancement of information communication technology, by analyzing a
mathematical model of population dynamics in a continuous opinion space. We
adopted the interaction kernel approach to model enhancement of people's
information gathering ability and introduced a generalized non-local gradient
as individuals' perception kernel. We found that the characteristic distance
between population peaks becomes greater as the wider range of opinions becomes
available to individuals or the greater attention is attracted to opinions
distant from theirs. These findings may provide a possible explanation for why
disagreement is growing in today's increasingly interconnected society, without
attributing its cause only to specific individuals or events.


  The use of network based approaches to model and analyse large datasets is
currently a growing research field. For instance in biology and medicine,
networks are used to model interactions among biological molecules as well as
relations among patients. Similarly, data coming from social networks can be
trivially modelled by using graphs. More recently, the use of dual networks
gained the attention of researchers. A dual network model uses a pair of graphs
to model a scenario in which one of the two graphs is usually unweighted (a
network representing physical associations among nodes) while the other one is
edge-weighted (a network representing conceptual associations among nodes). In
this paper we focus on the problem of finding the Densest Connected sub-graph
(DCS) having the largest density in the conceptual network which is also
connected in the physical network. The problem is relevant but also
computationally hard, therefore the need for introducing of novel algorithms
arises. We formalise the problem and then we map DCS into a graph alignment
problem. Then we propose a possible solution. A set of experiments is also
presented to support our approach.


  DNA and proteins are chiral: their three-dimensional structure cannot be
superimposed with its mirror image. Circular dichroism spectroscopy is widely
used to characterize chiral compounds, but data interpretation is difficult in
the case of mixtures. We recorded for the first time the electronic circular
dichroism spectra of DNA helices separated in a mass spectrometer. We
electrosprayed guanine-rich strands having various secondary structures as
negative ions, irradiated them with a laser, and measured the difference in
electron photodetachment efficiency between left and right circularly polarized
light. The reconstructed circular dichroism ion spectra resemble the solution
ones, thereby allowing us to assign the DNA helical topology. The ability to
measure circular dichroism directly on biomolecular ions expands the
capabilities of mass spectrometry for structural analysis.


  The emergence of R, a freely available data analysis environment, brought to
the researcher in any science field a set of well-concerted instruments of
immense power and low cost. In botany and zoology, these instruments could be
used, for example, to speed up work in two distant but related fields: analysis
of DNA markers and preparation of natural history manuals. Both of these tasks
require a significant amount of monotonous work, which could be automated with
software.
  I developed "Ripeline and "Rmanual," two highly customizable R-based
applications, designed with a goal of simplicity, reproducibility, and
effectiveness. Ripeline is a pipeline that allows for a continuously updated
analysis of multiple DNA markers. Rmanual is a "living book" which allows the
creation and continuous update of manuals and checklists.
  Comparing with more traditional ways of DNA marker analysis and manual
preparation, Ripeline and Rmanual allow for a significant reduction of time,
which is usually spent doing repetitive tasks. They also provide tools which
can be used in a broad spectrum of further applications.


  Electronically-active organic molecules have demonstrated great promise as
novel soft materials for energy harvesting and transport. Self-assembled
nanoaggregates formed from $\pi$-conjugated oligopeptides composed of an
aromatic core flanked by oligopeptide wings offer emergent optoelectronic
properties within a water soluble and biocompatible substrate. Nanoaggregate
properties can be controlled by tuning core chemistry and peptide composition,
but the sequence-structure-function relations remain poorly characterized. In
this work, we employ coarse-grained molecular dynamics simulations within an
active learning protocol employing deep representational learning and Bayesian
optimization to efficiently identify molecules capable of assembling pseudo-1D
nanoaggregates with good stacking of the electronically-active $\pi$-cores. We
consider the DXXX-OPV3-XXXD oligopeptide family, where D is an Asp residue and
OPV3 is an oligophenylene vinylene oligomer (1,4-distyrylbenzene), to identify
the top performing XXX tripeptides within all 20$^3$ = 8,000 possible
sequences. By direct simulation of only 2.3% of this space, we identify
molecules predicted to exhibit superior assembly relative to those reported in
prior work. Spectral clustering of the top candidates reveals new design rules
governing assembly. This work establishes new understanding of DXXX-OPV3-XXXD
assembly, identifies promising new candidates for experimental testing, and
presents a computational design platform that can be generically extended to
other peptide-based and peptide-like systems.


  Robustness and evolvability are essential properties to the evolution of
biological networks. To determine if a biological network is robust and/or
evolvable, it is required to compare its functions before and after mutations.
However, this sometimes takes a high computational cost as the network size
grows. Here we develop a predictive method to estimate the robustness and
evolvability of biological networks without an explicit comparison of
functions. We measure antifragility in Boolean network models of biological
systems and use this as the predictor. Antifragility occurs when a system
benefits from external perturbations. By means of the differences of
antifragility between the original and mutated biological networks, we train a
convolutional neural network (CNN) and test it to classify the properties of
robustness and evolvability. We found that our CNN model successfully
classified the properties. Thus, we conclude that our antifragility measure can
be used as a predictor of the robustness and evolvability of biological
networks.


  Predicting stochastic cellular dynamics as emerging from the mechanistic
models of molecular interactions is a long-standing challenge in systems
biology: low-level chemical reaction network (CRN) models give raise to a
highly-dimensional continuous-time Markov chain (CTMC) which is computationally
demanding and often prohibitive to analyse in practice. A recently proposed
abstraction method uses deep learning to replace this CTMC with a discrete-time
continuous-space process, by training a mixture density deep neural network
with traces sampled at regular time intervals (which can obtained either by
simulating a given CRN or as time-series data from experiment). The major
advantage of such abstraction is that it produces a computational model that is
dramatically cheaper to execute, while preserving the statistical features of
the training data. In general, the abstraction accuracy improves with the
amount of training data. However, depending on a CRN, the overall quality of
the method -- the efficiency gain and abstraction accuracy -- will also depend
on the choice of neural network architecture given by hyper-parameters such as
the layer types and connections between them. As a consequence, in practice,
the modeller would have to take care of finding the suitable architecture
manually, for each given CRN, through a tedious and time-consuming
trial-and-error cycle. In this paper, we propose to further automatise deep
abstractions for stochastic CRNs, through learning the optimal neural network
architecture along with learning the transition kernel of the abstract process.
Automated search of the architecture makes the method applicable directly to
any given CRN, which is time-saving for deep learning experts and crucial for
non-specialists. We implement the method and demonstrate its performance on a
number of representative CRNs with multi-modal emergent phenotypes.


  Lung cancer has a high rate of recurrence in early-stage patients. Predicting
the post-surgical recurrence in lung cancer patients has traditionally been
approached using single modality information of genomics or radiology images.
We investigate the potential of multimodal fusion for this task. By combining
computed tomography (CT) images and genomics, we demonstrate improved
prediction of recurrence using linear Cox proportional hazards models with
elastic net regularization. We work on a recent non-small cell lung cancer
(NSCLC) radiogenomics dataset of 130 patients and observe an increase in
concordance-index values of up to 10%. Employing non-linear methods from the
neural network literature, such as multi-layer perceptrons and visual-question
answering fusion modules, did not improve performance consistently. This
indicates the need for larger multimodal datasets and fusion techniques better
adapted to this biological setting.


  The neuromuscular junction represents a relevant substrate for revealing
important biophysical mechanisms of synaptic transmission. In this context,
calcium ions are important in the synapse machinery, providing the nervous
impulse transmission to the muscle fiber. In this work, we carefully
investigated whether intervals of spontaneous electrical activity, recorded in
seven different calcium concentrations, conform Newcomb-Benford law. Our
analysis revealed that electrical discharge of neuromuscular junction obeys the
expected values for Newcomb-Benford law for first and second digits, while
first-two digits do not perfectly follows the law. We next examined previous
theoretical studies, establishing a relation between the law and lognormal and
Weibull distributions. We showed that Weibull distribution is more appropriate
to fit the intervals as compared to lognormal distribution. Altogether, the
present findings strongly suggest that spontaneous activity is a base-scale
invariant phenomenon.


  State-space models (SSMs) are an important modeling framework for analyzing
ecological time series. These hierarchical models are commonly used to model
population dynamics, animal movement, and capture-recapture data, and are now
increasingly being used to model other ecological processes. SSMs are popular
because they are flexible and they model the natural variation in ecological
processes separately from observation error. Their flexibility allows
ecologists to model continuous, count, binary, and categorical data with linear
or nonlinear processes that evolve in discrete or continuous time. Modeling the
two sources of stochasticity separately allows researchers to differentiate
between biological variation (e.g., in birth processes) and imprecision in the
sampling methodology, and generally provides better estimates of the ecological
quantities of interest than if only one source of stochasticity is directly
modeled. Since the introduction of SSMs, a broad range of fitting procedures
have been proposed. However, the variety and complexity of these procedures can
limit the ability of ecologists to formulate and fit their own SSMs. We provide
the knowledge for ecologists to create SSMs that are robust to common, and
often hidden, estimation problems, and the model selection and validation tools
that can help them assess how well their models fit their data. In this paper,
we present a review of SSMs that will provide a strong foundation to ecologists
interested in learning about SSMs, introduce new tools to veteran SSM users,
and highlight promising research directions for statisticians interested in
ecological applications. The review is accompanied by an in-depth tutorial that
demonstrates how SSMs models can be fitted and validated in R. Together, the
review and tutorial present an introduction to SSMs that will help ecologists
to formulate, fit, and validate their models.


  This manuscript contains nothing new, but synthesizes known results: For the
theoretical population geneticist with a probabilistic background, we provide a
summary of some key results on stochastic differential equations. For the
evolutionary game theorist, we give a new perspective on the derivations of
results obtained when using discrete birth-death processes. For the theoretical
biologist familiar with deterministic modeling, we outline how to derive and
work with stochastic versions of classical ecological and evolutionary
processes.


  A novel technique was demonstrated that overcome important drawbacks to
crosslink cells by irradiation with ultrashort UV laser pulses
(L-crosslinking). To use this technique coupled to Chromatin
ImmunoPrecipitation (ChIP) in a high throughput context, a pre-screening fast
method needs to be implemented to set up suitable irradiation conditions of the
cell sample for efficient L-crosslinking with no final and long ChIP analysis.
Here a fast method is reported where living human cells have been first
transfected with a vector coding for Estrogen Receptor {\alpha} (ER{\alpha}),
linked to Green Florescent protein (ER{\alpha}-GFP), so that the well-known
interaction between the Estrogen Receptor Elements (ERE) region of the cell DNA
and the ER{\alpha} protein can be detected by studying the fluorometric
response of the irradiated cells. The damage induced to cells by UV irradiation
is characterized by looking at DNA integrity, proteins stability and cellular
viability. A second novel approach is presented to analyze or re-visit DNA and
RNA sequences and their molecular configurations. This approach is based on
methods derived from Chern-Simons super-gravity adapted to describe mutations
in DNA/RNA strings, as well as interactions between nucleic acids. As a
preliminary case we analyze the KRAS human gene sequence and some of its
mutations. Interestingly, our model shows how the Chern-Simons current are
capable to characterize the mutations within a sequence, in particular giving a
quantitative indication of the mutation likelihood.


  The monitoring, analysis and prediction of epidemic spread in the region
require the construction of mathematical model, big data processing and
visualization because the amount of population and the size of the region could
be huge. One of the important steps is refinement of mathematical model, i.e.
determination of initial data and coefficients of system of differential
equations which describe the epidemiology processes. We analyze numerical
method for solving inverse problem of epidemiology based on genetic algorithm
and traditional optimization ideas. Numerical results are applied to analysis
and prediction of epidemic situation in regions of Russian Federation, Republic
of Kazakhstan and People's Republic of China. Due to a great amount of data we
use a special Geo-information system for visualization of epidemic process,
i.e. a special software named Digital Earth.


  Genomics is changing our understanding of humans, evolution, diseases, and
medicines to name but a few. As sequencing technology is developed collecting
DNA sequences takes less time thereby generating more genetic data every day.
Today the rate of generating genetic data is outpacing the rate of computation
power growth. Current sequencing machines can sequence 50 humans genome per
day; however, aligning the read sequences against a reference genome and
assembling the genome will take 1300 CPU hours. The main step in constructing
the genome is aligning the reads against a reference genome. Numerous
accelerators have been proposed to accelerate the DNA alignment process.
Providing massive parallelism, FPGA-based accelerators have shown great
performance in accelerating DNA alignment algorithms. Additionally, FPGA-based
accelerators provide better energy efficiency than general-purpose processors.
In this survey, we introduce three main DNA alignment algorithms and FPGA-based
implementation of these algorithms to accelerate the DNA alignment. We also,
compare these three alignment categories and show how accelerators are
developing during the time.


  How to model the 2019 CoronaVirus (2019-nCov) spread in China is one of the
most urgent and interesting problems in applied mathematics. In this paper, we
propose a novel time delay dynamic system with external source to describe the
trend of local outbreak for the 2019-nCoV. The external source is introduced in
the newly proposed dynamic system, which can be considered as the suspected
people travel to different areas. The numerical simulations exhibit the dynamic
system with the external source is more reliable than the one without it, and
the rate of isolation is extremely important for controlling the increase of
cumulative confirmed people of 2019-nCoV. Based on our numerical simulation
results with the public data, we suggest that the local government should have
some more strict measures to maintain the rate of isolation. Otherwise the
local cumulative confirmed people of 2019-nCoV might be out of control.


  Resource support between individuals is of particular importance in
controlling or mitigating epidemic spreading, especially during pandemics.
Whereas there remains the question of how we can protect ourselves from being
infected while helping others by donating resources in fighting against the
epidemic. To answer the question, we propose a novel resource allocation model
by considering the awareness of self-protection of individuals. In the model, a
tuning parameter is introduced to quantify the reaction strength of individuals
when they are aware of the disease. And then, a coupled model of resource
allocation and disease spreading is proposed to study the impact of
self-awareness on resource allocation and, its impact on the dynamics of
epidemic spreading. Through theoretical analysis and extensive Monte Carlo
simulations, we find that in the stationary state, the system converges to two
states: the whole healthy or the completely infected, which indicates an abrupt
increase in the prevalence when there is a shortage of resources. More
importantly, we find that too cautious and too selfless for the people during
the outbreak of an epidemic are both not suitable for disease control. Through
extensive simulations, we find the optimal point, at which there is a maximum
value of the epidemic threshold, and an outbreak can be delayed to the greatest
extent. At last, we study further the effects of network structure on the
coupled dynamics. We find that the degree heterogeneity promotes the outbreak
of disease, and the network structure does not alter the optimal phenomenon in
behavior response.


  The affinity of antibodies (Abs) produced in vivo for their target antigens
(Ags) is typically well below the maximum affinity possible. Nearly 25 years
ago, Foote and Eisen explained how an 'affinity ceiling' could arise from
constraints associated with the acquisition of soluble antigen by B cells.
However, recent studies have shown that B cells in germinal centers (where Ab
affinity maturation occurs) acquire Ag not in soluble form but presented as
receptor-bound immune complexes on follicular dendritic cells (FDCs). How the
affinity ceiling arises in such a scenario is unclear. Here, we argue that the
ceiling arises from the weakest link of the chain of protein complexes that
bridges B cells and FDCs and is broken during Ag acquisition. This hypothesis
explains the affinity ceiling realized in vivo and suggests that strengthening
the weakest link could raise the ceiling and improve Ab responses.


  Climate change threatens biodiversity directly by influencing biophysical
variables that drive species' geographic distributions and indirectly through
socio-economic changes that influence land use patterns, driven by global
consumption, production and climate. To date, no detailed analyses have been
produced that assess the relative importance of, or interaction between, these
direct and indirect climate change impacts on biodiversity at large scales.
Here, we apply a new integrated modelling framework to quantify the relative
influence of biophysical and socio-economically mediated impacts on avian
species in Vietnam and Australia. We find that socio-economically mediated
impacts on suitable ranges are largely outweighed by biophysical impacts, but
global shifts of production are likely to result in adverse impacts on habitats
worldwide. By translating economic futures and shocks into spatially explicit
predictions of biodiversity change, we now have the power to analyse in a
consistent way outcomes for nature and people of any change to policy,
regulation, trading conditions or consumption trend at any scale from
sub-national to global.


  In this paper, we establish a reaction-diffusion-advection partial
differential equation model to describe the growth of algae depending on both
nutrients and inorganic carbon in a poorly mixed water column. Nutrients from
the water bottom and inorganic carbon from the water surface form an asymmetric
resource supply mechanism on the algal growth. The existence and stability of
semi-trivial steady state and coexistence steady state of the model are proved,
and a threshold condition for the regime shift from extinction to survival of
algae is established. The influence of environmental parameters on the vertical
distribution of algae is investigated in the water column. It is shown that the
vertical distribution of algae can exhibit many different profiles under the
joint limitation of nutrients and inorganic carbon.


  As the infection of 2019-nCoV coronavirus is quickly developing into a global
pneumonia epidemic, careful analysis of its transmission and cellular
mechanisms is sorely needed. In this report, we re-analyzed the computational
approaches and findings presented in two recent manuscripts by Ji et al.
(https://doi.org/10.1002/jmv.25682) and by Pradhan et al.
(https://doi.org/10.1101/2020.01.30.927871), which concluded that snakes are
the intermediate hosts of 2019-nCoV and that the 2019-nCoV spike protein
insertions shared a unique similarity to HIV-1. Results from our
re-implementation of the analyses, built on larger-scale datasets using
state-of-the-art bioinformatics methods and databases, do not support the
conclusions proposed by these manuscripts. Based on our analyses and existing
data of coronaviruses, we concluded that the intermediate hosts of 2019-nCoV
are more likely to be mammals and birds than snakes, and that the "novel
insertions" observed in the spike protein are naturally evolved from bat
coronaviruses.


  One of the chief paradoxes of molecular oxygen (O$_2$) is that it is an
essential requirement for multicellular eukaryotes on Earth while
simultaneously posing a threat to their survival via the formation of reactive
oxygen species. In this paper, the constraints imposed by O$_2$ on Earth-like
complex life are invoked to explore whether worlds with abiotic O$_2$
inventories can harbor such organisms. By taking the major O$_2$ sources and
sinks of Earth-like planets into account using a simple model, it is suggested
that worlds that receive time-averaged X-ray and extreme ultraviolet fluxes
that are $\gtrsim 10$ times higher than Earth might not be capable of hosting
complex lifeforms because the photolysis of molecules such as water may lead to
significant O$_2$ buildup. Methods for testing this hypothesis by searching for
anticorrelations between biosignatures and indicators of abiotic O$_2$
atmospheres are described. In the event, however, that life successfully adapts
to high-oxygen environments, these worlds could permit the evolution of large
and complex organisms.


  The novel coronavirus (2019-nCoV) is a recently emerged human pathogen that
has spread widely since January 2020. Initially, the basic reproductive number,
R0, was estimated to be 2.2 to 2.7. Here we provide a new estimate of this
quantity. We collected extensive individual case reports and estimated key
epidemiology parameters, including the incubation period. Integrating these
estimates and high-resolution real-time human travel and infection data with
mathematical models, we estimated that the number of infected individuals
during early epidemic double every 2.4 days, and the R0 value is likely to be
between 4.7 and 6.6. We further show that quarantine and contact tracing of
symptomatic individuals alone may not be effective and early, strong control
measures are needed to stop transmission of the virus.


  Active ion transport is very critical for living cells to maintain and
regular internal and external environment. It is known that adenosine
5'-triphosphate (ATP) is a general energy source that applies bond energy for
pumps to overcome ion concentration gradient. In this study, we introduce a
novel thermodynamic model for active ion transport, which allows the pump to
act as a 'Maxwell's demon', and also conforms to the second law of
thermodynamics. Transport against the gradient is caused by the thermodynamic
fluctuation with information recording in ATP stream. Besides, exhaustive
experiments about Na-K ATPase performed in the red cell ghost system are
reviewed to verify this model, and all the results can support the solvable
theory about the working mechanism of demon pumps. Our findings indicate the
possible role of ATP as an information carrier but not energy currency. The
high-energy phosphate bonds can improve the efficiency of information recording
in relevance to ion transport and may mainly convert to heat.


  The analysis of coevolution of residues in homologous proteins is a powerful
tool to predict their native conformation. The standard framework in which
coevolutionary analysis is usually worked out is that of equilibrium Potts
models, assuming that proteins have evolved for enough time to reach
thermodynamic equilibrium in sequence space. Here we propose a model to
describe correlations in sequences based on an explicit description of the
evolutionary kinetics of proteins. We show that this procedure improves the
correct prediction of native contacts with respect to equilibrium--based
models.


  The complete characterization of enzymatic activities between molecules
remains incomplete, hindering biological engineering and limiting biological
discovery. We develop in this work a technique, Enzymatic Link Prediction
(ELP), for predicting the likelihood of an enzymatic transformation between two
molecules. ELP models enzymatic reactions catalogued in the KEGG database as a
graph. ELP is innovative over prior works in using graph embedding to learn
molecular representations that capture not only molecular and enzymatic
attributes but also graph connectivity.
  We explore both transductive (test nodes included in the training graph) and
inductive (test nodes not part of the training graph) learning models. We show
that ELP achieves high AUC when learning node embeddings using both graph
connectivity and node attributes. Further, we show that graph embedding for
predicting enzymatic links improves link prediction by 24% over
fingerprint-similarity-based approaches. To emphasize the importance of graph
embedding in the context of biochemical networks, we illustrate how graph
embedding can also guide visualization.
  The code and datasets are available through
https://github.com/HassounLab/ELP.


  We present the findings of "The Alzheimer's Disease Prediction Of
Longitudinal Evolution" (TADPOLE) Challenge, which compared the performance of
92 algorithms from 33 international teams at predicting the future trajectory
of 219 individuals at risk of Alzheimer's disease. Challenge participants were
required to make a prediction, for each month of a 5-year future time period,
of three key outcomes: clinical diagnosis, Alzheimer's Disease Assessment Scale
Cognitive Subdomain (ADAS-Cog13), and total volume of the ventricles. The
methods used by challenge participants included multivariate linear regression,
machine learning methods such as support vector machines and deep neural
networks, as well as disease progression models. No single submission was best
at predicting all three outcomes. For clinical diagnosis and ventricle volume
prediction, the best algorithms strongly outperform simple baselines in
predictive ability. However, for ADAS-Cog13 no single submitted prediction
method was significantly better than random guesswork. Two ensemble methods
based on taking the mean and median over all predictions, obtained top scores
on almost all tasks. Better than average performance at diagnosis prediction
was generally associated with the additional inclusion of features from
cerebrospinal fluid (CSF) samples and diffusion tensor imaging (DTI). On the
other hand, better performance at ventricle volume prediction was associated
with inclusion of summary statistics, such as the slope or maxima/minima of
biomarkers. TADPOLE's unique results suggest that current prediction algorithms
provide sufficient accuracy to exploit biomarkers related to clinical diagnosis
and ventricle volume, for cohort refinement in clinical trials for Alzheimer's
disease. However, results call into question the usage of cognitive test scores
for patient selection and as a primary endpoint in clinical trials.


  Motivated by problems arising with the symbolic analysis of steady state
ideals in Chemical Reaction Network Theory, we consider the problem of testing
whether the points in a complex or real variety with non-zero coordinates form
a coset of a multiplicative group. That property corresponds to Shifted
Toricity, a recent generalization of toricity of the corresponding polynomial
ideal. The key idea is to take a geometric view on varieties rather than an
algebraic view on ideals. Recently, corresponding coset tests have been
proposed for complex and for real varieties. The former combine numerous
techniques from commutative algorithmic algebra with Gr\"obner bases as the
central algorithmic tool. The latter are based on interpreted first-order logic
in real closed fields with real quantifier elimination techniques on the
algorithmic side. Here we take a new logic approach to both theories, complex
and real, and beyond. Besides alternative algorithms, our approach provides a
unified view on theories of fields and helps to understand the relevance and
interconnection of the rich existing literature in the area, which has been
focusing on complex numbers, while from a scientific point of view the
(positive) real numbers are clearly the relevant domain in chemical reaction
network theory. We apply prototypical implementations of our new approach to a
set of 129 models from the BioModels repository.


  The temporal growth in the number of deaths in the COVID-19 epidemic is
subexponential. Here we show that a piecewise quadratic law provides an
excellent fit during the thirty days after the first three fatalities on
January 20 and later since the end of March 2020. There is also a brief
intermediate period of exponential growth. During the second quadratic growth
phase, the characteristic time of the growth is about eight times shorter than
in the beginning, which can be understood as the occurrence of separate
hotspots. Quadratic behavior can be motivated by peripheral growth when further
spreading occurs only on the outskirts of an infected region. We also study
numerical solutions of a simple epidemic model, where the spatial extend of the
system is taken into account. To model the delayed onset outside China together
with the early one in China within a single model with minimal assumptions, we
adopt an initial condition of several hotspots, of which one reaches saturation
much earlier than the others. At each site, quadratic growth commences when the
local number of infections has reached a certain saturation level. The total
number of deaths does then indeed follow a piecewise quadratic behavior.


  In the past decades microRNAs (miRNA) have much attracted the attention of
researchers at the interface between life and theoretical sciences for their
involvement in post-transcriptional regulation and related diseases. Thanks to
the always more sophisticated experimental techniques, the role of miRNAs as
"noise processing units" has been further elucidated and two main ways of miRNA
noise-control have emerged by combinations of theoretical and experimental
studies. While on one side miRNA were thought to buffer gene expression noise,
it has recently been suggested that miRNA could also increase the cell-to-cell
variability of their targets. In this Mini Review, we focus on the miRNA role
in noise processing and on the inference of the parameters defined by the
related theoretical modelling.


  Flagella of eukaryotic cells are transient long cylindrical protrusions. The
proteins needed to form and maintain flagella are synthesized in the cell body
and transported to the distal tips. What `rulers' or `timers' a specific type
of cells use to strike a balance between the outward and inward transport of
materials so as to maintain a particular length of its flagella in the steady
state is one of the open questions in cellular self-organization. Even more
curious is how the two flagella of biflagellates, like Chlamydomonas
Reinhardtii, communicate through their base to coordinate their lengths. In
this paper we develop a stochastic model for flagellar length control based on
a time-of-flight (ToF) mechanism. This ToF mechanism decides whether or not
structural proteins are to be loaded onto an intraflagellar transport (IFT)
train just before it begins its motorized journey from the base to the tip of
the flagellum. Because of the ongoing turnover, the structural proteins
released from the flagellar tip are transported back to the cell body also by
IFT trains. We represent the traffic of IFT trains as a totally asymmetric
simple exclusion process (TASEP). The ToF mechanism for each flagellum,
together with the TASEP-based description of the IFT trains, combined with a
scenario of sharing of a common pool of flagellar structural proteins in
biflagellates, can account for all key features of experimentally known
phenomena. These include ciliogenesis, resorption, deflagellation as well as
regeneration after selective amputation of one of the two flagella. We also
show that the experimental observations of Ishikawa and Marshall are consistent
with the ToF mechanism of length control if the effects of the mutual exclusion
of the IFT trains captured by the TASEP are taken into account. Moreover, we
make new predictions on the flagellar length fluctuations and the role of the
common pool.


  Since the end of 2019 an outbreak of a new strain of coronavirus, called
2019--nCoV, is reported from China and later other parts of the world. Since
January 21, WHO reports daily data on confirmed cases and deaths from both
China and other countries. In this work we present some discrete and continuous
models to discribe the disease dynamics in China and estimate the needed
epidemiological parameters. Good agreement with the current dynamics has be
found for both a discrete transmission model and a slightly modified SIR-model.


  The recent outbreak of a novel coronavirus (2019-nCoV) has quickly evolved
into a global health crisis. The transmission potential of 2019-nCoV has been
modelled and studied in several recent research works. The key factors such as
the basic reproductive number, $R_{0}$, of the virus have been identified by
fitting contagious disease spreading models to aggregated data. The data
include the reported cases both within China and in closely connected cities
over the world. In this paper, we study the transmission potential of 2019-nCoV
from the perspective of the robustness of the statistical estimation, in light
of varying data quality and timeliness in the initial stage of the outbreak.
Sample consensus algorithm has been adopted to improve model fitting when
outliers are present. The robust estimation enables us to identify two clusters
of transmission models, both are of substantial concern, one with
$R_0:8\sim14$, comparable to that of measles and the other dictates a large
initial infected group.


  While bacterial chromosomes were long thought to be amorphous, recent
experiments reveal pronounced organizational features. However, the extent of
bacterial chromosome organization remains unclear. Here, we develop a fully
data-driven maximum entropy approach to extract the distribution of single-cell
chromosome conformations from experimental normalized Hi-C data. We apply this
inference to the model organism Caulobacter crescentus. On small genomic scales
of 104-105 basepairs, our model reveals a pattern of local chromosome
extensions that correlates with transcriptional and DNA loop extrusion
activity. On larger genomic scales, we find that chromosome structure is
predominantly present along the long cell axis: chromosomal loci not only have
well-defined axial positions, they also exhibit long-ranged correlations due
interacting large emergent genomic clusters, termed Super Domains. Finally, our
model reveals information contained in chromosome structure that can guide
cellular processes. Our approach can be generalized to other species, providing
a principled way of analyzing spatial chromosome organization.


  The basic reproductive number -- $R_0$ -- is one of the most common and most
commonly misapplied numbers in public health. Although often used to compare
outbreaks and forecast pandemic risk, this single number belies the complexity
that two different pathogens can exhibit, even when they have the same $R_0$.
Here, we show how to predict outbreak size using estimates of the distribution
of secondary infections, leveraging both its average $R_0$ and the underlying
heterogeneity. To do so, we reformulate and extend a classic result from random
network theory that relies on contact tracing data to simultaneously determine
the first moment ($R_0$) and the higher moments (representing the
heterogeneity) in the distribution of secondary infections. Further, we show
the different ways in which this framework can be implemented in the
data-scarce reality of emerging pathogens. Lastly, we demonstrate that without
data on the heterogeneity in secondary infections for emerging infectious
diseases like COVID-19, the uncertainty in outbreak size ranges dramatically.
Taken together, our work highlights the critical need for contact tracing
during emerging infectious disease outbreaks and the need to look beyond $R_0$
when predicting epidemic size.


  The metabolic wiring of patient cells is altered drastically in many
diseases, including cancer. Understanding the nature of such changes may pave
the way for new therapeutic opportunities, as well as the development of
personalized treatment strategies for patients. In this paper, we propose an
algorithm, Metabolitics, which allows systems-level analysis of changes in the
biochemical network of cells in disease states. It enables the study of a
disease at both reaction- and pathway-level granularities for detailed and
summarized view of disease etiology. Metabolitics employs flux variability
analysis with an objective function which is dynamically built based on the
biofluid metabolomics measurements in a personalized manner. Moreover,
Metabolitics builds classification models for each disease to diagnose patients
and predict their subgroups based on the computed metabolic network changes. We
demonstrate the use of Metabolitics on three distinct diseases, namely, breast
cancer, Crohn's disease, and colorectal cancer. Our results show that the
constructed supervised learning models successfully diagnose patients by
f1-score of over 90% on the average. Besides, in addition to the confirmation
of previously reported breast cancer associated pathways, we discovered that
Butanoate metabolism experiences significantly decreased activity, while
Arginine and Proline Metabolism is subject to significant increase in activity,
which have not been reported before. Metabolitics is made available as a Python
package in pypi.


  We show how highly-diverse ecological communities may display persistent
abundance fluctuations, when interacting through resource competition and
subjected to migration from a species pool. This turns out to be closely
related to the ratio of realized species diversity to the number of resources.
This ratio is set by competition, through the balance between species being
pushed out and invading. When this ratio is smaller than one, dynamics will
reach stable equilibria. When this ratio is larger than one, fixed-points are
either unstable or marginally stable, as expected by the competitive exclusion
principle. If they are unstable, the system is repelled from fixed points, and
abundances forever fluctuate. While marginally-stable fixed points are in
principle allowed and predicted by some models, they become structurally
unstable at high diversity. This means that even small changes to the model,
such as non-linearities in how resources combine to generate species' growth,
will result in persistent abundance fluctuations.


  We have used Langevin dynamics to simulate the forced translocation of linked
polymer rings through a narrow pore. For fixed size (i.e. fixed number of
monomers) the translocation time depends on the link type and on whether the
rings are knotted or unknotted. For links with two unknotted rings, the
crossings between the rings can slow down the translocation and are responsible
for a delay as the crossings pass through the pore. The results fall on a set
of relatively smooth curves for different link families with the translocation
time not always increasing with crossings number within the same family. When
one ring is knotted the results depend on whether the link is prime or
composite and, for the composite case, they depend on whether the knotted or
unknotted ring enters the pore first. We find a similar situation for
3-component links where the results depend on whether the link is prime or
composite. These results contribute to our understanding of how the
entanglement complexity between filaments impacts their translocation dynamics
and should be useful for extending nanopore-sensing techniques to probe the
topological properties of these systems.


  At the moment of writing (12 February, 2020), the future evolution of the
2019-nCoV virus is unclear. Predictions of the further course of the epidemic
are decisive to deploy targeted disease control measures. We consider a
network-based model to describe the 2019-nCoV epidemic in the Hubei province.
The network is composed of the cities in Hubei and their interactions (e.g.,
traffic flow). However, the precise interactions between cities is unknown and
must be inferred from observing the epidemic. We propose a network-based method
to predict the future prevalence of the 2019-nCoV virus in every city. Our
results indicate that network-based modelling is beneficial for an accurate
forecast of the epidemic outbreak.


  Contemporary Niche Theory is a useful framework for understanding how
organisms interact with each other and with their shared environment. Its
graphical representation, popularized by Tilman's Resource Ratio Hypothesis,
facilitates the analysis of the equilibrium structure of complex dynamical
models including species coexistence. This theory has been applied primarily to
resource competition since its early beginnings. Here, we integrate mutualism
into niche theory by expanding Tilman's graphical representation to the
analysis of consumer-resource dynamics of plant-pollinator networks. We
graphically explain the qualitative phenomena previously found by numerical
simulations, including the effects on community dynamics of nestedness,
adaptive foraging, and pollinator invasions. Our graphical approach promotes
the unification of niche and network theories, and deepens the synthesis of
different types of interactions within a consumer-resource framework.


  A metric phylogenetic tree relating a collection of taxa induces weighted
rooted triples and weighted quartets for all subsets of three and four taxa,
respectively. New intertaxon distances are defined that can be calculated from
these weights, and shown to exactly fit the same tree topology, but with edge
weights rescaled by certain factors dependent on the associated split size.
These distances are analogs for metric trees of similar ones recently
introduced for topological trees that are based on induced unweighted rooted
triples and quartets. The distances introduced here lead to new statistically
consistent methods of inferring a metric species tree from a collection of
topological gene trees generated under the multispecies coalescent model of
incomplete lineage sorting. Simulations provide insight into their potential.


  The outbreak of novel coronavirus (2019-nCoV) in China has caused a viral
epidemic affecting tens of thousands of persons. Though the danger of this
epidemic is evident, the statistical analysis of data offered in this paper
indicates that the increase in new cases in China will stabilize in the coming
days or weeks. Our forecast could serve to evaluate risks and control the
evolution of this disease.


  Since the SARS outbreak in 2003, a lot of predictive epidemiological models
have been proposed. At the end of 2019, a novel coronavirus, termed as
2019-nCoV, has broken out and is propagating in China and the world. Here we
propose a multi-model ordinary differential equation set neural network
(MMODEs-NN) and model-free methods to predict the interprovincial transmissions
in mainland China, especially those from Hubei Province. Compared with the
previously proposed epidemiological models, the proposed network can simulate
the transportations with the ODEs activation method, while the model-free
methods based on the sigmoid function, Gaussian function, and Poisson
distribution are linear and fast to generate reasonable predictions. According
to the numerical experiments and the realities, the special policies for
controlling the disease are successful in some provinces, and the transmission
of the epidemic, whose outbreak time is close to the beginning of China Spring
Festival travel rush, is more likely to decelerate before February 18 and to
end before April 2020. The proposed mathematical and artificial intelligence
methods can give consistent and reasonable predictions of the 2019-nCoV ending.
We anticipate our work to be a starting point for comprehensive prediction
researches of the 2019-nCoV.


  In a case-control study aimed at locating autosomal disease variants for a
disease of interest, association between markers and the disease status is
often tested by comparing the marker minor allele frequencies (MAFs) between
cases and controls. For most common allele-based tests the statistical power is
highly dependent on the actual values of these MAFs, where associated markers
with low MAFs have less power to be detected compared to associated markers
with high MAFs. Therefore, the popular strategy of selecting markers for
follow-up studies based primarily on their p-values is likely to preferentially
select markers with high MAFs. We propose a new test which does not favor
markers with high MAFs and improves the power for markers with low to moderate
MAFs without sacrificing performance for markers with high MAFs and is
therefore superior to most existing tests in this regard. An explicit formula
for the asymptotic power function of the proposed test is derived
theoretically, which allows for fast and easy computation of the corresponding
p-values. The performance of the proposed test is compared with several
existing tests both in the asymptotic and the finite sample size settings.


  The initial cluster of severe pneumonia cases that triggered the 2019-nCoV
epidemic was identified in Wuhan, China in December 2019. While early cases of
the disease were linked to a wet market, human-to-human transmission has driven
the rapid spread of the virus throughout China. The ongoing outbreak presents a
challenge for modelers, as limited data are available on the early growth
trajectory, and the epidemiological characteristics of the novel coronavirus
are yet to be fully elucidated. We provide timely short-term forecasts of the
cumulative number of confirmed reported cases in Hubei province, the epicenter
of the epidemic, and for the overall trajectory in China, excluding the
province of Hubei. We collect daily reported cumulative case data for the
2019-nCoV outbreak for each Chinese province from the National Health
Commission of China. Here, we provide 5, 10, and 15 day forecasts for five
consecutive days, February 5th through February 9th, with quantified
uncertainty based on a generalized logistic growth model, the Richards growth
model, and a sub-epidemic wave model. Our most recent forecasts reported here
based on data up until February 9, 2020, largely agree across the three models
presented and suggest an average range of 7,409-7,496 additional cases in Hubei
and 1,128-1,929 additional cases in other provinces within the next five days.
Models also predict an average total cumulative case count between 37,415 -
38,028 in Hubei and 11,588 - 13,499 in other provinces by February 24, 2020.
Mean estimates and uncertainty bounds for both Hubei and other provinces have
remained relatively stable in the last three reporting dates (February 7th -
9th). Our forecasts suggest that the containment strategies implemented in
China are successfully reducing transmission and that the epidemic growth has
slowed in recent days.


  Global transport and communication networks enable information, ideas and
infectious diseases now to spread at speeds far beyond what has historically
been possible. To effectively monitor, design, or intervene in such
epidemic-like processes, there is a need to predict the speed of a particular
contagion in a particular network, and to distinguish between nodes that are
more likely to become infected sooner or later during an outbreak. Here, we
study these quantities using a message-passing approach to derive simple and
effective predictions which are validated against epidemic simulations on a
variety of real-world networks with good agreement. In addition to
individualized predictions for different nodes, we find an overall sudden
transition from low density to almost full network saturation as the contagion
develops in time. Our theory is developed and explained in the setting of
simple contagions on tree-like networks, but we are also able to show how the
method extends remarkably well to complex contagions and highly clustered
networks.


  Cooperation is the backbone of modern human societies, making it a priority
to understand how successful cooperation-sustaining mechanisms operate. Cyclic
dominance, a non-transitive setup comprising at least three strategies wherein
the first strategy overrules the second which overrules the third which, in
turn, overrules the first strategy, is known to maintain bio-diversity, drive
competition between bacterial strains, and preserve cooperation in social
dilemmas. Here, we present a novel route to cyclic dominance in voluntary
social dilemmas by adding to the traditional mix of cooperators, defectors, and
loners, a fourth player type, risk-averse hedgers, who enact tit-for-tat upon
paying a hedging cost to avoid being exploited. When this cost is sufficiently
small, cooperators, defectors, and hedgers enter a loop of cyclic dominance
that preserves cooperation even under the most adverse conditions. In contrast,
when the hedging cost is large, hedgers disappear, consequently reverting to
the traditional interplay of cooperators, defectors, and loners. In the interim
region of hedging costs, complex evolutionary dynamics ensues, prompting
transitions between states with two, three, or four competing strategies. Our
results thus reveal that voluntary participation is but one pathway to
sustained cooperation via cyclic dominance.


  Pairwise sequence alignment is one of the most computationally intensive
kernels in genomic data analysis, accounting for more than 90% of the runtime
for key bioinformatics applications. This method is particularly expensive for
third-generation sequences due to the high computational cost of analyzing
sequences of length between 1Kb and 1Mb. Given the quadratic overhead of exact
pairwise algorithms for long alignments, the community primarily relies on
approximate algorithms that search only for high-quality alignments and stop
early when one is not found. In this work, we present the first GPU
optimization of the popular X-drop alignment algorithm, that we named LOGAN.
Results show that our high-performance multi-GPU implementation achieves up to
181.6 GCUPS and speed-ups up to 6.6x and 30.7x using 1 and 6 NVIDIA Tesla V100,
respectively, over the state-of-the-art software running on two IBM Power9
processors using 168 CPU threads, with equivalent accuracy. We also demonstrate
a 2.3x LOGAN speed-up versus ksw2, a state-of-art vectorized algorithm for
sequence alignment implemented in minimap2, a long-read mapping software. To
highlight the impact of our work on a real-world application, we couple LOGAN
with a many-to-many long-read alignment software called BELLA, and demonstrate
that our implementation improves the overall BELLA runtime by up to 10.6x.
Finally, we adapt the Roofline model for LOGAN and demonstrate that our
implementation is near-optimal on the NVIDIA Tesla V100s.


  A fully atomistic modelling of biological macromolecules at relevant length-
and time-scales is often cumbersome or not even desirable, both in terms of
computational effort required and it a posteriori analysis. This difficulty can
be overcome with the use of multi-resolution models, in which different regions
of the same system are concurrently described at different levels of detail. In
enzymes, computationally expensive atomistic detail is crucial in the modelling
of the active site in order to capture e.g. the chemically subtle process of
ligand binding. In contrast, important yet more collective properties of the
remainder of the protein can be reproduced with a coarser description. In the
present work, we demonstrate the effectiveness of this approach through the
calculation of the binding free energy of hen egg white lysozyme (HEWL) with
the inhibitor di-N-acetylchitotriose. Particular attention is posed to the
impact of the mapping, i.e. the selection of atomistic and coarse-grained
residues, on the binding free energy. It is shown that, in spite of small
variations of the binding free energy with respect to the active site
resolution, the separate contributions coming from different energetic terms
(such as electrostatic and van der Waals interactions) manifest a stronger
dependence on the mapping, thus pointing to the existence of an optimal level
of intermediate resolution.


  In December 2019, a novel coronavirus (SARS-CoV-2) has been identified to
cause acute respiratory disease in humans. An outbreak of this disease has been
reported in mainland China with the city of Wuhan as the recognized epicenter.
The disease has also been exported to other countries, including the
Philippines, but the level of spread is still under control (as of 08 February
2020). To describe and predict the dynamics of the disease, several preliminary
mathematical models are formulated by various international study groups. Here,
the insights that can be drawn from these models are discussed, especially as
inputs for designing strategies to control the epidemics. Proposed model-based
strategies on how to prevent the spread of the disease in local setting, such
as during large social gatherings, are also presented. The model shows that the
exposure time is a significant factor in spreading the disease. With a basic
reproduction number equal to 2, and 14-day infectious period, an infected
person staying more than 9 hours in the event could infect other people.
Assuming the exposure time is 18 hours, the model recommends that attendees of
the social gathering should have a protection with more than 70 percent
effectiveness.


  The question of whether a population will persist or go extinct is of key
interest throughout ecology and biology. Various mathematical techniques allow
us to generate knowledge regarding individual behaviour, which can be analysed
to obtain predictions about the ultimate survival or extinction of the
population. A common model employed to describe population dynamics is the
lattice-based random walk model with crowding (exclusion). This model can
incorporate behaviour such as birth, death and movement, while including
natural phenomena such as finite size effects. Performing sufficiently many
realisations of the random walk model to extract representative population
behaviour is computationally intensive. Therefore, continuum approximations of
random walk models are routinely employed. However, standard continuum
approximations are notoriously incapable of making accurate predictions about
population extinction. Here, we develop a new continuum approximation, the
state space diffusion approximation, which explicitly accounts for population
extinction. Predictions from our approximation faithfully capture the behaviour
in the random walk model, and provides additional information compared to
standard approximations. We examine the influence of the number of lattice
sites and initial number of individuals on the long-term population behaviour,
and demonstrate the reduction in computation time between the random walk model
and our approximation.


  Data taking values on discrete sample spaces are the embodiment of modern
biological research. "Omics" experiments produce millions of symbolic outcomes
in the form of reads (i.e., DNA sequences of a few dozens to a few hundred
nucleotides). Unfortunately, these intrinsically non-numerical datasets are
often highly contaminated, and the possible sources of contamination are
usually poorly characterized. This contrasts with numerical datasets where
Gaussian-type noise is often well-justified. To overcome this hurdle, we
introduce the notion of latent weight, which measures the largest expected
fraction of samples from a contaminated probabilistic source that conform to a
model in a well-structured class of desired models. We examine various
properties of latent weights, which we specialize to the class of exchangeable
probability distributions. As proof of concept, we analyze DNA methylation data
from the 22 human autosome pairs. Contrary to what it is usually assumed, we
provide strong evidence that highly specific methylation patterns are
overrepresented at some genomic locations when contamination is taken into
account.


  The knowledge of potentially druggable binding sites on proteins is an
important preliminary step towards the discovery of novel drugs. The
computational prediction of such areas can be boosted by following the recent
major advances in the deep learning field and by exploiting the increasing
availability of proper data. In this paper, a novel computational method for
the prediction of potential binding sites is proposed, called DeepSurf.
DeepSurf combines a surface-based representation, where a number of 3D
voxelized grids are placed on the protein's surface, with state-of-the-art deep
learning architectures. After being trained on the large database of scPDB,
DeepSurf demonstrates superior results on three diverse testing datasets, by
surpassing all its main deep learning-based competitors, while attaining
competitive performance to a set of traditional non-data-driven approaches.


  Patterned vegetation is a characteristic feature of many dryland ecosystems.
While plant densities on the ecosystem-wide scale are typically low, a spatial
self-organisation principle leads to the occurrence of alternating patches of
high biomass and patches of bare soil. Nevertheless, intraspecific competition
dynamics other than competition for water over long spatial scales are commonly
ignored in mathematical models for vegetation patterns. In this chapter, I
address the impact of local intraspecific competition on a modelling framework
for banded vegetation patterns. Firstly, I show that in the context of a
single-species model, neglecting local intraspecific competition leads to an
overestimation of a patterned ecosystem's resilience to increases in aridity.
Secondly, in the context of a multispecies model, I argue that local
intraspecific competition is a key element in the successful capture of species
coexistence in model solutions representing a vegetation pattern. For both
models, a detailed bifurcation analysis is presented to analyse the onset,
existence and stability of patterns. Besides the strengths of local
intraspecific competition, also the the difference between two species has a
significant impact on the bifurcation structure, providing crucial insights
into the complex ecosystem dynamics. Predictions on future ecosystem dynamics
presented in this chapter, especially on pattern onset and pattern stability,
can aid the development of conservation programs.


  Querying new information from knowledge sources, in general, and published
literature, in particular, aims to provide precise and quick answers to
questions raised about a system under study. In this paper, we present
ACCORDION (Automated Clustering Conditional On Relating Data of Interactions tO
a Network), a novel tool and a methodology to enable efficient answering of
biological questions by automatically assembling new, or expanding existing
models using published literature. Our approach integrates information
extraction and clustering with simulation and formal analysis to allow for an
automated iterative process that includes assembling, testing and selecting the
most relevant models, given a set of desired system properties. We applied our
methodology to a model of the circuitry that con-trols T cell differentiation.
To evaluate our approach, we compare the model that we obtained, using our
automated model extension approach, with the previously published manually
extended T cell differentiation model. Besides demonstrating automated and
rapid reconstruction of a model that was previously built manually, ACCORDION
can assemble multiple models that satisfy desired properties. As such, it
replaces large number of tedious or even imprac-tical manual experiments and
guides alternative hypotheses and interventions in biological systems.


  By using the public data from Jan. 20 to Feb. 11, 2020, we perform
data-driven analysis and forecasting on the COVID-19 epidemic in mainland
China, especially Hubei province. Our results show that the turning points of
the daily infections are predicted to be Feb. 6 and Feb. 1, 2020, for Hubei and
China other than Hubei, respectively. The epidemic in China is predicted to end
up after Mar. 10, 2020, and the number of the total infections are predicted to
be 51600. The data trends reveal that quick and active strategies taken by
China to reduce human exposure have already had a good impact on the control of
the epidemic.