  The history of southern Africa involved interactions between indigenous
hunter-gatherers and a range of populations that moved into the region. Here we
use genome-wide genetic data to show that there are at least two admixture
events in the history of Khoisan populations (southern African hunter-gatherers
and pastoralists who speak non-Bantu languages with click consonants). One
involved populations related to Niger-Congo-speaking African populations, and
the other introduced ancestry most closely related to west Eurasian (European
or Middle Eastern) populations. We date this latter admixture event to
approximately 900-1,800 years ago, and show that it had the largest demographic
impact in Khoisan populations that speak Khoe-Kwadi languages. A similar signal
of west Eurasian ancestry is present throughout eastern Africa. In particular,
we also find evidence for two admixture events in the history of Kenyan,
Tanzanian, and Ethiopian populations, the earlier of which involved populations
related to west Eurasians and which we date to approximately 2,700 - 3,300
years ago. We reconstruct the allele frequencies of the putative west Eurasian
population in eastern Africa, and show that this population is a good proxy for
the west Eurasian ancestry in southern Africa. The most parsimonious
explanation for these findings is that west Eurasian ancestry entered southern
Africa indirectly through eastern Africa.


  The de novo assembly of large, complex genomes is a significant challenge
with currently available DNA sequencing technology. While many de novo assembly
software packages are available, comparatively little attention has been paid
to assisting the user with the assembly. This paper addresses the practical
aspects of de novo assembly by introducing new ways to perform quality
assessment on a collection of DNA sequence reads. The software implementation
calculates per-base error rates, paired-end fragment size histograms and
coverage metrics in the absence of a reference genome. Additionally, the
software will estimate characteristics of the sequenced genome, such as repeat
content and heterozygosity, that are key determinants of assembly difficulty.
The software described is freely available and open source under the GNU Public
License.


  Mississippi River nutrient loads and water stratification on the
Louisiana-Texas shelf contribute to an annually recurring, short-lived hypoxic
bottom layer in areas of the northern Gulf of Mexico comprising less than 2% of
the total Gulf of Mexico bottom area. This paper observes that the NOAA and
LUMCON have now published errant predictions of possible record size areas of
temporary bottom water hypoxia ("dead zones") three times since 2005, in 2008,
2011, and 2013 and that the LUMCON predictions of the area of hypoxic bottom
water average 31% higher than the actual measured hypoxic areas from 2006 to
2014. These systematically high predictions depend on the assumption that the
susceptibility of the Gulf of Mexico to forming hypoxic areas in response to
nutrient loading has been relatively constant since 2001, though the
susceptibility has been occasionally adjusted upward in different models. It
has been previously suggested that tropical storms in a given year that occur
on the Louisiana-Texas shelf between the peak nutrient loading in spring and
formation of the hypoxic zone in summer can mitigate the size of the hypoxic
zone that year through mixing of stratified well oxygenated lighter and warmer
surface layers and oxygen depleted heavier and cooler bottom layers. This paper
suggests several reasons why the Louisiana-Texas shelf may be systematically
growing less susceptible to a given level of nutrient loading over time so that
predictions based on the measured area of temporary bottom water hypoxia prior
to 2006 tend to be too big in recent years.


  Based on the measurements of noise in gene expression performed during the
last decade, it has become customary to think of gene regulation in terms of a
two-state model, where the promoter of a gene can stochastically switch between
an ON and an OFF state. As experiments are becoming increasingly precise and
the deviations from the two-state model start to be observable, we ask about
the experimental signatures of complex multi-state promoters, as well as the
functional consequences of this additional complexity. In detail, we (i) extend
the calculations for noise in gene expression to promoters described by state
transition diagrams with multiple states, (ii) systematically compute the
experimentally accessible noise characteristics for these complex promoters,
and (iii) use information theory to evaluate the channel capacities of complex
promoter architectures and compare them to the baseline provided by the
two-state model. We find that adding internal states to the promoter
generically decreases channel capacity, except in certain cases, three of which
(cooperativity, dual-role regulation, promoter cycling) we analyze in detail.


  Intrahost simian immunodeficiency virus (SIV) evolution is marked by repeated
viral escape from cytotoxic T-lymphocyte (CTLs) response. Typically, the first
such CTL escape occurs in a matter of days, starting around the time of peak
viral load. Many authors have developed methods to quantify the strength of CTL
response by measuring the rate at which CTL escape occurs, but such methods
usually depend on sampling at two or more timepoints, while many datasets
capture the dynamics of the first CTL escape at only a single timepoint. Here,
we develop inference methods for CTL escape rates applicable to single
timepoint datasets. Through a model of early infection dynamics, we construct
confidence intervals for escape rates, but since early infection dynamics are
not completely understood, we also develop a one-sided confidence interval
serving as a lower bound for escape rates over a collection of early infection
models. We apply our methods to two SIV datasets, using our lower bounds and
existing methods to show that escape rates are relatively high during the
initial days of the first CTL escape and then drop to lower levels as the
escape proceeds. We also compare escape in the lymph nodes and the rectal
mucosa, showing that escape in the lymph nodes is initially faster, but as the
first escape proceeds, the rate of escape in the lymph nodes drops below the
rate seen in the rectal mucosa.


  Asynchronous replication of the genome has been associated with different
rates of point mutation and copy number variation (CNV) in human populations.
Here, we explored if the bias in the generation of CNV that is associated to
DNA replication timing might have conditioned the birth of new protein-coding
genes during evolution. We show that genes that were duplicated during primate
evolution are more commonly found among the human genes located in
late-replicating CNV regions. We traced the relationship between replication
timing and the evolutionary age of duplicated genes. Strikingly, we found that
there is a significant enrichment of evolutionary younger duplicates in late
replicating regions of the human and mouse genome. Indeed, the presence of
duplicates in late replicating regions gradually decreases as the evolutionary
time since duplication extends. Our results suggest that the accumulation of
recent duplications in late replicating CNV regions is an active process
influencing genome evolution.


  We study a stochastic spatial model of biological competition in which two
species have the same birth and death rates, but different diffusion constants.
In the absence of this difference, the model can be considered as an
off-lattice version of the Voter model and presents similar coarsening
properties. We show that even a relative difference in diffusivity on the order
of a few percent may lead to a strong bias in the coarsening process favoring
the more agile species. We theoretically quantify this selective advantage and
present analytical formulas for the average growth of the fastest species and
its fixation probability.


  Although there has been much interest in estimating divergence and admixture
from genomic data, it has proven difficult to distinguish gene flow after
divergence from alternative histories involving structure in the ancestral
population. The lack of a formal test to distinguish these scenarios has
sparked recent controversy about the possibility of interbreeding between
Neandertals and modern humans in Eurasia. We derive the probability of
mutational configurations in non-recombining sequence blocks under alternative
histories of divergence with admixture and ancestral structure. Dividing the
genome into short blocks makes it possible to compute maximum likelihood
estimates of parameters under both models. We apply this method to triplets of
human Neandertal genomes and quantify the relative support for models of
long-term population structure in the ancestral African popuation and admixture
from Neandertals into Eurasian populations after their expansion out of Africa.
Our analysis allows us -- for the first time -- to formally reject a history of
ancestral population structure and instead reveals strong support for admixture
from Neandertals into Eurasian populations at a higher rate (3.4%-7.9%) than
suggested previously.


  The statistical and dynamical properties of ions in the selectivity filter of
the KcsA ion channel are considered on the basis of molecular dynamics (MD)
simulations of the KcsA protein embedded in a lipid membrane surrounded by an
ionic solution. A new approach to the derivation of a Brownian dynamics (BD)
model of ion permeation through the filter is discussed, based on unbiased MD
simulations. It is shown that depending on additional assumptions, ion's
dynamics can be described either by under-damped Langevin equation with
constant damping and white noise or by Langevin equation with a fractional
memory kernel. A comparison of the potential of the mean force derived from
unbiased MD simulations with the potential produced by the umbrella sampling
method demonstrates significant differences in these potentials. The origin of
these differences is an open question that requires further clarifications.


  With the advent of Next-Generation (NG) sequencing, it has become possible to
sequence an entire genome quickly and inexpensively. However, in some
experiments one only needs to extract and assembly a portion of the sequence
reads, for example when performing transcriptome studies, sequencing
mitochondrial genomes, or characterizing exomes. With the raw DNA-library of a
complete genome it would appear to be a trivial problem to identify reads of
interest. But it is not always easy to incorporate well-known tools such as
BLAST, BLAT, Bowtie, and SOAP directly into a bioinformatics pipelines before
the assembly stage, either due to in- compatibility with the assembler's file
inputs, or because it is desirable to incorporate information that must be
extracted separately. For example, in order to incorporate flowgrams from a
Roche 454 sequencer into the Newbler assembler it is necessary to first extract
them from the original SFF files. We present SlopMap, a bioinformatics software
utility which allows rapid identification similar to provided target sequences
from either Roche 454 or Illumnia DNA library. With a simple and intuitive
command- line interface along with file output formats compatible with assembly
programs, SlopMap can be directly embedded in biological data processing
pipeline without any additional programming work. In addition, SlopMap
preserves flowgram information needed for Roche 454 assembler.


  Study sample sizes in human genetics are growing rapidly, and in due course
it will become routine to analyze samples with hundreds of thousands if not
millions of individuals. In addition to posing computational challenges, such
large sample sizes call for carefully re-examining the theoretical foundation
underlying commonly-used analytical tools. Here, we study the accuracy of the
coalescent, a central model for studying the ancestry of a sample of
individuals. The coalescent arises as a limit of a large class of random mating
models and it is an accurate approximation to the original model provided that
the population size is sufficiently larger than the sample size. We develop a
method for performing exact computation in the discrete-time Wright-Fisher
(DTWF) model and compare several key genealogical quantities of interest with
the coalescent predictions. For realistic demographic scenarios, we find that
there are a significant number of multiple- and simultaneous-merger events
under the DTWF model, which are absent in the coalescent by construction.
Furthermore, for large sample sizes, there are noticeable differences in the
expected number of rare variants between the coalescent and the DTWF model. To
balance the tradeoff between accuracy and computational efficiency, we propose
a hybrid algorithm that utilizes the DTWF model for the recent past and the
coalescent for the more distant past. Our results demonstrate that the hybrid
method with only a handful of generations of the DTWF model leads to a
frequency spectrum that is quite close to the prediction of the full DTWF
model.


  The amount of sequence data available today highly facilitates the access to
genes from many gene families. Universal primers amplifying the desired genes
over a range of species are readily obtained by aligning conserved gene
regions, and laborious gene isolation procedures can often be replaced by
quicker PCR-based approaches. However, in case of multigene families, PCR-based
approaches bear the risk of incomplete isolation of family members. This
problem is most prominent in gene families with highly variable and thus
unpredictable number of gene copies among species, such as in the major
histocompatibility complex (MHC). In the present study we (i) report new
primers for the isolation of the MHC class IIB (MHCIIB) gene family in birds,
and (ii) share our experience with isolating MHCIIB genes from an unprecedented
number of avian species from all over the avian phylogeny. We report important
and usually underappreciated problems encountered during PCR-based multigene
family isolation, and provide a collection of measures that may help to
significantly improve the chance of successfully isolating complete multigene
families using PCR-based approaches.


  Natural selection drives populations towards higher fitness, but crossing
fitness valleys or plateaus may facilitate progress up a rugged fitness
landscape involving epistasis. We investigate quantitatively the effect of
subdividing an asexual population on the time it takes to cross a fitness
valley or plateau. We focus on a generic and minimal model that includes only
population subdivision into equivalent demes connected by global migration, and
does not require significant size changes of the demes, environmental
heterogeneity or specific geographic structure. We determine the optimal
speedup of valley or plateau crossing that can be gained by subdivision, if the
process is driven by the deme that crosses fastest. We show that isolated demes
have to be in the sequential fixation regime for subdivision to significantly
accelerate crossing. Using Markov chain theory, we obtain analytical
expressions for the conditions under which optimal speedup is achieved: valley
or plateau crossing by the subdivided population is then as fast as that of its
fastest deme. We verify our analytical predictions through stochastic
simulations. We demonstrate that subdivision can substantially accelerate the
crossing of fitness valleys and plateaus in a wide range of parameters
extending beyond the optimal window. We study the effect of varying the degree
of subdivision of a population, and investigate the trade-off between the
magnitude of the optimal speedup and the width of the parameter range over
which it occurs. Our results also hold for weakly beneficial intermediate
mutations. We extend our work to the case of a population connected by
migration to one or several smaller islands. Our results demonstrate that
subdivision with migration alone can significantly accelerate the crossing of
fitness valleys and plateaus, and shed light onto the quantitative conditions
necessary for this to occur.


  Genetic regulatory networks are defined by their topology and by a multitude
of continuously adjustable parameters. Here we present a class of simple models
within which the relative importance of topology vs. interaction strengths
becomes a well-posed problem. We find that complexity - the ability of the
network to adopt multiple stable states - is dominated by the adjustable
parameters. We comment on the implications for real networks and their
evolution.


  Most non-synonymous mutations are thought to be deleterious because of their
effect on protein sequence. These polymorphisms are expected to be removed or
kept at low frequency by the action of natural selection, and rare deleterious
variants have been implicated as a possible explanation for the "missing
heritability" seen in many studies of complex traits. Nonetheless, the effect
of positive selection on linked sites or drift in small or inbred populations
may also impact the evolution of deleterious alleles. Here, we made use of
genome-wide genotyping data to characterize deleterious variants in a large
panel of maize inbred lines. We show that, in spite of small effective
population sizes and inbreeding, most putatively deleterious SNPs are indeed at
low frequencies within individual genetic groups. We find that genes showing
associations with a number of complex traits are enriched for deleterious
variants. Together these data are consistent with the dominance model of
heterosis, in which complementation of numerous low frequency, weak deleterious
variants contribute to hybrid vigor.


  The survival and proliferation of cells and organisms require a highly
coordinated allocation of cellular resources to ensure the efficient synthesis
of cellular components. In particular, the total enzymatic capacity for
cellular metabolism is limited by finite resources that are shared between all
enzymes, such as cytosolic space, energy expenditure for amino-acid synthesis,
or micro-nutrients. While extensive work has been done to study constrained
optimization problems based only on stoichiometric information, mathematical
results that characterize the optimal flux in kinetic metabolic networks are
still scarce. Here, we study constrained enzyme allocation problems with
general kinetics, using the theory of oriented matroids. We give a rigorous
proof for the fact that optimal solutions of the non-linear optimization
problem are elementary flux modes. This finding has significant consequences
for our understanding of optimality in metabolic networks as well as for the
identification of metabolic switches and the computation of optimal flux
distributions in kinetic metabolic networks.


  We study the average asymptotic growth rate of cells in randomly fluctuating
environments, with multiple viable phenotypes per environment. We show that any
information processing strategy has an asymptotic growth rate, which is the sum
of: (i) the maximal growth rate at the worst possible distribution of
environments, (ii) relative information between the actual distribution of
environments to the worst one, and (iii) information utilization rate, which is
the information rate of the sensory devices minus the "information dissipation
rate", the amount of information not utilized by the cell for growth. In
non-stationary environments, we find that the optimal phenotypic switching
times equally partition the information dissipation rate between consecutive
switching intervals.


  The Maximum Entropy Theory of Ecology (METE) predicts a universal
species-area relationship (SAR) that can be fully characterized using only the
total abundance (N) and species richness (S) at a single spatial scale. This
theory has shown promise for characterizing scale dependence in the SAR.
However, there are currently four different approaches to applying METE to
predict the SAR and it is unclear which approach should be used due to a lack
of empirical evaluation. Specifically, METE can be applied recursively or a
non-recursively and can use either a theoretical or observed species-abundance
distribution (SAD). We compared the four different combinations of approaches
using empirical data from 16 datasets containing over 1000 species and 300,000
individual trees and herbs. In general, METE accurately downscaled the SAR
(R^2> 0.94), but the recursive approach consistently under-predicted richness,
and METEs accuracy did not depend strongly on using the observed or predicted
SAD. This suggests that best approach to scaling diversity using METE is to use
a combination of non-recursive scaling and the theoretical abundance
distribution, which allows predictions to be made across a broad range of
spatial scales with only knowledge of the species richness and total abundance
at a single scale.


  Social animals have to make collective decisions on a daily basis. In most
instances, these decisions are taken by consensus, when the group does what the
majority of individuals want. Individuals have to base these decisions on the
information they perceive from their socioecological landscape. The perception
mechanisms they use can influence the cost of collective decisions. Here I show
that when group-living individuals perceive their environment concurrently for
the same decisions, a quantum-like perception entanglement process can confer
less costly collective decisions than when individuals collect their
information independently. This highlights a mechanism that can help explain
what may seem to be irrational group-living behavior and opens avenues to
develop empirical tests for quantum decision theory.


  The commentators have brought a wealth of new perspectives to the question of
how culture evolves. Each of their diverse disciplines--ranging from psychology
to biology to anthropology to economics to engineering--has a valuable
contribution to make to our understanding of this complex, multifaceted topic.
Though the vast majority of their comments were supportive of my approach, it
is natural that a reply such as this focus on points where my views differ from
that of the commentators. ... I conclude by saying that I am grateful to the
commentators for their diverse perspectives and insights, their overall support
for the project, and provocative ideas for where to go from here. Clearly there
are many fascinating avenues to explore as we move forward on our quest to
understand how culture evolves.


  The spatial structure of the cell is highly organized at all levels: from
small complexes and assemblies, to local nano- and micro-clusters, to global,
micrometer scales across and between cells. We suggest that this multiscale
spatial cell organization also organizes signaling and coordinates cellular
behavior. We propose a new view of the spatial structure of cell signaling
systems. This new view describes cell signaling in terms of dynamic allosteric
interactions within and among distinct, spatially organized transient clusters.
The clusters vary over time and space and are on length scales from nanometers
to micrometers. When considered across these length-scales, primary factors in
the spatial organization are cell membrane domains and the actin cytoskeleton,
both also highly dynamic. A key challenge is to understand the interplay across
these multiple scales, link it to the physicochemical basis of the
conformational behavior of single molecules, and ultimately relate it to
cellular function. Overall, our premise is that at these scales, cell signaling
should be thought of not primarily as a sequence of diffusion-controlled
molecular collisions, but instead transient, allostery-driven cluster
re-forming interactions.


  The reproduction number of deterministic models is an essential quantity to
predict whether an epidemic will spread or die out. Thresholds for disease
extinction contribute crucial knowledge on disease control, elimination, and
mitigation of infectious diseases. Relationships between the basic reproduction
numbers of two network-based ordinary differential equation vector-host models,
and extinction thresholds of corresponding continuous-time Markov chain models
are derived under some assumptions. Numerical simulation results for malaria
and Rift Valley fever transmission on heterogeneous networks are in agreement
with analytical results without any assumptions, reinforcing the relationships
may always exist and proposing a mathematical problem of proving their
existences in general. Moreover, numerical simulations show that the
reproduction number is not monotonically increasing or decreasing with the
extinction threshold. Key parameters in predicting uncertainty of extinction
thresholds are identified using Latin Hypercube Sampling/Partial Rank
Correlation Coefficient. Consistent trends of extinction probability observed
through numerical simulations provide novel insights into mitigation strategies
to increase the disease extinction probability. Research findings may improve
understandings of thresholds for disease persistence in order to control
vector-borne diseases.


  What makes economic and ecological networks so unlike other highly skewed
networks in their tendency toward turbulence and collapse? Here, we explore the
consequences of a defining feature of these networks: their nodes are tied
together by flow. We show that flow networks tend to the power law degree
distribution (PLDD) due to a self-reinforcing process involving position within
the global network structure, and thus present the first random graph model for
PLDDs that does not depend on a rich-get-richer function of nodal degree. We
also show that in contrast to non-flow networks, PLDD flow networks are
dramatically more vulnerable to catastrophic failure than non-PLDD flow
networks, a finding with potential explanatory power in our age of resource-
and financial-interdependence and turbulence.


  The Maximum Entropy Theory of Ecology (METE) is a unified theory of
biodiversity that predicts a large number of macroecological patterns using
only information on the species richness, total abundance, and total metabolic
rate of the community. We evaluated four major predictions of METE
simultaneously at an unprecedented scale using data from 60 globally
distributed forest communities including over 300,000 individuals and nearly
2000 species. METE successfully captured 96% and 89% of the variation in the
species abundance distribution and the individual size distribution, but
performed poorly when characterizing the size-density relationship and
intraspecific distribution of individual size. Specifically, METE predicted a
negative correlation between size and species abundance, which is weak in
natural communities. By evaluating multiple predictions with large quantities
of data, our study not only identifies a mismatch between abundance and body
size in METE, but also demonstrates the importance of conducting strong tests
of ecological theories.


  Very low levels of genetic diversity have been reported in vertebrates with
large genomes, notably salamanders and lungfish [1-3]. Interpreting differences
in heterozygosity, which reflects genetic diversity in a population, is
complicated because levels of heterozygosity vary widely between conspecific
populations, and correlate with many different physiological and demographic
variables such as body size and effective population size. Here we return to
the question of genetic variability in salamanders and report on the
relationship between evolutionary rates and genome sizes in five different
salamander families. We found that rates of evolution are exceptionally low in
salamanders as a group. Evolutionary rates are as low as those reported for
cartilaginous fish, which have the slowest rates recorded so far in vertebrates
[4]. We also found that, independent of life history, salamanders with the
smallest genomes (14 pg) are evolving at rates two to three times faster than
salamanders with the largest genomes (>50 pg). After accounting for
evolutionary duration, we conclude that more recently evolved species have
correspondingly smaller genomes compared to older taxa and concomitantly higher
rates of mutation and evolution.


  Due to the increasing use of nanometric cerium oxide in applications,
concerns about the toxicity of these particles have been raised and have
resulted in a large number of investigations. We report here on the
interactions between 7 nm anionically charged cerium oxide particles and living
mammalian cells. By a modification of the particle coating including
low-molecular weight ligands and polymers, two generic behaviors are compared:
particles coated with citrate ions that precipitate in biofluids and particles
coated with poly(acrylic acid) that are stable and remain nanometric. We find
that nanoceria covered with both coating agents are taken up by mouse
fibroblasts and localized into membrane-bound compartments. However, flow
cytometry and electron microscopy reveal that as a result of their
precipitation, citrate-coated particles interact more strongly with cells. At
cerium concentration above 1 mM, only citrate-coated nanoceria (and not
particles coated with poly(acrylic acid)) display toxicity and moderate
genotoxicity. The results demonstrate that the control of the surface chemistry
of the particles and its ability to prevent aggregation can affect the toxicity
of nanomaterials.


  We study the stationary state of a population evolving under the action of
random genetic drift, selection and recombination in which both deleterious and
reverse beneficial mutations can occur. We find that the equilibrium fraction
of deleterious mutations decreases as the population size is increased. We
calculate exactly the steady state frequency in a nonrecombining population
when population size is infinite and for a neutral finite population, and
obtain bounds on the fraction of deleterious mutations. We also find that for
small and very large populations, the number of deleterious mutations depends
weakly on recombination, but for moderately large populations, recombination
alleviates the effect of deleterious mutations. An analytical argument shows
that recombination decreases disadvantageous mutations appreciably when
beneficial mutations are rare as is the case in adapting microbial populations,
whereas it has a moderate effect on codon bias where the mutation rates between
the preferred and unpreferred codons are comparable.


  A five-species predator-prey model is studied on a square lattice where each
species has two prey and two predators on the analogy to the
Rock-Paper-Scissors-Lizard-Spock game. The evolution of the spatial
distribution of species is governed by site exchange and invasion between the
neighboring predator-prey pairs, where the cyclic symmetry can be characterized
by two different invasion rates. The mean-field analysis has indicated periodic
oscillations in the species densities with a frequency becoming zero for a
specific ratio of invasion rates. When varying the ratio of invasion rates, the
appearance of this zero-eigenvalue mode is accompanied by neutrality between
the species associations. Monte Carlo simulations of the spatial system reveal
diverging fluctuations at a specific invasion rate, which can be related to the
vanishing dominance between all pairs of species associations.


  Metabolic substrates, such as oxygen and glucose, are rapidly delivered to
the cell through filtration across microvessels walls. Modelling this important
process is complicated by the coupling between flow and transport equations,
which are linked through the osmotic pressure induced by the colloidal plasma
proteins. Microvessel wall is a composite media with the internal glycocalyx
layer exerting a remarkable sieving effect on macromolecules, with respect to
the external layer composed by endothelial cells. The physiological structure
of microvessel is represented as the superimposition of two membranes with
different properties; the inner membrane represents the glycocalyx, while the
outer membrane represents the surrounding endothelial cells. Application of
mass conservation principle and thermodynamic considerations lead to a model
composed by two coupled second-order partial differential equations in the
hydrostatic and osmotic pressures, one expressing volumetric mass conservation
and the other, which is non-linear, expressing macromolecules mass
conservation. Despite system complexity, the assumption that the properties of
the layers are piece-wise constant allows us to obtain analytical solutions for
the two pressures. This solution is in agreement with experimental
observations, which contrary to common belief, show that flow reversal cannot
occur in steady state unless luminal hydrostatic pressure drops below
physiologically plausible values. The observed variations of volumetric and
solute mass flux in case of a significant reduction of luminal hydrostatic
pressure are in qualitative agreement with observed variations during
experiments reported in the literature. On the other hand, homogenising
microvessel wall into a single-layer membrane with equivalent properties leads
to a different distribution of pressure across microvessel wall, not consistent
with observations.


  Forward-time models of diversification (i.e., speciation and extinction)
produce phylogenetic trees that grow "vertically" as time goes by. Pruning the
extinct lineages out of such trees leads to natural models for reconstructed
trees (i.e., phylogenies of extant species). Alternatively, reconstructed trees
can be modelled by coalescent point processes (CPP), where trees grow
"horizontally" by the sequential addition of vertical edges. Each new edge
starts at some random speciation time and ends at the present time; speciation
times are drawn from the same distribution independently. CPP lead to extremely
fast computation of tree likelihoods and simulation of reconstructed trees.
Their topology always follows the uniform distribution on ranked tree shapes
(URT). We characterize which forward-time models lead to URT reconstructed
trees and among these, which lead to CPP reconstructed trees. We show that for
any "asymmetric" diversification model in which speciation rates only depend on
time and extinction rates only depend on time and on a non-heritable trait
(e.g., age), the reconstructed tree is CPP, even if extant species are
incompletely sampled. If rates additionally depend on the number of species,
the reconstructed tree is (only) URT (but not CPP). We characterize the common
distribution of speciation times in the CPP description, and discuss incomplete
species sampling as well as three special model cases in detail: 1) extinction
rate does not depend on a trait; 2) rates do not depend on time; 3) mass
extinctions may happen additionally at certain points in the past.


  Cytotoxic T-lymphocytes (CTLs) recognize viral protein fragments displayed by
major histocompatibility complex (MHC) molecules on the surface of virally
infected cells and generate an anti-viral response that can kill the infected
cells. Virus variants whose protein fragments are not efficiently presented on
infected cells or whose fragments are presented but not recognized by CTLs
therefore have a competitive advantage and spread rapidly through the
population. We present a method that allows a more robust estimation of these
escape rates from serially sampled sequence data. The proposed method accounts
for competition between multiple escapes by explicitly modeling the
accumulation of escape mutations and the stochastic effects of rare multiple
mutants. Applying our method to serially sampled HIV sequence data, we estimate
rates of HIV escape that are substantially larger than those previously
reported. The method can be extended to complex escapes that require
compensatory mutations. We expect our method to be applicable in other contexts
such as cancer evolution where time series data is also available.


  Multiple Sclerosis (MS) is a disorder that usually appears in adults in their
thirties. It has a prevalence that ranges between 2 and 150 per 100 000.
Epidemiological studies of MS have provided hints on possible causes for the
disease ranging from genetic, environmental and infectious factors to other
factors of vascular origin. Despite the tremendous effort spent in the last few
years, none of the hypotheses formulated so far has gained wide acceptance and
the causes of the disease remain unknown. From a clinical point of view, a high
correlation has been recently observed between MS and Chronic Cerebro-Spinal
Venous Insufficiency (CCSVI) in a statistically significant number of patients.
In this pathological situation CCSVI may induce alterations of blood pressure
in brain microvessels, thereby perturbing the exchange of small hydrophilic
molecules between the blood and the external cells. In the presence of large
pressure alterations it cannot be excluded also the leakage of macromolecules
that otherwise would not cross the vessel wall. All these disorders may trigger
immune defenses with the destruction of myelin as a side effect. In the present
work we investigate the role of perturbed blood pressure in brain microvessels
as driving force for an altered exchange of small hydrophilic solutes and
leakage of macromolecules into the interstitial fluid. With a simplified, yet
realistic, model we obtain closed-form steady-state solutions for fluid flow
and solute transport across the microvessel wall. Finally, we use these results
(i) to interpret experimental data available in the literature and (ii) to
carry out a preliminary analysis of the disorder in the exchange processes
triggered by an increase of blood pressure, thereby relating our preliminary
results to the hypothesised vascular connection to MS.


  Most genome assemblers construct point estimates, choosing a genome sequence
from among many alternative hypotheses that are supported by the data. We
present a Markov Chain Monte Carlo approach to sequence assembly that instead
generates distributions of assembly hypotheses with posterior probabilities,
providing an explicit statistical framework for evaluating alternative
hypotheses and assessing assembly uncertainty. We implement this approach in a
prototype assembler and illustrate its application to the bacteriophage
PhiX174.


  We introduce and analyze a within-host dynamical model of the coevolution
between rapidly mutating pathogens and the adaptive immune response. Pathogen
mutation and a homeostatic constraint on lymphocytes both play a role in
allowing the development of chronic infection, rather than quick pathogen
clearance. The dynamics of these chronic infections display emergent structure,
including branching patterns corresponding to asexual pathogen speciation,
which is fundamentally driven by the coevolutionary interaction. Over time,
continued branching creates an increasingly fragile immune system, and leads to
the eventual catastrophic loss of immune control.


  Living cells deploy many resources to sense their environments, including
receptors, downstream signaling molecules, time and fuel. However, it is not
known which resources fundamentally limit the precision of sensing, like weak
links in a chain, and which can compensate each other, leading to trade-offs
between them. We show by modeling that in equilibrium systems the precision is
limited by the number of receptors; the downstream network can never increase
precision. This limit arises from a trade-off between the removal of extrinsic
noise in the receptor and intrinsic noise in the downstream network.
Non-equilibrium systems can lift this trade-off by storing the receptor state
over time in chemical modification states of downstream molecules. As we
quantify for a push-pull network, this requires i) time and receptors; ii)
downstream molecules; iii) energy (fuel turnover) to drive modification. These
three resource classes cannot compensate each other, and it is the limiting
class which sets the fundamental sensing limit. Within each class, trade-offs
are possible. Energy allows a power-speed trade-off, while time can be traded
against receptors.


  Measuring natural selection on genomic elements involved in the
cis-regulation of gene expression -- such as transcriptional enhancers and
promoters -- is critical for understanding the evolution of genomes, yet it
remains a major challenge. Many studies have attempted to detect positive or
negative selection in these noncoding elements by searching for those with the
fastest or slowest rates of evolution, but this can be problematic. Here we
introduce a new approach to this issue, and demonstrate its utility on three
mammalian transcriptional enhancers. Using results from saturation mutagenesis
studies of these enhancers, we classified all possible point mutations as
up-regulating, down-regulating, or silent, and determined which of these
mutations have occurred on each branch of a phylogeny. Applying a framework
analogous to Ka/Ks in protein-coding genes, we measured the strength of
selection on up-regulating and down-regulating mutations, in specific branches
as well as entire phylogenies. We discovered distinct modes of selection acting
on different enhancers: while all three have experienced negative selection
against down-regulating mutations, the selection pressures on up-regulating
mutations vary. In one case we detected positive selection for up-regulation,
while the other two had no detectable selection on up-regulating mutations. Our
methodology is applicable to the growing number of saturation mutagenesis data
sets, and provides a detailed picture of the mode and strength of natural
selection acting on cis-regulatory elements.


  Despite recent advances in our ability to detect adaptive evolution involving
the cis-regulation of gene expression, our knowledge of the molecular
mechanisms underlying these adaptations has lagged far behind. Across all model
organisms the causal mutations have been discovered for only a handful of gene
expression adaptations, and even for these, mechanistic details (e.g. the
trans-regulatory factors involved) have not been determined. We previously
reported a polygenic gene expression adaptation involving down-regulation of
the ergosterol biosynthesis pathway in the budding yeast Saccharomyces
cerevisiae. Here we investigate the molecular mechanism of a cis-acting
mutation affecting a member of this pathway, ERG28. We show that the causal
mutation is a two-base deletion in the promoter of ERG28 that strongly reduces
the binding of two transcription factors, Sok2 and Mot3, thus abolishing their
regulation of ERG28. This down-regulation increases resistance to a widely used
antifungal drug targeting ergosterol, similar to mutations disrupting this
pathway in clinical yeast isolates. The identification of the causal genetic
variant revealed that the selection likely occurred after the deletion was
already present at high frequency in the population, rather than when it was a
new mutation. These results provide a detailed view of the molecular mechanism
of a cis-regulatory adaptation, and underscore the importance of this view to
our understanding of evolution at the molecular level.


  We study the spatial evolutionary dynamics of solid tumors as they obtain
additional driver mutations. We start with a cancer clone that expands
uniformly in three dimensions giving rise to a spherical shape. We assume that
cell division occurs on the surface of the growing tumor. Each cell division
has a chance to give rise to a mutation that activates an additional driver
gene. The resulting clone has an enhanced growth rate, which generates a local
ensemble of faster growing cells, thereby distorting the spherical shape of the
tumor. We derive analytic formulas for the geometric boundary that separates
the original cancer clone from the new mutant as well as the expanding frontier
of the new mutant. The total number of original cancer cells converges to a
constant as time goes to infinity, because this clone becomes enveloped by
mutants. We derive formulas for the abundance and diversity of additional
driver mutations as function of time. Our model is semi-deterministic: the
spatial growth of the various cancer clones follows deterministic equations,
but the arrival of a new mutant is a stochastic event.


  Upon infection of their bacterial hosts temperate phages must chose between
lysogenic and lytic developmental strategies. Here we apply the game-theoretic
bet-hedging strategy introduced by Kelly to derive the optimal lysogenic
fraction of the total population of phages as a function of frequency and
intensity of environmental downturns affecting the lytic subpopulation.
"Well-temperate" phage from our title is characterized by the best long-term
population growth rate. We show that it is realized when the lysogenization
frequency is approximately equal to the probability of lytic population
collapse. We further predict the existence of sharp boundaries in system's
environmental, ecological, and biophysical parameters separating the regions
where this temperate strategy is optimal from those dominated by purely
virulent or} dormant (purely lysogenic) strategies. We show that the virulent
strategy works best for phages with large diversity of hosts, and access to
multiple independent environments reachable by diffusion. Conversely,
progressively more temperate or even dormant strategies are favored in the
environments, that are subject to frequent and severe temporal downturns.


  Regular vegetation patterns in semiarid ecosystems are believed to arise from
the interplay between long-range competition and facilitation processes acting
at smaller distances. We show that, under rather general conditions, long-range
competition alone may be enough to shape these patterns. To this end we propose
a simple, general model for the dynamics of vegetation, which includes only
long-range competition between plants. Competition is introduced through a
nonlocal term, where the kernel function quantifies the intensity of the
interaction. We recover the full spectrum of spatial structures typical of
vegetation models that also account for facilitation in addition to
competition.


  We developed an online database, ExprTargetDB, of human miRNA targets
predicted by an approach that integrates gene expression profiling into a
broader framework involving important features of miRNA target site
predictions.


  Protein translation is one of the most important processes in cell life but,
despite being well understood biochemically, the implications of its intrinsic
stochastic nature have not been fully elucidated. In this paper we develop a
microscopic and stochastic model which describes a crucial step in protein
translation, namely the binding of the tRNA to the ribosome. Our model
explicitly takes into consideration tRNA recharging dynamics, spatial
inhomogeneity and stochastic fluctuations in the number of charged tRNAs around
the ribosome. By analyzing this non-equilibrium system we are able to derive
the statistical distribution of the times needed by the tRNAs to bind to the
ribosome, and to show that it deviates from an exponential due to the coupling
between the fluctuations of charged and uncharged populations of tRNA.


  We pointed out that a substantial number of CpG probes on the Illumina 450K
array could be mapped to multiple loci across the human genome. These CpGs need
to be considered when interpreting results using this platform.


  Do cancer cells escape their confinement of their original habitat in the
primary tumor or are they forced out by ecological changes in their home niche?
Describing metastasis in terms of a simple one-way migration of cells from the
primary to target organs is an insufficient concept to cover the nuances of
cancer spread. A diaspora is the scattering of people away from an established
homeland. To date, diaspora has been a uniquely human term utilized by social
scientists, however, the application of the diaspora concept to metastasis may
yield new biological insights as well as therapeutic paradigms. The diaspora
paradigm takes into account and models several variables: the quality of the
primary tumor microenvironment, the fitness of individual cancer cell migrants
as well as migrant populations, the rate of bidirectional migration of cancer
and host cells between cancer sites, and the quality of the target
microenvironments to establish metastatic sites. Ecological scientific
principles can be applied to the cancer diaspora to develop new therapeutic
strategies. For example, ecological traps, habitats that lead to the extinction
of a species, can be developed to attract cancer cells to a place where they
can be better exposed to treatments or to cells of the immune system for
improved antigen presentation. Merging the social science concept of diaspora
with ecological and population sciences concepts can inform the cancer field to
understand the biology of tumorigenesis and metastasis and inspire new ideas
for therapy.


  Besides the structure of interactions within networks, also the interactions
between networks are of the outmost importance. We therefore study the outcome
of the public goods game on two interdependent networks that are connected by
means of a utility function, which determines how payoffs on both networks
jointly influence the success of players in each individual network. We show
that an unbiased coupling allows the spontaneous emergence of interdependent
network reciprocity, which is capable to maintain healthy levels of public
cooperation even in extremely adverse conditions. The mechanism, however,
requires simultaneous formation of correlated cooperator clusters on both
networks. If this does not emerge or if the coordination process is disturbed,
network reciprocity fails, resulting in the total collapse of cooperation.
Network interdependence can thus be exploited effectively to promote
cooperation past the limits imposed by isolated networks, but only if the
coordination between the interdependent networks is not disturbed.


  Motivation. Protein contact map describes the pairwise spatial and functional
relationship of residues in a protein and contains key information for protein
3D structure prediction. Although studied extensively, it remains very
challenging to predict contact map using only sequence information. Most
existing methods predict the contact map matrix element-by-element, ignoring
correlation among contacts and physical feasibility of the whole contact map. A
couple of recent methods predict contact map based upon residue co-evolution,
taking into consideration contact correlation and enforcing a sparsity
restraint, but these methods require a very large number of sequence homologs
for the protein under consideration and the resultant contact map may be still
physically unfavorable.
  Results. This paper presents a novel method PhyCMAP for contact map
prediction, integrating both evolutionary and physical restraints by machine
learning and integer linear programming (ILP). The evolutionary restraints
include sequence profile, residue co-evolution and context-specific statistical
potential. The physical restraints specify more concrete relationship among
contacts than the sparsity restraint. As such, our method greatly reduces the
solution space of the contact map matrix and thus, significantly improves
prediction accuracy. Experimental results confirm that PhyCMAP outperforms
currently popular methods no matter how many sequence homologs are available
for the protein under consideration. PhyCMAP can predict contacts within
minutes after PSIBLAST search for sequence homologs is done, much faster than
the two recent methods PSICOV and EvFold.
  See http://raptorx.uchicago.edu for the web server.


  Linkage disequilibrium (LD) is often summarized using the "LD curve," which
relates the LD between pairs of sites to the distance that separates them along
the chromosome. This paper shows how the LD curve responds to changes in
population size. An expansion of population size generates an LD curve that
declines steeply, especially if that expansion has followed a bottleneck. A
reduction in size generates an LD curve that is high but relatively flat. In
European data, the curve is steep, suggesting a history of population
expansion.
  These conclusions emerge from the study of $\sigma_d^2$, a measure of LD that
has never played a central role. It has been seen merely as an approximation to
another measure, $r^2$. Yet $\sigma_d^2$ has different dynamical behavior and
provides deeper time depth. Furthermore, it is easily estimated from data and
can be predicted from population history using a fast, deterministic algorithm.


  Eukaryotic DNA replication follows a specific temporal program, with some
genomic regions consistently replicating earlier than others, yet what
determines this program is largely unknown. Highly transcribed regions have
been observed to replicate in early S-phase in all plant and animal species
studied to date, but this relationship is thought to be absent from both
budding yeast and fission yeast. No association between cell-cycle regulated
transcription and replication timing has been reported for any species. Here I
show that in budding yeast, fission yeast, and human, the genes most highly
transcribed during S-phase replicate early, whereas those repressed in S-phase
replicate late. Transcription during other cell-cycle phases shows either the
opposite correlation with replication timing, or no relation. The relationship
is strongest near late-firing origins of replication, which is not consistent
with a previously proposed model -- that replication timing may affect
transcription -- and instead suggests a potential mechanism involving the
recruitment of limiting replication initiation factors during S-phase. These
results suggest that S-phase transcription may be an important determinant of
DNA replication timing across eukaryotes, which may explain the
well-established association between transcription and replication timing.


  The mechanical stretching of single poly-proteins is an emerging tool for the
study of protein (un)folding, chemical catalysis and polymer physics at the
single molecule level. The observed processes i.e unfolding or reduction
events, are typically considered to be stochastic and by its nature are
susceptible to be censored by the finite duration of the experiment. Here we
develop a formal analytical and experimental description on the number of
observed events under various conditions of practical interest. We provide a
rule of thumb to define the experiment protocol duration. Finally we provide a
methodology to accurately estimate the number of stretched molecules based on
the number of observed unfolding events. Using this analysis on experimental
data we conclude for the first time that poly-ubiquitin binds at a random
position both to the substrate and to the pulling probe and that observing all
the existing modules is the less likely event.


  Lung transplantation remains the only viable treatment option for the
majority of patients with advanced lung diseases. However, 5-year
post-transplant survival rates remain low primarily secondary to chronic
rejection. Novel insights from global gene expression profiles may provide
molecular phenotypes and therapeutic targets to improve outcomes after lung
transplantation. We showed the presence of a significant number of dysregulated
genes, particularly those genes involved in pathways and biological processes
such as immune response and defense, in the PBMCs derived from a cohort of
patients after lung transplantation. The contribution of miRNAs in regulating
these differential genes was also demonstrated.


  The recurrent infectious diseases and their increasing impact on the society
has promoted the study of strategies to slow down the epidemic spreading. In
this review we outline the applications of percolation theory to describe
strategies against epidemic spreading on complex networks. We give a general
outlook of the relation between link percolation and the
susceptible-infected-recovered model, and introduce the node void percolation
process to describe the dilution of the network composed by healthy individual,
$i.e$, the network that sustain the functionality of a society. Then, we survey
two strategies: the quenched disorder strategy where an heterogeneous
distribution of contact intensities is induced in society, and the intermittent
social distancing strategy where health individuals are persuaded to avoid
contact with their neighbors for intermittent periods of time. Using
percolation tools, we show that both strategies may halt the epidemic
spreading. Finally, we discuss the role of the transmissibility, $i.e$, the
effective probability to transmit a disease, on the performance of the
strategies to slow down the epidemic spreading.


  Background: With the fast development of next generation sequencing
technologies, increasing numbers of genomes are being de novo sequenced and
assembled. However, most are in fragmental and incomplete draft status, and
thus it is often difficult to know the accurate genome size and repeat content.
Furthermore, many genomes are highly repetitive or heterozygous, posing
problems to current assemblers utilizing short reads. Therefore, it is
necessary to develop efficient assembly-independent methods for accurate
estimation of these genomic characteristics. Results: Here we present a
framework for modeling the distribution of k-mer frequency from sequencing data
and estimating the genomic characteristics such as genome size, repeat
structure and heterozygous rate. By introducing novel techniques of k-mer
individuals, float precision estimation, and proper treatment of sequencing
error and coverage bias, the estimation accuracy of our method is significantly
improved over existing methods. We also studied how the various genomic and
sequencing characteristics affect the estimation accuracy using simulated
sequencing data, and discussed the limitations on applying our method to real
sequencing data. Conclusion: Based on this research, we show that the k-mer
frequency analysis can be used as a general and assembly-independent method for
estimating genomic characteristics, which can improve our understanding of a
species genome, help design the sequencing strategy of genome projects, and
guide the development of assembly algorithms. The programs developed in this
research are written using C/C++, and freely accessible at Github URL
(https://github.com/fanagislab/GCE) or BGI ftp (
ftp://ftp.genomics.org.cn/pub/gce).


  Many organs of higher organisms are heavily branched structures and arise by
an at first sight similar process of branching morphogenesis. Yet the
regulatory components and local interactions that have been identified differ
greatly in these organs. It is an open question whether the regulatory
processes work according to a common principle and in how far physical and
geometric constraints determine the branching process. Here we review the known
regulatory factors and physical constraints in lung, kidney, pancreas,
prostate, mammary and salivary gland branching morphogenesis, and describe the
models that have been formulated to analyse their impacts.


  Background: Recent assays for individual-specific genome-wide DNA methylation
profiles have enabled epigenome-wide association studies to identify specific
CpG sites associated with a phenotype. Computational prediction of CpG
site-specific methylation levels is important, but current approaches tackle
average methylation within a genomic locus and are often limited to specific
genomic regions. Results: We characterize genome-wide DNA methylation patterns,
and show that correlation among CpG sites decays rapidly, making predictions
solely based on neighboring sites challenging. We built a random forest
classifier to predict CpG site methylation levels using as features neighboring
CpG site methylation levels and genomic distance, and co-localization with
coding regions, CGIs, and regulatory elements from the ENCODE project, among
others. Our approach achieves 91% -- 94% prediction accuracy of genome-wide
methylation levels at single CpG site precision. The accuracy increases to 98%
when restricted to CpG sites within CGIs. Our classifier outperforms
state-of-the-art methylation classifiers and identifies features that
contribute to prediction accuracy: neighboring CpG site methylation status, CpG
island status, co-localized DNase I hypersensitive sites, and specific
transcription factor binding sites were found to be most predictive of
methylation levels. Conclusions: Our observations of DNA methylation patterns
led us to develop a classifier to predict site-specific methylation levels that
achieves the best DNA methylation predictive accuracy to date. Furthermore, our
method identified genomic features that interact with DNA methylation,
elucidating mechanisms involved in DNA methylation modification and regulation,
and linking different epigenetic processes.


  The role of epistatic interactions among loci is a central question in
evolutionary biology and is increasingly relevant in the genomic age. While the
population genetics of compensatory substitution have received considerable
attention, most studies have focused on the case when natural selection is very
strong against deleterious intermediates. In the biologically-plausible
scenario of weak to moderate selection there exist two alternate pathways for
compensatory substitution. In one pathway, a deleterious mutation becomes fixed
prior to occurrence of the compensatory mutation. In the other, the two loci
are simultaneously polymorphic. The rates of compensatory substitution along
these two pathways and their relative probabilities are functions of the
population size, selection strength, mutation rate, and recombination rate. In
this paper these rates and path probabilities are derived analytically and
verified using population genetic simulations. The expected time durations of
these two paths are similar when selection is moderate, but not when selection
is weak. The effect of recombination on the dynamics of the substitution
process are explored using simulation. Using the derived rates, a phylogenetic
substitution model of the compensatory evolution process is presented that
could be used for inference of population genetic parameters from interspecific
data.


  Bacterial genomes and large-scale computer software projects both consist of
a large number of components (genes or software packages) connected via a
network of mutual dependencies. Components can be easily added or removed from
individual systems and their usage frequencies vary over many orders of
magnitude. We study this frequency distribution in genomes of ~500 bacterial
species and in over 2 million of Linux computers and find that in both cases it
is described by the same scale-free power law distribution with an additional
peak near the tail of the distribution corresponding to nearly universal
components. We argue that this is a general property of any modular system with
a multi-layered dependency network. We demonstrate that the frequency of a
component is positively correlated with its dependency degree given by the
total number of upstream components whose operation directly or indirectly
depends on the selected component. The observed frequency/dependency degree
distributions are reproduced in a simple mathematically tractable model
introduced and analyzed in this study.


  Genome wide association studies directly assay 10^6 single nucleotide
polymorphisms (SNPs) across a study cohort. Probabilistic estimation of
additional sites by genotype imputation can increase this set of variants by
10- to 40-fold. Even with modest sample sizes (10^3-10^4), these resulting
imputed datasets, containing 10^10-10^11 double-precision values, are
incompatible with simultaneous lossless storage in RAM using standard methods.
Existing solutions for this problem require compromises in either genotype
accuracy or complexity of permissible statistical methods. Here, we present a
C/C++ library that dynamically compresses probabilistic genotype data as they
are loaded into memory. This method uses a customization of the DEFLATE (gzip)
algorithm, and maintains constant-time access to any SNP. Average compression
ratios of more than 9-fold are observed in test data.


  Borrelia burgdorferi sensu stricto is a bacterial spirochete prevalent in the
Northeastern United States that causes Lyme disease. Lyme disease is the most
common arthropod-borne disease in the United States; affecting mice, deer,
humans and other mammals. The disease is spread by Ixodes Scapularis, a species
of tick whose primary food source are deer and mice. Reducing the population of
ticks feeding on both large and small mammals below some critical threshold can
decrease the prevalence of Lyme disease among humans. A simplified,
six-dimensional Susceptible-Infected, SI, model is used to capture the
mice-deer-tick dynamics while considering the impact of varying
population-specific death rates on infected population size. We analyzed the
stability of the models two equilibria, the unstable disease free equilibrium
and the endemic equilibrium. Static forward sensitivity analysis is conducted
on the basic reproduction number and the endemic equilibrium. A dynamic
approach was explored to observe change in the sensitivity of the death rates
over time. These analyses were conducted to determine the efficacy of changing
death rates in order to reduce prevalence of Lyme disease.


  We investigate evolutionary dynamics of two-strategy matrix games with
zealots in finite populations. Zealots are assumed to take either strategy
regardless of the fitness. When the strategy selected by the zealots is the
same, the fixation of the strategy selected by the zealots is a trivial
outcome. We study fixation time in this scenario. We show that the fixation
time is divided into three main regimes, in one of which the fixation time is
short, and in the other two the fixation time is exponentially long in terms of
the population size. Different from the case without zealots, there is a
threshold selection intensity below which the fixation is fast for an arbitrary
payoff matrix. We illustrate our results with examples of various social
dilemma games.


  The thesis examines in detail the folding and unfolding processes of a number
of proteins including hbSBD, DDLNF4, single and multi Ubiquitin. Using
simplified coarse-grained off-lattice Go model and CD experiments we have shown
the two-state behavior of hbSBD protein. It was shown that refolding pathways
of single Ubiquitin depend on what end is anchored to the surface. Namely, the
fixation of the N-terminal changes refolding pathways but anchoring the
C-terminal leaves them unchanged. Interestingly, the end fixation has no effect
on multi-domain Ubiquitin. Using the Go modeling and all-atom models with
explicit water, we have studied the mechanical unfolding mechanism of DDFLN4 in
detail. We predict that, contrary to the AFM experiments, an additional
unfolding peak would occur at the end-to-end $\Delta R \approx 1.5 $nm in the
force-extension curve. Our study reveals the important role of non-native
interactions which are responsible for a peak located at $\Delta R \approx 22
$nm. This peak can not be encountered by the Go models in which the non-native
interactions are neglected. Our finding may stimulate further experimental and
theoretical studies on this protein. The cooperativity of denaturation
transition was investigated using lattice and off-latice models. Our studies
reveal that the sharpness of the transition enhances as the number of amino
acids grows. The corresponding scaling behavior is governed by an universal
critical exponent. A new force replica exchange (RE) method was developed for
efficient configuration sampling of biomolecules pulled by an external
mechanical force. Contrary to the standard temperature RE, the exchange is
carried out between different forces (replicas). Our method was successfully
applied to study thermodynamics of a three-domain Ubiquitin.


  We describe a combination of all-atom simulations with CABS, a
well-established coarse-grained protein modeling tool, into a single multiscale
protocol. The simulation method has been tested on the C-terminal beta hairpin
of protein G, a model system of protein folding. After reconstructing atomistic
details, conformations derived from the CABS simulation were subjected to
replica-exchange molecular dynamics simulations with OPLS-AA and AMBER99sb
force fields in explicit solvent. Such a combination accelerates system
convergence several times in comparison with all-atom simulations starting from
the extended chain conformation, demonstrated by the analysis of melting
curves, the number of native-like conformations as a function of time and
secondary structure propagation. The results strongly suggest that the proposed
multiscale method could be an efficient and accurate tool for high-resolution
studies of protein folding dynamics in larger systems.


  1. Species distribution models (SDM) are tools used to determine
environmental features that influence the geographic distribution of species'
abundance and have been used to analyze presence-only records. Analysis of
presence-only records may require correction for nondetection sampling bias to
yield reliable conclusions. In addition, individuals of some species of animals
may be highly aggregated and standard SDMs ignore environmental features that
may influence aggregation behavior. 2. We contend that nondetection sampling
bias can be treated as missing data. Statistical theory and corrective methods
are well developed for missing data, but have been ignored in the literature on
SDMs. We developed a marked inhomogeneous Poisson point process model that
accounted for nondetection and aggregation behavior in animals and tested our
methods on simulated data. 3. Correcting for nondetection sampling bias
requires estimates of the probability of detection which must be obtained from
auxiliary data, as presence-only data do not contain information about the
detection mechanism. Weighted likelihood methods can be used to correct for
nondetection if estimates of the probability of detection are available. We
used an inhomogeneous Poisson point process model to model group abundance, a
zero-truncated generalized linear model to model group size, and combined these
two models to describe the distribution of abundance. Our methods performed
well on simulated data when nondetection was accounted for and poorly when
detection was ignored. 4. We recommend researchers consider the effects of
nondetection sampling bias when modeling species distributions using
presence-only data. If information about the detection process is available, we
recommend researchers explore the effects of nondetection and, when warranted,
correct the bias using our methods.


  A general framework for age-structured predator-prey systems is introduced.
Individuals are distinguished into two classes, juveniles and adults, and
several possible interactions are considered. The initial system of partial
differential equations is reduced to a system of (neutral) delay differential
equations with one or two delays. Thanks to this approach, physically correct
models for predator-prey with delay are provided. Previous models are
considered and analysed in view of the above results. A Rosenzweig-MacArthur
model with delay is presented as an example.


  A classical result in distance based tree-reconstruction characterizes when
for a distance $D$ on some finite set $X$ there exist a uniquely determined
dendrogram on $X$ (essentially a rooted tree $T=(V,E)$ with leaf set $X$ and no
degree two vertices but possibly the root and an edge weighting $\omega:E\to
\mathbb R_{\geq 0}$) such that the distance $D_{(T,\omega)}$ induced by
$(T,\omega)$ on $X$ is $D$. Moreover, algorithms that quickly reconstruct
$(T,\omega)$ from $D$ in this case are known. However in many areas where
dendrograms are being constructed such as Computational Biology not all
distances on $X$ are always available implying that the sought after dendrogram
need not be uniquely determined anymore by the available distances with regards
to topology of the underlying tree, edge-weighting, or both. To better
understand the structural properties a set $\cL\subseteq {X\choose 2}$ has to
satisfy to overcome this problem, various types of lassos have been introduced.
Here, we focus on the question of when a lasso uniquely determines the topology
of a dendrogram's underlying tree, that is, it is a topological lasso for that
tree. We show that any set-inclusion minimal topological lasso for such a tree
$T$ can be transformed into a 'distinguished' minimal topological lasso $\cL$
for $T$, that is, the graph $(X,\cL)$ is a claw-free block graph. Furthermore,
we characterize such lassos in terms of the novel concept of a cluster marker
map for $T$ and present results concerning the heritability of such lassos in
the context of the subtree and supertree problems.


  This paper studies the mechanisms, implications, and potential applications
of the recently discovered class of Zero Determinant (ZD) strategies in
iterated 2x2 games. These strategies were reported to successfully extort pure
economic maximizers, and to mischievously determine the set of feasible
long-term payoffs in iterated Prisoners' Dilemma by enforcing linear
constraints on both players' expected average scores.
  These results are generalized for all symmetric 2x2 games and a general
Battle of the Sexes, exemplified by four common games. Additionally, a
comparison to conventional strategies is made and typical ZD gameplay
simulations are analyzed along with convergence speeds. Several response
strategies are discussed, including a glance on how time preferences change
previous results. Furthermore, a possibility of retaliation is presented: when
maximin scores exceed the minimum symmetric payoff, it is possible to extort
the extortioner.
  Finally, a summary of findings from evolutionary game theory shows that
mischief is limited by its own malice. Nevertheless, this does not challenge
the result that mindless economic maximization is subject to extortion: the
study of ZD strategies reveals exciting new perspectives and opportunities in
game theory, both evolutionary and classic.


  Molecular motors are specialized proteins which perform active, directed
transport of cellular cargoes on cytoskeletal filaments. In many cases, cargo
motion powered by motor proteins is found to be bidirectional, and may be
viewed as a biased random walk with fast unidirectional runs interspersed with
slow `tug-of-war' states. The statistical properties of this walk are not known
in detail, and here, we study memory and bias, as well as directional
correlations between successive runs in bidirectional transport. We show, based
on a study of the direction reversal probabilities of the cargo using a purely
stochastic (tug-of-war) model, that bidirectional motion of cellular cargoes
is, in general, a correlated random walk. In particular, while the motion of a
cargo driven by two oppositely pulling motors is a Markovian random walk,
memory of direction appears when multiple motors haul the cargo in one or both
directions. In the latter case, the Markovian nature of the underlying single
motor processes is hidden by internal transitions between degenerate run and
pause states of the cargo. Interestingly, memory is found to be a non-monotonic
function of the number of motors. Stochastic numerical simulations of the
tug-of-war model support our mathematical results and extend them to
biologically relevant situations.


  Protein folding produces characteristic and functional three-dimensional
structures from unfolded polypeptides or disordered coils. The emergence of
extraordinary complexity in the protein folding process poses astonishing
challenges to theoretical modeling and computer simulations. The present work
introduces molecular nonlinear dynamics (MND), or molecular chaotic dynamics,
as a theoretical framework for describing and analyzing protein folding. We
unveil the existence of intrinsically low dimensional manifolds (ILDMs) in the
chaotic dynamics of folded proteins. Additionally, we reveal that the
transition from disordered to ordered conformations in protein folding
increases the transverse stability of the ILDM. Stated differently, protein
folding reduces the chaoticity of the nonlinear dynamical system, and a folded
protein has the best ability to tame chaos. Additionally, we bring to light the
connection between the ILDM stability and the thermodynamic stability, which
enables us to quantify the disorderliness and relative energies of folded,
misfolded and unfolded protein states. Finally, we exploit chaos for protein
flexibility analysis and develop a robust chaotic algorithm for the prediction
of Debye-Waller factors, or temperature factors, of protein structures.


  An ongoing debate in ecology concerns the impacts of ecological drift and
selection on community assembly. Here, we show that there is a sharp phase
transition in diverse ecological communities between a selection dominated
regime (the niche phase) and a drift dominated regime (the neutral phase).
Simulations and analytic arguments show that the niche phase is favored in
communities with large population sizes and relatively constant environments,
whereas the neutral phase is favored in communities with small population sizes
and fluctuating environments. Our results demonstrate how apparently neutral
populations may arise even in communities inhabited by species with varying
traits.


  Stern-judging is one of the best-known assessment rules in indirect
reciprocity. Indirect reciprocity is a fundamental mechanism for the evolution
of cooperation. It relies on mutual monitoring and assessments, i.e.,
individuals judge, following their own assessment rules, whether other
individuals are "good" or "bad" according to information on their past
behaviors. Among many assessment rules, stern-judging is known to provide
stable cooperation in a population, as observed when all members in the
population know all about others' behaviors (public information case) and when
the members never commit an assessment error. In this paper, the effect of
assessment error and private information on stern-judging is investigated. By
analyzing the image matrix, which describes who is good in the eyes of whom in
the population, we analytically show that private information and assessment
error cause the collapse of stern-judging: all individuals assess other
individuals as "good" at random with a probability of 1/2.


  Quantum mechanical calculations are performed on the proteins that constitute
the ubiquitin-Dsk2 complex whose atomic structure has been experimentally
determined by NMR spectroscopy (PDB id 1WR1). The results indicate that the
dipole moment vectors of the two proteins are aligned in a head-to-tail
orientation while forming and angle of ~130{\deg}. Hence, attractive
dipole-dipole interactions not only stabilize the protein-protein complex but
they are likely to favor the correct orientation of the proteins during the
formation of the complex.


  We discuss the use of a structure based C$\alpha$-Go model and Langevin
dynamics to study in detail the mechanical properties and unfolding pathway of
the titin I27 domain. We show that a simple Go-model does detect correctly the
origin of the mechanical stability of this domain. The unfolding free energy
landscape parameters $x_u$ and $\Delta G^{\ddagger}$, extracted from
dependencies of unfolding forces on pulling speeds, are found to agree
reasonably well with experiments. We predict that above $v=10^4$ nm/s the
additional force-induced intermediate state is populated at an end-to-end
extension of about $75 \mathring{A}$. The force-induced switch in the unfolding
pathway occurs at the critical pulling speed $v_{crit} \approx 10^6-10^7$ nm/s.
We argue that this critical pulling speed is an upper limit of the interval
where Bell's theory works. However, our results suggest that the Go-model fails
to reproduce the experimentally observed mechanical unfolding pathway properly,
yielding an incomplete picture of the free energy landscape. Surprisingly, the
experimentally observed intermediate state with the A strand detached is not
populated in Go-model simulations over a wide range of pulling speeds. The
discrepancy between simulation and experiment is clearly seen from the early
stage of the unfolding process which shows the limitation of the Go model in
reproducing unfolding pathways and deciphering the complete picture of the free
energy landscape.


  Stochasticity is an indispensable aspect of biochemical processes at the
cellular level. Studies on how the noise enters and propagates in biochemical
systems provided us with nontrivial insights into the origins of stochasticity,
in total however they constitute a patchwork of different theoretical analyses.
Here we present a flexible and generally applicable noise decomposition tool,
that allows us to calculate contributions of individual reactions to the total
variability of a system's output. With the package it is therefore possible to
quantify how the noise enters and propagates in biochemical systems. We also
demonstrate and exemplify using the JAK-STAT signalling pathway that it is
possible to infer noise contributions resulting from individual reactions
directly from experimental data. This is the first computational tool that
allows to decompose noise into contributions resulting from individual
reactions.


  If two species exhibit different nonlinear responses to a single shared
resource, and if each species modifies the resource dynamics such that this
favors its competitor, they may stably coexist. This coexistence mechanism,
known as relative nonlinearity of competition, is well understood
theoretically, but less is known about its evolutionary properties and its
prevalence in real communities. We address this challenge by using adaptive
dynamics theory and individual-based simulations to compare community
stabilization and evolutionary stability of species that coexist by relative
nonlinearity. In our analysis, evolution operates on the species'
density-compensation strategies, and we consider a trade-off between population
growth rates at high and low resource availability. We confirm previous
findings that, irrespective of the particular model of density dependence,
there are many combinations of overcompensating and undercompensating
density-compensation strategies that allow stable coexistence by relative
nonlinearity. However, our analysis also shows that most of these strategy
combinations are not evolutionarily stable and will be outcompeted by an
intermediate density-compensation strategy. Only very specific trade-offs lead
to evolutionarily stable coexistence by relative nonlinearity. As we find no
reason why these particular trade-offs should be common in nature, we conclude
that the sympatric evolution and evolutionary stability of relative
nonlinearity, while possible in principle, seems rather unlikely. We speculate
that this may, at least in part, explain why empirical demonstrations of this
coexistence mechanism are rare, noting, however, that the difficulty to detect
relative nonlinearity in the field [...]


  In line with the importance of RNA-seq, the bioinformatics community has
produced numerous data analysis tools incorporating methods to correct
sample-specific biases. However, few advanced simulation tools exist to enable
benchmarking of competing correction methods. We introduce the first framework
to reproduce the properties of individual RNA-seq runs and, by applying it on
several datasets, we demonstrate the importance of accounting for
sample-specificity in realistic simulations.


  The ominous warnings of a `data deluge' in the life sciences from
high-throughput DNA sequencing data are being supplanted by a second deluge, of
cliches bemoaning our collective scientific fate unless we address the genomic
data `tsunami'. It is imperative that we explore the many facets of the genome,
not just sequence but also transcriptional and epigenetic variability,
integrating these observations in order to attain a genuine understanding of
how genes function, towards a goal of genomics-based personalized medicine.
Determining any individual's genomic properties requires comparison to many
others, sifting out the specific from the trends, requiring access to the many
in order to yield information relevant to the few. This is the central big data
challenge in genomics that still requires some sort of resolution. Is there a
practical, feasible way of directly connecting the scientific community to this
data universe? The best answer could be in the stars overhead.


  Since 1927, until recently, models describing the spread of disease have
mostly been of the SIR-compartmental type, based on the assumption that
populations are homogeneous and well-mixed. The focus of these models have
typically been on large-scale analysis of scenarios such as cities, nations or
even world scale. SIR models are appealing because of their simplicity, but
their parameters, especially the transmission rate, are complex and depend on a
number of factors, which makes it hard to predict how a change of a single
environmental, demographic, or epidemiological factor will affect the
population. Therefore, in this contribution we start to unpick the
transmission-rate parameter. Analysing the implications that arise when taking
crowd behaviour explicitly into account, we show how both the rate of infection
as well as the walking speed depend on the local crowd density around an
infected individual. The combined effect is that the rate of infection at a
population scale has an analytically tractable non-linear dependency on crowd
density. We model the spread of a hypothetical disease in a corridor and
compare our new model with a typical compartmental model, which highlights the
regime in which current models may not produce credible results.


  Tumors bearing Ras mutations are notoriously difficult to treat. Drug
combinations targeting the Ras protein or its pathway have also not met with
success. Pathway drug cocktails, which are combinations aiming at parallel
pathways, appear more promising; however, to be usefully exploited, a
repertoire of classified pathway combinations is desirable. This challenge
would be facilitated by the availability of the structural network of signaling
pathways. When integrated with functional and systems-level clinical data they
can be powerful in advancing novel therapeutic platforms. Based on structural
knowledge, drug cocktails may tear into multiple cellular processes that drive
tumorigenesis, and help in deciphering the interrelationship between Ras
mutations and the rewired Ras network. The pathway drug cocktail paradigm can
be applied to other signaling protein targets.


  More often than not, bad decisions are bad regardless of where and when they
are made. Information sharing might thus be utilized to mitigate them. Here we
show that sharing the information about strategy choice between players
residing on two different networks reinforces the evolution of cooperation. In
evolutionary games the strategy reflects the action of each individual that
warrants the highest utility in a competitive setting. We therefore assume that
identical strategies on the two networks reinforce themselves by lessening
their propensity to change. Besides network reciprocity working in favour of
cooperation on each individual network, we observe the spontaneous emerge of
correlated behaviour between the two networks, which further deters defection.
If information is shared not just between individuals but also between groups,
the positive effect is even stronger, and this despite the fact that
information sharing is implemented without any assumptions with regards to
content.


  Global environmental change is affecting species distribution and their
interactions with other species. In particular, the main drivers of
environmental change strongly affect the strength of interspecific interactions
with considerable consequences to biodiversity. However, extrapolating the
effects observed on pair-wise interactions to entire ecological networks is
challenging. Here we propose a framework to estimate the tolerance to changes
in the strength of mutualistic interaction that species in mutualistic networks
can sustain before becoming extinct. We identify the scenarios where generalist
species can be the least tolerant. We show that the least tolerant species
across different scenarios do not appear to have uniquely common
characteristics. Species tolerance is extremely sensitive to the direction of
change in the strength of mutualistic interaction, as well as to the observed
mutualistic trade-offs between the number of partners and the strength of the
interactions.


  In this paper, new techniques that allow conditional entropy to estimate the
combinatorics of symbols are applied to animal communication studies to
estimate the communication's repertoire size. By using the conditional entropy
estimates at multiple orders, the paper estimates the total repertoire sizes
for animal communication across bottlenose dolphins, humpback whales, and
several species of birds for N-grams length one to three. In addition to
discussing the impact of this method on studies of animal communication
complexity, the reliability of these estimates is compared to other methods
through simulation. While entropy does undercount the total repertoire size due
to rare N-grams, it gives a more accurate picture of the most frequently used
repertoire than just repertoire size alone.


  RNA-seq has rapidly become the de facto technique to measure gene expression.
However, the time required for analysis has not kept up with the pace of data
generation. Here we introduce Sailfish, a novel computational method for
quantifying the abundance of previously annotated RNA isoforms from RNA-seq
data. Sailfish entirely avoids mapping reads, which is a time-consuming step in
all current methods. Sailfish provides quantification estimates much faster
than existing approaches (typically 20-times faster) without loss of accuracy.


  Quantum mechanical calculations are performed on 116 conformers of the
protein ubiquitin (Lange et al., Science 2008, 320, 1471-1475). The results
indicate that the heat of formation (HOF), dipole moment, energy of the
frontier orbitals HOMO and LUMO, and HOMO-LUMO gap fluctuate within their
corresponding ranges. This study thus provides a link between the
conformational dynamics of a protein and its electronic structure.


  To simulate long time and length scale processes involving DNA it is
necessary to use a coarse-grained description. Here we provide an overview of
different approaches to such coarse graining, focussing on those at the
nucleotide level that allow the self-assembly processes associated with DNA
nanotechnology to be studied. OxDNA, our recently-developed coarse-grained DNA
model, is particularly suited to this task, and has opened up this field to
systematic study by simulations. We illustrate some of the range of DNA
nanotechnology systems to which the model is being applied, as well as the
insights it can provide into fundamental biophysical properties of DNA.


  The sustainable use of multicomponent treatments such as combination
therapies, combination vaccines/chemicals, and plants carrying multigenic
resistance requires an understanding of how their population-wide deployment
affects the speed of the pathogen adaptation. Here, we develop a stochastic
model describing the emergence of a mutant pathogen and its dynamics in a
heterogeneous host population split into various types by the management
strategy. Based on a multi-type Markov birth and death process, the model can
be used to provide a basic understanding of how the life-cycle parameters of
the pathogen population, and the controllable parameters of a management
strategy affect the speed at which a pathogen adapts to a multicomponent
treatment. Our results reveal the importance of coupling stochastic mutation
and migration processes, and illustrate how their stochasticity can alter our
view of the principles of managing pathogen adaptive dynamics at the population
level. In particular, we identify the growth and migration rates that allow
pathogens to adapt to a multicomponent treatment even if it is deployed on only
small proportions of the host. In contrast to the accepted view, our model
suggests that treatment durability should not systematically be identified with
mutation cost. We show also that associating a multicomponent treatment with
defeated monocomponent treatments can be more durable than associating it with
intermediate treatments including only some of the components. We conclude that
the explicit modelling of stochastic processes underlying evolutionary dynamics
could help to elucidate the principles of the sustainable use of multicomponent
treatments in population-wide management strategies intended to impede the
evolution of harmful populations.


  A novel self-assembly strategy for polypeptide nanostructure design was
presented in [Design of a single-chain polypeptide tetrahedron assembled from
coiled-coil segments, Nature Chemical Biology 9 (2013) 362--366]. The first
mathematical model (polypeptide nanostructure can naturally be presented as a
skeleton graph of a polyhedron) from [Stable traces as a model for
self-assembly of polypeptide nanoscale polyhedrons, MATCH Commun. Math. Comput.
Chem. 70 (2013) 317-330] introduced stable traces as the appropriate
mathematical description, yet we find them deficient in modeling graphs with
either very small (less or equal to 2) or large (greater or equal to 6) degree
vertices. We introduce strong traces which remedy both of the above mentioned
drawbacks. We show that every connected graph admits a strong trace by studying
a connection between strong traces and graph embeddings. Further we also
characterize graphs which admit parallel (resp. antiparallel) strong traces.


  Major unresolved questions in evolutionary genetics include determining the
contributions of different mutational sources to the total pool of genetic
variation in a species, and understanding how these different forms of genetic
variation interact with natural selection. Recent work has shown that
structural variants (insertions, deletions, inversions and transpositions) are
a major source of genetic variation, often out-numbering single nucleotide
variants in terms of total bases affected. Despite the near ubiquity of
structural variants, major questions about their interaction with natural
selection remain. For example, how does the allele frequency spectrum of
structural variants differ when compared to single nucleotide variants? How
often do structural variants affect genes, and what are the consequences? To
begin to address these questions, we have systematically identified and
characterized a large set submicroscopic insertion and deletion (indel)
variants (between 1 kb to 200 kb in length) among ten individuals from a single
natural population of the plant species Mimulus guttatus. After extensive
computational filtering, we focused on a set of 4,142 high-confidence indels
that showed an experimental validation rate of 73%. All but one of these indels
were < 200 kb. While the largest were generally at lower frequencies in the
population, a surprising number of large indels are at intermediate
frequencies. While indels overlapping with genes were much rarer than expected
by chance, nearly 600 genes were affected by an indel. NBS-LRR defense response
genes were the most enriched among the gene families affected. Most indels
associated with genes were rare and appeared to be under purifying selection,
though we do find four high-frequency derived insertion alleles that show
signatures of recent positive selection.


  Genome-wide association studies (GWAS) have identified hundreds of loci at
very stringent levels of statistical significance across many different human
traits. However, it is now clear that very large samples (n~10^4-10^5) are
needed to find the majority of genetic variants underlying risk for most human
diseases. Therefore, the field has engaged itself in a race to increase study
sample sizes with some studies yielding very successful results but also
studies which provide little or no new insights. This project started early on
in this new wave of studies and I decided to use an alternative approach that
uses prior biological knowledge to improve both interpretation and power of
GWAS. The project aimed to a) implement and develop new gene-based methods to
derive gene-level statistics to use GWAS in well established system biology
tools; b) use of these gene-level statistics in networks and gene-set analyses
of GWAS data; c) mine GWAS of neuropsychiatric disorders using gene, gene-sets
and integrative biology analyses with gene-expression studies; and d) explore
the ability of these methods to improve the analysis GWAS on disease
sub-phenotypes which usually suffer of very small sample sizes.


  Van der Waals density functional theory is integrated with analysis of a
non-redundant set of protein-DNA crystal structures from the Nucleic Acid
Database to study the stacking energetics of CG:CG base-pair steps,
specifically the role of cytosine 5-methylation. Principal component analysis
of the steps reveals the dominant collective motions to correspond to a tensile
'opening' mode and two shear 'sliding' and 'tearing' modes in the orthogonal
plane. The stacking interactions of the methyl groups globally inhibit CG:CG
step overtwisting while simultaneously softening the modes locally via
potential energy modulations that create metastable states. Additionally, the
indirect effects of the methyl groups on possible base-pair steps neighboring
CG:CG are observed to be of comparable importance to their direct effects on
CG:CG. The results have implications for the epigenetic control of DNA
mechanics.