  Sequence alignment has had an enormous impact on our understanding of
biology, evolution, and disease. The alignment of biological {\em networks}
holds similar promise. Biological networks generally model interactions between
biomolecules such as proteins, genes, metabolites, or mRNAs. There is strong
evidence that the network topology -- the "structure" of the network -- is
correlated with the functions performed, so that network topology can be used
to help predict or understand function. However, unlike sequence comparison and
alignment -- which is an essentially solved problem -- network comparison and
alignment is an NP-complete problem for which heuristic algorithms must be
used.
  Here we introduce SANA, the {\it Simulated Annealing Network Aligner}. SANA
is one of many algorithms proposed for the arena of biological network
alignment. In the context of global network alignment, SANA stands out for its
speed, memory efficiency, ease-of-use, and flexibility in the arena of
producing alignments between 2 or more networks. SANA produces better
alignments in minutes on a laptop than most other algorithms can produce in
hours or days of CPU time on large server-class machines. We walk the user
through how to use SANA for several types of biomolecular networks.
  Availability: https://github.com/waynebhayes/SANA


  The end-Permian mass extinction is the most severe known from the fossil
record. The most likely cause is massive volcanic activity associated with the
formation of the Permo-Triassic Siberian flood basalts. A proposed mechanism
for extinction due to this volcanic activity is depletion of stratospheric
ozone, leading to increased penetration of biologically damaging Solar
ultraviolet-B (UVB) radiation to Earth's surface. Previous work has modeled the
atmospheric chemistry effects of volcanic emission at the end-Permian. Here we
use those results as input for detailed radiative transfer simulations to
investigate changes in surface-level Solar irradiance in the ultraviolet-B,
ultraviolet-A and photosynthetically available (visible light) wave bands. We
then evaluate the potential biological effects using biological weighting
functions. In addition to changes in ozone column density we also include
gaseous sulfur dioxide (SO2) and sulfate aerosols. Ours is the first such study
to include these factors and we find they have a significant impact on
transmission of Solar radiation through the atmosphere. Inclusion of SO2 and
aerosols greatly reduces the transmission of radiation across the ultraviolet
and visible wavelengths, with subsequent reduction in biological impacts by
UVB. We conclude that claims of a UVB mechanism for this extinction are likely
overstated.


  Characterization of the differences between biological and random networks
can reveal the design principles that enable the robust realization of crucial
biological functions including the establishment of different cell types.
Previous studies, focusing on identifying topological features that are present
in biological networks but not in random networks, have, however, provided few
functional insights. We use a Boolean modeling framework and ideas from spin
glass literature to identify functional differences between five real
biological networks and random networks with similar topological features. We
show that minimal frustration is a fundamental property that allows biological
networks to robustly establish cell types and regulate cell fate choice, and
this property can emerge in complex networks via Darwinian evolution. The study
also provides clues regarding how the regulation of cell fate choice can go
awry in a disease like cancer and lead to the emergence of aberrant cell types.


  This paper investigates a dynamical predator-prey interaction model that
incorporates: (a) hunting cooperation among predators; (b) Allee effect in
prey. We show all possible boundary and interior solutions. In order to analyze
the stability of the solution, we make use of the Jacobian matrix and the
resultant characteristic polynomial. Particularly, the sign of the eigenvalue
is used to determine the stability of a solution. We then provide proof for
stability of the interior solution. Finally, we verify our results numerically
in MATLAB by plotting: (1) predator-prey intersection graphs; (2) prey-predator
vs hunting cooperation graphs; (3) initial condition trajectory for equilibrium
solution. It is interesting to notice that hunting cooperation can switch the
stability of coexistence equilibrium solutions. Through numerical simulations,
it was verified that increasing the hunting cooperation could lead to the
extinction of both prey and predator population for alpha greater than 0.96,
given our choice of parameters.


  Drug combinations are frequently used for the treatment of cancer patients in
order to increase efficacy, decrease adverse side effects, or overcome drug
resistance. Given the enormous number of drug combinations, it is cost- and
time-consuming to screen all possible drug pairs experimentally. Currently, it
has not been fully explored to integrate multiple networks to predict
synergistic drug combinations using recently developed deep learning
technologies. In this study, we proposed a Graph Convolutional Network (GCN)
model to predict synergistic drug combinations in particular cancer cell lines.
Specifically, the GCN method used a convolutional neural network model to do
heterogeneous graph embedding, and thus solved a link prediction task. The
graph in this study was a multimodal graph, which was constructed by
integrating the drug-drug combination, drug-protein interaction, and
protein-protein interaction networks. We found that the GCN model was able to
correctly predict cell line-specific synergistic drug combinations from a large
heterogonous network. The majority (30) of the 39 cell line-specific models
show an area under the receiver operational characteristic curve (AUC) larger
than 0.80, resulting in a mean AUC of 0.84. Moreover, we conducted an in-depth
literature survey to investigate the top predicted drug combinations in
specific cancer cell lines and found that many of them have been found to show
synergistic antitumor activity against the same or other cancers in vitro or in
vivo. Taken together, the results indicate that our study provides a promising
way to better predict and optimize synergistic drug pairs in silico.


  Social categorizations divide people into "us" and "them," often along
continuous attributes such as political ideology or skin color. This division
results in both positive consequences, such as a sense of community, and
negative ones, such as group conflict. Further, individuals in the middle of
the spectrum can fall through the cracks of this categorization process and are
seen as out-group by individuals on either side of the spectrum, becoming
inbetweeners. Here, we propose a quantitative, dynamical-system model that
studies the joint influence of cognitive and social processes. We model where
two social groups draw the boundaries between "us" and "them" on a continuous
attribute. Our model predicts that both groups tend to draw a more restrictive
boundary than the middle of the spectrum. As a result, each group sees the
individuals in the middle of the attribute space as an out-group. We test this
prediction using U.S. political survey data on how political independents are
perceived by registered party members as well as existing experiments on the
perception of racially ambiguous faces, and find support.


  The zinc finger structure where a Zn2+ ion binds to 4 cysteine or histidine
amino acids in a tetrahedral structure is very common motif of nucleic acid
binding proteins. The corresponding interaction model is present in 3% of the
genes of human genome. As a result, zinc finger has been shown to be extremely
useful in various therapeutic and research capacities, as well as in
biotechnology. In stable configuration, the cysteine amino acids are
deprotonated and become negatively charged. This means the Zn2+ ion is
overscreened by 4 cysteine charges (overcharged). It is question of whether
this overcharged configuration is also stable when such negatively charged zinc
finger binds to negatively charged DNA molecule. Using all atom molecular
dynamics simulation up to microsecond range of an androgen receptor protein
dimer, we investigate how the deprotonated state of cysteine influences its
structure, dynamics, and function in binding o DNA molecules. Our results show
that the deprotonated state of cysteine residues are essential for mechanical
stabilization of the functional, folded conformation. Not only this state
stabilizes the protein structure, it also stabilizes the protein-DNA binding
complex. The differences in structural and energetic properties of the two
(sequence-identical) monomers are also investigated showing the strong
influence of DNA on the structure of zinc fingers upon complexation. Our result
has potential impact on better molecular understanding of one of the most
common classes of zinc fingers


  In the prospect of ecology and biology, studying insect-plant predation will
considerably contribute to pest control, benefit agriculture and afforestation,
and also help people to better understand insect-plant co-evolution. Therefore,
we are motivated to do two work in this study. The first part is to cluster the
insect-plant predation, in such manner, unobserved predation could be
estimated. The second part is to explore the connection between predation and
bio-taxonomy, and we find insects get more divergence than plants during the
insect-plant co-evolution.


  The savanna biome is characterised by a continuous vegetation cover,
comprised of herbaceous and woody plants. The coexistence of species in arid
savannas, where water availability is the main limiting resource for plant
growth, provides an apparent contradiction to the classical principle of
competitive exclusion. Previous theoretical work using nonspatial models has
focussed on the development of an understanding of coexistence mechanisms
through the consideration of resource niche separation and ecosystem
disturbances. In this paper, we propose that a spatial self-organisation
principle, caused by a positive feedback between local vegetation growth and
water redistribution, is sufficient for species coexistence in savanna
ecosystems. We propose a spatiotemporal ecohydrological model of partial
differential equations, based on the Klausmeier reaction-advection-diffusion
model for vegetation patterns, to investigate the effects of spatial
interactions on species coexistence on sloped terrain. Our results suggest that
species coexistence is a possible model outcome, if a balance is kept between
the species' average fitness (a measure of a species' competitive abilities in
a spatially uniform setting) and their colonisation abilities. Spatial
heterogeneities in resource availability are utilised by the superior coloniser
(grasses), before it is outcompeted by the species of higher average fitness
(trees). A stability analysis of the spatially nonuniform coexistence solutions
further suggests that grasses act as ecosystem engineers and facilitate the
formation of a continuous tree cover for precipitation levels unable to support
a uniform tree density in the absence of a grass species.


  Patterns of vegetation are a characteristic feature of many semi-arid
regions. The limiting resource in these ecosystems is water, which is added to
the system through short and intense rainfall events that cause a pulse of
biological processes such as plant growth and seed dispersal. We propose an
impulsive model based on the Klausmeier reaction-advection-diffusion system,
analytically investigate the effects of rainfall intermittency on the onset of
patterns, and augment our results by numerical simulations of model extensions.
Our investigation focuses on the parameter region in which a transition between
uniform and patterned vegetation occurs. Results show that decay-type processes
associated with a low frequency of precipitation pulses inhibit the onset of
patterns and that under intermittent rainfall regimes, a spatially uniform
solution is sustained at lower total precipitation volumes than under
continuous rainfall, if plant species are unable to efficiently use low soil
moisture levels. Unlike in the classical setting of a reaction-diffusion model,
patterns are not caused by a diffusion-driven instability but by a combination
of sufficiently long periods of droughts between precipitation pulses and water
diffusion. Our results further indicate that the introduction of pulse-type
seed dispersal weakens the effects of changes to width and shape of the plant
dispersal kernel on the onset of patterns.


  Vegetation patterns are a characteristic feature of semi-deserts occurring on
all continents except Antarctica. In some semi-arid regions, the climate is
characterised by seasonality, which yields a synchronisation of seed dispersal
with the dry season or the beginning of the wet season. We reformulate the
Klausmeier model, a reaction-advection-diffusion system that describes the
plant-water dynamics in semi-arid environments, as an integrodifference model
to account for the temporal separation of plant growth processes during the wet
season and seed dispersal processes during the dry season. The model further
accounts for nonlocal processes involved in the dispersal of seeds. Our
analysis focusses on the onset of spatial patterns. The Klausmeier partial
differential equations (PDE) model is linked to the integrodifference model in
an appropriate limit, which yields a control parameter for the temporal
separation of seed dispersal events. We find that the conditions for pattern
onset in the integrodifference model are equivalent to those for the continuous
PDE model and hence independent of the time between seed dispersal events. We
thus conclude that in the context of seed dispersal, a PDE model provides a
sufficiently accurate description, even if the environment is seasonal. This
emphasises the validity of results that have previously been obtained for the
PDE model. Further, we numerically investigate the effects of changes to seed
dispersal behaviour on the onset of patterns. We find that long-range seed
dispersal inhibits the formation of spatial patterns and that the seed
dispersal kernel's decay at infinity is a significant regulator of patterning.


  Vegetation patterns are a ubiquitous feature of water-deprived ecosystems.
Despite the competition for the same limiting resource, coexistence of several
plant species is commonly observed. We propose a two-species reaction-diffusion
model based on the single-species Klausmeier model, to analytically investigate
the existence of states in which both species coexist. Ecologically, the study
finds that coexistence is supported if there is a small difference in the plant
species' average fitness, measured by the ratio of a species' capabilities to
convert water into new biomass to its mortality rate. Mathematically,
coexistence is not a stable solution of the system, but both spatially uniform
and patterned coexistence states occur as metastable states. In this context, a
metastable solution in which both species coexist corresponds to a long
transient (exceeding $10^3$ years in dimensional parameters) to a stable
one-species state. This behaviour is characterised by the small size of a
positive eigenvalue which has the same order of magnitude as the average
fitness difference between the two species. Two mechanisms causing the
occurrence of metastable solutions are established: a spatially uniform
unstable equilibrium and a stable one-species pattern which is unstable to the
introduction of a competitor. We further discuss effects of asymmetric
interspecific competition (e.g. shading) on the metastability property.


  Vegetation patterns are a characteristic feature of semi-arid regions. On
hillsides these patterns occur as stripes running parallel to the contours. The
Klausmeier model, a coupled reaction-advection-diffusion system, is a
deliberately simple model describing the phenomenon. In this paper, we replace
the diffusion term describing plant dispersal by a more realistic nonlocal
convolution integral to account for the possibility of long-range dispersal of
seeds. Our analysis focuses on the rainfall level at which there is a
transition between uniform vegetation and pattern formation. We obtain results,
valid to leading order in the large parameter comparing the rate of water flow
downhill to the rate of plant dispersal, for a negative exponential dispersal
kernel. Our results indicate that both a wider dispersal of seeds and an
increase in dispersal rate inhibit the formation of patterns. Assuming an
evolutionary trade-off between these two quantities, mathematically motivated
by the limiting behaviour of the convolution term, allows us to make
comparisons to existing results for the original reaction-advection-diffusion
system. These comparisons show that the nonlocal model always predicts a larger
parameter region supporting pattern formation. We then numerically extend the
results to other dispersal kernels, showing that the tendency to form patterns
depends on the type of decay of the kernel.


  In noisy environments such as the cell, many processes involve target sites
that are often hidden or inactive, and thus not always available for reaction
with diffusing entities. To understand reaction kinetics in these situations,
we study the first hitting time statistics of a Brownian particle searching for
a target site that switches stochastically between visible and hidden phases.
At high crypticity, an unexpected rate limited power-law regime emerges for the
first hitting time density, which markedly differs from the classic $t^{-3/2}$
scaling for steady targets. Our problem admits an asymptotic mapping onto a
mixed, or Robin, boundary condition. Similar results are obtained with
non-Markov targets and particles diffusing anomalously.


  Complex I (NADH:ubiquinone oxidoreductase) is a redox-driven proton pump that
powers synthesis of adenosine triphosphate (ATP) and active transport in most
organisms. This gigantic enzyme reduces quinone (Q) to quinol (QH2) in its
hydrophilic domain, and transduces the released free energy into pumping of
protons across its membrane domain, up to ca. 200 {\AA} away from its active
Q-reduction site. Recently resolved molecular structures of complex I from
several species have made it possible for the first time to address the
energetics and dynamics of the complete complex I using multi-scale methods of
computational biochemistry. Here it is described how molecular simulations can
provide important mechanistic insights into the function of the remarkable
pumping machinery in complex I and stimulate new experiments.


  Despite the recent progress in genome sequencing and assembly, many of the
currently available assembled genomes come in a draft form. Such draft genomes
consist of a large number of genomic fragments (scaffolds), whose order and/or
orientation (i.e., strand) in the genome are unknown. There exist various
scaffold assembly methods, which attempt to determine the order and orientation
of scaffolds along the genome chromosomes. Some of these methods (e.g., based
on FISH physical mapping, chromatin conformation capture, etc.) can infer the
order of scaffolds, but not necessarily their orientation. This leads to a
special case of the scaffold orientation problem (i.e., deducing the
orientation of each scaffold) with a known order of the scaffolds.
  We address the problem of orientating ordered scaffolds as an optimization
problem based on given weighted orientations of scaffolds and their pairs
(e.g., coming from pair-end sequencing reads, long reads, or homologous
relations). We formalize this problem using notion of a scaffold graph (i.e., a
graph, where vertices correspond to the assembled contigs or scaffolds and
edges represent connections between them). We prove that this problem is
NP-hard, and present a polynomial-time algorithm for solving its special case,
where orientation of each scaffold is imposed relatively to at most two other
scaffolds. We further develop an FPT algorithm for the general case of the OOS
problem.


  The human-associated microbiome is closely tied to human health and is of
substantial clinical interest. Metagenomics-based tools are emerging for
clinical diagnostics, tracking the spread of diseases, and surveillance of
potential pathogens. In some cases, these tools are overcoming limitations of
traditional clinical approaches. Metagenomics has limitations barring the tools
from clinical validation. Once these hurdles are overcome, clinical
metagenomics will inform doctors of the best, targeted treatment for their
patients and provide early detection of disease. Here we present an overview of
metagenomics methods with a discussion of computational challenges and
limitations.


  Finding tumour genetic markers is essential to biomedicine due to their
relevance for cancer detection and therapy development. In this paper, we
explore a recently released dataset of chromosome rearrangements in 2,586
cancer patients, where different sorts of alterations have been detected. Using
a Random Forest classifier, we evaluate the relevance of several features (some
directly available in the original data, some engineered by us) related to
chromosome rearrangements. This evaluation results in a set of potential tumour
genetic markers, some of which are validated in the bibliography, while others
are potentially novel.


  Predicting the three-dimensional (3D) functional structures of proteins
remains an important computational milestone in molecular biology to be
achieved. This feat is hinged on a clear understanding of the mechanism which
proteins use to fold into their native structures. Since Levinthal's paradox,
there has been a lot of progress in understanding this mechanism. Most of the
earlier attempts were caught between assigning either hydrophobic interactions
or hydrogen bonding as the dominant folding force. However, a consensus now
seems to be emerging about hydrogen bonding being a stronger force.
Interestingly, a view from chaperone action may further throw some light on the
nature of the folding mechanism. Thus the very mechanisms which prevent protein
aggregation and misfolding, could help us have a better understanding of the
folding mechanism itself.


  Many disease models focus on characterizing the underlying transmission
mechanism but make simple, possibly naive assumptions about how infections are
reported. In this note, we use a simple deterministic
Susceptible-Infected-Removed (SIR) model to compare two common assumptions
about disease incidence reports: individuals can report their infection as soon
as they become infected or as soon as they recover. We show that incorrect
assumptions about the underlying observation processes can bias estimates of
the basic reproduction number and lead to overly narrow confidence intervals.


  Background: Information on influenza-associated mortality in Russia is
limited.
  Methods: Using previously developed methodology (Goldstein et al.,
Epidemiology 2012), we regressed the monthly rates of mortality for respiratory
causes, as well as circulatory causes linearly against the monthly proxies for
the incidence of influenza A/H3N2, A/H1N1 and B (obtained using data from the
Smorodintsev Research Institute of Influenza (RII)), adjusting for the baseline
rates of mortality not associated with influenza circulation and temporal
trends.
  Results: For the 2013/14 through the 2018/19 seasons, influenza circulation
was associated with an average annual 17636 (95% CI (9482,25790)) deaths for
circulatory causes and 4179 (3250,5109) deaths for respiratory causes, with the
largest number of deaths (32298 (18071,46525) for circulatory causes and 6689
(5019,8359) for respiratory causes) estimated during the 2014/15 influenza
season. The biggest contributor to both circulatory and respiratory
influenza-associated deaths was influenza A/H3N2, followed by influenza B and
A/H1N1. Compared to the 2013/14 through the 2015/16 seasons, during the 2016/17
through the 2018/19 seasons (when levels of influenza vaccination were
significantly higher), the number of influenza-associated deaths declined by
about 16.1%, or 3809 annual respiratory and circulatory deaths.
  Conclusions: Influenza circulation is associated with a substantial mortality
burden in Russia, particularly for circulatory deaths, with some reduction
observed following the major increase in vaccination coverage. Those results
support the potential utility of further extending the levels of influenza
vaccination, the use of quadrivalent influenza vaccines, and extra efforts for
protecting individuals with circulatory disease in Russia, including
vaccination and the use of antiviral medications.


  Much traditional statistical modelling assumes that the outcome variables of
interest are independent of each other when conditioned on the explanatory
variables. This assumption is strongly violated in the case of infectious
diseases, particularly in close-contact settings such as households, where each
individual's probability of infection is strongly influenced by whether other
household members experience infection. On the other hand, general multi-type
transmission models of household epidemics quickly become unidentifiable from
data as the number of types increases. This has led to a situation where it is
has not been possible to draw consistent conclusions from household studies of
infectious diseases, for example in the event of an influenza pandemic. Here,
we present a generalised linear modelling framework for binary outcomes in
sub-units that can (i) capture the effects of non-independence arising from a
transmission process and (ii) adjust estimates of disease risk and severity for
differences in study population characteristics. This model allows for
computationally fast estimation, uncertainty quantification, covariate choice
and model selection. In application to real pandemic influenza household data,
we show that it is formally favoured over existing modelling approaches.


  Increasing rates of global trade and travel, as well as changing climatic
patterns, have led to more frequent outbreaks of plant disease epidemics
worldwide. Mathematical modelling is a key tool in predicting where and how
these new threats will spread, as well as in assessing how damaging they might
be. Models can also be used to inform disease management, providing a rational
methodology for comparing the performance of possible control strategies
against one another. For emerging epidemics, in which new pathogens or pathogen
strains are actively spreading into new regions, the spatial component of
spread becomes particularly important, both to make predictions and to optimise
disease control. In this chapter we illustrate how the spatial spread of
emerging plant diseases can be modelled at the landscape scale via spatially
explicit compartmental models. Our particular focus is on the crucial role of
the dispersal kernel-which parameterises the probability of pathogen spread
from an infected host to susceptible hosts at any given distance-in determining
outcomes of epidemics. We add disease management to our model by testing
performance of a simple ''one off'' form of reactive disease control, in which
sites within a particular distance of locations detected to contain infection
are removed in a single round of disease management. We use this simplified
model to show how ostensibly arcane decisions made by the modeller-most notably
whether or not the underpinning disease model allows for stochasticity (i.e.
randomness)-can greatly impact on disease management recommendations. Our
chapter is accompanied by example code in the programming language R available
via an online repository, allowing the reader to run the models we present for
him/herself.


  Ecosystems are governed by dynamic processes such as competition for
resources, reproduction and dispersal. These shape their biodiversity and how
the system responds to change. Current approaches to modelling ecosystems,
especially plants, focus on either describing fine-scale processes for
individual species or broad-scale patterns for limited groups of plant
functional types.
  Digitisation of herbarium and other plant records has unlocked a wealth of
information that can be used to drive models of plant communities and make
predictions for their future under different scenarios of climate change. The
advent of increased computational capacity and fast, high level programming
languages allows for simulation of such landscapes at unprecedented scales.
  Here, we demonstrate a tool for Ecosystem Simulation through Integrated
Species Trait-Environment Modelling (EcoSISTEM), which models plant species
across multiple ecosystem sizes, from patches and small islands to regions and
entire continents. These simulated ecosystems support the ability to generate
many different types of habitat, as well as reproducing different disturbance
scenarios such as climate change, habitat loss and invasion. EcoSISTEM also
reproduces examples of real-world species distributions by integrating plant
occurrence records and global climate reconstructions to simulate plant species
throughout the continent of Africa for the past century.
  EcoSISTEM allows us to flexibly explore the dynamics of tens of thousands of
species interacting across a continent. The code parallelises efficiently
across multiple nodes on high performance computing platforms, and has been
scaled up to run on over 1000 cores. It allows us to study the impact of
changes to climate, resources and habitat and investigate real-life mechanisms
surrounding climate change and biodiversity loss.


  Background: Haplotypes, the ordered lists of single nucleotide variations
that distinguish chromosomal sequences from their homologous pairs, may reveal
an individual's susceptibility to hereditary and complex diseases and affect
how our bodies respond to therapeutic drugs. Reconstructing haplotypes of an
individual from short sequencing reads is an NP-hard problem that becomes even
more challenging in the case of polyploids. While increasing lengths of
sequencing reads and insert sizes {\color{black} helps improve accuracy of
reconstruction}, it also exacerbates computational complexity of the haplotype
assembly task. This has motivated the pursuit of algorithmic frameworks capable
of accurate yet efficient assembly of haplotypes from high-throughput
sequencing data.
  Results: We propose a novel graphical representation of sequencing reads and
pose the haplotype assembly problem as an instance of community detection on a
spatial random graph. To this end, we construct a graph where each read is a
node with an unknown community label associating the read with the haplotype it
samples. Haplotype reconstruction can then be thought of as a two-step
procedure: first, one recovers the community labels on the nodes (i.e., the
reads), and then uses the estimated labels to assemble the haplotypes. Based on
this observation, we propose ComHapDet - a novel assembly algorithm for diploid
and ployploid haplotypes which allows both bialleleic and multi-allelic
variants.
  Conclusions: Performance of the proposed algorithm is benchmarked on
simulated as well as experimental data obtained by sequencing Chromosome $5$ of
tetraploid biallelic \emph{Solanum-Tuberosum} (Potato). The results demonstrate
the efficacy of the proposed method and that it compares favorably with the
existing techniques.


  The stationary sampling distribution of a neutral decoupled Moran or
Wright-Fisher diffusion with neutral mutations is known to first order for a
general rate matrix with small but otherwise unconstrained mutation rates.
Using this distribution as a starting point we derive results for maximum
likelihood estimates of scaled mutation rates from site frequency data under
three model assumptions: a twelve-parameter general rate matrix, a
nine-parameter reversible rate matrix, and a six-parameter strand-symmetric
rate matrix. The site frequency spectrum is assumed to be sampled from a fixed
size population in equilibrium, and to consist of allele frequency data at a
large number of unlinked sites evolving with a common mutation rate matrix
without selective bias. We correct an error in a previous treatment of the same
problem (Burden and Tang, 2017) affecting the estimators for the general and
strand-symmetric rate matrices. The method is applied to a biological dataset
consisting of a site frequency spectrum extracted from short autosomal introns
in a sample of Drosophila melanogaster individuals.


  The coupling between population growth and transport accessibility has been
an elusive problem for more than 60 years now. Due to the lack of theoretical
foundations, most of the studies that considered how the evolution of
transportation networks impacts the population growth are based on regression
analysis in order to identify relevant variables. The recent availability of
large amounts of data allows us to envision the construction of new approaches
for understanding this coupling between transport and population growth. Here,
we use a detailed dataset for about 36000 municipalities in France from 1968
until now. In the case of large urban areas such as Paris, we show that growth
rate statistical variations decay in time and display a trend towards
homogeneization where local aspects are less relevant. We also show that growth
rate differences due to accessibility are very small and can mostly be observed
for cities that experienced very large accessibility variations. This suggests
that the relevant variable for explaining growth rate variations is not the
accessibility but its temporal variation. We propose a model that integrates
the stochastic internal variation of the municipalitie's population and an
inter-urban migration term that we show to be proportional to the accessibility
variation and has a limited time duration. This model provides a simple
theoretical framework that allows to go beyond econometric studies and sheds a
new light on the impact of transportation modes on city growth.


  The majority of cancer treatments end in failure due to Intra-Tumor
Heterogeneity (ITH). ITH in cancer is represented by clonal evolution where
different sub-clones compete with each other for resources under conditions of
Darwinian natural selection. Predicting the growth of these sub-clones within a
tumour is among the key challenges of modern cancer research. Predicting tumor
behavior enables the creation of risk profiles for patients and the
optimisation of their treatment by therapeutically targeting sub-clones more
likely to grow. Current research efforts in this space are focused on
mathematical modelling of population genetics to quantify the selective
advantage of sub-clones, thus enabling predictions of which sub-clones are more
likely to grow. These tumor evolution models are based on assumptions which are
not valid for real-world tumor micro-environment. Furthermore, these models are
often fit on a single instance of a tumor, and hence prediction models cannot
be validated. This paper presents an alternative approach for predicting cancer
evolution using a data-driven machine learning method. Our proposed method is
based on the intuition that if we can capture the true characteristics of
sub-clones within a tumor and represent it in the form of features, a
sophisticated machine learning algorithm can be trained to predict its
behavior. The work presented here provides a novel approach to predicting
cancer evolution, utilizing a data-driver approach. We strongly believe that
the accumulation of data from microbiologists, oncologists and machine learning
researchers could be used to encapsulate the true essence of tumor sub-clones,
and can play a vital role in selecting the best cancer treatments for patients.


  In this paper, we consider a SIRS model with general nonmonotone and
saturated incidence rate and perform stability and bifurcation analysis. We
show that the system has saddle-node bifurcation and displays bistable
behavior. We obtain the critical thresholds that characterize the dynamical
behaviors of the model. We find with surprise that the system always admits a
disease free equilibrium E0 which is always asymptotically stable. Numerical
simulations are carried out to verify our results.


  We propose a new method to estimate plant biodiversity with R{\'e}nyi and Rao
indexes through the so called High Order Singular Value Decomposition (HOSVD)
of tensors. Starting from NASA multispectral images we evaluate biodiversity
and we compare original biodiversity estimates with those realised via the
HOSVD compression methods for big data. Our strategy turns out to be extremely
powerful in terms of storage memory and precision of the outcome. The obtained
results are so promising that we can support the efficiency of our method in
the ecological framework.


  Local ancestry inference (LAI) allows identification of the ancestry of all
chromosomal segments in admixed individuals, and it is a critical step in the
analysis of human genomes with applications from pharmacogenomics and precision
medicine to genome-wide association studies. In recent years, many LAI
techniques have been developed in both industry and academic research. However,
these methods require large training data sets of human genomic sequences from
the ancestries of interest. Such reference data sets are usually limited,
proprietary, protected by privacy restrictions, or otherwise not accessible to
the public. Techniques to generate training samples that resemble real haploid
sequences from ancestries of interest can be useful tools in such scenarios,
since a generalized model can often be shared, but the unique human sample
sequences cannot. In this work we present a class-conditional VAE-GAN to
generate new human genomic sequences that can be used to train local ancestry
inference (LAI) algorithms. We evaluate the quality of our generated data by
comparing the performance of a state-of-the-art LAI method when trained with
generated versus real data.


  Background: Information on influenza-associated mortality in Russia is
limited.
  Methods: Using previously developed methodology (Goldstein et al.,
Epidemiology 2012), we regressed the monthly rates of mortality for respiratory
causes, as well as circulatory causes linearly against the monthly proxies for
the incidence of influenza A/H3N2, A/H1N1 and B (obtained using data from the
Smorodintsev Research Institute of Influenza (RII)), adjusting for the baseline
rates of mortality not associated with influenza circulation and temporal
trends.
  Results: For the 2013/14 through the 2018/19 seasons, influenza circulation
was associated with an average annual 17636 (95% CI (9482,25790)) deaths for
circulatory causes and 4179 (3250,5109) deaths for respiratory causes, with the
largest number of deaths (32298 (18071,46525) for circulatory causes and 6689
(5019,8359) for respiratory causes) estimated during the 2014/15 influenza
season. The biggest contributor to both circulatory and respiratory
influenza-associated deaths was influenza A/H3N2, followed by influenza B and
A/H1N1. Compared to the 2013/14 through the 2015/16 seasons, during the 2016/17
through the 2018/19 seasons (when levels of influenza vaccination were
significantly higher), the volume of influenza-associated mortality declined by
about 16.1%, or 3809 annual respiratory and circulatory deaths.
  Conclusions: Influenza circulation is associated with a substantial mortality
burden in Russia, particularly for circulatory deaths, with some reduction
observed following the major increase in vaccination coverage. Those results
support the potential utility of further extending the levels of influenza
vaccination, the use of quadrivalent influenza vaccines, and extra efforts for
protecting individuals with circulatory disease in Russia, including
vaccination and the use of antiviral medications.


  In this study, we present Flatsomatic - a Variational Auto Encoder (VAE)
optimized to compress somatic mutations that allow for unbiased data
compression whilst maintaining the signal. We compared two different neural
network architectures for the VAE: Multilayer Perceptron (MLP) and
bidirectional LSTM. The somatic profiles we used to train our models consisted
of 8,062 Pan-Cancer patients from The Cancer Genome Atlas and 989 cell lines
from the COSMIC cell line project. The profiles for each patient were
represented by the genomic loci where somatic mutations occurred and, to reduce
sparsity, the locations with a frequency <5 were removed. We enhanced the VAE
performance by changing its evidence lower bound, and devised an F1-score based
loss showing that it helps the VAE learn better than with binary cross-entropy.
We also employed beta-VAE to weight the variational regularisation term in the
loss function and showed the best performance through a preliminary function to
increase the weight of the regularisation term with each epoch. We assessed the
reconstruction ability of the VAE using the micro F1-score metric and showed
that our best performing model was a 2-layer deep MLP VAE. Our analysis also
showed that the size of the latent space did not have a significant effect on
the VAE learning ability. We compared the Flatsomatic embeddings created to a
lower dimension version of the data from principal component analysis, showing
superior performance of Flatsomatic, and performed K-means clustering on both
datasets to draw comparisons to known cancer types of each profile. Finally, we
present results that confirm that the Flatsomatic representations of 64
dimensions maintain the same predictive power as the original 8,298 dimensions
vector, through prediction of drug response.


  According to the features of drug addiction, this paper constructs an
SEIR-based SUC model to describe and predict the spread of drug addiction.
Predictions are that the number of drug addictions will continue to fluctuate
with reduced amplitude and eventually stabilize. To seek the fountainhead of
heroin, we identified the most likely origins of drugs in Philadelphia, PA,
Cuyahoga and Hamilton, OH, Jefferson, KY, Kanawha, WV, and Bedford, VA. Based
on the facts, advised concentration includes the spread of Oxycodone,
Hydrocodone, Heroin, and Buprenorphine. In other words, drug transmission in
the two states of Ohio and Pennsylvania require awareness. According to the
propagation curve predicted by our model, the transfer of KY state is still in
its early stage, while that of VA, WV is in the middle point, and OH, PA in its
latter ones. As a result of this, the number of drug addictions in KY, OH, and
VA is projected to increase in three years. For methodology, with the Principal
component analysis technique, 22 variables in socio-economic data related to
the continuous use of Opioid drugs was filtered, where the 'Relationship' Part
deserves a highlight.
  Based on them, by using the K-means algorithm, 464 counties were categorized
into three baskets. To combat the opioid crisis, a specific action will discuss
in the sensitivity analysis section. After modeling and analytics, innovation
is required to control addicts and advocate anti-drug news campaigns. This part
also verified the effectiveness of model when $d_1<0.2; r_1,r_2,r_3<0.3;
15<\beta_1,\beta_2,\beta_3<25$. In other words, if such boundary exceeded, the
number of drug addictions may rocket and peak in a short period.


  The cornerstone of computational drug design is the calculation of binding
affinity between two biological counterparts, especially a chemical compound,
i.e., a ligand, and a protein. Predicting the strength of protein-ligand
binding with reasonable accuracy is critical for drug discovery. In this paper,
we propose a data-driven framework named DeepAtom to accurately predict the
protein-ligand binding affinity. With 3D Convolutional Neural Network (3D-CNN)
architecture, DeepAtom could automatically extract binding related atomic
interaction patterns from the voxelized complex structure. Compared with the
other CNN based approaches, our light-weight model design effectively improves
the model representational capacity, even with the limited available training
data. With validation experiments on the PDBbind v.2016 benchmark and the
independent Astex Diverse Set, we demonstrate that the less feature engineering
dependent DeepAtom approach consistently outperforms the other state-of-the-art
scoring methods. We also compile and propose a new benchmark dataset to further
improve the model performances. With the new dataset as training input,
DeepAtom achieves Pearson's R=0.83 and RMSE=1.23 pK units on the PDBbind v.2016
core set. The promising results demonstrate that DeepAtom models can be
potentially adopted in computational drug development protocols such as
molecular docking and virtual screening.


  We study the stochastic dynamics of a system of interacting species in a
stochastic environment by means of a continuous-time Markov chain with
transition rates depending on the state of the environment. Models of gene
regulation in systems biology take this form. We characterise the finite-time
distribution of the Markov chain, provide conditions for ergodicity, and
characterise the stationary distribution (when it exists) as a mixture of
Poisson distributions. The mixture measure is uniquely identified as the law of
a fixed point of a stochastic recurrence equation. This recursion is crucial
for statistical computation of moments and other distributional features.


  Most of the formalizations of the Inclusive Fitness and Kin Selection
concepts contain two highly simplifying assumptions (which actually complicate
the problem): first, is that the models are described in terms of "fitness", an
abstract parameter vaguely related to the population growth rate or expected
reproductive success. The second is that they ignore the division between
Donors and Receivers of altruism and the distribution of those roles in the
population. In this paper the Inclusive Fitness and Kin Selection approaches
will be expressed by explicit demographic parameters describing the probability
of death during focal interaction when cooperative trait can be exhibited. In
addition, new framework will take into account the distribution of the roles of
Donor or the Receiver in the population. This description will be used for the
derivation of the population growth model describing the competition between
Cooperative and Noncooperative strategies. The obtained approach will be
sufficient for description of the cases when the roles are independently drawn
during each interaction. However, we can imagine situations when survival and
the change of the role are not correlated (such as food support for the
infected individual, which keeps him alive, but cannot cure). To cover those
cases, new model is extended to the case with explicit dynamics of the role
distributions among carriers of different strategies, driven by some general
mechanisms. In effect it is shown that even in the case when fluxes between
roles are driven by selectively neutral mechanisms (acting in the same way on
all strategies), the differences in the mortalities in the focal interaction
will lead to different distributions of roles for different strategies. This
leads to more complex rules for cooperation than the classical Hamilton's rule.


  The paper contains the attempt to integration of the classical evolutionary
game theory based on replicator dynamics and the state based approach of
Houston and Mcnamara. In the new approach, individuals have different heritable
strategies, however the individuals carrying the same strategy can differ on
the state, role or situation in which they act. Thus, the classical replicator
dynamics is completed by the additional subsystem of differential equations
describing the dynamics of transitions between different states. In effect the
interactions described by game structure, in addition to the demographic
payoffs (constituted by births and deaths) can lead to the change of state of
the competing individuals. The special cases of the new framework of stage
structured models where the state changes describe developmental steps or aging
are derived. New approach is illustrated by the example of Owner-Intruder game
with explicit dynamics of the role changes. New model is the generalization of
the demographic version of the Hawk-Dove game, the difference is that opponents
in the game are drawn from two separate subpopulations consisting of Owners and
Intruders. Intruders check random nest sites, and play the Hawk-Dove game with
the Owner if they are occupied. Interesting feedback mechanism is produced by
fluxes of individuals between subpopulations. Owners produce newborns which
become Intruders, since they should find a free nest site to reproduce.


  The explosion of data throughout the biomedical sciences provides
unprecedented opportunities to learn about the dynamics of evolution and
disease progression, but harnessing these large and diverse datasets remains
challenging. Here, we describe a highly generalisable statistical platform to
infer the dynamic pathways by which many, potentially interacting, discrete
traits are acquired or lost over time in biomedical systems. The platform uses
HyperTraPS (hypercubic transition path sampling) to learn progression pathways
from cross-sectional, longitudinal, or phylogenetically-linked data with
unprecedented efficiency, readily distinguishing multiple competing pathways,
and identifying the most parsimonious mechanisms underlying given observations.
Its Bayesian structure quantifies uncertainty in pathway structure and allows
interpretable predictions of behaviours, such as which symptom a patient will
acquire next. We exploit the model's topology to provide visualisation tools
for intuitive assessment of multiple, variable pathways. We apply the method to
ovarian cancer progression and the evolution of multidrug resistance in
tuberculosis, demonstrating its power to reveal previously undetected dynamic
pathways.


  The domestic dog is known to have evolved from gray wolves, about 15,000
years ago. They majorly exist as free-ranging populations across the world.
They are typically scavengers and well adapted to living among humans. Most
canids living in and around urban habitats tend to avoid humans and show
crepuscular activity peaks. In this study, we carried out a detailed
population-level survey on free-ranging dogs in West Bengal, India, to
understand the activity patterns of free-ranging dogs in relation to human
activity. Using 5669 sightings of dogs, over a period of 1 year, covering the
24 hours of the day, we carried out an analysis of the time-activity budget of
free-ranging dogs to conclude that they are generalists in their habit. They
remain active when humans are active. Their activity levels are affected
significantly by age class and time of the day. Multivariate analysis revealed
the presence of certain behavioural clusters on the basis of time of the day
and energy expenditure in the behaviours. In addition, we provide a detailed
ethogram of free-ranging dogs. This, to our knowledge, is the first study of
this kind, which might be used to further study the eco-ethology of these dogs.


  This article is a presentation of specific recent results describing scaling
limits of individual-based models. Thanks to them, we wish to relate the
time-scales typical of demographic dynamics and natural selection to the
parameters of the individual-based models. Although these results are by no
means exhaustive, both on the mathematical and the biological level, they
complement each other. Indeed, they provide a viewpoint for many classical
time-scales. Namely, they encompass the timescale typical of the
life-expectancy of a single individual, the longer one wherein a population can
be characterized through its demographic dynamics, and at least four
interconnected ones wherein selection occurs. The limiting behavior is
generally deterministic. Yet, since there are selective effects on randomness
in the history of lineages, probability theory is shown to be a key factor in
understanding the results. Besides, randomness can be maintained in the
limiting dynamics, for instance to model rare mutations fixing in the
population.


  Evidence of human activities during the Middle to Upper Palaeolithic
transition is well represented from rock_shelters, caves and open_air sites
across Italy. Over the past decade, both the revision of taphonomic processes
affecting archaeological faunal assemblages and new zooarchaeological studies
have allowed archaeologists to better understand subsistence strategies and
cultural behaviors attributed to groups of Neandertal and modern humans living
in the region. This work presents the preliminary results of a 5 years research
programme (ERC n. 724046_SUCCESS) and offers a state_of_the_art synthesis of
archaeological faunal assemblages including mammals and birds uncovered in
Italy between 50 and 35 ky ago. The present data were recovered in primary Late
Mousterian, Uluzzian, and Protoaurignacian stratigraphic contexts from Northern
Italy (Grotta di Fumane, Riparo del Broion, Grotta Maggiore di San Bernardino,
Grotta del Rio Secco, Riparo Bombrini), and Southern Italy (Grotta di
Castelcivita, Grotta della Cala, Grotta del Cavallo, and Riparo
l'Oscurusciuto). The available Number of Identified Specimens (NISP) is
analysed through intra- and inter-site comparisons at a regional scale, while
aoristic analysis is applied to the sequence documented at Grotta di Fumane.
Results of qualitative comparisons suggest an increase in the number of hunted
taxa since the end of the Middle Palaeolithic, and a marked change in
ecological settings beginning with the Protoaurignacian, with a shift to lower
temperatures and humidity. The distribution of carnivore remains and taphonomic
analyses hint at a possible change in faunal exploitation and butchering
processing between the Middle and Upper Palaeolithic. A preliminary comparison
between bone frequencies and the distribution of burned bones poses interesting
questions concerning the management of fire. Eventually, the combined use
of...(continue)


  The prediction of critical transitions, such as extinction events, is vitally
important to preserving vulnerable populations in the face of a rapidly
changing climate and continuously increasing human resource usage. Predicting
such events in spatially distributed populations is challenging because of the
high dimensionality of the system and the complexity of the system dynamics.
Here, we reduce the dimensionality of the problem by quantifying spatial
patterns via Betti numbers ($\beta_0$ and $\beta_1$), which count particular
topological features in a topological space. Spatial patterns representing
regions occupied by the population are analyzed in a coupled patch population
model with Ricker map growth and nearest-neighbors dispersal on a
two-dimensional lattice. We illustrate how Betti numbers can be used to
characterize spatial patterns by type, which in turn may be used to track
spatiotemporal changes via Betti number time series and characterize asymptotic
dynamics of the model parameter space. En route to a global extinction event,
we find that the Betti number time series of a population exhibits
characteristic changes. We hope these preliminary results will be used to aide
in the prediction of critical transitions in spatially extended systems.
Additional applications of this technique include analysis of spatial data
(e.g., GIS) and model validation.


  Defining the processes involved in the technical/cultural shifts from the
Late Middle to the Early Upper Palaeolithic in Europe (~50-39 thousand years
BP) is one of the most important tasks facing prehistoric studies. In this
debate Italy plays a pivotal role, due to its geographical position between
eastern and western Mediterranean Europe as well as to it being the location of
several sites showing Late Mousterian, Uluzzian and Protoaurignacian evidence
distributed across the Peninsula. Our study aims to provide a synthesis of the
available lithic evidence from this key area through a review of the evidence
collected from a number of reference sites. The main technical features of the
Late Mousterian, the Uluzzian and the Protoaurignacian traditions are examined
from a diachronic and spatial perspective. Our overview allows the
identification of major differences in the technological behaviour of these
populations, making it possible to propose a number of specific working
hypotheses on the basis of which further studies can be carried out. This study
presents a detailed comparative study of the whole corpus of the lithic
production strategies documented during this interval, and crucial element thus
emerge.These data are of primary importance in order to assess the nature of
the "transition" phenomenon in Italy, thus contributing to the larger debate
about the disappearance of Neandertals and the arrival of early Modern Humans
in Europe.


  Most of the DNA that composes a complex organism is non-coding and defined as
junk. Even the coding part is composed of genes that affect the phenotype
differently. Therefore, a random mutation has an effect on the specimen fitness
that strongly depends on the DNA region where it occurs. Such heterogeneous
composition should be linked to the evolutionary process. However, the way is
still unknown.
  Here, we study a minimal model for the evolution of an ecosystem where two
antagonist species struggle for survival on a lattice. Each specimen possesses
a toy genome, encoding for its phenotype. The gene pool of populations changes
in time due to the effect of random mutations on genes (entropic force) and of
interactions with the environment and between individuals (natural selection).
  We prove that the relevance of each gene in the manifestation of the
phenotype is a key feature for evolution. In the presence of a uniform gene
relevance, a mutational meltdown is observed. Natural selection acts quenching
the ecosystem in a non-equilibriumstate that slowly drifts, decreasing the
fitness and leading to the extinction of the species. Conversely, if a specimen
is provided with a heterogeneous gene relevance, natural selection wins against
entropic forces, and the species evolves increasing its fitness. We finally
show that heterogeneity together with spatial correlations is responsible for
spontaneous sympatric speciation.


  A blood cell lineage consists of several consecutive developmental stages
from the pluripotent or multipotent stem cell to a particular stage of
terminally differentiated cells. There is considerable interest in identifying
the key regulatory genes that govern blood cell development from the gene
expression data without considering the underlying network between
transcription factors (TFs) and their target genes. In this study, we introduce
a novel expression pattern that key regulators expose along the differentiation
path. We deploy this pattern to identify the cell-specific key regulators
responsible for the development. As proof of concept, we consider this approach
to data on six developmental stages from mouse embryonic stem cells to
terminally differentiated macrophages.


  An important goal in microbial computational genomics is to identify crucial
events in the evolution of a gene that severely alter the duplication, loss and
mobilization patterns of the gene within the genomes in which it disseminates.
In this paper, we formalize this microbiological goal as a new pattern-matching
problem in the domain of Gene tree and Species tree reconciliation, denoted
"Reconciliation-Scenario Altering Mutation (RSAM) Discovery". We propose an
$O(m\cdot n\cdot k)$ time algorithm to solve this new problem, where $m$ and
$n$ are the number of vertices of the input Gene tree and Species tree,
respectively, and $k$ is a user-specified parameter that bounds from above the
number of optimal solutions of interest. The algorithm first constructs a
hypergraph representing the $k$ highest scoring reconciliation scenarios
between the given Gene tree and Species tree, and then interrogates this
hypergraph for subtrees matching a pre-specified RSAM Pattern. Our algorithm is
optimal in the sense that the number of hypernodes in the hypergraph can be
lower bounded by $\Omega(m\cdot n\cdot k)$. We implement the new algorithm as a
tool, called RSAM-finder, and demonstrate its application to -the
identification of RSAMs in toxins and drug resistance elements across a dataset
spanning hundreds of species.


  The genomic profile underlying an individual tumor can be highly informative
in the creation of a personalized cancer treatment strategy for a given
patient; a practice known as precision oncology. This involves next generation
sequencing of a tumor sample and the subsequent identification of genomic
aberrations, such as somatic mutations, to provide potential candidates of
targeted therapy. The identification of these aberrations from sequencing noise
and germline variant background poses a classic classification-style problem.
This has been previously broached with many different supervised machine
learning methods, including deep-learning neural networks. However, these
neural networks have thus far not been tailored to give any indication of
confidence in the mutation call, meaning an oncologist could be targeting a
mutation with a low probability of being true. To address this, we present here
a deep bayesian recurrent neural network for cancer variant calling, which
shows no degradation in performance compared to standard neural networks. This
approach enables greater flexibility through different priors to avoid
overfitting to a single dataset. We will be incorporating this approach into
software for oncologists to obtain safe, robust, and statistically confident
somatic mutation calls for precision oncology treatment choices.


  Recent theoretical research proposes that computational complexity can be
seen as an ultimate constraint that allows for open-ended biological evolution
on finite static fitness landscapes. Whereas on easy fitness landscapes,
evolution will quickly converge to a local fitness peaks, on hard fitness
landscapes this computational constraints prevents evolution from reaching any
local fitness peak in polynomial time. Valued constraint satisfaction problems
(VCSPs) can be used to represent both easy and hard fitness landscapes. Thus
VCSPS can be seen as a natural way of linking the theory of evolution with
notions of computer science to better understand the features that make
landscapes hard. However, there are currently no simulators that study
VCSP-structured fitness landscapes.
  This report describes the design and build of an evolution simulator for
VCSP-structured fitness landscapes. The platform is used for simulating various
instances of easy and hard fitness landscapes. In particular, we look at
evolution under more realistic assumptions than fittest mutant strong-selection
weak mutation dynamics on the winding semismooth fitness landscape. The results
obtained match with the theoretical expectations, while also providing new
information about the limits of evolution. The last part of the report
introduces a mathematical model for smooth fitness landscapes and uses it to
better understand why these landscapes are easy.


  Modelling the substitution of nucleotides along a phylogenetic tree is
usually done by a hidden Markov process. This allows to define a distribution
of characters at the leaves of the trees and one might be able to obtain
polynomial relationships among the probabilities of different characters. The
study of these polynomials and the geometry of the algebraic varieties defined
by them can be used to reconstruct phylogenetic trees. However, not all points
in these algebraic varieties have biological sense. In this paper, we explore
the extent to which adding semi-algebraic conditions arising from the
restriction to parameters with statistical meaning can improve existing methods
of phylogenetic reconstruction. To this end, our aim is to compute the distance
of data points to algebraic varieties and to the stochastic part of these
varieties. Computing these distances involves optimization by nonlinear
programming algorithms. We use analytical methods to find some of these
distances for quartet trees evolving under the Kimura 3-parameter or the
Jukes-Cantor models. Numerical algebraic geometry and computational algebra
play also a fundamental role in this paper.


  To date, twelve complete genomes representing eleven species belonging to six
genera have been sequenced in salmonids. For the genus Salvelinus, it was
supposed to sequence the genome of Arctic char, one of the most variable
species of vertebrate animals. Sequencing was carried out (Christensen et al.,
2018) using the tissues of the female IW2-2015 obtained from the company
engaged in industrial aquaculture of chars - Icy Waters Ltd. The company
exploits two of its own broodstocks - NL and TR, originating from the chars
from the Nauyuk Lake and the Tree River (Nunavut, Canada). Since the complete
mitochondrial genome of the female IW2-2015 was absent in the published
assembly ASM291031v2, we determined its type and complete sequence from the
sequence read archives taken from Genbank. It was found that the female's
mitogenome belongs to the BERING haplogroup, which is characteristic of
Northern Dolly Varden S. malma malma. Analysis of other unlinked diagnostic
loci encoded by nuclear DNA (ITS1, RAG1, SFO-12, SFO-18, SMM-21) also revealed
distinctive characters of Northern Dolly Varden in female IW2-2015. It was
concluded that the genomic assembly ASM291031v2 was obtained not from an
individual of Arctic char S. alpinus, but from an individual of a related
species - Northern Dolly Varden S. malma malma. The identical to the IW2-2015
female characteristics of diagnostic loci were found in other individuals from
the broodstock TR. Apparently, the broodstock TR is entirely a strain derived
from Northern Dolly Varden. Since assembly ASM291031v2 was obtained from a
specimen originated from the marginal population of Northern Dolly Varden (Tree
R.) isolated from the main range of the species and with some traces of
introgressive hybridization, this assembly can hardly be considered as a
description of a typical genome of S. malma malma.


  There are many experimental data indicating the correlations of the changes
in high level of organization of chromatin in human cells and changes in the
state of the whole organism related to disease, state of tiredness or aging. In
our previous work: arXiv.org-2018 (1812.00186) we analyzed the publications on
the topic up to 2017. In this work we focused on works upon the problem of
connection of the state of chromatin with human diseases published in
2018-2019. In the modern literature the most attention is paid to problem of
chromatin transformations in different forms of cancer, Alzheimer'r disease,
and hereditary diseases. Summing up, the tendency of scientific research of
noncommunicable diseases is shifting towards investigation of aspects of
nuclear regulation of disease origin, connected with conformation of chromatin.


  There are many proteins or protein complexes which have multiple DNA binding
domains. This allows them to bind to multiple points on a DNA molecule (or
chromatin fibre) at the same time. There are also many proteins which have been
found to be able to compact DNA in vitro, and many others have been observed in
foci or puncta when fluorescently labelled and imaged in vivo. In this work we
study, using coarse-grained Langevin dynamics simulations, the compaction of
polymers by simple model proteins and a phenomenon known as the
"bridging-induced attraction". The latter is a mechanism observed in previous
simulations [Brackley et al., Proc. Natl. Acad. Sci. USA 110 (2013)], where
proteins modelled as spheres form clusters via their multivalent interactions
with a polymer, even in the absence of any explicit protein-protein attractive
interactions. Here we extend this concept to consider more detailed model
proteins, represented as simple "patchy particles" interacting with a
semi-flexible bead-and-spring polymer. We find that both the compacting ability
and the effect of the bridging-induced attraction depend on the valence of the
model proteins. These effects also depend on the shape of the protein, which
determines its ability to form bridges.


  Biological and cultural inspired optimization algorithms are nowadays part of
the basic toolkit of a great many research domains. By mimicking processes in
nature and animal societies, these general-purpose search algorithms promise to
deliver optimal or near-optimal solutions using hardly any information on the
optimization problems they are set to tackle. Here we study the performances of
a cultural-inspired algorithm -- the imitative learning search -- as well as of
asexual and sexual variants of evolutionary algorithms in finding the global
maxima of NK-fitness landscapes. The main performance measure is the total
number of agent updates required by the algorithms to find those global maxima
and the baseline performance, which establishes the effectiveness of the
cooperative algorithms, is set by the blind search in which the agents explore
the problem space (binary strings) by flipping bits at random. We find that
even for smooth landscapes that exhibit a single maximum, the evolutionary
algorithms do not perform much better than the blind search due to the
stochastic effects of the genetic roulette. The imitative learning is immune to
this effect thanks to the deterministic choice of the fittest string in the
population, which is used as a model for imitation. The tradeoff is that for
rugged landscapes the imitative learning search is more prone to be trapped in
local maxima than the evolutionary algorithms. In fact, in the case of rugged
landscapes with a mild density of local maxima, the blind search either beats
or matches the cooperative algorithms regardless of whether the task is to find
the global maximum or to find the fittest state within a given runtime.


  In this article, we propose a covariance based method for combining partial
data sets in the genotype to phenotype spectrum. In particular, an
expectation-maximization algorithm that can be used to combine partially
overlapping relationship/covariance matrices is introduced. Combining data this
way, based on relationship matrices, can be contrasted with a feature
imputation based approach. We used several public genomic data sets to explore
the accuracy of combining genomic relationship matrices. We have also used the
heterogeneous genotype/phenotype data sets in the https://triticeaetoolbox.org/
to illustrate how this new method can be used in genomic prediction, phenomics,
and graphical modeling.


  Evolutionary algorithms, inspired by natural evolution, aim to optimize
difficult objective functions without computing derivatives. Here we detail the
relationship between population genetics and evolutionary optimization and
formulate a new evolutionary algorithm. Optimization of a continuous objective
function is analogous to searching for high fitness phenotypes on a fitness
landscape. We summarize how natural selection moves a population along the
non-euclidean gradient that is induced by the population on the fitness
landscape (the natural gradient). Under normal approximations common in
quantitative genetics, we show how selection is related to Newton's method in
optimization. We find that intermediate selection is most informative of the
fitness landscape. We describe the generation of new phenotypes and introduce
an operator that recombines the whole population to generate variants that
preserve normal statistics. Finally, we introduce a proof-of-principle
algorithm that combines natural selection, our recombination operator, and an
adaptive method to increase selection. Our algorithm is similar to covariance
matrix adaptation and natural evolutionary strategies in optimization, and has
similar performance. The algorithm is extremely simple in implementation with
no matrix inversion or factorization, does not require storing a covariance
matrix, and may form the basis of more general model-based optimization
algorithms with natural gradient updates.


  The generalized Lotka-Volterra (gLV) equations are a mathematical proxy for
ecological dynamics. We focus on a gLV model of the gut microbiome, in which
the evolution of the gut microbial state is determined in part by pairwise
inter-species interaction parameters that encode environmentally-mediated
resource competition between microbes. We develop an in silico method that
controls the steady-state outcome of the system by adjusting these interaction
parameters. This approach is confined to a bistable region of the gLV model.
The two steady states of interest are idealized as either a "healthy" or
"diseased" steady state of the gut microbiome. In this method, a dimensionality
reduction technique called steady-state reduction (SSR) is first used to
generate a two-dimensional (2D) gLV model that approximates the
high-dimensional dynamics on the 2D subspace spanned by the two steady states.
Then a bifurcation analysis of the 2D model analytically determines parameter
modifications that drive a disease-prone initial condition to the healthy
steady state. This parameter modification of the reduced 2D model guides
parameter modifications of the original high-dimensional model, resulting in a
change of steady-state outcome in the high-dimensional model. This control
method, called SPARC (SSR-guided parameter change), bypasses the computational
challenge of directly determining parameter modifications in the original
high-dimensional system. SPARC could guide the development of indirect
bacteriotherapies, which seek to change microbial compositions by deliberately
modifying gut environmental variables such as gut acidity or macronutrient
availability.


  The time it takes the fastest searcher out of $N\gg1$ searchers to find a
target determines the timescale of many physical, chemical, and biological
processes. This time is called an extreme first passage time (FPT) and is
typically much faster than the FPT of a single searcher. Extreme FPTs of
diffusion have been studied for decades, but little is known for other types of
stochastic processes. In this paper, we study the distribution of extreme FPTs
of piecewise deterministic Markov processes (PDMPs). PDMPs are a broad class of
stochastic processes that evolve deterministically between random events. Using
classical extreme value theory, we prove general theorems which yield the
distribution and moments of extreme FPTs in the limit of many searchers based
on the short time distribution of the FPT of a single searcher. We then apply
these theorems to some canonical PDMPs, including run and tumble searchers in
one, two, and three space dimensions. We discuss our results in the context of
some biological systems and show how our approach accounts for an unphysical
property of diffusion which can be problematic for extreme statistics.


  In static monitoring cameras, useful contextual information can stretch far
beyond the few seconds typical video understanding models might see: subjects
may exhibit similar behavior over multiple days, and background objects remain
static. Due to power and storage constraints, sampling frequencies are low,
often no faster than one frame per second, and sometimes are irregular due to
the use of a motion trigger. In order to perform well in this setting, models
must be robust to irregular sampling rates. In this paper we propose a method
that leverages temporal context from the unlabeled frames of a novel camera to
improve performance at that camera. Specifically, we propose an attention-based
approach that allows our model, Context R-CNN, to index into a long term memory
bank constructed on a per-camera basis and aggregate contextual features from
other frames to boost object detection performance on the current frame.
  We apply Context R-CNN to two settings: (1) species detection using camera
traps, and (2) vehicle detection in traffic cameras, showing in both settings
that Context R-CNN leads to performance gains over strong baselines. Moreover,
we show that increasing the contextual time horizon leads to improved results.
When applied to camera trap data from the Snapshot Serengeti dataset, Context
R-CNN with context from up to a month of images outperforms a single-frame
baseline by 17.9% mAP, and outperforms S3D (a 3d convolution based baseline) by
11.2% mAP.


  An explicit bound is given for the Kolmogorov distance between a mixture of
normal distributions and a normal distribution with properly chosen parameter
values. A random variable X has a mixture of normal distributions if its
conditional distribution given some sigma-algebra is normal. The bound depends
only on the first two moments of the first two conditional moments of X given
this sigma-algebra. As an application, the Yule-Ornstein-Uhlenbeck model, used
in the field of phylogenetic comparative methods, is considered. A bound is
derived for the Kolmogorov distance between the distribution of the average
value of a phenotypic trait over n related species and a normal distribution.
The bound goes to 0 as n goes to infinity, extending earlier limit theorems by
Bartoszek and Sagitov.


  The pace and unpredictability of evolution are critically relevant in a
variety of modern challenges: combating drug resistance in pathogens and
cancer, understanding how species respond to environmental perturbations like
climate change, and developing artificial selection approaches for agriculture.
Great progress has been made in quantitative modeling of evolution using
fitness landscapes, allowing a degree of prediction for future evolutionary
histories. Yet fine-grained control of the speed and the distributions of these
trajectories remains elusive. We propose an approach to achieve this using
ideas originally developed in a completely different context: counterdiabatic
driving to control the behavior of quantum states for applications like quantum
computing and manipulating ultra-cold atoms. Implementing these ideas for the
first time in a biological context, we show how a set of external control
parameters (i.e. varying drug concentrations / types, temperature, nutrients)
can guide the probability distribution of genotypes in a population along a
specified path and time interval. This level of control, allowing empirical
optimization of evolutionary speed and trajectories, has myriad potential
applications, from enhancing adaptive therapies for diseases, to the
development of thermotolerant crops in preparation for climate change, to
accelerating bioengineering methods built on evolutionary models, like directed
evolution of biomolecules.


  During phloem unloading, multiple cell-to-cell transport events move organic
substances to the root meristem. Although the primary unloading event from the
sieve elements to the phloem pole pericycle has been characterized to some
extent, little is known about post-sieve element unloading. Here, we report a
novel gene, PHLOEM UNLOADING MODULATOR (PLM), in the absence of which
plasmodesmata-mediated symplastic transport through the phloem pole
pericycle--endodermis interface is specifically enhanced. Increased unloading
is attributable to a defect in the formation of the endoplasmic
reticulum--plasma membrane tethers during plasmodesmal morphogenesis, resulting
in the majority of pores lacking a visible cytoplasmic sleeve. PLM encodes a
putative enzyme required for the biosynthesis of sphingolipids with
very-long-chain fatty acid. Taken together, our results indicate that
post-sieve element unloading involves sphingolipid metabolism, which affects
plasmodesmal ultrastructure. They also raise the question of how and why
plasmodesmata with no cytoplasmic sleeve facilitate molecular trafficking.


  Phylogenetic networks are important for the study of evolution. The number of
methods to find such networks is increasing, but most such methods can only
reconstruct small networks. To find bigger networks, one can attempt to combine
small networks. In this paper, we study the {\sc Network Hybridization}
problem, a problem of combining networks into another network with low
complexity. We characterize this complexity via a restricted problem, {\sc
Tree-child Network Hybridization}, and we present an FPT algorithm to
efficiently solve this restricted problem.


  In Diffusive Molecular Communication (DMC), information is transmitted by
diffusing molecules. Synaptic signaling is a natural implementation of this
paradigm. It is responsible for relaying information from one neuron to
another, but also provides support for complex functionalities, such as
learning and memory. Many of its features are not yet understood, some are,
however, known to be critical for robust, reliable neural communication. In
particular, some synapses feature a re-uptake mechanism at the presynaptic
neuron, which provides a means for removing neurotransmitters from the synaptic
cleft and for recycling them for future reuse. In this paper, we develop a
comprehensive channel model for synaptic DMC encompassing a spatial model of
the synaptic cleft, molecule re-uptake at the presynaptic neuron, and
reversible binding to individual receptors at the postsynaptic neuron. Based on
this model, we derive an analytical time domain expression for the channel
impulse response (CIR) of the synaptic DMC system. Our model explicitly
incorporates macroscopic physical channel parameters and can be used to
evaluate the impact of re-uptake, receptor density, and channel width on the
CIR of the synaptic DMC system. Furthermore, we provide results from
particlebased computer simulation, which validate the analytical model. The
proposed comprehensive channel model for synaptic DMC systems can be exploited
for the investigation of challenging problems, like the quantification of the
inter-symbol interference between successive synaptic signals and the design of
synthetic neural communication systems.


  Defining and identifying causal intervention effects for transmissible
infectious disease outcomes is challenging because a treatment -- such as a
vaccine -- given to one individual may affect the infection outcomes of others.
Epidemiologists have proposed causal estimands to quantify effects of
interventions under contagion using a two-person partnership model. These
simple conceptual models have helped researchers develop causal estimands
relevant to clinical evaluation of vaccine effects. However, many of these
partnership models are formulated under structural assumptions that preclude
realistic infectious disease transmission dynamics, limiting their conceptual
usefulness in defining and identifying causal treatment effects in empirical
intervention trials. In this paper, we propose causal intervention effects in
two-person partnerships under arbitrary infectious disease transmission
dynamics, and give nonparametric identification results showing how effects can
be estimated in empirical trials using time-to-infection or binary outcome
data. The key insight is that contagion is a causal phenomenon that induces
conditional independencies on infection outcomes that can be exploited for the
identification of clinically meaningful causal estimands. These new estimands
are compared to existing quantities, and results are illustrated using a
realistic simulation of an HIV vaccine trial.


  Plant-soil feedback is recognized as a causal mechanism for the emergence of
vegetation patterns of the same species especially when water is not a limiting
resource (e.g. humid environments). Nevertheless, in the field, plants rarely
grow in monoculture but compete with other plant species. In these cases,
plant-soil feedback was shown to play a key role in plant-species coexistence.
Using a mathematical model consisting of four PDEs, we investigate mechanisms
of inter- and intra-specific plant-soil feedback on the coexistence of two
competing plant species. In particular, the model takes into account both
negative and positive feedback influencing the growth of the same and the other
plant species. Both the coexistence of the plant species and the dominance of a
particular plant species is examined with respect to all model parameters
together with the emergence of spatial vegetation patterns.


  The emerging field of precision oncology relies on the accurate pinpointing
of alterations in the molecular profile of a tumor to provide personalized
targeted treatments. Current methodologies in the field commonly include the
application of next generation sequencing technologies to a tumor sample,
followed by the identification of mutations in the DNA known as somatic
variants. The differentiation of these variants from sequencing error poses a
classic classification problem, which has traditionally been approached with
Bayesian statistics, and more recently with supervised machine learning methods
such as neural networks. Although these methods provide greater accuracy,
classic neural networks lack the ability to indicate the confidence of a
variant call. In this paper, we explore the performance of deep Bayesian neural
networks on next generation sequencing data, and their ability to give
probability estimates for somatic variant calls. In addition to demonstrating
similar performance in comparison to standard neural networks, we show that the
resultant output probabilities make these better suited to the disparate and
highly-variable sequencing data-sets these models are likely to encounter in
the real world. We aim to deliver algorithms to oncologists for which model
certainty better reflects accuracy, for improved clinical application. By
moving away from point estimates to reliable confidence intervals, we expect
the resultant clinical and treatment decisions to be more robust and more
informed by the underlying reality of the tumor molecular profile.


  The bulk of stochastic gene expression models in the literature do not have
an explicit description of the age of a cell within a generation and hence they
cannot capture events such as cell division and DNA replication. Instead, many
models incorporate cell cycle implicitly by assuming that dilution due to cell
division can be described by an effective decay reaction with first-order
kinetics. If it is further assumed that protein production occurs in bursts
then the stationary protein distribution is a negative binomial. Here we seek
to understand how accurate these implicit models are when compared with more
detailed models of stochastic gene expression. We derive the exact stationary
solution of the chemical master equation describing bursty protein dynamics,
binomial partitioning at mitosis, age-dependent transcription dynamics
including replication, and random interdivision times sampled from Erlang or
more general distributions; the solution is different for single lineage and
population snapshot settings. We show that protein distributions are well
approximated by the solution of implicit models (a negative binomial) when the
mean number of mRNAs produced per cycle is low and the cell cycle length
variability is large. When these conditions are not met, the distributions are
either almost bimodal or else display very flat regions near the mode and
cannot be described by implicit models. We also show that for genes with low
transcription rates, the size of protein noise has a strong dependence on the
replication time, it is almost independent of cell cycle variability for
lineage measurements and increases with cell cycle variability for population
snapshot measurements. In contrast for large transcription rates, the size of
protein noise is independent of replication time and increases with cell cycle
variability for both lineage and population measurements.


  The aim of this study was to investigate human erythrocytes as a carrier for
targeted drug delivery of primaquine (PQ). The process of PQ loading in human
erythrocytes, as well as the effect of PQ loading on the oxidative status of
erythrocytes, was also studied. At PQ concentrations of 2, 4, 6, and 8 mg/mL
and an incubation time of 2 h, the ratios of the concentrations of PQ entrapped
in erythrocytes to that in the incubation medium were 0.515, 0.688, 0.697 and
0.788, respectively. The maximal decline of erythrocyte reduced glutathione
content was observed at 8 mg/mL of PQ compared with native erythrocytes p <
0.001. In contrast, malondialdehyde and protein carbonyl were significantly
increased in cells loaded with PQ (p < 0.001). Furthermore, osmotic fragility
of PQ carrier erythrocytes was increased in comparison with unloaded cells.
Electron microscopy revealed spherocyte formation with PQ carrier erythrocytes.
PQ-loaded cells showed sustained drug release over a 48 h period. Erythrocytes
were loaded with PQ successfully, but there were some biochemical as well as
physiological changes that resulted from the effect of PQ on the oxidative
status of drug-loaded erythrocytes. These changes may result in favorable
targeting of PQ-loaded cells to reticulo-endothelial organs. The relative
impact of these changes remains to be explored in ongoing animal studies.


  Stochasticity in gene expression can result in fluctuations in gene product
levels. Recent experiments indicated that feedback regulation plays an
important role in controlling the noise in gene expression. A quantitative
understanding of the feedback effect on gene expression requires analysis of
the corresponding stochastic model. However, for stochastic models of gene
expression with general regulation functions, exact analytical results for gene
product distributions have not been given so far. Here, we propose a technique
to solve a generalized ON-OFF model of stochastic gene expression with
arbitrary (positive or negative, linear or nonlinear) feedbacks including
posttranscriptional or posttranslational regulation. The obtained results,
which generalize results obtained previously, provide new insights into the
role of feedback in regulating gene expression. The proposed analytical
framework can easily be extended to analysis of more complex models of
stochastic gene expression.


  Recently, machine learning (ML) has established itself in various worldwide
benchmarking competitions in computational biology, including Critical
Assessment of Structure Prediction (CASP) and Drug Design Data Resource (D3R)
Grand Challenges. However, the intricate structural complexity and high ML
dimensionality of biomolecular datasets obstruct the efficient application of
ML algorithms in the field. In addition to data and algorithm, an efficient ML
machinery for biomolecular predictions must include structural representation
as an indispensable component. Mathematical representations that simplify the
biomolecular structural complexity and reduce ML dimensionality have emerged as
a prime winner in D3R Grand Challenges. This review is devoted to the recent
advances in developing low-dimensional and scalable mathematical
representations of biomolecules in our laboratory. We discuss three classes of
mathematical approaches, including algebraic topology, differential geometry,
and graph theory. We elucidate how the physical and biological challenges have
guided the evolution and development of these mathematical apparatuses for
massive and diverse biomolecular data. We focus the performance analysis on the
protein-ligand binding predictions in this review although these methods have
had tremendous success in many other applications, such as protein
classification, virtual screening, and the predictions of solubility, solvation
free energy, toxicity, partition coefficient, protein folding stability changes
upon mutation, etc.


  There are many ways to represent a molecule as input to a machine learning
model and each is associated with loss and retention of certain kinds of
information. In the interest of preserving three-dimensional spatial
information, including bond angles and torsions, we have developed libmolgrid,
a general-purpose library for representing three-dimensional molecules using
multidimensional arrays. This library also provides functionality for composing
batches of data suited to machine learning workflows, including data
augmentation, class balancing, and example stratification according to a
regression variable or data subgroup, and it further supports temporal and
spatial recurrences over that data to facilitate work with recurrent neural
networks, dynamical data, and size extensive modeling. It was designed for
seamless integration with popular deep learning frameworks, including Caffe,
PyTorch, and Keras, providing good performance by leveraging graphical
processing units (GPUs) for computationally-intensive tasks and efficient
memory usage through the use of memory views over preallocated buffers.
libmolgrid is a free and open source project that is actively supported,
serving the growing need in the molecular modeling community for tools that
streamline the process of data ingestion, representation construction, and
principled machine learning model development.


  The dynamics of a population exhibiting exponential growth can be modelled as
a birth-death process, which naturally captures the stochastic variation in
population size over time. In this article, we consider a supercritical
birth-death process, started at a random time in the past, and conditioned to
have n sampled individuals at the present. The genealogy of individuals sampled
at the present time is then described by the reversed reconstructed process
(RRP), which traces the ancestry of the sample backwards from the present. We
show that a simple, analytic, time rescaling of the RRP provides a
straightforward way to derive its inter-event times. The same rescaling
characterises other distributions underlying this process, obtained elsewhere
in the literature via more cumbersome calculations. We also consider the case
of incomplete sampling of the population, in which each leaf of the genealogy
is retained with an independent Bernoulli trial with probability $\psi$, and we
show that corresponding results for Bernoulli-sampled RRPs can be derived using
time rescaling, for any values of the underlying parameters. A central result
is the derivation of a scaling limit as $\psi$ approaches 0, corresponding to
the underlying population growing to infinity, using the time rescaling
formalism. We show that in this setting, after a linear time rescaling, the
event times are the order statistics of $n$ logistic random variables with mode
$\log(1/\psi)$; moreover, we show that the inter-event times are approximately
exponentially distributed.


  The underlying structure of the canonical amino acid substitution matrix
(aaSM) is examined by considering stepwise improvements in the differential
recognition of amino acids according to their chemical properties during the
branching history of the two aminoacyl-tRNA synthetase (aaRS) superfamilies.
The evolutionary expansion of the genetic code is described by a simple
parameterization of the aaSM, in which (i) the number of distinguishable amino
acid types, (ii) the matrix dimension, and (iii) the number of parameters, each
increases by one for each bifurcation in an aaRS phylogeny. Parameterized
matrices corresponding to trees in which the size of an amino acid sidechain is
the only discernible property behind its categorization as a substrate,
exclusively for a Class I or II aaRS, provide a significantly better fit to
empirically determined aaSM than trees with random bifurcation patterns. A
second split between polar and nonpolar amino acids in each Class effects a
vastly greater further improvement. The earliest Class-separated epochs in the
phylogenies of the aaRS reflect these enzymes' capability to distinguish tRNAs
through the recognition of acceptor stem identity elements via the minor (Class
I) and major (Class II) helical grooves, which is how the ancient Operational
Code functioned. The advent of tRNA recognition using the anticodon loop
supports the evolution of the optimal map of amino acid chemistry found in the
later Genetic Code, an essentially digital categorization, in which polarity is
the major functional property, compensating for the unrefined, haphazard
differentiation of amino acids achieved by the Operational Code.


  Neopetrosia exigua has received great attention in natural product chemistry.
The diversity of N. exigua constituents has been demonstrated by the continued
discovery of novel bioactive metabolites such as antimicrobial metabolites. In
this study, in order to localise the active component of N. exigua biomass
according to the polarity, a sequential gradient partition with different
solvents (nhexane, carbon tetrachloride, dichloromethane, n-butanol, and water)
was performed to obtain fractions containing metabolites distributed according
to their polarity. The antimicrobial activities of N. exigua fractions were
then evaluated using disc diffusion and microdilution methods (influence on the
growth curve, Minimum Inhibitory Concentration (MIC) and Minimum Bactericidal
Concentration (MBC)). The results showed that the active metabolites were
present in n-hexane, CH2Cl2, nBuOH, and water fractions. n-hexane, CH2Cl2, and
n-BuOH fractions were the most effective fractions. Among microbes tested,
Staphylococcus aureus was the most susceptible microbe evaluated. The obtained
results are considered sufficient for further study to isolate the compounds
represent the antimicrobial activity.


  In this chapter, an application of Mathematical Epidemiology to crop
vector-borne diseases is presented to investigate the interactions between
crops, vectors, and virus. The main illustrative example is the cassava mosaic
disease (CMD). The CMD virus has two routes of infection: through vectors and
also through infected crops. In the field, the main tool to control CMD
spreading is roguing. The presented biological model is sufficiently generic
and the same methodology can be adapted to other crops or crop vector-borne
diseases. After an introduction where a brief history of crop diseases and
useful information on Cassava and CMD is given, we develop and study a
compartmental temporal model, taking into account the crop growth and the
vector dynamics. A brief qualitative analysis of the model is provided,i.e.,
existence and uniqueness of a solution,existence of a disease-free equilibrium
and existence of an endemic equilibrium. We also provide conditions for local
(global) asymptotic stability and show that a Hopf Bifurcation may occur, for
instance, when diseased plants are removed. Numerical simulations are provided
to illustrate all possible behaviors. Finally, we discuss the theoretical and
numerical outputs in terms of crop protection.