  We study the effect of spatial structure, genetic drift, mutation, and
selective pressure on the evolutionary dynamics in a simplified model of
asexual organisms colonizing a new territory. Under an appropriate
coarse-graining, the evolutionary dynamics is related to the directed
percolation processes that arise in voter models, the Domany-Kinzel (DK) model,
contact process, etc. We explore the differences between linear (flat front)
expansions and the much less familiar radial (curved front) range expansions.
For the radial expansion, we develop a generalized, off-lattice DK model that
minimizes otherwise persistent lattice artifacts. With both simulations and
analytical techniques, we study the survival probability of advantageous
mutants, the spatial correlations between domains of neutral strains, and the
dynamics of populations with deleterious mutations. "Inflation" at the frontier
leads to striking differences between radial and linear expansions. For a
colony with initial radius $R_0$ expanding at velocity $v$, significant genetic
demixing, caused by local genetic drift, only occurs up to a finite time $t^* =
R_0/v$, after which portions of the colony become causally disconnected due to
the inflating perimeter of the expanding front. As a result, the effect of a
selective advantage is amplified relative to genetic drift, increasing the
survival probability of advantageous mutants. Inflation also modifies the
underlying directed percolation transition, introducing novel scaling functions
and modifications similar to a finite size effect. Finally, we consider radial
range expansions with deflating perimeters, as might arise from colonization
initiated along the shores of an island.


  DNA samples are often pooled, either by experimental design, or because the
sample itself is a mixture. For example, when population allele frequencies are
of primary interest, individual samples may be pooled together to lower the
cost of sequencing. Alternatively, the sample itself may be a mixture of
multiple species or strains (e.g. bacterial species comprising a microbiome, or
pathogen strains in a blood sample). We present an expectation-maximization
(EM) algorithm for estimating haplotype frequencies in a pooled sample directly
from mapped sequence reads, in the case where the possible haplotypes are
known. This method is relevant to the analysis of pooled sequencing data from
selection experiments, as well as the calculation of proportions of different
strains within a metagenomics sample. Our method outperforms existing methods
based on single- site allele frequencies, as well as simple approaches using
sequence read data. We have implemented the method in a freely available
open-source software tool.


  Homologous recombination is an important operator in the evolution of
biological organisms. However, there is still no clear, generally accepted
understanding of why it exists and under what circumstances it is useful. In
this paper we consider its utility in the context of an infinite population
haploid model with selection and homologous recombination. We define utility in
terms of two metrics - the increase in frequency of fit genotypes, and the
increase in average population fitness, relative to those associated with
selection only. Explicitly, we exhaustively explore the eight-dimensional
parameter space of a two-locus two-allele system, showing, as a function of the
landscape and the initial population, that recombination is beneficial in terms
of our metrics in two distinct regimes: a landscape independent regime - the
"search" regime - where recombination aids in the search for a fit genotype
that is absent or at low frequency in the population; and the "modular" regime,
associated with quasi-additive fitness landscapes with low epistasis, where
recombination allows for the juxtaposition of fit "modules" or Building Blocks.
Thus, we conclude that the ubiquity and utility of recombination is intimately
associated with the existence of modularity in biological fitness landscapes.


  We present simulations that reveal a surprisingly large effect of
hydrodynamic coupling on the speed of thermal ratchet motors. The model that we
use considers particles performing thermal ratchet motion in a hydrodynamic
solvent. Using particle-based, mesoscopic simulations that maintain local
momentum conservation, we analyze quantitatively how the coupling to the
surrounding fluid affects ratchet motion. We find that coupling can increase
the mean velocity of the moving particles by almost two orders of magnitude,
precisely because ratchet motion has both a diffusive and a deterministic
component. The resulting coupling also leads to the formation of aggregates at
longer times. The correlated motion that we describe increases the efficiency
of motor-delivered cargo transport and we speculate that the mechanism that we
have uncovered may play a key role in speeding up molecular motor-driven
intracellular transport.


  Alzheimer's disease is a human brain disease that affects a significant
fraction of the population by causing problems with short-term memory,
thinking, spatial orientation and behavior, memory loss and other intellectual
abilities. Up to date there is no singular test that can definitively diagnose
Alzheimer's disease, although imaging technology designed to detect Alzheimer's
plaques and tangles is rapidly becoming more powerful and precise. In this
paper we introduce a decision-making model, based on the combination of
mitochondrial hypothesis-dynamics with the role of electromagnetic influences
of the metal ions into the inner mitochondrial membrane and the quantitative
analysis of mitochondrial population. While there are few disappointing
clinical-trial results for drug treatments in patients with Alzheimer's
disease, scientific community need alternative diagnostic tools rather
investing mainly in amyloid-targeting drugs.


  Hybridization between species is an important mechanism for the origin of
novel lineages and adaptation to new environments. Increased allelic variation
and modification of the transcriptional network are the two recognized forces
currently deemed to be responsible for the phenotypic properties seen in
hybrids. However, since the majority of the biological functions in a cell are
carried out by protein complexes, inter-specific protein assemblies therefore
represent another important source of natural variation upon which evolutionary
forces can act. Here we studied the composition of six protein complexes in two
different Saccharomyces "sensu strictu" hybrids, to understand whether chimeric
interactions can be freely formed in the cell in spite of species-specific
co-evolutionary forces, and whether the different types of complexes cause a
change in hybrid fitness. The protein assemblies were isolated from the hybrids
via affinity chromatography and identified via mass spectrometry. We found
evidence of spontaneous chimericity for four of the six protein assemblies
tested and we showed that different types of complexes can cause a variety of
phenotypes in selected environments. In the case of TRP2/TRP3 complex, the
effect of such chimeric formation resulted in the fitness advantage of the
hybrid in an environment lacking tryptophan, while only one type of parental
combination of the MBF complex could confer viability to the hybrid under
respiratory conditions. This study provides empirical evidence that chimeric
protein complexes can freely assemble in cells and reveals a new mechanism to
generate phenotypic novelty and plasticity in hybrids to complement the genomic
innovation resulting from gene duplication. The ability to exchange orthologous
members has also important implications for the adaptation and subsequent
genome evolution of the hybrids in terms of pattern of gene loss.


  Influenza virus contains two highly variable envelope glycoproteins,
hemagglutinin (HA) and neuraminidase (NA). The structure and properties of HA,
which is responsible for binding the virus to the cell that is being infected,
change significantly when the virus is transmitted from avian or swine species
to humans. Previously we identified much smaller human individual evolutionary
amino acid mutational changes in NA, which cleaves sialic acid groups and is
required for influenza virus replication. We showed that these smaller changes
can be monitored very accurately across many Uniprot and NCBI strains using
hydropathicity scales to quantify the roughness of water film packages, which
increases gradually due to migration, but decreases abruptly under large-scale
vaccination pressures. Here we show that, while HA evolution is much more
complex, it still shows abrupt punctuation changes linked to those of NA. HA
exhibits proteinquakes, which resemble earthquakes and are related to
hydropathic shifting of sialic acid binding regions. HA proteinquakes based on
sialic acid interactions are required for optimal balance between the
receptor-binding and receptor-destroying activities of HA and NA for efficient
virus replication. Our comprehensive results present an historical (1945-2011)
panorama of HA evolution over thousands of strains, and are consistent with
many studies of HA and NA interactions based on a few mutations of a few
strains. While the common influenza virus discussed here has been rendered
almost harmless by decades of vaccination programs, the sequential decoding
lessons learned here are applicable to other viruses that are emerging as
powerful weapons for controlling and even curing common organ cancers. Those
engineered oncolytic drugs will be discussed in future papers.


  Based on the oligomer-world hypothesis we propose an abstract model where the
molecular recognition among oligomers is described in the shape space. The
origin of life in the oligomer world is regarded as the establishment of a
metabolic cycle in a primitive cell. The cycle is sustained by the molecular
recognition. If an original cell acquires the ability of the replication of
oligomers, the relationship among oligomers changes due to the poor fidelity of
the replication. This change leads to the diversification of metabolic cycles.
The selection among diverse cycles is the basis of the evolution. The
evolvability is one of the essential characters of life. We demonstrate the
origin and diversification of the metabolic cycle by the computer simulation of
our model. Such a simulation is expected to be the simplified demonstration of
what actually occurred in the primordial soup. Our model describes an analog
era preceding the digital era based on the genetic code.


  Our previously presented method for high throughput computational screening
of mutant activity (Hediger et al., arXiv:1203.2950) is benchmarked against
experimentally measured amidase activity for 22 mutants of Candida antarctica
lipase B (CalB). Using an appropriate cutoff criterion for the computed
barriers, the qualitative activity of 15 out of 22 mutants is correctly
predicted. The method identifies four of the six most active mutants with
>=3-fold wild type activity and seven out of the eight least active mutants
with <=0.5-fold wild type activity. The method is further used to screen all
sterically possible (386) double-, triple- and quadruple-mutants constructed
from the most active single mutants. Based on the benchmark test at least 20
new promising mutants are identified.


  It has been known for some time that human autoimmune diseases can be
triggered by viral infections. Several possible mechanisms of interactions
between a virus and immune system have been analysed, with a prevailing opinion
being that the onset of autoimmunity can in many cases be attributed to
"molecular mimicry", where linear peptide epitopes, processed from viral
proteins, mimic normal host self proteins, thus leading to a cross-reaction of
immune response against virus with host cells. In this paper we present a
mathematical model for the dynamics of an immune response to a viral infection
and autoimmunity, which takes into account T cells with different activation
thresholds. We show how the infection can be cleared by the immune system, as
well as how it can lead to a chronic infection or recurrent infection with
relapses and remissions. Numerical simulations of the model are performed to
illustrate various dynamical regimes, as well as to analyse the potential
impact of treatment of autoimmune disease in the chronic and recurrent states.
The results provide good qualitative agreement with available data on immune
responses to viral infections and progression of autoimmune diseases.


  We study stochastic copying schemes in which discrimination between a right
and a wrong match is achieved via different kinetic barriers or different
binding energies of the two matches. We demonstrate that, in single-step
reactions, the two discrimination mechanisms are strictly alternative and can
not be mixed to further reduce the error fraction. Close to the lowest error
limit, kinetic discrimination results in a diverging copying velocity and
dissipation per copied bit. On the opposite, energetic discrimination reaches
its lowest error limit in an adiabatic regime where dissipation and velocity
vanish. By analyzing experimentally measured kinetic rates of two DNA
polymerases, T7 and Pol{\gamma}, we argue that one of them operates in the
kinetic and the other in the energetic regime. Finally, we show how the two
mechanisms can be combined in copying schemes implementing error correction
through a proofreading pathway


  How do mammalian cells that share the same genome exist in notably distinct
phenotypes, exhibiting differences in morphology, gene expression patterns, and
epigenetic chromatin statuses? Furthermore how do cells of different phenotypes
differentiate reproducibly from a single fertilized egg? These fundamental
questions are closely related to a deeply rooted paradigm in developmental
biology that cell differentiation is irreversible. Yet, recently a growing body
of research suggests the possibility of cell reprogramming, which offers the
potential for us to convert one type of cell into another. Despite the
significance of quantitative understandings of cell reprogramming, theoretical
efforts often suffer from the complexity of large circuits maintaining cell
phenotypes coupled at many different epigenetic and gene regulation levels. To
capture the global architecture of cell phenotypes, we propose an "epigenetic
state network" approach that translates the classical concept of an epigenetic
landscape into a simple-yet-predictive mathematical model. As a testing case,
we apply the approach to the reprogramming of fibroblasts (FB) to
cardiomyocytes (CM). The epigenetic state network for this case predicts three
major pathways of reprogramming. One pathway goes by way of induced pluripotent
stem cells (iPSC) and continues on to the normal pathway of cardiomyocyte
differentiation. The other two pathways involve transdifferentiation (TD)
either indirectly through cardiac progenitor (CP) cells or directly from
fibroblast to cardiomyocyte. Numerous experimental observations support the
predicted states and pathways.


  Centromeres are essential for chromosome segregation, yet their DNA sequences
evolve rapidly. In most animals and plants that have been studied, centromeres
contain megabase-scale arrays of tandem repeats. Despite their importance, very
little is known about the degree to which centromere tandem repeats share
common properties between different species across different phyla. We used
bioinformatic methods to identify high-copy tandem repeats from 282 species
using publicly available genomic sequence and our own data. The assumption that
the most abundant tandem repeat is the centromere DNA was true for most species
whose centromeres have been previously characterized, suggesting this is a
general property of genomes. Our methods are compatible with all current
sequencing technologies. Long Pacific Biosciences sequence reads allowed us to
find tandem repeat monomers up to 1,419 bp. High-copy centromere tandem repeats
were found in almost all animal and plant genomes, but repeat monomers were
highly variable in sequence composition and in length. Furthermore,
phylogenetic analysis of sequence homology showed little evidence of sequence
conservation beyond ~50 million years of divergence. We find that despite an
overall lack of sequence conservation, centromere tandem repeats from diverse
species showed similar modes of evolution, including the appearance of higher
order repeat structures in which several polymorphic monomers make up a larger
repeating unit. While centromere position in most eukaryotes is epigenetically
determined, our results indicate that tandem repeats are highly prevalent at
centromeres of both animals and plants. This suggests a functional role for
such repeats, perhaps in promoting concerted evolution of centromere DNA across
chromosomes.


  Plasmodium falciparum aquaglyceroporin (PfAQP) is a multifunctional membrane
protein in the plasma membrane of P. falciparum, the parasite that causes the
most severe form of malaria. The current literature has established the science
of PfAQP's structure, functions, and hydrogen-bonding interactions but left
unanswered the following fundamental question: Does glycerol modulate water
permeation through aquaglyceroporin that conducts both glycerol and water? This
paper provides an affirmative answer to this question of essential importance
to the protein's functions. On the basis of the chemical-potential profile of
glycerol from the extracellular bulk region, throughout PfAQP's conducting
channel, to the cytoplasmic bulk region, this study shows the existence of a
bound state of glycerol inside aquaglyceroporin's permeation pore, from which
the dissociation constant is approximately 14 microM. A glycerol molecule
occupying the bound state occludes the conducting pore through which permeating
molecules line up in single file by hydrogen-bonding with one another and with
the luminal residues of aquaglyceroporin. In this way, glycerol inhibits
permeation of water and other permeants through aquaglyceroporin. The
biological implications of this theory are discussed and shown to agree with
the existent in vitro data. It turns out that the structure of aquaglyceroporin
is perfect for the van der Waals interactions between the protein and glycerol
to cause the existence of the bound state deep inside the conducting pore and,
thus, to play an unexpected but significant role in aquaglyceroporin's
functions.


  One way to understand the role history plays on evolutionary trajectories is
by giving ancient life a second opportunity to evolve. Our ability to
empirically perform such an experiment, however, is limited by current
experimental designs. Combining ancestral sequence reconstruction with
synthetic biology allows us to resurrect the past within a modern context and
has expanded our understanding of protein functionality within a historical
context. Experimental evolution, on the other hand, provides us with the
ability to study evolution in action, under controlled conditions in the
laboratory. Here we describe a novel experimental setup that integrates two
disparate fields - ancestral sequence reconstruction and experimental
evolution. This allows us to rewind and replay the evolutionary history of
ancient biomolecules in the laboratory. We anticipate that our combination will
provide a deeper understanding of the underlying roles that contingency and
determinism play in shaping evolutionary processes.


  We consider a birth-death process with the birth rates $i\lambda$ and death
rates $i\mu +i(i-1)\theta$, where $i$ is the current state of the process. A
positive competition rate $\theta$ is assumed to be small. In the supercritical
case when $\lambda>\mu$ this process can be viewed as a demographic model for a
population with a high carrying capacity around $\lambda-\mu\over\theta$.
  The article reports in a self-contained manner on the asymptotic properties
of the time to extinction for this logistic branching process as $\theta\to0$.
All three reproduction regimes $\lambda>\mu$, $\lambda<\mu$, and $\lambda=\mu$
are studied.


  The linear noise approximation (LNA) offers a simple means by which one can
study intrinsic noise in monostable biochemical networks. Using simple physical
arguments, we have recently introduced the slow-scale LNA (ssLNA) which is a
reduced version of the LNA under conditions of timescale separation. In this
paper, we present the first rigorous derivation of the ssLNA using the
projection operator technique and show that the ssLNA follows uniquely from the
standard LNA under the same conditions of timescale separation as those
required for the deterministic quasi-steady state approximation. We also show
that the large molecule number limit of several common stochastic model
reduction techniques under timescale separation conditions constitutes a
special case of the ssLNA.


  Patterns of isolation-by-distance arise when population differentiation
increases with increasing geographic distances. Patterns of
isolation-by-distance are usually caused by local spatial dispersal, which
explains why differences of allele frequencies between populations accumulate
with distance. However, spatial variations of demographic parameters such as
migration rate or population density can generate non-stationary patterns of
isolation-by-distance where the rate at which genetic differentiation
accumulates varies across space. To characterize non-stationary patterns of
isolation-by-distance, we infer local genetic differentiation based on Bayesian
kriging. Local genetic differentiation for a sampled population is defined as
the average genetic differentiation between the sampled population and fictive
neighboring populations. To avoid defining populations in advance, the method
can also be applied at the scale of individuals making it relevant for
landscape genetics. Inference of local genetic differentiation relies on a
matrix of pairwise similarity or dissimilarity between populations or
individuals such as matrices of FST between pairs of populations. Simulation
studies show that maps of local genetic differentiation can reveal barriers to
gene flow but also other patterns such as continuous variations of gene flow
across habitat. The potential of the method is illustrated with 2 data sets:
genome-wide SNP data for human Swedish populations and AFLP markers for alpine
plant species. The software LocalDiff implementing the method is available at
http://membres-timc.imag.fr/Michael.Blum/LocalDiff.html


  The necks of the sauropod dinosaurs reached 15 m in length: six times longer
than that of the world record giraffe and five times longer than those of all
other terrestrial animals. Several anatomical features enabled this extreme
elongation, including: absolutely large body size and quadrupedal stance
providing a stable platform for a long neck; a small, light head that did not
orally process food; cervical vertebrae that were both numerous and
individually elongate; an efficient air-sac-based respiratory system; and
distinctive cervical architecture. Relevant features of sauropod cervical
vertebrae include: pneumatic chambers that enabled the bone to be positioned in
a mechanically efficient way within the envelope; and muscular attachments of
varying importance to the neural spines, epipophyses and cervical ribs. Other
long-necked tetrapods lacked important features of sauropods, preventing the
evolution of longer necks: for example, giraffes have relatively small torsos
and large, heavy heads, share the usual mammalian constraint of only seven
cervical vertebrae, and lack an air-sac system and pneumatic bones. Among
non-sauropods, their saurischian relatives the theropod dinosaurs seem to have
been best placed to evolve long necks, and indeed they probably surpassed those
of giraffes. But 150 million years of evolution did not suffice for them to
exceed a relatively modest 2.5 m.


  Maize Abnormal chromosome 10 (Ab10) contains a classic meiotic drive system
that exploits asymmetry of meiosis to preferentially transmit itself and other
chromosomes containing specialized heterochromatic regions called knobs. The
structure and diversity of the Ab10 meiotic drive haplotype is poorly
understood. We developed a BAC library from an Ab10 line and used the data to
develop sequence-based markers, focusing on the proximal portion of the
haplotype that shows partial homology to normal chromosome 10. These molecular
and additional cytological data demonstrate that two previously identified Ab10
variants (Ab10-I and Ab10-II) share a common origin. Dominant PCR markers were
used with FISH to assay 160 diverse teosinte and maize landrace populations
from across the Americas, resulting in the identification of a previously
unknown but prevalent form of Ab10 (Ab10-III). We find that Ab10 occurs in at
least 75% of teosinte populations at a mean frequency of 15%. Ab10 was also
found in 13% of the maize landraces, but does not appear to be fixed in any
wild or cultivated population. Quantitative analyses suggest that the abundance
and distribution of Ab10 is governed by a complex combination of intrinsic
fitness effects as well as extrinsic environmental variability.


  The emergence of cross species interactions at protein level is a part of
molecular mechanisms that lead to parasitic diseases. Comprehensive modelling
can capture such interactions and could be useful to understand their
pathophysiology and assist in identifying novel drug targets. Using combination
of databases, text minig and predictive methods, we expanded the sparse
information space of protein-protein interactions in three parasitic diseases
namely, malaria, sleeping sickness and cattle east coast fever. These network
models revealed significant similarities in molecular mechanisms underlaying
host's invasion, immuno-modulation and energy metabolism. The models also
suggested new possible pathways in the inter-species protein interaction maps.
Enrichment of these maps with drug-target informations showed a plethora of
drug space to be explored and led to proposal of two new targets for each
malaria and sleeping sickness.


  Biochemical oscillation systems often consist of negative feedback loops with
repressive transcription regulation. Such systems have distinctive
characteristics in comparison with ordinary chemical systems: i) the numbers of
molecules involved are small, ii) there are typically only a couple of genes in
a cell with a finite regulation time scale. Due to the fluctuations caused by
these features, the system behavior can be quite different from the one
obtained by rate equations, because the rate equations ignore molecular
fluctuations and thus are exact only in the infinite molecular number limit.
The molecular fluctuations on a free-running circadian system have been studied
by Gonze et al. (2002) by introducing a scale parameter $\Omega$ for the system
size. They consider, however, only the first effect, assuming that the gene
process is fast enough for the second effect to be ignored, but this has not
been examined systematically yet. In this work, we study fluctuation effects
due to the finite gene regulation time by introducing a new scale parameter
$\tau$, which we take as the unbinding time of a nuclear protein from the gene.
We focus on the case where the fluctuations due to small molecular numbers can
be ignored. In simulations on the same system studied by Gonze et al., we find
the system is sensitive to the fluctuation in the transcription regulation; the
period of oscillation fluctuates about 30 min even when the regulation time
scale $\tau$ is around 30 s, that is even smaller than 1/1000 of its circadian
period. We also demonstrate that the distribution width for the oscillation
period and the amplitude scales with $\sqrt\tau$, and the correlation time of
the oscillation scales with $1/\tau$ in the small $\tau$ regime. The relative
fluctuations for the period are about half of that for the amplitude, namely,
the periodicity is more stable than the amplitude.


  Background: Biological effects of extra-low-frequency (ELF) magnetic fields
(MF) have lacked a credible mechanism of interaction between MFs and living
material. Objectives: Examine the effect of ELF-MFs on cancer cells. Methods:
Five cancer cell lines were exposed to ELF-MFs within the range of 0.025 to 5
microT, and the cells were examined for karyotype changes after 6 days.
Results: All cancer cells lines lost chromosomes from MF exposure, with a
mostly flat dose-response. Constant MF exposures for three weeks allow a rising
return to the baseline, unperturbed karyotypes. From this point, small MF
increases or decreases are again capable of inducing karyotype contractions.
Our data suggests that the karyotype contractions are caused by MF interference
with mitochondria's ATP synthase (ATPS), compensated by the action of
AMP-activated Protein Kinase (AMPK). The effects of MFs are similar to those of
the ATPS inhibitor oligomycin. They are amplified by metformin, an AMPK
stimulator, and attenuated by resistin, an AMPK inhibitor. Over environmental
MFs, karyotype contractions of various cancer cell lines show exceptionally
wide and flat dose-responses, except for those of erythro-leukemia cells, which
display a progressive rise from 0.025 to 0.4 microT. Conclusions: The
biological effects of MFs are connected to an alteration in the structure of
water that impedes the flux of protons in ATPS channels. These results may be
environmentally important, in view of the central roles played in human
physiology by ATPS and AMPK, particularly in their links to diabetes, cancer
and longevity.


  Proteins encoded by genes containing regions of variable number tandem
repeats (VNTRs) are known to be polymorphic within species but the influence of
their instability in molecular interactions remains unclear. VNTRs are
overrepresented in encoding sequence of particular functional groups where
their presence could influence protein interactions. Using human consensus
coding sequence, this work examines if genomic instability, determined by
regions of VNTRs, influences the number of protein interactions. Findings
reveal that, in relation to protein function, the frequency of unique
interactions in human proteins increase with the number of repeated regions.
This supports experimental evidence that repeat expansion may lead to an
increase in molecular interactions. Genetic diversity, estimated by Ka/Ks,
appeared to decrease as the number of protein-protein interactions increased.
Additionally, G+C and CpG content were negatively correlated with increasing
occurrence of VNTRs. This may indicate that nucleotide composition along with
selective processes can increase genomic stability and thereby restrict the
expansion of repeated regions. Proteins involved in acetylation are associated
with a high number of repeated regions and interactions but a low G+C and CpG
content. While in contrast, less interactive membrane proteins contain a lower
number of repeated regions but higher levels of C+G and CpGs. This work
provides further evidence that VNTRs may provide the genetic variability to
generate unique interactions between proteins.


  Background: Microbiota have evolved to acclimate themselves to many
environments. Humanity is become ever increasingly medicated and many of those
medications are antibiotics. Sadly, Microbiota are adapting to medication and
with each passing generation they become more difficult to subdue. The 16S
small subunit of bacterial ribosomal rRNA provides a wealth of information for
classifying the species level taxonomy of bacteria. Methodology/Principal
Findings: Experiments were collected utilizing broad and narrow spectrum
antibiotics, which act primarily on DNA. In each experiment a statistically
significant, unique and predictable pattern of sequential and thermodynamic
stability or instability was found to correlate to antibiotic resistance.
Conclusions/Significance: Classification of antibiotic resistance is possible
for some species and antibiotic combinations using the 16S rRNA sequential and
thermodynamic properties.


  We use a recently developed coarse-grained model to simulate the
overstretching of duplex DNA. Overstretching at 23C occurs at 74 pN in the
model, about 6-7 pN higher than the experimental value at equivalent salt
conditions. Furthermore, the model reproduces the temperature dependence of the
overstretching force well. The mechanism of overstretching is always
force-induced melting by unpeeling from the free ends. That we never see S-DNA
(overstretched duplex DNA), even though there is clear experimental evidence
for this mode of overstretching under certain conditions, suggests that S-DNA
is not simply an unstacked but hydrogen-bonded duplex, but instead probably has
a more exotic structure.


  With a simple phenomenological metapopulation model, which characterizes the
invasion process of an influenza pandemic from a source to a subpopulation at
risk, we compare the efficiency of inter- and intra-population interventions in
delaying the arrival of an influenza pandemic. We take travel restriction and
patient isolation as examples, since in reality they are typical control
measures implemented at the inter- and intra-population levels, respectively.
We find that the intra-population interventions, e.g., patient isolation,
perform better than the inter-population strategies such as travel restriction
if the response time is small. However, intra-population strategies are
sensitive to the increase of the response time, which might be inevitable due
to socioeconomic reasons in practice and will largely discount the efficiency.


  We examine the distribution of heterozygous sites in nine European and nine
Yoruban individuals whose genomic sequences were made publicly available by
Complete Genomics. We show that it is possible to obtain detailed information
about inbreeding when a relatively small set of whole-genome sequences is
available. Rather than focus on testing for deviations from Hardy-Weinberg
genotype frequencies at each site, we analyze the entire distribution of
heterozygotes conditioned on the number of copies of the derived
(non-chimpanzee) allele. Using Levene's exact test, we reject Hardy-Weinberg in
both populations. We generalized Levene's distribution to obtain the exact
distribution of the number of heterozygous individuals given that every
individual has the same inbreeding coefficient, F. We estimated F to be 0.0026
in Europeans and 0.0005 in Yorubans, but we could also reject the hypothesis
that F was the same in each individual. We used a composite likelihood method
to estimate F in each individual and within each chromosome. Variation in F
across chromosomes within individuals was too large to be consistent with
sampling effects alone. Furthermore, estimates of F for each chromosome in
different populations were not correlated. Our results show how detailed
comparisons of population genomic data can be made to theoretical predictions.
The application of methods to the Complete Genomics data set shows that the
extent of apparent inbreeding varies across chromosomes and across individuals,
and estimates of inbreeding coefficients are subject to unexpected levels of
variation which might be partly accounted for by selection.


  Nested canalizing Boolean (NCF) functions play an important role in
biological motivated regulative networks and in signal processing, in
particular describing stack filters. It has been conjectured that NCFs have a
stabilizing effect on the network dynamics. It is well known that the average
sensitivity plays a central role for the stability of (random) Boolean
networks. Here we provide a tight upper bound on the average sensitivity for
NCFs as a function of the number of relevant input variables. As conjectured in
literature this bound is smaller than 4/3 This shows that a large number of
functions appearing in biological networks belong to a class that has very low
average sensitivity, which is even close to a tight lower bound.


  Sensory hair cells in auditory and vestibular organs rely on active
mechanisms to achieve high sensitivity and frequency selectivity. Recent
experimental studies have documented self-sustained oscillations in hair cells
of lower vertebrates on two distinct levels. First, the hair bundle can undergo
spontaneous mechanical oscillations. Second, somatic electric voltage
oscillations across the baso-lateral membrane of the hair cell have been
observed. We develop a biophysical model of the bullfrog's saccular hair cell
consisting of two compartments, mechanical and electrical, to study how the
mechanical and the voltage oscillations interact to produce coherent
self-sustained oscillations and how this interaction contributes to the overall
sensitivity and selectivity of the hair cell. The model incorporates nonlinear
mechanical stochastic hair bundle system coupled bi-directionally to a
Hodgkin-Huxley type system describing somatic ionic currents. We isolate
regions of coherent spontaneous oscillations in the parameter space of the
model and then study how coupling between compartments affect sensitivity of
the hair cell to external mechanical perturbations. We show that spontaneous
electrical oscillations may enhance significantly the sensitivity and
selectivity of the hair cell.


  Apoptosis is one of the most basic biological processes. In apoptosis, tens
of species are involved in many biochemical reactions with times scales of
widely differing orders of magnitude. By the law of mass action, the process is
mathematically described with a large and stiff system of ODEs (ordinary
differential equations). The goal of this work is to simplify such systems of
ODEs with the PEA (partial equilibrium approximation) method. In doing so, we
propose a general framework of the PEA method together with some conditions,
under which the PEA method can be justified rigorously. The main condition is
the principle of detailed balance for fast reactions as a whole. With the
justified method as a tool, we made many attempts via numerical tests to
simplify the Fas-signaling pathway model due to Hua et al. (2005) and found
that nine of reactions therein can be well regarded as relatively fast. This
paper reports our simplification of Hua at el.'s model with the PEA method
based on the fastness of the nine reactions, together with numerical results
which confirm the reliability of our simplified model.


  Mass action systems capture chemical reaction networks in homogeneous and
dilute solutions. We suggest a notion of generalized mass action systems that
admits arbitrary nonnegative power-law rate functions and serves as a more
realistic model for reaction networks in intracellular environments. In
addition to the chemical complexes and the related stoichiometric subspace, we
introduce corresponding kinetic complexes, which represent the nonnegative
exponents in the rate functions and determine the kinetic-order subspace. We
show that several results of Chemical Reaction Network Theory carry over to the
case of generalized mass action kinetics. Our main result essentially states
that, if the sign vectors of the stoichiometric and the kinetic-order subspace
coincide, there exists a unique positive complex balancing equilibrium in every
stoichiometric compatibility class. However, in contrast to classical mass
action systems, multiple complex balancing equilibria in one stoichiometric
compatibility class are possible in general.


  The H-NS chromosome-organizing protein in E. coli can stabilize genomic DNA
loops, and form oligomeric structures connected to repression of gene
expression. Motivated by the link between chromosome organization, protein
binding and gene expression, we analyzed publicly available genomic data sets
of various origins, from genome-wide protein binding profiles to evolutionary
information, exploring the connections between chromosomal organization,
genesilencing, pseudo-gene localization and horizontal gene transfer. We report
the existence of transcriptionally silent contiguous areas corresponding to
large regions of H-NS protein binding along the genome, their position
indicates a possible relationship with the known large-scale features of
chromosome organization.


  A population has two types of individuals, each occupying an island. One of
those, where individuals of type 1 live, offers a variable environment. Type 2
individuals dwell on the other island, in a constant environment. Only one-way
migration (1->2) is possible. We study the asymptotics of the survival
probability in critical and subcritical cases.


  We investigate extinction dynamics in the paradigmatic model of two competing
species A and B that reproduce (A-->2A, B-->2B), self-regulate by annihilation
(2A-->0, 2B-->0), and compete (A+B-->A, A+B-->B). For a finite system that is
in the well-mixed limit, a quasi-stationary state arises which describes
coexistence of the two species. Because of discrete noise, both species
eventually become extinct in time that is exponentially long in the
quasi-stationary population size. For a sizable range of asymmetries in the
growth and competition rates, the paradoxical situation arises in which the
numerically disadvantaged species according to the deterministic rate equations
survives much longer.


  The Gene Ontology (GO) provides biologists with a controlled terminology that
describes how genes are associated with functions and how functional terms are
related to each other. These term-term relationships encode how scientists
conceive the organization of biological functions, and they take the form of a
directed acyclic graph (DAG). Here, we propose that the network structure of
gene-term annotations made using GO can be employed to establish an alternate
natural way to group the functional terms which is different from the
hierarchical structure established in the GO DAG. Instead of relying on an
externally defined organization for biological functions, our method connects
biological functions together if they are performed by the same genes, as
indicated in a compendium of gene annotation data from numerous different
experiments. We show that grouping terms by this alternate scheme is distinct
from term relationships defined in the ontological structure and provides a new
framework with which to describe and predict the functions of experimentally
identified sets of genes.


  While the concepts involved in Self-Organized Criticality have stimulated
thousands of theoretical models, only recently have these models addressed
problems of biological and clinical importance. Here we outline how SOC can be
used to engineer hybrid viral proteins whose properties, extrapolated from
those of known strains, may be sufficiently effective to cure cancer.


  Martincorena et al. estimated synonymous diversity ({\theta}s = 2N{\mu})
across 2,930 orthologous gene alignments from 34 Escherichia coli genomes, and
found substantial variation among genes in the density of synonymous
polymorphisms. They argue that this pattern reflects variation in the mutation
rate per nucleotide ({\mu}) among genes. However, the effective population size
(N) is not necessarily constant across the genome. In particular, different
genes may have different histories of horizontal gene transfer (HGT), whereas
Martincorena et al. used a model with random recombination to calculate
{\theta}s. They did filter alignments in an effort to minimize the effects of
HGT, but we doubt that any procedure can completely eliminate HGT among closely
related genomes, such as E. coli living in the complex gut community.
  Here we show that there is no significant variation among genes in rates of
synonymous substitutions in a long-term evolution experiment with E. coli and
that the per-gene rates are not correlated with {\theta}s estimates from genome
comparisons. However, there is a significant association between {\theta}s and
HGT events. Together, these findings imply that {\theta}s variation reflects
different histories of HGT, not local optimization of mutation rates to reduce
the risk of deleterious mutations as proposed by Martincorena et al.


  Over the past decade, advances in super-resolution microscopy and
particle-based modeling have driven an intense interest in investigating
spatial heterogeneity at the level of single molecules in cells. Remarkably, it
is becoming clear that spatiotemporal correlations between just a few molecules
can have profound effects on the signaling behavior of the entire cell. While
such correlations are often explicitly imposed by molecular structures such as
rafts, clusters, or scaffolds, they also arise intrinsically, due strictly to
the small numbers of molecules involved, the finite speed of diffusion, and the
effects of macromolecular crowding. In this chapter we review examples of both
explicitly imposed and intrinsic correlations, focusing on the mechanisms by
which microscopic heterogeneity is amplified to macroscopic effect.


  Ray-finned fishes constitute the dominant radiation of vertebrates with over
30,000 species. Although molecular phylogenetics has begun to disentangle major
evolutionary relationships within this vast section of the Tree of Life, there
is no widely available approach for efficiently collecting phylogenomic data
within fishes, leaving much of the enormous potential of massively parallel
sequencing technologies for resolving major radiations in ray-finned fishes
unrealized. Here, we provide a genomic perspective on longstanding questions
regarding the diversification of major groups of ray-finned fishes through
targeted enrichment of ultraconserved nuclear DNA elements (UCEs) and their
flanking sequence. Our workflow efficiently and economically generates data
sets that are orders of magnitude larger than those produced by traditional
approaches and is well-suited to working with museum specimens. Analysis of the
UCE data set recovers a well-supported phylogeny at both shallow and deep
time-scales that supports a monophyletic relationship between Amia and
Lepisosteus (Holostei) and reveals elopomorphs and then osteoglossomorphs to be
the earliest diverging teleost lineages. Divergence time estimation based upon
14 fossil calibrations reveals that crown teleosts appeared ~270 Ma at the end
of the Permian and that elopomorphs, osteoglossomorphs, ostarioclupeomorphs,
and euteleosts diverged from one another by 205 Ma during the Triassic. Our
approach additionally reveals that sequence capture of UCE regions and their
flanking sequence offers enormous potential for resolving phylogenetic
relationships within ray-finned fishes.


  All cells must keep time to consistently perform vital biological functions.
To that end, the coupling and interrelatedness of diverse subsecond events in
the complex cellular environment, such as protein folding or translation rates,
cannot simply result from the chance convergence of the inherent chemical
properties of these phenomena, but may instead be synchronized through a
cell-wide pacemaking mechanism. Picosecond vibrations of lipid membranes may
play a role in such a mechanism.


  We develop an algorithm based on three basic DNA editing operations suggested
by a model for ciliate micronuclear decryption, to transform a given
permutation into another. The number of ciliate operations performed by our
algorithm during such a transformation is taken to be the distance between two
such permutations. Applying well-known clustering methods to such distance
functions enables one to determine phylogenies among the items to which the
distance functions apply. As an application of these ideas we explore the
relationships among the chromosomes of eight fruitfly (drosophila) species,
using the well-known UPGMA algorithm on the distance function provided by our
algorithm.


  The Cancer Genome Atlas (TCGA) provides researchers with clinicopathological
data and genomic characterizations of various carcinomas. These data sets
include expression microarrays for genes and microRNAs -- short, non-coding
strands of RNA that downregulate gene expression through RNA interference -- as
well as days_to_death and days_to_last_followup fields for each tumor sample.
Our aim is to develop a software tool that screens TCGA data sets for
genes/miRNAs with functional involvement in specific cancers. Furthermore, our
computational pipeline is intended to produce a set of visualizations, or
profiles, that place our screened outputs in a pathway-centric context. We
accomplish our 'screening' by ranking genes/miRNAs by the correlation of their
expression misregulation with differential patient survival. In other words, if
a gene/miRNA is consistently misregulated in patients with poor survival rates
and, on the other hand, is expressed more 'normally' in patients with longer
survival rates, then it is ranked highly; if its misregulation has no such
correlation with good/bad survival in patients, then its rank is low. Our
pathway profiling pipeline produces several outputs, which allow us to examine
the functional roles played by highly ranked genes discovered by our screening.
Running the OV (ovarian serous cystadenocarcinoma) data set through our
analysis pipeline, we find that several highly ranked pathways and functional
groups of genes (VEGF, Jun, Fos, etc.) have already been shown to play some
part in the development of epithelial ovarian carcinomas. We also observe that
the dysfunction of the Wnt signaling pathway, which regulates cell-fate
specification and progenitor cell differentiation, has a disproportionate
impact on the survival of ovarian cancer patients.


  The Jacobian matrix of a dynamic system and its principal minors play a
prominent role in the study of qualitative dynamics and bifurcation analysis.
When interpreting the Jacobian as an adjacency matrix of an interaction graph,
its principal minors correspond to sets of disjoint cycles in this graph and
conditions for various dynamic behaviors can be inferred from its cycle
structure. For chemical reaction systems, more fine-grained analyses are
possible by studying a bipartite species-reaction graph. Several results on
injectivity, multistationarity, and bifurcations of a chemical reaction system
have been derived by using various definitions of such bipartite graph. Here,
we present a new definition of the species-reaction graph that more directly
connects the cycle structure with determinant expansion terms, principal
minors, and the coefficients of the characteristic polynomial and encompasses
previous graph constructions as special cases. This graph has a direct relation
to the interaction graph, and properties of cycles and sub-graphs can be
translated in both directions. A simple equivalence relation enables to
decompose determinant expansions more directly and allows simpler and more
direct proofs of previous results.


  Despite considerable progress in genome- and proteome-based high-throughput
screening methods and in rational drug design, the increase in approved drugs
in the past decade did not match the increase of drug development costs.
Network description and analysis not only give a systems-level understanding of
drug action and disease complexity, but can also help to improve the efficiency
of drug design. We give a comprehensive assessment of the analytical tools of
network topology and dynamics. The state-of-the-art use of chemical similarity,
protein structure, protein-protein interaction, signaling, genetic interaction
and metabolic networks in the discovery of drug targets is summarized. We
propose that network targeting follows two basic strategies. The central hit
strategy selectively targets central nodes/edges of the flexible networks of
infectious agents or cancer cells to kill them. The network influence strategy
works against other diseases, where an efficient reconfiguration of rigid
networks needs to be achieved by targeting the neighbors of central nodes or
edges. It is shown how network techniques can help in the identification of
single-target, edgetic, multi-target and allo-network drug target candidates.
We review the recent boom in network methods helping hit identification, lead
selection optimizing drug efficacy, as well as minimizing side-effects and drug
toxicity. Successful network-based drug development strategies are shown
through the examples of infections, cancer, metabolic diseases,
neurodegenerative diseases and aging. Summarizing more than 1200 references we
suggest an optimized protocol of network-aided drug development, and provide a
list of systems-level hallmarks of drug quality. Finally, we highlight
network-related drug development trends helping to achieve these hallmarks by a
cohesive, global approach.


  From the analysis of sizes of approximately 130 small icosahedral viruses we
find that there is a typical structural capsid protein, having a mean diameter
of 5 nm and a mean thickness of 3 nm, with more than two thirds of the analyzed
capsid proteins having thicknesses between 2 nm and 4 nm. To investigate
whether, in addition to the fairly conserved geometry, capsid proteins show
similarities in the way they interact with one another, we examined the shapes
of the capsids in detail. We classified them numerically according to their
similarity to sphere and icosahedron and an interpolating set of shapes in
between, all of them obtained from the theory of elasticity of shells. In order
to make a unique and straightforward connection between an idealized,
numerically calculated shape of an elastic shell and a capsid, we devised a
special shape fitting procedure, the outcome of which is the idealized elastic
shape fitting the capsid best. Using such a procedure we performed statistical
analysis of a series of virus shapes and we found similarities between the
capsid elastic properties of even very different viruses. As we explain in the
paper, there are both structural and functional reasons for the convergence of
protein sizes and capsid elastic properties. Our work presents a specific
quantitative scheme to estimate relatedness between different proteins based on
the details of the (quaternary) shape they form (capsid). As such, it may
provide an information complementary to the one obtained from the studies of
other types of protein similarity, such as the overall composition of
structural elements, topology of the folded protein backbone, and sequence
similarity.


  Phylogenetic networks generalize evolutionary trees, and are commonly used to
represent evolutionary histories of species that undergo reticulate
evolutionary processes such as hybridization, recombination and lateral gene
transfer. Recently, there has been great interest in trying to develop methods
to construct rooted phylogenetic networks from triplets, that is rooted trees
on three species. However, although triplets determine or encode rooted
phylogenetic trees, they do not in general encode rooted phylogenetic networks,
which is a potential issue for any such method. Motivated by this fact, Huber
and Moulton recently introduced trinets as a natural extension of rooted
triplets to networks. In particular, they showed that level-1 phylogenetic
networks are encoded by their trinets, and also conjectured that all
"recoverable" rooted phylogenetic networks are encoded by their trinets. Here
we prove that recoverable binary level-2 networks and binary tree-child
networks are also encoded by their trinets. To do this we prove two
decomposition theorems based on trinets which hold for all recoverable binary
rooted phylogenetic networks. Our results provide some additional evidence in
support of the conjecture that trinets encode all recoverable rooted
phylogenetic networks, and could also lead to new approaches to construct
phylogenetic networks from trinets.


  "Pull moves" is a popular move set for lattice polymer model simulations. We
show that the proof given for its reversibility earlier is flawed, and some
moves are irreversible, which leads to biases in the parameters estimated from
the simulations. We show how to make the move set fully reversible.


  Cooperation is beneficial for the species as a whole, but, at the level of an
individual, defection pays off. Natural selection is then expected to favor
defectors and eliminate cooperation. This prediction is in stark contrast with
the abundance of cooperation at all levels of biological systems: from cells
cooperating to form a biofilm or an organism to ecosystems and human societies.
Several explanations have been proposed to resolve this paradox, including
direct reciprocity, kin, and group selection. However, our work builds upon an
observation that selection on cooperators might depend both on their relative
frequency in the population and on the population density. We find that this
feedback between the population and evolutionary dynamics can substantially
increase the frequency of cooperators at the front of an expanding population,
and can even lead to a splitting of cooperators from defectors. After
splitting, only cooperators colonize new territories, while defectors slowly
invade them from behind. Since range expansions are very common in nature, our
work provides a new explanation of the maintenance of cooperation.


  The spatial propagation of many livestock infectious diseases critically
depends on the animal movements among premises; so the knowledge of movement
data may help us to detect, manage and control an outbreak. The identification
of robust spreading features of the system is however hampered by the temporal
dimension characterizing population interactions through movements. Traditional
centrality measures do not provide relevant information as results strongly
fluctuate in time and outbreak properties heavily depend on geotemporal initial
conditions. By focusing on the case study of cattle displacements in Italy, we
aim at characterizing livestock epidemics in terms of robust features useful
for planning and control, to deal with temporal fluctuations, sensitivity to
initial conditions and missing information during an outbreak. Through spatial
disease simulations, we detect spreading paths that are stable across different
initial conditions, allowing the clustering of the seeds and reducing the
epidemic variability. Paths also allow us to identify premises, called
sentinels, having a large probability of being infected and providing critical
information on the outbreak origin, as encoded in the clusters. This novel
procedure provides a general framework that can be applied to specific
diseases, for aiding risk assessment analysis and informing the design of
optimal surveillance systems.


  We study by Small Angle Neutron Scattering (SANS) the structure of Hyaluronan
-Lysozyme complexes. Hyaluronan (HA) is a polysaccharide of 9 nm intrinsic
persistence length that bears one negative charge per disaccharide monomer
(Mmol = 401.3 g/mol); two molecular weights, Mw = 6000 and 500 000 Da were
used. The pH was adjusted at 4.7 and 7.4 so that lysozyme has a global charge
of +10 and + 8 respectively. The lysozyme concentration was varied from 3 to 40
g/L, at constant HA concentration (10 g/L). At low protein concentration,
samples are monophasic and SANS experiments reveal only fluctuations of
concentration although, at high protein concentration, clusters are observed by
SANS in the dense phase of the diphasic samples. In between, close to the onset
of the phase separation, a distinct original scattering is observed. It is
characteristic of a rod-like shape, which could characterize "single" complexes
involving one or a few polymer chains. For the large molecular weight (500 000)
the rodlike rigid domains extend to much larger length scale than the
persistence length of the HA chain alone in solution and the range of the SANS
investigation. They can be described as a necklace of proteins attached along a
backbone of diameter one or a few HA chains. For the short chains (Mw ~ 6000),
the rod length of the complexes is close to the chain contour length (~ 15 nm).


  In just the last decade, a multitude of bio-technologies and software
pipelines have emerged to revolutionize genomics. To further their central
goal, they aim to accelerate and improve the quality of de novo whole-genome
assembly starting from short DNA reads. However, the performance of each of
these tools is contingent on the length and quality of the sequencing data, the
structure and complexity of the genome sequence, and the resolution and quality
of long-range information. Furthermore, in the absence of any metric that
captures the most fundamental "features" of a high-quality assembly, there is
no obvious recipe for users to select the most desirable assembler/assembly.
International competitions such as Assemblathons or GAGE tried to identify the
best assembler(s) and their features. Some what circuitously, the only
available approach to gauge de novo assemblies and assemblers relies solely on
the availability of a high-quality fully assembled reference genome sequence.
Still worse, reference-guided evaluations are often both difficult to analyze,
leading to conclusions that are difficult to interpret. In this paper, we
circumvent many of these issues by relying upon a tool, dubbed FRCbam, which is
capable of evaluating de novo assemblies from the read-layouts even when no
reference exists. We extend the FRCurve approach to cases where lay-out
information may have been obscured, as is true in many deBruijn-graph-based
algorithms. As a by-product, FRCurve now expands its applicability to a much
wider class of assemblers -- thus, identifying higher-quality members of this
group, their inter-relations as well as sensitivity to carefully selected
features, with or without the support of a reference sequence or layout for the
reads. The paper concludes by reevaluating several recently conducted assembly
competitions and the datasets that have resulted from them.


  We study a model of species survival recently proposed by Michael and Volkov.
We interpret it as a variant of empirical processes, in which the sample size
is random and when decreasing, samples of smallest numerical values are
removed. Micheal and Volkov proved that the empirical distributions converge to
the sample distribution conditioned not to be below a certain threshold. We
prove a functional central limit theorem for the fluctuations. There exists a
threshold above which the limit process is Gaussian with variance bounded below
by a positive constant, while at the threshold it is half-Gaussian.


  Genomic and proteomic methods yield networks of biological regulatory
interactions but do not provide direct insight into how those interactions are
organized into functional modules, or how information flows from one module to
another. In this work we introduce an approach that provides this complementary
information and apply it to the bacterium Caulobacter crescentus, a paradigm
for cell-cycle control. Operationally, we use an inducible promoter to express
the essential transcriptional regulatory gene ctrA in a periodic, pulsed
fashion. This chemical perturbation causes the population of cells to divide
synchronously, and we use the resulting advance or delay of the division times
of single cells to construct a phase resetting curve. We find that delay is
strongly favored over advance. This finding is surprising since it does not
follow from the temporal expression profile of CtrA and, in turn, simulations
of existing network models. We propose a phenomenological model that suggests
that the cell-cycle network comprises two distinct functional modules that
oscillate autonomously and couple in a highly asymmetric fashion. These
features collectively provide a new mechanism for tight temporal control of the
cell cycle in C. crescentus. We discuss how the procedure can serve as the
basis for a general approach for probing network dynamics, which we term
chemical perturbation spectroscopy (CPS).


  Early warning signals have been proposed to forecast the possibility of a
critical transition, such as the eutrophication of a lake, the collapse of a
coral reef, or the end of a glacial period. Because such transitions often
unfold on temporal and spatial scales that can be difficult to approach by
experimental manipulation, research has often relied on historical observations
as a source of natural experiments. Here we examine a critical difference
between selecting systems for study based on the fact that we have observed a
critical transition and those systems for which we wish to forecast the
approach of a transition. This difference arises by conditionally selecting
systems known to experience a transition of some sort and failing to account
for the bias this introduces -- a statistical error often known as the
Prosecutor's Fallacy. By analysing simulated systems that have experienced
transitions purely by chance, we reveal an elevated rate of false positives in
common warning signal statistics. We further demonstrate a model-based approach
that is less subject to this bias than these more commonly used summary
statistics. We note that experimental studies with replicates avoid this
pitfall entirely.


  The methylation of DNA regulates gene expression. On cell division the
methylation state of the DNA is typically inherited from parent to daughter
cells. While the chemical bond between the methyl group and the DNA is very
strong, changes to the methylation state do occur and are observed to occur
rapidly in response to external stimulus. The loss of methylation can be active
where enzyme physically breaks the bond, or passive where on cell division the
newly constructed strand of DNA is not properly inherited.
  Here we present a mathematical model of single locus passive demethylation
for a dividing population of cells. The model describes the heterogenity in the
population expected from passive mechanisms. We see that even when the site
specific probabilities of passive demethylation are independent, conservation
of methylation on the inherited strand gives rise to site-site correlations of
the methylation state. We then extend the model to incorporate correlations
between sites in the locus for demethylation rates. Biologically, correlations
in demethylation rates might correspond to locus wide changes such as the
inability of methyltransferase to access the locus. We also look at the effects
of selection on the multicellular population.
  The model of passive demethylation not only provides a tool for measurement
of parameters in loci-specific cases where passive demethylation is the
dominant mechanism, but also provides a baseline in the search for active
mechanisms. The model tells us that there are states of methylation
inaccessible by passive mechanisms. Observation of these states constitutes
evidence of active mechanisms, either de novo methylation or enzymatic
demethylation. We also see that selection and passive demethylation combined
can give rise to a stable heterogeneous distribution of gene methylation states
in a population.


  Gene networks exhibiting oscillatory dynamics are widespread in biology. The
minimal regulatory designs giving rise to oscillations have been implemented
synthetically and studied by mathematical modeling. However, most of the
available analyses generally neglect the coupling of regulatory circuits with
the cellular "chassis" in which the circuits are embedded. For example, the
intracellular macromolecular composition of fast-growing bacteria changes with
growth rate. As a consequence, important parameters of gene expression, such as
ribosome concentration or cell volume, are growth-rate dependent, ultimately
coupling the dynamics of genetic circuits with cell physiology. This work
addresses the effects of growth rate on the dynamics of a paradigmatic example
of genetic oscillator, the repressilator. Making use of empirical growth-rate
dependences of parameters in bacteria, we show that the repressilator dynamics
can switch between oscillations and convergence to a fixed point depending on
the cellular state of growth, and thus on the nutrients it is fed. The physical
support of the circuit (type of plasmid or gene positions on the chromosome)
also plays an important role in determining the oscillation stability and the
growth-rate dependence of period and amplitude. This analysis has potential
application in the field of synthetic biology, and suggests that the coupling
between endogenous genetic oscillators and cell physiology can have substantial
consequences for their functionality.


  We review modeling attempts for the paradigmatic contact process (or SIS
model) on adaptive networks. Elaborating on one particular proposed mechanism
of topology change (rewiring) and its mean field analysis, we obtain a
coarse-grained view of coevolving network topology in the stationary active
phase of the system. Introducing an alternative framework applicable to a wide
class of adaptive networks, active stationary states are detected, and an
extended description of the resulting steady-state statistics is given for
three different rewiring schemes. We find that slight modifications of the
standard rewiring rule can result in either minuscule or drastic change of
steady-state network topologies.


  The Escalator Boxcar Train (EBT) is a numerical method that is widely used in
theoretical biology to investigate the dynamics of physiologically structured
population models, i.e., models in which individuals differ by size or other
physiological characteristics. The method was developed more than two decades
ago, but has so far resisted attempts to give a formal proof of convergence.
Using a modern framework of measure-valued solutions, we investigate the EBT
method and show that the sequence of approximating solution measures generated
by the EBT method converges weakly to the true solution measure under weak
conditions on the growth rate, birth rate, and mortality rate. In rigorously
establishing the convergence of the EBT method, our results pave the way for
wider acceptance of the EBT method beyond theoretical biology and constitutes
an important step towards integration with established numerical schemes.


  We introduce the epidemic quasimetric on graphs and study its behavior with
respect to clustering techniques. In particular we compare its behavior to
known objects such as the graph distance, effective resistance, and modulus of
path families.


  Summary: Biospectrogam is an open-source software for the spectral analysis
of DNA and protein sequences. The software can fetch (from NCBI server), import
and manage biological data. One can analyze the data using Digital Signal
Processing (DSP) techniques since the software allows the user to convert the
symbolic data into numerical data using 23 popular encodings and then apply
popular transformations such as Fast Fourier Transform (FFT) etc. and export
it. The ability of exporting (both encoding files and transform files) as a
MATLAB .m file gives the user an option to apply variety of techniques of DSP.
User can also do window analysis (both sliding in forward and backward
directions and stagnant) with different size windows and search for meaningful
spectral pattern with the help of exported MATLAB file in a dynamic manner by
choosing time delay in the plot using Biospectrogram. Random encodings and user
choice encoding allows software to search for many possibilities in spectral
space.
  Availability: Biospectrogam is written in Java and is available to download
freely from http://www.guptalab.org/biospectrogram. Software has been optimized
to run on Windows, Mac OSX and Linux. User manual and you-tube (product demo)
tutorial is also available on the website. We are in the process of acquiring
open source license for it.


  Evolutionary relationships among birds in Neoaves, the clade comprising the
vast majority of avian diversity, have vexed systematists due to the ancient,
rapid radiation of numerous lineages. We applied a new phylogenomic approach to
resolve relationships in Neoaves using target enrichment (sequence capture) and
high-throughput sequencing of ultraconserved elements (UCEs) in avian genomes.
We collected sequence data from UCE loci for 32 members of Neoaves and one
outgroup (chicken) and analyzed data sets that differed in their amount of
missing data. An alignment of 1,541 loci that allowed missing data was 87%
complete and resulted in a highly resolved phylogeny with broad agreement
between the Bayesian and maximum-likelihood (ML) trees. Although results from
the 100% complete matrix of 416 UCE loci were similar, the Bayesian and ML
trees differed to a greater extent in this analysis, suggesting that increasing
from 416 to 1,541 loci led to increased stability and resolution of the tree.
Novel results of our study include surprisingly close relationships between
phenotypically divergent bird families, such as tropicbirds (Phaethontidae) and
the sunbittern (Eurypygidae) as well as between bustards (Otididae) and turacos
(Musophagidae). This phylogeny bolsters support for monophyletic waterbird and
landbird clades and also strongly supports controversial results from previous
studies, including the sister relationship between passerines and parrots and
the non-monophyly of raptorial birds in the hawk and falcon families. Although
significant challenges remain to fully resolving some of the deep relationships
in Neoaves, especially among lineages outside the waterbirds and landbirds,
this study suggests that increased data will yield an increasingly resolved
avian phylogeny.


  In both population genetics and forensic genetics it is important to know how
haplotypes are distributed in a population. Simulation of population dynamics
helps facilitating research on the distribution of haplotypes. In forensic
genetics, the haplotypes can for example consist of lineage markers such as
short tandem repeat loci on the Y chromosome (Y-STR). A dominating model for
describing population dynamics is the simple, yet powerful, Fisher-Wright
model. We describe an efficient algorithm for exact forward simulation of exact
Fisher-Wright populations (and not approximative such as the coalescent model).
The efficiency comes from convenient data structures by changing the
traditional view from individuals to haplotypes. The algorithm is implemented
in the open-source R package 'fwsim' and is able to simulate very large
populations. We focus on a haploid model and assume stochastic population size
with flexible growth specification, no selection, a neutral single step
mutation process, and self-reproducing individuals. These assumptions make the
algorithm ideal for studying lineage markers such as Y-STR.


  Predictability of undesired events is a question of great interest in many
scientific disciplines including seismology, economy, and epidemiology. Here,
we focus on the predictability of invasion of a broad class of epidemics caused
by diseases that lead to permanent immunity of infected hosts after recovery or
death. We approach the problem from the perspective of the science of
complexity by proposing and testing several strategies for the estimation of
important characteristics of epidemics, such as the probability of invasion.
Our results suggest that parsimonious approximate methodologies may lead to the
most reliable and robust predictions. The proposed methodologies are first
applied to analysis of experimentally observed epidemics: invasion of the
fungal plant pathogen \emph{Rhizoctonia solani} in replicated host microcosms.
We then consider numerical experiments of the SIR
(susceptible-infected-removed) model to investigate the performance of the
proposed methods in further detail. The suggested framework can be used as a
valuable tool for quick assessment of epidemic threat at the stage when
epidemics only start developing. Moreover, our work amplifies the significance
of the small-scale and finite-time microcosm realizations of epidemics
revealing their predictive power.


  Edwards et al question aspects of the methods used in two of our published
papers that report results showing Levy walk like and Levy flight movement
patterns of marine predators.The criticisms are focused on the applicability of
some statistical methodologies used to detect power law distributions.We reply
to the principal criticisms levelled at each of these papers in turn including
our own reanalysis of specific datasets and find that neither of our papers
conclusions are overturned in any part by the issues raised.Indeed, in addition
to the findings of our research reported in these papers there is strong
evidence accumulating from studies worldwide that organisms show movements and
behaviour consistent with scale invariant patterns such as Levy flights.


  Antagonistic coevolution between hosts and parasites can have a major impact
on host population structures, and hence on the evolution of social traits.
Using stochastic modelling techniques in the context of bacteria-virus
interactions, we investigate the impact of coevolution across a continuum of
host-parasite genetic specificity (specifically, where genotypes have the same
infectivity/resistance ranges (matching alleles, MA) to highly variable ranges
(gene-for-gene, GFG)) on population genetic structure, and on the social
behaviour of the host. We find that host cooperation is more likely to be
maintained towards the MA end of the continuum, as the more frequent
bottlenecks associated with an MA-like interaction can prevent defector
invasion, and can even allow migrant cooperators to invade populations of
defectors.


  MicroRNAs (miRNAs) are small RNA molecules, about 22 nucleotide long, which
post-transcriptionally regulate their target messenger RNAs (mRNAs). They
accomplish key roles in gene regulatory networks, ranging from signaling
pathways to tissue morphogenesis, and their aberrant behavior is often
associated with the development of various diseases. Recently it has been shown
that, in analogy with the better understood case of small RNAs in bacteria, the
way miRNAs interact with their targets can be described in terms of a titration
mechanism characterized by threshold effects, hypersensitivity of the system
near the threshold, and prioritized cross-talk among targets. The latter
characteristic has been lately identified as competing endogenous RNA (ceRNA)
effect to mark those indirect interactions among targets of a common pool of
miRNAs they are in competition for. Here we analyze the equilibrium and
out-of-equilibrium properties of a general stochastic model of $M$ miRNAs
interacting with $N$ mRNA targets. In particular we are able to describe in
details the peculiar equilibrium and non-equilibrium phenomena that the system
displays around the threshold: (i) maximal cross-talk and correlation between
targets, (ii) robustness of ceRNA effect with respect to the model's parameters
and in particular to the catalyticity of the miRNA-mRNA interaction, and (iii)
anomalous response-time to external perturbations.


  It has recently been suggested that the competition for a finite pool of
microRNAs (miRNA) gives rise to effective interactions among their common
targets (competing endogenous RNAs or ceRNAs) that could prove to be crucial
for post-transcriptional regulation (PTR). We have studied a minimal model of
PTR where the emergence and the nature of such interactions can be
characterized in detail at steady state. Sensitivity analysis shows that
binding free energies and repression mechanisms are the key ingredients for the
cross-talk between ceRNAs to arise. Interactions emerge in specific ranges of
repression values, can be symmetrical (one ceRNA influences another and
vice-versa) or asymmetrical (one ceRNA influences another but not the reverse)
and may be highly selective, while possibly limited by noise. In addition, we
show that non-trivial correlations among ceRNAs can emerge in experimental
readouts due to transcriptional fluctuations even in absence of miRNA-mediated
cross-talk.


  Comprehensive discovery of structural variation (SV) in human genomes from
DNA sequencing requires the integration of multiple alignment signals including
read-pair, split-read and read-depth. However, owing to inherent technical
challenges, most existing SV discovery approaches utilize only one signal and
consequently suffer from reduced sensitivity, especially at low sequence
coverage and for smaller SVs. We present a novel and extremely flexible
probabilistic SV discovery framework that is capable of integrating any number
of SV detection signals including those generated from read alignments or prior
evidence. We demonstrate improved sensitivity over extant methods by combining
paired-end and split-read alignments and emphasize the utility of our framework
for comprehensive studies of structural variation in heterogeneous tumor
genomes. We further discuss the broader utility of this approach for
probabilistic integration of diverse genomic interval datasets.


  High-throughput pooled resequencing offers significant potential for whole
genome population sequencing. However, its main drawback is the loss of
haplotype information. In order to regain some of this information, we present
LDx, a computational tool for estimating linkage disequilibrium (LD) from
pooled resequencing data. LDx uses an approximate maximum likelihood approach
to estimate LD (r2) between pairs of SNPs that can be observed within and among
single reads. LDx also reports r2 estimates derived solely from observed
genotype counts. We demonstrate that the LDx estimates are highly correlated
with r2 estimated from individually resequenced strains. We discuss the
performance of LDx using more stringent quality conditions and infer via
simulation the degree to which performance can improve based on read depth.
Finally we demonstrate two possible uses of LDx with real and simulated pooled
resequencing data. First, we use LDx to infer genomewide patterns of decay of
LD with physical distance in D. melanogaster population resequencing data.
Second, we demonstrate that r2 estimates from LDx are capable of distinguishing
alternative demographic models representing plausible demographic histories of
D. melanogaster.


  Accurately determining and classifying the structure of complex networks is
the focus of much current research. One class of network of particular interest
are metabolic pathways, which have previously been studied from a graph
theoretical viewpoint in a number of ways. Metabolic networks describe the
chemical reactions within cells and are thus of prime importance from a
biological perspective.
  Here we analyse metabolic networks from a section of microorganisms, using a
range of metrics and attempt to address anomalies between the observed metrics
and current descriptions of the graphical structure. We propose that the growth
of the network may in some way be regulated by network size and attempt to
reproduce networks with similar metrics to the metabolic pathways using a
generative approach. We provide some hypotheses as to why biological networks
may evolve according to these model criteria.


  In this paper, by using a stochastic reaction-diffusion-taxis model, we
analyze the picophytoplankton dynamics in the basin of the Mediterranean Sea,
characterized by poorly mixed waters. The model includes intraspecific
competition of picophytoplankton for light and nutrients. The multiplicative
noise sources present in the model account for random fluctuations of
environmental variables. Phytoplankton distributions obtained from the model
show a good agreement with experimental data sampled in two different sites of
the Sicily Channel. The results could be extended to analyze data collected in
different sites of the Mediterranean Sea and to devise predictive models for
phytoplankton dynamics in oligotrophic waters.


  We present a new method for inferring species trees from multi-copy gene
trees. Our method is based on a generalization of the Robinson-Foulds (RF)
distance to multi-labeled trees (mul-trees), i.e., gene trees in which multiple
leaves can have the same label. Unlike most previous phylogenetic methods using
gene trees, this method does not assume that gene tree incongruence is caused
by a single, specific biological process, such as gene duplication and loss,
deep coalescence, or lateral gene transfer. We prove that it is NP-hard to
compute the RF distance between two mul-trees, but it is easy to calculate the
generalized RF distance between a mul-tree and a singly-labeled tree. Motivated
by this observation, we formulate the RF supertree problem for mul-trees
(MulRF), which takes a collection of mul-trees and constructs a species tree
that minimizes the total RF distance from the input mul-trees. We present a
fast heuristic algorithm for the MulRF supertree problem. Simulation
experiments demonstrate that the MulRF method produces more accurate species
trees than gene tree parsimony methods when incongruence is caused by gene tree
error, duplications and losses, and/or lateral gene transfer. Furthermore, the
MulRF heuristic runs quickly on data sets containing hundreds of trees with up
to a hundred taxa.


  We investigate kinetic pathways of the DNA melting transition using
variable-range versions of the Poland-Scheraga (PS) and Peyrard-Dauxois-Bishop
(PDB) models of DNA. In the PS model, we construct a phi^4-field theory to
calculate the critical droplet profile, the initial growth modes, and the
exponent characterizing the divergence of the susceptibility near the spinodal.
In the PDB model, we use a mean field analysis to calculate susceptibility
exponent. We compare these theoretical results with Monte Carlo and Brownian
dynamic simulations on the PS and PDB models, respectively. We find that by
increasing the range of interaction, the system can be brought close to a
pseudospinodal, and that in this region the nucleating droplet is diffuse in
contrast to the compact droplets predicted by classical nucleation theory.


  RNA-Seq technology allows for studying the transcriptional state of the cell
at an unprecedented level of detail. Beyond quantification of whole-gene
expression, it is now possible to disentangle the abundance of individual
alternatively spliced transcript isoforms of a gene. A central question is to
understand the regulatory processes that lead to differences in relative
abundance variation due to external and genetic factors. Here, we present a
mixed model approach that allows for (i) joint analysis and genetic mapping of
multiple transcript isoforms and (ii) mapping of isoform-specific effects.
Central to our approach is to comprehensively model the causes of variation and
correlation between transcript isoforms, including the genomic background and
technical quantification uncertainty. As a result, our method allows to
accurately test for shared as well as transcript-specific genetic regulation of
transcript isoforms and achieves substantially improved calibration of these
statistical tests. Experiments on genotype and RNA-Seq data from 126 human
HapMap individuals demonstrate that our model can help to obtain a more
fine-grained picture of the genetic basis of gene expression variation.


  Complex networks have been shown to be robust against random structural
perturbations, but vulnerable against targeted attacks. Robustness analysis
usually simulates the removal of individual or sets of nodes, followed by the
assessment of the inflicted damage. For complex metabolic networks, it has been
suggested that evolutionary pressure may favor robustness against reaction
removal. However, the removal of a reaction and its impact on the network may
as well be interpreted as selective regulation of pathway activities,
suggesting a tradeoff between the efficiency of regulation and vulnerability.
Here, we employ a cascading failure algorithm to simulate the removal of single
and pairs of reactions from the metabolic networks of two organisms, and
estimate the significance of the results using two different null models:
degree preserving and mass-balanced randomization. Our analysis suggests that
evolutionary pressure promotes larger cascades of non-viable reactions, and
thus favors the ability of efficient metabolic regulation at the expense of
robustness.


  Many biological, physical, and social interactions have a particular
dependence on where they take place. In living cells, protein movement between
the nucleus and cytoplasm affects cellular response (i.e., proteins must be
present in the nucleus to regulate their target genes). Here we use recent
developments from dynamical systems and chemical reaction network theory to
identify and characterize the key-role of the spatial organization of
eukaryotic cells in cellular information processing. In particular the
existence of distinct compartments plays a pivotal role in whether a system is
capable of multistationarity (multiple response states), and is thus directly
linked to the amount of information that the signaling molecules can represent
in the nucleus. Multistationarity provides a mechanism for switching between
different response states in cell signaling systems and enables multiple
outcomes for cellular-decision making. We find that introducing species
localization can alter the capacity for multistationarity and mathematically
demonstrate that shuttling confers flexibility for and greater control of the
emergence of an all-or-none response.


  In spite of the abundance and commercial importance of these two species,
there is little published weight-length data for the gafftopsail catfish (Bagre
marinus) and hardhead catfish (Ariopsis felis). For this study 84 catfish were
caught (hook and line) from the Calcasieu Estuary in Southwest Louisiana near
the Gulf of Mexico and estuaries and near shore waters close to bayou
Lafourche. Using least squares regression, best fit curves were determined for
weight (W) vs. total length (L) relationships in gafftopsail catfish, W(L) =
1000(L/484.73)^3.2660, while the best-fit equation for the hardhead was W(L) =
1000(L/469.53)^3.0188, where W is the weight in grams and L is the total length
in millimeters. Results showed that, when compared to gafftopsail and hardhead
catfish caught in Florida, Louisiana gafftopsail catfish tend to weigh less at
similar lengths; whereas, the Louisiana hardhead catfish tend to weigh about
the same. Results also show that the total length (TL) can be related to the
fork length (FL) as TL = 1.184 FL in gafftopsail catfish with R^2 = 0.9922 and
TL = 1.140 FL with R^2 = 0.9667 in hardhead catfish.


  When modeling coexpression networks from high-throughput time course data,
Pearson Correlation Coefficient (PCC) is one of the most effective and popular
similarity functions. However, its reliability is limited since it cannot
capture non-linear interactions and time shifts. Here we propose to overcome
these two issues by employing a novel similarity function, Dynamic Time Warping
Maximal Information Coefficient (DTW-MIC), combining a measure taking care of
functional interactions of signals (MIC) and a measure identifying horizontal
displacements (DTW). By using the Hamming-Ipsen-Mikhailov (HIM) metric to
quantify network differences, the effectiveness of the DTW-MIC approach is
demonstrated on both synthetic and transcriptomic datasets.


  HIV integrase is a 32 kDa protein produced from the C-terminal portion of the
Pol gene product, and is an attractive target for new anti-HIV drugs. Integrase
is an enzyme produced by a retrovirus (such as HIV) that enables its genetic
material to be integrated into the DNA of the infected cell. Raltegravir and
Elvitegravir are two important drugs against integrase.


  We analyze the dynamics of structural transitions between normal right-handed
B form and unusual left-handed Z form for a linear DNA molecule. The dynamics
under the external torque in physiological buffer is modeled by a Langevin
equation, with the potential term given by the authors previously [Phys. Rev. E
84, 021926 (2011)]. With this model, we first simulate the relaxation processes
around B-form structure after sudden changes of the external torques, where
slow relaxation $\sim t^{-1/2}$ as a function of the elapsed time $t$ is
observed. Then, the dynamics of structural transition from Z form to B form is
computed under various external torque strength. For small external torques,
the transition proceeds via nucleation and the growth, while for higher
torques, Z-form structure becomes unstable, and the transition mechanism is
switched to a spinodal-like process. These numerical results are qualitatively
understood by simple phenomenological arguments.


  We construct an energy function that describes the crystallographic structure
of spermwhale myoglobin backbone. As a model in our construction, we use the
Protein Data Bank entry 1ABS that has been measured at liquid helium
temperature. Consequently the thermal B-factor fluctuations are very small,
which is an advantage in our construction. The energy function that we utilize
resembles that of the discrete non-linear Schrodinger equation. Likewise, ours
supports solitons as local minimum energy configurations. We describe the 1ABS
backbone in terms of solitons with a precision that deviates from 1ABS by an
average root-mean-square distance, which is less than the experimentally
observed Debye-Waller B-factor fluctuation distance. We then subject the
multisoliton solution to extensive numerical heating and cooling experiments,
over a very wide range of temperatures. We concentrate in particular to
temperatures above 300K and below the theta-point unfolding temperature, which
is around 348K. We confirm that the behavior of the multisoliton is fully
consistent with Anfinsen's principle, up to very high temperatures. We observe
that the structure responds to an increase of temperature consistently in a
very similar manner. This enables us to characterize the onset of thermally
induced conformational changes in terms of three distinct backbone ligand
gates. One of the gates is made of the helix F and the helix E. This is a
pathway that is presumed to have a major role in ligand migration between the
heme and the exterior. The two other gates are chosen similarly, when open they
provide a direct access route for a ligand to reach the heme. We find that out
of the three gates we investigate, the one which is formed by helices B and G
is the most sensitive one to thermally induced conformational changes. Our
approach provides a novel perspective to the important problem of ligand
migration.


  Given two rooted phylogenetic trees on the same set of taxa X, the Maximum
Agreement Forest problem (MAF) asks to find a forest that is, in a certain
sense, common to both trees and has a minimum number of components. The Maximum
Acyclic Agreement Forest problem (MAAF) has the additional restriction that the
components of the forest cannot have conflicting ancestral relations in the
input trees. There has been considerable interest in the special cases of these
problems in which the input trees are required to be binary. However, in
practice, phylogenetic trees are rarely binary, due to uncertainty about the
precise order of speciation events. Here, we show that the general, nonbinary
version of MAF has a polynomial-time 4-approximation and a fixed-parameter
tractable (exact) algorithm that runs in O(4^k poly(n)) time, where n = |X| and
k is the number of components of the agreement forest minus one. Moreover, we
show that a c-approximation algorithm for nonbinary MAF and a d-approximation
algorithm for the classical problem Directed Feedback Vertex Set (DFVS) can be
combined to yield a d(c+3)-approximation for nonbinary MAAF. The algorithms for
MAF have been implemented and made publicly available.


  Genetic variants in cis-regulatory elements or trans-acting regulators
commonly influence the quantity and spatiotemporal distribution of gene
transcription. Recent interest in expression quantitative trait locus (eQTL)
mapping has paralleled the adoption of genome-wide association studies (GWAS)
for the analysis of complex traits and disease in humans. Under the hypothesis
that many GWAS associations tag non-coding SNPs with small effects, and that
these SNPs exert phenotypic control by modifying gene expression, it has become
common to interpret GWAS associations using eQTL data. To exploit the
mechanistic interpretability of eQTL-GWAS comparisons, an improved
understanding of the genetic architecture and cell type specificity of eQTLs is
required. We address this need by performing an eQTL analysis in four parts:
first we identified eQTLs from eleven studies on seven cell types; next we
quantified cell type specific eQTLs across the studies; then we integrated eQTL
data with cis-regulatory element (CRE) data sets from the ENCODE project;
finally we built a classifier to predict cell type specific eQTLs. Consistent
with prior studies, we demonstrate that allelic heterogeneity is pervasive at
cis-eQTLs and that cis-eQTLs are often cell type specific. Within and between
cell type eQTL replication is associated with eQTL SNP overlap with hundreds of
cell type specific CRE element classes, including enhancer, promoter, and
repressive chromatin marks, regions of open chromatin, and many classes of DNA
binding proteins. Using a random forest classifier including 526 CRE data sets
as features, we successfully predict the cell type specificity of eQTL SNPs in
the absence of gene expression data from the cell type of interest. We
anticipate that such integrative, predictive modeling will improve our ability
to understand the mechanistic basis of human complex phenotypic variation.


  A problem of substantial interest is to systematically map variation in
chromatin structure to gene expression regulation across conditions,
environments, or differentiated cell types. We developed and applied a
quantitative framework for determining the existence, strength, and type of
relationship between high-resolution chromatin structure in terms of DNaseI
hypersensitivity (DHS) and genome-wide gene expression levels in 20 diverse
human cell lines. We show that ~25% of genes show cell-type specific expression
explained by alterations in chromatin structure. We find that distal regions of
chromatin structure (e.g., +/- 200kb) capture more genes with this relationship
than local regions (e.g., +/- 2.5kb), yet the local regions show a more
pronounced effect. By exploiting variation across cell-types, we were capable
of pinpointing the most likely hypersensitive sites related to cell-type
specific expression, which we show have a range of contextual usages. This
quantitative framework is likely applicable to other settings aimed at relating
continuous genomic measurements to gene expression variation.