  Nowadays, studying introns is becoming a very promising field in the
genomics. Even though they play a role in the dynamic regulation of gene and in
the organism's evolution, introns have not attracted enough attention like
exons did; especially of digital signal processing researchers. Thus, we focus
on analysis of the C.elegans introns. In this paper, we propose the complex
Morlet wavelet analysis to investigate introns' characterization in the
C.elegans genes. However, catching the change in frequency response with
respect to time of the gene sequences is hindered by their presence in the form
of strings of characters. This can only be counteracted by assigning numerical
values to each of the DNA characters. This operation defines the so called "DNA
coding approach". In this context, we propose a new coding technique based on
the Frequency Chaos Game Representation (FCGR) that we name the "Frequency
Chaos Game Signal" (FCGS). Results of the complex Morlet wavelet Analysis
applied to the Celegans FCGS are showing a very distinguished texture. The
visual interpretation of the colour scalograms is proved to be an efficient
tool for revealing significant information about intronic sequences.


  The complex correlation structure of a collection of orthologous DNA
sequences is uniquely captured by the "ancestral recombination graph" (ARG), a
complete record of coalescence and recombination events in the history of the
sample. However, existing methods for ARG inference are computationally
intensive, highly approximate, or limited to small numbers of sequences, and,
as a consequence, explicit ARG inference is rarely used in applied population
genomics. Here, we introduce a new algorithm for ARG inference that is
efficient enough to apply to dozens of complete mammalian genomes. The key idea
of our approach is to sample an ARG of n chromosomes conditional on an ARG of
n-1 chromosomes, an operation we call "threading." Using techniques based on
hidden Markov models, we can perform this threading operation exactly, up to
the assumptions of the sequentially Markov coalescent and a discretization of
time. An extension allows for threading of subtrees instead of individual
sequences. Repeated application of these threading operations results in highly
efficient Markov chain Monte Carlo samplers for ARGs. We have implemented these
methods in a computer program called ARGweaver. Experiments with simulated data
indicate that ARGweaver converges rapidly to the true posterior distribution
and is effective in recovering various features of the ARG for dozens of
sequences generated under realistic parameters for human populations. In
applications of ARGweaver to 54 human genome sequences from Complete Genomics,
we find clear signatures of natural selection, including regions of unusually
ancient ancestry associated with balancing selection and reductions in allele
age in sites under directional selection. Preliminary results also indicate
that our methods can be used to gain insight into complex features of human
population structure, even with a noninformative prior distribution.


  A large, annually recurring region of hypoxia in the northern Gulf of Mexico
has been attributed to water stratification and nutrient loading of nitrogen
and phosphorus delivered by the Mississippi and Atchafalaya rivers. This
nutrient loading increased nearly 300% since 1950, primarily due to increased
use of agricultural fertilizers. Over this same time period, the red snapper
(Lutjanus campechanus) population in the Gulf of Mexico has shifted strongly
from being dominated by the eastern Gulf of Mexico to being dominated by the
northern and western Gulf of Mexico, with the bulk of the current population in
the same regions with significant nutrient loading from the Mississippi and
Atchafalaya rivers and in or near areas with development of mid-summer hypoxic
zones. The population decline of red snapper in the eastern Gulf is almost
certainly attributable to overfishing, but the cause of the population increase
in the northern and western Gulf is subject to broad debate, with the impact of
artificial reefs (primarily oil platforms which have incr eased greatly since
the 1960s) being the most contentious point. Nutrient loading has been shown to
positively impact secondary production of fish in many marine systems. The
present paper offers the hypothesis that increased nutrient loading has
contributed significantly to increased red snapper population in the northern
and western Gulf of Mexico. Nutrient loading may be working in synergy with the
abundant oil platforms both increasing primary production and providing
structure encouraging red snapper to feed throughout the water column.


  The entire T7 bacteriophage genome contains 39937 base pairs (Database NCBI
RefSeq N1001604). Here, electrostatic potential distribution around double
helical T7 DNA was calculated by Coulomb method using the computer program of
Sorokin A.A. Electrostatic profiles of 17 promoters recognized by T7 phage
specific RNA polymerase were analyzed. It was shown that electrostatic profiles
of all T7 RNA polymerase specific promoters can be characterized by distinctive
motifs which are specific for each promoter class. Comparative analysis of
electrostatic profiles of native T7 promoters of different classes demonstrates
that T7 RNA polymerase can differentiate them due to their electrostatic
features.


  Population genetic studies have found evidence for dramatic population growth
in recent human history. It is unclear how this recent population growth,
combined with the effects of negative natural selection, has affected patterns
of deleterious variation, as well as the number, frequencies, and effect sizes
of mutations that contribute risk to complex traits. Here I use simulations
under population genetic models where a proportion of the heritability of the
trait is accounted for by mutations in a subset of the exome. I show that
recent population growth increases the proportion of nonsynonymous variants
segregating in the population, but does not affect the genetic load relative to
that in a population that did not expand. Under a model where a mutation's
effect on a trait is correlated with its effect on fitness, rare variants
explain a greater portion of the additive genetic variance of the trait in a
population that has recently expanded than in a population that did not
recently expand. Further, when using a single-marker test, for a given
false-positive rate and sample size, recent population growth decreases the
expected number of significant association with the trait relative to the
number detected in a population that did not expand. However, in a model where
there is no correlation between a mutation's effect on fitness and the effect
on the trait, common variants account for much of the additive genetic
variance, regardless of demography. Moreover, here demography does not affect
the number of significant association detected. These finding suggest recent
population history may be an important factor influencing the power of
association tests in accounting for the missing heritability of certain complex
traits.


  Temperature affects both the timing and outcome of animal development, but
the detailed effects of temperature on the progress of early development have
been poorly characterized. To determine the impact of temperature on the order
and timing of events during Drosophila melanogaster embryogenesis, we used
time-lapse imaging to track the progress of embryos from shortly after egg
laying through hatching at seven precisely maintained temperatures between
17.5C and 32.5C. We employed a combination of automated and manual annotation
to determine when 36 milestones occurred in each embryo. D. melanogaster
embryogenesis takes ~33 hours at 17.5C, and accelerates with increasing
temperature to 16 hours at 27.5C, above which embryogenesis slows slightly.
Remarkably, while the total time of embryogenesis varies over two fold, the
relative timing of events from cellularization through hatching is constant
across temperatures. To further explore the relationship between temperature
and embryogenesis, we expanded our analysis to cover ten additional Drosophila
species of varying climatic origins. Six of these species, like D.
melanogaster, are of tropical origin, and embryogenesis time at different
temperatures was similar for them all. D. mojavensis, a sub-tropical fly,
develops slower than the tropical species at lower temperatures, while D.
virilis, a temperate fly, exhibits slower development at all temperatures. The
alpine sister species D. persimilis and D. pseudoobscura develop as rapidly as
tropical flies at cooler temperatures, but exhibit diminished acceleration
above 22.5C and have drastically slowed development by 30C. Despite ranging
from 13 hours for D. erecta at 30C to 46 hours for D. virilis at 17.5C, the
relative timing of events from cellularization through hatching is constant
across all species and temperatures, suggesting the existence of a timer
controlling embryogenesis.


  The conformational complexity of linear polymers far exceeds that of
point-like atoms and molecules. Polymers can bend, twist, even become knotted.
Thus they may also display a much richer phase structure than point particles.
But it is not very easy to characterize the phase of a polymer. Essentially,
the only attribute is the radius of gyration. The way how it changes when the
degree of polymerization becomes different, and how it evolves when the ambient
temperature and solvent properties change, discloses the phase of the polymer.
Moreover, in any finite length chain there are corrections to scaling, that
complicate the detailed analysis of the phase structure. Here we introduce a
quantity that we call the folding angle, a novel tool to identify and
scrutinize the phases of polymers. We argue for a mean-field relationship
between its values and those of the scaling exponent in the radius of gyration.
But unlike in the case of the radius of gyration, the value of the folding
angle can be evaluated from a single structure. As an example we estimate the
value of the folding angle in the case of crystallographic {\alpha}-helical
protein structures in the Protein Data Bank (PDB). We also show how the value
can be numerically computed using a theoretical model of {\alpha}-helical
chiral homopolymers.


  The digital nature of genes combined with the associated low copy numbers of
proteins regulating them is a significant source of stochasticity, which
affects the phase of biochemical oscillations. We provide a theoretical
framework for understanding the dephasing evolution of genetic oscillations by
combining the phenomenological stochastic limit cycle dynamics and the discrete
Markov state models that describe the genetic oscillations. Through simulations
of the realistic model of the NF\kappa B/I\kappa B network we illustrate the
dephasing phenomena which are important for reconciling single cell and
population based experiments on this system.


  Mississippi River nutrient loads and water stratification on the
Louisiana-Texas shelf contribute to an annually recurring, short-lived hypoxic
bottom layer in areas of the northern Gulf of Mexico comprising less than 2% of
the total Gulf of Mexico bottom area. Many publications demonstrate increases
in biomass and fisheries production attributed to nutrient loading from river
plumes. Decreases in fisheries production when nutrient loads are decreased are
also well documented. However, the National Oceanic and Atmospheric
Administration (NOAA) persists in describing the area adjacent to the
Mississippi River discharge as a "dead zone" and predicting dire consequences
if nutrient loads are not reduced. In reality, these areas teem with aquatic
life and provide 70-80% of the Gulf of Mexico fishery production. On June 18,
2013, NOAA published a misleading figure purporting to show the "dead zone" in
an article predicting a possible record dead zone area for 2013
(http://www.noaanews.noaa.gov/stories2013/20130618_deadzone.html). This area is
not a region of hypoxic bottom water at all nor is it related directly to 2013
predicted hypoxia. This figure appeared as early as 2004 in a National
Aeronautics and Space Administration (NASA) article
(http://www.nasa.gov/vision/earth/environment/dead_zone.html) as a satellite
image where the red area represents turbidity and is much larger than the
short-lived areas of hypoxic bottom water documented in actual NOAA
measurements. Thus, it is misleading for NOAA to characterize the red area in
that image as a "dead zone." The NOAA has also published other misleading and
exaggerated descriptions of the consequences of nutrient loading.


  The prevalence of sexual reproduction ("sex") in eukaryotes is an enigma of
evolutionary biology. Sex increases genetic variation only tells its long-term
superiority in essence. The accumulation of harmful mutations causes an
immediate and ubiquitous pressure for organisms. Contrary to the common sense,
our theoretical model suggests that reproductive rate can influence
initiatively the accumulation of harmful mutations. The interaction of
reproductive rate and the integrated harm of mutations causes a critical
reproductive rate R*. A population will become irreversibly extinct once the
reproductive rate reduces to lower than R*. A sexual population has a R* lower
than 1 and an asexual population has a R* higher than 1. The mean reproductive
rate of a population reached to the carrying capacity has to reduce to 1. That
explains the widespread sex as well as the persistence of facultative and
asexual organisms. Computer simulations support significantly our conclusion.


  This paper uses a recently presented abstract, tuneable Boolean regulatory
network model extended to consider aspects of mobile DNA, such as transposons.
The significant role of mobile DNA in the evolution of natural systems is
becoming increasingly clear. This paper shows how dynamically controlling
network node connectivity and function via transposon-inspired mechanisms can
be selected for in computational intelligence tasks to give improved
performance. The designs of dynamical networks intended for implementation
within the slime mould Physarum polycephalum and for the distributed control of
a smart surface are considered.


  Summary: We develop TaSer (TabAnno and SeqMiner), a toolkit for annotating
and querying next generation sequence (NGS) dataset in tab-delimited files.
TabAnno is a powerful and efficient command-line tool designed to pre-process
sequence data, annotate variations and generate an indexed feature-enriched
project file that can integrate multiple sources of information. Using the
project file generated by TabAnno, complex queries to the sequence dataset can
be performed using SeqMiner, an R-package designed to efficiently access large
datasets. Extracted information can be conveniently viewed and analyzed by
tools in R. TaSer is optimized and computationally more efficient than software
using database systems. It enables annotating and querying NGS dataset using
moderate computing resource. Availability and implementation: TabAnno can be
downloaded from github (zhanxw.github.io/anno/). SeqMiner is distributed on
CRAN (cran.r-project.org/web/packages/seqminer). Contact: X.Z.
(zhanxw@umich.edu) D.J.L (dajiang@umich.edu)


  Modern ecology has re-emphasized the need for a quantitative understanding of
the original 'survival of the fittest theme' based on analyzis of the intricate
trade-offs between competing evolutionary strategies that characterize the
evolution of life. This is key to the understanding of species coexistence and
ecosystem diversity under the omnipresent constraint of limited resources. In
this work we propose an agent based model replicating a community of
interacting individuals, e.g. plants in a forest, where all are competing for
the same finite amount of resources and each competitor is characterized by a
specific growth-reproduction strategy. We show that such an evolution dynamics
drives the system towards a stationary state characterized by an emergent
optimal strategy, which in turn depends on the amount of available resources
the ecosystem can rely on. We find that the share of resources used by
individuals is power-law distributed with an exponent directly related to the
optimal strategy. The model can be further generalized to devise optimal
strategies in social and economical interacting systems dynamics.


  The increasing interest in photoactivated proteins as natural replacement of
standard inorganic materials in photocells drives to the compared analysis of
bacteriorhodopsin and proteorhodopsin, two widely diffused proteins belonging
to the family of \textit{type-1} opsins. These proteins share similar
behaviours but exhibit relevant differences in the sequential chain of the
amino acids constituting their tertiary structure. The use of an impedance
network analogue to model the protein main features provides a microscopic
interpretation of a set of experiments on their photoconductance properties. In
particular, this model links the protein electrical responses to the tertiary
structure and to the interactions among neighbouring amino acids. The same
model is also used to predict the small-signal response in terms of the Nyquist
plot. Interesting enough, these rhodopsins are found to behave like a wide gap
semiconductor with intrinsic conductivities of the order of $10^{-7}$ S/cm.


  Bacterial plant pathogens rely on a battalion of transcription factors to
fine-tune their response to changing environmental conditions and marshal the
genetic resources required for successful pathogenesis. Prediction of
transcription factor binding sites represents an important tool for elucidating
regulatory networks, and has been conducted in multiple genera of plant
pathogenic bacteria for the purpose of better understanding mechanisms of
survival and pathogenesis. The major categories of transcription factor binding
sites that have been characterized are reviewed here with emphasis on in silico
methods used for site identification and challenges therein, their
applicability to different types of sequence datasets, and insights into
mechanisms of virulence and survival that have been gained through binding site
mapping. An improved strategy for establishing E value cutoffs when using
existing models to screen uncharacterized genomes is also discussed.


  Robustness of ecological flow networks under random failure of arcs is
considered with respect to two different functionalities: coherence and
circulation. In our previous work, we showed that each functionality is
associated with a natural path notion: lateral path for the former and directed
path for the latter. Robustness of a network is measured in terms of the size
of the giant laterally connected arc component and that of the giant strongly
connected arc component, respectively. We study how realistic structures of
ecological flow networks affect the robustness with respect to each
functionality. To quantify the impact of realistic network structures, two null
models are considered for a given real ecological flow network: one is random
networks with the same degree distribution and the other is those with the same
average degree. Robustness of the null models is calculated by theoretically
solving the size of giant components for the configuration model. We show that
realistic network structures have positive effect on robustness for coherence,
whereas they have negative effect on robustness for circulation.


  Meloidogyne root knot nematodes (RKN) can infect most of the world's
agricultural crop species and are among the most important of all plant
pathogens. As yet however we have little understanding of their origins or the
genomic basis of their extreme polyphagy. The most damaging pathogens reproduce
by mitotic parthenogenesis and are suggested to originate by interspecific
hybridizations between unknown parental taxa. We sequenced the genome of the
diploid meiotic parthenogen Meloidogyne floridensis, and use a comparative
genomic approach to test the hypothesis that it was involved in the hybrid
origin of the tropical mitotic parthenogen M. incognita. Phylogenomic analysis
of gene families from M. floridensis, M. incognita and an outgroup species M.
hapla was used to trace the evolutionary history of these species' genomes,
demonstrating that M. floridensis was one of the parental species in the hybrid
origins of M. incognita. Analysis of the M. floridensis genome revealed many
gene loci present in divergent copies, as they are in M. incognita, indicating
that it too had a hybrid origin. The triploid M. incognita is shown to be a
complex double-hybrid between M. floridensis and a third, unidentified parent.
The agriculturally important RKN have very complex origins involving the mixing
of several parental genomes by hybridization and their extreme polyphagy and
agricultural success may be related to this hybridization, producing
transgressive variation on which natural selection acts. Studying RKN variation
via individual marker loci may fail due to the species' convoluted origins, and
multi-species population genomics is essential to understand the hybrid
diversity and adaptive variation of this important species complex. This
comparative genomic analysis provides a compelling example of the importance
and complexity of hybridization in generating animal species diversity more
generally.


  We model a spatially detailed, two-sex population dynamics, to study the cost
of ecological restoration. We assume that cost is proportional to the number of
individuals introduced into a large habitat. We treat dispersal as homogeneous
diffusion. The local population dynamics depends on sex ratio at birth, and
allows mortality rates to differ between sexes. Furthermore, local density
dependence induces a strong Allee effect, implying that the initial population
must be sufficiently large to avert rapid extinction. We address three
different initial spatial distributions for the introduced individuals; for
each we minimize the associated cost, constrained by the requirement that the
species must be restored throughout the habitat. First, we consider spatially
inhomogeneous, unstable stationary solutions of the model's equations as
plausible candidates for small restoration cost. Second, we use numerical
simulations to find the smallest cluster size, enclosing a spatially
homogeneous population density, that minimizes the cost of assured restoration.
Finally, by employing simulated annealing, we minimize restoration cost among
all possible initial spatial distributions of females and males. For biased sex
ratios, or for a significant between-sex difference in mortality, we find that
sex-specific spatial distributions minimize the cost. But as long as the sex
ratio maximizes the local equilibrium density for given mortality rates, a
common homogeneous distribution for both sexes that spans a critical distance
yields a similarly low cost.


  Metagenomics provides a powerful new tool set for investigating evolutionary
interactions with the environment. However, an absence of model-based
statistical methods means that researchers are often not able to make full use
of this complex information. We present a Bayesian method for inferring the
phylogenetic relationship among related organisms found within metagenomic
samples. Our approach exploits variation in the frequency of taxa among samples
to simultaneously infer each lineage haplotype, the phylogenetic tree
connecting them, and their frequency within each sample. Applications of the
algorithm to simulated data show that our method can recover a substantial
fraction of the phylogenetic structure even in the presence of strong mixing
among samples. We provide examples of the method applied to data from green
sulfur bacteria recovered from an Antarctic lake, plastids from mixed
Plasmodium falciparum infections, and virulent Neisseria meningitidis samples.


  The availability of data from many different sources and fields of science
has made it possible to map out an increasing number of networks of contacts
and interactions. However, quantifying how reliable these data are remains an
open problem. From Biology to Sociology and Economy, the identification of
false and missing positives has become a problem that calls for a solution. In
this work we extend one of newest, best performing models -due to Guimera and
Sales-Pardo in 2009- to directed networks. The new methodology is able to
identify missing and spurious directed interactions, which renders it
particularly useful to analyze data reliability in systems like trophic webs,
gene regulatory networks, communication patterns and social systems. We also
show, using real-world networks, how the method can be employed to help
searching for new interactions in an efficient way.


  L\'evy flights represent the best strategy to randomly search for a target in
an unknown environment, and have been widely observed in many animal species.
Here, we inspect and discuss recent results concerning human behavior and
cognition. Different studies have shown that human mobility can be described in
terms of L\'evy flights, while fresh evidence indicates that the same pattern
accounts for human mental searches in online gambling sites. Thus, L\'evy
flights emerge as a unifying concept with broad cross-disciplinary
implications. We argue that the ubiquity of such a pattern, both in behavior
and cognition, suggests that the brain regions responsible for this behavior
are likely to be evolutionarily old (i.e. no frontal cortex is involved), and
that fMRI techniques might help to confirm this hypothesis.


  Background.
  Models of ancestral gene order reconstruction have progressively integrated
different evolutionary patterns and processes such as unequal gene content,
gene duplications, and implicitly sequence evolution via reconciled gene trees.
In unicellular organisms, these models have so far ignored lateral gene
transfer, even though it can have an important confounding effect on such
models, as well as a rich source of information on the function of genes
through the detection of transfers of entire clusters of genes.
  Result.
  We report an algorithm together with its implementation, DeCoLT, that
reconstructs ancestral genome organization based on reconciled gene trees which
summarize information on sequence evolution, gene origination, duplication,
loss, and lateral transfer. DeCoLT finds in polynomial time the minimum number
of rearrangements, computed as the number of gains and breakages of adjacencies
between pairs of genes. We apply DeCoLT to 1099 gene families from 36
cyanobacteria genomes.
  Conclusion.
  DeCoLT is able to reconstruct adjacencies in 35 ancestral bacterial genomes
with a thousand genes families in a few hours, and detects clusters of
co-transferred genes. As there is no constraint on genome organization,
adjacencies can be generalized to any relationship between genes to reconstruct
ancestral interactions, functions or complexes with the same framework.


  We analyze the force-extension curve for a general class of systems, which
are described at the mesoscopic level by a free energy depending on the
extension of its components. Similarly to what is done in real experiments, the
total length of the system is the controlled parameter. This imposes a global
constraint in the minimization procedure leading to the equilibrium values of
the extensions. As a consequence, the force-extension curve has multiple
branches in a certain range of forces. The stability of these branches is
governed by the free energy: there are a series of first-order phase
transitions at certain values of the total length, in which the free energy
itself is continuous but its first derivative, the force, has a finite jump.
This behavior is completely similar to the one observed in real experiments
with biomolecules like proteins, and other complex systems.


  Single virus epidemics over complete networks are widely explored in the
literature as the fraction of infected nodes is, under appropriate microscopic
modeling of the virus infection, a Markov process. With non-complete networks,
this macroscopic variable is no longer Markov. In this paper, we study virus
diffusion, in particular, multi-virus epidemics, over non-complete stochastic
networks. We focus on multipartite networks. In companying work
http://arxiv.org/abs/1306.6198, we show that the peer-to-peer local random
rules of virus infection lead, in the limit of large multipartite networks, to
the emergence of structured dynamics at the macroscale. The exact fluid limit
evolution of the fraction of nodes infected by each virus strain across islands
obeys a set of nonlinear coupled differential equations, see
http://arxiv.org/abs/1306.6198. In this paper, we develop methods to analyze
the qualitative behavior of these limiting dynamics, establishing conditions on
the virus micro characteristics and network structure under which a virus
persists or a natural selection phenomenon is observed.


  The appearance of cancer in a tissue is thought to be the result of two or
more successive mutations. We propose a stochastic model that allows for an
exact computation of the distribution of the waiting time for a second
mutation. This models the time of appearance of the first cancerous cell in a
tissue. Our model is an alternative to the Moran model with mutations.


  The control of complex networks is of paramount importance in areas as
diverse as ecosystem management, emergency response, and cell reprogramming. A
fundamental property of networks is that perturbations to one node can affect
other nodes, potentially causing the entire system to change behavior or fail.
Here, we show that it is possible to exploit the same principle to control
network behavior. Our approach accounts for the nonlinear dynamics inherent to
real systems, and allows bringing the system to a desired target state even
when this state is not directly accessible due to constraints that limit the
allowed interventions. Applications show that this framework permits
reprogramming a network to a desired task as well as rescuing networks from the
brink of failure---which we illustrate through the mitigation of cascading
failures in a power-grid network and the identification of potential drug
targets in a signaling network of human cancer.


  The human genome remains incomplete, with multi-megabase sized gaps
representing the endogenous centromeres and other heterochromatic regions.
These regions are commonly enriched with long arrays of near-identical tandem
repeats, known as satellite DNAs, that offer a limited number of variant sites
to differentiate individual repeat copies across millions of bases. This
substantial sequence homogeneity challenges available assembly strategies, and
as a result, centromeric regions are omitted from ongoing genomic studies. To
address this problem, we present a locally ordered assembly across two haploid
human satellite arrays on chromosomes X and Y, resulting in an initial linear
representation of 3.83 Mb of centromeric DNA within an individual genome. To
further expand the utility of each centromeric reference sequence, we evaluate
sites within the arrays for short-read mappability and chromosome specificity.
As satellite DNAs evolve in a concerted manner, we use these centromeric
assemblies to assess the extent of sequence variation among 372 individuals
from distinct human populations. In doing so, we identify two ancient satellite
array variants in both X and Y centromeres as determined by array length and
sequence composition. This study provides an initial linear representation and
comprehensive sequence characterization of a regional centromere and
establishes a foundation to extend genomic characterization to these sites as
well as to other repeat-rich regions within complex genomes.


  We suggest to discriminate single DNA bases via transverse ionic transport,
namely by detecting the ionic current that flows in a channel while a
single-stranded DNA is driven through an intersecting nanochannel. Our all-atom
molecular dynamics simulations indeed show that the ionic currents of the four
bases are statistically distinct, thus offering another possible approach to
sequence DNA.


  We present a detailed dynamical model of the behavior of
transcription-translation circuits in vitro that makes explicit the roles
played by essential molecular resources. A set of simple two-gene test circuits
operating in a cell-free biochemical 'breadboard' validate this model and
highlight the consequences of limited resource availability. In particular, we
are able to confirm the existence of biomolecular 'crosstalk' and isolate its
individual sources. The implications of crosstalk for biomolecular circuit
design and function are discussed.


  This paper presents a novel DNA sequences alignment method based on inverted
index. Now most large scale information retrieval system are all use inverted
index as the basic data structure. But its application in DNA sequence
alignment is still not found. This paper just discuss such applications. Three
main problems, DNA segmenting, long DNA query search, DNA search ranking
algorithm and evaluation method are detailed respectively. This research
presents a new avenue to build more effective DNA alignment methods.


  Epidemiological models may give some basic guidelines for public health
practitioners, allowing to analyze issues that can influence the strategies to
prevent and fight a disease. To be used in decision-making, however, a
mathematical model must be carefully parameterized and validated with
epidemiological and entomological data. Here a SIR (S for susceptible, I for
infectious, R for recovered individuals) and ASI (A for the aquatic phase of
the mosquito, S for susceptible and I for infectious mosquitoes)
epidemiological model describing a dengue disease is presented, as well as the
associated basic reproduction number. A sensitivity analysis of the
epidemiological model is performed in order to determine the relative
importance of the model parameters to the disease transmission.


  We studied the structure and mechanical properties of DNA i-motif nanowires
by means of molecular dynamics computer simulations. We built up to 230 nm long
nanowires, based on a repeated TC5 sequence from crystallographic data, fully
relaxed and equilibrated in water. The unusual stacked C*C+ stacked structure,
formed by four ssDNA strands arranged in an intercalated tetramer, is here
fully characterized both statically and dynamically. By applying stretching,
compression and bending deformation with the steered molecular dynamics and
umbrella sampling methods, we extract the apparent Young's and bending moduli
of the nanowire, as wel as estimates for the tensile strength and persistence
length. According to our results, the i-motif nanowire shares similarities with
structural proteins, as far as its tensile stiffness, but is closer to nucleic
acids and flexible proteins, as far as its bending rigidity is concerned.
Furthermore, thanks to its very thin cross section, the apparent tensile
toughness is close to that of a metal. Besides their yet to be clarified
biological significance, i-motif nanowires may qualify as interesting
candidates for nanotechnology templates, due to such outstanding mechanical
properties.


  Kidney branching morphogenesis has been studied extensively, but the
mechanism that defines the branch points is still elusive. Here we obtained a
2D movie of kidney branching morphogenesis in culture to test different models
of branching morphogenesis with physiological growth dynamics. We carried out
image segmentation and calculated the displacement fields between the frames.
The models were subsequently solved on the 2D domain, that was extracted from
the movie. We find that Turing patterns are sensitive to the initial conditions
when solved on the epithelial shapes. A previously proposed diffusion-dependent
geometry effect allowed us to reproduce the growth fields reasonably well, both
for an inhibitor of branching that was produced in the epithelium, and for an
inducer of branching that was produced in the mesenchyme. The latter could be
represented by Glial-derived neurotrophic factor (GDNF), which is expressed in
the mesenchyme and induces outgrowth of ureteric branches. Considering that the
Turing model represents the interaction between the GDNF and its receptor RET
very well and that the model reproduces the relevant expression patterns in
developing wildtype and mutant kidneys, it is well possible that a combination
of the Turing mechanism and the geometry effect control branching
morphogenesis.


  Tests of the neutral evolution hypothesis are usually built on the standard
null model which assumes that mutations are neutral and population size remains
constant over time. However, it is unclear how such tests are affected if the
last assumption is dropped. Here, we extend the unifying framework for tests
based on the site frequency spectrum, introduced by Achaz and Ferretti, to
populations of varying size. A key ingredient is to specify the first two
moments of the frequency spectrum. We show that these moments can be determined
analytically if a population has experienced two instantaneous size changes in
the past. We apply our method to data from ten human populations gathered in
the 1000 genomes project, estimate their demographies and define
demography-adjusted versions of Tajima's $D$, Fay & Wu's $H$, and Zeng's $E$.
The adjusted test statistics facilitate the direct comparison between
populations and they show that most of the differences among populations seen
in the original tests can be explained by demography. We carried out whole
genome screens for deviation from neutrality and identified candidate regions
of recent positive selection. We provide track files with values of the
adjusted and original tests for upload to the UCSC genome browser.


  Branching processes are classical growth models in cell kinetics. In their
construction, it is usually assumed that cell lifetimes are independent random
variables, which has been proved false in experiments. Models of dependent
lifetimes are considered here, in particular bifurcating Markov chains. Under
hypotheses of stationarity and multiplicative ergodicity, the corresponding
branching process is proved to have the same type of asymptotics as its classic
counterpart in the i.i.d. supercritical case: the cell population grows
exponentially, the growth rate being related to the exponent of multiplicative
ergodicity, in a similar way as to the Laplace transform of lifetimes in the
i.i.d. case. An identifiable model for which the multiplicative ergodicity
coefficients and the growth rate can be explicitly computed is proposed.


  Monolayer spontaneous curvatures for cholesterol, DOPE, POPE, DOPC, DPPC,
DSPC, POPC, SOPC, and egg sphingomyelin were obtained using small-angle X-ray
scattering (SAXS) on inverted hexagonal phases (HII). Spontaneous curvatures of
bilayer forming lipids were estimated by adding controlled amounts to a HII
forming template following previously established protocols. Spontanous
curvatures of both phosphatidylethanolamines and cholesterol were found to be
at least a factor of two more negative than those of phosphatidylcholines,
whose J0 are closer to zero. Interestingly, a significant positive J0 value
(+0.1 1/nm) was retrieved for DPPC at 25 {\deg}C. We further determined the
temperature dependence of the spontaneous curvatures J0(T) in the range from 15
to 55 \degC, resulting in a quite narrow distribution of -1 to -3 * 10^-3
1/nm{\deg}C for most investigated lipids. The data allowed us to estimate the
monolayer spontaneous curvatures of ternary lipid mixtures showing liquid
ordered / liquid disordered phase coexistence. We report spontaneous curvature
phase diagrams for DSPC/DOPC/Chol, DPPC/DOPC/Chol and SM/POPC/Chol and discuss
effects on protein insertion and line tension.


  Reciprocity is a vital feature of social networks, but relatively little is
known about its temporal structure or the mechanisms underlying its persistence
in real world behavior. In pursuit of these two questions, we study the
stationary and dynamical signals of reciprocity in a network of manioc beer
(Spanish: chicha; Tsimane': shocdye') drinking events in a Tsimane' village in
lowland Bolivia. At the stationary level, our analysis reveals that social
exchange within the community is heterogeneously patterned according to kinship
and spatial proximity. A positive relationship between the frequencies at which
two families host each other, controlling for kinship and proximity, provides
evidence for stationary reciprocity. Our analysis of the dynamical structure of
this network presents a novel method for the study of conditional, or
non-stationary, reciprocity effects. We find evidence that short-timescale
reciprocity (within three days) is present among non- and distant-kin pairs;
conversely, we find that levels of cooperation among close kin can be accounted
for on the stationary hypothesis alone.


  Processes that proceed reliably from a variety of initial conditions to a
unique final form, regardless of moderately changing conditions, are of obvious
importance in biophysics. Protein folding is a case in point. We show that the
action principle can be applied directly to study the stability of biological
processes. The action principle in classical physics starts with the first
variation of the action and leads immediately to the equations of motion. The
second variation of the action leads in a natural way to powerful theorems that
provide quantitative treatment of stability and focusing and also explain how
some very complex processes can behave as though some seemingly important
forces drop out. We first apply these ideas to the non-equilibrium states
involved in two-state folding. We treat torsional waves and use the action
principle to talk about critical points in the dynamics. For some proteins the
theory resembles TST. We reach several quantitative and qualitative
conclusions. Besides giving an explanation of why TST often works in folding,
we find that the apparent smoothness of the energy funnel is a natural
consequence of the putative critical points in the dynamics. These ideas also
explain why biological proteins fold to unique states and random polymers do
not. The insensitivity to perturbations which follows from the presence of
critical points explains how folding to a unique shape occurs in the presence
of dilute denaturing agents in spite of the fact that those agents disrupt the
folded structure of the native state. This paper contributes to the theoretical
armamentarium by directing attention to the logical progression from first
physical principles to the stability theorems related to catastrophe theory as
applied to folding. This can potentially have the same success in biophysics as
it has enjoyed in optics.


  In the past several years, the problem of genome assembly has received
considerable attention from both biologists and computer scientists. An
important component of current assembly methods is the scaffolding process.
This process involves building ordered and oriented linear collections of
contigs (continuous overlapping sequence reads) called scaffolds and relies on
the use of mate pair data. A mate pair is a set of two reads that are sequenced
from the ends of a single fragment of DNA, and therefore have opposite mutual
orientations. When two reads of a mate-pair are placed into two different
contigs, one can infer the mutual orientation of these contigs. While several
orientation algorithms exist as part of assembly programs, all encounter
challenges while solving the orientation problem due to errors from
mis-assemblies in contigs or errors in read placements. In this paper we
present an algorithm based on hierarchical clustering that independently solves
the orientation problem and is robust to errors. We show that our algorithm can
correctly solve the orientation problem for both faux (generated) assembly data
and real assembly data for {\em R. sphaeroides bacteria}. We demonstrate that
our algorithm is stable to both changes in the initial orientations as well as
noise in the data, making it advantageous compared to traditional approaches.


  Development of resistance limits transferability of most anticancer therapies
into curative treatment and understanding mechanisms beyond it remains a big
challenge. Many high resolution experimental observations show enormous
intratumor heterogeneity at molecular, genetic and cellular levels which is
made responsible for emerging resistance to therapy. Therefore, researchers
search techniques to influence development of intratumor heterogeneity, which
requires understanding its role within the context of integrative, logically
consistent, framework, such as evolutionary theory. Although it is agreed that
intratumor heterogeneity increases probability of the emergence of therapy
resistant clones, more instructive role of its structure in the process of
cancer dynamics and metastasis is needed. In the paper, intratumor
heterogeneity is viewed as a product of two, in general stochastic, processes,
evolutionary optimization and changing environment, respectively. In
evolutionary theory, common risk-diversifying strategy displayed by isogenic
populations in unpredictably changing environments is bet-hedging. We suggest,
that the structure of intratumor heterogeneity is evolutionary trait evolving
to maximize the clonal fitness in changing (or uncertain) environment and that
its structure corresponds to bet-hedging strategy. We advocate our view by
reviewing and combining important cancer relevant concepts.


  Dual trap optical tweezers are often used in high-resolution measurements in
single-molecule biophysics. Such measurements can be hindered by the presence
of extraneous noise sources, the most prominent of which is the coupling of
fluctuations along different spatial directions, which may affect any optical
tweezers setup. In this paper we analyze, both from the theoretical and the
experimental points of view, the most common source for these couplings in Dual
Trap Optical Tweezers setups: the misalignment of traps and tether. We give
criteria to distinguish different kinds of misalignment, to estimate their
quantitative relevance and to include them in the data analysis. The
experimental data is obtained in a novel dual trap optical tweezers setup which
directly measures forces. In the case in which misalignment is negligible we
provide a method to measure the stiffness of traps and tether based on variance
analysis. This method can be seen as a calibration technique valid beyond the
linear trap region. Our analysis is then employed to measure the persistence
length of ds-DNA tethers of three different lengths spanning two orders of
magnitude. The effective persistence length of such tethers is shown to
decrease with the contour length, in accordance with previous studies.


  The effect of magnesium ion Mg2+ on the dielectric relaxation of semidilute
DNA aqueous solutions has been studied by means of dielectric spectroscopy in
the 100 Hz-100 MHz frequency range. De Gennes-Pfeuty-Dobrynin semidilute
solution correlation length is the pertinent fundamental length scale for
sufficiently low concentration of added salt, describing the collective
properties of Mg-DNA solutions. No relaxation fingerprint of the DNA
denaturation bubbles, leading to exposed hydrophobic core scaling, was detected
at low DNA concentrations, thus indicating an increased stability of the
double-stranded conformation in Mg-DNA solutions as compared to the case of
Na-DNA solutions. Some changes are detected in the behavior of the fundamental
length scale pertaining to the single molecule DNA properties, reflecting
modified electrostatic screening effects of the Odijk-Skolnick-Fixman type. All
results consistently demonstrate that Mg2+ ions interact with DNA in a similar
way as Na1+ ions do, their effect being mostly describable through an enhanced
screening.


  In large populations, multiple beneficial mutations may be simultaneously
spreading. In asexual populations, these mutations must either arise on the
same background or compete against each other. In sexual populations,
recombination can bring together beneficial alleles from different backgrounds,
but tightly linked alleles may still greatly interfere with each other. We show
for well-mixed populations that when this interference is strong, the genome
can be seen as consisting of many effectively asexual stretches linked
together. The rate at which beneficial alleles fix is thus roughly proportional
to the rate of recombination, and depends only logarithmically on the mutation
supply and the strength of selection. Our scaling arguments also allow to
predict, with reasonable accuracy, the distribution of effects of fixed
mutations when new mutations have broadly-distributed effects. We focus on the
regime in which crossovers occur more frequently than beneficial mutations, as
is likely to be the case for many natural populations.


  We present a dual-trap optical tweezers setup which directly measures forces
using linear momentum conservation. The setup uses a counter-propagating
geometry, which allows momentum measurement on each beam separately. The
experimental advantages of this setup include low drift due to all-optical
manipulation, and a robust calibration (independent of the features of the
trapped object or buffer medium) due to the force measurement method. Although
this design does not attain the high-resolution of some co-propagating setups,
we show that it can be used to perform different single molecule measurements:
fluctuation-based molecular stiffness characterization at different forces and
hopping experiments on molecular hairpins. Remarkably, in our setup it is
possible to manipulate very short tethers (such as molecular hairpins with
short handles) down to the limit where beads are almost in contact. The setup
is used to illustrate a novel method for measuring the stiffness of optical
traps and tethers on the basis of equilibrium force fluctuations, i.e. without
the need of measuring the force vs molecular extension curve. This method is of
general interest for dual trap optical tweezers setups and can be extended to
setups which do not directly measure forces.


  Many biological systems, such as metabolic pathways, exhibit bistability
behavior: these biological systems exhibit two distinct stable states with
switching between the two stable states controlled by certain conditions. Since
understanding bistability is key for understanding these biological systems,
mathematical modeling of the bistability phenomenon has been at the focus of
researches in quantitative and system biology. Recent study shows that Boolean
networks offer relative simple mathematical models that are capable of
capturing these essential information. Thus a better understanding of the
Boolean networks with bistability property is desirable for both theoretical
and application purposes. In this paper, we describe an algebraic condition for
the number of stable states (fixed points) of a Boolean network based on its
polynomial representation, and derive algorithms for a Boolean network to have
a single stable state or two stable states. As an example, we also construct a
Boolean network with exactly two stable states for the lac operon's
$\beta$-galactosidase regulatory pathway when glucose is absent based on a
delay differential equation model


  With the astonishing rate that the genomic and metagenomic sequence data sets
are accumulating, there are many reasons to constrain the data analyses. One
approach to such constrained analyses is to focus on select subsets of gene
families that are particularly well suited for the tasks at hand. Such gene
families have generally been referred to as marker genes. We are particularly
interested in identifying and using such marker genes for phylogenetic and
phylogeny-driven ecological studies of microbes and their communities. We
therefore refer to these as PhyEco (for phylogenetic and phylogenetic ecology)
markers. The dual use of these PhyEco markers means that we needed to develop
and apply a set of somewhat novel criteria for identification of the best
candidates for such markers. The criteria we focused on included universality
across the taxa of interest, ability to be used to produce robust phylogenetic
trees that reflect as much as possible the evolution of the species from which
the genes come, and low variation in copy number across taxa. We describe here
an automated protocol for identifying potential PhyEco markers from a set of
complete genome sequences. The protocol combines rapid searching, clustering
and phylogenetic tree building algorithms to generate protein families that
meet the criteria listed above. We report here the identification of PhyEco
markers for different taxonomic levels including 40 for all bacteria and
archaea, 114 for all bacteria, and much more for some of the individual phyla
of bacteria. This new list of PhyEco markers should allow much more detailed
automated phylogenetic and phylogenetic ecology analyses of these groups than
possible previously.


  Large sets of genotypes give rise to the same phenotype because phenotypic
expression is highly redundant. Accordingly, a population can accept mutations
without altering its phenotype, as long as thegenotype mutates into another one
on the same set. By linking every pair of genotypes that are mutually
accessible through mutation, genotypes organize themselves into neutral
networks (NN). These networks are known to be heterogeneous and assortative,
and these properties affect the evolutionary dynamics of the population. By
studying the dynamics of populations on NN with arbitrary topology we analyze
the effect of assortativity, of NN (phenotype) fitness, and of network size. We
find that the probability that the population leaves the network is smaller the
longer the time spent on it. This progressive "phenotypic entrapment" entails a
systematic increase in the overdispersion of the process with time and an
acceleration in the fixation rate of neutral mutations. We also quantify the
variation of these effects with the size of the phenotype and with its fitness
relative to that of neighbouring alternatives.


  Gene expression is a central process to any form of life. It involves
multiple temporal and functional scales that extend from specific protein-DNA
interactions to the coordinated regulation of multiple genes in response to
intracellular and extracellular changes. This diversity in scales poses
fundamental challenges among traditional approaches to fully understand even
the simplest gene expression systems. Recent advances in computational systems
biophysics have provided promising avenues to reliably integrate the molecular
detail of biophysical process into the system behavior. Here, we review recent
advances in the description of gene regulation as a system of biophysical
processes that extend from specific protein-DNA interactions to the
combinatorial assembly of nucleoprotein complexes. There is now basic
mechanistic understanding on how promoters controlled by multiple, local and
distal, DNA binding sites for transcription factors can actively control
transcriptional noise, cell-to-cell variability, and other properties of gene
regulation, including precision and flexibility of the transcriptional
responses.


  We model and discuss advantages of pan-editing, the complex way of expressing
mitochondrial genes in kinetoplastids. The rapid spread and preservation of
pan-editing seems to be due to its concomitant fragment dispersal. Such
dispersal prevents losing temporarily non expressed mitochondrial genes upon
intense intraspecific competition, by linking non expressed fragments to parts
which are still needed. We mathematically modelled protection against gene
loss, due to the absence of selection, by this kind of fragment association.
This gives ranges of values for parameters like scrambling extent, population
size, and number of generations still retaining full genomes despite limited
selection. Values obtained seem consistent with those observed. We find a
quasi-linear correlation between dispersal and number of generations after
which populations lose genes, showing that pan-editing can be selected to
effectively limit gene loss under relaxed selective pressure.


  A methodology is proposed to automatically detect significant symbol
associations in genomic databases. A new statistical test is proposed to assess
the significance of a group of symbols when found in several genesets of a
given database. Applied to symbol pairs, the thresholded p-values of the test
define a graph structure on the set of symbols. The cliques of that graph are
significant symbol associations, linked to a set of genesets where they can be
found. The method can be applied to any database, and is illustrated MSigDB C2
database. Many of the symbol associations detected in C2 or in non-specific
selections did correspond to already known interactions. On more specific
selections of C2, many previously unkown symbol associations have been
detected. These associations unveal new candidates for gene or protein
interactions, needing further investigation for biological evidence.


  MicroRNAs (miRNAs) are small non-coding RNAs which regulate gene expression
by binding to the 3' UTR of the corresponding messenger RNAs. We construct
miRNA co-target networks for 22 different species using a target prediction
database, MicroCosm Tagets. The miRNA pairs of individual species having one or
more common target genes are connected and the number of co-targets are
assigned as the weight of these links. We show that the link-weight
distributions of all the species collapse remarkably onto each other when
scaled suitably. It turns out that the scale-factor is a measure of complexity
of the species. A simple model, where targets are chosen randomly by miRNAs,
could provide the correct scaling function and explain the universality.


  Molecular signaling networks are believed to determine cancer robustness.
Although cancer patient survivability was reported to correlate with the
heterogeneous connectivity of the signaling networks inspired by theoretical
studies on the increase of network robustness due to the heterogeneous
connectivity, other theoretical and data analytic studies suggest an
alternative explanation: the impact of modular organization of networks on
biological robustness or adaptation to changing environments. In this study,
thus, we evaluate whether the modularity--robustness hypothesis is applicable
to cancer using network analysis. We focus on 14 specific cancer types whose
molecular signaling networks are available in databases, and show that modular
organization of cancer signaling networks is associated with the patient
survival rate. In particular, the cancers with less modular signaling networks
are more curable. This result is consistent with a prediction from the
modularity--robustness hypothesis. Furthermore, we show that the network
modularity is a better descriptor of the patient survival rate than the
heterogeneous connectivity. However, these results do not contradict the
importance of the heterogeneous connectivity. Rather, they provide new and
different insights into the relationship between cellular networks and cancer
behaviors. Despite several limitations of data analysis, these findings enhance
our understanding of adaptive and evolutionary mechanisms of cancer cells.


  The primary structure of proteins, that is their sequence, represents one of
the most abundant set of experimental data concerning biomolecules. The study
of correlations in families of co--evolving proteins by means of an inverse
Ising--model approach allows to obtain information on their native
conformation. Following up on a recent development along this line, we optimize
the algorithm to calculate effective energies between the residues, validating
the approach both back-calculating interaction energies in a model system, and
predicting the free energies associated to mutations in real systems. Making
use of these effective energies, we study the networks of interactions which
stabilizes the native conformation of some well--studied proteins, showing that
it display different properties than the associated contact network.


  I propose an approach to identify, among several strategies of phylogenetic
analysis, those producing the most accurate results. This approach is based on
the hypothesis that the more a result is reproduced from independent data, the
more it reflects the historical signal common to the analysed data. Under this
hypothesis, the capacity of an analytical strategy to extract historical signal
should correlate positively with the coherence of the obtained results. I apply
this approach to a series of analyses on empirical data, basing the coherence
measure on the Robinson-Foulds distances between the obtained trees. At first
approximation, the analytical strategies most suitable for the data produce the
most coherent results. However, risks of false positives and false negatives
are identified, which are difficult to rule out.


  A new mathematical model for evolutionary games on graphs is proposed to
extend the classical replicator equation to finite populations of players
organized on a network with generic topology. Classical results from game
theory, evolutionary game theory and graph theory are used. More specifically,
each player is placed in a vertex of the graph and he is seen as an infinite
population of replicators which replicate within the vertex. At each time
instant, a game is played by two replicators belonging to different connected
vertices, and the outcome of the game influences their ability of producing
offspring. Then, the behavior of a vertex player is determined by the
distribution of strategies used by the internal replicators. Under suitable
hypotheses, the proposed model is equivalent to the classical replicator
equation. Extended simulations are performed to show the dynamical behavior of
the solutions and the potentialities of the developed model.


  Interactions among multiple infectious agents are increasingly recognized as
a fundamental issue in the understanding of key questions in public health,
regarding pathogen emergence, maintenance, and evolution. The full description
of host-multipathogen systems is however challenged by the multiplicity of
factors affecting the interaction dynamics and the resulting competition that
may occur at different scales, from the within-host scale to the spatial
structure and mobility of the host population. Here we study the dynamics of
two competing pathogens in a structured host population and assess the impact
of the mobility pattern of hosts on the pathogen competition. We model the
spatial structure of the host population in terms of a metapopulation network
and focus on two strains imported locally in the system and having the same
transmission potential but different infectious periods. We find different
scenarios leading to competitive success of either one of the strain or to the
codominance of both strains in the system. The dominance of the strain
characterized by the shorter or longer infectious period depends exclusively on
the structure of the population and on the the mobility of hosts across
patches. The proposed modeling framework allows the integration of other
relevant epidemiological, environmental and demographic factors opening the
path to further mathematical and computational studies of the dynamics of
multipathogen systems.


  Umbrella sampling is an efficient method for the calculation of free energy
changes of a system along well-defined reaction coordinates. However, when
multiple parallel channels along the reaction coordinate or hidden barriers in
directions perpendicular to the reaction coordinate exist, it is difficult for
conventional umbrella sampling methods to generate sufficient sampling within
limited simulation time. Here we propose an efficient approach to combine
umbrella sampling with the integrated tempering sampling method. The umbrella
sampling method is applied to conformational degrees of freedom which possess
significant barriers and are chemically more relevant. The integrated tempering
sampling method is employed to facilitate the sampling of other degrees of
freedom in which statistically non-negligible barriers may exist. The combined
method is applied to two model systems and show significantly improved sampling
efficiencies as compared to standalone conventional umbrella sampling or
integrated tempering sampling approaches. Therefore, the combined approach will
become a very efficient method in the simulation of biomolecular processes
which often involve sampling of complex rugged energy landscapes.


  A major revision of the geological timescale was published in 2012. We
re-examine our past finding of a 27 Myr periodicity in marine extinction rates
by re-assigning dates to the extinction data used previously. We find that the
spectral power at this period is somewhat increased, and persists at a narrow
bandwidth, which supports our previous contention that the Nemesis hypothesis
is untenable as an explanation for the periodicity that was first noted by Raup
and Sepkoski in the 1980s. We enumerate a number of problems in a recent study
comparing extinction rates with time series models.


  We scrutinize the temporally-resolved speed of active cargo transport in
living cells, and show intermittent bursting motions. These nonlinear
fluctuations follow a scaling law over several decades of time and space, the
statistical regularities displaying a time-averaged shape that we interpret to
reflect stress buildup followed by rapid release. The power law of scaling is
the same as seen in driven jammed colloids, granular, and magnetic systems. The
implied regulation of active transport with environmental obstruction extends
the classical notion of molecular crowding.


  Stabilizing the dynamics of complex, non-linear systems is a major concern
across several scientific disciplines including ecology and conservation
biology. Unfortunately, most methods proposed to reduce the fluctuations in
chaotic systems are not applicable for real, biological populations. This is
because such methods typically require detailed knowledge of system specific
parameters and the ability to manipulate them in real time; conditions often
not met by most real populations. Moreover, real populations are often noisy
and extinction-prone, which can sometimes render such methods ineffective. Here
we investigate a control strategy, which works by perturbing the population
size, and is robust to reasonable amounts of noise and extinction probability.
This strategy, called the Adaptive Limiter Control (ALC), has been previously
shown to increase constancy and persistence of laboratory populations and
metapopulations of Drosophila melanogaster. Here we present a detailed
numerical investigation of the effects of ALC on the fluctuations and
persistence of metapopulations. We show that at high migration rates,
application of ALC does not require a priori information about the population
growth rates. We also show that ALC can stabilize metapopulations even when
applied to as low as one-tenth of the total number of subpopulations. Moreover,
ALC is effective even when the subpopulations have high extinction rates:
conditions under which one other control algorithm has previously failed to
attain stability. Importantly, ALC not only reduces the fluctuation in
metapopulation sizes, but also the global extinction probability. Finally, the
method is robust to moderate levels of noise in the dynamics and the carrying
capacity of the environment. These results, coupled with our earlier empirical
findings, establish ALC to be a strong candidate for stabilizing real
biological metapopulations.


  It has recently been noted that the relative prevalence of the various kinds
of epistasis varies along an adaptive walk. This has been explained as a result
of mean regression in NK model fitness landscapes. Here we show that this
phenomenon occurs quite generally in fitness landscapes. We propose a simple
and general explanation for this phenomemon, confirming the role of mean
regression. We provide support for this explanation with simulations, and
discuss the empirical relevance of our findings.


  It has been well accepted that the RNA secondary structures of most
functional non-coding RNAs (ncRNAs) are closely related to their functions and
are conserved during evolution. Hence, prediction of conserved secondary
structures from evolutionarily related sequences is one important task in RNA
bioinformatics; the methods are useful not only to further functional analyses
of ncRNAs but also to improve the accuracy of secondary structure predictions
and to find novel functional RNAs from the genome. In this review, I focus on
common secondary structure prediction from a given aligned RNA sequence, in
which one secondary structure whose length is equal to that of the input
alignment is predicted. I systematically review and classify existing tools and
algorithms for the problem, by utilizing the information employed in the tools
and by adopting a unified viewpoint based on maximum expected gain (MEG)
estimators. I believe that this classification will allow a deeper
understanding of each tool and provide users with useful information for
selecting tools for common secondary structure predictions.


  We study the planar matching problem, defined by a symmetric random matrix
with independent identically distributed entries, taking values 0 and 1. We
show that the existence of a perfect planar matching structure is possible only
above a certain critical density, $p_{c}$, of allowed contacts (i.e. of '1').
Using a formulation of the problem in terms of Dyck paths and a matrix model of
planar contact structures, we provide an analytical estimation for the value of
the transition point, $p_{c}$, in the thermodynamic limit. This estimation is
close to the critical value, $p_{c} \approx 0.379$, obtained in numerical
simulations based on an exact dynamical programming algorithm. We characterize
the corresponding critical behavior of the model and discuss the relation of
the perfect-imperfect matching transition to the known molten-glass transition
in the context of random RNA secondary structure's formation. In particular, we
provide strong evidence supporting the conjecture that the molten-glass
transition at T=0 occurs at $p_{c}$.


  The detection of genomic structural variations (SV) remains a difficult
challenge in analyzing sequencing data, and the growing size and number of
sequenced genomes have rendered SV detection a bona fide big data problem.
MapReduce is a proven, scalable solution for distributed computing on huge data
sets. We describe a conceptual framework for SV detection algorithms in
MapReduce based on computing local genomic features, and use it to develop a
deletion and insertion detection algorithm, Cloudbreak. On simulated and real
data sets, Cloudbreak achieves accuracy improvements over popular SV detection
algorithms, and genotypes variants from diploid samples. It provides
dramatically shorter runtimes and the ability to scale to big data volumes on
large compute clusters. Cloudbreak includes tools to set up and configure
MapReduce (Hadoop) clusters on cloud services, enabling on-demand cluster
computing. Our implementation and source code are available at
http://github.com/cwhelan/cloudbreak .


  Modeling epidemic dynamics plays an important role in studying how diseases
spread, predicting their future course, and designing strategies to control
them. In this letter, we introduce a model of SIR
(susceptible-infected-removed) type which explicitly incorporates the effect of
{\it cooperative coinfection}. More precisely, each individual can get infected
by two different diseases, and an individual already infected with one disease
has an increased probability to get infected by the other. Depending on the
amount of this increase, we observe different threshold scenarios. Apart from
the standard continuous phase transition for single disease outbreaks, we
observe continuous transitions where both diseases must coexist, but also
discontinuous transitions are observed, where a finite fraction of the
population is already affected by both diseases at the threshold. All our
results are obtained in a mean field model using rate equations, but we argue
that they should hold also in more general frameworks.


  Motivated by modeling the dynamics of a population living in a flowing medium
where the environmental factors are random in space, we have studied an
asymmetric variant of the one-dimensional contact process, where the quenched
random reproduction rates are systematically greater in one direction than in
the opposite one. The spatial disorder turns out to be a relevant perturbation
but, according to results of Monte Carlo simulations, the behavior of the model
at the extinction transition is different from the (infinite randomness)
critical behavior of the disordered, symmetric contact process. Depending on
the strength $a$ of the asymmetry, the critical population drifts either with a
finite velocity or with an asymptotically vanishing velocity as $x(t)\sim
t^{\mu(a)}$, where $\mu(a)<1$. Dynamical quantities are non-self-averaging at
the extinction transition; the survival probability, for instance, shows
multiscaling, i.e. it is characterized by a broad spectrum of effective
exponents. For a sufficiently weak asymmetry, a Griffiths phase appears below
the extinction transition, where the survival probability decays as a
non-universal power of the time while, above the transition, another extended
phase emerges, where the front of the population advances anomalously with a
diffusion exponent continuously varying with the control parameter.


  Sequencing of pools of individuals (Pool-Seq) represents a reliable and cost-
effective approach for estimating genome-wide SNP and transposable element
insertion frequencies. However, Pool-Seq does not provide direct information on
haplotypes so that for example obtaining inversion frequencies has not been
possible until now. Here, we have developed a new set of diagnostic marker SNPs
for 7 cosmopolitan inversions in Drosophila melanogaster that can be used to
infer inversion frequencies from Pool-Seq data. We applied our novel marker set
to Pool-Seq data from an experimental evolution study and from North American
and Australian latitudinal clines. In the experimental evolution data, we find
evidence that positive selection has driven the frequencies of In(3R)C and
In(3R)Mo to increase over time. In the clinal data, we confirm the existence of
frequency clines for In(2L)t, In(3L)P and In(3R)Payne in both North America and
Australia and detect a previously unknown latitudinal cline for In(3R)Mo in
North America. The inversion markers developed here provide a versatile and
robust tool for characterizing inversion frequencies and their dynamics in
Pool- Seq data from diverse D. melanogaster populations.


  Adaptive landscape has been a fundamental concept in many branches of modern
biology since Wright's first proposition in 1932. Meanwhile, the general
existence of landscape remains controversial. The causes include the mixed uses
of different landscape definitions with their own different aims and
advantages. Sometimes the difficulty and the impossibility of the landscape
construction for complex models are also equated. To clarify these confusions,
based on a recent formulation of Wright's theory, the current authors construct
generalized adaptive landscape in a two-loci population model with non-gradient
dynamics, where the conventional gradient landscape does not exist. On the
generalized landscape, a population moves along an evolutionary trajectory
which always increases or conserves adaptiveness but does not necessarily
follow the steepest gradient direction. Comparisons of different aspects of
various landscapes lead to a conclusion that the generalized landscape is a
possible direction to continue the exploration of Wright's theory for complex
dynamics.


  Free energy landscapes decisively determine the progress of enzymatically
catalyzed reactions[1]. Time-resolved macromolecular crystallography unifies
transient-state kinetics with structure determination [2-4] because both can be
determined from the same set of X-ray data. We demonstrate here how barriers of
activation can be determined solely from five-dimensional crystallography [5].
Directly linking molecular structures with barriers of activation between them
allows for gaining insight into the structural nature of the barrier. We
analyze comprehensive time series of crystal-lographic data at 14 different
temperature settings and determine entropy and enthalpy contributions to the
barriers of activation. 100 years after the discovery of X-ray scattering, we
advance X-ray structure determination to a new frontier, the determination of
energy landscapes.


  In the present work we study and compare unfolding of a small protein,
chicken villin headpiece (HP-36) in two different aqueous binary mixtures,
namely water-ethanol (EtOH) and water-dimethyl sulphoxide (DMSO). In both the
binary mixtures, HP-36 is found to unfold (fully or partially, depending on the
mixture) under ambient conditions, that otherwise requires temperature as high
as ~600 K to denature in pure aqueous solvent. In all the cases, first step of
unfolding is found to be similar, i.e. separation of the cluster formed by
three hydrophobic (phenylalanine) residues, namely Phe-7, Phe-11 and Phe-18,
which constitute the hydrophobic core, thereby initiating melting of helix-2 of
the protein. Subsequent unfolding steps follow different paths in different
chemical environments. As both water-DMSO and water-ethanol show composition
dependent anomalies, so do the details of unfolding dynamics. With an increase
of co-solvent concentration different partially unfolded intermediates are
found to be formed in both the cases. This is reflected in a remarkable
non-monotonic composition dependence of several order parameters, including
fraction of native contacts and protein-solvent interaction energy. The
emergence of such partially unfolded states is particularly attributed to the
preferential solvation of the hydrophobic residues by the ethyl groups of
ethanol and methyl groups of DMSO. While in DMSO the protein gradually attains
a completely unfolded state at xDMSO=0.30, unfolding in water-ethanol appears
to be more complex and sensitive to solvent composition.


  A universal feature of the biochemistry of any living system is that all the
molecules and catalysts that are required for reactions of the system can be
built up from an available food source by repeated application of reactions
from within that system. RAF (reflexively autocatalytic and food-generated)
theory provides a formal way to study such processes. Beginning with Kauffman's
notion of "collectively autocatalytic sets", this theory has been further
developed over the last decade with the discovery of efficient algorithms and
new mathematical analysis. In this paper, we study how the behaviour of a
simple binary polymer model can be extended to models where the pattern of
catalysis more precisely reflects the ligation and cleavage reactions involved.
We find that certain properties of these models are similar to, and can be
accurately predicted from, the simple binary polymer model; however, other
properties lead to slightly different estimates. We also establish a number of
new results concerning the structure of RAFs in these systems.


  I present a stochastic population model that combines cooperative
interactions of the type often used in physics with the process of reproduction
and death familiar to biology, and I refer to reasons why such interlocking may
be of interest to both fields.


  Cooperativity arising from local interactions in equilibrium receptor systems
provides gain, but does not increase sensory performance, as measured by the
signal-to-noise ratio (SNR) due to a fundamental tradeoff between gain and
intrinsic noise. Here we allow sensing to be a nonequilibrium process and show
that energy dissipation cannot circumvent the fundamental tradeoff, so that SNR
is still optimal for independent receptors. For systems requiring high gain,
nonequilibrium 2D-coupled receptors maximize SNR, revealing a new design
principle for biological sensors.


  Synthetic lipid membranes in the absence of proteins can display quantized
conduction events for ions that are virtually indistinguishable from those of
protein channel. By indistinguishable we mean that one cannot decide based on
the current trace alone whether conductance events originate from a membrane,
which does or does not contain channel proteins. Additional evidence is
required to distinguish between the two cases, and it is not always certain
that such evidence can be provided. The phenomenological similarities are
striking and span a wide range of phenomena: The typical conductances are of
equal order and both lifetime distributions and current histograms are similar.
One finds conduction bursts, flickering, and multistep-conductance. Lipid
channels can be gated by voltage, and can be blocked by drugs. They respond to
changes in lateral membrane tension and temperature. Thus, they behave like
voltage-gated, temperature-gated and mechano-sensitive protein channels, or
like receptors. Lipid channels are remarkably under-appreciated. However, the
similarity between lipid and protein channels poses an eminent problem for the
interpretation of protein channel data. For instance, the Hodgkin-Huxley theory
for nerve pulse conduction requires a selective mechanism for the conduction of
sodium and potassium ions. To this end, the lipid membrane must act both as a
capacitor and as an insulator. Non-selective ion conductance by mechanisms
other than the gated protein-channels challenges the proposed mechanism for
pulse propagation. ... Some important questions arise: Are lipid and protein
channels similar due a common mechanism, or are these similarities fortuitous?
Is it possible that both phenomena are different aspects of the same
phenomenon? Are lipid and protein channels different at all? ... (abbreviated)


  Polypeptides can self-assemble into hierarchically organized fibrils
consisting of a stack of individually folded polypeptides driven together by
hydrophobic interaction. Using a coarse grained model, we systematically
studied this self-assembly as a function of temperature and hydrophobicity of
the residues on the outside of the building block. We find the self-assembly
can occur via two different pathways - a random aggregation-folding route, and
a templated-folding process - thus indicating a strong coupling between folding
and assembly. The simulation results can explain experimental evidence that
assembly through stacking of folded building blocks is rarely observed, at the
experimental concentrations. The model thus provides a generic picture of
hierarchical fibril formation.


  When four species compete stochastically in a cyclic way, the formation of
two teams of mutually neutral partners is observed. In this paper we study
through numerical simulations the extinction processes that can take place in
this system both in the well mixed case as well as on different types of
lattices. The different routes to extinction are revealed by the probability
distribution of the domination time, i.e. the time needed for one team to fully
occupy the system. If swapping is allowed between neutral partners, then the
probability distribution is dominated by very long-lived states where a few
very large domains persist, each domain being occupied by a mix of individuals
from species that form one of the teams. Many aspects of the possible
extinction scenarios are lost when only considering averaged quantities as for
example the mean domination time.


  Two synthetic gene circuits -- the genetic toggle switch and the
repressilator -- are analyzed quantitatively and discussed in the context of an
educational module on gene circuits and feedback that constitutes the final
topic of a year-long introductory physics sequence, aimed at biology and
premedical undergraduate students. The genetic toggle switch consists of two
genes, each of whose protein product represses the other's expression, while
the repressilator consists of three genes, each of whose protein product
represses the next gene's expression. Analytic, numerical, and electronic
treatments of the genetic toggle switch shows that this gene circuit realizes
bistability. A simplified treatment of the repressilator reveals that this
circuit can realize sustained oscillations. In both cases, a "phase diagram" is
obtained, that specifies the region of parameter space in which bistability or
oscillatory behavior, respectively, occurs.


  Protein interaction networks (PINs) are often used to "learn" new biological
function from their topology. Since current PINs are noisy, their computational
de-noising via link prediction (LP) could improve the learning accuracy. LP
uses the existing PIN topology to predict missing and spurious links. Many of
existing LP methods rely on shared immediate neighborhoods of the nodes to be
linked. As such, they have limitations. Thus, in order to comprehensively study
what are the topological properties of nodes in PINs that dictate whether the
nodes should be linked, we had to introduce novel sensitive LP measures that
overcome the limitations of the existing methods.
  We systematically evaluate the new and existing LP measures by introducing
"synthetic" noise to PINs and measuring how well the different measures
reconstruct the original PINs. Our main findings are: 1) LP measures that favor
nodes which are both "topologically similar" and have large shared extended
neighborhoods are superior; 2) using more network topology often though not
always improves LP accuracy; and 3) our new LP measures are superior to the
existing measures. After evaluating the different methods, we use them to
de-noise PINs. Importantly, we manage to improve biological correctness of the
PINs by de-noising them, with respect to "enrichment" of the predicted
interactions in Gene Ontology terms. Furthermore, we validate a statistically
significant portion of the predicted interactions in independent, external PIN
data sources.
  Software executables are freely available upon request.


  It is well known that, under suitable conditions, microRNAs are able to fine
tune the relative concentration of their targets to any desired value. We show
that this function is particularly effective when one of the targets is a
Transcription Factor (TF) which regulates the other targets. This combination
defines a new class of feed-forward loops (FFLs) in which the microRNA plays
the role of master regulator. Using both deterministic and stochastic equations
we show that these FFLs are indeed able not only to fine-tune the TF/target
ratio to any desired value as a function of the miRNA concentration but also,
thanks to the peculiar topology of the circuit, to ensures the stability of
this ratio against stochastic fluctuations. These two effects are due to the
interplay between the direct transcriptional regulation and the indirect
TF/Target interaction due to competition of TF and target for miRNA binding
(the so called "sponge effect"). We then perform a genome wide search of these
FFLs in the human regulatory network and show that they are characterizedby a
very peculiar enrichment pattern. In particular they are strongly enriched in
all the situations in which the TF and its target have to be precisely kept at
the same concentration notwithstanding the environmental noise. As an example
we discuss the FFL involving E2F1 as Transcription Factor, RB1 as target and
miR-17 family as master regulator. These FFLs ensure a tight control of the
E2F/RB ratio which in turns ensures the stability of the transition from the
G0/G1 to the S phase in quiescent cells.


  Motivation: Since susceptibility to diseases increases with age, studying
aging gains importance. Analyses of gene expression or sequence data, which
have been indispensable for investigating aging, have been limited to studying
genes and their protein products in isolation, ignoring their connectivities.
However, proteins function by interacting with other proteins, and this is
exactly what biological networks (BNs) model. Thus, analyzing the proteins' BN
topologies could contribute to understanding of aging. Current methods for
analyzing systems-level BNs deal with their static representations, even though
cells are dynamic. For this reason, and because different data types can give
complementary biological insights, we integrate current static BNs with
aging-related gene expression data to construct dynamic, age-specific BNs.
Then, we apply sensitive measures of topology to the dynamic BNs to study
cellular changes with age.
  Results: While global BN topologies do not significantly change with age,
local topologies of a number of genes do. We predict such genes as
aging-related. We demonstrate credibility of our predictions by: 1) observing
significant overlap between our predicted aging-related genes and "ground
truth" aging-related genes; 2) showing that our aging-related predictions group
by functions and diseases that are different than functions and diseases of
genes that are not predicted as aging-related; 3) observing significant overlap
between functions and diseases that are enriched in our aging-related
predictions and those that are enriched in "ground truth" aging-related data;
4) providing evidence that diseases which are enriched in our aging-related
predictions are linked to human aging; and 5) validating all of our
high-scoring novel predictions via manual literature search.


  We report SInC (SNV, Indel and CNV) simulator and read generator, an
open-source tool capable of simulating biological variants taking into account
a platform-specific error model. SInC is capable of simulating and generating
single- and paired-end reads with user-defined insert size with high efficiency
compared to the other existing tools. SInC, due to its multi-threaded
capability during read generation, has a low time footprint. SInC is currently
optimised to work in limited infrastructure setup and can efficiently exploit
the commonly used quad-core desktop architecture to simulate short sequence
reads with deep coverage for large genomes. Sinc can be downloaded from
https://sourceforge.net/projects/sincsimulator/.


  We present a continuum level analytical model of a droplet of active
contractile fluid consisting of filaments and motors. We calculate the steady
state flows that result from a splayed polarisation of the filaments. We
account for the interaction with an arbitrary external medium by imposing a
viscous friction at the fixed droplet boundary. We then show that the droplet
has non-zero force dipole and quadrupole moments, the latter of which is
essential for self-propelled motion of the droplet at low Reynolds' number.
Therefore, this calculation describes a simple mechanism for the motility of a
droplet of active contractile fluid embedded in a 3D environment, which is
relevant to cell migration in confinement (for example, embedded within a gel
or tissue). Our analytical results predict how the system depends on various
parameters such as the effective friction coefficient, the phenomenological
activity parameter and the splay of the imposed polarisation.


  Chronic diseases and cancer account for over 75 percent of healthcare costs
in the US. Increased prevention services and improved primary care are thought
to decrease costs. Current models for detecting changes in the health of
populations are cumbersome and expensive, and are not sensitive in the short
term. In this paper we model population health as a dynamical system to predict
the time evolution of the new diagnosis of chronic diseases and cancer. This
provides a reliable forecasting tool and a means of measuring short-term
changes in the health status of the population resulting from preventive care
programs. Twelve month forecasts of cancer and chronic populations were
accurate with errors lying between 3 percent and 6 percent. We confirmed what
other studies have demonstrated that diabetes patients are at increased cancer
risk but, interestingly, we also discovered that all of the studied chronic
conditions increased cancer risk just as diabetes did, and by a similar amount.
The model(i)yields a new metric for measuring performance of preventive and
clinical care programs that can provide timely feedback for quality improvement
programs;(ii)helps understand "savings" in the context of preventive care
programs and explains how they can be calculated in the short term, even though
they materialize only in the long term and(iii)provides an analytic tool and
metrics to infer correlations and derive insights on the effect of changes in
socio-economic factors affecting population health on improving health and
lowering costs of populations.


  Plasmodium falciparum 70 kDa heat shock proteins (PfHsp70s) are expressed at
all stages of the pathogenic erythrocytic phase of the malaria parasite
lifecycle. There are six PfHsp70s, all of which have orthologues in other
plasmodial species, except for PfHsp70-x which is unique to P. falciparum. This
paper highlights a number of original results obtained by a detailed
bioinformatics analysis of the protein. Large scale sequence analysis indicated
the presence of an extended transit peptide sequence of PfHsp70-x which
potentially directs it to the endoplasmic reticulum (ER). Further analysis
showed that PfHsp70-x does not have an ER-retention sequence, suggesting that
the protein transits through the ER and is secreted into the parasitophorous
vacuole (PV) or beyond into the erythrocyte cytosol. These results are
consistent with experimental findings. Next, possible interactions between
PfHsp70-x and exported P. falciparum Hsp40s or host erythrocyte DnaJs were
interrogated by modeling and docking. Docking results indicated that
interaction between PfHsp70-x and each of the Hsp40s, regardless of biological
feasibility, seems equally likely. This suggests that J domain might not
provide the specificity in the formation of unique Hsp70-Hsp40 complexes, but
that the specificity might be provided by other domains of Hsp40s. By studying
different structural conformations of PfHsp70-x, it was shown that Hsp40s can
only bind when PfHsp70-x is in a certain conformation. Additionally, this work
highlighted the possible dependence of the substrate binding domain residues on
the orientation of the {\alpha}-helical lid for formation of the substrate
binding pocket.


  Motivation: Illumina Sequencing data can provide high coverage of a genome by
relatively short (100 bp150 bp) reads at a low cost. Our goal is to produce
trimmed and error-corrected reads to improve genome assemblies. Our error
correction procedure aims at producing a set of error-corrected reads (1)
minimizing the number of distinct false k-mers, i.e. that are not present in
the genome, in the set of reads and (2) maximizing the number that are true,
i.e. that are present in the genome. Because coverage of a genome by Illumina
reads varies greatly from point to point, we cannot simply eliminate k-mers
that occur rarely.
  Results: Our software, called QuorUM, provides reasonably accurate correction
and is suitable for large data sets (1 billion bases checked and corrected per
day per core).
  Availability: QuorUM is distributed as an independent software package and as
a module of the MaSuRCA assembly software. Both are available under the GPL
open source license at http://www.genome.umd.edu.
  Contact: gmarcais@umd.edu


  The Kruskal-Wallis test is a popular non-parametric statistical test for
identifying expression quantitative trait loci (eQTLs) from genome-wide data
due to its robustness against variations in the underlying genetic model and
expression trait distribution, but testing billions of marker-trait
combinations one-by-one can become computationally prohibitive. We developed
kruX, an algorithm implemented in Matlab, Python and R that uses matrix
multiplications to simultaneously calculate the Kruskal-Wallis test statistic
for several millions of marker-trait combinations at once. KruX is more than
ten thousand times faster than computing associations one-by-one on a typical
human dataset. We used kruX and a dataset of more than 500k SNPs and 20k
expression traits measured in 102 human blood samples to compare eQTLs detected
by the Kruskal-Wallis test to eQTLs detected by the parametric ANOVA and linear
model methods. We found that the Kruskal-Wallis test is more robust against
data outliers and heterogeneous genotype group sizes and detects a higher
proportion of non-linear associations, but is more conservative for calling
additive linear associations. In summary, kruX enables the use of robust
non-parametric methods for massive eQTL mapping without the need for a
high-performance computing infrastructure.


  Disease-causing mutations usually change the interacting partners of mutant
proteins. In this article, we propose that the biological consequences of
mutation are directly related to the alteration of corresponding protein
protein interaction networks (PPIN). Mutation of Huntingtin (HTT) which causes
Huntington's disease (HD) and mutations to TP53 which is associated with
different cancers are studied as two example cases. We construct the PPIN of
wild type and mutant proteins separately and identify the structural modules of
each of the networks. The functional role of these modules are then assessed by
Gene Ontology (GO) enrichment analysis for biological processes (BPs). We find
that a large number of significantly enriched (p<0.0001) GO terms in mutant
PPIN were absent in the wild type PPIN indicating the gain of BPs due to
mutation. Similarly some of the GO terms enriched in wild type PPIN cease to
exist in the modules of mutant PPIN, representing the loss. GO terms common in
modules of mutant and wild type networks indicate both loss and gain of BPs. We
further assign relevant biological function(s) to each module by classifying
the enriched GO terms associated with it. It turns out that most of these
biological functions in HTT networks are already known to be altered in HD and
those of TP53 networks are altered in cancers. We argue that gain of BPs, and
the corresponding biological functions, are due to new interacting partners
acquired by mutant proteins. The methodology we adopt here could be applied to
genetic diseases where mutations alter the ability of the protein to interact
with other proteins.