  Reconstruction of gene regulatory networks or 'reverse-engineering' is a
process of identifying gene interaction networks from experimental microarray
gene expression profile through computation techniques. In this paper, we tried
to reconstruct cancer-specific gene regulatory network using information
theoretic approach - mutual information. The considered microarray data
consists of large number of genes with 20 samples - 12 samples from colon
cancer patient and 8 from normal cell. The data has been preprocessed and
normalized. A t-test statistics has been applied to filter differentially
expressed genes. The interaction between filtered genes has been computed using
mutual information and ten different networks has been constructed with varying
number of interactions ranging from 30 to 500. We performed the topological
analysis of the reconstructed network, revealing a large number of interactions
in colon cancer. Finally, validation of the inferred results has been done with
available biological databases and literature.


  We propose a novel combination of methods that (i) portrays quantitative
characteristics of a DNA sequence as an image, (ii) computes distances between
these images, and (iii) uses these distances to output a map wherein each
sequence is a point in a common Euclidean space. In the resulting "Molecular
Distance Map" each point signifies a DNA sequence, and the geometric distance
between any two points reflects the degree of relatedness between the
corresponding sequences and species.
  Molecular Distance Maps present compelling visual representations of
relationships between species and could be used for taxonomic clarifications,
for species identification, and for studies of evolutionary history. One of the
advantages of this method is its general applicability since, as sequence
alignment is not required, the DNA sequences chosen for comparison can be
completely different regions in different genomes. In fact, this method can be
used to compare any two DNA sequences. For example, in our dataset of 3,176
mitochondrial DNA sequences, it correctly finds the mtDNA sequences most
closely related to that of the anatomically modern human (the Neanderthal, the
Denisovan, and the chimp), and it finds that the sequence most different from
it belongs to a cucumber. Furthermore, our method can be used to compare real
sequences to artificial, computer-generated, DNA sequences. For example, it is
used to determine that the distances between a Homo sapiens sapiens mtDNA and
artificial sequences of the same length and same trinucleotide frequencies can
be larger than the distance between the same human mtDNA and the mtDNA of a
fruit-fly.
  We demonstrate this method's promising potential for taxonomical
clarifications by applying it to a diverse variety of cases that have been
historically controversial, such as the genus Polypterus, the family Tarsiidae,
and the vast (super)kingdom Protista.


  Protein complexes conserved across species indicate processes that are core
to cellular machinery (e.g. cell-cycle or DNA damage-repair complexes conserved
across human and yeast). While numerous computational methods have been devised
to identify complexes from the protein interaction (PPI) networks of individual
species, these are severely limited by noise and errors (false positives) in
currently available datasets. Our analysis using human and yeast PPI networks
revealed that these methods missed several important complexes including those
conserved between the two species (e.g. the MLH1-MSH2-PMS2-PCNA mismatch-repair
complex). Here, we note that much of the functionalities of yeast complexes
have been conserved in human complexes not only through sequence conservation
of proteins but also of critical functional domains. Therefore, integrating
information of domain conservation might throw further light on conservation
patterns between yeast and human complexes.


  A mechanism of the translation in oligomer world is proposed. The translation
is carried out by a minimum cycle, which is sustained by adaptors and
helicases, and the first information processing in oligomer world. We expect
that such a cycle actually worked in a primitive cell and can be constructed in
vitro. By computer simulation we have shown that a proofreading is achieved by
the fluctuation in the cell. It is rather paradoxical that the proofreading is
effective for the system consisting of molecular machines with low efficiency.


  To function as gene regulatory elements in response to environmental signals,
riboswitches must adopt specific secondary structures on appropriate time
scales. We employ kinetic Monte Carlo simulation to model the time-dependent
folding during transcription of TPP riboswitch expression platforms. According
to our simulations, riboswitch transcriptional terminators, which must adopt a
specific hairpin configuration by the time they have been transcribed, fold
with higher efficiency than Shine-Dalgarno sequesterers, whose proper structure
is required only at the time of ribosomal binding. Our findings suggest both
that riboswitch transcriptional terminator sequences have been naturally
selected for high folding efficiency, and that sequesterers can maintain their
function even in the presence of significant misfolding.


  Infinite Leslie matrices, introduced by Demetrius forty years ago are
mathematical models of age-structured populations defined by a countable
infinite number of age classes. This article is concerned with determining
solutions of the discrete dynamical system in finite time. We address this
problem by appealing to the concept of kneading matrices and kneading
determinants. Our analysis is applicable not only to populations models, but to
models of self-reproducing machines and self-reproducing computer programs. The
dynamics of theses systems can also be described in terms of infinite Leslie
matrices.uation.


  In this work we examine how volume exclusion caused by regions of high
chromatin density might influence the time required for proteins to find
specific DNA binding sites. The spatial variation of chromatin density within
mouse olfactory sensory neurons is determined from soft X-ray tomography
reconstructions of five nuclei. We show that there is a division of the nuclear
space into regions of low-density euchromatin and high-density heterochromatin.
Volume exclusion experienced by a diffusing protein caused by this varying
density of chromatin is modeled by a repulsive potential. The value of the
potential at a given point in space is chosen to be proportional to the density
of chromatin at that location. The constant of proportionality, called the
volume exclusivity, provides a model parameter that determines the strength of
volume exclusion. Numerical simulations demonstrate that the mean time for a
protein to locate a binding site localized in euchromatin is minimized for a
finite, non-zero volume exclusivity. For binding sites in heterochromatin, the
mean time is minimized when the volume exclusivity is zero (the protein
experiences no volume exclusion). An analytical theory is developed to explain
these results. The theory suggests that for binding sites in euchromatin there
is an optimal level of volume exclusivity that balances a reduction in the
volume searched in finding the binding site, with the height of effective
potential barriers the protein must cross during the search process.


  A steady influx of a single deleterious multilocus genotype will impose
genetic load on the resident population and leave multiple descendants carrying
various numbers of the foreign alleles. Provided that the foreign types are
rare at equilibrium, and that all immigrant genes will eventually be eliminated
by selection, the population structure can be inferred explicitly from the
deterministic branching process taking place within a single immigrant lineage.
Unless the migration and recombination rates were high, this simple method was
a very close approximation to the simulated migration-selection balance with
all possible multilocus genotypes considered.


  The shift from outcrossing to self-fertilization is among the most common
transitions in plants. Until recently, however, a genome-wide view of this
transition has been obscured by a dearth of appropriate data and the lack of
appropriate population genomic methods to interpret such data. Here, we present
novel analyses detailing the origin of the selfing species, Capsella rubella,
which recently split from its outcrossing sister, Capsella grandiflora. Due to
the recency of the split, most variation within C. rubella is found within C.
grandiflora. We can therefore identify genomic regions where two C. rubella
individuals have inherited the same or different segments of ancestral
diversity (i.e. founding haplotypes) present in C. rubella's founder(s). Based
on this analysis, we show that C. rubella was founded by multiple individuals
drawn from a diverse ancestral population closely related to extant C.
grandiflora, that drift and selection have rapidly homogenized most of this
ancestral variation since C. rubella's founding, and that little novel
variation has accumulated within this time. Despite the extensive loss of
ancestral variation, the approximately 25% of the genome for which two C.
rubella individuals have inherited different founding haplotypes makes up
roughly 90% of the genetic variation between them. To extend these findings, we
develop a coalescent model that utilizes the inferred frequency of founding
haplotypes and variation within founding haplotypes to estimate that C. rubella
was founded by a potentially large number of individuals 50-100 kya, and has
subsequently experienced a 20X reduction in its effective population size. As
population genomic data from an increasing number of outcrossing/selfing pairs
are generated, analyses like this here will facilitate a fine-scaled view of
the evolutionary and demographic impact of the transition to
self-fertilization.


  While much effort has focused on detecting positive and negative directional
selection in the human genome, relatively little work has been devoted to
balancing selection. This lack of attention is likely due to the paucity of
sophisticated methods for identifying sites under balancing selection. Here we
develop two composite likelihood ratio tests for detecting balancing selection.
Using simulations, we show that these methods outperform competing methods
under a variety of assumptions and demographic models. We apply the new methods
to whole-genome human data, and find a number of previously-identified loci
with strong evidence of balancing selection, including several HLA genes.
Additionally, we find evidence for many novel candidates, the strongest of
which is FANK1, an imprinted gene that suppresses apoptosis, is expressed
during meiosis in males, and displays marginal signs of segregation distortion.
We hypothesize that balancing selection acts on this locus to stabilize the
segregation distortion and negative fitness effects of the distorter allele.
Thus, our methods are able to reproduce many previously-hypothesized signals of
balancing selection, as well as discover novel interesting candidates.


  In social situations with which evolutionary game is concerned, individuals
are considered to be heterogeneous in various aspects. In particular, they may
differently perceive the same outcome of the game owing to heterogeneity in
idiosyncratic preferences, fighting abilities, and positions in a social
network. In such a population, an individual may imitate successful and similar
others, where similarity refers to that in the idiosyncratic fitness function.
I propose an evolutionary game model with two subpopulations on the basis of
multipopulation replicator dynamics to describe such a situation. In the
proposed model, pairs of players are involved in a two-person game as a
well-mixed population, and imitation occurs within subpopulations in each of
which players have the same payoff matrix. It is shown that the model does not
allow any internal equilibrium such that the dynamics differs from that of
other related models such as the bimatrix game. In particular, even a slight
difference in the payoff matrix in the two subpopulations can make the opposite
strategies to be stably selected in the two subpopulations in the snowdrift and
coordination games.


  Neither of the two prevalent theories, namely thermodynamic stability and
kinetic stability, provides a comprehensive understanding of protein folding.
The thermodynamic theory is misleading because it assumes that free energy is
the exclusive dominant mechanism of protein folding, and attributes the
structural transition from one characteristic state to another to energy
barriers. Conversely, the concept of kinetic stability overemphasizes dominant
mechanisms that are related to kinetic factors. This article explores the
stability condition of protein structures from the viewpoint of meso-science,
paying attention to the compromise in the competition between minimum free
energy and other dominant mechanisms. Based on our study of complex systems, we
propose that protein folding is a meso-scale, dissipative, nonlinear and
non-equilibrium process that is dominated by the compromise between free energy
and other dominant mechanisms such as environmental factors. Consequently, a
protein shows dynamic structures, featuring characteristic states that appear
alternately and dynamically, only one of which is the state with minimum free
energy. To provide evidence for this concept, we analyzed the time series of
energetic and structural changes of three simulations of protein
folding/unfolding. Our results indicate that thorough consideration of the
multiple dynamic characteristic structures generated by multiple mechanisms may
be the key to understanding protein folding.


  In a recent paper, we applied time series analysis methods to study the
possible influence of the solar motion through the Galaxy on terrestrial
extinction (Feng & Bailer-Jones 2013). We drew conclusions about the relative
probabilities of how well different models explain the geological extinction
record. We found no strong evidence for either a periodicity in the extinction
record or for a link between the solar orbit and this extinction record. In
response to this, Melott & Bambach (2013) dedicate the majority of their
article to criticizing our analysis. Their main objection is to evidence-based
model comparison. Here we deal with their criticisms one by one, and show that
they are largely misplaced, leaving our conclusions unaffected. The bottom line
is that although one can get a periodic model at some period to fit the
extinction data, there are other models which explain the data better.


  Addressing the functionality of genomes is one of the most important and
challenging tasks of today's biology. In particular the ability to link
genotypes to corresponding phenotypes is of interest in the reconstruction and
biotechnological manipulation of metabolic pathways. Over the last years, the
OmniLogTM Phenotype Microarray (PM) technology has been used to address many
specific issues related to the metabolic functionality of microorganisms.
However, computational tools that could directly link PM data with the gene(s)
of interest followed by the extraction of information on genephenotype
correlation are still missing. Here we present DuctApe, a suite that allows the
analysis of both genomic sequences and PM data, to find metabolic differences
among PM experiments and to correlate them with KEGG pathways and gene
presence/absence patterns. As example, an application of the program to four
bacterial datasets is presented. The source code and tutorials are available at
http://combogenomics.github.io/DuctApe/.


  Genes are not located randomly along genomes. Synteny, the conservation of
their relative positions in genomes of different species, reflects fundamental
constraints on natural evolution. We present approaches to infer pairs of
co-localized genes from multiple genomes, describe their organization, and
study their evolutionary history. In bacterial genomes, we thus identify
synteny units, or "syntons", which are clusters of proximal genes that
encompass and extend operons. The size distribution of these syntons divide
them into large syntons, which correspond to fundamental macro-molecular
complexes of bacteria, and smaller ones, which display a remarkable exponential
distribution of sizes. This distribution is "universal" in two respects: it
holds for vastly different genomes, and for functionally distinct genes.
Similar statistical laws have been reported previously in studies of bacterial
genomes, and generally attributed to purifying selection or neutral processes.
Here, we perform a new analysis based on the concept of parsimony, and find
that the prevailing evolutionary mechanism behind the formation of small
syntons is a selective process of gene aggregation. Altogether, our results
imply a common evolutionary process that selectively shapes the organization
and diversity of bacterial genomes.


  Computational modeling of Wnt signaling pathway has gained prominence for its
use as computer aided diagnostic tool to develop therapeutic cancer target
drugs and predict of test samples as cancerous and non cancerous. This
manuscript focuses on development of simple static bayesian network models of
varying complexity that encompasses prior partially available biological
knowledge about intra and extra cellular factors affecting the Wnt pathway and
incorporates epigenetic information like methylation and histone modification
of a few genes known to have inhibitory affect on Wnt pathway. It might be
expected that such models not only increase cancer prediction accuracies and
also form basis for understanding Wnt signaling activity in different states of
tumorigenesis. Initial results in human colorectal cancer cases indicate that
incorporation of epigenetic information increases prediction accuracy of test
samples as being tumorous or normal. Receiver Operator Curves (ROC) and their
respective area under the curve (AUC) measurements, obtained from predictions
of state of test sample and corresponding predictions of the state of
activation of transcription complex of the Wnt pathway for the test sample,
indicate that there is significant difference between the Wnt pathway being on
(off) and its association with the sample being tumorous (normal). Two sample
Kolmogorov-Smirnov test confirm the statistical deviation between the
distributions of these predictions. At a preliminary stage, use of these models
may help in understanding the yet unknown effect of certain factors like DKK2,
DKK3-1 and SFRP-2/3/5 on {\beta}-catenin transcription complex.


  Empirical evidence suggesting that living systems might operate in the
vicinity of critical points, at the borderline between order and disorder, has
proliferated in recent years, with examples ranging from spontaneous brain
activity to flock dynamics. However, a well-founded theory for understanding
how and why interacting living systems could dynamically tune themselves to be
poised in the vicinity of a critical point is lacking. Here we employ tools
from statistical mechanics and information theory to show that complex adaptive
or evolutionary systems can be much more efficient in coping with diverse
heterogeneous environmental conditions when operating at criticality.
Analytical as well as computational evolutionary and adaptive models vividly
illustrate that a community of such systems dynamically self-tunes close to a
critical state as the complexity of the environment increases while they remain
non-critical for simple and predictable environments. A more robust convergence
to criticality emerges in co-evolutionary and co-adaptive set-ups in which
individuals aim to represent other agents in the community with fidelity,
thereby creating a collective critical ensemble and providing the best possible
trade-off between accuracy and flexibility. Our approach provides a
parsimonious and general mechanism for the emergence of critical-like behavior
in living systems needing to cope with complex environments or trying to
efficiently coordinate themselves as an ensemble.


  In contrast to the neutral population cycles of the deterministic mean-field
Lotka--Volterra rate equations, including spatial structure and stochastic
noise in models for predator-prey interactions yields complex spatio-temporal
structures associated with long-lived erratic population oscillations.
Environmental variability in the form of quenched spatial randomness in the
predation rates results in more localized activity patches. Population
fluctuations in rare favorable regions in turn cause a remarkable increase in
the asymptotic densities of both predators and prey. Very intriguing features
are found when variable interaction rates are affixed to individual particles
rather than lattice sites. Stochastic dynamics with demographic variability in
conjunction with inheritable predation efficiencies generate non-trivial time
evolution for the predation rate distributions, yet with overall essentially
neutral optimization.


  A model of telomere dynamics is proposed and examined. Our model, which
extends a previously introduced two-compartment model that incorporates stem
cells as progenitors of new cells, imposes the Hayflick Limit, the maximum
number of cell divisions that are possible. This new model leads to cell
populations for which the average telomere length is not necessarily a
monotonically decreasing function of time, in contrast to previously published
models. We provide a phase diagram indicating where such results would be
expected. In addition, qualitatively different results are obtained for the
evolution of the total cell population. Last, in comparison to available
leukocyte baboon data, this new model is shown to provide a better fit to
biological data.


  We study the combined influence of selection and random fluctuations on the
evolutionary dynamics of two-strategy ("cooperation" and "defection") games in
populations comprising cooperation facilitators. The latter are individuals
that support cooperation by enhancing the reproductive potential of cooperators
relative to the fitness of defectors. By computing the fixation probability of
a single cooperator in finite and well-mixed populations that include a fixed
number of facilitators, and by using mean field analysis, we determine when
selection promotes cooperation in the important classes of prisoner's dilemma,
snowdrift and stag-hunt games. In particular, we identify the circumstances
under which selection favors the replacement and invasion of defection by
cooperation. Our findings, corroborated by stochastic simulations, show that
the spread of cooperation can be promoted through various scenarios when the
density of facilitators exceeds a critical value whose dependence on the
population size and selection strength is analyzed. We also determine under
which conditions cooperation is more likely to replace defection than vice
versa.
  Keywords: Evolutionary games; dynamics of cooperation; social dilemmas;
fixation; population dynamics.


  A reply to Drake (2013) "Early warning signals of stochastic switching"
http://dx.doi.org/10.1098/rspb.2013.0686


  Under random mating, a progeny's alleles are independently sampled from the
parental gene pools. Here is a new proof which avoids the usual algebraic
complexity, based on a restated Mendel's First Law. Another simplified proof
along the old approach led to the discovery that allelic independence can hold
under random mating and fertility selection. The theoretical existence and
number of such selection coefficients are established.


  Clathrin mediated endocytosis (CME) is an ubiquitous cellular pathway that
regulates central aspects of cell physiology such as nutrient uptake,
modulation of signal transduction, synaptic transmission and membrane
turn-over. Endocytic vesicle formation depends on the timed production of
specific phosphoinositides and their interactions with various endocytic
proteins. Recently, it has been found that
phosphatidylinositol-3,4-bisphosphate (PI(3,4)P2) produced by the class II
phosphatidylinositol 3-kinase C2alpha plays a key role in the recruitment of
the PX-BAR domain protein SNX9, which is proposed to play a role in the
constriction of the endocytic vesicle neck [Posor et al, Nature 499, p233
(2013)]. Interestingly, SNX9 and its close paralog SNX18 are not fully specific
to PI(3,4)P2 but can also bind other phospholipids, in particular to PI(4,5)P2,
an abundant plasma membrane lipid required for the recruitment of many
endocytic proteins. In order to understand the dynamical interplay between
phospholipids and endocytic proteins, we developed a computational model of the
temporal changes in the population of the phosphoinositide-associated endocytic
proteins and their spatial distribution at a clathrin-coated pit (CCP). The
model resolves single molecules in time and space, and incorporates the complex
interplay of proteins and lipids, as well as their movement within the CCP. We
find that the comparably small differences in lipid binding affinities of
endocytic proteins are amplified by competition among them, allowing for the
selective enrichment of SNX9 at late stage CCPs as a result of timed PI(3,4)P2
production.


  Despite the enormous relevance of zoonotic infections to world- wide public
health, and despite much effort in modeling individual zoonoses, a fundamental
understanding of the disease dynamics and the nature of outbreaks arising in
such systems is still lacking. We introduce a simple stochastic model of
susceptible-infected- recovered dynamics in a coupled animal-human
metapopulation, and solve analytically for several important properties of the
cou- pled outbreaks. At early timescales, we solve for the probability and time
of spillover, and the disease prevalence in the animal population at spillover
as a function of model parameters. At long times, we characterize the
distribution of outbreak sizes and the critical threshold for a large human
outbreak, both of which show a strong dependence on the basic reproduction
number in the animal population. The coupling of animal and human infection
dynamics has several crucial implications, most importantly al- lowing for the
possibility of large human outbreaks even when human-to-human transmission is
subcritical.


  Most available antimicrobial peptides (AMP) prediction methods use common
approach for different classes of AMP. Contrary to available approaches, we
suggest, that a strategy of prediction should be based on the fact, that there
are several kinds of AMP which are vary in mechanisms of action, structure,
mode of interaction with membrane etc. According to our suggestion for each
kind of AMP a particular approach has to be developed in order to get high
efficacy. Consequently in this paper a particular but the biggest class of AMP
- linear cationic antimicrobial peptides (LCAP) - has been considered and a
newly developed simple method of LCAP prediction described. The aim of this
study is the development of a simple method of discrimination of AMP from
non-AMP, the efficiency of which will be determined by efficiencies of selected
descriptors only and comparison the results of the discrimination procedure
with the results obtained by more complicated discriminative methods. As
descriptors the physicochemical characteristics responsible for capability of
the peptide to interact with an anionic membrane were considered. The following
characteristics such as hydrophobicity, amphiphaticity, location of the peptide
in relation to membrane, charge, propensity to disordered structure were
studied. On the basis of these characteristics a new simple algorithm of
prediction is developed and evaluation of efficacies of the characteristics as
descriptors performed. The results show that three descriptors: hydrophobic
moment, charge and location of the peptide along the membranes can be used as
discriminators of LCAPs. For the training set our method gives the same level
of accuracy as more complicated machine learning approaches offered as CAMP
database service tools. For the test set sensitivity obtained by our method
gives even higher value than the one obtained by CAMP prediction tools.


  Clinical adoption of human genome sequencing requires methods with known
accuracy of genotype calls at millions or billions of positions across a
genome. Previous work showing discordance amongst sequencing methods and
algorithms has made clear the need for a highly accurate set of genotypes
across a whole genome that could be used as a benchmark. We present methods to
make highly confident SNP, indel, and homozygous reference genotype calls for
NA12878, the pilot genome for the Genome in a Bottle Consortium. We minimize
bias towards any method by integrating and arbitrating between 14 datasets from
5 sequencing technologies, 7 mappers, and 3 variant callers. Regions for which
no confident genotype call could be made are identified as uncertain, and
classified into different reasons for uncertainty. Our highly confident
genotype calls are publicly available on the Genome Comparison and Analytic
Testing (GCAT) website to enable real-time benchmarking of any method.


  Understanding the causes and effects of network structural features is a key
task in deciphering complex systems. In this context, the property of network
nestedness has aroused a fair amount of interest as regards ecological
networks. Indeed, Bastolla et al. introduced a simple measure of network
nestedness which opened the door to analytical understanding, allowing them to
conclude that biodiversity is strongly enhanced in highly nested mutualistic
networks. Here, we suggest a slightly refined version of such a measure and go
on to study how it is influenced by the most basic structural properties of
networks, such as degree distribution and degree-degree correlations (i.e.
assortativity). We find that heterogeneity in the degree has a very strong
influence on nestedness. Once such an influence has been discounted, we find
that nestedness is strongly correlated with disassortativity and hence, as
random (neutral) networks have been recently found to be naturally
disassortative, they tend to be naturally nested just as the result of chance.


  Epithelial to mesenchymal transition (EMT) plays important roles in embryonic
development, tissue regeneration and cancer metastasis. While several feedback
loops have been shown to regulate EMT, it remains elusive how they coordinately
modulate EMT response to TGF-\beta treatment. We construct a mathematical model
for the core regulatory network controlling TGF-\beta-induced EMT. Through
deterministic analyses and stochastic simulations, we show that EMT is a
sequential two-step program that an epithelial cell first transits to partial
EMT then to the mesenchymal state, depending on the strength and duration of
TGF-\beta stimulation. Mechanistically the system is governed by coupled
reversible and irreversible bistable switches. The SNAIL1/miR-34 double
negative feedback loop is responsible for the reversible switch and regulates
the initiation of EMT, while the ZEB/miR-200 feedback loop is accountable for
the irreversible switch and controls the establishment of the mesenchymal
state. Furthermore, an autocrine TGF-\beta/miR-200 feedback loop makes the
second switch irreversible, modulating the maintenance of EMT. Such coupled
bistable switches are robust to parameter variation and molecular noise. We
provide a mechanistic explanation on multiple experimental observations. The
model makes several explicit predictions on hysteretic dynamic behaviors,
system response to pulsed stimulation and various perturbations, which can be
straightforwardly tested.


  Evolutionary graph theory has grown to be an area of intense study. Despite
the amount of interest in the field, it seems to have grown separate from other
subfields of population genetics and evolution. In the current work I introduce
the concept of Fisher's (1930) reproductive value into the study of evolution
on graphs. Reproductive value is a measure of the expected genetic contribution
of an individual to a distant future generation. In a heterogeneous
graph-structured population, differences in the number of connections among
individuals translates into differences in the expected number of offspring,
even if all individuals have the same fecundity. These differences are
accounted for by reproductive value. The introduction of reproductive value
permits the calculation of the fixation probability of a mutant in a neutral
evolutionary process in any graph-structured population for either the moran
birth-death or death-birth process.


  Many kinase inhibitors have been approved as cancer therapies. Recently,
libraries of kinase inhibitors have been extensively profiled, thus providing a
map of the strength of action of each compound on a large number of its
targets. These profiled libraries define drug-kinase networks that can predict
the effectiveness of new untested drugs and elucidate the role played by
specific kinases in different cellular systems. Predictions of drug
effectiveness based on a comprehensive network model of cellular signalling are
difficult, due to our partial knowledge of the complex biological processes
downstream of the targeted kinases. We have developed the Kinase Inhibitors
Elastic Net (KIEN) method, which integrates information contained in
drug-kinase networks with in vitro screening. The method uses the in vitro cell
response of single drugs and drug pair combinations as a training set to build
linear and nonlinear regression models. Besides predicting the effectiveness of
untested drugs, the method identifies sets of kinases that are statistically
associated to drug sensitivity in a given cell line. We compare different
versions of the method, which is based on a regression technique known as
elastic net. Data from two-drug combinations leads to predictive models, and
predictivity can be improved by applying logarithmic transformation to the
data. The method is applied to the A549 lung cancer cell line. A pathway
enrichment analysis of the set of kinases identified by the method shows that
axon guidance, activation of Rac, and semaphorin interactions pathways are
associated to a selective response to therapeutic intervention in this cell
line.


  In our previous study, we showed that the branching process approximation is
useful for estimating metabolic robustness, measured using the impact degree.
By applying a theory of random family forests, we here extend the branching
process approximation to consider the knockout of {\it multiple} reactions,
inspired by the importance of multiple knockouts reported by recent
computational and experimental studies. In addition, we propose a better
definition of the number of offspring of each reaction node, allowing for an
improved estimation of the impact degree distribution obtained as a result of a
single knockout. Importantly, our proposed approach is also applicable to
multiple knockouts. The comparisons between theoretical predictions and
numerical results using real-world metabolic networks demonstrate the validity
of the modeling based on random family forests for estimating the impact degree
distributions resulting from the knockout of multiple reactions.


  Standing genetic variation provides a rich reservoir of potentially useful
mutations facilitating the adaptation to novel environments. Experimental
evolution studies have demonstrated that rapid and strong phenotypic responses
to selection can also be obtained in the laboratory. When combined with the
Next Generation Sequencing technology, these experiments promise to identify
the individual loci contributing to adaption. Nevertheless, until now, very
little is known about the design of such evolve and resequencing (E&R) studies.
Here, we use forward simulations of entire genomes to evaluate different
experimental designs that aim to maximize the power to detect selected
variants. We show that low linkage disequilibrium in the starting population,
population size, duration of the experiment and the number of replicates are
the key factors in determining the power and accuracy of E&R studies.
Furthermore, replication of E&R is more important for detecting the targets of
selection than increasing the population size. Using an optimized design
beneficial loci with a selective advantage as low as s=0.005 can be identified
at the nucleotide level. Even when a large number of loci are selected
simultaneously, up to 56% can be reliably detected without incurring large
numbers of false positives. Our computer simulations suggest that, with an
adequate experimental design, E&R studies are a powerful tool to identify
adaptive mutations from standing genetic variation and thereby provide an
excellent means to analyze the trajectories of selected alleles in evolving
populations.


  We prove that the seasonally-forced SIR model with a T-periodic forcing has a
periodic solution with period T whenever the basic reproductive number R0>1.
The proof uses the Leray-Schauder degree theory. We also describe some
numerical results in which we compute the T-periodic solution, where in order
to obtain the T-periodic solution when the behavior of the system is
subharmonic or chaotic, we use a Galerkin scheme.


  The reconstruction of a central tendency `species tree' from a large number
of conflicting gene trees is a central problem in systematic biology. Moreover,
it becomes particularly problematic when taxon coverage is patchy, so that not
all taxa are present in every gene tree. Here, we list four desirable
properties that a method for estimating a species tree from gene trees should
have. We show that while these can be achieved when taxon coverage is complete
(by the Adams consensus method), they cannot all be satisfied in the more
general setting of partial taxon coverage.


  We use the worm-like chain model to study supercoiling of DNA under tension
and torque. The model reproduces experimental data for a broad range of forces,
salt concentrations and contour lengths. We find a plane of first order phase
transitions ending in a smeared out line of critical points, the
multi-plectoneme phase, which is characterized by a fast twist mediated
diffusion of plectonemes and a torque that rises after plectoneme formation
with increasing linking number. The discovery of this new phase at the same
time resolves the discrepancies between existing models and experiment.


  A method for detecting a trend change in cross-border epidemic transmission
is developed for a standard epidemiological SIR compartment model and a
meta-population network model. The method is applicable to investigating the
efficacy of the implemented public health intervention in managing infectious
travelers across borders from a time series of the number of new cases reported
in multiple geographical regions. It is found that the change point of the
probability of travel movements was one week after the WHO worldwide alert on
the SARS outbreak in 2003. The alert was effective in managing infectious
travelers. On the other hand, it is found that the probability of travel
movements did not change at all for the flu pandemic in 2009. The pandemic did
not affect potential travelers despite the WHO alert.


  We study the mean field approximation to a simple spatial host-pathogen model
that has been shown to display interesting evolutionary properties. We show
that previous derivations of the mean field equations for this model are
actually only low-density approximations to the true mean field limit. We
derive the correct equations and the corresponding equations including
pair-correlations. The process of invasion by a mutant type of pathogen is also
discussed. [This article was published as Physical Review E 67, 047102 (2003).
Errata for the published version are corrected here and explicitly listed at
the end of this document.]


  We performed a stability analysis of 119 quantitative food webs which were
compiled using a standard methodology to build Ecopath mass-balance models. Our
analysis reveals that classic descriptors of complexity do not affect stability
in natural food webs. Food web structure, which is non-random in real
communities, reflects another form of complexity that we found influences
dramatically the stability of real communities. We conclude that the occurrence
of complex communities in nature is possible owing to their trophic structure.


  Urea, a polar molecule with a large dipole moment, not only destabilizes the
folded RNA structures, but can also enhance the folding rates of large
ribozymes. Unlike the mechanism of urea-induced unfolding of proteins, which is
well understood, the action of urea on RNA has barely been explored. We
performed extensive all atom molecular dynamics (MD) simulations to determine
the molecular underpinnings of urea-induced RNA denaturation. Urea displays its
denaturing power in both secondary and tertiary motifs of the riboswitch (RS)
structure. Our simulations reveal that the denaturation of RNA structures is
mainly driven by the hydrogen bonds and stacking interactions of urea with the
bases. Through detailed studies of the simulation trajectories, we found that
geminate pairs between urea and bases due to hydrogen bonds and stacks persist
only ~ (0.1-1) ns, which suggests that urea-base interaction is highly dynamic.
Most importantly, the early stage of base pair disruption is triggered by
penetration of water molecules into the hydrophobic domain between the RNA
bases. The infiltration of water into the narrow space between base pairs is
critical in increasing the accessibility of urea to transiently disrupted
bases, thus allowing urea to displace inter base hydrogen bonds. This
mechanism, water-induced disruption of base-pairs resulting in the formation of
a "wet" destabilized RNA followed by solvation by urea, is the exact opposite
of the two-stage denaturation of proteins by urea. In the latter case, initial
urea penetration creates a dry-globule, which is subsequently solvated by water
penetration leading to global protein unfolding. Our work shows that the
ability to interact with both water and polar, non-polar components of
nucleotides makes urea a powerful chemical denaturant for nucleic acids.


  Flocks of birds exhibit a remarkable degree of coordination and collective
response. It is not just that thousands of individuals fly, on average, in the
same direction and at the same speed, but that even the fluctuations around the
mean velocity are correlated over long distances. Quantitative measurements on
flocks of starlings, in particular, show that these fluctuations are
scale-free, with effective correlation lengths proportional to the linear size
of the flock. Here we construct models for the joint distribution of velocities
in the flock that reproduce the observed local correlations between individuals
and their neighbors, as well as the variance of flight speeds across
individuals, but otherwise have as little structure as possible. These
minimally structured, or maximum entropy models provide quantitative,
parameter-free predictions for the spread of correlations throughout the flock,
and these are in excellent agreement with the data. These models are
mathematically equivalent to statistical physics models for ordering in
magnets, and the correct prediction of scale-free correlations arises because
the parameters - completely determined by the data - are in the critical
regime. In biological terms, criticality allows the flock to achieve maximal
correlation across long distances with limited speed fluctuations.


  RNA/protein interactions play crucial roles in controlling gene expression.
They are becoming important targets for pharmaceutical applications. Due to RNA
flexibility and to the strength of electrostatic interactions, standard docking
methods are insufficient. We here present a computational method which allows
studying the binding of RNA molecules and charged peptides with atomistic,
explicit-solvent molecular dynamics. In our method, a suitable estimate of the
electrostatic interaction is used as an order parameter (collective variable)
which is then accelerated using bi-directional pulling simulations. Since the
electrostatic interaction is only used to enhance the sampling, the
approximations used to compute it do not affect the final accuracy. The method
is employed to characterize the binding of TAR RNA from HIV-1 and a small
cyclic peptide. Our simulation protocol allows blindly predicting the binding
pocket and pose as well as the binding affinity. The method is general and
could be applied to study other electrostatics-driven binding events.


  We consider the continuous-time presentation of the strand symmetric
phylogenetic substitution model (in which rate parameters are unchanged under
nucleotide permutations given by Watson-Crick base conjugation). Algebraic
analysis of the model's underlying structure as a matrix group leads to a
change of basis where the rate generator matrix is given by a two-part block
decomposition. We apply representation theoretic techniques and, for any
(fixed) number of phylogenetic taxa $L$ and polynomial degree $D$ of interest,
provide the means to classify and enumerate the associated Markov invariants.
In particular, in the quadratic and cubic cases we prove there are precisely
1/3$(3^L+(-1)^L)$ and $6^{L-1}$ linearly independent Markov invariants,
respectively. Additionally, we give the explicit polynomial forms of the Markov
invariants for (i) the quadratic case with any number of taxa $L$, and (ii) the
cubic case in the special case of a three-taxa phylogenetic tree. We close by
showing our results are of practical interest since the quadratic Markov
invariants provide independent estimates of phylogenetic distances based on (i)
substitution rates within Watson-Crick conjugate pairs, and (ii) substitution
rates across conjugate base pairs.


  Collective behaviour is a widespread phenomenon in biology, cutting through a
huge span of scales, from cell colonies up to bird flocks and fish schools. The
most prominent trait of collective behaviour is the emergence of global order:
individuals synchronize their states, giving the stunning impression that the
group behaves as one. In many biological systems, though, it is unclear whether
global order is present. A paradigmatic case is that of insect swarms, whose
erratic movements seem to suggest that group formation is a mere epiphenomenon
of the independent interaction of each individual with an external landmark. In
these cases, whether or not the group behaves truly collectively is debated.
Here, we experimentally study swarms of midges in the field and measure how
much the change of direction of one midge affects that of other individuals. We
discover that, despite the lack of collective order, swarms display very strong
correlations, totally incompatible with models of noninteracting particles. We
find that correlation increases sharply with the swarm's density, indicating
that the interaction between midges is based on a metric perception mechanism.
By means of numerical simulations we demonstrate that such growing correlation
is typical of a system close to an ordering transition. Our findings suggest
that correlation, rather than order, is the true hallmark of collective
behaviour in biological systems.


  We have shown that commonly used model-fitting methods fail when the
underlying ecological system is chaotic. Importantly, Hartig & Dormann
[arXiv:1305.3544] suggest an alternative method that is able to overcome this
problem for the logistic model when the "true model" is known. This is
encouraging, however the more relevant test is of the method when the "true
model" is unknown, as it always is for real ecosystems. By making a small
modification to the "true model," we find that the suggested method provides
inaccurate forecasts and severely biased estimates of the ecological dynamics.
In contrast, a model-free approach is able to provide accurate forecasts and a
useful foundation for ecological understanding without the burdensome
assumption of a "true model."


  A major challenge in single particle reconstruction methods using
cryo-electron microscopy is to attain a resolution sufficient to interpret fine
details in three-dimensional (3D) macromolecular structures. Obtaining high
resolution 3D reconstructions is difficult due to unknown orientations and
positions of the imaged particles, possible incomplete coverage of the viewing
directions, high level of noise in the projection images, and limiting effects
of the contrast transfer function of the electron microscope. In this paper, we
focus on the 3D reconstruction problem from projection images assuming an
existing estimate for their orientations and positions. We propose a fast and
accurate Fourier-based Iterative Reconstruction Method (FIRM) that exploits the
Toeplitz structure of the operator ${\bf A}^{*}{\bf A}$, where $\bf A$ is the
forward projector and ${\bf A}^{*}$ is the back projector. The operator ${\bf
A}^{*}{\bf A}$ is equivalent to a convolution with a kernel. The kernel is
pre-computed using the non-uniform Fast Fourier Transform and is efficiently
applied in each iteration step. The iterations by FIRM are therefore
considerably faster than those of traditional iterative algebraic approaches,
while maintaining the same accuracy even when the viewing directions are
unevenly distributed. The time complexity of FIRM is comparable to the direct
Fourier inversion method. Moreover, FIRM combines images from different defocus
groups simultaneously and can handle a wide range of regularization terms. We
provide experimental results on simulated data that demonstrate the speed and
accuracy of FIRM in comparison with current methods.


  We analyze experimental observations of microtubules undergoing small
fluctuations about a "balance point" when mixed in solution of two different
kinesin motor proteins, KLP61F and Ncd. It has been proposed that the
microtubule movement is due to stochastic variations in the densities of the
two species of motor proteins. We test this hypothesis here by showing how it
maps onto a one-dimensional random walk in a random environment. Our estimate
of the amplitude of the fluctuations agrees with experimental observations. We
point out that there is an initial transient in the position of the microtubule
where it will typically move of order its own length. We compare the physics of
this gliding assay to a recent theory of the role of antagonistic motors on
restricting interpolar microtubule sliding of a cell's mitotic spindle during
prometaphase. It is concluded that randomly positioned antagonistic motors can
restrict relative movement of microtubules, however they do so imperfectly. A
variation in motor concentrations is also analyzed and shown to lead to greater
control of spindle length.


  Distance-based phylogenetic algorithms attempt to solve the NP-hard least
squares phylogeny problem by mapping an arbitrary dissimilarity map
representing biological data to a tree metric. The set of all dissimilarity
maps is a Euclidean space properly containing the space of all tree metrics as
a polyhedral fan. Outputs of distance-based tree reconstruction algorithms such
as UPGMA and Neighbor-Joining are points in the maximal cones in the fan. Tree
metrics with polytomies lie at the intersections of maximal cones.
  A phylogenetic algorithm divides the space of all dissimilarity maps into
regions based upon which combinatorial tree is reconstructed by the algorithm.
Comparison of phylogenetic methods can be done by comparing the geometry of
these regions. We use polyhedral geometry to compare the local nature of the
subdivisions induced by least squares phylogeny, UPGMA, and Neighbor-Joining.
Our results suggest that in some circumstances, UPGMA and Neighbor-Joining
poorly match least squares phylogeny when the true tree has a polytomy.


  The GC-content is very variable in different genome regions and species but
although many hypothesis we still do not know the reason why. Here we show that
a relationship exists with the mutation rate, in particular we noticed a new
recurrence in the amino acids coding table. Moreover we analyze recombination
frequency taking into account Single Nucleotide Polymorphisms hourglass
distribution.


  Models for chemotaxis are based on gradient sensing of individual organisms.
The key contribution of Keller and Segel is showing that erratic movements of
individuals may result in an accurate chemotaxis phenomenon as a group. In this
paper we provide another option to understand chemotactic behavior when
individuals do not sense the gradient of chemical concentration by any means.
We show that, if individuals increase their motility to find food when they are
hungry, an accurate chemotactic behavior may obtained without sensing the
gradient. Such a random dispersal has been suggested by Cho and Kim and is
called starvation driven diffusion. This model is surprisingly similar to the
original derivation of Keller-Segel model. A comprehensive picture of traveling
band and front solutions is provided with numerical simulations.


  The attack of novel avian influenza (H7N9) in east China caused a serious
health crisis and public panic. In this paper, we empirically analyze the onset
patterns of human cases of the novel avian influenza and observe several
spatial and temporal properties that are similar to other infective diseases.
More deeply, using the empirical analysis and modeling studies, we find that
the spatio-temporal network that connects the cities with human cases along the
order of outbreak timing emerges two-section-power-law edge-length
distribution, indicating the picture that several islands with higher and
heterogeneous risk straggle in east China. The proposed method is applicable to
the analysis on the spreading situation in early stage of disease outbreak
using quite limited dataset.


  Development of the new antimicrobial agents against antibiotic resistance
pathogens is the nowadays challenge. Antimicrobial peptides (AMP) occur as
important defence agents in many organisms and offer a viable alternative to
conventional antibiotics. Therefore they have become increasingly recognized in
current research as templates for prospective antibiotic agents. The efficient
designing of the new antimicrobials on the basis of antimicrobial peptides
requires comprehensive knowledge on those general physical-chemical
characteristics which allow to differ antimicrobial peptides from non-active
against microbs ones. According to supposed mechanisms of action, AMP interact
with and physically disrupt the bacterial membranes. Consequently,
hydrophobicity, amphiphilicity and intrinsic aggregation propensities are
considered as such major characteristics of the peptide, which determine the
results of peptide-membrane interactions. For some kind of peptides such
characteristics as hydrophobicity, amphiphilicity and aggregation bias
determines their ability to compose transmembrane domain of the membrane
protein, whilst for others the same properties are respond for their
antimicrpobial activity, i.e. give them ability of membrane permeability and
its damage. In this review we analyze the data about hydrophobicity,
amphiphilicity and intrinsic aggregation propensities available in literature
in order to compare antimicrobial and transmembrane peptides and show what is
the common and what is the difference in this respect between them.


  Establishing the conditions that guarantee the spreading or the sustenance of
altruistic traits in a population is the main goal of intergroup selection
models. Of particular interest is the balance of the parameters associated to
group size, migration and group survival against the selective advantage of the
non-altruistic individuals. Here we use Kimura's diffusion model of intergroup
selection to determine those conditions in the case the group survival
probability is a nonlinear non-decreasing function of the proportion of
altruists in a group. In the case this function is linear, there are two
possible steady states which correspond to the non-altruistic and the
altruistic phases. At the discontinuous transition line separating these phases
there is a non-ergodic coexistence phase. For a continuous concave survival
function, we find an ergodic coexistence phase that occupies a finite region of
the parameter space in between the altruistic and the non-altruistic phases,
and is separated from these phases by continuous transition lines. For a convex
survival function, the coexistence phase disappears altogether but a bistable
phase appears for which the choice of the initial condition determines whether
the evolutionary dynamics leads to the altruistic or the non-altruistic steady
state.


  In the past decade, transcriptome data have become an important component of
many phylogenetic studies. Phylogenetic studies now regularly include genes
from newly sequenced transcriptomes, as well as publicly available
transcriptomes and genomes. Implementing such a phylogenomic study, however, is
computationally intensive, requires the coordinated use of many complex
software tools, and includes multiple steps for which no published tools exist.
Phylogenomic studies have therefore been manual or semiautomated. In addition
to taking considerable user time, this makes phylogenomic analyses difficult to
reproduce, compare, and extend. In addition, methodological improvements made
in the context of one study often cannot be easily applied and evaluated in the
context of other studies. We present Agalma, an automated tool that conducts
phylogenomic analyses. The user provides raw Illumina transcriptome data, and
Agalma produces annotated assemblies, aligned gene sequence matrices, a
preliminary phylogeny, and detailed diagnostics that allow the investigator to
make extensive assessments of intermediate analysis steps and the final
results. Sequences from other sources, such as externally assembled genomes and
transcriptomes, can also be incorporated in the analyses. Agalma tracks
provenance, profiles processor and memory use, records diagnostics, manages
metadata, and enables rich HTML reports for all stages of the analysis. Agalma
includes a test data set and a built-in test analysis of these data. In
addition to describing Agalma, we here present a sample analysis of a larger
seven-taxon data set. Agalma is available for download at
https://bitbucket.org/caseywdunn/agalma. Agalma allows complex phylogenomic
analyses to be implemented and described unambiguously as a series of
high-level commands. This will enable phylogenomic studies to be readily
reproduced, modified, and extended.


  We investigate what structural aspects of a collection of twelve empirical
temporal networks of human contacts are important to disease spreading. We scan
the entire parameter spaces of the two canonical models of infectious disease
epidemiology -- the Susceptible-Infectious-Susceptible (SIS) and
Susceptible-Infectious-Removed (SIR) models. The results from these simulations
are compared to reference data where we eliminate structures in the interevent
intervals, the time to the first contact in the data, or the time from the last
contact to the end of the sampling. The picture we find is that the birth and
death of links, and the total number of contacts over a link, are essential to
predict outbreaks. On the other hand, the exact times of contacts between the
beginning and end, or the interevent interval distribution, do not matter much.
In other words, a simplified picture of these empirical data sets that suffices
for epidemiological purposes is that links are born, is active with some
intensity, and die.


  Moment control of a simple quadratic reaction network describing a
dimerization process is addressed. It is shown that the moment closure problem
can be circumvented without invoking any moment closure technique. Local
stabilization and convergence of the average dimer population to any desired
reference value is ensured using a pure integral control law. Explicit bounds
on the controller gain are provided and shown to be valid for any reference
value. As a byproduct, an explicit upper-bound of the variance of the monomer
species, acting on the system as unknown input due to the moment openness, is
obtained. The obtained results are illustrated by an example relying on the
simulation of a cell population using stochastic simulation algorithms.


  Negatively charged DNA can be compacted by positively charged dendrimers and
the degree of compaction is a delicate balance between the strength of the
electrostatic interaction and the elasticity of DNA. We report various elastic
properties of short double stranded DNA (dsDNA) and the effect of dendrimer
binding using fully atomistic molecular dynamics and numerical simulations. In
equilibrium at room temperature, the contour length distribution P(L) and
end-to-end distance distribution P(R) are nearly Gaussian, the former gives an
estimate of the stretch modulus {\gamma}_1 of dsDNA in quantitative agreement
with the literature value. The bend angle distribution P({\theta}) of the dsDNA
also has a Gaussian form and allows to extract a persistence length, L_p of 43
nm. When the dsDNA is compacted by positively charged dendrimer, the stretch
modulus stays invariant but the effective bending rigidity estimated from the
end-to-end distance distribution decreases dramatically due to backbone charge
neutralization of dsDNA by dendrimer. We support our observations with
numerical solutions of the worm-like-chain (WLC) model as well as using
non-equilibrium dsDNA stretching simulations. These results are helpful in
understanding the dsDNA elasticity at short length scales as well as how the
elasticity is modulated when dsDNA binds to a charged object such as a
dendrimer or protein.


  Non-processive molecular motors have to work together in ensembles in order
to generate appreciable levels of force or movement. In skeletal muscle, for
example, hundreds of myosin II molecules cooperate in thick filaments. In
non-muscle cells, by contrast, small groups with few tens of non-muscle myosin
II motors contribute to essential cellular processes such as transport, shape
changes or mechanosensing. Here we introduce a detailed and analytically
tractable model for this important situation. Using a three-state crossbridge
model for the myosin II motor cycle and exploiting the assumptions of fast
power stroke kinetics and equal load sharing between motors in equivalent
states, we reduce the stochastic reaction network to a one-step master equation
for the binding and unbinding dynamics (parallel cluster model) and derive the
rules for ensemble movement. We find that for constant external load, ensemble
dynamics is strongly shaped by the catch bond character of myosin II, which
leads to an increase of the fraction of bound motors under load and thus to
firm attachment even for small ensembles. This adaptation to load results in a
concave force-velocity relation described by a Hill relation. For external load
provided by a linear spring, myosin II ensembles dynamically adjust themselves
towards an isometric state with constant average position and load. The
dynamics of the ensembles is now determined mainly by the distribution of
motors over the different kinds of bound states. For increasing stiffness of
the external spring, there is a sharp transition beyond which myosin II can no
longer perform the power stroke. Slow unbinding from the pre-power-stroke state
protects the ensembles against detachment.


  Evolutionary forces shape patterns of genetic diversity within populations
and contribute to phenotypic variation. In particular, recurrent positive
selection has attracted significant interest in both theoretical and empirical
studies. However, most existing theoretical models of recurrent positive
selection cannot easily incorporate realistic confounding effects such as
interference between selected sites, arbitrary selection schemes, and
complicated demographic processes. It is possible to quantify the effects of
arbitrarily complex evolutionary models by performing forward population
genetic simulations, but forward simulations can be computationally prohibitive
for large population sizes ($> 10^5$). A common approach for overcoming these
computational limitations is rescaling of the most computationally expensive
parameters, especially population size. Here, we show that ad hoc approaches to
parameter rescaling under the recurrent hitchhiking model do not always provide
sufficiently accurate dynamics, potentially skewing patterns of diversity in
simulated DNA sequences. We derive an extension of the recurrent hitchhiking
model that is appropriate for strong selection in small population sizes, and
use it to develop a method for parameter rescaling that provides the best
possible computational performance for a given error tolerance. We perform a
detailed theoretical analysis of the robustness of rescaling across the
parameter space. Finally, we apply our rescaling algorithms to parameters that
were previously inferred for Drosophila, and discuss practical considerations
such as interference between selected sites.


  Using an Ising-like model of protein mechanical unfolding, we introduce a
diffusive dynamics on its exactly known free energy profile, reducing the
nonequilibrium dynamics of the model to a biased random walk. As an
illustration, the model is then applied to the protein translocation
phenomenon, taking inspiration from a recent experiment on the green
fluorescent protein pulled by a molecular motor. The average translocation time
is evaluated exactly, and the analysis of single trajectories shows that
translocation proceeds through an intermediate state, similar to that observed
in the experiment.


  We show that the generation time -- a notion usually described in a
biological context -- can be defined in a general way as a return time in a
conveniently constructed finite Markov chain. The simple formula we obtain
agrees with previous results derived for structured populations projected in
discrete time, and allows to define the generation time of any process
described by a weighted directed graph whose matrix is primitive.


  Coevolution is expected to follow two alternative dynamics, often called
trench warfare and arms races in plant-pathogen systems. Trench warfare
situations are stable cycles of allele frequencies at the coevolving loci of
both host and parasite, and it is predicted that the loci will show molecular
evolutionary signatures of balancing selection, while arms races involve
successive selective sweeps at the interacting loci. We study a haploid
gene-for-gene model that includes mutation and genetic drift due to finite
population size. We study the outcomes under different coevolutionary
parameters to quantify the frequency of fixation of alleles, i.e. occurrence of
an arms race dynamics. We find that contrary to the conventional wisdom, trench
warfare situations do not imply larger numbers of coevolutionary cycles per
unit time than arms races. Therefore, one cannot infer the nature of the
dynamics in such systems based on the speed of coevolution estimated from cycle
times. We subsequently perform coalescent simulations to generate sequences at
the host and parasite loci. We ask whether the signatures expected under
balancing selection or selective sweeps (unexpectedly high or low diversity,
and high or low Tajima D values, respectively) are likely to be observable in
genomic data. Genomic footprints of recurrent selective sweeps are often found,
whereas trench warfare yields signatures of balancing selection only in
parasite sequences, and only in a limited parameter space with high effective
population sizes and long-term selection. Therefore, the existence of a
deterministic polymorphic equilibrium does not imply long-term trench warfare
necessary for the signature of balancing selection to be observed in the
coevolving genes sequence. Our results suggest that to search for signatures of
coevolution via population genomics, it is best to study pathogen rather than
host genomes.


  Energy landscape theory describes how a full-length protein can attain its
native fold after sampling only a tiny fraction of all possible structures.
Although protein folding is now understood to be concomitant with synthesis on
the ribosome there have been few attempts to modify energy landscape theory by
accounting for cotranslational folding. This paper introduces a model for
cotranslational folding that leads to a natural definition of a nested energy
landscape. By applying concepts drawn from submanifold differential geometry
the dynamics of protein folding on the ribosome can be explored in a
quantitative manner and conditions on the nested potential energy landscapes
for a good cotranslational folder are obtained. A generalisation of diffusion
rate theory using van Kampen's technique of composite stochastic processes is
then used to account for entropic contributions and the effects of variable
translation rates on cotranslational folding. This stochastic approach agrees
well with experimental results and Hamiltionian formalism in the deterministic
limit.


  Many DNA sequence variants influence phenotypes by altering gene expression.
Our understanding of these variants is limited by sample sizes of current
studies and by measurements of mRNA rather than protein abundance. We developed
a powerful method for identifying genetic loci that influence protein
expression in very large populations of the yeast Saccharomyes cerevisiae. The
method measures single-cell protein abundance through the use of
green-fluorescent-protein tags. We applied this method to 160 genes and
detected many more loci per gene than previous studies. We also observed closer
correspondence between loci that influence protein abundance and loci that
influence mRNA abundance of a given gene. Most loci cluster at hotspot
locations that influence multiple proteins - in some cases, more than half of
those examined. The variants that underlie these hotspots have profound effects
on the gene regulatory network and provide insights into genetic variation in
cell physiology between yeast strains.


  Background: Analysing tumour architecture for metastatic potential usually
focuses on phenotypic differences due to cellular morphology or specific
genetic mutations, but often ignore the cell's position within the
heterogeneous substructure. Similar disregard for local neighborhood structure
is common in mathematical models.
  Methods: We view the dynamics of disease progression as an evolutionary game
between cellular phenotypes. A typical assumption in this modeling paradigm is
that the probability of a given phenotypic strategy interacting with another
depends exclusively on the abundance of those strategies without regard local
heterogeneities. We address this limitation by using the Ohtsuki-Nowak
transform to introduce spatial structure to the go vs. grow game.
  Results: We show that spatial structure can promote the invasive (go)
strategy. By considering the change in neighbourhood size at a static boundary
-- such as a blood-vessel, organ capsule, or basement membrane -- we show an
edge effect that allows a tumour without invasive phenotypes in the bulk to
have a polyclonal boundary with invasive cells. We present an example of this
promotion of invasive (EMT positive) cells in a metastatic colony of prostate
adenocarcinoma in bone marrow.
  Interpretation: Pathologic analyses that do not distinguish between cells in
the bulk and cells at a static edge of a tumour can underestimate the number of
invasive cells. We expect our approach to extend to other evolutionary game
models where interaction neighborhoods change at fixed system boundaries.


  Background: Systems biology projects and omics technologies have led to a
growing number of biochemical pathway reconstructions. However, mathematical
models are still most often created de novo, based on reading the literature
and processing pathway data manually. Results: To increase the efficiency with
which such models can be created, we automatically generated mathematical
models from pathway representations using a suite of freely available software.
We produced models that combine data from KEGG PATHWAY, BioCarta, MetaCyc and
SABIO-RK; According to the source data, three types of models are provided:
kinetic, logical and constraint-based. All models are encoded using SBML Core
and Qual packages, and available through BioModels Database. Each model
contains the list of participants, the interactions, and the relevant
mathematical constructs, but, in most cases, no meaningful parameter values.
Most models are also available as easy to understand graphical SBGN maps.
Conclusions: to date, the project has resulted in more than 140000 models
freely available. We believe this resource can tremendously accelerate the
development of mathematical models by providing initial starting points ready
for parametrization.


  The mechanisms controlling the transport of proteins across the Golgi stack
of mammalian and plant cells is the subject of intense debate, with two models,
cisternal progression and inter-cisternal exchange, emerging as major
contenders. A variety of transport experiments have claimed support for each of
these models. We reevaluate these experiments using a single quantitative
coarse-grained framework of intra-Golgi transport that accounts for both
transport models and their many variants. Our analysis makes a definitive case
for the existence of inter-cisternal exchange both for small membrane proteins
(VSVG) and large protein complexes (procollagen) -- this implies that membrane
structures larger than the typical protein-coated vesicles must be involved in
transport. Notwithstanding, we find that current observations on protein
transport cannot rule out cisternal progression as contributing significantly
to the transport process. To discriminate between the different models of
intra-Golgi transport, we suggest experiments and an analysis based on our
extended theoretical framework that compare the dynamics of transiting and
resident proteins.


  We discuss several models of the dynamics of interacting populations. The
models are constructed by nonlinear differential equations and have two sets of
parameters: growth rates and coefficients of interaction between populations.
We assume that the parameters depend on the densities of the populations. In
addition the parameters can be influenced by different factors of the
environment. This influence is modelled by noise terms in the equations for the
growth rates and interaction coefficients. Thus the model differential
equations become stochastic. In some particular cases these equations can be
reduced to a Foker-Plancnk equation for the probability density function of the
densities of the interacting populations.


  Understanding how genetic variants influence cellular-level processes is an
important step towards understanding how they influence important
organismal-level traits, or "phenotypes", including human disease
susceptibility. To this end scientists are undertaking large-scale genetic
association studies that aim to identify genetic variants associated with
molecular and cellular phenotypes, such as gene expression, transcription
factor binding, or chromatin accessibility. These studies use high-throughput
sequencing assays (e.g. RNA-seq, ChIP-seq, DNase-seq) to obtain high-resolution
data on how the traits vary along the genome in each sample. However, typical
association analyses fail to exploit these high-resolution measurements,
instead aggregating the data at coarser resolutions, such as genes, or windows
of fixed length. Here we develop and apply statistical methods that better
exploit the high-resolution data. The key idea is to treat the sequence data as
measuring an underlying "function" that varies along the genome, and then,
building on wavelet-based methods for functional data analysis, test for
association between genetic variants and the underlying function. Applying
these methods to identify genetic variants associated with chromatin
accessibility (dsQTLs) we find that they identify substantially more
associations than a simpler window-based analysis, and in total we identify 772
novel dsQTLs not identified by the original analysis.


  Modern maize breeding relies upon selection in inbreeding populations to
improve the performance of cross-population hybrids. The United States
Department of Agriculture - Agricultural Research Service reciprocal recurrent
selection experiment between the Iowa Stiff Stalk Synthetic (BSSS) and the Iowa
Corn Borer Synthetic No. 1 (BSCB1) populations represents one of the longest
standing models of selection for hybrid performance. To investigate the genomic
impact of this selection program, we used the Illumina MaizeSNP50 high-density
SNP array to determine genotypes of progenitor lines and over 600 individuals
across multiple cycles of selection. Consistent with previous research (Messmer
et al., 1991; Labate et al., 1997; Hagdorn et al., 2003; Hinze et al., 2005),
we found that genetic diversity within each population steadily decreases, with
a corresponding increase in population structure. High marker density also
enabled the first view of haplotype ancestry, fixation and recombination within
this historic maize experiment. Extensive regions of haplotype fixation within
each population are visible in the pericentromeric regions, where large blocks
trace back to single founder inbreds. Simulation attributes most of the
observed reduction in genetic diversity to genetic drift. Signatures of
selection were difficult to observe in the background of this strong genetic
drift, but heterozygosity in each population has fallen more than expected.
Regions of haplotype fixation represent the most likely targets of selection,
but as observed in other germplasm selected for hybrid performance (Feng et
al., 2006), there is no overlap between the most likely targets of selection in
the two populations. We discuss how this pattern is likely to occur during
selection for hybrid performance, and how it poses challenges for dissecting
the impacts of modern breeding and selection on the maize genome.


  We explore the effect of an attractive interaction between parallel-aligned
polymers, which are perpendicularly grafted on a substrate. Such an attractive
interaction could be due to, e.g., reversible cross-links. The competition
between permanent grafting favoring a homogeneous state of the polymer brush
and the attraction, which tends to induce in-plane collapse of the aligned
polymers, gives rise to an instability of the homogeneous phase to a bundled
state. In this latter state the in-plane translational symmetry is
spontaneously broken and the density is modulated with a finite wavelength,
which is set by the length scale of transverse fluctuations of the grafted
polymers. We analyze the instability for two models of aligned polymers:
directed polymers with a line tension and weakly bending chains with a bending
stiffness.


  Drug target identification is of significant commercial interest to
pharmaceutical companies, and there is a vast amount of research done related
to the topic of therapeutic target identification. Interdisciplinary research
in this area involves both the biological network community and the graph
algorithms community. Key steps of a typical therapeutic target identification
problem include synthesizing or inferring the complex network of interactions
relevant to the disease, connecting this network to the disease-specific
behavior, and predicting which components are key mediators of the behavior.
All of these steps involve graph theoretical or graph algorithmic aspects. In
this perspective, we provide modelling and algorithmic perspectives for
therapeutic target identification and highlight a number of algorithmic
advances, which have gotten relatively little attention so far, with the hope
of strengthening the ties between these two research communities.


  The genome of a 650 year old Yersinia pestis bacteria, responsible for the
medieval Black Death, was recently sequenced and assembled into 2,105 contigs
from the main chromosome. According to the point mutation record, the medieval
bacteria could be an ancestor of most Yersinia pestis extant species, which
opens the way to reconstructing the organization of these contigs using a
comparative approach. We show that recent computational paleogenomics methods,
aiming at reconstructing the organization of ancestral genomes from the
comparison of extant genomes, can be used to correct, order and complete the
contig set of the Black Death agent genome, providing a full chromosome
sequence, at the nucleotide scale, of this ancient bacteria. This sequence
suggests that a burst of mobile elements insertions predated the Black Death,
leading to an exceptional genome plasticity and increase in rearrangement rate.


  In plants, a subset of genes exhibit imprinting in endosperm tissue such that
expression is primarily from the maternal or paternal allele. Imprinting may
arise as a consequence of mechanisms for silencing of transposons during
reproduction, and in some cases imprinted expression of particular genes may
provide a selective advantage such that it is conserved across species.
Separate mechanisms for the origin of imprinted expression patterns and
maintenance of these patterns may result in substantial variation in the
targets of imprinting in different species. Here we present deep sequencing of
RNAs isolated from reciprocal crosses of four diverse maize genotypes,
providing a comprehensive analysis of imprinting in maize that allows
evaluation of imprinting at more than 95% of endosperm-expressed genes. We find
that over 500 genes exhibit statistically significant parent-of-origin effects
in maize endosperm tissue, but focused our analyses on a subset of these genes
that had >90% expression from the maternal allele (69 genes) or from the
paternal allele (108 genes) in at least one reciprocal cross. Over 10% of
imprinted genes show evidence of allelic variation for imprinting. A comparison
of imprinting in maize and rice reveals that only 13% of genes with syntenic
orthologs in both species exhibit conserved imprinting. Genes that exhibit
conserved imprinting in maize relative to rice have elevated dN/dS ratios
compared to other imprinted genes, suggesting a history of more rapid
evolution. Together, these data suggest that imprinting only has functional
relevance at a subset of loci that currently exhibit imprinting in maize.


  In a previous paper Sharkey et al. [13] proved the exactness of closures at
the level of triples for Markovian SIR (susceptible-infected-removed) dynamics
on tree-like networks. This resulted in a deterministic representation of the
epidemic dynamics on the network that can be numerically evaluated. In this
paper, we extend this modelling framework to certain classes of networks
exhibiting loops. We show that closures where the loops are kept intact are
exact, and lead to a simplified and numerically solvable system of ODEs
(ordinary-differential-equations). The findings of the paper lead us to a
generalisation of closures that are based on partitioning the network around
nodes that are cut-vertices (i.e. the removal of such a node leads to the
network breaking down into at least two disjointed components or subnetworks).
Exploiting this structural property of the network yields some natural
closures, where the evolution of a particular state can typically be exactly
given in terms of the corresponding or projected sates on the subnetworks and
the cut-vertex. A byproduct of this analysis is an alternative probabilistic
proof of the exactness of the closures for tree-like networks presented in
Sharkey et al. [13]. In this paper we also elaborate on how the main result can
be applied to more realistic networks, for which we write down the ODEs
explicitly and compare output from these to results from simulation.
Furthermore, we give a general, recipe-like method of how to apply the
reduction by closures technique for arbitrary networks, and give an upper bound
on the maximum number of equations needed for an exact representation.


  The Wright-Fisher process with selection is an important tool in population
genetics theory. Traditional analysis of this process relies on the diffusion
approximation. The diffusion approximation is usually studied in a partial
differential equations framework. In this paper, I introduce a path integral
formalism to study the Wright-Fisher process with selection and use that
formalism to obtain a simple perturbation series to approximate the transition
density. The perturbation series can be understood in terms of Feynman
diagrams, which have a simple probabilistic interpretation in terms of
selective events. The perturbation series proves to be an accurate
approximation of the transition density for weak selection and is shown to be
arbitrarily accurate for any selection coefficient.


  Adaptation in response to selection on polygenic phenotypes may occur via
subtle allele frequencies shifts at many loci. Current population genomic
techniques are not well posed to identify such signals. In the past decade,
detailed knowledge about the specific loci underlying polygenic traits has
begun to emerge from genome-wide association studies (GWAS). Here we combine
this knowledge from GWAS with robust population genetic modeling to identify
traits that may have been influenced by local adaptation. We exploit the fact
that GWAS provide an estimate of the additive effect size of many loci to
estimate the mean additive genetic value for a given phenotype across many
populations as simple weighted sums of allele frequencies. We first describe a
general model of neutral genetic value drift for an arbitrary number of
populations with an arbitrary relatedness structure. Based on this model we
develop methods for detecting unusually strong correlations between genetic
values and specific environmental variables, as well as a generalization of
$Q_{ST}/F_{ST}$ comparisons to test for over-dispersion of genetic values among
populations. Finally we lay out a framework to identify the individual
populations or groups of populations that contribute to the signal of
overdispersion. These tests have considerably greater power than their single
locus equivalents due to the fact that they look for positive covariance
between like effect alleles, and also significantly outperform methods that do
not account for population structure. We apply our tests to the Human Genome
Diversity Panel (HGDP) dataset using GWAS data for height, skin pigmentation,
type 2 diabetes, body mass index, and two inflammatory bowel disease datasets.
This analysis uncovers a number of putative signals of local adaptation, and we
discuss the biological interpretation and caveats of these results.


  Knowing the location of a protein within the cell is important for
understanding its function, role in biological processes, and potential use as
a drug target. Much progress has been made in developing computational methods
that predict single locations for proteins, assuming that proteins localize to
a single location. However, it has been shown that proteins localize to
multiple locations. While a few recent systems have attempted to predict
multiple locations of proteins, they typically treat locations as independent
or capture inter-dependencies by treating each locations-combination present in
the training set as an individual location-class. We present a new method and a
preliminary system we have developed that directly incorporates
inter-dependencies among locations into the multiple-location-prediction
process, using a collection of Bayesian network classifiers. We evaluate our
system on a dataset of single- and multi-localized proteins. Our results,
obtained by incorporating inter-dependencies are significantly higher than
those obtained by classifiers that do not use inter-dependencies. The
performance of our system on multi-localized proteins is comparable to a top
performing system (YLoc+), without restricting predictions to be based only on
location-combinations present in the training set.


  The sheer amounts of biological data that are generated in recent years have
driven the development of network analysis tools to facilitate the
interpretation and representation of these data. A fundamental challenge in
this domain is the reconstruction of a protein-protein subnetwork that
underlies a process of interest from a genome-wide screen of associated genes.
Despite intense work in this area, current algorithmic approaches are largely
limited to analyzing a single screen and are, thus, unable to account for
information on condition-specific genes, or reveal the dynamics (over time or
condition) of the process in question. Here we propose a novel formulation for
network reconstruction from multiple-condition data and devise an efficient
integer program solution for it. We apply our algorithm to analyze the response
to influenza infection in humans over time as well as to analyze a pair of ER
export related screens in humans. By comparing to an extant, single-condition
tool we demonstrate the power of our new approach in integrating data from
multiple conditions in a compact and coherent manner, capturing the dynamics of
the underlying processes.


  In order to overcome the limitations imposed by DNA barcoding when
multiplexing a large number of samples in the current generation of
high-throughput sequencing instruments, we have recently proposed a new
protocol that leverages advances in combinatorial pooling design (group
testing) doi:10.1371/journal.pcbi.1003010. We have also demonstrated how this
new protocol would enable de novo selective sequencing and assembly of large,
highly-repetitive genomes. Here we address the problem of decoding pooled
sequenced data obtained from such a protocol. Our algorithm employs a
synergistic combination of ideas from compressed sensing and the decoding of
error-correcting codes. Experimental results on synthetic data for the rice
genome and real data for the barley genome show that our novel decoding
algorithm enables significantly higher quality assemblies than the previous
approach.


  Large RNA molecules often carry multiple functional domains whose spatial
arrangement is an important determinant of their function. Pre-mRNA splicing,
furthermore, relies on the spatial proximity of the splice junctions that can
be separated by very long introns. Similar effects appear in the processing of
RNA virus genomes. Albeit a crude measure, the distribution of spatial
distances in thermodynamic equilibrium therefore provides useful information on
the overall shape of the molecule can provide insights into the interplay of
its functional domains. Spatial distance can be approximated by the
graph-distance in RNA secondary structure. We show here that the equilibrium
distribution of graph-distances between arbitrary nucleotides can be computed
in polynomial time by means of dynamic programming. A naive implementation
would yield recursions with a very high time complexity of O(n^11). Although we
were able to reduce this to O(n^6) for many practical applications a further
reduction seems difficult. We conclude, therefore, that sampling approaches,
which are much easier to implement, are also theoretically favorable for most
real-life applications, in particular since these primarily concern long-range
interactions in very large RNA molecules.


  Evolution is a process that is influenced by various environmental factors,
e.g. the interactions between different species, genes, and biogeographical
properties. Hence, it is interesting to study the combined evolutionary history
of multiple species, their genes, and the environment they live in. A common
approach to address this research problem is to describe each individual
evolution as a phylogenetic tree and construct a tree reconciliation which is
parsimonious with respect to a given event model. Unfortunately, most of the
previous approaches are designed only either for host-parasite systems, for
gene tree/species tree reconciliation, or biogeography. Hence, a method is
desirable, which addresses the general problem of mapping phylogenetic trees
and covering all varieties of coevolving systems, including e.g., predator-prey
and symbiotic relationships. To overcome this gap, we introduce a generalized
cophylogenetic event model considering the combinatorial complete set of local
coevolutionary events. We give a dynamic programming based heuristic for
solving the maximum parsimony reconciliation problem in time O(n^2), for two
phylogenies each with at most n leaves. Furthermore, we present an exact
branch-and-bound algorithm which uses the results from the dynamic programming
heuristic for discarding partial reconciliations. The approach has been
implemented as a Java application which is freely available from
http://pacosy.informatik.uni-leipzig.de/coresym.


  We propose a simple tractable pair hidden Markov model for pairwise sequence
alignment that accounts for the presence of short tandem repeats. Using the
framework of gain functions, we design several optimization criteria for
decoding this model and describe the resulting decoding algorithms, ranging
from the traditional Viterbi and posterior decoding to block-based decoding
algorithms specialized for our model. We compare the accuracy of individual
decoding algorithms on simulated data and find our approach superior to the
classical three-state pair HMM in simulations.


  Recent chromosome conformation capture experiments have led to the discovery
of dense, contiguous, megabase-sized topological domains that are similar
across cell types and conserved across species. These domains are strongly
correlated with a number of chromatin markers and have since been included in a
number of analyses. However, functionally-relevant domains may exist at
multiple length scales. We introduce a new and efficient algorithm that is able
to capture persistent domains across various resolutions by adjusting a single
scale parameter. The identified novel domains are substantially different from
domains reported previously and are highly enriched for insulating factor CTCF
binding and histone modfications at the boundaries.


  We study a stochastic matrix model to understand the mechanics of
risk-spreading (or bet-hedging) by dispersion. Such model has been mostly dealt
numerically except for well-mixed case, so far. Here, we present an analytical
result, which shows that optimal dispersion leads to Zipf's law. Moreover, we
found that the arithmetic ensemble average of the total growth rate converges
to the geometric one, because the sample size is finite.


  Genome assembly using high throughput data with short reads, arguably,
remains an unresolvable task in repetitive genomes, since when the length of a
repeat exceeds the read length, it becomes difficult to unambiguously connect
the flanking regions. The emergence of third generation sequencing (Pacific
Biosciences) with long reads enables the opportunity to resolve complicated
repeats that could not be resolved by the short read data. However, these long
reads have high error rate and it is an uphill task to assemble the genome
without using additional high quality short reads. Recently, Koren et al. 2012
proposed an approach to use high quality short reads data to correct these long
reads and, thus, make the assembly from long reads possible. However, due to
the large size of both dataset (short and long reads), error-correction of
these long reads requires excessively high computational resources, even on
small bacterial genomes. In this work, instead of error correction of long
reads, we first assemble the short reads and later map these long reads on the
assembly graph to resolve repeats.
  Contribution: We present a hybrid assembly approach that is both
computationally effective and produces high quality assemblies. Our algorithm
first operates with a simplified version of the assembly graph consisting only
of long contigs and gradually improves the assembly by adding smaller contigs
in each iteration. In contrast to the state-of-the-art long reads error
correction technique, which requires high computational resources and long
running time on a supercomputer even for bacterial genome datasets, our
software can produce comparable assembly using only a standard desktop in a
short running time.


  Comparing strains within the same microbial species has proven effective in
the identification of genes and genomic regions responsible for virulence, as
well as in the diagnosis and treatment of infectious diseases. In this paper,
we present Sibelia, a tool for finding synteny blocks in multiple closely
related microbial genomes using iterative de Bruijn graphs. Unlike most other
tools, Sibelia can find synteny blocks that are repeated within genomes as well
as blocks shared by multiple genomes. It represents synteny blocks in a
hierarchy structure with multiple layers, each of which representing a
different granularity level. Sibelia has been designed to work efficiently with
a large number of microbial genomes; it finds synteny blocks in 31 S. aureus
genomes within 31 minutes and in 59 E.coli genomes within 107 minutes on a
standard desktop. Sibelia software is distributed under the GNU GPL v2 license
and is available at: https://github.com/bioinf/Sibelia Sibelia's web-server is
available at: http://etool.me/software/sibelia


  Presented here is a simple method for cross-validated genome-wide association
studies (cvGWAS). Focusing on phenotype prediction, the method is able to
reveal a significant amount of missing heritability by properly selecting a
small number of loci with implicit predictive ability. The results provide new
insights into the missing heritability problem and the underlying genetic
architecture of complex traits.