  In this paper, we obtain an asymptotic formula for the persistence
probability in the positive real line of a random polynomial arising from
evolutionary game theory. It corresponds to the probability that a multi-player
two-strategy random evolutionary game has no internal equilibria. The key
ingredient is to approximate the sequence of random polynomials indexed by
their degrees by an appropriate centered stationary Gaussian process.


  Background: Since the invention of next-generation RNA sequencing (RNA-seq)
technologies, they have become a powerful tool to study the presence and
quantity of RNA molecules in biological samples and have revolutionized
transcriptomic studies. The analysis of RNA-seq data at four different levels
(samples, genes, transcripts, and exons) involve multiple statistical and
computational questions, some of which remain challenging up to date.
  Results: We review RNA-seq analysis tools at the sample, gene, transcript,
and exon levels from a statistical perspective. We also highlight the
biological and statistical questions of most practical considerations.
  Conclusion: The development of statistical and computational methods for
analyzing RNA- seq data has made significant advances in the past decade.
However, methods developed to answer the same biological question often rely on
diverse statical models and exhibit different performance under different
scenarios. This review discusses and compares multiple commonly used
statistical models regarding their assumptions, in the hope of helping users
select appropriate methods as needed, as well as assisting developers for
future method development.


  The aetiology of polygenic obesity is multifactorial, which indicates that
life-style and environmental factors may influence multiples genes to aggravate
this disorder. Several low-risk single nucleotide polymorphisms (SNPs) have
been associated with BMI. However, identified loci only explain a small
proportion of the variation ob-served for this phenotype. The linear nature of
genome wide association studies (GWAS) used to identify associations between
genetic variants and the phenotype have had limited success in explaining the
heritability variation of BMI and shown low predictive capacity in
classification studies. GWAS ignores the epistatic interactions that less
significant variants have on the phenotypic outcome. In this paper we utilise a
novel deep learning-based methodology to reduce the high dimensional space in
GWAS and find epistatic interactions between SNPs for classification purposes.
SNPs were filtered based on the effects associations have with BMI. Since
Bonferroni adjustment for multiple testing is highly conservative, an important
proportion of SNPs involved in SNP-SNP interactions are ignored. Therefore,
only SNPs with p-values < 1x10-2 were considered for subsequent epistasis
analysis using stacked auto encoders (SAE). This allows the nonlinearity
present in SNP-SNP interactions to be discovered through progressively smaller
hidden layer units and to initialise a multi-layer feedforward artificial
neural network (ANN) classifier. The classifier is fine-tuned to classify
extremely obese and non-obese individuals. The best results were obtained with
2000 compressed units (SE=0.949153, SP=0.933014, Gini=0.949936,
Lo-gloss=0.1956, AUC=0.97497 and MSE=0.054057). Using 50 compressed units it
was possible to achieve (SE=0.785311, SP=0.799043, Gini=0.703566,
Logloss=0.476864, AUC=0.85178 and MSE=0.156315).


  MOTIVATION: Proteins fold into complex structures that are crucial for their
biological functions. Experimental determination of protein structures is
costly and therefore limited to a small fraction of all known proteins. Hence,
different computational structure prediction methods are necessary for the
modelling of the vast majority of all proteins. In most structure prediction
pipelines, the last step is to select the best available model and to estimate
its accuracy. This model quality estimation problem has been growing in
importance during the last decade, and progress is believed to be important for
large scale modelling of proteins. The current generation of model quality
estimation programs performs well at separating incorrect and good models, but
fails to consistently identify the best possible model. State-of-the-art model
quality assessment methods use a combination of features that describe a model
and the agreement of the model with features predicted from the protein
sequence.
  RESULTS: We first introduce a deep neural network architecture to predict
model quality using significantly fewer input features than state-of-the-art
methods. Thereafter, we propose a methodology to train the deep network that
leverages the comparative structure of the problem. We also show the
possibility of applying transfer learning on databases of known protein
structures. We demonstrate its viability by reaching state-of-the-art
performance using only a reduced set of input features and a coarse description
of the models.
  AVAILABILITY: The code will be freely available for download at
github.com/ElofssonLab/ProQ4.


  Bayesian network models are finding success in characterizing
enzyme-catalyzed reactions, slow conformational changes, predicting enzyme
inhibition, and genomics. In this work, we apply them to statistical modeling
of peptides by simultaneously identifying amino acid sequence motifs and using
a motif-based model to clarify the role motifs may play in antimicrobial
activity. We construct models of increasing sophistication, demonstrating how
chemical knowledge of a peptide system may be embedded without requiring new
derivation of model fitting equations after changing model structure. These
models are used to construct classifiers with good performance (94% accuracy,
Matthews correlation coefficient of 0.87) at predicting antimicrobial activity
in peptides, while at the same time being built of interpretable parameters. We
demonstrate use of these models to identify peptides that are potentially both
antimicrobial and antifouling, and show that the background distribution of
amino acids could play a greater role in activity than sequence motifs do. This
provides an advancement in the type of peptide activity modeling that can be
done and the ease in which models can be constructed.


  In mechanical manipulation experiments, a single DNA molecule overwound at
constant force undergoes a discontinuous drop in extension as it buckles and
forms a superhelical loop (a plectoneme). Further overwinding the DNA, we
observe an unanticipated cascade of highly regular discontinuous extension
changes associated with stepwise plectoneme lengthening. This phenomenon is
consistent with a model in which the force-extended DNA forms barriers to
plectoneme lengthening caused by topological writhe. Furthermore, accounting
for writhe in a fluctuating solenoid gives an improved description of the
measured force-dependent effective torsional modulus of DNA, providing a
reliable formula to estimate DNA torque. Our data and model thus provide
context for further measurements and theories that capture the structures and
mechanics of supercoiled biopolymers.


  In this work, we aim to answer the question: what triggers cooperative
behavior in the thermodynamic limit by taking recourse to the Public goods
game. Using the idea of mapping the 1D Ising model Hamiltonian with nearest
neighbor coupling to payoffs in the game theory we calculate the Magnetisation
of the game in the thermodynamic limit. We see a phase transition in the
thermodynamic limit of the two player Public goods game. We observe that
punishment acts as an external field for the two player Public goods game
triggering cooperation or provide strategy, while cost can be a trigger for
suppressing cooperation or free riding. Finally, reward also acts as a trigger
for providing while the role of inverse temperature (fluctuations in choices)
is to introduce randomness in strategic choices.


  Zhou and Smith (2017) investigate different multi-species harvesting
scenarios using a simple Holling-Tanner model. Among these scenarios are two
methods for implementing balanced harvesting, where fishing is distributed
across trophic levels in accordance with their productivity. This note examines
the effects of a different quantitative implementation of balanced harvesting,
where the fishing mortality rate is proportional to the total production rate
of each trophic level. The results show that setting fishing mortality rate to
be proportional to total production rate, rather than to productivity per unit
biomass, better preserves trophic structure and provides a crucial safeguard
for rare and threatened ecological groups. This is a key ingredient of balanced
harvesting if it is to meet its objective of preserving biodiversity.


  Evolutionary dynamics is often viewed as a subtle process of change
accumulation that causes a divergence among organisms and their genomes.
However, this interpretation is an inheritance of a gradualistic view that has
been challenged at the macroevolutionary, ecological, and molecular level.
Actually, when the complex architecture of genotype spaces is taken into
account, the evolutionary dynamics of molecular populations becomes
intrinsically non-uniform, sharing deep qualitative and quantitative
similarities with slowly driven physical systems: non-linear responses
analogous to critical transitions, sudden state changes, or hysteresis, among
others. Furthermore, the phenotypic plasticity inherent to genotypes transforms
classical fitness landscapes into multiscapes where adaptation in response to
an environmental change may be very fast. The quantitative nature of adaptive
molecular processes is deeply dependent on a networks-of-networks multilayered
structure of the map from genotype to function that we begin to unveil.


  Recently, it has become feasible to generate large-scale, multi-tissue gene
expression data, where expression profiles are obtained from multiple tissues
or organs sampled from dozens to hundreds of individuals. When traditional
clustering methods are applied to this type of data, important information is
lost, because they either require all tissues to be analyzed independently,
ignoring dependencies and similarities between tissues, or to merge tissues in
a single, monolithic dataset, ignoring individual characteristics of tissues.
We developed a Bayesian model-based multi-tissue clustering algorithm, revamp,
which can incorporate prior information on physiological tissue similarity, and
which results in a set of clusters, each consisting of a core set of genes
conserved across tissues as well as differential sets of genes specific to one
or more subsets of tissues. Using data from seven vascular and metabolic
tissues from over 100 individuals in the STockholm Atherosclerosis Gene
Expression (STAGE) study, we demonstrate that multi-tissue clusters inferred by
revamp are more enriched for tissue-dependent protein-protein interactions
compared to alternative approaches. We further demonstrate that revamp results
in easily interpretable multi-tissue gene expression associations to key
coronary artery disease processes and clinical phenotypes in the STAGE
individuals. Revamp is implemented in the Lemon-Tree software, available at
https://github.com/eb00/lemon-tree


  Cooperation among self-interested players in a social dilemma is fragile and
easily interrupted by mistakes. In this work, we study the repeated $n$-person
public-goods game and search for a strategy that forms a cooperative Nash
equilibrium in the presence of implementation error with a guarantee that the
resulting payoff will be no less than any of the co-players'. By enumerating
strategic possibilities for $n=3$, we show that such a strategy indeed exists
when its memory length $m$ equals three. It means that a deterministic strategy
can be publicly employed to stabilize cooperation against error with avoiding
the risk of being exploited. We furthermore show that, for general $n$-person
public-goods game, $m \geq n$ is necessary to satisfy the above criteria.


  The distributions of species lifetimes and species in space are related,
since species with good local survival chances have more time to colonize new
habitats and species inhabiting large areas have higher chances to survive
local disturbances. Yet, both distributions have been discussed in mostly
separate communities. Here, we study both patterns simultaneously using a
spatially explicit, evolutionary community assembly approach. We present and
investigate a metacommunity model, consisting of a grid of patches, where each
patch contains a local food web. Species survival depends on predation and
competition interactions, which in turn depend on species body masses as the
key traits. The system evolves due to the migration of species to neighboring
patches, the addition of new species as modifications of existing species, and
local extinction events. The structure of each local food web thus emerges in a
self-organized manner as the highly non-trivial outcome of the relative time
scales of these processes. Our model generates a large variety of complex,
multi-trophic networks and therefore serves as a powerful tool to investigate
ecosystems on long temporal and large spatial scales. We find that the observed
lifetime distributions and species-area relations resemble power laws over
appropriately chosen parameter ranges and thus agree qualitatively with
empirical findings. Moreover, we observe strong finite-size effects, and a
dependence of the relationships on the trophic level of the species. By
comparing our results to simple neutral models found in the literature, we
identify the features that are responsible for the values of the exponents.


  Perfect phylogenies are fundamental in the study of evolutionary trees
because they capture the situation when each evolutionary trait emerges only
once in history; if such events are believed to be rare, then by Occam's Razor
such parsimonious trees are preferable as a hypothesis of evolution. A
classical result states that 2-state characters permit a perfect phylogeny
precisely if each subset of 2 characters permits one. More recently, it was
shown that for 3-state characters the same property holds but for size-3
subsets. A long-standing open problem asked whether such a constant exists for
each number of states. More precisely, it has been conjectured that for any
fixed integer $r$, there exists a constant $f(r)$ such that a set of $r$-state
characters $C$ has a perfect phylogeny if and only if every subset of at most
$f(r)$ characters has a perfect phylogeny. In this paper, we show that this
conjecture is false. In particular, we show that for any constant $t$, there
exists a set $C$ of $8$-state characters such that $C$ has no perfect
phylogeny, but there exists a perfect phylogeny for every subset of $t$
characters. This negative result complements the two negative results
("strikes") of Bodlaender et al. We reflect on the consequences of this third
strike, pointing out that while it does close off some routes for efficient
algorithm development, many others remain open.


  A fifteen state kinetic model for Na+ channel gating that describes the
coupling between three activation sensors, a two-stage fast inactivation
process and slow inactivated states, may be reduced to equations for a six
state system by application of the method of multiple scales. By expressing the
occupation probabilities for closed states and the open state in terms of
activation and fast inactivation variables, and assuming that activation has a
faster relaxation than inactivation and that the activation sensors are
mutually independent, the kinetic equations may be further reduced to rate
equations for activation, and coupled fast and slow inactivation that describe
spike frequency adaptation, a repetitive bursting oscillation in the neural
membrane, and a cardiac action potential with a plateau oscillation. The fast
inactivation rate function is, in general, dependent on the activation variable
m(t) but may be approximated by a voltage-dependent function, and the rate
function for entry into the slow inactivated state is dependent on the fast
inactivation variable.


  Managing blood lipid levels is important for the treatment and prevention of
diabetes, cardiovascular disease, and obesity. An easy-to-use, portable lipid
blood test will accelerate more frequent testing by patients and at-risk
populations. We used smartphone systems that are already familiar to many
people. Because smartphone systems can be carried around everywhere, blood can
be measured easily and frequently. We compared the results of lipid tests with
those of existing clinical diagnostic laboratory methods. We found that
smartphone-based point-of-care lipid blood tests are as accurate as
hospital-grade laboratory tests. Our system will be useful for those who need
to manage blood lipid levels to motivate them to track and control their
behavior.


  Humans interact with each other on a daily basis by developing and
maintaining various social norms and it is critical to form a deeper
understanding of how such norms develop, how they change, and how fast they
change. In this work, we develop an evolutionary game-theoretic model based on
research in cultural psychology that shows that humans in various cultures
differ in their tendencies to conform with those around them. Using this model,
we analyze the evolutionary relationships between the tendency to conform and
how quickly a population reacts when conditions make a change in norm
desirable. Our analysis identifies conditions when a tipping point is reached
in a population, causing norms to change rapidly.


  The yeast Saccharomyces cerevisiae is one of the best characterized
eukaryotic models. The secretory pathway was the first trafficking pathway
clearly understood mainly thanks to the work done in the laboratory of Randy
Schekman in the 1980s. They have isolated yeast sec mutants unable to secrete
an extracellular enzyme and these SEC genes were identified as encoding key
effectors of the secretory machinery. For this work, the 2013 Nobel Prize in
Physiology and Medicine has been awarded to Randy Schekman; the prize is shared
with James Rothman and Thomas S{\"u}dhof. Here, we present the different
trafficking pathways of yeast S. cerevisiae. At the Golgi apparatus newly
synthesized proteins are sorted between those transported to the plasma
membrane (PM), or the external medium, via the exocytosis or secretory pathway
(SEC), and those targeted to the vacuole either through endosomes (vacuolar
protein sorting or VPS pathway) or directly (alkaline phosphatase or ALP
pathway). Plasma membrane proteins can be internalized by endocytosis (END) and
transported to endosomes where they are sorted between those targeted for
vacuolar degradation and those redirected to the Golgi (recycling or RCY
pathway). Studies in yeast S. cerevisiae allowed the identification of most of
the known effectors, protein complexes, and trafficking pathways in eukaryotic
cells, and most of them are conserved among eukaryotes.


  We present the bond-calculus, a process algebra for modelling biological and
chemical systems featuring nonlinear dynamics, multiway interactions, and
dynamic bonding of agents. Mathematical models based on differential equations
have been instrumental in modelling and understanding the dynamics of
biological systems. Quantitative process algebras aim to build higher level
descriptions of biological systems, capturing the agents and interactions
underlying their behaviour, and can be compiled down to a range of lower level
mathematical models. The bond-calculus builds upon the work of Kwiatkowski,
Banks, and Stark's continuous pi-calculus by adding a flexible multiway
communication operation based on affinity patterns and general kinetic laws. We
develop a compositional semantics based on vector fields and linear operators,
which we use to define the time evolution of this system. This enables
simulation and analysis via differential equation generation or stochastic
simulation. Finally, we apply our framework to an existing biological model:
Kuznetsov's classic model of tumour immune interactions.


  Evolution has fascinated quantitative and physical scientists for decades:
how can the random process of mutation, recombination, and duplication of
genetic information generate the diversity of life? What determines the rate of
evolution? Are there quantitative laws that govern and constrain evolution? Is
evolution repeatable or predictable? Historically, the study of evolution
involved classifying and comparing species, typically based on morphology. In
addition to phenotypes on the organismal and molecular scales, we now use
whole-genome sequencing to uncover not only the differences between species but
also to characterize genetic diversity within-species in unprecedented detail.
This diversity can be compared to predictions of quantitative models of
evolutionary dynamics. Here, we review key theoretical models of population
genetics and evolution along with examples of data from lab evolution
experiments, longitudinal sampling of viral populations, microbial communities
and the studies of immune repertoires. In all these systems, evolution is
shaped by often variable biological and physical environments. While these
variable environments can be modeled implicitly in cases such as host-pathogen
co-evolution, the dynamic environment, and emerging ecology often cannot be
ignored. Integrating dynamics on different scales, both in terms of observation
and theoretical models, remains a major challenge towards a better
understanding of evolution.


  Many important analgesics relieve pain by binding to the $\mu$-Opioid
Receptor ($\mu$OR), which makes the $\mu$OR among the most clinically relevant
proteins of the G Protein Coupled Receptor (GPCR) family. Despite previous
studies on the activation pathways of the GPCRs, the mechanism of opiate
binding and the selectivity of $\mu$OR are largely unknown. We performed
extensive molecular dynamics (MD) simulation and analysis to find the selective
allosteric binding sites of the $\mu$OR and the path opiates take to bind to
the orthosteric site. In this study, we predicted that the allosteric site is
responsible for the attraction and selection of opiates. Using Markov state
models and machine learning, we traced the pathway of opiates in binding to the
orthosteric site, the main binding pocket. Our results have important
implications in designing novel analgesics.


  Mutations in the MTM1 gene, encoding the phosphoinositide phosphatase
myotubularin, are responsible for the X-linked centronuclear myopathy (XLCNM)
or X-linked myotubular myopathy (XLMTM). The MTM1 gene was first identified in
1996 and its function as a PtdIns3P and PtdIns(,5)P2 phosphatase was discovered
in 2000. In recent years, very important progress has been made to set up good
models to study MTM1 and the XLCNM disease such as knockout or knockin mice,
the Labrador Retriever dog, the zebrafish and the yeast Saccharomyces
cerevisiae. These helped to better understand the cellular function of MTM1 and
of its four conserved domains: PH-GRAM (Pleckstrin
Homology-Glucosyltransferase, Rab-like GTPase Activator and Myotubularin), RID
(Rac1-Induced recruitment Domain), PTP/DSP (Protein Tyrosine
Phosphatase/Dual-Specificity Phosphatase) and SID (SET-protein Interaction
Domain). This review presents the cellular function of human myotubularin MTM1
and its yeast homolog yeast protein Ymr1, and the role of MTM1 in the
centronuclear myopathy (CNM) disease.


  Pathogen transmission and virulence are main evolutionary variables broadly
assumed to be linked through trade-offs. In well-mixed populations, these
trade-offs are often ascribed to physiological restrictions, while populations
with spatial self-structuring might evolve emergent trade-offs. Here, we
reexamine a model of the latter kind proposed by Ballegooijen and Boerlijst
with the aim of characterising the mechanisms causing the emergence of the
trade-off and its structural robustness. Using invadability criteria, we
establish the conditions under which an evolutionary feedback between
transmission and virulence mediated by pattern formation can poise the system
to a critical boundary separating a disordered state (without emergent
trade-off) from a self-structured phase (where the trade-off emerges), and
analytically calculate the functional shape of the boundary in a certain
approximation. Beyond evolutionary parameters, the success of an invasion
depends on the size and spatial structure of the invading and invaded
populations. Spatial self-structuring is often destroyed when hosts are mobile,
changing the evolutionary dynamics to those of a well-mixed population. In a
metapopulation scenario, the systematic extinction of the pathogen in the
disordered phase may counteract the disruptive effect of host mobility, favour
pattern formation and therefore recover the emergent trade-off.


  Spreading processes have been largely studied in the literature, both
analytically and by means of large-scale numerical simulations. These processes
mainly include the propagation of diseases, rumors and information on top of a
given population. In the last two decades, with the advent of modern network
science, we have witnessed significant advances in this field of research. Here
we review the main theoretical and numerical methods developed for the study of
spreading processes on complex networked systems. Specifically, we formally
define epidemic processes on single and multilayer networks and discuss in
detail the main methods used to perform numerical simulations. Throughout the
review, we classify spreading processes (disease and rumor models) into two
classes according to the nature of time: (i) continuous-time and (ii) cellular
automata approach, where the second one can be further divided into synchronous
and asynchronous updating schemes. Our revision includes the heterogeneous
mean-field, the quenched-mean field, and the pair quenched mean field
approaches, as well as their respective simulation techniques, emphasizing
similarities and differences among the different techniques. The content
presented here offers a whole suite of methods to study epidemic-like processes
in complex networks, both for researchers without previous experience in the
subject and for experts.


  The United States has an opioid epidemic, and Pennsylvania's Allegheny County
is among the worst. This motivates a deeper exploration of what characterizes
the epidemic, such as what are risk factors for people who ultimately overdose
and die due to opioids. We show that some interesting trends and factors can be
identified from openly available autopsy data, and demonstrate the power of
building an interactive data exploration tool for policy makers. However, there
is still a pressing need to incorporate further demographic factors. We show
this by using synthetic Electronic Medical Record (EMR) data to simulate the
predictive power of random forests and neural networks when given additional
loosely correlated features. In addition, we give examples of useful feature
extraction that enable model enhancement without sacrificing privacy.


  A detailed algorithmic explanation is required for how a network of chemical
reactions can generate the sophisticated behavior displayed by living cells.
Though several previous works have shown that reaction networks are
computationally universal and can in principle implement any algorithm, there
is scope for constructions that map well onto biological reality, make
efficient use of the computational potential of the native dynamics of reaction
networks, and make contact with statistical mechanics. We describe a new
reaction network scheme for solving a large class of statistical problems
including the problem of how a cell would infer its environment from
receptor-ligand bindings. Specifically we show how reaction networks can
implement information projection, and consequently a generalized
Expectation-Maximization algorithm, to solve maximum likelihood estimation
problems in partially-observed exponential families on categorical data. Our
scheme can be thought of as an algorithmic interpretation of E. T. Jaynes's
vision of statistical mechanics as statistical inference.


  Asymmetric segregation of key proteins at cell division -- be it a beneficial
or deleterious protein -- is ubiquitous in unicellular organisms and often
considered as an evolved trait to increase fitness in a stressed environment.
Here, we provide a general framework to describe the evolutionary origin of
this asymmetric segregation. We compute the population fitness as a function of
the protein segregation asymmetry $a$, and show that the value of $a$ which
optimizes the population growth manifests a phase transition between symmetric
and asymmetric partitioning phases. Surprisingly, the nature of phase
transition is different for the case of beneficial proteins as opposed to
proteins which decrease the single-cell growth rate. Our study elucidates the
optimization problem faced by evolution in the context of protein segregation,
and motivates further investigation of asymmetric protein segregation in
biological systems.


  We study competition of two non-motile bacterial strains in a
three-dimensional channel numerically, and analyze how their configuration
evolves in space and time. We construct a lattice model that takes into account
self-replication, mutation, and killing of bacteria. When mutation is not
significant, the two strains segregate and form stripe patterns along the
channel. The formed lanes are gradually rearranged, with increasing length
scales in the two-dimensional cross-sectional plane. We characterize it in
terms of coarsening and phase ordering in statistical physics. In particular,
for the simple model without mutation and killing, we find logarithmically slow
coarsening, which is characteristic of the two-dimensional voter model. With
mutation and killing, we find a phase transition from a monopolistic phase, in
which lanes are formed and coarsened until the system is eventually dominated
by one of the two strains, to an equally mixed and disordered phase without
lane structure. Critical behavior at the transition point is also studied and
compared with the generalized voter class and the Ising class. These results
are accounted for by continuum equations, obtained by applying a mean field
approximation along the channel axis. Our findings indicate relevance of
critical coarsening of two-dimensional systems in the problem of bacterial
competition within anisotropic three-dimensional geometry.


  Protein-protein interaction (PPI) network alignment is a canonical operation
to transfer biological knowledge among species. The alignment of PPI-networks
has many applications, such as the prediction of protein function, detection of
conserved network motifs, and the reconstruction of species' phylogenetic
relationships. A good multiple-network alignment (MNA), by considering the data
related to several species, provides a deep understanding of biological
networks and system-level cellular processes. With the massive amounts of
available PPI data and the increasing number of known PPI networks, the problem
of MNA is gaining more attention in the systems-biology studies.
  In this paper, we introduce a new scalable and accurate algorithm, called
MPGM, for aligning multiple networks. The MPGM algorithm has two main steps:
(i) SEEDGENERATION and (ii) MULTIPLEPERCOLATION. In the first step, to generate
an initial set of seed tuples, the SEEDGENERATION algorithm uses only protein
sequence similarities. In the second step, to align remaining unmatched nodes,
the MULTIPLEPERCOLATION algorithm uses network structures and the seed tuples
generated from the first step. We show that, with respect to different
evaluation criteria, MPGM outperforms the other state-of-the-art algorithms. In
addition, we guarantee the performance of MPGM under certain classes of network
models. We introduce a sampling-based stochastic model for generating k
correlated networks. We prove that for this model, if a sufficient number of
seed tuples are available, the MULTIPLEPERCOLATION algorithm correctly aligns
almost all the nodes. Our theoretical results are supported by experimental
evaluations over synthetic networks.


  All living systems can function only far away from equilibrium, and for this
reason chemical kinetic methods are critically important for uncovering the
mechanisms of biological processes. Here we present a new theoretical method of
investigating dynamics of protein-DNA interactions, which govern all major
biological processes. It is based on a first-passage analysis of biochemical
and biophysical transitions, and it provides a fully analytic description of
the processes. Our approach is explained for the case of a single protein
searching for a specific binding site on DNA. In addition, the application of
the method to investigations of the effect of DNA sequence heterogeneity, and
the role multiple targets and traps in the protein search dynamics are
discussed.


  Yuri Manin's approach to Zipf's law (Kolmogorov complexity as energy) is
applied to investigation of biological evolution. Model of constructive
statistical mechanics where complexity is a contribution to energy is proposed
to model genomics. Scaling laws in genomics are discussed in relation to Zipf's
law. This gives a model of Eugene Koonin's Third Evolutionary Synthesis
(physical model which should describe scaling in genomics).


  Telomeres, the nucleoprotein complexes at the termini of linear chromosomes,
are essential for the processes of end replication, end protection, and
chromatin segregation. The Mre11 complex is involved in multiple cellular roles
in DNA repair and structure in the regulation and function of telomere size
homeostasis. In this study, we characterize yeast telomere chromatin structure,
phenotypic heritability, and chromatin segregation in both wild-type [MRE11]
and A470 motif alleles. MRE11 strains confer a telomere size of 300 base pairs
of G+T irregular simple sequence repeats. This DNA and a portion of
subtelomeric DNA is embedded in a telosome: a MNase-resistant non-nucleosomal
particle. Chromatin immunoprecipitation shows a three to four-fold lower
occupancy of Mre11A470T proteins than wild-type proteins in telosomes.
Telosomes containing the Mre11A470T protein confer a greater resistance to
MNase digestion than wild-type telosomes. The integration of a wild-type MRE11
allele into an ectopic locus in the genome of an mre11A470T mutant and the
introduction of an mre11A470T allele at an ectopic site in a wild-type strain
lead to unexpectedly differing results. In each case, the replicated sister
chromatids inherit telosomes containing only the protein encoded by the genomic
mre11 locus, even in the presence of protein encoded by the opposing ectopic
allele. We hypothesize that the telosome segregates by a conservative
mechanism. These data support a mechanism for the linkage between sister
chromatid replication and maintenance of either identical mutant or identical
wild-type telosomes after replication of sister chromatids. These data suggest
the presence of an active mechanism for chromatin segregation in yeast.


  Advanced mathematics, such as multiscale weighted colored graph and element
specific persistent homology, and machine learning including deep neural
networks were integrated to construct mathematical deep learning models for
pose and binding affinity prediction and ranking in the last two D3R grand
challenges in computer-aided drug design and discovery. D3R Grand Challenge 2
(GC2) focused on the pose prediction and binding affinity ranking and free
energy prediction for Farnesoid X receptor ligands. Our models obtained the top
place in absolute free energy prediction for free energy Set 1 in Stage 2. The
latest competition, D3R Grand Challenge 3 (GC3), is considered as the most
difficult challenge so far. It has 5 subchallenges involving Cathepsin S and
five other kinase targets, namely VEGFR2, JAK2, p38-$\alpha$, TIE2, and ABL1.
There is a total of 26 official competitive tasks for GC3. Our predictions were
ranked 1st in 10 out of 26 official competitive tasks.


  The rapid advancement in high-throughput techniques has fueled the generation
of large volume of biological data rapidly with low cost. Some of these
techniques are microarray and next generation sequencing which provides genome
level insight of living cells. As a result, the size of most of the biological
databases, such as NCBI-GEO, NCBI-SRA, is exponentially growing. These
biological data are analyzed using computational techniques for knowledge
discovery - which is one of the objectives of bioinformatics research. Gene
regulatory network (GRN) is a gene-gene interaction network which plays pivotal
role in understanding gene regulation process and disease studies. From the
last couple of decades, the researchers are interested in developing
computational algorithms for GRN inference (GRNI) using high-throughput
experimental data. Several computational approaches have been applied for
inferring GRN from gene expression data including statistical techniques
(correlation coefficient), information theory (mutual information), regression
based approaches, probabilistic approaches (Bayesian networks, naive byes),
artificial neural networks, and fuzzy logic. The fuzzy logic, along with its
hybridization with other intelligent approach, is well studied in GRNI due to
its several advantages. In this paper, we present a consolidated review on
fuzzy logic and its hybrid approaches for GRNI developed during last two
decades.


  Virus capsids in interchromatin corrals of a cell nucleus are experimentally
known to exhibit anomalous diffusion as well as normal diffusion, leading to
the Gaussian distribution of the diffusion-exponent fluctuations over the
corrals. Here, the sojourn-time distribution of the virus capsid in local areas
of the corral, i.e., probability distribution of the sojourn time
characterizing diffusion in the local areas, is examined. Such an area is
regarded as a virtual cubic block, the diffusion property in which is normal or
anomalous. The distribution, in which the Gaussian fluctuation is incorporated,
is shown to tend to slowly decay. Then, the block-size dependence of average
sojourn time is discussed. A comment is also made on (non-)Markovianity of the
process of moving through the blocks.


  The genetic material of a eukaryotic cell comprises both nuclear DNA (ncDNA)
and mitochondrial DNA (mtDNA). These differ markedly in several aspects but
nevertheless must encode proteins that are compatible with one another. Here we
introduce a network model of the hypothetical co-evolution of the two most
common modes of cellular division for reproduction: by mitosis (supporting
asexual reproduction) and by meiosis (supporting sexual reproduction). Our
model is based on a random hypergraph, with two nodes for each possible
genotype, each encompassing both ncDNA and mtDNA. One of the nodes is
necessarily generated by mitosis occurring at a parent genotype, the other by
meiosis occurring at two parent genotypes. A genotype's fitness depends on the
compatibility of its ncDNA and mtDNA. The model has two probability parameters,
$p$ and $r$, the former accounting for the diversification of ncDNA during
meiosis, the latter for the diversification of mtDNA accompanying both meiosis
and mitosis. Another parameter, $\lambda$, is used to regulate the relative
rate at which mitosis- and meiosis-generated genotypes are produced. We have
found that, even though $p$ and $r$ do affect the existence of evolutionary
pathways in the network, the crucial parameter regulating the coexistence of
the two modes of cellular division is $\lambda$. Depending on genotype size,
$\lambda$ can be valued so that either mode of cellular division prevails. Our
study is closely related to a recent hypothesis that views the appearance of
cellular division by meiosis, as opposed to division by mitosis, as an
evolutionary strategy for boosting ncDNA diversification to keep up with that
of mtDNA. Our results indicate that this may well have been the case, thus
lending support to the first hypothesis in the field to take into account the
role of such ubiquitous and essential organelles as mitochondria.


  Probabilistic modeling is fundamental to the statistical analysis of complex
data. In addition to forming a coherent description of the data-generating
process, probabilistic models enable parameter inference about given data sets.
This procedure is well-developed in the Bayesian perspective, in which one
infers probability distributions describing to what extent various possible
parameters agree with the data. In this paper we motivate and review
probabilistic modeling for adaptive immune receptor repertoire data then
describe progress and prospects for future work, from germline haplotyping to
adaptive immune system deployment across tissues. The relevant quantities in
immune sequence analysis include not only continuous parameters such as gene
use frequency, but also discrete objects such as B cell clusters and lineages.
Throughout this review, we unravel the many opportunities for probabilistic
modeling in adaptive immune receptor analysis, including settings for which the
Bayesian approach holds substantial promise (especially if one is optimistic
about new computational methods). From our perspective the greatest prospects
for progress in probabilistic modeling for repertoires concern ancestral
sequence estimation for B cell receptor lineages, including uncertainty from
germline genotype, rearrangement, and lineage development.


  An important task for many if not all the scientific domains is efficient
knowledge integration, testing and codification. It is often solved with model
construction in a controllable computational environment. In spite of that, the
throughput of in-silico simulation-based observations become similarly
intractable for thorough analysis. This is especially the case in molecular
biology, which served as a subject for this study. In this project, we aimed to
test some approaches developed to deal with the curse of dimensionality. Among
these we found dimension reduction techniques especially appealing. They can be
used to identify irrelevant variability and help to understand critical
processes underlying high-dimensional datasets. Additionally, we subjected our
data sets to nonlinear time series analysis, as those are well established
methods for results comparison. To investigate the usefulness of dimension
reduction methods, we decided to base our study on a concrete sample set. The
example was taken from the domain of systems biology concerning dynamic
evolution of sub-cellular signaling. Particularly, the dataset relates to the
yeast pheromone pathway and is studied in-silico with a stochastic model. The
model reconstructs signal propagation stimulated by a mating pheromone. In the
paper, we elaborate on the reason of multidimensional analysis problem in the
context of molecular signaling, and next, we introduce the model of choice,
simulation details and obtained time series dynamics. A description of used
methods followed by a discussion of results and their biological interpretation
finalize the paper.


  A large number of recent genome-wide association studies (GWASs) for complex
phenotypes confirm the early conjecture for polygenicity, suggesting the
presence of large number of variants with only tiny or moderate effects.
However, due to the limited sample size of a single GWAS, many associated
genetic variants are too weak to achieve the genome-wide significance. These
undiscovered variants further limit the prediction capability of GWAS.
Restricted access to the individual-level data and the increasing availability
of the published GWAS results motivate the development of methods integrating
both the individual-level and summary-level data. How to build the connection
between the individual-level and summary-level data determines the efficiency
of using the existing abundant summary-level resources with limited
individual-level data, and this issue inspires more efforts in the existing
area.
  In this study, we propose a novel statistical approach, LEP, which provides a
novel way of modeling the connection between the individual-level data and
summary-level data. LEP integrates both types of data by \underline{LE}veraing
\underline{P}leiotropy to increase the statistical power of risk variants
identification and the accuracy of risk prediction. The algorithm for parameter
estimation is developed to handle genome-wide-scale data. Through comprehensive
simulation studies, we demonstrated the advantages of LEP over the existing
methods. We further applied LEP to perform integrative analysis of Crohn's
disease from WTCCC and summary statistics from GWAS of some other diseases,
such as Type 1 diabetes, Ulcerative colitis and Primary biliary cirrhosis. LEP
was able to significantly increase the statistical power of identifying risk
variants and improve the risk prediction accuracy from 63.39\% ($\pm$ 0.58\%)
to 68.33\% ($\pm$ 0.32\%) using about 195,000 variants.


  Targeting the mitochondrial enzyme FoF1-ATP synthase and modulating its
catalytic activities with small molecules is a promising new approach for
treatment of autoimmune diseases. The immuno-modulatory compound Bz-423 is such
a drug that binds to subunit OSCP of the mitochondrial FoF1-ATP synthase and
induces apoptosis via increased reactive oxygen production in coupled, actively
respiring mitochondria. Here we review the experimental progress to reveal the
binding of Bz-423 to the mitochondrial target and discuss how subunit rotation
of FoF1-ATP synthase is affected by Bz-423. Briefly, we report how F\"orster
resonance energy transfer (FRET) can be employed to colocalize the enzyme and
the fluorescently tagged Bz-423 within the mitochondria of living cells with
nanometer resolution.


  A matrix Lie algebra is a linear space of matrices closed under the operation
$ [A, B] = AB-BA $. The "Lie closure" of a set of matrices is the smallest
matrix Lie algebra which contains the set. In the context of Markov chain
theory, if a set of rate matrices form a Lie algebra, their corresponding
Markov matrices are closed under matrix multiplication; this has been found to
be a useful property in phylogenetics. Inspired by previous research involving
Lie closures of DNA models, it was hypothesised that finding the Lie closure of
a codon model could help to solve the problem of mis-estimation of the
non-synonymous/synonymous rate ratio, $ \omega $. We propose two different
methods of finding a linear space from a model: the first is the \emph{linear
closure} which is the smallest linear space which contains the model, and the
second is the \emph{linear version} which changes multiplicative constraints in
the model to additive ones. For each of these linear spaces we then find the
Lie closures of them. Under both methods, it was found that closed codon models
would require thousands of parameters, and that any partial solution to this
problem that was of a reasonable size violated stochasticity. Investigation of
toy models indicated that finding the Lie closure of matrix linear spaces which
deviated only slightly from a simple model resulted in a Lie closure that was
close to having the maximum number of parameters possible. Given that Lie
closures are not practical, we propose further consideration of the two
variants of linearly closed models.


  Aligning millions of short DNA or RNA reads, of 75 to 250 base pairs each, to
a reference genome is a significant computation problem in bioinformatics. We
present a flexible and fast FPGA-based short read alignment tool. Our aligner
makes use of the processing power of FPGAs in conjunction with the greater host
memory bandwidth and flexibility of software to improve performance and achieve
a high level of configurability. This flexible design supports a variety of
reference genome sizes without the performance degradation suffered by other
software and FPGA-based aligners. It is also better able to support the
features of new alignment algorithms, which frequently crop up in the rapidly
evolving field of bioinformatics. We demonstrate these advantages in a case
study where we align RNA-Seq data from a hypothesized mouse / human xenograft.
In this case study, our aligner provides a speedup of 5.6x over BWA-SW with
energy savings of 21%, while also reducing incorrect short read classification
by 29%. To demonstrate the flexibility of our system we show that the speedup
can be substantially improved while retaining most of the accuracy gains over
BWA-SW. The speedup can be increased to 71.3x, while still enjoying a 28%
incorrect classification improvement and 52% improvement in unaligned reads.


  In this study, we investigate a porous medium-type flux limited
reaction--diffusion equation that arises in morphogenesis modeling. This
nonlinear partial differential equation is an extension of the generalized
Fisher--Kolmogorov--Petrovsky--Piskunov (Fisher-KPP) equation in
one-dimensional space. The approximate analytical traveling wave solution is
found by using a perturbation method. We show that the morphogen concentration
propagates as a sharp wave front where the wave speed has a saturated value.
The numerical solutions of this equation are also provided to compare them with
the analytical predictions. Finally, we qualitatively compare our theoretical
results with those obtained in experimental studies.


  This paper is based on the complete classification of evolutionary scenarios
for the Moran process with two strategies given by Taylor et al. (B. Math.
Biol. 66(6): 1621--1644, 2004). Their classification is based on whether each
strategy is a Nash equilibrium and whether the fixation probability for a
single individual of each strategy is larger or smaller than its value for
neutral evolution. We improve on this analysis by showing that each
evolutionary scenario is characterized by a definite graph shape for the
fixation probability function. A second class of results deals with the
behavior of the fixation probability when the population size tends to
infinity. We develop asymptotic formulae that approximate the fixation
probability in this limit and conclude that some of the evolutionary scenarios
cannot exist when the population size is large.


  In phylogenetic studies, biologists often wish to estimate the ancestral
discrete character state at an interior vertex $v$ of an evolutionary tree $T$
from the states that are observed at the leaves of the tree. A simple and fast
estimation method --- maximum parsimony --- takes the ancestral state at $v$ to
be any state that minimises the number of state changes in $T$ required to
explain its evolution on $T$. In this paper, we investigate the reconstruction
accuracy of this estimation method further, under a simple symmetric model of
state change, and obtain a number of new results, both for 2-state characters,
and $r$--state characters ($r>2$). Our results rely on establishing new
identities and inequalities, based on a coupling argument that involves a
simpler `coin toss' approach to ancestral state reconstruction.


  Colletotrichum represent a genus of fungal species primarily known as plant
pathogens with severe economic impacts in temperate, subtropical and tropical
climates Consensus taxonomy and classification systems for Colletotrichum
species have been undergoing revision as high resolution genomic data becomes
available. Here we propose an alternative annotation that provides a complete
sequence for a Colletotrichum YPT1 gene homolog using the whole genome shotgun
sequence of Colletotrichum incanum isolated from soybean crops in Illinois,
USA.


  Some population is made of n individuals that can be of p possible species
(or types). The update of the species abundance occupancies is from a Moran
mutational model designed by Karlin and McGregor in 1967. We first study the
equilibrium species counts as a function of n, p and the total mutation
probability $\nu$ before considering various asymptotic regimes on n, p and
$\nu$. Running title: KMG Model with Mutations.


  This study explores the validity of chain effects of clean water, which are
known as the "Mills-Reincke phenomenon," in early twentieth-century Japan.
Recent studies have reported that water purifications systems are responsible
for huge contributions to human capital. Although some studies have
investigated the instantaneous effects of water-supply systems in pre-war
Japan, little is known about the chain effects of these systems. By analyzing
city-level cause-specific mortality data from 1922-1940, we find that a decline
in typhoid deaths by one per 1,000 people decreased the risk of death due to
non-waterborne diseases such as tuberculosis and pneumonia by 0.742-2.942 per
1,000 people. Our finding suggests that the observed Mills-Reincke phenomenon
could have resulted in the relatively rapid decline in the mortality rate in
early twentieth-century Japan.


  Modelling gene regulatory networks not only requires a thorough understanding
of the biological system depicted but also the ability to accurately represent
this system from a mathematical perspective. Throughout this chapter, we aim to
familiarise the reader with the biological processes and molecular factors at
play in the process of gene expression regulation.We first describe the
different interactions controlling each step of the expression process, from
transcription to mRNA and protein decay. In the second section, we provide
statistical tools to accurately represent this biological complexity in the
form of mathematical models. Amongst other considerations, we discuss the
topological properties of biological networks, the application of deterministic
and stochastic frameworks and the quantitative modelling of regulation. We
particularly focus on the use of such models for the simulation of expression
data that can serve as a benchmark for the testing of network inference
algorithms.


  The Colless index is one of the most popular and natural balance indices for
bifurcating phylogenetic trees, but it makes no sense for multifurcating trees.
In this paper we propose a family of Colless-like balance indices
$\mathfrak{C}_{D,f}$, which depend on a dissimilarity $D$ and a function
$f:\mathbb{N}\to \mathbb{R}_{\geq 0}$, that generalize the Colless index to
multifurcating phylogenetic trees. We provide two functions $f$ such that the
most balanced phylogenetic trees according to the corresponding indices
$\mathfrak{C}_{D,f}$ are exactly the fully symmetric ones. Next, for each one
of these two functions $f$ and for three popular dissimilarities $D$ (the
variance, the standard deviation, and the mean deviation from the median), we
determine the range of values of $\mathfrak{C}_{D,f}$ on the sets of
phylogenetic trees with a given number $n$ of leaves. We end the paper by
assessing the performance of one of these indices on TreeBASE and using it to
show that the trees in this database do not seem to follow either the uniform
model for multifurcating trees or the $\alpha$-$\gamma$-model, for any values
of $\alpha$ and $\gamma$.


  The objective of this article is to create a framework to study asymptotic
equilibria in human populations with a special focus on immigration. We present
a new model, based on Resource Dependent Branching Processes, which is now
broad enough to cope with the goal of finding equilibrium criteria under
reasonable hypotheses. Our equations are expressed in terms of natality rates,
mean productivity and mean consumption of the home-population and the immigrant
population as well as policies of the Society to distribute resources among
individuals. We also study the impact of integration of one sub-population into
the other one, and in a third model, the additional influence of a continuous
stream of new immigrants. Proofs of the results are based on classical limit
theorems, on Borel-Cantelli type arguments, on the Theorem of envelopment of
Bruss and Duerinckx (2015), on a maximum inequality of Bruss and Robertson
(1991) and, in particular, on an extension of J.M. Steele (2016) of the latter.
Conditions for the existence of an equilibrium often prove to be severe, and
sometimes surprisingly sensitive. This underlines how demanding the real world
of immigration can be for politicians trying to make sound decisions. Our main
objective is to provide help through insights from an adequate theory. Another
objective of the present study is to learn which of the possible control
measures are best for combining feasibility and efficiency to reach an
equilibrium, and to recognise the corresponding steps one has to take towards
controls. We also make preliminary suggestions to envisage ways to optimal
control.
  As far as the author is aware, all results are new.


  We propose a stochastic SICA epidemic model for HIV transmission, described
by stochastic ordinary differential equations, and discuss its perturbation by
environmental white noise. Existence and uniqueness of the global positive
solution to the stochastic HIV system is proven, and conditions under which
extinction and persistence in mean hold, are given. The theoretical results are
illustrated via numerical simulations.


  The complex dynamics of gene expression in living cells can be
well-approximated using Boolean networks. The average sensitivity is a natural
measure of stability in these systems: values below one indicate typically
stable dynamics associated with an ordered phase, whereas values above one
indicate chaotic dynamics. This yields a theoretically motivated adaptive
advantage to being near the critical value of one, at the boundary between
order and chaos. Here, we measure average sensitivity for 66 publicly available
Boolean network models describing the function of gene regulatory circuits
across diverse living processes. We find the average sensitivity values for
these networks are clustered around unity, indicating they are near critical.
In many types of random networks, mean connectivity <K> and the average
activity bias of the logic functions <p> have been found to be the most
important network properties in determining average sensitivity, and by
extension a network's criticality. Surprisingly, many of these gene regulatory
networks achieve the near-critical state with <K> and <p> far from that
predicted for critical systems: randomized networks sharing the local causal
structure and local logic of biological networks better reproduce their
critical behavior than controlling for macroscale properties such as <K> and
<p> alone. This suggests the local properties of genes interacting within
regulatory networks are selected to collectively be near-critical, and this
non-local property of gene regulatory network dynamics cannot be predicted
using the density of interactions alone.


  Reconstruction of gene regulatory networks is the process of identifying gene
dependency from gene expression profile through some computation techniques. In
our human body, though all cells pose similar genetic material but the
activation state may vary. This variation in the activation of genes helps
researchers to understand more about the function of the cells. Researchers get
insight about diseases like mental illness, infectious disease, cancer disease
and heart disease from microarray technology, etc. In this study, a
cancer-specific gene regulatory network has been constructed using a simple and
novel machine learning approach. In First Step, linear regression algorithm
provided us the significant genes those expressed themselves differently. Next,
regulatory relationships between the identified genes has been computed using
Pearson correlation coefficient. Finally, the obtained results have been
validated with the available databases and literatures. We can identify the hub
genes and can be targeted for the cancer diagnosis.


  How to prevent the spread of human diseases is a great challenge for the
scientific community and so far there are many studies in which immunization
strategies have been developed. However, these kind of strategies usually do
not consider that medical institutes may have limited vaccine resources
available. In this manuscript, we explore the Susceptible-Infected-Recovered
(SIR) model with local dynamic vaccination, and considering limited vaccines.
In this model, susceptibles in contact with an infected individual, are
vaccinated -with probability $\omega$- and then get infected -with probability
$\beta$. However, when the fraction of immunized individuals reaches a
threshold $V_L$, the vaccination stops, after which only the infection is
possible. In the steady state, besides the critical points $\beta_c$ and
$\omega_c$ that separate a non-epidemic from an epidemic phase, we find for a
range of $V_L$ another transition points, $\beta^*>\beta_c$ and
$\omega^*<\omega_c$, which correspond to a novel discontinuous phase
transition. This critical value separates a phase where the amount of vaccines
is sufficient, from a phase where the disease is strong enough to exhaust all
the vaccination units. For a disease with fixed $\beta$, the vaccination
probability $\omega$ can be controlled in order to drastically reduce the
number of infected individuals, using efficiently the available vaccines.
  Furthermore, the temporal evolution of the system close to $\beta^*$ or
$\omega^*$, shows that after a peak of infection the system enters into a
quasi-stationary state, with only a few infected cases. But if there are no
more vaccines, these few infected individuals could originate a second
outbreak, represented by a second peak of infection. This state of apparent
calm, could be dangerous since it may lead to misleading conclusions and to an
abandon of the strategies to control the disease.


  I review the evolutionary history of human populations in Europe with an
emphasis on what has been learned in recent years through the study of ancient
DNA. Human populations in Europe ~430-39kya (archaic Europeans) included
Neandertals and their ancestors, who were genetically differentiated from other
archaic Eurasians (such as the Denisovans of Siberia), as well as modern
humans. Modern humans arrived to Europe by ~45kya, and are first genetically
attested by ~39kya when they were still mixing with Neandertals. The first
Europeans who were recognizably genetically related to modern ones appeared in
the genetic record shortly thereafter at ~37kya. At ~15kya a largely
homogeneous set of hunter-gatherers became dominant in most of Europe, but with
some admixture from Siberian hunter-gatherers in the eastern part of the
continent. These hunter-gatherers were joined by migrants from the Near East
beginning at ~8kya: Anatolian farmers settled most of mainland Europe, and
migrants from the Caucasus reached eastern Europe, forming steppe populations.
After ~5kya there was migration from the steppe into mainland Europe and vice
versa. Present-day Europeans (ignoring the long-distance migrations of the
modern era) are largely the product of this Bronze Age collision of steppe
pastoralists with Neolithic farmers.


  Biological networks are a very convenient modelling and visualisation tool to
discover knowledge from modern high-throughput genomics and postgenomics data
sets. Indeed, biological entities are not isolated, but are components of
complex multi-level systems. We go one step further and advocate for the
consideration of causal representations of the interactions in living
systems.We present the causal formalism and bring it out in the context of
biological networks, when the data is observational. We also discuss its
ability to decipher the causal information flow as observed in gene expression.
We also illustrate our exploration by experiments on small simulated networks
as well as on a real biological data set.


  The meaningful comparison of ion mobility (IM) results and of collision cross
section (CCS) values on different platforms is a prerequisite for using CCS for
identification or structural assignment. The amount of internal energy imparted
to the ions prior to the ion mobility cell is a source of experimental
variation. Here we investigated the effects of virtually all tuning parameters
of the Agilent 6560 IM-Q-TOF on the arrival time distributions of Ubiquitin7+,
and found conditions in which the native state prevails. We will discuss the
effects of solvent evaporation conditions in the source, in the entire pre-IM
DC voltage gradient, and with the funnel RF amplitudes, and will also report on
ubiquitin7+ conformations in different solvents, including native supercharging
conditions. Collision-induced unfolding (CIU) can be conveniently provoked in
two distinct regions: behind the source capillary (by changing the fragmentor
voltage) and in the trapping funnel (by changing the trap entrance grid delta
voltage). The softness of the instrumental conditions were then optimized with
the benchmark DNA G-quadruplex [(dG4T4G4)2.(NH4+)3-8H]5-, for which ion
activation results in ammonia loss. To reduce the ion internal energy and
obtain the intact 3-NH4+ complex, we reduced the post-IM voltage gradient, but
this resulted in a lower IM resolving power due to increased diffusion behind
the drift tube. The article thus describes the various trade-offs between ion
activation, ion transmission, and ion mobility performance for native MS of
very fragile structures.


  We explore the emergence of persistent infection in a closed region where the
disease progression of the individuals is given by the SIRS model, with an
individual becoming infected on contact with another infected individual within
a given range. We focus on the role of synchronization in the persistence of
contagion. Our key result is that higher degree of synchronization, both
globally in the population and locally in the neighborhoods, hinders
persistence of infection. Importantly, we find that early short-time asynchrony
appears to be a consistent precursor to future persistence of infection, and
can potentially provide valuable early warnings for sustained contagion in a
population patch. Thus transient synchronization can help anticipate the
long-term persistence of infection. Further we demonstrate that when the range
of influence of an infected individual is wider, one obtains lower persistent
infection. This counter-intuitive observation can also be understood through
the relation of synchronization to infection burn-out.


  The bacterial SOS response is a cellular reaction to DNA damage, that, among
other actions, triggers the expression of colicin - toxic bacteriocins in
Escherichia coli that are released to kill close relatives competing for
resources. However, it is largely unknown, how the complex network regulating
toxin expression controls the time-point of toxin release to prevent premature
release of inefficient protein concentrations. Here, we study how different
regulatory mechanisms affect production and release of the bacteriocin
ColicinE2 in Escherichia coli. Combining experimental and theoretical
approaches, we demonstrate that the global carbon storage regulator CsrA
controls the duration of the delay between toxin production and release and
emphasize the importance of CsrA sequestering elements for the timing of
ColicinE2 release. In particular, we show that ssDNA originating from
rolling-circle replication of the toxin-producing plasmid represents a yet
unknown additional CsrA sequestering element, which is essential in the
ColicinE2-producing strain to enable toxin release by reducing the amount of
free CsrA molecules in the bacterial cell. Taken together, our findings show
that CsrA times ColicinE2 release and reveal a dual function for CsrA as an
ssDNA and mRNA-binding protein, introducing ssDNA as an important
post-transcriptional gene regulatory element.


  In genetic drift of small population, it is well known that even when the
ratio of alleles is 0.5, specific genes are fixed in or disappear from the
population. It seems the reason why inbreeding is avoided. On the other hand,
this phenomenon suggests an interesting possibility. The mutant gene does not
increase the number of genes at once in a large population. A gene is partially
fixed by increasing the number within a small population because of inbreeding,
and the gene increases in a large group by Darwin's natural selection. It would
be more reasonable to think in this way. We studied this mathematically based
on the concept of genetic drift. This suggested that inbreeding could be useful
as a trigger for fixation of mutation.


  The aim of this work is to provide a rigorous mathematical analysis of a
stochastic concatenation model presented by Sobottka and Hart (2011) which
allows approximation of the first-order stochastic structure in bacterial DNA
by means of a stationary Markov chain. Two probabilistic constructions that
rigorously formalize the model are presented. Necessary and sufficient
conditions for a Markov chain to be generated by the model are given, as well
as the theoretical background needed for designing new algorithms for
statistical analyses of real bacterial genomes. It is shown that the model
encompasses the Markov chains satisfying intra-strand parity, a property
observed in most DNA sequences.


  Uncovering modular structure in networks is fundamental for systems in
biology, physics, and engineering. Community detection identifies candidate
modules as hypotheses, which then need to be validated through experiments,
such as mutagenesis in a biological laboratory. Only a few communities can
typically be validated, and it is thus important to prioritize which
communities to select for downstream experimentation. Here we develop CRank, a
mathematically principled approach for prioritizing network communities. CRank
efficiently evaluates robustness and magnitude of structural features of each
community and then combines these features into the community prioritization.
CRank can be used with any community detection method. It needs only
information provided by the network structure and does not require any
additional metadata or labels. However, when available, CRank can incorporate
domain-specific information to further boost performance. Experiments on many
large networks show that CRank effectively prioritizes communities, yielding a
nearly 50-fold improvement in community prioritization.


  GeneVis is a web-based tool to visualize complementary data sets of different
disciplines within the field of genetics. It overlays gene-cluster information,
gene-interaction data and gene-disease association data by means of web-based
interactive graph visualizations. This allows an intuitive and quick assessment
of possible relations between the different datasets. By starting from a
high-level graph abstraction based on gene clusters, which can be selected for
detailed inspection at the gene-interaction level in a separate window, GeneVis
circumvents the common visual clutter problem when using gene datasets with a
high number of gene entries.


  The production of mechanical stresses in living organisms largely relies on
localized, force-generating active units embedded in filamentous matrices.
Numerical simulations of discrete fiber networks with fixed boundaries have
shown that buckling in the matrix dramatically amplifies the resulting active
stresses. Here we extend this result to a bucklable continuum elastic medium
subjected to an arbitrary external stress, and derive analytical expressions
for the active, nonlinear constitutive relations characterizing the full active
medium. Inserting these relations into popular "active gel" descriptions of
living tissues and the cytoskeleton will enable investigations into nonlinear
regimes previously inaccessible due to the phenomenological nature of these
theories.


  Prosociality is fundamental to human social life, and, accordingly, much
research has attempted to explain human prosocial behavior. Capraro and Rand
(Judgment and Decision Making, 13, 99-111, 2018) recently provided experimental
evidence that prosociality in anonymous, one-shot interactions (such as
Prisoner's Dilemma and Dictator Game experiments) is not driven by
outcome-based social preferences - as classically assumed - but by a
generalized morality preference for "doing the right thing". Here we argue that
the key experiments reported in Capraro and Rand (2018) comprise prominent
methodological confounds and open questions that bear on influential
psychological theory. Specifically, their design confounds: (i) preferences for
efficiency with self-interest; and (ii) preferences for action with preferences
for morality. Furthermore, their design fails to dissociate the preference to
do "good" from the preference to avoid doing "bad". We thus designed and
conducted a preregistered, refined and extended test of the morality preference
hypothesis (N=801). Consistent with this hypothesis, our findings indicate that
prosociality in the anonymous, one-shot Dictator Game is driven by preferences
for doing the morally right thing. Inconsistent with influential psychological
theory, however, our results suggest the preference to do "good" was as potent
as the preference to avoid doing "bad" in this case.


  We report a procedure that, in one step from continuous data with minimal
preparation, recovers the graph found by Sachs et al. \cite{sachs2005causal},
with only a few edges different. The algorithm, Fast Adjacency Skewness (FASK),
relies on a mixture of linear reasoning and reasoning from the skewness of
variables; the Sachs data is a good candidate for this procedure since the
skewness of the variables is quite pronounced. We review the ground truth model
from Sachs et al. as well as some of the fluctuations seen in the protein
abundances in the system, give the Sachs model and the FASK model, and perform
a detailed comparison. Some variation in hyper-parameters is explored, though
the main result uses values at or near the defaults learned from work modeling
fMRI data.


  We study the problem of identifying the source of emerging large-scale
outbreaks of foodborne disease. To solve the source identification problem we
formulate a probabilistic model of the contamination diffusion process as a
random walk on a network and derive the maximum likelihood estimator for the
source location. By modeling the transmission process as a random walk, we are
able to develop a novel, computationally tractable solution that accounts for
all possible paths of travel through the network. This is in contrast to
existing approaches to network source identification, which assume that the
contamination travels along either the shortest or highest probability paths.
We demonstrate the benefits of the multiple-paths approach through application
to different network topologies, including stylized models of food supply
network structure and real data from the 2011 EHEC outbreak in Germany. We show
significant improvements in accuracy and reliability compared with the relevant
state-of-the-art approach to source identification. Beyond foodborne disease,
these methods should find application in identifying the source of spread in
network-based diffusion processes more generally, including in networks not
well approximated by tree-like structure.


  Reaching the full potential of precision medicine depends on the quality of
personalized genome interpretation. In order to facilitate precision medicine
in regions of the Middle East and North Africa (MENA), a population-specific
reference genome for the indigenous Arab popula-tion of Qatar (QTRG) was
constructed by incorporating allele frequency data from sequencing of 1,161
Qataris, representing 0.4% of the population. A total of 20.9 million SNP and
3.1 million indels were observed in Qatar, including an average of 1.79% novel
variants per individual ge-nome. Replacement of the GRCh37 standard reference
with QTRG in a best practices genome analysis workflow resulted in an average
of 7* deeper coverage depth (an improvement of 23%), and 756,671 fewer variants
on average, a reduction of 16% that is attributed to common Qatari alleles
being present in the QTRG reference. The benefit for using QTRG varies across
ances-tries, a factor that should be taken into consideration when selecting an
appropriate reference for analysis.


  Infection by many viruses begins with fusion of viral and cellular lipid
membranes, followed by entry of viral contents into the target cell and
ultimately, after many biochemical steps, integration of viral DNA into that of
the host cell. The early steps of membrane fusion and viral capsid entry are
mediated by adsorption to the cell surface, and receptor and coreceptor
binding. HIV-1 specifically targets CD4+ helper T-cells of the human immune
system and binds to the receptor CD4 and coreceptor CCR5 before fusion is
initiated. Previous experiments have been performed using a cell line
(293-Affinofile) in which the expression of CD4 and CCR5 concentration were
independently controlled. After exposure to HIV-1 of various strains, the
resulting infectivity was measured through the fraction of infected cells. To
design and evaluate the effectiveness of drug therapies that target the
inhibition of the entry processes, an accurate functional relationship between
the CD4/CCR5 concentrations and infectivity is desired in order to more
quantitatively analyze experimental data. We propose three kinetic models
describing the possible mechanistic processes involved in HIV entry and fit
their predictions to infectivity measurements, contrasting and comparing
different outcomes. Our approach allows interpretation of the clustering of
infectivity of different strains of HIV-1 in the space of mechanistic kinetic
parameters. Our model fitting also allows inference of nontrivial
stoichiometries of receptor and coreceptor binding and provides a framework
through which to quantitatively investigate the effectiveness of fusion
inhibitors and neutralizing antibodies.


  Networks are ubiquitous in biology where they encode connectivity patterns at
all scales of organization, from molecular to the biome. However, biological
networks are noisy due to the limitations of measurement technology and
inherent natural variation, which can hamper discovery of network patterns and
dynamics. We propose Network Enhancement (NE), a method for improving the
signal-to-noise ratio of undirected, weighted networks. NE uses a doubly
stochastic matrix operator that induces sparsity and provides a closed-form
solution that increases spectral eigengap of the input network. As a result, NE
removes weak edges, enhances real connections, and leads to better downstream
performance. Experiments show that NE improves gene function prediction by
denoising tissue-specific interaction networks, alleviates interpretation of
noisy Hi-C contact maps from the human genome, and boosts fine-grained
identification accuracy of species. Our results indicate that NE is widely
applicable for denoising biological networks.


  This article is aimed at proposing a generalization of the Ross-Macdonald
model for the transmission of Vector-borne diseases in which human-to-human
contagions are also considered. We first present this generalized model by
formulating a mean field theory, checking its validity by comparing to
numerical simulations. To make the premises of our model more realistic, we
adapt the mean field equations to the case in which human contacts are
described by a complex network. In this case we are also able to derive an
analytical expression for the epidemic threshold. In both the mean-field and
network-based models, we estimate the value of the epidemic threshold which
corresponds to the boundary between the disease-free and epidemic regimes. The
expression of this threshold allows us to discuss the impact that
human-to-human contagions have on the spread of vector-borne diseases.


  Empirical observations show that ecological communities can have a huge
number of coexisting species, also with few or limited number of resources.
These ecosystems are characterized by multiple type of interactions, in
particular displaying cooperative behaviors. However, standard modeling of
population dynamics based on Lotka-Volterra type of equations predicts that
ecosystem stability should decrease as the number of species in the community
increases and that cooperative systems are less stable than communities with
only competitive and/or exploitative interactions. Here we propose a stochastic
model of population dynamics, which includes exploitative interactions as well
as cooperative interactions induced by cross-feeding. The model is exactly
solved and we obtain results for relevant macro-ecological patterns, such as
species abundance distributions and correlation functions. In the large system
size limit, any number of species can coexist for a very general class of
interaction networks and stability increases as the number of species grows.
For pure mutualistic/commensalistic interactions we determine the topological
properties of the network that guarantee species coexistence. We also show that
the stationary state is globally stable and that inferring species interactions
through species abundance correlation analysis may be misleading. Our
theoretical approach thus show that appropriate models of cooperation naturally
leads to a solution of the long-standing question about complexity-stability
paradox and on how highly biodiverse communities can coexist.


  Supertree construction is the process by which a set of phylogenetic trees,
each on a subset of the overall set X of species, is combined into a tree on
the full set S. The traditional use of supertree methods is the assembly of a
large species tree from previously computed smaller species trees; however,
supertree methods are also used to address large-scale tree estimation using
divide-and-conquer (i.e., a dataset is divided into overlapping subsets, trees
are constructed on the subsets, and then combined using the supertree method).
Because most supertree methods are heuristics for NP-hard optimization
problems, the use of supertree estimation on large datasets is challenging,
both in terms of scalability and accuracy. In this paper, we describe the
current state of the art in supertree construction and the use of supertree
methods in divide-and-conquer strategies. Finally, we identify directions where
future research could lead to improved supertree methods.


  High quality gene models are necessary to expand the molecular and genetic
tools available for a target organism, but these are available for only a
handful of model organisms that have undergone extensive curation and
experimental validation over the course of many years. The majority of gene
models present in biological databases today have been identified in draft
genome assemblies using automated annotation pipelines that are frequently
based on orthologs from distantly related model organisms. Manual curation is
time consuming and often requires substantial expertise, but is instrumental in
improving gene model structure and identification. Manual annotation may seem
to be a daunting and cost-prohibitive task for small research communities but
involving undergraduates in community genome annotation consortiums can be
mutually beneficial for both education and improved genomic resources. We
outline a workflow for efficient manual annotation driven by a team of
primarily undergraduate annotators. This model can be scaled to large teams and
includes quality control processes through incremental evaluation. Moreover, it
gives students an opportunity to increase their understanding of genome biology
and to participate in scientific research in collaboration with peers and
senior researchers at multiple institutions.


  The response of proteins to chemical reactions or impulsive excitation that
occurs within the molecule has fascinated chemists for decades. In recent years
ultrafast X-ray studies have provided ever more detailed information about the
evolution of protein structural change following ligand photolysis, and
time-resolved IR and Raman techniques, e.g., have provided detailed pictures of
the nature and rate of energy transport in peptides and proteins, including
recent advances in identifying transport through individual amino acids of
several heme proteins. Computational tools to locate energy transport pathways
in proteins have also been advancing. Energy transport pathways in proteins
have since some time been identified by molecular dynamics (MD) simulations,
and more recent efforts have focused on the development of coarse graining
approaches, some of which have exploited analogies to thermal transport in
other molecular materials. With the identification of pathways in proteins and
protein complexes, network analysis has been applied to locate residues that
control protein dynamics and possibly allostery, where chemical reactions at
one binding site mediate reactions at distance sites of the protein. In this
chapter we review approaches for locating computationally energy transport
networks in proteins. We present background into energy and thermal transport
in condensed phase and macromolecules that underlies the approaches we discuss
before turning to a description of the approaches themselves. We also
illustrate the application of the computational methods for locating energy
transport networks and simulating energy dynamics in proteins with several
examples.


  We present ensemble methods in a machine learning (ML) framework combining
predictions from five known motif/binding site exploration algorithms. For a
given TF the ensemble starts with position weight matrices (PWM's) for the
motif, collected from the component algorithms. Using dimension reduction, we
identify significant PWM-based subspaces for analysis. Within each subspace a
machine classifier is built for identifying the TF's gene (promoter) targets
(Problem 1). These PWM-based subspaces form an ML-based sequence analysis tool.
Problem 2 (finding binding motifs) is solved by agglomerating k-mer (string)
feature PWM-based subspaces that stand out in identifying gene targets. We
approach Problem 3 (binding sites) with a novel machine learning approach that
uses promoter string features and ML importance scores in a classification
algorithm locating binding sites across the genome. For target gene
identification this method improves performance (measured by the F1 score) by
about 10 percentage points over the (a) motif scanning method and (b) the
coexpression-based association method. Top motif outperformed 5 component
algorithms as well as two other common algorithms (BEST and DEME). For
identifying individual binding sites on a benchmark cross species database
(Tompa et al., 2005) we match the best performer without much human
intervention. It also improved the performance on mammalian TFs.
  The ensemble can integrate orthogonal information from different weak
learners (potentially using entirely different types of features) into a
machine learner that can perform consistently better for more TFs. The TF gene
target identification component (problem 1 above) is useful in constructing a
transcriptional regulatory network from known TF-target associations. The
ensemble is easily extendable to include more tools as well as future PWM-based
information.


  Many biological activities are induced by cellular chemical reactions of
diffusing reactants. The dynamics of such systems can be captured by stochastic
reaction networks. A recent numerical study has shown that diffusion can
significantly enhance the fluctuations in gene regulatory networks. However,
the universal relation between diffusion and stochastic system dynamics remains
veiled. Within the approximation of reaction-diffusion master equation (RDME),
we find general relation that the steady-state distribution in complex balanced
networks is diffusion-independent. Here, complex balance is the nonequilibrium
generalization of detailed balance. We also find that for a diffusion-included
network with a Poisson-like steady-state distribution, the diffusion can be
ignored at steady state. We then derive a necessary and sufficient condition
for networks holding such steady-state distributions. Moreover, we show that
for linear reaction networks the RDME reduces to the chemical master equation,
which implies that the stochastic dynamics of networks is unaffected by
diffusion at any arbitrary time. Our findings shed light on the fundamental
question of when diffusion can be neglected, or (if nonnegligible) its effects
on the stochastic dynamics of the reaction network.


  The Alzheimer's Disease Prediction Of Longitudinal Evolution (TADPOLE)
Challenge compares the performance of algorithms at predicting future evolution
of individuals at risk of Alzheimer's disease. TADPOLE Challenge participants
train their models and algorithms on historical data from the Alzheimer's
Disease Neuroimaging Initiative (ADNI) study or any other datasets to which
they have access. Participants are then required to make monthly forecasts over
a period of 5 years from January 2018, of three key outcomes for ADNI-3
rollover participants: clinical diagnosis, Alzheimer's Disease Assessment Scale
Cognitive Subdomain (ADAS-Cog13), and total volume of the ventricles. These
individual forecasts are later compared with the corresponding future
measurements in ADNI-3 (obtained after the TADPOLE submission deadline). The
first submission phase of TADPOLE was open for prize-eligible submissions
between 15 June and 15 November 2017. The submission system remains open via
the website: https://tadpole.grand-challenge.org, although since 15 November
2017 submissions are not eligible for the first round of prizes. This paper
describes the design of the TADPOLE Challenge.


  By providing a framework of accounting for the shared ancestry inherent to
all life, phylogenetics is becoming the statistical foundation of biology. The
importance of model choice continues to grow as phylogenetic models continue to
increase in complexity to better capture micro and macroevolutionary processes.
In a Bayesian framework, the marginal likelihood is how data update our prior
beliefs about models, which gives us an intuitive measure of comparing model
fit that is grounded in probability theory. Given the rapid increase in the
number and complexity of phylogenetic models, methods for approximating
marginal likelihoods are increasingly important. Here we try to provide an
intuitive description of marginal likelihoods and why they are important in
Bayesian model testing. We also categorize and review methods for estimating
marginal likelihoods of phylogenetic models, highlighting several recent
methods that provide well-behaved estimates. Furthermore, we review some
empirical studies that demonstrate how marginal likelihoods can be used to
learn about models of evolution from biological data. We discuss promising
alternatives that can complement marginal likelihoods for Bayesian model
choice, including posterior-predictive methods. Using simulations, we find one
alternative method based on approximate-Bayesian computation (ABC) to be
biased. We conclude by discussing the challenges of Bayesian model choice and
future directions that promise to improve the approximation of marginal
likelihoods and Bayesian phylogenetics as a whole.


  The mainstream of research in genetics, epigenetics and imaging data analysis
focuses on statistical association or exploring statistical dependence between
variables. Despite their significant progresses in genetic research,
understanding the etiology and mechanism of complex phenotypes remains elusive.
Using association analysis as a major analytical platform for the complex data
analysis is a key issue that hampers the theoretic development of genomic
science and its application in practice. Causal inference is an essential
component for the discovery of mechanical relationships among complex
phenotypes. Many researchers suggest making the transition from association to
causation. Despite its fundamental role in science, engineering and
biomedicine, the traditional methods for causal inference require at least
three variables. However, quantitative genetic analysis such as QTL, eQTL,
mQTL, and genomic-imaging data analysis requires exploring the causal
relationships between two variables. This paper will focus on bivariate causal
discovery. We will introduce independence of cause and mechanism (ICM) as a
basic principle for causal inference, algorithmic information theory and
additive noise model (ANM) as major tools for bivariate causal discovery.
Large-scale simulations will be performed to evaluate the feasibility of the
ANM for bivariate causal discovery. To further evaluate their performance for
causal inference, the ANM will be applied to the construction of gene
regulatory networks. Also, the ANM will be applied to trait-imaging data
analysis to illustrate three scenarios: presence of both causation and
association, presence of association while absence of causation, and presence
of causation, while lack of association between two variables.


  We give a combinatorial description of the toric ideal of invariants of the
Cavender-Farris-Neyman model with a molecular clock (CFN-MC) on a rooted binary
phylogenetic tree and prove results about the polytope associated to this toric
ideal. Key results about the polyhedral structure include that the number of
vertices of this polytope is a Fibonacci number, the facets of the polytope can
be described using the combinatorial "cluster" structure of the underlying
rooted tree, and the volume is equal to an Euler zig-zag number. The toric
ideal of invariants of the CFN-MC model has a quadratic Groebner basis with
squarefree initial terms. Finally, we show that the Ehrhart polynomial of these
polytopes, and therefore the Hilbert series of the ideals, depends only on the
number of leaves of the underlying binary tree, and not on the topology of the
tree itself. These results are analogous to classic results for the
Cavender-Farris-Neyman model without a molecular clock. However, new techniques
are required because the molecular clock assumption destroys the toric fiber
product structure that governs group-based models without the molecular clock.


  Multi-strain competition on networks is observed in many contexts, including
infectious disease ecology, information dissemination or behavioral adaptation
to epidemics. Despite a substantial body of research has been developed
considering static, time-aggregated networks, it remains a challenge to
understand the transmission of concurrent strains when links of the network are
created and destroyed over time. Here we analyze how network dynamics shapes
the outcome of the competition between an initially endemic strain and an
emerging one, when both strains follow a susceptible-infected-susceptible
dynamics, and spread at time scales comparable with the network evolution one.
Using time-resolved data of close-proximity interactions between patients
admitted to a hospital and medical health care workers, we analyze the impact
of temporal patterns and initial conditions on the dominance diagram and
coexistence time. We find that strong variations in activity volume cause the
probability that the emerging strain replaces the endemic one to be highly
sensitive to the time of emergence. The temporal structure of the network
shapes the dominance diagram, with significant variations in the replacement
probability (for a given set of epidemiological parameters) observed from the
empirical network and a randomized version of it. Our work contributes towards
the description of the complex interplay between competing pathogens on
temporal networks.


  Many organisms repartition their proteome in a circadian fashion in response
to the daily nutrient changes in their environment. A striking example is
provided by cyanobacteria, which perform photosynthesis during the day to fix
carbon. These organisms not only face the challenge of rewiring their proteome
every 12 hours, but also the necessity of storing the fixed carbon in the form
of glycogen to fuel processes during the night. In this manuscript, we extend
the framework developed by Hwa and coworkers (Scott et al., Science 330, 1099
(2010)) for quantifying the relatinship between growth and proteome composition
to circadian metabolism. We then apply this framework to investigate the
circadian metabolism of the cyanobacterium Cyanothece, which not only fixes
carbon during the day, but also nitrogen during the night, storing it in the
polymer cyanophycin. Our analysis reveals that the need to store carbon and
nitrogen tends to generate an extreme growth strategy, in which the cells
predominantly grow during the day, as observed experimentally. This strategy
maximizes the growth rate over 24 hours, and can be quantitatively understood
by the bacterial growth laws. Our analysis also shows that the slow relaxation
of the proteome, arising from the slow growth rate, puts a severe constraint on
implementing this optimal strategy. Yet, the capacity to estimate the time of
the day, enabled by the circadian clock, makes it possible to anticipate the
daily changes in the environment and mount a response ahead of time. This
significantly enhances the growth rate by counteracting the detrimental effects
of the slow proteome relaxation.


  Boolean tensor decomposition approximates data of multi-way binary
relationships as product of interpretable low-rank binary factors, following
the rules of Boolean algebra. Here, we present its first probabilistic
treatment. We facilitate scalable sampling-based posterior inference by
exploitation of the combinatorial structure of the factor conditionals. Maximum
a posteriori decompositions feature higher accuracies than existing techniques
throughout a wide range of simulated conditions. Moreover, the probabilistic
approach facilitates the treatment of missing data and enables model selection
with much greater accuracy. We investigate three real-world data-sets. First,
temporal interaction networks in a hospital ward and behavioural data of
university students demonstrate the inference of instructive latent patterns.
Next, we decompose a tensor with more than 10 billion data points, indicating
relations of gene expression in cancer patients. Not only does this demonstrate
scalability, it also provides an entirely novel perspective on relational
properties of continuous data and, in the present example, on the molecular
heterogeneity of cancer. Our implementation is available on GitHub:
https://github.com/TammoR/LogicalFactorisationMachines.