  To date, most genetic analyses of phenotypes have focused on analyzing single
traits or, analyzing each phenotype independently. However, joint epistasis
analysis of multiple complementary traits will increase statistical power, and
hold the key to understanding the complicated genetic structure of the complex
diseases. Despite their importance in uncovering the genetic structure of
complex traits, the statistical methods for identifying epistasis in multiple
phenotypes remains fundamentally unexplored. To fill this gap, we formulate a
test for interaction between two gens in multiple quantitative trait analysis
as a multiple functional regression (MFRG) in which the genotype functions
(genetic variant profiles) are defined as a function of the genomic position of
the genetic variants. We use large scale simulations to calculate its type I
error rates for testing interaction between two genes with multiple phenotypes
and to compare its power with multivariate pair-wise interaction analysis and
single trait interaction analysis by a single variate functional regression
model. To further evaluate its performance, the MFRG for epistasis analysis is
applied to five phenotypes and exome sequence data from the NHLBI Exome
Sequencing Project (ESP) to detect pleiotropic epistasis. A total of 136 pairs
of genes that formed a genetic interaction network showed significant evidence
of epistasis influencing five traits. The results demonstrate that the joint
interaction analysis of multiple phenotypes has much higher power to detect
interaction than the interaction analysis of single trait and may open a new
direction to fully uncovering the genetic structure of multiple phenotypes.


  We investigate Bayesian non-parametric inference of the $\Lambda$-measure of
$\Lambda$-coalescent processes with recurrent mutation, parametrised by
probability measures on the unit interval. We give verifiable criteria on the
prior for posterior consistency when observations form a time series, and prove
that any non-trivial prior is inconsistent when all observations are
contemporaneous. We then show that the likelihood given a data set of size $n
\in \mathbb{N}$ is constant across $\Lambda$-measures whose leading $n - 2$
moments agree, and focus on inferring truncated sequences of moments. We
provide a large class of functionals which can be extremised using finite
computation given a credible region of posterior truncated moment sequences,
and a pseudo-marginal Metropolis-Hastings algorithm for sampling the posterior.
Finally, we compare the efficiency of the exact and noisy pseudo-marginal
algorithms with and without delayed acceptance acceleration using a simulation
study.


  Routine public health surveillance of notifiable infectious diseases gives
rise to weekly counts of reported cases -- possibly stratified by region and/or
age group. We investigate how an age-structured social contact matrix can be
incorporated into a spatio-temporal endemic-epidemic model for infectious
disease counts. To illustrate the approach, we analyze the spread of norovirus
gastroenteritis over 6 age groups within the 12 districts of Berlin, 2011-2015,
using contact data from the POLYMOD study. The proposed age-structured model
outperforms alternative scenarios with homogeneous or no mixing between age
groups. An extended contact model suggests a power transformation of the
survey-based contact matrix towards more within-group transmission.


  Network inference is a rapidly advancing field, with new methods being
proposed on a regular basis. Understanding the advantages and limitations of
different network inference methods is key to their effective application in
different circumstances. The common structural properties shared by diverse
networks naturally pose a challenge when it comes to devising accurate
inference methods, but surprisingly, there is a paucity of comparison and
evaluation methods. Historically, every new methodology has only been tested
against \textit{gold standard} (true values) purpose-designed synthetic and
real-world (validated) biological networks. In this paper we aim to assess the
impact of taking into consideration aspects of topological and information
content in the evaluation of the final accuracy of an inference procedure.
Specifically, we will compare the best inference methods, in both
graph-theoretic and information-theoretic terms, for preserving topological
properties and the original information content of synthetic and biological
networks. New methods for performance comparison are introduced by borrowing
ideas from gene set enrichment analysis and by applying concepts from
algorithmic complexity. Experimental results show that no individual algorithm
outperforms all others in all cases, and that the challenging and non-trivial
nature of network inference is evident in the struggle of some of the
algorithms to turn in a performance that is superior to random guesswork.
Therefore special care should be taken to suit the method to the purpose at
hand. Finally, we show that evaluations from data generated using different
underlying topologies have different signatures that can be used to better
choose a network reconstruction method.


  Proteins fold using a two-state or multi-state kinetic mechanisms, but up to
now there isn't a first-principle model to explain this different behaviour. We
exploit the network properties of protein structures by introducing novel
observables to address the problem of classifying the different types of
folding kinetics. These observables display a plain physical meaning, in terms
of vibrational modes, possible configurations compatible with the native
protein structure, and folding cooperativity. The relevance of these
observables is supported by a classification performance up to 90%, even with
simple classifiers such as disciminant analysis.


  Established populations often exhibit oscillations in their sizes. If a
population is isolated, intrinsic stochasticity of elemental processes can
ultimately bring it to extinction. Here we study extinction of oscillating
populations in a stochastic version of the Rosenzweig-MacArthur predator-prey
model. To this end we extend a WKB approximation (after Wentzel, Kramers and
Brillouin) of solving the master equation to the case of extinction from a
limit cycle in the space of population sizes. We evaluate the extinction rates
and find the most probable paths to extinction by applying Floquet theory to
the dynamics of an effective WKB Hamiltonian. We show that the entropic
barriers to extinction change in a non-analytic way as the system passes
through the Hopf bifurcation. We also study the subleading pre-exponential
factors of the WKB approximation.


  Comprehensive and predictive simulation of coupled reaction networks has long
been a goal of biology and other fields. Currently, metabolic network models
that utilize enzyme mass action kinetics have predictive power but are limited
in scope and application by the fact that the determination of enzyme rate
constants is laborious and low throughput. We present a statistical
thermodynamic formulation of the law of mass action for coupled reactions at
both steady states and non-stationary states. The formulation is based on a
fluctuation theorem for coupled reactions and uses chemical potentials instead
of rate constants. When used to model deterministic systems, the theory
corresponds to a rescaling of the time dependent reactions in such a way that
steady states can be reached on the same time scale but with significantly
fewer computational steps. The significance for applications in systems biology
is discussed.


  A better characterization of the early growth dynamics of an epidemic is
needed to dissect the important drivers of disease transmission. We introduce a
2-parameter generalized-growth model to characterize the ascending phase of an
outbreak and capture epidemic profiles ranging from sub-exponential to
exponential growth. We test the model against empirical outbreak data
representing a variety of viral pathogens and provide simulations highlighting
the importance of sub-exponential growth for forecasting purposes. We applied
the generalized-growth model to 20 infectious disease outbreaks representing a
range of transmission routes. We uncovered epidemic profiles ranging from very
slow growth (p=0.14 for the Ebola outbreak in Bomi, Liberia (2014)) to near
exponential (p>0.9 for the smallpox outbreak in Khulna (1972), and the 1918
pandemic influenza in San Francisco). The foot-and-mouth disease outbreak in
Uruguay displayed a profile of slower growth while the growth pattern of the
HIV/AIDS epidemic in Japan was approximately linear. The West African Ebola
epidemic provided a unique opportunity to explore how growth profiles vary by
geography; analysis of the largest district-level outbreaks revealed
substantial growth variations (mean p=0.59, range: 0.14-0.97). Our findings
reveal significant variation in epidemic growth patterns across different
infectious disease outbreaks and highlights that sub-exponential growth is a
common phenomenon. Sub-exponential growth profiles may result from
heterogeneity in contact structures or risk groups, reactive behavior changes,
or the early onset of interventions strategies, and consideration of
"deceleration parameters" may be useful to refine existing mathematical
transmission models and improve disease forecasts.


  Game theoretic tools are utilized to analyze a one-locus continuous selection
model of sex-specific meiotic drive by considering nonequivalence of the
viabilities of reciprocal heterozygotes that might be noticed at an imprinted
locus. The model draws attention to the role of viability selections of
different types to examine the stable nature of polymorphic equilibrium. A
bridge between population genetics and evolutionary game theory has been built
up by applying the concept of the Fundamental Theorem of Natural Selection. In
addition to pointing out the influences of male and female segregation ratios
on selection, configuration structure reveals some noted results, e.g.,
Hardy-Weinberg frequencies hold in replicator dynamics, occurrence of faster
evolution at the maximized variance fitness, existence of mixed Evolutionarily
Stable Strategy ($ESS$) in asymmetric games, the tending evolution to follow
not only a $1:1$ sex ratio but also a $1:1$ different alleles ratio at
particular gene locus. Through construction of replicator dynamics in the group
selection framework, our selection model introduces a redefining bases of game
theory to incorporate non-random mating where a mating parameter associated
with population structure is dependent on the social structure. Also, the model
exposes the fact that the number of polymorphic equilibria will depend on the
algebraic expression of population structure.


  Spiking activity from populations of neurons display causal interactions and
memory effects. Therefore, they are expected to show some degree of
irreversibility in time. Motivated by the spike train statistics, in this paper
we build a framework to quantify the degree of irreversibility of any maximum
entropy distribution. Our approach is based on the transfer matrix technique,
which enables us to find an homogeneous irreducible Markov chain that shares
the same maximum entropy measure. We provide relevant examples in the context
of spike train statistics


  We propose, develop, and analyze a mathematical model of intergroup attitudes
in a community that is divided between two distinct social groups (which may be
distinguished by religion, ethnicity, or some other socially distinguishing
factor). The model is based on very simple premises that are both intuitive and
justified by sociological research. We investigate the behavior of the model in
various special cases, for various model configurations. We discuss the
stability of the model, and the continuous or discontinuous dependence of model
behavior on various parameters. Finally, we discuss possible implications for
strategies to improve intergroup affinity, and to defuse tension and prevent
deterioration of intergroup relationships.


  Designing algorithms that generate networks with a given degree sequence
while varying both subgraph composition and distribution of subgraphs around
nodes is an important but challenging research problem. Current algorithms lack
control of key network parameters, the ability to specify to what subgraphs a
node belongs to, come at a considerable complexity cost or, critically, sample
from a limited ensemble of networks. To enable controlled investigations of the
impact and role of subgraphs, especially for epidemics, neuronal activity or
complex contagion, it is essential that the generation process be versatile and
the generated networks as diverse as possible. In this paper, we present two
new network generation algorithms that use subgraphs as building blocks to
construct networks preserving a given degree sequence. Additionally, these
algorithms provide control over clustering both at node and global level. In
both cases, we show that, despite being constrained by a degree sequence and
global clustering, generated networks have markedly different topologies as
evidenced by both subgraph prevalence and distribution around nodes, and
large-scale network structure metrics such as path length and betweenness
measures. Simulations of standard epidemic and complex contagion models on
those networks reveal that degree distribution and global clustering do not
always accurately predict the outcome of dynamical processes taking place on
them. We conclude by discussing the benefits and limitations of both methods.


  Modeling of biomolecular systems plays an essential role in understanding
biological processes, such as ionic flow across channels, protein modification
or interaction, and cell signaling. The continuum model described by the
Poisson-Boltzmann (PB)/Poisson-Nernst-Planck (PNP) equations has made great
contributions towards simulation of these processes. However, the model has
shortcomings in its commonly used form and cannot capture (or cannot accurately
capture) some important physical properties of biological systems. Considerable
efforts have been made to improve the continuum model to account for discrete
particle interactions and to make progress in numerical methods to provide
accurate and efficient simulation. This review will summarize recent main
improvements in continuum modeling for biomolecular systems, with focus on the
size-modified models, the coupling of the classical density functional theory
and PNP equations, the coupling of polar and nonpolar interactions, and
numerical progress.


  Photosymbiosis is one of the most important evolutionary trajectories,
resulting in the chloroplast and the subsequent development of all complex
photosynthetic organisms. The ciliate Paramecium bursaria and the alga
Chlorella have a well established and well studied light dependent
endosymbiotic relationship. Despite its prominence there remain many unanswered
questions regarding the exact mechanisms of the photosymbiosis. Of particular
interest is how a host maintains and manages its symbiont load in response to
the allocation of nutrients between itself and its symbionts. Here we construct
a detailed mathematical model, parameterised from the literature, that
explicitly incorporates nutrient trading within a deterministic model of both
partners. The model demonstrates how the symbiotic relationship can manifest as
parasitism of the host by the symbionts, mutualism, wherein both partners
benefit, or exploitation of the symbionts by the hosts. We show that the
precise nature of the photosymbiosis is determined by both environmental
conditions (how much light is available for photosynthesis) and the level of
control a host has over its symbiont load. Our model provides a framework
within which it is possible to pose detailed questions regarding the
evolutionary behaviour of this important example of an established light
dependent endosymbiosis; we focus on one question in particular, namely the
evolution of host control, and show using an adaptive dynamics approach that a
moderate level of host control may evolve provided the associated costs are not
prohibitive.


  A resource selection function is a model of the likelihood that an available
spatial unit will be used by an animal, given its resource value. But how do we
appropriately define availability? Step-selection analysis deals with this
problem at the scale of the observed positional data, by matching each used
step (connecting two consecutive observed positions of the animal) with a set
of available steps randomly sampled from a distribution of observed steps or
their characteristics. Here we present a simple extension to this approach,
termed integrated step-selection analysis (iSSA), which relaxes the implicit
assumption that observed movement attributes (i.e. velocities and their
temporal autocorrelations) are independent of resource selection. Instead, iSSA
relies on simultaneously estimating movement and resource-selection parameters,
thus allowing simple likelihood-based inference of resource selection within a
mechanistic movement model. We provide theoretical underpinning of iSSA, as
well as practical guidelines to its implementation. Using computer simulations,
we evaluate the inferential and predictive capacity of iSSA compared to
currently used methods. Our work demonstrates the utility of iSSA as a general,
flexible and user-friendly approach for both evaluating a variety of ecological
hypotheses, and predicting future ecological patterns.


  Identifying genetic variants that regulate quantitative traits, or QTLs, is
the primary focus of the field of statistical genetics. Most current methods
are limited to identifying mean effects, or associations between genotype and
the mean value of a quantitative trait. It is possible, however, that a genetic
variant may affect the variance of the quantitative trait in lieu of, or in
addition to, affecting the trait mean. Here, we develop a general
methodological approach to identifying covariates with variance effects on a
quantitative trait using a Bayesian heteroskedastic linear regression model. We
show that our Bayesian test for heteroskedasticity (BTH) outperforms classical
tests for differences in variation across a large range of simulations drawn
from scenarios common to the analysis of quantitative traits. We apply BTH to
methylation QTL study data and expression QTL study data to identify variance
QTLs. When compared with three tests for heteroskedasticity used in genomics,
we illustrate the benefits of using our approach, including avoiding
overfitting by incorporating uncertainty and flexibly identifying
heteroskedastic effects.


  Locating recombination hotspots in genomic data is an important but difficult
task. Current methods frequently rely on estimating complicated models at high
computational cost. In this paper we develop an extremely fast, scalable method
for inferring recombination hot spots in a population of genomic sequences that
is based on the singular value decomposition. Our method performs well in
several synthetic data scenarios. We also apply our technique to a real data
investigation of the evolution of drug therapy resistance in a population of
HIV genomic sequences. Finally, we compare our method both on real and
simulated data to a state of the art algorithm.


  Social norms have traditionally been difficult to quantify. In any particular
society, their sheer number and complex interdependencies often limit a
system-level analysis. One exception is that of the network of norms that
sustain the online Wikipedia community. We study the fifteen-year evolution of
this network using the interconnected set of pages that establish, describe,
and interpret the community's norms. Despite Wikipedia's reputation for
\textit{ad hoc} governance, we find that its normative evolution is highly
conservative. The earliest users create norms that both dominate the network
and persist over time. These core norms govern both content and interpersonal
interactions using abstract principles such as neutrality, verifiability, and
assume good faith. As the network grows, norm neighborhoods decouple
topologically from each other, while increasing in semantic coherence. Taken
together, these results suggest that the evolution of Wikipedia's norm network
is akin to bureaucratic systems that predate the information age.


  The sequence of amino acid monomers in the primary structure of a protein is
decided by the corresponding sequence of codons (triplets of nucleic acid
monomers) on the template messenger RNA (mRNA). The polymerization of a
protein, by incorporation of the successive amino acid monomers, is carried out
by a molecular machine called ribosome. We develop a stochastic kinetic model
that captures the possibilities of mis-reading of mRNA codon and prior
mis-charging of a tRNA. By a combination of analytical and numerical methods we
obtain the distribution of the times taken for incorporation of the successive
amino acids in the growing protein in this mathematical model. The
corresponding exact analytical expression for the average rate of elongation of
a nascent protein is a `biologically motivated' generalization of the {\it
Michaelis-Menten formula} for the average rate of enzymatic reactions. This
generalized Michaelis-Menten-like formula (and the exact analytical expressions
for a few other quantities) that we report here display the interplay of four
different branched pathways corresponding to selection of four different types
of tRNA.


  Motivation: Single Molecule Real-Time (SMRT) sequencing technology and Oxford
Nanopore technologies (ONT) produce reads over 10kbp in length, which have
enabled high-quality genome assembly at an affordable cost. However, at
present, long reads have an error rate as high as 10-15%. Complex and
computationally intensive pipelines are required to assemble such reads.
  Results: We present a new mapper, minimap, and a de novo assembler, miniasm,
for efficiently mapping and assembling SMRT and ONT reads without an error
correction stage. They can often assemble a sequencing run of bacterial data
into a single contig in a few minutes, and assemble 45-fold C. elegans data in
9 minutes, orders of magnitude faster than the existing pipelines. We also
introduce a pairwise read mapping format (PAF) and a graphical fragment
assembly format (GFA), and demonstrate the interoperability between ours and
current tools.
  Availability and implementation: https://github.com/lh3/minimap and
https://github.com/lh3/miniasm
  Contact: hengli@broadinstitute.org


  Faithful segregation of genetic material during cell division requires
alignment of chromosomes between two spindle poles and attachment of their
kinetochores to each of the poles. Failure of these complex dynamical processes
leads to chromosomal instability (CIN), a characteristic feature of several
diseases including cancer. While a multitude of biological factors regulating
chromosome congression and bi-orientation have been identified, it is still
unclear how they are integrated so that coherent chromosome motion emerges from
a large collection of random and deterministic processes. Here we address this
issue by a three dimensional computational model of motor-driven chromosome
congression and bi-orientation during mitosis. Our model reveals that
successful cell division requires control of the total number of microtubules:
if this number is too small bi-orientation fails, while if it is too large not
all the chromosomes are able to congress. The optimal number of microtubules
predicted by our model compares well with early observations in mammalian cell
spindles. Our results shed new light on the origin of several pathological
conditions related to chromosomal instability.


  We present a short introduction into the framework of piecewise deterministic
Markov processes. We illustrate the abstract mathematical setting with a series
of examples related to dispersal of biological systems, cell cycle models, gene
expression, physiologically structured populations, as well as neural activity.
General results concerning asymptotic properties of stochastic semigroups
induced by such Markov processes are applied to specific examples.


  Living cells can enhance their fitness by anticipating environmental change.
We study how accurately linear signaling networks in cells can predict future
signals. We find that maximal predictive power results from a combination of
input-noise suppression, linear extrapolation, and selective readout of
correlated past signal values. Single-layer networks generate exponential
response kernels, which suffice to predict Markovian signals optimally.
Multilayer networks allow oscillatory kernels that can optimally predict
non-Markovian signals. At low noise, these kernels exploit the signal
derivative for extrapolation, while at high noise, they capitalize on signal
values in the past that are strongly correlated with the future signal. We show
how the common motifs of negative feedback and incoherent feed-forward can
implement these optimal response functions. Simulations reveal that E. coli can
reliably predict concentration changes for chemotaxis, and that the integration
time of its response kernel arises from a trade-off between rapid response and
noise suppression.


  Genome-wide association studies have proven to be essential for understanding
the genetic basis of disease. However, many complex traits---personality
traits, facial features, disease subtyping---are inherently high-dimensional,
impeding simple approaches to association mapping. We developed a nonparametric
Bayesian reduced rank regression model for multi-SNP, multi-trait association
mapping that does not require the rank of the linear subspace to be specified.
We show in simulations and real data that our model shares strength over SNPs
and over correlated traits, improving statistical power to identify genetic
associations with an interpretable, SNP-supervised low-dimensional linear
projection of the high-dimensional phenotype. On the HapMap phase 3 gene
expression QTL study data, we identify pleiotropic expression QTLs that
classical univariate tests are underpowered to find and that two step
approaches cannot recover. Our Python software, BERRRI, is publicly available
at GitHub: https://github.com/ashlee1031/BERRRI.


  Shannon information has, in the past, been applied to quantify the genetic
diversity of many natural populations. Here, we apply the Shannon concept to
consecutive generations of alleles as they evolve over time. We suppose a
genetic system analogous to the discrete noisy channel of Shannon, where the
signal emitted by the input (mother population) is a number of alleles that
will form the next generation (offspring). The alleles received at a given
generation are conditional upon the previous generation. Knowledge of this
conditional probability law allows us to track the evolution of the allele
entropies and mutual information values from one generation to the next. We
apply these laws to numerical computer simulations and to real data
(Stryphnodendron adstringens). We find that, due to the genetic sampling
process, in the absence of new mutations the mutual information increases
between generations toward a maximum value. Lastly, after sufficient
generations the system has a level of mutual information equal to the entropy
(diversity) that it had at the beginning of the process (mother population).
This implies no increase in genetic diversity in the absence of new mutations.
Now, obviously, mutations are essential to the evolution of species. In
addition, we observe that when a population shows at least a low level of
genetic diversity, the highest values of mutual information between the
generations occurs when the system is neither too orderly nor too disorderly.
We also find that mutual information is a valid measure of allele fixation.


  A tree-based network on a set $X$ of $n$ leaves is said to be universal if
any rooted binary phylogenetic tree on $X$ can be its base tree. Francis and
Steel showed that there is a universal tree-based network on $X$ in the case of
$n=3$, and asked whether such a network exists in general. We settle this
problem by proving that there are infinitely many universal tree-based networks
for any $n>1$.


  The dynamics of many microbial ecosystems are driven by cross-feeding
interactions, in which metabolites excreted by some species are metabolised
further by others. The population dynamics of such ecosystems are governed by
frequency-dependent selection, which allows for stable coexistence of two or
more species. We have analysed a model of cross-feeding based on the replicator
equation, with the aim of establishing criteria for coexistence in ecosystems
containing three species, given the information of the three species' ability
to coexist in their three separate pairs, i.e. the long term dynamics in the
three two-species component systems. The triple-system is studied statistically
and the probability of coexistence in the species triplet is computed for two
models of species interactions. The interaction parameters are modelled either
as stochastically independent or organised in a hierarchy where any derived
metabolite carries less energy than previous nutrients in the metabolic chain.
We differentiate between different modes of coexistence with respect to the
pair-wise dynamics of the species, and find that the probability of coexistence
is close to $\frac12$ for triplet systems with three pair-wise coexistent pairs
and for so-called intransitive systems. Systems with two and one pair-wise
coexistent pairs are more likely to exist for random interaction parameters,
but are on the other hand much less likely to exhibit triplet coexistence.
Hence we conclude that certain species triplets are, from a statistical point
of view, rare, but if allowed to interact are likely to coexist. This knowledge
might be helpful when constructing synthetic microbial communities for
industrial purposes.


  Reanalysis of data (Bruice and York 1961) for the pH dependences of the
calculated apparent second-order rate constants (k2) for the reaction of
tris-(hydroxymethyl)-aminomethane (TRIS) with phenyl esters reveals that the pH
profiles are consistently much better represented as multiphasic, i.e. as a
series of straight lines separated by discontinuities in the form of sharp
breaks or noncontiguities (jumps), than by the continuous function used by the
authors. For p-nitrophenyl acetate, m-nitrophenyl acetate, p-chlorophenyl
acetate and phenyl acetate, the profiles are best represented as multiphasic in
plots of log k2 vs. pH, whereas the profile for p-methylphenyl acetate is only
multiphasic in a plot of k2 vs. pH. Data for the pseudo-first-order rate
constants (kobs) give multiphasic pH profiles with patterns identical or very
similar to those for k2. Extensive reanalyses of data for a wide variety of
enzymatic and non-enzymatic reactions show, similarly, that pH profiles are in
general better represented as multiphasic than as curvilinear.


  Mathematical and computational modelling of biochemical networks is often
done in terms of either the concentrations of molecular species or the fluxes
of biochemical reactions. When is mathematical modelling from either
perspective equivalent to the other? Mathematical duality translates concepts,
theorems or mathematical structures into other concepts, theorems or
structures, in a one-to-one manner. We present a novel stoichiometric condition
that is necessary and sufficient for duality between unidirectional fluxes and
concentrations. Our numerical experiments, with computational models derived
from a range of genome-scale biochemical networks, suggest that this
flux-concentration duality is a pervasive property of biochemical networks. We
also provide a combinatorial characterisation that is sufficient to ensure
flux-concentration duality. That is, for every two disjoint sets of molecular
species, there is at least one reaction complex that involves species from only
one of the two sets. When unidirectional fluxes and molecular species
concentrations are dual vectors, this implies that the behaviour of the
corresponding biochemical network can be described entirely in terms of either
concentrations or unidirectional fluxes.


  The evolution of species habitat range is an important topic over a wide
range of research fields. In higher organisms, habitat range evolution is
generally associated with genetic events such as gene duplication. However, the
specific factors that determine habitat variability remain unclear at higher
levels of biological organization (e.g., biochemical networks). One widely
accepted hypothesis developed from both theoretical and empirical analyses is
that habitat variability promotes network modularity; however, this
relationship has not yet been directly tested in higher organisms. Therefore, I
investigated the relationship between habitat variability and metabolic network
modularity using compound and enzymatic networks in flies and mammals. Contrary
to expectation, there was no clear positive correlation between habitat
variability and network modularity. As an exception, the network modularity
increased with habitat variability in the enzymatic networks of flies. However,
the observed association was likely an artifact, and the frequency of gene
duplication appears to be the main factor contributing to network modularity.
These findings raise the question of whether or not there is a general
mechanism for habitat range expansion at a higher level (i.e., above the gene
scale). This study suggests that the currently widely accepted hypothesis for
habitat variability should be reconsidered.


  An analytical expression for the optimal control of an Ebola problem is
obtained. The analytical solution is found as a first-order approximation to
the Pontryagin Maximum Principle via the Euler-Lagrange equation. An
implementation of the method is given using the computer algebra system Maple.
Our analytical solutions confirm the results recently reported in the
literature using numerical methods.


  The number of molecules involved in a cell or subcellular structure is
sometimes rather small. In this situation, ordinary macroscopic-level
fluctuations can be overwhelmed by non-negligible large fluctuations, which
results in drastic changes in chemical-reaction dynamics and statistics
compared to those observed under a macroscopic system (i.e., with a large
number of molecules). In order to understand how salient changes emerge from
fluctuations in molecular number, we here quantitatively define small-number
effect by focusing on a `mesoscopic' level, in which the concentration
distribution is distinguishable both from micro- and macroscopic ones, and
propose a criterion for determining whether or not such an effect can emerge in
a given chemical reaction network. Using the proposed criterion, we
systematically derive a list of motifs of chemical reaction networks that can
show small-number effects, which includes motifs showing emergence of the power
law and the bimodal distribution observable in a mesoscopic regime with respect
to molecule number. The list of motifs provided herein is helpful in the search
for candidates of biochemical reactions with a small-number effect for possible
biological functions, as well as for designing a reaction system whose behavior
can change drastically depending on molecule number, rather than concentration.


  Dangerous damage to mitochondrial DNA (mtDNA) can be ameliorated during
mammalian development through a highly debated mechanism called the mtDNA
bottleneck. Uncertainty surrounding this process limits our ability to address
inherited mtDNA diseases. We produce a new, physically motivated, generalisable
theoretical model for mtDNA populations during development, allowing the first
statistical comparison of proposed bottleneck mechanisms. Using approximate
Bayesian computation and mouse data, we find most statistical support for a
combination of binomial partitioning of mtDNAs at cell divisions and random
mtDNA turnover, meaning that the debated exact magnitude of mtDNA copy number
depletion is flexible. New experimental measurements from a wild-derived mtDNA
pairing in mice confirm the theoretical predictions of this model. We
analytically solve a mathematical description of this mechanism, computing
probabilities of mtDNA disease onset, efficacy of clinical sampling strategies,
and effects of potential dynamic interventions, thus developing a quantitative
and experimentally-supported stochastic theory of the bottleneck.


  In this work we formalize a method to compute the degree of importance of
biological agents that participates on the dynamics of a biological phenomenon
build upon a complex network. We call this new procedure by theoretical
knock-out (KO). To devise this method, we make two approaches: algebraically
and algorithmically. In both cases we compute a vector on an asymptotic state,
called flux vector. The flux is given by a random walk on a directed graph that
represents a biological phenomenon. This vector gives us the information about
the relative flux of walkers on a vertex which represents a biological agent.
With two vector of this kind, we can calculate the relative mean error between
them by averaging over its coefficients. This quantity allows us to assess the
degree of importance of each vertex of a complex network that evolves in time
and has experimental background. We find out that this procedure can be applied
in any sort of biological phenomena in which we can know the role and
interrelationships of its agents. These results also provide experimental
biologists to predict the order of importance of biological agents on a mounted
complex network.


  Spatial stochastic molecular simulations in biology are limited by the
intense computation required to track molecules in space either in a discrete
time or discrete space framework, meaning that the serial limit has already
been reached in sub-cellular models. This calls for parallel simulations that
can take advantage of the power of modern supercomputers; however exact methods
are known to be inherently serial. We introduce an operator splitting
implementation for irregular grids with a novel method to improve accuracy, and
demonstrate potential for scalable parallel simulations in an initial MPI
version. We foresee that this groundwork will enable larger scale, whole-cell
stochastic simulations in the near future.


  It is well established that amyloid fibril solubility is protein specific,
but how solubility depends on the interactions between the fibril building
blocks is not clear. Here we use a simple protein model and perform Monte Carlo
simulations to directly measure the solubility of amyloid fibrils as a function
of the interaction between the fibril building blocks. Our simulations confirms
that the fibril solubility depends on the fibril thickness and that the
relationship between the interactions and the solubility can be described by a
simple analytical formula. The results presented in this study reveal general
rules how side-chain side-chain interactions, backbone hydrogen bonding and
temperature affect amyloid fibril solubility, which might prove a powerful tool
to design protein fibrils with desired solubility and aggregation properties in
general.


  World population and the number of cultural artifacts are growing
exponentially or faster, while cultural interaction approaches the fidelity of
a global nervous system. Every day hundreds of millions of images are loaded
into social networks by users all over the world. As this myriad of new
artifacts veils the view into the past, like city lights covering the night
sky, it is easy to forget that there is more than one Starry Night, the
painting by Van Gogh. Like in ecology, where saving rare species may help us in
treating disease, art and architectural history can reveal insights into the
past, which may hold keys to our own future. With humanism under threat, facing
the challenge of understanding the structure and dynamics of art and culture,
both qualitatively and quantitatively, is more crucial now than it ever was.
The purpose of this article is to provide perspective in the aim of figuring
out the process of art history - not art history as a discipline, but the
actual history of all made things, in the spirit of George Kubler and Marcel
Duchamp. In other words, this article deals with the grand challenge of
developing a systematic science of art and culture, no matter what, and no
matter how.


  Recent advances in next-generation sequencing have revolutionized genomic
research. 16S rRNA amplicon sequencing using paired-end sequencing on the MiSeq
platform from Illumina, Inc., is being used to characterize the composition and
dynamics of extremely complex/diverse microbial communities. For this analysis
on the Illumina platform, merging and quality filtering of paired-end reads are
essential first steps in data analysis to ensure the accuracy and reliability
of downstream analysis. We have developed the Merging and Filtering Tool
(MeFiT) to combine these pre-processing steps into one simple, intuitive
pipeline. MeFiT provides an open-source solution that permits users to merge
and filter paired end illumina reads based on user-selected quality parameters.
The tool has been implemented in python and the source-code is freely available
at https://github.com/nisheth/MeFiT.


  Next generation sequencing technology rapidly produces massive volume of data
and quality control of this sequencing data is essential to any genomic
analysis. Here we present MEEPTOOLS, which is a collection of open-source tools
based on maximum expected error as a percentage of read length (MEEP score) to
filter, trim, truncate and assess next generation DNA sequencing data in FASTQ
file format. MEEPTOOLS provides a non-traditional approach towards read
filtering/trimming based on maximum error probabilities of the bases in the
read on a non-logarithmic scale. This method simultaneously retains more
reliable bases and removes more unreliable bases than the traditional quality
filtering strategies.


  The spread of Ebola virus in 2014 is unprecedented. The epidemic is still
affecting West Africa, exacerbated by extraordinary socioeconomic disadvantages
and health system inadequacies. With the aim of understanding, predicting, and
control the propagation of the virus in the populations of affected countries,
it is crucial to model the dynamics of the virus and study several strategies
to control it. In this paper, we present a very simple mathematical model that
describes quite well the spread of Ebola. Then, we discuss several strategies
for the control of the propagation of this lethal virus into populations, in
order to predict the impact of vaccine programmes, treatment, and the impact of
educational campaigns.


  We present three major transitions that occur on the way to the elaborate and
diverse societies of the modern era. Our account links the worlds of social
animals such as pigtail macaques and monk parakeets to examples from human
history, including 18th Century London and the contemporary online phenomenon
of Wikipedia. From the first awareness and use of group-level social facts to
the emergence of norms and their self-assembly into normative bundles, each
transition represents a new relationship between the individual and the group.
At the center of this relationship is the use of coarse-grained information
gained via lossy compression. The role of top-down causation in the origin of
society parallels that conjectured to occur in the origin and evolution of life
itself.


  Paper provides biological description of both males and females of Arctic
lamprey, Lethenteron camtschaticum ammocoetes from one of Sakhalin rivers,
Malaya Khuzi. It briefly describes stages of development of the larvae gonads
and cases of infantilism and transverse.


  A method was developed for Bayesian inference of species phylogeny using the
multi-species coalescent model. To improve the mixing properties of the Markov
chain Monte Carlo (MCMC) algorithm that traverses the space of species trees,
we implement two efficient MCMC proposals: the first is based on the Subtree
Pruning and Regrafting (SPR) algorithm and the second is based on a novel
node-slider algorithm. Like the Nearest-Neighbor Interchange (NNI) algorithm we
implemented previously, both algorithms propose changes to the species tree,
while simultaneously altering the gene trees at multiple genetic loci to
automatically avoid conflicts with the newly-proposed species tree. The method
integrates over gene trees, naturally taking account of the uncertainty of gene
tree topology and branch lengths given the sequence data. A simulation study
was performed to examine the statistical properties of the new method. We found
that it has excellent statistical performance, inferring the correct species
tree with near certainty when analyzing 10 loci. The prior on species trees has
some impact, particularly for small numbers of loci. An empirical dataset (for
rattlesnakes) was reanalyzed. While the 18 nuclear loci and one mitochondrial
locus support largely consistent species trees under the multi-species
coalescent model estimates of parameters suggest drastically different
evolutionary dynamics between the nuclear and mitochondrial loci.


  In 1938, H. Selye proposed the notion of adaptation energy and published
"Experimental evidence supporting the conception of adaptation energy".
Adaptation of an animal to different factors appears as the spending of one
resource. Adaptation energy is a hypothetical extensive quantity spent for
adaptation. This term causes much debate when one takes it literally, as a
physical quantity, i.e. a sort of energy. The controversial points of view
impede the systematic use of the notion of adaptation energy despite
experimental evidence. Nevertheless, the response to many harmful factors often
has general non-specific form and we suggest that the mechanisms of
physiological adaptation admit a very general and nonspecific description.
  We aim to demonstrate that Selye's adaptation energy is the cornerstone of
the top-down approach to modelling of non-specific adaptation processes. We
analyse Selye's axioms of adaptation energy together with Goldstone's
modifications and propose a series of models for interpretation of these
axioms. {\em Adaptation energy is considered as an internal coordinate on the
`dominant path' in the model of adaptation}. The phenomena of `oscillating
death' and `oscillating remission' are predicted on the base of the dynamical
models of adaptation. Natural selection plays a key role in the evolution of
mechanisms of physiological adaptation. We use the fitness optimization
approach to study of the distribution of resources for neutralization of
harmful factors, during adaptation to a multifactor environment, and analyse
the optimal strategies for different systems of factors.


  Over the past century, nonlinear difference and differential equations have
been used to understand conditions for species coexistence. However, these
models fail to account for random fluctuations due to demographic and
environmental stochasticity which are experienced by all populations. I review
some recent mathematical results about persistence and coexistence for models
accounting for each of these forms of stochasticity. Demographic stochasticity
stems from populations and communities consisting of a finite number of
interacting individuals, and often are represented by Markovian models with a
countable number of states. For closed populations in a bounded world,
extinction occurs in finite time but may be preceded by long-term transients.
Quasi-stationary distributions (QSDs) of these Makov models characterize this
meta-stable behavior. For sufficiently large "habitat sizes", QSDs are shown to
concentrate on the positive attractors of deterministic models. Moreover, the
probability extinction decreases exponentially with habitat size.
Alternatively, environmental stochasticity stems from fluctuations in
environmental conditions which influence survival, growth, and reproduction.
Stochastic difference equations can be used to model the effects of
environmental stochasticity on population and community dynamics. For these
models, stochastic persistence corresponds to empirical measures placing
arbitrarily little weight on arbitrarily low population densities. Sufficient
and necessary conditions for stochastic persistence are reviewed. These
conditions involve weighted combinations of Lyapunov exponents corresponding to
"average" per-capita growth rates of rare species. The results are illustrated
with Bay checkerspot butterflies, coupled sink populations, the storage effect,
and stochastic rock-paper-scissor communities. Open problems and conjectures
are presented.


  Stochasticity plays a fundamental role in various biochemical processes, such
as cell regulatory networks and enzyme cascades. Isothermal, well-mixed systems
can be modelled as Markov processes, typically simulated using the Gillespie
Stochastic Simulation Algorithm (SSA). While easy to implement and exact, the
computational cost of using the Gillespie SSA to simulate such systems can
become prohibitive as the frequency of reaction events increases. This has
motivated numerous coarse-grained schemes, where the "fast" reactions are
approximated either using Langevin dynamics or deterministically. While such
approaches provide a good approximation when all reactants are abundant, the
approximation breaks down when one or more species exist only in small
concentrations and the fluctuations arising from the discrete nature of the
reactions becomes significant. This is particularly problematic when using such
methods to compute statistics of extinction times for chemical species, as well
as simulating non-equilibrium systems such as cell-cycle models in which a
single species can cycle between abundance and scarcity. In this paper, a
hybrid jump-diffusion model for simulating well- mixed stochastic kinetics is
derived. It acts as a bridge between the Gillespie SSA and the chemical
Langevin equation. For low reactant reactions the underlying behaviour is
purely discrete, while purely diffusive when the concentrations of all species
is large, with the two different behaviours coexisting in the intermediate
region. A bound on the weak error in the classical large volume scaling limit
is obtained, and three different numerical discretizations of the
jump-diffusion model are described. The benefits of such a formalism are
illustrated using computational examples.


  We describe a robust, fast, and memory-efficient procedure that can cluster
millions of structures derived from molecular dynamics simulations. The essence
of the method is based on a peak-picking algorithm applied to three- and
five-dimensional distributions of the principal components derived from the
trajectory and automatically supports both Cartesian and dihedral PCA-based
clustering. The density threshold required for identifying isolated peaks
(which correspond to discrete clusters) is determined through the application
of a variance-explained criterion which allows for a completely automated
clustering procedure with no user intervention. In this communication we
describe the algorithm and present some of the results obtained from the
application of the method as implemented in the molecular dynamics analysis
programs carma, grcarma. and cluster5D. We conclude with a discussion of the
limitations and possible pitfalls of this method.


  It is well stablished that danger information can be transmitted by ants
through relatively small distances, provoking either a state of alarm when they
move away from potentially dangerous stimulus, or charge toward it
aggressively. There is almost no knowledge if danger information can be
transmitted along large distances. In this paper, we perturb leaf cutting ants
of the species Atta insularis while they forage in their natural evioronment at
a certain point of the foraging line, so ants make a "U" turn to escape from
the danger zone and go back to the nest. Our results strongly suggest that
those ants do not transmit "danger information" to other nestmates marching
towards the danger area. The individualistic behavior of the ants returning
from the danger zone results in a depression of the foraging activity due to
the systematic sacrifice of non-informed individuals.


  What is the boundary between a vigorous argument and a breakdown of
relations? What drives a group of individuals across it? Taking Wikipedia as a
test case, we use a hidden Markov model to approximate the computational
structure and social grammar of more than a decade of cooperation and conflict
among its editors. Across a wide range of pages, we discover a bursty war/peace
structure where the systems can become trapped, sometimes for months, in a
computational subspace associated with significantly higher levels of
conflict-tracking "revert" actions. Distinct patterns of behavior characterize
the lower-conflict subspace, including tit-for-tat reversion. While a fraction
of the transitions between these subspaces are associated with top-down actions
taken by administrators, the effects are weak. Surprisingly, we find no
statistical signal that transitions are associated with the appearance of
particularly anti-social users, and only weak association with significant news
events outside the system. These findings are consistent with transitions being
driven by decentralized processes with no clear locus of control. Models of
belief revision in the presence of a common resource for information-sharing
predict the existence of two distinct phases: a disordered high-conflict phase,
and a frozen phase with spontaneously-broken symmetry. The bistability we
observe empirically may be a consequence of editor turn-over, which drives the
system to a critical point between them.


  A fundamental issue discussed in evolutionary biology is the transition from
unicellular to multicellular organisms. Here we develop non-robust models
provided in [1] and attempt to get robust models investigated how
differentiation of types and energy constraints influence on the optimal
behavior of colonies with different size. Constructed models show that each
large-sized colony with high initial costs of reproduction tends to full
specialization, no matter are all cells in this colony identical or are there
cells with different types in this colony. The level of type's diversity
determines the number of cells specialized in soma. In small - sized colonies
with low initial costs of reproduction, when type's diversity is week, an
unspecialized state may bring colony some benefits. However, these benefits may
be only local and in optimum in the colony some cells would be specialized,
others - unspecialized. The amount of specialized cells in small - sized colony
depends on the level of type's diversity in this colony. Adding energy
constraint, we may receive robust models even in convex case. In optimum, the
colony with different types of cells and energy restriction may be indifferent
between some optimal patterns of states. Arbitrary chosen cell may be soma or
germ in some states or may be unspecialized in other. Moreover, in each optimal
state levels of fecundity and viability of each cell lie in limited ranges.
This result reflects the fact that some cell in the colony may lose the
potential ability to achieve, for example, high level of fecundity, but does
not lose the possibility to perform a reproductive function at all. It means
that provided model can describe organisms, which represents the intermediate
between unspecialized colonies and full-specialized multicellular organisms.


  We present here a system of self-propelled particles that follow a very
simple motion law in continuous space in a deterministic and asynchronous way.
This system of particles is capable of producing, depending on the particle
density in the habitat, several spatio-temporal patterns emerging from an
initial randomized spatial configuration. We found that those structures show
specific population dynamics which arise from death (decay) and growth
(self-replication) of those structures, thus we call the system Primordial
Particle System (PPS), as the model can be interpreted as a simplistic model of
emergence of self-replicating chemical structures from initially chaotic mixed
components in the "primordial soup" at the beginning of life. We describe the
observed dynamics, show the emerging spatio-temporal structures and present a
macroscopic top-down model as well as a probabilistic microscopic bottom-up
model of the system.


  Larvae of the blackspotted pliers support beetle, Rhagium mordax, express
antifreeze proteins in their haemolymph during temperate climate winter. It is
believed that they also express antifreeze proteins in their cuticle as a means
of preventing inoculative freezing. Larvae of Rhagium mordax were collected
during winter (March) and summer (May) of 2011. Larvae were fixated, embedded
in paraffin wax, sectioned on a microtome, incubated with custom made anti-AFP
antibodies and isualised on a fluorescence microscope. The larvae of both
winter and summer showed AFP activity in their cuticle, gut lumen and
-epithelium. Due to the long synthesis process of AFPs, the larvae contain them
all year round. The distribution of these AFPs change during summer, possibly
relocating to vesicles in the cuticle and gut lumen/epithelium.


  A cell-molecular based evolutionary model of tumor development driven by a
stochastic Moran birth-death process is developed, where each cell carries
molecular information represented by a four-digit binary string, used to
differentiate cells into 16 molecular types. The binary string value determines
cell fitness, with lower fit cells (e.g. 0000) defined as healthy phenotypes,
and higher fit cells (e.g. 1111) defined as malignant phenotypes. At each step
of the birth-death process, the two phenotypic sub-populations compete in a
prisoner's dilemma evolutionary game with healthy cells (cooperators) competing
with cancer cells (defectors). Fitness and birth-death rates are defined via
the prisoner's dilemma payoff matrix. Cells are able undergo two types of
stochastic point mutations passed to the daughter cell's binary string during
birth: passenger mutations (conferring no fitness advantage) and driver
mutations (increasing cell fitness). Dynamic phylogenetic trees show clonal
expansions of cancer cell sub-populations from an initial malignant cell. The
tumor growth equation states that the growth rate is proportional to the
logarithm of cellular heterogeneity, here measured using the Shannon entropy of
the distribution of binary sequences in the tumor cell population. Nonconstant
tumor growth rates, (exponential growth during sub-clinical range of the tumor
and subsequent slowed growth during tumor saturation) are associated with a
Gompertzian growth curve, an emergent feature of the model explained here using
simple statistical mechanics principles related to the degree of functional
coupling of the cell states. Dosing strategies at early stage development,
mid-stage (clinical stage), and late stage development of the tumor are
compared, showing therapy is most effective during the sub-clinical stage,
before the cancer subpopulation is selected for growth.


  Tumor development is an evolutionary process in which a heterogeneous
population of cells with differential growth capabilities compete for resources
in order to gain a proliferative advantage. What are the minimal ingredients
needed to recreate some of the emergent features of such a developing complex
ecosystem? What is a tumor doing before we can detect it? We outline a
mathematical model, driven by a stochastic Moran process, in which cancer cells
and healthy cells compete for dominance in the population. Each are assigned
payoffs according to a Prisoner's Dilemma evolutionary game where the healthy
cells are the cooperators and the cancer cells are the defectors. With point
mutational dynamics, heredity, and a fitness landscape controlling birth and
death rates, natural selection acts on the cell population and simulated
"cancer-like" features emerge, such as Gompertzian tumor growth driven by
heterogeneity, the log-kill law which (linearly) relates therapeutic dose
density to the (log) probability of cancer cell survival, and the Norton-Simon
hypothesis which (linearly) relates tumor regression rates to tumor growth
rates. We highlight the utility, clarity, and power that such models provide,
despite (and because of) their simplicity and built-in assumptions.


  Biology offers a boundless source of adaptation, innovation, and inspiration.
A wide range of photosynthetic organisms exist that are capable of harvesting
solar light in an exceptionally efficient way, using abundant and low-cost
materials. These natural light-harvesting complexes consist of proteins that
strongly bind a high density of chromophores to capture solar photons and
rapidly transfer the excitation energy to the photochemical reaction centre.
The amount of harvested light is also delicately tuned to the level of solar
radiation to maintain a constant energy throughput at the reaction centre and
avoid the accumulation of the products of charge separation. In this Review,
recent developments in the understanding of light harvesting by plants will be
discussed, based on results obtained from single molecule spectroscopy studies.
Three design principles of the main light-harvesting antenna of plants will be
highlighted: (a) fine, photoactive control over the intrinsic protein disorder
to efficiently use intrinsically available thermal energy dissipation
mechanisms; (b) the design of the protein microenvironment of a low-energy
chromophore dimer to control the amount of shade absorption; (c) the design of
the exciton manifold to ensure efficient funneling of the harvested light to
the terminal emitter cluster.


  Certain dissipative physical systems closely resemble Hamiltonian systems in
$\mathbb{R}^{2n}$, but with the canonical equation for one of the variables in
each conjugate pair rescaled by a real parameter. To generalise these dynamical
systems to symplectic manifolds in this paper we introduce and study the
properties of deformed Hamiltonian vector fields on Lagrangian fibrations. We
describe why these objects have some interesting applications to symplectic
geometry and discuss how their physical interpretation motivates new problems
in mathematics.


  In population genetic studies, the allele frequency spectrum (AFS)
efficiently summarizes genome-wide polymorphism data and shapes a variety of
allele frequency-based summary statistics. While existing theory typically
features equilibrium conditions, emerging methodology requires an analytical
understanding of the build-up of the allele frequencies over time. In this
work, we use the framework of Poisson random fields to derive new
representations of the non-equilibrium AFS for the case of a Wright-Fisher
population model with selection. In our approach, the AFS is a scaling-limit of
the expectation of a Poisson stochastic integral and the representation of the
non-equilibrium AFS arises in terms of a fixation time probability
distribution. The known duality between the Wright-Fisher diffusion process and
a birth and death process generalizing Kingman's coalescent yields an
additional representation. The results carry over to the setting of a random
sample drawn from the population and provide the non-equilibrium behavior of
sample statistics. Our findings are consistent with and extend a previous
approach where the non-equilibrium AFS solves a partial differential forward
equation with a non-traditional boundary condition. Moreover, we provide a
bridge to previous coalescent-based work, and hence tie several frameworks
together. Since frequency-based summary statistics are widely used in
population genetics, for example, to identify candidate loci of adaptive
evolution, to infer the demographic history of a population, or to improve our
understanding of the underlying mechanics of speciation events, the presented
results are potentially useful for a broad range of topics.


  Phase diagram of the ground states of DNA in a bad solvent is studied for a
semi-flexible polymer model with a generalized local elastic bending potential
characterized by a nonlinearity parameter $x$ and effective self-attraction
promoting compaction. $x=1$ corresponds to the worm-like chain model.
Surprisingly, the phase diagram as well as the transition lines between the
ground states are found to be a function of $x$. The model provides a simple
explanation for the results of prior experimental and computational studies and
makes predictions for the specific geometries of the ground states. The results
underscore the impact of the form of the microscopic bending energy at
macroscopic observable scales.


  Integrating expression data with gene interactions in a network is essential
for understanding the functional organization of the cells. Consequently,
knowledge of interaction types in biological networks is important for data
interpretation. Signing of Regulatory Networks (SIREN) plugin for Cytoscape is
an open-source Java tool for discrimination of interaction type (activatory or
inhibitory) in gene regulatory networks. Utilizing an information theory based
concept, SIREN seeks to identify the interaction type of pairs of genes by
examining their corresponding gene expression profiles. We introduce SIREN, a
fast and memory efficient tool with low computational complexity, that allows
the user to easily consider it as a complementary approach for many network
reconstruction methods. SIREN allows biologists to use independent expression
data to predict interaction types for known gene regulatory networks where
reconstruction methods do not provide any information about the nature of their
interaction types. The SIREN Cytoscape plugin is implemented in Java and is
freely available at http://baderlab.org/Software/SIRENplugin and via the
Cytoscape app manager.


  A eukaryotic gene consists of multiple exons (protein coding regions) and
introns (non-coding regions), and a splice junction refers to the boundary
between a pair of exon and intron. Precise identification of spice junctions on
a gene is important for deciphering its primary structure, function, and
interaction. Experimental techniques for determining exon/intron boundaries
include RNA-seq, which is often accompanied by computational approaches.
Canonical splicing signals are known, but computational junction prediction
still remains challenging because of a large number of false positives and
other complications. In this paper, we exploit deep recurrent neural networks
(RNNs) to model DNA sequences and to detect splice junctions thereon. We test
various RNN units and architectures including long short-term memory units,
gated recurrent units, and recently proposed iRNN for in-depth design space
exploration. According to our experimental results, the proposed approach
significantly outperforms not only conventional machine learning-based methods
but also a recent state-of-the-art deep belief network-based technique in terms
of prediction accuracy.


  The Tangled Nature Model of biological and cultural evolution features
interacting agents which compete for limited resources and reproduce in an
error prone fashion and at a rate depending on the `tangle' of interactions
they maintain with others.
  The set of interactions linking a TNM individual to others is key to its
reproductive success and arguably constitutes its most important property. Yet,
in many studies, the interactions of an individual and those of its mutated
off-spring are unrelated, a rather unrealistic feature corresponding to a point
mutation turning a giraffe into an elephant. To bring out the structural and
dynamical effects of trait inheritance , we introduce and numerically analyze a
family of TNM models where a positive integer $K$ parametrises correlations
between the interactions of an agent and those of its mutated offspring. For
$K=1$ a single point mutation randomizes all the interactions, while increasing
$K$ up to the length of the genome ensures an increasing level of trait
inheritance. We show that the distribution of the interactions generated by our
rule is nearly independent of the value of $K$. Changing $K$ strengthens the
core structure of the ecology, leads to population abundance distributions
which are better approximated by log-normal probability densities and increases
the probability that a species extant at time $t_{\rm w}$ is also extant at a
later time $t$. In particular, survival probabilities are shown to decay as
powers of the ratio $t/t_{\rm w}$, similarity to the pure aging behaviour
approximately describing glassy systems of physical origin. Increasing the
value of $K$ decreases the numerical value of the decay exponent of the power
law, which is a clear quantitative dynamical effect of trait inheritance.


  The knowledge of cell molecular mechanisms implicated in human diseases is
expanding and should be converted into guidelines for deciphering pathological
cell signaling and suggesting appropriate treatment. The basic assumption is
that during a pathological transformation, the cell does not create new
signaling mechanisms, but rather it hijacks the existing molecular programs.
This affects not only intracellular functions, but also a crosstalk between
different cell types resulting in a new, yet pathological status of the system.
There is a certain combination of molecular characteristics dictating specific
cell signaling states that sustains the pathological disease status.
Identifying and manipulating the key molecular players controlling these cell
signaling states, and shifting the pathological status toward the desired
healthy phenotype, are the major challenge for molecular biology of human
diseases.


  Quartet trees displayed by larger phylogenetic trees have long been used as
inputs for species tree and supertree reconstruction. Computational constraints
prevent the use of all displayed quartets in many practical problems due to the
number of taxa. We introduce the notion of an Efficient Quartet System (EQS) to
represent a phylogenetic tree with a subset of the quartets displayed by the
tree. We show mathematically that the set of quartets obtained from a tree via
an EQS contains all of the combinatorial information of the tree itself. Using
performance tests on simulated datasets, we also demonstrate that using an EQS
to reduce the number of quartets in pipelines for summary methods of species
tree inference and supertree inference results in only small reductions in
accuracy.


  Interaction between proteins is a fundamental mechanism that underlies
virtually all biological processes. Many important interactions are conserved
across a large variety of species. The need to maintain interaction leads to a
high degree of co-evolution between residues in the interface between partner
proteins. The inference of protein-protein interaction networks from the
rapidly growing sequence databases is one of the most formidable tasks in
systems biology today. We propose here a novel approach based on the
Direct-Coupling Analysis of the co-evolution between inter-protein residue
pairs. We use ribosomal and trp operon proteins as test cases: For the small
resp. large ribosomal subunit our approach predicts protein-interaction
partners at a true-positive rate of 70% resp. 90% within the first 10
predictions, with areas of 0.69 resp. 0.81 under the ROC curves for all
predictions. In the trp operon, it assigns the two largest interaction scores
to the only two interactions experimentally known. On the level of residue
interactions we show that for both the small and the large ribosomal subunit
our approach predicts interacting residues in the system with a true positive
rate of 60% and 85% in the first 20 predictions. We use artificial data to show
that the performance of our approach depends crucially on the size of the joint
multiple sequence alignments and analyze how many sequences would be necessary
for a perfect prediction if the sequences were sampled from the same model that
we use for prediction. Given the performance of our approach on the test data
we speculate that it can be used to detect new interactions, especially in the
light of the rapid growth of available sequence data.


  Age dependent population dynamics are frequently modeled with generalizations
of the classic McKendrick-von Foerster equation. These are deterministic
systems, and a stochastic generalization was recently reported in [1,2]. Here
we develop a fully stochastic theory for age-structured populations via quantum
field theoretical Doi-Peliti techniques. This results in a path integral
formulation where birth and death events correspond to cubic and quadratic
interaction terms. This formalism allows us to efficiently recapitulate the
results in [1,2], exemplifying the utility of Doi-Peliti methods. Furthermore,
we find that the path integral formulation for age-structured moments has an
exact perturbative expansion that explicitly relates to the hereditary
structure between correlated individuals. These methods are then generalized
with a binary fission model of cell division.


  We study mathematical models describing the evolution of stochastic
age-structured populations. After reviewing existing approaches, we present a
full kinetic framework for age-structured interacting populations undergoing
birth, death and fission processes, in spatially dependent environments. We
define the complete probability density for the population-size-age-chart and
find results under specific conditions. Connections with more classical models
are also explicitly derived. In particular, we show that factorial moments for
non-interacting processes are described by a natural generalization of the
McKendrick-von Foerster equation, which describes mean-field deterministic
behaviour. Our approach utilizes mixed type, multi-dimensional probability
distributions similar to those employed in the study of gas kinetics, with
terms that satisfy BBGKY-like equation hierarchies.


  Genome-wide experiments to map the DNA-binding locations of
transcription-associated factors (TFs) have shown that the number of genes
bound by a TF far exceeds the number of possible direct target genes.
Distinguishing functional from non-functional binding is therefore a major
challenge in the study of transcriptional regulation. We hypothesized that
functional targets can be discovered by correlating binding and expression
profiles across multiple experimental conditions. To test this hypothesis, we
obtained ChIP-seq and RNA-seq data from matching cell types from the human
ENCODE resource, considered promoter-proximal and distal cumulative regulatory
models to map binding sites to genes, and used a combination of linear and
non-linear measures to correlate binding and expression data. We found that a
high degree of correlation between a gene's TF-binding and expression profiles
was significantly more predictive of the gene being differentially expressed
upon knockdown of that TF, compared to using binding sites in the cell type of
interest only. Remarkably, TF targets predicted from correlation across a
compendium of cell types were also predictive of functional targets in other
cell types. Finally, correlation across a time course of ChIP-seq and RNA-seq
experiments was also predictive of functional TF targets in that tissue.


  Genetic differences between individuals associated to quantitative phenotypic
traits, including disease states, are usually found in non-coding genomic
regions. These genetic variants are often also associated to differences in
expression levels of nearby genes (they are "expression quantitative trait
loci" or eQTLs for short) and presumably play a gene regulatory role, affecting
the status of molecular networks of interacting genes, proteins and
metabolites. Computational systems biology approaches to reconstruct causal
gene networks from large-scale omics data have therefore become essential to
understand the structure of networks controlled by eQTLs together with other
regulatory genes, and to generate detailed hypotheses about the molecular
mechanisms that lead from genotype to phenotype. Here we review the main
analytical methods and softwares to identify eQTLs and their associated genes,
to reconstruct co-expression networks and modules, to reconstruct causal
Bayesian gene and module networks, and to validate predicted networks in
silico.


  The Moran process, as studied by Lieberman, Hauert and Nowak, is a randomised
algorithm modelling the spread of genetic mutations in populations. The
algorithm runs on an underlying graph where individuals correspond to vertices.
Initially, one vertex (chosen u.a.r.) possesses a mutation, with fitness r>1.
All other individuals have fitness 1. During each step of the algorithm, an
individual is chosen with probability proportional to its fitness, and its
state (mutant or non-mutant) is passed on to an out-neighbour which is chosen
u.a.r. If the underlying graph is strongly connected then the algorithm will
eventually reach fixation, in which all individuals are mutants, or extinction,
in which no individuals are mutants. An infinite family of directed graphs is
said to be strongly amplifying if, for every r>1, the extinction probability
tends to 0 as the number of vertices increases. Lieberman et al. proposed two
potentially strongly-amplifying families - superstars and metafunnels.
Heuristic arguments have been published, arguing that there are infinite
families of superstars that are strongly amplifying. The same has been claimed
for metafunnels. In this paper, we give the first rigorous proof that there is
an infinite family of directed graphs that is strongly amplifying. We call the
graphs in the family "megastars". When the algorithm is run on an n-vertex
graph in this family, starting with a uniformly-chosen mutant, the extinction
probability is roughly $n^{-1/2}$ (up to logarithmic factors). We prove that
all infinite families of superstars and metafunnels have larger extinction
probabilities (as a function of n). Finally, we prove that our analysis of
megastars is fairly tight - there is no infinite family of megastars such that
the Moran algorithm gives a smaller extinction probability (up to logarithmic
factors).


  Evolutionary scenarios displaying reticulation events are often represented
by rooted phylogenetic networks. Due to biological reasons, those events occur
very rarely, and, thus, networks containing a minimum number of such events,
so-called minimum hybridization networks, are of particular interest for
research. Moreover, to study reticulate evolution, biologist need not only a
subset but all of those networks. To achieve this goal, the less complex
concept of rooted phylogenetic trees can be used as building block. Here, as a
first important step, the trees are disjoint into common parts, so-called
maximum acyclic agreement forests, which can then be turned into minimum
hybridization networks by applying further network building algorithms. In this
paper, we present two modifications of the first non-naive algorithm --- called
allMAAFs --- computing all maximum acyclic agreement forests for two rooted
binary phylogenetic trees on the same set of taxa. By a simulation study, we
indicate that through these modifications the algorithm is on average 8 times
faster than the original algorithm making this algorithm accessible to larger
input trees and, thus, to a wider range of biological problems.


  There exist several methods dealing with the reconstruction of rooted
phylogenetic networks explaining different evolutionary histories given by
rooted binary phylogenetic trees. In practice, however, due to insufficient
information of the underlying data, phylogenetic trees are in general not
completely resolved and, thus, those methods can often not be applied to
biological data. In this work, we make a first important step to approach this
goal by presenting the first algorithm --- called allMulMAAFs --- that enables
the computation of all relevant nonbinary maximum acyclic agreement forests for
two rooted (nonbinary) phylogenetic trees on the same set of taxa. Notice that
our algorithm is part of the freely available software Hybroscale computing
minimum hybridization networks for a set of rooted (nonbinary) phylogenetic
trees on an overlapping set of taxa.


  With the rise of network science old topics in ecology and economics are
resurfacing. One such topic is structural stability (often referred to as
qualitative stability or sign stability). A system is deemed structurally
stable if the system remains stable for all possible parameter variations so
long as the parameters do not change sign. This type of stability analysis is
appealing when studying real systems as the underlying stability result only
requires the scientist or engineer to know the sign of the parameters in the
model and not the specific values. The necessary and sufficient conditions for
qualitative stability however are opaque. In order to shed light on those
conditions root locus analysis is employed. This technique allows us to
illustrate the necessary conditions for qualitative stability.


  Although RNA-Seq data provide unprecedented isoform-level expression
information, detection of alternative isoform regulation (AIR) remains
difficult, particularly when working with an incomplete transcript annotation.
We introduce JunctionSeq, a new method that builds on the statistical
techniques used by the well-established DEXSeq package to detect differential
usage of both exonic regions and splice junctions. In particular, JunctionSeq
is capable of detecting differentials in novel splice junctions without the
need for an additional isoform assembly step, greatly improving performance
when the available transcript annotation is flawed or incomplete. JunctionSeq
also provides a powerful and streamlined visualization toolset that allows
bioinformaticians to quickly and intuitively interpret their results. We tested
our method on publicly available data from several experiments performed on the
rat pineal gland and Toxoplasma gondii, successfully detecting known and
previously validated AIR genes in 19 out of 19 gene-level hypothesis tests. Due
to its ability to query novel splice sites, JunctionSeq is still able to detect
these differentials even when all alternative isoforms for these genes were not
included in the transcript annotation. JunctionSeq thus provides a powerful
method for detecting alternative isoform regulation even with low-quality
annotations. An implementation of JunctionSeq is available as an R/Bioconductor
package.


  The Asian chestnut gall wasp \emph{Dryocosmus kuriphilus}, native of China,
has become a pest when it appeared in Japan, Korea, and the United States. In
Europe it was first found in Italy, in 2002. In 1982 the host-specific
parasitoid \emph{Torymus sinensis} was introduced in Japan, in an attempt to
achieve a biological control of the pest. After an apparent initial success,
the two species seem to have locked in predator-prey cycles of decadal length.
We have developed a spatially explicit mathematical model that describes the
seasonal time evolution of the adult insect populations, and the competition
for finding egg deposition sites. In a spatially homogeneous situation the
model reduces to an iterated map for the egg density of the two species. While
the map would suggest, for realistic parameters, that both species should
become locally extinct (somewhat corroborating the hypothesis of biological
control), the full model, for the same parameters, shows that the introduction
of \emph{T. sinensis} sparks a traveling wave of the parasitoid population that
destroys the pest on its passage. Depending on the value of the diffusion
coefficients of the two species, the pest can later be able to re-colonize the
empty area left behind the wave. When this occurs the two populations do not
seem to attain a state of spatial homogeneity, but produce an ever-changing
pattern of traveling waves.


  The vertebrate adaptive immune system provides a flexible and diverse set of
molecules to neutralize pathogens. Yet, viruses such as HIV can cause chronic
infections by evolving as quickly as the adaptive immune system, forming an
evolutionary arms race. Here we introduce a mathematical framework to study the
coevolutionary dynamics of antibodies with antigens within a host. We focus on
changes in the binding interactions between the antibody and antigen
populations, which result from the underlying stochastic evolution of genotype
frequencies driven by mutation, selection, and drift. We identify the critical
viral and immune parameters that determine the distribution of antibody-antigen
binding affinities. We also identify definitive signatures of coevolution that
measure the reciprocal response between antibodies and viruses, and we
introduce experimentally measurable quantities that quantify the extent of
adaptation during continual coevolution of the two opposing populations. Using
this analytical framework, we infer rates of viral and immune adaptation based
on time-shifted neutralization assays in two HIV-infected patients. Finally, we
analyze competition between clonal lineages of antibodies and characterize the
fate of a given lineage in terms of the state of the antibody and viral
populations. In particular, we derive the conditions that favor the emergence
of broadly neutralizing antibodies, which may be useful in designing a vaccine
against HIV.


  Much effort has been directed towards using mathematical models to understand
and predict contagious disease, in particular Ebola outbreaks. Classical SIR
(susceptible-infected-recovered) compartmental models capture well the dynamics
of the outbreak in certain communities, and accurately describe the differences
between them based on a variety of parameters. However, repeated resurgence of
Ebola contagions suggests that there are components of the global disease
dynamics that we don't yet fully understand and can't effectively control.
  In order to understand the dynamics of a more widespread contagion, we placed
SIR models within the framework of dynamic networks, with the communities at
risk of contracting the virus acting as nonlinear systems, coupled based on a
connectivity graph. We study how the effects of the disease (measured as the
outbreak impact and duration) change with respect to local parameters, but also
with changes in both short-range and long-range connectivity patterns in the
graph. We discuss the implications of optimizing both these measures in
increasingly realistic models of coupled communities.


  Entropy, under a variety of names, has long been used as a measure of
diversity in ecology, as well as in genetics, economics and other fields. There
is a spectrum of viewpoints on diversity, indexed by a real parameter q giving
greater or lesser importance to rare species. Leinster and Cobbold proposed a
one-parameter family of diversity measures taking into account both this
variation and the varying similarities between species. Because of this latter
feature, diversity is not maximized by the uniform distribution on species. So
it is natural to ask: which distributions maximize diversity, and what is its
maximum value?
  In principle, both answers depend on q, but our main theorem is that neither
does. Thus, there is a single distribution that maximizes diversity from all
viewpoints simultaneously, and any list of species has an unambiguous maximum
diversity value. Furthermore, the maximizing distribution(s) can be computed in
finite time, and any distribution maximizing diversity from some particular
viewpoint q > 0 actually maximizes diversity for all q.
  Although we phrase our results in ecological terms, they apply very widely,
with applications in graph theory and metric geometry.


  A central task in analyzing complex dynamics is to determine the loci of
information storage and the communication topology of information flows within
a system. Over the last decade and a half, diagnostics for the latter have come
to be dominated by the transfer entropy. Via straightforward examples, we show
that it and a derivative quantity, the causation entropy, do not, in fact,
quantify the flow of information. At one and the same time they can
overestimate flow or underestimate influence. We isolate why this is the case
and propose several avenues to alternate measures for information flow. We also
address an auxiliary consequence: The proliferation of networks as a now-common
theoretical model for large-scale systems, in concert with the use of
transfer-like entropies, has shoehorned dyadic relationships into our
structural interpretation of the organization and behavior of complex systems.
This interpretation thus fails to include the effects of polyadic dependencies.
The net result is that much of the sophisticated organization of complex
systems may go undetected.


  We consider the Moran process with two populations competing under an
iterated Prisoners' Dilemma in the presence of mutation, and concentrate on the
case where there are multiple Evolutionarily Stable Strategies. We perform a
complete bifurcation analysis of the deterministic system which arises in the
infinite population size. We also study the Master equation and obtain
asymptotics for the invariant distribution and metastable switching times for
the stochastic process in the case of large but finite population. We also show
that the stochastic system has asymmetries in the form of a skew for parameter
values where the deterministic limit is symmetric.


  This methods paper presents computational protocols for the identification of
non-coding RNA genes or RNA motifs within genomic sequences. An application to
bacterial small RNA is proposed.


  Several results about sign properties of Metzler matrices are obtained. It is
first established that checking the sign-stability of a Metzler sign-matrix can
be either characterized in terms of the Hurwitz stability of the unit
sign-matrix in the corresponding qualitative class, or in terms the negativity
of the diagonal elements of the Metzler sign-matrix and the acyclicity of the
associated directed graph. Similar results are obtained for the case of Metzler
block-matrices and Metzler mixed-matrices, the latter being a class of Metzler
matrices containing both sign- and real-type entries. The problem of assessing
the sign-stability of the convex hull of a finite and summable family of
Metzler matrices is also solved, and a necessary and sufficient condition for
the existence of common Lyapunov functions for all the matrices in the convex
hull is obtained. The concept of sign-stability is then generalized to the
concept of Ker$_+(B)$-sign-stability, a problem that arises in the analysis of
certain jump Markov processes. A sufficient condition for the
Ker$_+(B)$-sign-stability of Metzler sign-matrices is obtained and formulated
using inverses of sign-matrices and the concept of $L^+$-matrices. Several
applications of the results are discussed in the last section.


  Phylogenetic varieties related to equivariant substitution models have been
studied largely in the last years. One of the main objectives has been finding
a set of generators of the ideal of these varieties, but this has not yet been
achieved in some cases (for example, for the general Markov model this involves
the open "salmon conjecture") and it is not clear how to use all generators in
practice. Motivated by applications in biology, we tackle the problem from
another point of view. The elements of the ideal that could be useful for
applications in phylogenetics only need to describe the variety around certain
"points of no evolution". We produce a collection of explicit equations that
describe the variety on a Zariski open neighborhood of these points. Namely,
for any tree on any number of leaves (and any degrees at the interior nodes)
and for any equivariant model on any set of states, we compute the codimension
of the corresponding phylogenetic variety. We prove that this variety is smooth
at general points of no evolution, and provide an algorithm to produce a
complete intersection that describes the variety around these points.


  For group-living animals, reaching consensus to stay cohesive is crucial for
their fitness, particularly when collective motion starts and stops.
Understanding the decision-making at individual and collective levels upon
sudden disturbances is central in the study of collective animal behavior, and
concerns the broader question of how information is distributed and evaluated
in groups. Despite the relevance of the problem, well-controlled experimental
studies that quantify the collective response of groups facing disruptive
events are lacking. Here we study the behavior of groups of uninformed
individuals subject to the departure and stop of a trained conspecific within
small-sized groups. We find that the groups reach an effective consensus:
either all uninformed individuals follow the trained one (and collective motion
occurs) or none does it. Combining experiments and a simple mathematical model
we show that the observed phenomena results from the interplay between simple
mimetic rules and the characteristic duration of the stimulus, here, the time
the trained individual is moving away. The proposed mechanism strongly depends
on group size, as observed in the experiments, and though group splitting can
occur, the most likely outcome is always a coherent collective group response
(consensus). The prevalence of a consensus is expected even if the groups of
naives face conflicting information, e.g. if groups contain two subgroups of
trained individuals, one trained to stay and one trained to leave. Our results
indicate that collective decision-making and consensus in (small) animal groups
are likely to be self-organized phenomena that do not involve concertation or
even communication among the group members.


  Can we quantify the change of complexity throughout evolutionary processes?
We attempt to address this question through an empirical approach. In very
general terms, we simulate two simple organisms on a computer that compete over
limited available resources. We implement Global Rules that determine the
interaction between two Elementary Cellular Automata on the same grid. Global
Rules change the complexity of the state evolution output which suggests that
some complexity is intrinsic to the interaction rules themselves. The largest
increases in complexity occurred when the interacting elementary rules had very
little complexity, suggesting that they are able to accept complexity through
interaction only. We also found that some Class 3 or 4 CA rules are more
fragile than others to Global Rules, while others are more robust, hence
suggesting some intrinsic properties of the rules independent of the Global
Rule choice. We provide statistical mappings of Elementary Cellular Automata
exposed to Global Rules and different initial conditions onto different
complexity classes.


  In phylogenetics, distances are often used to measure the incongruence
between a pair of phylogenetic trees that are reconstructed by different
methods or using different regions of genome. Motivated by the maximum
parsimony principle in tree inference, we recently introduced the maximum
parsimony (MP) distance, which enjoys various attractive properties due to its
connection with several other well-known tree distances, such as TBR and SPR.
Here we show that computing the MP distance between two trees, a NP-hard
problem in general, is fixed parameter tractable in terms of the TBR distance
between the tree pair. Our approach is based on two reduction rules--the chain
reduction and the subtree reduction--that are widely used in computing TBR and
SPR distances. More precisely, we show that reducing chains to length 4 (but
not shorter) preserves the MP distance. In addition, we describe a
generalization of the subtree reduction which allows the pendant subtrees to be
rooted in different places, and show that this still preserves the MP distance.
On a slightly different note we also show that Monadic Second Order Logic
(MSOL), posited over an auxiliary graph structure known as the display graph
(obtained by merging the two trees at their leaves), can be used to obtain an
alternative proof that computation of MP distance is fixed parameter tractable
in terms of TBR-distance. We conclude with an extended discussion in which we
focus on similarities and differences between MP distance and TBR distance and
present a number of open problems. One particularly intriguing question,
emerging from the MSOL formulation, is whether two trees with bounded MP
distance induce display graphs of bounded treewidth.


  A commonly used approach to study stability in a complex system is by
analyzing the Jacobian matrix at an equilibrium point of a dynamical system.
The equilibrium point is stable if all eigenvalues have negative real parts.
Here, by obtaining eigenvalue bounds of the Jacobian, we show that stable
complex systems will favor mutualistic and competitive relationships that are
asymmetrical (non-reciprocative) and trophic relationships that are symmetrical
(reciprocative). Additionally, we define a measure called the interdependence
diversity that quantifies how distributed the dependencies are between the
dynamical variables in the system. We find that increasing interdependence
diversity has a destabilizing effect on the equilibrium point, and the effect
is greater for trophic relationships than for mutualistic and competitive
relationships. These predictions are consistent with empirical observations in
ecology. More importantly, our findings suggest stabilization algorithms that
can apply very generally to a variety of complex systems.