  The margins within the geographic range of species are often specific in
terms of ecological and evolutionary processes, and can strongly influence the
species' reaction to climate change. One of the frequently observed features at
range margins is fragmentation, caused internally by population dynamics or
externally by the limited availability of suitable habitat sites. We study both
causes, and describe the transition from a connected to a fragmented state
across space by means of a gradient metapopulation model. The main features of
our approach are the following. 1) Inhomogeneities can occur at two spatial
scales: there is a broad-scale gradient, which can be patterned by fine-scale
heterogeneities. The latter is implemented by dispersing a variable number of
small obstacles over the terrain, which can be penetrable or unpenetrable by
the spreading species. 2) We study the occupancy of this terrain in a
steady-state on two temporal scales: in snapshots and by long-term averages.
The simulations reveal some general scaling laws that are applicable in various
environments, independently of the mechanism of fragmentation. The edge of the
connected region (the hull) is a fractal with dimension 7/4. Its width and
length changes with the gradient according to universal scaling laws, that are
characteristic for percolation transitions. The results suggest that
percolation theory is a powerful tool for understanding the structure of range
margins in a broad variety of real-life scenarios, including those in which the
environmental gradient is combined with fine-scale heterogeneity. This provides
a new method for comparing the range margins of different species in various
geographic regions, and monitoring range shifts under climate change.


  Networks are widely used to model the contact structure within a population
and in the resulting models of disease spread. While networks provide a high
degree of realism, the analysis of the exact model is out of reach and even
numerical methods fail for modest network size. Hence, mean-field models (e.g.
pairwise) focusing on describing the evolution of some summary statistics from
the exact model gained a lot of traction over the last few decades. In this
paper we revisit the problem of deriving triple closures for pairwise models
and we investigate in detail the assumptions behind some of the well-known
closures as well as their validity. Using a top-down approach we start at the
level of the entire graph and work down to the level of triples and combine
this with information around nodes and pairs. We use our approach to derive
many of the existing closures and propose new ones and theoretically connect
the two well-studied models of multinomial link and Poisson link selection. The
theoretical work is backed up by numerical examples to highlight where the
commonly used assumptions may fail and provide some recommendations for how to
choose the most appropriate closure when using graphs with no or modest degree
heterogeneity.


  Identifying directed interactions between species from time series of their
population densities has many uses in ecology. This key statistical task is
equivalent to causal time series inference, which connects to the Granger
causality (GC) concept: $x$ causes $y$ if $x$ improves the prediction of $y$ in
a dynamic model. However, the entangled nature of nonlinear ecological systems
has led to question the appropriateness of Granger causality, especially in its
classical linear Multivariate AutoRegressive (MAR) model form. Convergent-cross
mapping (CCM), a nonparametric method developed for deterministic dynamical
systems, has been suggested as an alternative. Here, we show that linear GC and
CCM are able to uncover interactions with surprisingly similar performance, for
predator-prey cycles, 2-species deterministic (chaotic) or stochastic
competition, as well as 10- and 20-species interaction networks. There is no
correspondence between the degree of nonlinearity of the dynamics and which
method performs best. Our results therefore imply that Granger causality, even
in its linear MAR($p$) formulation, is a valid method for inferring
interactions in nonlinear ecological networks; using GC or CCM (or both) can
instead be decided based on the aims and specifics of the analysis.


  In the past decade, Genome-Wide Association Studies (GWAS) have delivered an
increasingly broad view of the genetic basis of human phenotypic variation. One
of the major developments from GWAS is polygenic scores, a genetic predictor of
an individual's genetic predisposition towards a trait constructed from GWAS.
The success of GWAS and polygenic scores seems to suggest that we will soon be
able to settle debates about whether phenotypic differences among groups are
driven in part by genetics. However, answering these questions is more
complicated than it seems at first glance and touches on many old issues about
the interpretation of human genetic variation. In this perspective piece, I
outline the ways in which issues of causality, stratification,
gene-by-environment interactions, and divergence among groups all complicate
the interpretation of among-population polygenic score differences.


  The TAZ1 domain of CREB binding protein is crucial for transcriptional
regulation and recognizes multiple targets. The interactions between TAZ1 and
its specific targets are related to the cellular hypoxic negative feedback
regulation. Previous experiments reported that one of the TAZ1 targets CITED2
is an efficient competitor of another target HIF-1{\alpha}. Here by developing
the structure-based models of TAZ1 complexes we have uncovered the underlying
mechanisms of the competitions between HIF-1{\alpha} and CITED2 binding to
TAZ1. Our results are consistent with the experimental hypothesis on the
competition mechanisms and the apparent affinity. In addition, the simulations
prove the dominant position of forming TAZ1-CITED2 complex in both
thermodynamics and kinetics. For thermodynamics, TAZ1-CITED2 is the lowest
basin located on the free energy surface of binding in the ternary system. For
kinetics, the results suggest that CITED2 binds to TAZ1 faster than
HIF-1{\alpha}. Besides, the analysis of contact map and f values in this study
will be helpful for further experiments on TAZ1 systems.


  Because of the recent advances of genome sequences, a large number of human
genome sequences are available for the study of human genetics. Genome-wide
association studies typically focus on associations between single-nucleotide
polymorphisms and traits such as major human diseases. However, the statistical
power of classical single-marker association analysis for rare variants is
limited. To address the challenge, rare and low-frequency variants are often
grouped into a gene or pathway level, and the effects of multiple variants
evaluated based on collapsing methods. The sequential kernel association test
(SKAT) is one of the most effective collapsing methods. SKAT utilizes the
kernel matrix. The size of the kernel matrix is O(n^2), where the sample size
is n, so that the calculation of the data using the kernel method requires a
long time. As the sample sizes of human genetic studies increase, the
computational time is getting more and more problematic. In the present paper,
the genotype value decomposition method is proposed for the handling the
sequential kernel in a short period of time. The method can be referred to as
genotype value decomposition. In the present paper, it is shown that the
genetic relationship matrix and Identity by State (IBS) matrix can be obtained
using the genotype value vectors. By using this method, the SKAT can be
conducted with time complexity O(n). The proposed method enables to conduct
SKAT for samples of human genetics.


  In this paper we discuss a one-dimensional Abelian Higgs model with
Chern-Simons interaction as an effective theory of one-dimensional curves
embedded in three-dimensional space. We demonstrate how this effective model is
compatible with the geometry of protein molecules. Using standard field theory
techniques we analyze phenomenologically interesting static configurations of
the model and discuss their stability. This simple model predicts some
characteristic relations for the geometry of secondary structure motifs of
proteins, and we show how this is consistent with the experimental data. After
using the data to universally fix basic local geometric parameters, such as the
curvature and torsion of the helical motifs, we are left with a single free
parameter. We explain how this parameter controls the abundance and shape of
the principal motifs (alpha helices, beta strands and loops connecting them).


  Recent electronic transport experiments using metallic contacts attached to
proteins identified some 'stylized facts' which contradict conventional wisdom
that increasing either the spatial distance between the electrodes or the
temperature suppresses conductance exponentially. These include nearly
temperature independent conductance over the protein in the 30-300K range,
distance independent conductance within a single-protein in the 1-10 nm range
and an anomalously large conductance in the 0.1-10 nS range. In this paper we
develop a generalization of the low temperature Landauer formula which can
account for the joint effects of tunneling and decoherence and can explain
these new experimental findings. We use novel approximations which greatly
simplify the mathematical treatment and allow us to calculate the conductance
in terms of a handful macroscopic parameters instead of the myriads of
microscopic parameters describing the details of an atomic level quantum
chemical computation. The new approach makes it possible to get predictions for
the outcomes of new experiments without relying solely on high performance
computing and can distinguish important and unimportant details of the protein
structures from the point of view of transport properties.


  There has been tremendous development of linear controllability of complex
networks. Real-world systems are fundamentally nonlinear. Is linear
controllability relevant to nonlinear dynamical networks? We identify a common
trait underlying both types of control: the nodal "importance." For nonlinear
and linear control, the importance is determined, respectively, by
physical/biological considerations and the probability for a node to be in the
minimum driver set. We study empirical mutualistic networks and a gene
regulatory network, for which the nonlinear nodal importance can be quantified
by the ability of individual nodes to restore the system from the aftermath of
a tipping-point transition. We find that the nodal importance ranking for
nonlinear and linear control exhibits opposite trends: for the former
large-degree nodes are more important but for the latter, the importance scale
is tilted towards the small-degree nodes, suggesting strongly irrelevance of
linear controllability to these systems. The recent claim of successful
application of linear controllability to C. elegans connectome is examined and
discussed.


  We study unzipping of a complementary RNA-DNA helix applied to an external
force, focusing on the force-force correlations. While at the microscopic level
these are given by the sequence, the experiment measures effective, macroscopic
correlations. The latter are sequence-independent, i.e. universal, and
constitute the central object of the underlying field theory of disordered
systems. Comparing field-theory predictions and the exact solution of a 1-d toy
model with the experiment, we find an excellent agreement, confirming
fundamental theoretical concepts via a biologically inspired experiments.


  Engineering the entire genome of an organism enables large-scale changes in
organization, function, and external interactions, with significant
implications for industry, medicine, and the environment. Improvements to DNA
synthesis and organism engineering are already enabling substantial changes to
organisms with megabase genomes, such as Escherichia coli and Saccharomyces
cerevisiae. Simultaneously, recent advances in genome-scale modeling are
increasingly informing the design of metabolic networks. However, major
challenges remain for integrating these and other relevant technologies into
workflows that can scale to the engineering of gigabase genomes.
  In particular, we find that a major under-recognized challenge is
coordinating the flow of models, designs, constructs, and measurements across
the large teams and complex technological systems that will likely be required
for gigabase genome engineering. We recommend that the community address these
challenges by 1) adopting and extending existing standards and technologies for
representing and exchanging information at the gigabase genomic scale, 2)
developing new technologies to address major open questions around data
curation and quality control, 3) conducting fundamental research on the
integration of modeling and design at the genomic scale, and 4) developing new
legal and contractual infrastructure to better enable collaboration across
multiple institutions.


  This paper proposes small tree-width graph decomposition computational
protein design CFN instances defined according to the model [1] with protocol
defined by Simononcini et al [2] . The proteins used in the benchmark have been
selected in the PDB (not on their biological interest) to explore the
efficiency of global search method based on tree-width decomposition. The
instances are bigger than those previously proposed in the paper [2] with one
backbone relaxation and the aka Beta November 2016 Rosetta force-field [3]. The
benchmark includes 21 proteins selected with a low level of sequences identity
(40%) . Those instances have been selected on the basis of 3D criteria by
applying a decreasing average coarse volume occupancy filter by Amino Acid
(-i.e. by CFN variable) . The instances characteristic (see Table 1) contain
from 130 up to n = 282 variables with a maximum domain size from 383 to 438,
and between 1706 and 6208 cost functions. The min-fill tree-width ranges from
21 to 68, and from 0.16 to 0.34 for a normalized tree width. Those instances
have been used for UDGVNS search algorithm[4] benchmarking. This approach is
suitable for evaluation of search methods that exploit the notion of graph
decomposition.


  FRET experiments can yield state-specific structural information on complex
dynamic biomolecular assemblies. However, FRET experiments need to be combined
with computer simulations to overcome their sparsity. We introduce (i) an
automated FRET experiment design tool determining optimal FRET pairs for
structural modeling, (ii) a protocol for efficient FRET-assisted computational
structural modeling at multiple scales, and (iii) a quantitative quality
estimate for judging the accuracy of determined structures. We tested against
simulated and experimental data.


  I discuss the so-called stochastic individual based model of adaptive
dynamics and in particular how different scaling limits can be obtained by
taking limits of large populations, small mutation rate, and small effect of
single mutations together with appropriate time rescaling. In particular, one
derives the trait substitution sequence, polymorphic evolution sequence, and
the canonical equation of adaptive dynamics. In addition, I show how the escape
from an evolutionary stable conditions can occur as a metastable transition.
This is a review paper that will appear in
  "Probabilistic Structures in Evolution", ed. by E. Baake and A. Wakolbinger.


  In our Brief Communication (DOI: 10.1038/s41592-018-0033-z), we presented the
method SAVER for recovering true gene expression levels in noisy single cell
RNA sequencing data. We evaluated the performance of SAVER, along with
comparable methods MAGIC and scImpute, in an RNA FISH validation experiment and
a data downsampling experiment. In a Comment [arXiv:1908.07084v1], Li & Li were
concerned with the use of the downsampled datasets, specifically focusing on
clustering results obtained from the Zeisel et al. data. Here, we will address
these comments and, furthermore, amend the data downsampling experiment to
demonstrate that the findings from the data downsampling experiment in our
Brief Communication are valid.


  Background: There is limited information on the volume of
influenza-associated mortality in Russia.
  Methods: Using previously developed methodology (Goldstein et al.,
Epidemiology 2012), we regressed the monthly rates of mortality for respiratory
causes, circulatory causes, and for certain infectious and parasitic diseases
(available from the Russian Federal State Statistics Service (Rosstat))
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) on levels of ILI/ARI consultations and the percent of
respiratory specimens testing positive for influenza A/H3N2, A/H1N1 and B),
adjusting for the baseline rates of mortality not associated with influenza
circulation and temporal trends.
  Results: For the 2010/11 through the 2014/15 seasons, influenza circulation
was associated with an average annual 11106 (95% CI (4321,17892)) deaths for
circulatory causes, 4552 (3744,5360) deaths for respiratory causes, and 343
(63,624) deaths for certain infectious and parasitic diseases, with influenza
A/H3N2 being responsible for the majority of influenza-associated respiratory
deaths and influenza B making a substantial contribution to
influenza-associated mortality for circulatory causes.
  Conclusions: Influenza circulation is associated with a substantial mortality
burden in Russia, particularly for circulatory deaths. Those results support
the potential utility of influenza vaccination (with the role played by
influenza B in circulatory mortality pointing to the benefit of quadrivalent
influenza vaccines), as well as of administration of antiviral drugs for
high-risk individuals during periods of active influenza circulation.


  We investigate gene interaction networks in various cancer cells by spectral
analysis of the adjacency matrices. We observe localization of the networks on
hub genes which have extraordinarily many links. The eigenvector centralities
take finite values only on special nodes when the hub degree exceeds a critical
value $d_c \simeq 40$. The degree correlation function shows the disassortative
behavior in the large degrees, and the nodes whose degrees $d \gtrsim 40$ have
tendencies to link to small degree nodes. The communities of the gene networks
centered at the hub genes are extracted by the amount of node degree
discrepancies between linked nodes. We verify the Wigner-Dyson distribution of
the nearest neighbor eigenvalues spacing distribution $P(s)$ in the small
degree discrepancy communities, and the Poisson $P(s)$ in the communities of
large degree discrepancies including the hubs.


  Magombedze and Mulder (2013) studied the gene regulatory system of
\textit{Mycobacterium Tuberculosis} (\textit{Mtb}) by partitioning this into
three subsystems based on putative gene function and role in dormancy/latency
development. Each subsystem, in the form of $S$-system, is represented by an
embedded chemical reaction network (CRN), defined by a species subset and a
reaction subset induced by the set of digraph vertices of the subsystem. Based
on the network decomposition theory initiated by Feinberg in 1987, we have
introduced the concept of incidence-independent and developed the theory of
$\mathscr{C}$- and $\mathscr{C}^*$-decompositions including their structure
theorems in terms of linkage classes. With the $S$-system CRN $\mathscr{N}$ of
Magombedze and Mulder's \textit{Mtb} model, its reaction set partition induced
decomposition of subnetworks that are not CRNs of $S$-system but constitute
independent decomposition of $\mathscr{N}$. We have also constructed a new
$S$-system CRN $\mathscr{N}^*$ for which the embedded networks are
$\mathscr{C}^*$-decomposition. We have shown that subnetworks of $\mathscr{N}$
and the embedded networks (subnetworks of $\mathscr{N}^*$) are digraph
homomorphisms. Lastly, we attempted to explore modularity in the context of
CRN.


  To understand the function of cortical circuits it is necessary to classify
their underlying cellular diversity. Traditional attempts based on comparing
anatomical or physiological features of neurons and glia, while productive,
have not resulted in a unified taxonomy of neural cell types. The recent
development of single-cell transcriptomics has enabled, for the first time,
systematic high-throughput profiling of large numbers of cortical cells and the
generation of datasets that hold the promise of being complete, accurate and
permanent. Statistical analyses of these data have revealed the existence of
clear clusters, many of which correspond to cell types defined by traditional
criteria, and which are conserved across cortical areas and species. To
capitalize on these innovations and advance the field, we, the Copenhagen
Convention Group, propose the community adopts a transcriptome-based taxonomy
of the cell types in the adult mammalian neocortex. This core classification
should be ontological, hierarchical and use a standardized nomenclature. It
should be configured to flexibly incorporate new data from multiple approaches,
developmental stages and a growing number of species, enabling improvement and
revision of the classification. This community-based strategy could serve as a
common foundation for future detailed analysis and reverse engineering of
cortical circuits and serve as an example for cell type classification in other
parts of the nervous system and other organs.


  I present a Hamiltonian model and a computational method suitable to evaluate
structural and thermodynamic properties of helical molecules embedded in
crowded environments which may confine the space available to the base pair
fluctuations. It is shown that, for the specific case of a short DNA fragment
in a nanochannel, the molecule is markedly over-twisted and stretched by
narrowing the width of the channel.


  The concept that catalytic enzymes can act as molecular machines transducing
chemical activity into motion has conceptual and experimental support, but much
of the claimed support comes from experimental conditions where the substrate
concentration is higher than biologically relevant and accordingly exceeds kM,
the Michaelis-Menten constant. Moreover, many of the enzymes studied
experimentally to date are oligomeric. Urease, a hexamer of subunits, has been
considered to be the gold standard demonstrating enhanced diffusion. Here we
show that urease and certain other oligomeric enzymes of high catalytic
activity above kM dissociate into their smaller subunit fragments that diffuse
more rapidly, thus providing a simple physical mechanism of enhanced diffusion
in this regime of concentrations. Mindful that this conclusion may be
controversial, our findings are sup-ported by four independent analytical
techniques, static light scattering, dynamic light scattering (DLS),
size-exclusion chroma-tography (SEC), and fluorescence correlation spectroscopy
(FCS). Data for urease are presented in the main text and the con-clusion is
validated for hexokinase and acetylcholinesterase with data presented in
supplementary information. For substrate concentration regimes below kM at
which these enzymes do not dissociate, our findings from both FCS and DLS
validate that enzymatic catalysis does lead to the enhanced diffusion
phenomenon. INTRODUCT


  Extrinsic noise-induced transitions to bimodal dynamics have been largely
investigated in a variety of chemical, physical, and biological systems. In the
standard approach in physical and chemical systems, the key properties that
make these systems mathematically tractable are that the noise appears linearly
in the dynamical equations, and it is assumed Gaussian and white. In biology,
the Gaussian approximation has been successful in specific systems, but the
relevant noise being usually non-Gaussian, non-white, and nonlinear poses
serious limitations to its general applicability. Here we revisit the
fundamental features of linear Gaussian noise, pinpoint its limitations, and
review recent new approaches based on nonlinear bounded noises, which highlight
novel mechanisms to account for transitions to bimodal behaviour. We do this by
considering a simple but fundamental gene expression model, the repressed gene,
which is characterized by linear and nonlinear dependencies on external
parameters. We then review a general methodology introduced recently, so-called
nonlinear noise filtering, which allows the investigation of linear, nonlinear,
Gaussian and non-Gaussian noises. We also present a derivation of it, which
highlights its dynamical origin. Testing the methodology on the repressed gene
confirms that the emergence of noise-induced transitions appears to be strongly
dependent on the type of noise adopted, and on the degree of nonlinearity
present in the system.


  Gene regulatory networks play a crucial role in controlling an organism's
biological processes, which is why there is significant interest in developing
computational methods that are able to extract their structure from
high-throughput genetic data. Many of these computational methods are designed
to infer individual regulatory relationships among genes from data on gene
expression. We propose a novel efficient Bayesian method for discovering local
causal relationships among triplets of (normally distributed) variables. In our
approach, we score covariance structures for each triplet in one go and
incorporate available background knowledge in the form of priors to derive
posterior probabilities over local causal structures. Our method is flexible in
the sense that it allows for different types of causal structures and
assumptions. We apply our approach to the task of learning causal regulatory
relationships among genes. We show that the proposed algorithm produces stable
and conservative posterior probability estimates over local causal structures
that can be used to derive an honest ranking of the most meaningful regulatory
relationships. We demonstrate the stability and efficacy of our method both on
simulated data and on real-world data from an experiment on yeast.


  The discovery of important biomarkers is a significant step towards
understanding the molecular mechanisms of carcinogenesis; enabling accurate
diagnosis for, and prognosis of, a certain cancer type. Before recommending any
diagnosis, genomics data such as gene expressions(GE) and clinical outcomes
need to be analyzed. However, complex nature, high dimensionality, and
heterogeneity in genomics data make the overall analysis challenging.
Convolutional neural networks(CNN) have shown tremendous success in solving
such problems. However, neural network models are perceived mostly as `black
box' methods because of their not well-understood internal functioning.
However, interpretability is important to provide insights on why a given
cancer case has a certain type. Besides, finding the most important biomarkers
can help in recommending more accurate treatments and drug repositioning. In
this paper, we propose a new approach called OncoNetExplainer to make
explainable predictions of cancer types based on GE data. We used genomics data
about 9,074 cancer patients covering 33 different cancer types from the
Pan-Cancer Atlas on which we trained CNN and VGG16 networks using
guided-gradient class activation maps++(GradCAM++). Further, we generate
class-specific heat maps to identify significant biomarkers and computed
feature importance in terms of mean absolute impact to rank top genes across
all the cancer types. Quantitative and qualitative analyses show that both
models exhibit high confidence at predicting the cancer types correctly giving
an average precision of 96.25%. To provide comparisons with the baselines, we
identified top genes, and cancer-specific driver genes using gradient boosted
trees and SHapley Additive exPlanations(SHAP). Finally, our findings were
validated with the annotations provided by the TumorPortal.


  Boolean networks model finite discrete dynamical systems with complex
behaviours. The state of each component is determined by a Boolean function of
the state of (a subset of) the components of the network. This paper addresses
the synthesis of these Boolean functions from constraints on their domain and
emerging dynamical properties of the resulting network. The dynamical
properties relate to the existence and absence of trajectories between
partially observed configurations, and to the stable behaviours (fixpoints and
cyclic attractors). The synthesis is expressed as a Boolean satisfiability
problem relying on Answer-Set Programming with a parametrized complexity, and
leads to a complete non-redundant characterization of the set of solutions.
Considered constraints are particularly suited to address the synthesis of
models of cellular differentiation processes, as illustrated on a case study.
The scalability of the approach is demonstrated on random networks with
scale-free structures up to 100 to 1,000 nodes depending on the type of
constraints.


  The quantitative description of marine systems is constrained by a major
issue of scale separation: most marine biochemical processes occur at
sub-centimeter scales, while the contribution to the Earth's biogeochemical
cycles is expressed at much larger scales, up to the planetary one. In spite of
vastly improved computing power and observational capabilities, the modeling
approach has remained anchored to an old view that sees the microscales as
unable to substantially affect larger ones. The lack of a widespread
theoretical appreciation of the interactions between vastly different scales
has led to the proliferation of numerical models with uncertain predictive
capabilities. We show that an enhanced Lagrangian modeling framework, allowing
for those interactions, can easily tackle puzzling problems such as the
phenology of phytoplankton blooms, or vertical variability in mixed layers.


  We mathematically model the uptake of phosphorus by a soil community
consisting of a plant and two bacterial groups: copiotrophs and oligotrophs.
Four equilibrium states emerge, one for each of the species monopolising the
resource and dominating the community and one with coexistence of all species.
We show that the dynamics are controlled by the ratio of chemical adsorption to
bacterial death permitting either oscillatory states or quasi-steady uptake. We
show how a steady state can emerge which has soil and plant nutrient content
unresponsive to increased fertilization. However, the additional fertilization
supports the copiotrophs leading to community reassembly. Our results
demonstrate the importance of time-series measurements in nutrient uptake
experiments.


  In this work, we introduce a computational model for the study of the
host-bacteria interaction and the influence of the intestinal microbiota on the
behavior and feeding pattern of an individual. The model is based on digital
entities, which we've called Digital Beings (DBs), modeled using dynamic
systems and genetic algorithms. We have successfully tested the use of the DBs
by reproducing observation in previously made studies using rats and humans.
Among these studies, we highlight those on how the bacteria in an individual's
stomach could influence their eating behavior and how a controlled and
continuous diet can affect the longevity of a certain population. Our results
point that the Digital Beings can be used as a tool for supporting the devising
of experiments and corroborating with theoretical hypotheses, reducing the
number of in vivo tests.


  The treatment of short-tandem-repeat (STR) loci on the Y chromosome presents
special problems in the forensic analysis of DNA mixtures, chiefly but not
exclusively relating to the linkage of Y-STR loci which precludes the use of
the `product rule' for estimating Y-haplotype match probabilities. In recent
paper, Andersen and Balding(2017) estimated, via a population simulation model,
the distribution of the number of haplotypes sharing a common profile over a
set of Y-STR loci, and argued for its use as an alternative to estimating
Y-haplotype match probabilities.
  In this paper we present a sub-critical branching process model that
approximates their population model, and show how to estimate the haplotype
number distribution numerically using multivariate probability generating
functions. It is shown that the approximation provides a good fit to their
simulations. The model is extended to propose a new framework for evaluating
the weight-of-evidence of Y-STR haplotype mixtures, and it is illustrated with
publicly available data of a three person DNA mixture.


  The inverse Potts problem to infer a Boltzmann distribution for homologous
protein sequences from their single-site and pairwise amino acid frequencies
recently attracts a great deal of attention in the studies of protein structure
and evolution. We study regularization and learning methods and how to tune
regularization parameters to correctly infer interactions in Boltzmann machine
learning. Using $L_2$ regularization for fields, group $L_1$ for couplings is
shown to be very effective for sparse couplings in comparison with $L_2$ and
$L_1$. Two regularization parameters are tuned to yield equal values for both
the sample and ensemble averages of evolutionary energy. Both averages smoothly
change and converge, but their learning profiles are very different between
learning methods. The Adam method is modified to make stepsize proportional to
the gradient for sparse couplings and to use a soft-thresholding function for
group $L_1$. It is shown by first inferring interactions from protein sequences
and then from Monte Carlo samples that the fields and couplings can be well
recovered, but that recovering the pairwise correlations in the resolution of a
total energy is harder for the natural proteins than for the protein-like
sequences. Selective temperature for folding/structural constrains in protein
evolution is also estimated.


  With the advent of deep generative models in computational chemistry, in
silico anticancer drug design has undergone an unprecedented transformation.
While state-of-the-art deep learning approaches have shown potential in
generating compounds with desired chemical properties, they disregard the
genetic profile and properties of the target disease. Here, we introduce the
first generative model capable of tailoring anticancer compounds for a specific
biomolecular profile. Using a RL framework, the transcriptomic profiles of
cancer cells are used as a context for the generation of candidate molecules.
Our molecule generator combines two separately pretrained variational
autoencoders (VAEs) - the first VAE encodes transcriptomic profiles into a
smooth, latent space which in turn is used to condition a second VAE to
generate novel molecular structures on the given transcriptomic profile. The
generative process is optimized through PaccMann, a previously developed drug
sensitivity prediction model to obtain effective anticancer compounds for the
given context (i.e., transcriptomic profile). We demonstrate how the molecule
generation can be biased towards compounds with high predicted inhibitory
effect against individual cell lines or specific cancer sites. We verify our
approach by investigating candidate drugs generated against specific cancer
types and find the highest structural similarity to existing compounds with
known efficacy against these cancer types. We envision our approach to
transform in silico anticancer drug design by leveraging the biomolecular
characteristics of the disease in order to increase success rates in lead
compound discovery.


  Evidence accumulated over the past decade provides support for liquid-liquid
phase separation as the mechanism underlying the formation of biomolecular
condensates, which include not only membraneless organelles such as nucleoli
and RNA granules, but additional assemblies involved in transcription,
translation and signaling. Understanding the molecular mechanisms of condensate
function requires knowledge of the structures of their constituents. Current
knowledge suggests that structures formed via multivalent domain-motif
interactions remain largely unchanged within condensates. Two different
viewpoints exist regarding structures of disordered low-complexity domains
within condensates; one argues that low-complexity domains remain largely
disordered in condensates and their multivalency is encoded in short motifs
called stickers, while the other argues that the sequences form cross-beta
structures resembling amyloid fibrils. We review these viewpoints and highlight
outstanding questions that will inform structure-function relationships for
biomolecular condensates.


  Fetal neuroinflammation and prenatal stress (PS) may contribute to lifelong
neurological disabilities. Astrocytes and microglia, among the brain's
non-neuronal glia cell populations, play a pivotal role in neurodevelopment,
predisposition to and initiation of disease throughout lifespan. One of the
most common neurodevelopmental disorders manifesting between 1-4 years of age
is autism spectrum disorder (ASD). A pathological glial-neuronal interplay is
thought to increase the risk for clinical manifestation of ASD in at-risk
children, but the mechanisms remain poorly understood and integrative,
multi-scale models are needed. We propose a model that integrates the data
across the scales of physiological organization, from genome to phenotype, and
provides a foundation to explain the disparate findings on the genomic level.
We hypothesize that via gene-environment interactions, fetal neuroinflammation
and PS may reprogram glial immunometabolic phenotypes that impact
neurodevelopment and neurobehavior. Drawing on genomic data from the recently
published series of ovine and rodent glial transcriptome analyses with fetuses
exposed to neuroinflammation or PS, we conduct an analysis on the Simons
Foundation Autism Research Initiative (SFARI) Gene database. We confirm 21 gene
hits. Using unsupervised statistical network analysis, we then identify six
clusters of probable protein-protein interactions mapping onto the
immunometabolic and stress response networks and epigenetic memory. These
findings support our hypothesis. We discuss the implications for ASD etiology,
early detection, and novel therapeutic approaches. We conclude with delineation
of the next steps to verify our model on the individual gene level in an
assumption-free manner.


  Conformation of a protein (CorA) is examined in a matrix with mobile solute
constituents as a function of solute-residue interaction strength (f) by a
coarse-grained model with a Monte Carlo simulation. Solute particles are found
to reach their targeted residue due to their unique interactions with the
residues. Degree of slowing down of the protein depends on the interaction
strength f. Unlike a predictable dependence of the radius of gyration of the
same protein on interaction in an effective medium, it does not show a
systematic dependence on interaction due to pinning caused by the solute
binding. Spread of the protein chain is quantified by estimating its effective
dimension (D) from scaling of the structure factor. Even with a lower
solute-residue interaction, the protein chain appears to conform to a
random-coil conformation (D ~ 2) in its native phase where it is globular in
absence of such solute environment. The structural spread at small length scale
differs from that at large scale in presence of stronger interactions: D ~ 2.3
at smaller length scale and D ~ 1.4 on larger scale with f = 3.5 while D ~ 1.4
at smaller length scale and D ~ 2.5 at larger length scales with f = 4.0.


  Tree congruence metrics are typically global indices that describe the
similarity or dissimilarity between dendrograms. This study principally focuses
on topological congruence metrics that quantify similarity between two
dendrograms and can give a normalised score between 0 and 1. Specifically, this
article describes and tests two metrics the Clade Retention Index (CRI) and the
MASTxCF which is derived from the combined information available from a maximum
agreement subtree and a strict consensus. The two metrics were developed to
study differences between evolutionary trees, but their applications are
multidisciplinary and can be used on hierarchical cluster diagrams derived from
analyses in science, technology, maths or social sciences disciplines. A
comprehensive, but non-exhaustive review of other tree congruence metrics is
provided and nine metrics are further analysed. 1,620 pairwise analyses of
simulated dendrograms (which could be derived from any type of analysis) were
conducted and are compared in Pac-man piechart matrices. Kendalls tau-b is used
to demonstrate the concordance of the different metrics and Spearmans rho
ranked correlations are used to support these findings. The results support the
use of the CRI and MASTxCF as part of a suite of metrics, but it is recommended
that permutation metrics such as SPR distances and weighted metrics are
disregarded for the specific purpose of measuring similarity.


  We show that very simple molecular systems, modeled as chemical reaction
networks, can have behaviors that exhibit dramatic phase transitions at certain
population thresholds. Moreover, the magnitudes of these thresholds can thwart
attempts to use simulation, model checking, or approximation by differential
equations to formally verify the behaviors of such systems at realistic
populations. We show how formal theorem provers can successfully verify some
such systems at populations where other verification methods fail.


  With small-scale quantum processors transitioning from experimental physics
labs to industrial products, these processors allow us to efficiently compute
important algorithms in various fields. In this paper, we propose a quantum
algorithm to address the challenging field of big data processing for genome
sequence reconstruction. This research describes an architecture-aware
implementation of a quantum algorithm for sub-sequence alignment. A new
algorithm named QiBAM (quantum indexed bidirectional associative memory) is
proposed, that uses approximate pattern-matching based on Hamming distances.
QiBAM extends the Grover's search algorithm in two ways to allow for: (1)
approximate matches needed for read errors in genomics, and (2) a distributed
search for multiple solutions over the quantum encoding of DNA sequences. This
approach gives a quadratic speedup over the classical algorithm. A full
implementation of the algorithm is provided and verified using the OpenQL
compiler and QX simulator framework. This represents a first exploration
towards a full-stack quantum accelerated genome sequencing pipeline design. The
open-source implementation can be found on
https://github.com/prince-ph0en1x/QAGS.


  There is compelling evidence that substrate stiffness affects cell adhesion
as well as cytoskeleton organization and contractile activity. This work was
designed to study the cytoskeletal contractile activity of cells plated on
microposts of different stiffness using a numerical model simulating the
intracellular tension of individual cells. We allowed cells to adhere onto
micropost substrates of various rigidities and used experimental traction force
data to infer cell contractility using a numerical model. The model
discriminates between the influence of substrate stiffness on cell tension and
shows that higher substrate stiffness leads to an increase in intracellular
tension. The strength of this model is its ability to calculate the mechanical
state of each cell in accordance to its individual cytoskeletal structure. This
is achieved by regenerating a numerical cytoskeleton based


  The shape that stem cells reach at the end of adhesion process influences
their differentiation. Rearrangement of cytoskeleton and modification of
intracellular tension may activate mechanotransduction pathways controlling
cell commitment. In the present study, the mechanical signals involved in cell
adhesion were computed in in vitro stem cells of different shapes using a
single cell model, the so-called Cytoskeleton Divided Medium (CDM) model. In
the CDM model, the filamentous cytoskeleton and nucleoskeleton networks were
represented as a mechanical system of multiple tensile and compressive
interactions between the nodes of a divided medium. The results showed that
intracellular tonus, focal adhesion forces as well as nuclear deformation
increased with cell spreading. The cell model was also implemented to simulate
the adhesion process of a cell that spreads on protein-coated substrate by
emitting filopodia and creating new distant focal adhesion points. As a result,
the cell model predicted cytoskeleton reorganisation and reinforcement during
cell spreading. The present model quantitatively computed the evolution of
certain elements of mechanotransduction and may be a powerful tool for
understanding cell mechanobiology and designing biomaterials with specific
surface properties to control cell adhesion and differentiation.


  Statistical averaging and asynchronous population dynamics as portfolio
mechanisms are considered as the most important processes with which
biodiversity contributes to ecosystem stability. However, portfolio theories
usually regard biodiversity as a fixed property, but overlook the dynamics of
biodiversity altered by other ecosystem components. Here, we proposed a new
mechanistic food chain model with nutrient-diversity feedback to investigate
how dynamics of phytoplankton species diversity determines ecosystem stability.
Our model focuses on nutrient, community biomass of phytoplankton and
zooplankton, and phytoplankton species richness. The model assumes diversity
effects of phytoplankton on trophic interaction strength along plankton food
chain: phytoplankton diversity influences nutrient uptake by phytoplankton and
zooplankton grazing on phytoplankton, which subsequently affects nutrient level
and community biomass of phytoplankton and zooplankton. The nutrient level in
turn affects phytoplankton diversity. These processes collectively form
feedbacks between phytoplankton diversity and dynamics of plankton and
nutrient. More importantly, nutrient-diversity feedback introduced additional
temporal variabilities in community biomass, which apparently implies a
destabilizing effect of phytoplankton diversity on ecosystem. However, the
variabilities made ecosystems more robust against extinction of plankton
because increasing phytoplankton diversity facilitates resource consumptions
when consumers prone to extinct; while, reducing diversity weakens
destabilizing dynamics caused by over-growth. Our results suggest the presence
of a novel stabilizing effect of biodiversity acting through nutrient-diversity
feedback, being independent of portfolio mechanisms.


  We define a novel quantitative strategy inspired by the ecological notion of
nestedness to single out the scale at which innovation complexity emerges from
the aggregation of specialized building blocks. Our analysis not only suggests
that the innovation space can be interpreted as a natural system in which
advantageous capabilities are selected by evolutionary pressure, but also that
the emerging structure of capabilities is not independent of the scale of
observation at which they are observed. Expanding on this insight allows us to
understand whether the capabilities characterizing the innovation space at a
given scale are compatible with a complex evolutionary dynamics or, rather, a
set of essentially independent activities allowing to reduce the system at that
scale to a set of disjoint non interacting sub-systems. This yields a measure
of the innovation complexity of the system, i.e. of the degree of
interdependence between the sets of capabilities underlying the system's
building blocks.


  Identifying the secondary structure of an RNA is crucial for understanding
its diverse regulatory functions. This paper focuses on how to enhance target
identification in a Boltzmann ensemble of structures via chemical probing data.
We employ an information-theoretic approach to solve the problem, via
considering a variant of the R\'{e}nyi-Ulam game. Our framework is centered
around the ensemble tree, a hierarchical bi-partition of the input ensemble,
that is constructed by recursively querying about whether or not a base pair of
maximum information entropy is contained in the target. These queries are
answered via relating local with global probing data, employing the modularity
in RNA secondary structures. We present that leaves of the tree are comprised
of sub-samples exhibiting a distinguished structure with high probability. In
particular, for a Boltzmann ensemble incorporating probing data, which is well
established in the literature, the probability of our framework correctly
identifying the target in the leaf is greater than $90\%$.


  The paper studies homeostatic ion channel trafficking in neurons. We derive a
nonlinear closed-loop model that captures active transport with degradation,
channel insertion, average membrane potential activity, and integral control.
We study the model via dominance theory and differential dissipativity to show
when steady regulation gives way to pathological oscillations. We provide
quantitative results on the robustness of the closed loop behavior to static
and dynamic uncertainties, which allows us to understand how cell growth
interacts with ion channel regulation.


  Computing the stationary distributions of a continuous-time Markov chain
(CTMC) involves solving a set of linear equations. In most cases of interest,
the number of equations is infinite or too large, and the equations cannot be
solved analytically or numerically. Several approximation schemes overcome this
issue by truncating the state space to a manageable size. In this review, we
first give a comprehensive theoretical account of the stationary distributions
and their relation to the long-term behaviour of CTMCs that is readily
accessible to non-experts and free of irreducibility assumptions made in
standard texts. We then review truncation-based approximation schemes for CTMCs
with infinite state spaces paying particular attention to the schemes'
convergence and the errors they introduce, and we illustrate their performance
with an example of a stochastic reaction network of relevance in biology and
chemistry. We conclude by discussing computational trade-offs associated with
error control and several open questions.


  We investigate the outcome of generalised Lotka-Volterra dynamics of
ecological communities with random interaction coefficients and non-linear
feedback. We show in simulations that the saturation of non-linear feedback
stabilises the dynamics. This is confirmed in an analytical
generating-functional approach to generalised Lotka-Volterra equations with
piecewise linear saturating response. For such systems we are able to derive
self-consistent relations governing the stable fixed-point phase, and to carry
out a linear stability analysis to predict the onset of unstable behaviour. We
investigate in detail the combined effects of the mean, variance and
co-variance of the random interaction coefficients, and the saturation value of
the non-linear response. We find that stability and diversity increases with
the introduction of non-linear feedback, where decreasing the saturation value
has a similar effect to decreasing the co-variance. We also find co-operation
to no longer have a detrimental effect on stability with non-linear feedback,
and the order parameters mean abundance and diversity to be less dependent on
the symmetry of interactions with stronger saturation.


  Collective decision making processes lie at the heart of many social,
political and economic challenges. The classical voter model is a
well-established conceptual model to study such processes. In this work, we
define a new form of adaptive (or co-evolutionary) voter model posed on a
simplicial complex, i.e., on a certain class of hypernetworks/hypergraphs. We
use the persuasion rule along edges of the classical voter model and the
recently studied re-wiring rule of edges towards like-minded nodes, and
introduce a new peer pressure rule applied to three nodes connected via a
2-simplex. This simplicial adaptive voter model is studied via numerical
simulation. We show that adding the effect of peer pressure to an adaptive
voter model leaves its fragmentation transition, i.e., the transition upon
varying the re-wiring rate from a single majority state into to a fragmented
state of two different opinion subgraphs, intact. Yet, above and below the
fragmentation transition, we observe that the peer pressure has substantial
quantitative effects. It accelerates the transition to a single-opinion state
below the transition and also speeds up the system dynamics towards
fragmentation above the transition. Furthermore, we quantify that there is a
multiscale hierarchy in the model leading to the depletion of 2-simplices,
before the depletion of active edges. This leads to the conjecture that many
other dynamic network models on simplicial complexes may show a similar
behaviour with respect to the sequential evolution of simplicies of different
dimensions.


  Grid diagrams with their relatively simple mathematical formalism provide a
convenient way to generate and model projections of various knots. It has been
an open question whether these 2D diagrams can be used to model a complex 3D
process such as the topoisomerase-mediated preferential unknotting of DNA
molecules. We model here topoisomerase-mediated passages of double-stranded DNA
segments through each other using the formalism of grid diagrams. We show that
this grid diagram-based modelling approach captures the essence of the
preferential unknotting mechanism, based on topoisomerase selectivity of hooked
DNA juxtapositions as the sites of intersegmental passages. We show that grid
diagram-based approach provide an important, new and computationally convenient
framework for investigating entanglement in biopolymers.


  Recent advances in high-throughput genomic technologies coupled with
exponential increases in computer processing and memory have allowed us to
interrogate the complex aberrant molecular underpinnings of human disease from
a genome-wide perspective. While the deluge of genomic information is expected
to increase, a bottleneck in conventional high-performance computing is rapidly
approaching. Inspired in part by recent advances in physical quantum
processors, we evaluated several unconventional machine learning (ML)
strategies on actual human tumor data. Here we show for the first time the
efficacy of multiple annealing-based ML algorithms for classification of
high-dimensional, multi-omics human cancer data from the Cancer Genome Atlas.
To assess algorithm performance, we compared these classifiers to a variety of
standard ML methods. Our results indicate the feasibility of using
annealing-based ML to provide competitive classification of human cancer types
and associated molecular subtypes and superior performance with smaller
training datasets, thus providing compelling empirical evidence for the
potential future application of unconventional computing architectures in the
biomedical sciences.


  Opinion formation is a process with strong implications for public policy. In
controversial debates with large consequences, the public opinion is often
trapped in a fifty-fifty stalemate, jeopardizing broadly accepted political
decisions. Emergent effects from millions of private discussions make it hard
to understand or influence this kind of opinion dynamics. Here we demonstrate
that repulsion from opinions favors fifty-fifty stalemates. We study a voter
model where agents can have two opinions or an undecided state in-between. In
pairwise discussions, undecided agents can be convinced or repelled from the
opinion expressed by another agent. If repulsion happens in at least one of
four cases, as in controversial debates, the frequencies of both opinions
equalize. Further we include transitions of decided agents to the undecided
state. If that happens often, the share of undecided agents becomes large, as
can be measured with the share of undecided answers in polls.


  In this letter, we mainly consider an MSM (men have sex with men) network to
analysis how monogamy behavior can influence the transmission of HIV. By
calculating and analyzing the basic reproductive number of that network, we
find the condition for when the monogamy rate can have a positive influence on
controlling the transmission of HIV. Numerical simulations are also done to
illustrate that monogamy can influence the transmission process of HIV.


  Almost all fields of science rely upon statistical inference to estimate
unknown parameters in theoretical and computational models. While the
performance of modern computer hardware continues to grow, the computational
requirements for the simulation of models are growing even faster. This is
largely due to the increase in model complexity, often including stochastic
dynamics, that is necessary to describe and characterize phenomena observed
using modern, high resolution, experimental techniques. Such models are rarely
analytically tractable, meaning that extremely large numbers of stochastic
simulations are required for parameter inference. In such cases, parameter
inference can be practically impossible. In this work, we present new
computational Bayesian techniques that accelerate inference for expensive
stochastic models by using computationally inexpensive approximations to inform
feasible regions in parameter space, and through learning transforms that
adjust the biased approximate inferences to closer represent the correct
inferences under the expensive stochastic model. Using topical examples from
ecology and cell biology, we demonstrate a speed improvement of an order of
magnitude without any loss in accuracy. This represents a substantial
improvement over current state-of-the-art methods for Bayesian computations
when appropriate model approximations are available.


  Modern digital electronics support remarkably reliable computing, especially
given the challenge of controlling nanoscale logical components that interact
in fluctuating environments. However, we demonstrate that the high-reliability
limit is subject to a fundamental error-energy-efficiency tradeoff that arises
from time-symmetric control: Requiring a low probability of error causes energy
consumption to diverge as logarithm of the inverse error rate for nonreciprocal
logical transitions. The reciprocity (self-invertibility) of a computation is a
stricter condition for thermodynamic efficiency than logical reversibility
(invertibility), the latter being the root of Landauer's work bound on erasing
information. Beyond engineered computation, the results identify a generic
error-dissipation tradeoff in steady-state transformations of genetic
information carried out by biological organisms. The lesson is that computation
under time-symmetric control cannot reach, and is often far above, the Landauer
limit. In this way, time-asymmetry becomes a design principle for
thermodynamically efficient computing.


  Objective: Determining the structure of mortality from causes and
sociodemographic characteristics, in the municipality of San Francisco,
Antioquia, 2001-2010. Methodology: Quantitative descriptive study with
retrospective longitudinal data obtained from secondary source of death events
through databases in electronic media supplied by the DANE. A description of
the sociodemographic variables was performed by groups of cause of death, the
life table by sex and years of potential life lost (APVP) for each year group
and cause of death was calculated. Results: External causes and assaults,
homicides, as the main cause of death occurs during the decade of study, and
especially in men, which had higher mortality rates, more likely to die and
less life expectancy during the period. On average men and external causes
showed a higher number of potential years of life lost for the years 2001-2010
life. Conclusions: Men have higher mortality rates over the decade, as external
causes, assaults, murders, and a greater proportion of young men. The causes of
death that bring more potential years of life lost are external causes. Life
expectancy at birth and throughout the decade is greater for women than it is
therefore essential for men, that the municipality analyze the current
situation of these causes of death, and to be carried out public policies that
contribute to its decline and hence to improve the life expectancy of the
population.


  Urea at sufficiently high concentration unfolds the secondary structure of
proteins leading to denaturation. In contrast, Choline Chloride (ChCl) and
urea, in 1:2 molar ratio form a deep eutectic mixture, a liquid at room
temperature and protect proteins from denaturation. In order to get a
microscopic picture of this phenomenon, we perform extensive all-atom molecular
dynamics simulations on a model protein HP-36. Based on our calculation of
Kirkwood-Buff integrals, we analyze the relative accumulation of these
osmolytes around the protein. Additional insights are drawn from the analyses
of the translational and rotational dynamics of solvent molecules and also from
the hydrogen bond auto-correlation functions. In the presence of urea, water
shows slow subdiffusive dynamics around the protein backbone as a consequence
of stronger interaction of water molecules with the backbone atoms. Urea also
shows subdiffusive motion. Addition of choline further slows down the dynamics
of urea restricting its inclusion around the protein backbone. Adding to this,
choline molecules in the first solvation shell of the protein shows the
strongest subdiffusive behavior. In other words, ChCl behaves as a nano-crowder
by excluding urea from the protein backbone and thereby slowing down the
dynamics of the water layer around the protein. This prevents the protein from
denaturation and makes it structurally rigid which in turn is supported by the
smaller radius of gyration and root mean square deviation values of HP-36 when
ChCl is present in the system.


  We analyze from basic physical considerations the Darwinian competition for
reproduction (evolutionary dynamics) of strategists in a Public Goods Game, the
archetype for $n$-agent (group) economical and biological interactions. In the
proposed setup, the population is organized into groups, being the individual
fitness linked to the group performance, while the evolutionary dynamics takes
place globally. Taking advantage of (groups) permutation symmetry, the
nonlinear analysis of the "mesoscale" Markov phase space for many competing
groups is feasible to a large extent, regarding the expected typicality of
evolutionary histories. These predictions are the basis for a sensible
understanding of the numerical simulation results of the agent (microscopic)
dynamics. Potential implications of these results on model-related issues as,
e.g. group selection, the role of "social norms", or sustainability of common
goods, are highlighted in concise terms, before the conclusion.


  We describe a new and computationally efficient Bayesian methodology for
inferring species trees and demographics from unlinked binary markers.
Likelihood calculations are carried out using diffusion models of allele
frequency dynamics combined with a new algorithm for numerically computing
likelihoods of quantitative traits. The diffusion approach allows for analysis
of datasets containing hundreds or thousands of individuals. The method, which
we call \snapper, has been implemented as part of the Beast2 package. We
introduce the models, the efficient algorithms, and report performance of
\snapper on simulated data sets and on SNP data from rattlesnakes and
freshwater turtles.


  Forest management relies on the evaluation of silviculture practices. The
increase in natural risk due to climate change makes it necessary to consider
evaluation criteria that take natural risk into account. Risk integration in
existing software requires advanced programming skills.We propose a
user-friendly software to simulate even-aged and monospecific forest at the
stand level, in order to evaluate and optimize forest management. The software
gives the possibility to run management scenarii with or without considering
the impact of natural risk. The control variables are the dates and rates of
thinning and the cutting age.The risk model is based on a Poisson processus.
The Faustmann approach, including tree damage risk, is used to evaluate future
benefits, economic or ecosystem services. It relies on the calculation of
expected values, for which a dedicated mathematical development has been done.
The optimized criteria used to evaluate the various scenarii are the Faustmann
value and the Averaged yield value.We illustrate the approach and the software
on two case studies: economic optimization of a beech stand and carbon
sequestration optimization of a pine stand.Software interface makes it easy for
users to write their own (growth-tree damage-economic) models without advanced
programming skills. The possibility to run management scenarii with/without
considering the impact of natural risk may contribute improving silviculture
guidelines and adapting them to climate change. We propose future lines of
research and improvement.


  Neurodegenerative diseases are characterized as the progressive loss of
neural cells, e.g. neurons, glial cells. Ageing, monogenic variations, viral
infections, and many other factors are determined and speculated as causes for
them. While many individual genes, such as APP for Alzheimer disease and HTT
for Huntington disease, and biological pathways are studied for
neurodegenerative diseases, system-wide pathogenesis studies are limited. In
this study, we carried out a meta-analysis of RNA-Seq studies for three
neurodegenerative diseases, namely Alzheimer's disease, Parkinson's disease and
Amyotrophic Lateral Sclerosis (ALS) to minimize the batch effect derived
differences and identify the similarly altered factors among studies. Our main
assumption is that these three diseases share some pathological pathway
pattern. For this purpose, we downloaded publicly available Alzheimer's disease
(84 patients + 33 controls = 117 individuals), Parkinson's disease (28 patients
+ 43 controls = 71 individuals) and ALS (2 studies: 46 patients + 25 control =
71 individuals) RNA-Seq data from Sequence Read Archive (SRA) database. The
significantly differentially expressed genes common to these studies were first
identified and analyzed for the patterns in their pathways and variations. Our
meta-analysis revealed the shared nature of differential gene expression and
mutation load of the cellular heat stress response and GABA synthesis in
neurodegenerative diseases. The downregulated GABA synthesis-related genes
(e.g. GAD1 and GAD2) and the upregulated cellular heat stress response-related
genes (e.g. DNAJB6 and HSP90AA1), in addition to their expression patterns,
contain unique variations in samples from patients with neurodegenerative
diseases. The significance of genes and pathways we identified in this study
corroborated by the recent literature on neurodegenerative diseases.


  Molecular chaperones are ATP-consuming biological machines, which facilitate
the folding of proteins and RNA molecules that are kinetically trapped in
misfolded states for long times. Unassisted folding occurs by the kinetic
partitioning mechanism according to which folding to the native state, with low
probability as well as misfolding to one of the many metastable states, with
high probability, occur rapidly on similar time scales. GroEL is an all-purpose
stochastic machine that assists misfolded substrate proteins (SPs) to fold. The
RNA chaperones (CYT-19) help the folding of ribozymes that readily misfold.
GroEL does not interact with the folded proteins but CYT-19 disrupts both the
folded and misfolded ribozymes. Despite this major difference, the Iterative
Annealing Mechanism (IAM) quantitatively explains all the available
experimental data for assisted folding of proteins and ribozymes. Driven by ATP
binding and hydrolysis and GroES binding, GroEL undergoes a catalytic cycle
during which it samples three allosteric states, referred to as T (apo), R (ATP
bound), and R'' (ADP bound). In accord with the IAM predictions, analyses of
the experimental data shows that the efficiency of the GroEL-GroES machinery
and mutants is determined by the resetting rate $k_{R''\rightarrow T}$, which
is largest for the wild type GroEL. Generalized IAM accurately predicts the
folding kinetics of Tetrahymena ribozyme and its variants. Chaperones maximize
the product of the folding rate and the steady state native state fold by
driving the substrates out of equilibrium. Neither the absolute yield nor the
folding rate is optimized.


  Most real-world networks are incompletely observed. Algorithms that can
accurately predict which links are missing can dramatically speedup the
collection of network data and improve the validity of network models. Many
algorithms now exist for predicting missing links, given a partially observed
network, but it has remained unknown whether a single best predictor exists,
how link predictability varies across methods and networks from different
domains, and how close to optimality current methods are. We answer these
questions by systematically evaluating 203 individual link predictor
algorithms, representing three popular families of methods, applied to a large
corpus of 548 structurally diverse networks from six scientific domains. We
first show that individual algorithms exhibit a broad diversity of prediction
errors, such that no one predictor or family is best, or worst, across all
realistic inputs. We then exploit this diversity via meta-learning to construct
a series of "stacked" models that combine predictors into a single algorithm.
Applied to a broad range of synthetic networks, for which we may analytically
calculate optimal performance, these stacked models achieve optimal or nearly
optimal levels of accuracy. Applied to real-world networks, stacked models are
also superior, but their accuracy varies strongly by domain, suggesting that
link prediction may be fundamentally easier in social networks than in
biological or technological networks. These results indicate that the
state-of-the-art for link prediction comes from combining individual
algorithms, which achieves nearly optimal predictions. We close with a brief
discussion of limitations and opportunities for further improvement of these
results.


  Levan is a bacterial homopolysaccharide, which consists of beta-2,6 linked
beta-D-fructose monomers. Because of its structural properties and its health
promoting effects, levan is a promising functional ingredient for the food,
cosmetic and pharma industry. The properties of levan have been reported to be
linked to its molecular weight. For a better understanding of how its molecular
weight determines its polymer conformation in aqueous solution, levan produced
by the food grade acetic acid bacterium Gluconobacter albidus TMW 2.1191 was
analysed over a broad molecular weight range using dynamic and static light
scattering and viscometry. Levan, with low molecular weight, exhibited a
compact random coil structure. As the molecular weight increased, the structure
transformed into a compact non-drained sphere. The density of the sphere
continued to increase with increasing molecular weight. This resulted in a
negative exponent in the Mark-Houwink-Sakurada Plot. For the first time, an
increase in molecular density with increasing molecular weight, as determined
by a negative Mark-Houwink-Sakurada exponent, could be shown for biopolymers.
Our results reveal the unique properties of high-molecular weight levan and
indicate the need of further systematic studies on the structure-function
relationship of levans for their targeted use in food, cosmetic and
pharmaceutical applications.


  We present the performances of our mathematical deep learning (MathDL) models
for D3R Grand Challenge 4 (GC4). This challenge involves pose prediction,
affinity ranking, and free energy estimation for beta secretase 1 (BACE) as
well as affinity ranking and free energy estimation for Cathepsin S (CatS). We
have developed advanced mathematics, namely differential geometry, algebraic
graph, and/or algebraic topology, to accurately and efficiently encode high
dimensional physical/chemical interactions into scalable low-dimensional
rotational and translational invariant representations. These representations
are integrated with deep learning models, such as generative adversarial
networks (GAN) and convolutional neural networks (CNN) for pose prediction and
energy evaluation, respectively. Overall, our MathDL models achieved the top
place in pose prediction for BACE ligands in Stage 1a. Moreover, our
submissions obtained the highest Spearman correlation coefficient on the
affinity ranking of 460 CatS compounds, and the smallest centered root mean
square error on the free energy set of 39 CatS molecules. It is worthy to
mention that our method for docking pose predictions has significantly improved
from our previous ones.


  Climate change impacts could cause progressive decrease of crop quality and
yield, up to harvest failures. In particular, heat waves and other climate
extremes can lead to localized food shortages and even threaten food security
of communities worldwide. In this study, we apply a deep learning architecture
for high resolution forecasting (300 m, 10 days) of the Leaf Area Index (LAI),
whose dynamics has been widely used to model the growth phase of crops and
impact of heat waves. LAI models can be computed at 0.1 degree spatial
resolution with an auto regressive component adjusted with weather conditions,
validated with remote sensing measurements. However model actionability is poor
in regions of varying terrain morphology at this scale (about 8 km at the Alps
latitude). Our deep learning model aims instead at forecasting LAI by training
multiscale multitemporal (MSMT) data from the Copernicus Global Land Service
(CGLS) project for all Europe at 300m resolution and medium-resolution
historical weather data. Further, the deep learning model inputs integrate
high-resolution land surface features, known to improve forecasts of
agricultural productivity. The historical weather data are then replaced with
forecast values to predict LAI values at 10 day horizon on Europe. We propose
the MSMT model to develop a high resolution crop-specific warning system for
mitigating damage due to heat waves and other extreme events.


  Environmental variation can play an important role in ecological competition
by influencing the relative advantage between competing species. Here, we
consider such effects by extending a classical, competitive Moran model to
incorporate an environment that fluctuates periodically in time. We adapt
methods from work on these classical models to investigate the effects of the
magnitude and frequency of environmental fluctuations on two important
population statistics: the probability of fixation and the mean time to
fixation. In particular, we find that for small frequencies, the system behaves
similar to a system with a constant fitness difference between the two species,
and for large frequencies, the system behaves similar to a neutrally
competitive model. Most interestingly, the system exhibits nontrivial behavior
for intermediate frequencies. We conclude by showing that our results agree
quite well with recent theoretical work on competitive models with a
stochastically changing environment, and discuss how the methods we develop
ease the mathematical analysis required to study such models.


  An ideal signature material and method, which can be used to prove the
authenticity of a physical item and against forgery, should be immune to the
fast developments in digital and engineering technologies. Herein, the design
of signature material combining cryptographic hash function and DNA
steganography is proposed. The encrypting materials are used to construct a
series of time-stamped records (blockchain) associated with published hash
values, while each DNA-encrypted block is associated with a set of DNA keys.
The decrypted DNA information, as digital keys, can be validated through a hash
function to compare with the published hash values. The blocks can also be
cross-referenced among different related signatures. While both digital
cryptography and DNA steganography can have large key size, automated brutal
force search is far more labor intensive and expensive for DNA steganography
with wet lab experiments, as compared to its digital counterpart. Moreover, the
time-stamped blockchain structure permits the incorporation of new
cryptographic functions and DNA steganographies over time, thus can evolve over
time without losing the continuous history line.


  We analyse the Hamiltonian structure of a system of first-order ordinary
differential equations used for modeling the interaction of an oncolytic virus
with a tumour cell population. The analysis is based on the existence of a
Jacobi Last Multiplier for the system and a time dependent first integral. For
suitable conditions on the model parameters this allows for the reduction of
the problem to a planar system of equations for which the time dependent
Hamiltonian flows are described. The geometry of the Hamiltonian flows are
finally investigated using the symplectic and cosymplectic methods.


  Recent studies classify the topology of proteins by analysing the
distribution of their projections using knotoids. The approximation of this
distribution depends on the number of projection directions that are sampled.
Here we investigate the relation between knotoids differing only by small
perturbations of the direction of projection. Since such knotoids are connected
by at most a single forbidden move, we characterise forbidden moves in terms of
equivariant band attachment between strongly invertible knots and of strand
passages between $\theta$-curves. This allows for the determination of the
optimal sample size needed to produce a well approximated knotoid distribution.
Based on that and on topological properties of the distribution, we propose a
numerical measure for the determination of deeply knotted proteins that does
not require the computationally expensive method of subchain analysis.


  Human traffic along roads can be a major vector for infectious diseases and
invasive species. Though most road traffic is local, a small number of
long-distance trips can suffice to move an invasion or disease front forward.
Therefore, understanding how many agents travel over long distances and which
routes they choose is key to successful management of diseases and invasions.
Stochastic gravity models have been used to estimate the distribution of trips
between origins and destinations of agents. However, in large-scale systems it
is hard to collect the data required to fit these models, as the number of
long-distance travellers is small, and origins and destinations can have
multiple access points. Therefore, gravity models often provide only relative
measures of the agent flow. Furthermore, gravity models yield no insights into
which roads agents use. We resolve these issues by combining a stochastic
gravity model with a stochastic route choice model. Our hybrid model can be
fitted to survey data collected at roads that are used by many long-distance
travellers. This decreases the sampling effort, allows us to obtain absolute
predictions of both vector pressure and pathways, and permits rigorous model
validation. After introducing our approach in general terms, we demonstrate its
benefits by applying it to the potential invasion of zebra and quagga mussels
(Dreissena spp.) to the Canadian province British Columbia. The model yields an
R-squared value of 0.73 for variance-corrected agent counts at survey
locations.


  The stability of social relationships is important to animals living in
groups, and social network analysis provides a powerful tool to help
characterize and understand their (in)stability and the consequences at the
group level. However, the use of dynamic social networks is still limited in
this context as it requires long-term social data and new analytical tools.
Here, we study the dynamic evolution of a group of 29 Guinea baboons using a
dataset of automatically collected cognitive tests comprising more than 16M
records collected over 3 years. We first built a monthly aggregated temporal
network describing the baboon's co-presence in the cognitive testing booths. We
then used a null model, considering the heterogeneity in the baboons' activity,
to define both positive (association) and negative (avoidance) monthly
networks. We tested social balance theory by combining these positive and
negative social networks. The results showed that the networks were
structurally balanced and that newly created edges also tended to preserve
social balance. We investigated several network metrics to gain insights into
the individual level and group level social networks long-term temporal
evolution. A measure of similarity between successive monthly networks was able
to pinpoint periods of stability and instability and to show how some baboons'
ego-networks remained stable while others changed radically. Our study confirms
the prediction of social balance theory but also shows that large fluctuations
in the numbers of triads may limit its applicability to study the dynamic
evolution of animal social networks. In contrast, the use of the similarity
measure proved to be very versatile and sensitive in detecting relationships'
(in)stabilities at different levels. The changes we identified can be linked,
at least in some cases, to females changing primary male, as observed in the
wild.


  State-space models (SSMs) are a popular tool for modeling animal abundances.
Inference difficulties for simple linear SSMs are well known, particularly in
relation to simultaneous estimation of process and observation variances.
Several remedies to overcome estimation problems have been studied for
relatively simple SSMs, but whether these challenges and proposed remedies
apply for nonlinear stage-structured SSMs, an important class of ecological
models, is less well understood. Here we identify improvements for inference
about nonlinear stage-structured SSMs fit with biased sequential life stage
data. Theoretical analyses indicate parameter identifiability requires
covariates in the state processes. Simulation studies show that plugging in
externally estimated observation variances, as opposed to jointly estimating
them with other parameters, reduces bias and standard error of estimates. In
contrast to previous results for simple linear SSMs, strong confounding between
jointly estimated process and observation variance parameters was not found in
the models explored here. However, when observation variance was also estimated
in the motivating case study, the resulting process variance estimates were
implausibly low (near-zero). As SSMs are used in increasingly complex ways,
understanding when inference can be expected to be successful, and what aids
it, becomes more important. Our study illustrates (i) the need for relevant
process covariates and (ii) the benefits of using externally estimated
observation variances for inference for nonlinear stage-structured SSMs.


  Biomolecular machines are protein complexes that convert between different
forms of free energy. They are utilized in nature to accomplish many cellular
tasks. As isothermal nonequilibrium stochastic objects at low Reynolds number,
they face a distinct set of challenges compared to more familiar
human-engineered macroscopic machines. Here we review central questions in
their performance as free energy transducers, outline theoretical and modeling
approaches to understand these questions, identify both physical limits on
their operational characteristics and design principles for improving
performance, and discuss emerging areas of research.


  Several processes in the cell, such as gene regulation, start when key
proteins recognise and bind to short DNA sequences. However, as these sequences
can be hundreds of million times shorter than the genome, they are hard to find
by simple diffusion: diffusion-limited association rates may underestimate
$in~vitro$ measurements up to several orders of magnitude. Moreover, the rates
increase if the DNA is coiled rather than straight. Here we model how this
works $in~vivo$ in mammalian cells. We use chromatin-chromatin contact data
from state-of-the-art Hi-C experiments to map the protein target-search onto a
network problem. The nodes represent a DNA segment and the weight of the links
is proportional to measured contact probabilities. We then put forward a master
equation for the density of searching protein that allows us to calculate the
association rates across the genome analytically. For segments where the rates
are high, we find that they are enriched with active genes and have high RNA
expression levels. This paper suggests that the DNA's 3D conformation is
important for protein search times $in~vivo$ and offers a method to interpret
protein-binding profiles in eukaryotes that cannot be explained by the DNA
sequence itself.


  In the past few decades, the frequency of pandemics has been increased due to
the growth of urbanization and mobility among countries. Since a disease
spreading in one country could become a pandemic with a potential worldwide
humanitarian and economic impact, it is important to develop models to estimate
the probability of a worldwide pandemic. In this paper, we propose a model of
disease spreading in a structural modular complex network (having communities)
and study how the number of bridge nodes $n$ that connect communities affects
disease spread. We find that our model can be described at a global scale as an
infectious transmission process between communities with global infectious and
recovery time distributions that depend on the internal structure of each
community and $n$. We find that near the critical point as $n$ increases, the
disease reaches most of the communities, but each community has only a small
fraction of recovered nodes. In addition, we obtain that in the limit $n \to
\infty$, the probability of a pandemic increases abruptly at the critical
point. This scenario could make the decision on whether to launch a pandemic
alert or not more difficult. Finally, we show that link percolation theory can
be used at a global scale to estimate the probability of a pandemic since the
global transmissibility between communities has a weak dependence on the global
recovery time.


  Noise-induced oscillations in individual cells are usually characterized by a
non-monotonic power spectrum with an oscillatory autocorrelation function. Here
we develop an analytical approach of stochastic oscillations in a minimal
stochastic gene expression model including promoter state switching, protein
synthesis and degradation, as well as a genetic feedback loop. The
autocorrelated function, power spectrum, and probability distribution of
protein concentration fluctuations are computed in closed form. Using the
solvable model, we illustrate oscillations as a stochastic circulation along a
hysteresis loop. A triphasic stochastic bifurcation upon the increasing
strength of negative feedback is observed, which reveals how stochastic bursts
evolve into stochastic oscillations. Translational bursting is found to enhance
the robustness and broaden the region of stochastic oscillations in our model.
These results provide deeper insights into R. Thomas' two conjectures for
single-cell stochastic gene expression kinetics.


  The timescales of many physical, chemical, and biological processes are
determined by first passage times (FPTs) of diffusion. The overwhelming
majority of FPT research studies the time it takes a single diffusive searcher
to find a target. However, the more relevant quantity in many systems is the
time it takes the fastest searcher to find a target from a large group of
searchers. This fastest FPT depends on extremely rare events and has a
drastically faster timescale than the FPT of a given single searcher. In this
work, we prove a simple explicit formula for every moment of the fastest FPT.
The formula is remarkably universal, as it holds for $d$-dimensional diffusion
processes (i) with general space-dependent diffusivities and force fields, (ii)
on Riemannian manifolds, (iii) in the presence of reflecting obstacles, and
(iv) with partially absorbing targets. Our results rigorously confirm,
generalize, correct, and unify various conjectures and heuristics about the
fastest FPT.


  Chen et al. solved the crystal structure of bovine DHX36 bound to a DNA with
a G-quadruplex (G4) and a single-stranded DNA segment. They believed that the
mechanism they proposed may represent a general model for describing how a
G4-unfolding helicase recognizes and unfolds G4 DNA. Their conclusion is
interesting, however, we noticed that their linear DNA substrate (DNAMyc) that
harbors a Myc-promoter-derived G4-forming sequence was directly used without
pre-folding. This raises the question whether the structure they obtained
really reflects DHX36-mediated G4 recognition and unfolding, or just only
represents a DHX36-binding-induced quasi-folded G4 structure. By a combination
of polymerase extension, DMS footprinting, stopped-flow, and smFRET assays, we
obtained clear evidences that do not support their ATP-independent one-base
translocation structural model. We further revealed that the oscillation of
FRET signal they observed should correspond to a repetitive G4 binding, but not
unfolding, by DHX36.


  Identifying important nodes for disease spreading is a central topic in
network epidemiology. We investigate how well the position of a node,
characterized by standard network measures, can predict its epidemiological
importance in any graph of a given number of nodes. This is in contrast to
other studies that deal with the easier prediction problem of ranking nodes by
their epidemic importance in given graphs. As a benchmark for epidemic
importance, we calculate the exact expected outbreak size given a node as the
source. We study exhaustively all graphs of a given size, so do not restrict
ourselves to certain generative models for graphs, nor to graph data sets. Due
to the large number of possible nonisomorphic graphs of a fixed size, we are
limited to 10-node graphs. We find that combinations of two or more
centralities are predictive ($R^2$ scores of 0.91 or higher) even for the most
difficult parameter values of the epidemic simulation. Typically, these
successful combinations include one normalized spectral centralities (such as
PageRank or Katz centrality) and one measure that is sensitive to the number of
edges in the graph.


  This paper starts from the Fisher-Kolmogorov-Petrovskii-Piskunov equation to
model diffusive populations. The main result, according to this model, is that
two connected patches in a system do not always contribute to each other.
Specifically, inserting a large fragment next to a small one is always positive
for life inside the small patch, while inserting a very small patch next to a
large one can be negative for life inside the large fragment. This result,
obtained to homogeneously fragmented regions, is possible from the general case
expression for the minimum sizes in a system of two patches. This expression by
itself is an interesting result, because it allows the study of other
characteristics not included in the present work.


  Prosocial behaviors are encountered in the donation game, the prisoner's
dilemma, relaxed social dilemmas, and public goods games. Many studies assume
that the population structure is homogeneous, meaning all individuals have the
same number of interaction partners, or that the social good is of one
particular type. Here, we explore general evolutionary dynamics for arbitrary
spatial structures and social goods. We find that heterogeneous networks,
wherein some individuals have many more interaction partners than others, can
enhance the evolution of prosocial behaviors. However, they often accumulate
most of the benefits in the hands of a few highly-connected individuals, while
many others receive low or negative payoff. Surprisingly, selection can favor
producers of social goods even if the total costs exceed the total benefits. In
summary, heterogeneous structures have the ability to strongly promote the
emergence of prosocial behaviors, but they also create the possibility of
generating large inequality.


  The distribution of living (rose Bengal stained) Reophax pyriformis
Yamashita, Sousa and Kaminski, n. sp., an agglutinated benthic foraminiferal
species, was analyzed in the area of the Campos Basin (southeastern Brazilian
continental margin). The study is based on 34 oceanographic stations (54
samples), located between 400 m and 3,000 m water depth. The distribution of
living Reophax pyriformis n.sp. density is compared to sedimentological
parameters, such as total organic carbon, total nitrogen, calcium carbonate,
phytopigment, lipids biomarkers (sterols, fatty acids and n-alcohols), total
lipids, and bacterial biomass, as well as the particulate organic matter flux
to the seafloor. This species was found in the range of 1,000--1,200 m water
depth, with an average living depth of 1.52 cm in the sediment. The
distribution of this species seems to be related to lipid biomarkers
(allochthone and/or continental derivative, zooplankton and/or fauna,
phytoplankton or primary producers) and total organic carbon under the
influence of Intermediate Western Boundary Current conditions.


  In the context of pharmacogenomics, whole genome sequencing provides a
powerful approach for identifying correlations between response variability to
specific drugs and genomic polymorphisms in a population, in an unbiased
manner. In this study, we employed whole genome sequencing of DNA samples from
patients showing extreme response (n=72) and non-response (n=27) to the
antiepileptic drug levetiracetam, in order to identify genomic variants that
underlie response to the drug. Although no common SNP (MAF>5%) crossed the
conventional genome-wide significance threshold of 5e-8, we found common
polymorphisms in genes SPNS3, HDC, MDGA2, NSG1 and RASGEF1C, which collectively
predict clinical response to levetiracetam in our cohort with ~91% predictive
accuracy. Among these genes, HDC, NSG1, MDGA2 and RASGEF1C are potentially
implicated in synaptic neurotransmission, while SPNS3 is an atypical solute
carrier transporter homologous to SV2A, the known molecular target of
levetiracetam. Furthermore, we performed gene- and pathway-based statistical
analysis on sets of rare and low-frequency variants (MAF<5%) and we identified
associations between the following genes or pathways and response to
levetiracetam: a) genes PRKCB and DLG2, which are involved in glutamatergic
neurotransmission, a known target of anticonvulsants, including levetiracetam;
b) genes FILIP1 and SEMA6D, which are involved in axon guidance and modelling
of neural connections; and c) pathways with a role in synaptic
neurotransmission, such as WNT5A-dependent internalization of FZD4 and
disinhibition of SNARE formation. In summary, our approach to utilise whole
genome sequencing on subjects with extreme response phenotypes is a feasible
route to generate plausible hypotheses for investigating the genetic factors
underlying drug response variability in cases of pharmaco-resistant epilepsy.


  Cells use genetic switches to shift between alternate stable gene expression
states, e.g., to adapt to new environments or to follow a developmental
pathway. Conceptually, these stable phenotypes can be considered as attractive
states on an epigenetic landscape with phenotypic changes being transitions
between states. Measuring these transitions is challenging because they are
both very rare in the absence of appropriate signals and very fast. As such, it
has proven difficult to experimentally map the epigenetic landscapes that are
widely believed to underly developmental networks. Here, we introduce a new
nonequilibrium perturbation method to help reconstruct a regulatory network's
epigenetic landscape. We derive the mathematical theory needed and then use the
method on simulated data to reconstruct the landscapes. Our results show that
with a relatively small number of perturbation experiments it is possible to
recover an accurate representation of the true epigenetic landscape. We propose
that our theory provides a general method by which epigenetic landscapes can be
studied. Finally, our theory suggests that the total perturbation impulse
required to induce a switch between metastable states is a fundamental quantity
in developmental dynamics.


  The analytic derivation of a more general model of survival-mortality and the
estimation of a parameter bx related to the Healthy Life Years Lost (HLYL) is
followed with the formulation of a computer program providing results similar
to those of the World Health Organization for the Healthy Life Expectancy
(HALE) and the corresponding HLYL estimates. This program is an extension of
the classical life table including more columns to estimate the cumulative
mortality, the average mortality, the person life years lost and finally the
HLYL parameter bx. Evenmore, a further extension of the Excel program based on
the Sullivan method provides estimates of the Healthy Life Expectancy at every
year of the lifespan for five different types of estimates that are the Direct,
WHO, Eurostat, Equal and Other. Estimates for several countries are presented.
It is also presented a methodology and a program to bridge the gap between the
World Health Organization (HALE) and Eurostat (HLE) healthy life expectancy
estimates. The latest version of this program (SKI-6) appear in the
Demographics2020 website.