  Evolutionary models for populations of constant size are frequently studied
using the Moran model, the Wright-Fisher model, or their diffusion limits. When
evolution is neutral, a random genealogy given through Kingman's coalescent is
used in order to understand basic properties of such models. Here, we address
the use of a genealogical perspective for models with weak frequency-dependent
selection, i.e. N s =: {\alpha} is small, and s is the fitness advantage of a
fit individual and N is the population size. When computing fixation
probabilities, this leads either to the approach proposed by Rousset (2003),
who argues how to use the Kingman's coalescent for weak selection, or to
extensions of the ancestral selection graph of Neuhauser and Krone (1997) and
Neuhauser (1999). As an application, we re-derive the one-third law of
evolutionary game theory (Nowak et al., 2004). In addition, we provide the
approximate distribution of the genealogical distance of two randomly sampled
individuals under linear frequency-dependence.


  We present force-clamp data on the collapse of ubiquitin polyproteins in
response to a quench in the force. These nonequilibrium trajectories are
analyzed using a general method based on a diffusive assumption of the
end-to-end length to reconstruct a downhill free energy profile at 5pN and an
energy plateau at 10pN with a slow diffusion coefficient on the order
of~100nm^2/s. The shape of the free energy and its linear scaling with the
protein length give validity to a physical model for the collapse. However, the
length independent diffusion coefficient suggests that internal rather than
viscous friction dominates and thermal noise is needed to capture the
variability in the measured times to collapse.


  Systems biology approaches combining theoretical modeling with experiments
have been singularly successful in uncovering novel features of cellular
phenomena. One such feature is that of binary gene expression in which the
expression level is either low or high, i.e., digital in nature. This gives
rise to two distinct subpopulations in a population of genetically identical
cells. The fraction of cells in the high expression state is raised as the
strength of the inducing signal is increased indicating that the response is
not graded. In this review, we discuss the possible origins of binary gene
expression with emphasis on three principal mechanisms: purely stochastic,
positive feedback-based and emergent bistability. In the latter case, two
stable expression states are obtained due to an autoregulatory positive
feedback loop in protein synthesis along with cell growth retardation by the
proteins synthesized. The theoretical foundations of the observed phenomena are
described in each case.


  Hemoglobin exhibits allosteric structural changes upon ligand binding due to
the dynamic interactions between the ligand binding sites, the amino acids
residues and some other solutes present under physiological conditions. In the
present study, the dynamical and quaternary structural changes occurring in two
unligated (deoxy-) T structures, and two fully ligated (oxy-) R, R2 structures
of adult human hemoglobin were investigated with molecular dynamics. It is
shown that, in the sub-microsecond time scale, there is no marked difference in
the global dynamics of the amino acids residues in both the oxy- and the deoxy-
forms of the individual structures. In addition, the R, R2 are relatively
stable and do not present quaternary conformational changes within the time
scale of our simulations while the T structure is dynamically more flexible and
exhibited the T\rightarrow R quaternary conformational transition, which is
propagated by the relative rotation of the residues at the {\alpha}1{\beta}2
and {\alpha}2{\beta}1 interface.


  The analysis of the evolutionary dynamics of a population with many
polymorphic loci is challenging since a large number of possible genotypes
needs to be tracked. In the absence of analytical solutions, forward computer
simulations are an important tool in multi-locus population genetics. The run
time of standard algorithms to simulate sexual populations increases as 8^L
with the number L of loci, or with the square of the population size N. We have
developed algorithms that allow to simulate large populations with a run-time
that scales as 3^L. The algorithm is based on an analog of the Fast-Fourier
Transform (FFT) and allows for arbitrary fitness functions (i.e. any epistasis)
and genetic maps. The algorithm is implemented as a collection of C++ classes
and a Python interface.


  Ion channels form pores across the lipid bilayer, selectively allowing
inorganic ions to cross the membrane down their electrochemical gradient. While
the study of ion desolvation free-energies have attracted much attention, the
role of water inside the pore is less clear. Here, molecular dynamics
simulations of a reduced model of the KcsA selectivity filter indicate that the
equilibrium position of Na+, but not of K+, is strongly influenced by confined
water. The latter forms a stable complex with Na+, moving the equilibrium
position of the ion to the plane of the backbone carbonyls. Almost at the
centre of the binding site, the water molecule is trapped by favorable
electrostatic interactions and backbone hydrogen-bonds. In the absence of
confined water the equilibrium position of both Na+ and K+ is identical. Our
observations strongly suggest a previously unnoticed active role of confined
water in the selectivity mechanism of ion channels.


  During the course of evolution, an organism's genome can undergo changes that
affect the large-scale structure of the genome. These changes include gene
gain, loss, duplication, chromosome fusion, fission, and rearrangement. When
gene gain and loss occurs in addition to other types of rearrangement,
breakpoints of rearrangement can exist that are only detectable by comparison
of three or more genomes. An arbitrarily large number of these "hidden"
breakpoints can exist among genomes that exhibit no rearrangements in pairwise
comparisons.
  We present an extension of the multichromosomal breakpoint median problem to
genomes that have undergone gene gain and loss. We then demonstrate that the
median distance among three genomes can be used to calculate a lower bound on
the number of hidden breakpoints present. We provide an implementation of this
calculation including the median distance, along with some practical
improvements on the time complexity of the underlying algorithm.
  We apply our approach to measure the abundance of hidden breakpoints in
simulated data sets under a wide range of evolutionary scenarios. We
demonstrate that in simulations the hidden breakpoint counts depend strongly on
relative rates of inversion and gene gain/loss. Finally we apply current
multiple genome aligners to the simulated genomes, and show that all aligners
introduce a high degree of error in hidden breakpoint counts, and that this
error grows with evolutionary distance in the simulation. Our results suggest
that hidden breakpoint error may be pervasive in genome alignments.


  ARISA (Automated Ribosomal Intergenic Spacer Analysis) is a low-cost
technique that allows for the rapid comparison of different microbial
environments. In this study, we asked if a set of ARISA profiles can
distinguish human microbial environments from one another with the same
accuracy as results generated from 454 high throughput DNA sequencing. Using a
set of human microbial communities where the sequencing results cluster by
subject, we tested how choices made during ARISA data processing influence
clustering. We found that choice of clustering methods had a profound effect
with Ward's clustering generating profiles the most similar to 454 sequencing.
Factors such as bin size, using presence or absence calls and technical
replicate manipulation had a negligible effect on clustering. In fact, no
established bin sizing method reported in the literature performed
significantly different results than simply picking bin intervals at random. We
conclude that in an analysis of ARISA data from an ecosystem of sufficient
complexity to saturate bins, a careful choice of clustering algorithm is
essential whereas differing strategies for choosing bins are likely to have a
much less pronounced effect on the outcome of the analysis. As a tool for
distinguishing complex microbial communities, ARISA closely approximates the
results obtained from DNA sequencing at a fraction of the cost; however ARISA
fails to reproduce the sequencing results perfectly.


  It has remained an open question for some time whether, given a set of not
necessarily binary (i.e. "nonbinary") trees T on a set of taxa X, it is
possible to determine in time f(r).poly(m) whether there exists a phylogenetic
network that displays all the trees in T, where r refers to the reticulation
number of the network and m=|X|+|T|. Here we show that this holds if one or
both of the following conditions holds: (1) |T| is bounded by a function of r;
(2) the maximum degree of the nodes in T is bounded by a function of r. These
sufficient conditions absorb and significantly extend known special cases,
namely when all the trees in T are binary, or T contains exactly two nonbinary
trees. We believe this result is an important step towards settling the issue
for an arbitrarily large and complex set of nonbinary trees. For completeness
we show that the problem is certainly solveable in polynomial time.


  The ability of information processing in biologically motivated Boolean
networks is of interest in recent information theoretic research. One measure
to quantify this ability is the well known mutual information. Using Fourier
analysis we show that canalizing functions maximize the mutual information
between an input variable and the outcome of the function. We proof our result
for Boolean functions with uniform distributed as well as product distributed
input variables.


  We propose a specialized string kernel for small bio-molecules, peptides and
pseudo-sequences of binding interfaces. The kernel incorporates
physico-chemical properties of amino acids and elegantly generalize eight
kernels, such as the Oligo, the Weighted Degree, the Blended Spectrum, and the
Radial Basis Function. We provide a low complexity dynamic programming
algorithm for the exact computation of the kernel and a linear time algorithm
for it's approximation. Combined with kernel ridge regression and SupCK, a
novel binding pocket kernel, the proposed kernel yields biologically relevant
and good prediction accuracy on the PepX database. For the first time, a
machine learning predictor is capable of accurately predicting the binding
affinity of any peptide to any protein. The method was also applied to both
single-target and pan-specific Major Histocompatibility Complex class II
benchmark datasets and three Quantitative Structure Affinity Model benchmark
datasets.
  On all benchmarks, our method significantly (p-value < 0.057) outperforms the
current state-of-the-art methods at predicting peptide-protein binding
affinities. The proposed approach is flexible and can be applied to predict any
quantitative biological activity. The method should be of value to a large
segment of the research community with the potential to accelerate
peptide-based drug and vaccine development.


  We introduce a lattice model of protein conformations which is able to
reproduce second structures of proteins (alpha--helices and beta--sheets). This
model is based on the following two main ideas. First, we model backbone parts
of amino acid residues in a peptide chain by edges in the cubic lattice which
are not parallel to the coordinate axes. Second, we describe possible contacts
of amino acid residues using a discrete model of the Ramachandran plot.
  This model allows to describe hydrogen bonds between the residues in the
backbone of the peptide chain. In particular the lattice secondary structures
have the correct structure of hydrogen bonds. We also take into account the
side chains of amino acid residues and their interaction.
  The expression for the energy of conformation of a lattice protein which
contains contributions from hydrogen bonds in the backbone of the peptide chain
and from interaction of the side chains is proposed. The lattice secondary
structures are local minima of the introduced energy.


  In comparative genomics, the rearrangement distance between two genomes
(equal the minimal number of genome rearrangements required to transform them
into a single genome) is often used for measuring their evolutionary
remoteness. Generalization of this measure to three genomes is known as the
median score (while a resulting genome is called median genome). In contrast to
the rearrangement distance between two genomes which can be computed in linear
time, computing the median score for three genomes is NP-hard. This inspires a
quest for simpler and faster approximations for the median score, the most
natural of which appears to be the halved sum of pairwise distances which in
fact represents a lower bound for the median score.
  In this work, we study relationship and interplay of pairwise distances
between three genomes and their median score under the model of
Double-Cut-and-Join (DCJ) rearrangements. Most remarkably we show that while a
rearrangement may change the sum of pairwise distances by at most 2 (and thus
change the lower bound by at most 1), even the most "powerful" rearrangements
in this respect that increase the lower bound by 1 (by moving one genome
farther away from each of the other two genomes), which we call strong, do not
necessarily affect the median score. This observation implies that the two
measures are not as well-correlated as one's intuition may suggest.
  We further prove that the median score attains the lower bound exactly on the
triples of genomes that can be obtained from a single genome with strong
rearrangements. While the sum of pairwise distances with the factor 2/3
represents an upper bound for the median score, its tightness remains unclear.
Nonetheless, we show that the difference of the median score and its lower
bound is not bounded by a constant.


  Long noncoding RNAs (lncRNAs) are a recently discovered class of non-protein
coding RNAs which have now increasingly been shown to be involved in a wide
variety of biological processes as regulatory molecules. Little is known
regarding the regulatory interactions between noncoding RNA classes. Recent
reports have suggested that lncRNAs could potentially interact with other
noncoding RNAs including miroRNAs (miRNAs) and modulate their regulatory role
through interactions. We hypothesized that long noncoding RNAs could
participate as a layer of regulatory interactions with miRNAs. The availability
of genome-scale datasets for argonaute targets across human transcriptome has
prompted us to reconstruct a genome-scale network of interactions between
miRNAs and lncRNAs.
  We used well characterized experimental
Photoactivatable-Ribonucleoside-Enhanced Crosslinking and Immunoprecipitation
(PAR-CLIP) datasets and the recent genome-wide annotations for lncRNAs in
public domain to construct a comprehensive transcriptome-wide map of miRNA
regulatory elements. Comparative analysis revealed many of the miRNAs could
target long noncoding RNAs, apart from the coding transcripts thus
participating in a novel layer of regulatory interactions between noncoding RNA
classes. We also find the miRNA regulatory elements have a positional
preference, clustering towards the 3' and 5' ends of the long noncoding
transcripts. We also further reconstruct a genome-wide map of miRNA
interactions with lncRNAs as well as messenger RNAs.
  This analysis suggests widespread regulatory interactions between noncoding
RNAs classes and suggests a novel functional role for lncRNAs. We also present
the first transcriptome scale study on lncRNA-miRNA interactions and the first
report of a genome-scale reconstruction of a noncoding RNA regulatory
interactome involving lncRNAs.


  Embryonic development is driven by spatial patterns of gene expression that
determine the fate of each cell in the embryo. While gene expression is often
highly erratic, embryonic development is usually exceedingly precise. In
particular, gene expression boundaries are robust not only against intrinsic
noise from gene expression and protein diffusion, but also against
embryo-to-embryo variations in the morphogen gradients, which provide
positional information to the differentiating cells. How development is robust
against intra- and inter-embryonic variations is not understood. A common motif
in the gene regulation networks that control embryonic development is mutual
repression between pairs of genes. To assess the role of mutual repression in
the robust formation of gene expression patterns, we have performed large-scale
stochastic simulations of a minimal model of two mutually repressing gap genes
in Drosophila, hunchback (hb) and knirps (kni). Our model includes not only
mutual repression between hb and kni, but also the stochastic and cooperative
activation of hb by the anterior morphogen Bicoid (Bcd) and of kni by the
posterior morphogen Caudal (Cad), as well as the diffusion of Hb and Kni. Our
analysis reveals that mutual repression can markedly increase the steepness and
precision of the gap gene expression boundaries. In contrast to spatial
averaging and cooperative gene activation, mutual repression thus allows for
gene-expression boundaries that are both steep and precise. Moreover, mutual
repression dramatically enhances their robustness against embryo-to-embryo
variations in the morphogen levels. Finally, our simulations reveal that gap
protein diffusion plays a critical role not only in reducing the width of gap
gene expression boundaries via spatial averaging, but also in repairing
patterning errors that could arise due to the bistability induced by mutual
repression.


  We report spin-selective tunneling of electrons along natural and artificial
double-stranded DNA (dsDNA) sandwiched by nonmagnetic leads. The results reveal
that the spin polarization strongly depends on the dsDNA sequence and is
dominated by its end segment. Both genomic and artificial dsDNA could be
efficient spin filters. The spin-filtering effects are sensitive to point
mutation which occurs in the end segment. These results are in good agreement
with recent experiments and are robust against various types of disorder, and
could help for designing DNA-based spintronic devices.


  Two methods of exploratory spatial data analysis (ESDA), analysis of spatial
heterogeneity and dependence (auto-correlation), - - were applied to the cover
patterns from ten paramo localities in the Central and Eastern cordilleras of
Colombia. Among the localities studied, the high montane region of the Serrania
de Perija, the paramo region of the Los Nevados National Park, and the paramo
region under management of CORPOGUAVIO showed a good state of conservation and
satisfactory level of connectivity among patches. Anthropic intervention,
expansion of potato cultivation and other changes in land use mainly in the
paramos of the Cordillera Oriental have alteredd the distribution patterns of
natural vegetation types and resulted in significant losses of connectively.
Under the circumstancesany attempt to restore the original conditions will be
difficult and expensive. This problem is most acute in the paramo localities of
La Rusia, Bel\'en, Guantiva (towards the eastern slopes) Telecom, Merchan
Tablazo, the northeastern sector of the paramo of Sumapaz and the eastern
slopes of the mountains of Carmen de Carupa and surrounding areas. For the
mountains of Guanacaste and Las Delicias (areas influenced by the volcano
Purac\'e) the problem is more acute towards the northwestern sector. The
results of the present study can be applied to ecological restoration and
conservation as they provide detailed information about the distribution of
vegetation types andhabitatconnectivity, both indicators of ecosystem health.


  We estimated the aerial biomass and stored carbon in twelve types of forests
in the department of C\'ordoba with annual rain fall ranging from 3000 mm
(super humid climates) to 1300 mm (semi-humid climate). Biomass was estimated
based on structural aspects of the vegetation (diameter at breast height, total
height, and wood specific weight). We tested nine allometric equations for
tropical forests available in the literature and selected those proposed by
Chave et al (2005) that are specific for humid and dry forests. Carbon content
in trees was measured in four tissues (stem, branch, bark, and leafs) through
an automated dry combustion method, which estimates the percentage of carbon in
a sample of known weight.


  In spring turnip rape (Brassica rapa L. spp. oleifera) the most promising F1
hybrid system would be the Ogu-INRA CMS/Rf system. A Kosena fertility restorer
gene Rfk1, homologue of the Ogura restorer gene Rfo, was successfully
transferred from oilseed rape into turnip rape and that restored the fertility
in female lines carrying Ogura cms. The trait was, however, unstable in
subsequent generations. The physical localization of the radish chromosomal
region carrying the Rfk1 gene was investigated using 8 GISH (genomic in situ
hybridization) and BAC-FISH (bacterial artificial chromosome fluorescence in
situ hybridization) methods. The metaphase chromosomes were hybridized using
radish DNA as the genomic probe and BAC64 probe, which is linked with Rfo gene.
Both probes showed a signal in the chromosome spreads of the restorer line
4021-2 Rfk of turnip rape but not in the negative control line 4021B. The GISH
analyses clearly showed that the turnip rape restorer plants were either
monosomic (2n=2x=20+1R) or disomic (2n=2x=20+2R) addition lines with one or two
copies of a single alien chromosome region originating from radish. In the
BAC-FISH analysis, double dot signals were detected in sub-terminal parts of
the radish chromosome arms showing that the fertility restorer gene Rfk1 was
located in this additional radish chromosome. Detected disomic addition lines
were found to be unstable for turnip rape hybrid production. Using the BAC-FISH
analysis, weak signals were sometimes visible in two chromosomes of turnip rape
and a homologous region of Rfk1 in chromosome 9 of the B. rapa A genome was
verified with BLAST analysis. In the future this homologous area in A genome
could be substituted with radish chromosome area carrying the Rfk1 gene.


  The evolution of cooperation often depends upon population structure, yet
nearly all models of cooperation implicitly assume that this structure remains
static. This is a simplifying assumption, because most organisms possess
genetic traits that affect their population structure to some degree. These
traits, such as a group size preference, affect the relatedness of interacting
individuals and hence the opportunity for kin or group selection. We argue that
models that do not explicitly consider their evolution cannot provide a
satisfactory account of the origin of cooperation, because they cannot explain
how the prerequisite population structures arise. Here, we consider the
concurrent evolution of genetic traits that affect population structure, with
those that affect social behavior. We show that not only does population
structure drive social evolution, as in previous models, but that the
opportunity for cooperation can in turn drive the creation of population
structures that support it. This occurs through the generation of linkage
disequilibrium between socio-behavioral and population-structuring traits, such
that direct kin selection on social behavior creates indirect selection
pressure on population structure. We illustrate our argument with a model of
the concurrent evolution of group size preference and social behavior.


  Selection on the level of loosely associated groups has been suggested as a
route towards the evolution of cooperation between individuals and the
subsequent formation of higher-level biological entities. Such group selection
explanations remain problematic, however, due to the narrow range of parameters
under which they can overturn within-group selection that favours selfish
behaviour. In principle, individual selection could act on such parameters so
as to strengthen the force of between-group selection and hence increase
cooperation and individual fitness, as illustrated in our previous work.
However, such a process cannot operate in parameter regions where group
selection effects are totally absent, since there would be no selective
gradient to follow. One key parameter, which when increased often rapidly
causes group selection effects to tend to zero, is initial group size, for when
groups are formed randomly then even moderately sized groups lack significant
variance in their composition. However, the consequent restriction of any group
selection effect to small sized groups is derived from models that assume
selfish types will competitively exclude their more cooperative counterparts at
within-group equilibrium. In such cases, diversity in the migrant pool can tend
to zero and accordingly variance in group composition cannot be generated. In
contrast, we show that if within-group dynamics lead to a stable coexistence of
selfish and cooperative types, then the range of group sizes showing some
effect of group selection is much larger.


  Multilevel selection and the evolution of cooperation are fundamental to the
formation of higher-level organisation and the evolution of biocomplexity, but
such notions are controversial and poorly understood in natural populations.
The theoretic principles of group selection are well developed in idealised
models where a population is neatly divided into multiple semi-isolated
sub-populations. But since such models can be explained by individual selection
given the localised frequency-dependent effects involved, some argue that the
group selection concepts offered are, even in the idealised case, redundant and
that in natural conditions where groups are not well-defined that a group
selection framework is entirely inapplicable. This does not necessarily mean,
however, that a natural population is not subject to some interesting localised
frequency-dependent effects -- but how could we formally quantify this under
realistic conditions? Here we focus on the presence of a Simpson's Paradox
where, although the local proportion of cooperators decreases at all locations,
the global proportion of cooperators increases. We illustrate this principle in
a simple individual-based model of bacterial biofilm growth and discuss various
complicating factors in moving from theory to practice of measuring group
selection.


  We investigate the dynamics of a broad class of stochastic copying processes
on a network that includes examples from population genetics
(spatially-structured Wright-Fisher models), ecology (Hubbell-type models),
linguistics (the utterance selection model) and opinion dynamics (the voter
model) as special cases. These models all have absorbing states of fixation
where all the nodes are in the same state. Earlier studies of these models
showed that the mean time when this occurs can be made to grow as different
powers of the network size by varying the the degree distribution of the
network. Here we demonstrate that this effect can also arise if one varies the
asymmetry of the copying dynamics whilst holding the degree distribution
constant. In particular, we show that the mean time to fixation can be
accelerated even on homogeneous networks when certain nodes are very much more
likely to be copied from than copied to. We further show that there is a
complex interplay between degree distribution and asymmetry when they may
co-vary; and that the results are robust to correlations in the network or the
initial condition.


  Transcriptional regulatory network inference methods have been studied for
years. Most of them relie on complex mathematical and algorithmic concepts,
making them hard to adapt, re-implement or integrate with other methods. To
address this problem, we introduce a novel method based on a minimal
statistical model for observing transcriptional regulatory interactions in
noisy expression data, which assumes that transcription factors (TFs) and their
targets are both differentially expressed in a gene-specific, critical sample
contrast, as measured by repeated two-way t-tests. This method is conceptually
simple and easy to implement and integrate in any statistical software
environment. Benchmarking on standard E. coli and yeast reference datasets
showed that it performs equally well as the best existing methods. Analysis of
the predicted interactions suggested that it works best to infer
context-specific TF-target interactions which only co-express locally. We
confirmed this hypothesis on a dataset of more than 1,000 normal human tissue
samples, where we found that our method predicts highly tissue-specific and
functionally relevant interactions, whereas a global co-expression method only
associates general TFs to non-specific biological processes.


  We analyze the social mechanisms that shape the popularity rise and fall of
the names given to newborn babies. During the initial stage, popularity
increases by imitation. As the people with the same name grow in number,
however, its usage is inhibited and eventually decays. This process mirrors the
dynamics of fashion fads. An activator-inhibitor dynamical model for the
interplay of the population bearing a name and the expecting couples wishing to
give it to their children provides a satisfactory explanation of historical
data from the Canadian province of Quebec during the twentieth century.


  The development of an aptamer-based impedimetric sensor for typing of
bacteria (AIST-B) is presented. Highly specific DNA aptamers to Salmonella
enteritidis were selected via Cell-SELEX technique. Twelve rounds of selection
were performed; each comprises a positive selection step against S. enteritidis
and a negative selection step against a mixture of related pathogens, including
Salmonella typhimurium, Escherichia coli, Staphylococcus aureus, Pseudomonas
aeruginosa, and Citrobacter freundii to ensure the species-specificity of the
selected aptamers. After sequencing of the pool showing the highest binding
affinity to S. enteritidis, a DNA aptamer sequence of high affinity to the
bacteria was integrated into an impedimetric sensor via self-assembly onto a
gold nanoparticles-modified screen-printed carbon electrode (GNPs-SPCE).
Remarkably, this aptasensor is highly selective and can successfully detect S.
enteritidis down to 600 CFU 1/mL (equivalent to 18 CFU in 30 \muL assay volume)
in 10 minutes and distinguish it from other salmonella species, including S.
typhimurium and S. choleraesuis. This report is envisaged to open a new venue
for the aptamer-based typing of a variety of microorganisms using a rapid,
economic, label-free electrochemical platform.


  Many biological characteristics of evolutionary interest are not scalar
variables but continuous functions. Here we use phylogenetic Gaussian process
regression to model the evolution of simulated function-valued traits. Given
function-valued data only from the tips of an evolutionary tree and utilising
independent principal component analysis (IPCA) as a method for dimension
reduction, we construct distributional estimates of ancestral function-valued
traits, and estimate parameters describing their evolutionary dynamics.


  Populations of widely distributed species often encounter and adapt to
specific environmental conditions. However, comprehensive characterization of
the genetic basis of adaptation is demanding, requiring genome-wide genotype
data, multiple sampled populations, and a good understanding of population
structure. We have used environmental and high-density genotype data to
describe the genetic basis of local adaptation in 21 populations of teosinte,
the wild ancestor of maize. We found that altitude, dispersal events and
admixture among subspecies formed a complex hierarchical genetic structure
within teosinte. Patterns of linkage disequilibrium revealed four mega-base
scale inversions that segregated among populations and had altitudinal clines.
Based on patterns of differentiation and correlation with environmental
variation, inversions and nongenic regions play an important role in local
adaptation of teosinte. Further, we note that strongly differentiated
individual populations can bias the identification of adaptive loci. The role
of inversions in local adaptation has been predicted by theory and requires
attention as genome-wide data become available for additional plant species.
These results also suggest a potentially important role for noncoding
variation, especially in large plant genomes in which the gene space represents
a fraction of the entire genome.


  It has been suggested that a Random Tree Puzzle (RTP) process leads to a
Yule-Harding (YH) distribution, when the number of taxa becomes large. In this
study, we formalize this conjecture, and we prove that the two tree
distributions converge for two particular properties, which suggests that the
conjecture may be true. However, we present evidence that, while the two
distributions are close, the RTP appears to converge on a different
distribution than does the YH.


  The ways in which natural selection can allow the proliferation of
cooperative behavior have long been seen as a central problem in evolutionary
biology. Most of the literature has focused on interactions between pairs of
individuals and on linear public goods games. This emphasis led to the
conclusion that even modest levels of migration would pose a serious problem to
the spread of altruism in group structured populations. Here we challenge this
conclusion, by analyzing evolution in a framework which allows for complex
group interactions and random migration among groups. We conclude that
contingent forms of strong altruism can spread when rare under realistic group
sizes and levels of migration. Our analysis combines group-centric and
gene-centric perspectives, allows for arbitrary strength of selection, and
leads to extensions of Hamilton's rule for the spread of altruistic alleles,
applicable under broad conditions.


  Motivated by questions in mass-action kinetics, we introduce the notion of
vertexical family of differential inclusions. Defined on open hypercubes, these
families are characterized by particular good behavior under projection maps.
The motivating examples are certain families of reaction networks -- including
reversible, weakly reversible, endotactic, and strongly endotactic reaction
networks -- that give rise to vertexical families of mass-action differential
inclusions. We prove that vertexical families are amenable to structural
induction. Consequently, a trajectory of a vertexical family approaches the
boundary if and only if either the trajectory approaches a vertex of the
hypercube, or a trajectory in a lower-dimensional member of the family
approaches the boundary. With this technology, we make progress on the global
attractor conjecture, a central open problem concerning mass-action kinetics
systems. Additionally, we phrase mass-action kinetics as a functor on reaction
networks with variable rates.


  We present a novel formulation for biochemical reaction networks in the
context of signal transduction. The model consists of input-output transfer
functions, which are derived from differential equations, using stable
equilibria. We select a set of 'source' species, which receive input signals.
Signals are transmitted to all other species in the system (the 'target'
species) with a specific delay and transmission strength. The delay is computed
as the maximal reaction time until a stable equilibrium for the target species
is reached, in the context of all other reactions in the system. The
transmission strength is the concentration change of the target species. The
computed input-output transfer functions can be stored in a matrix, fitted with
parameters, and recalled to build discrete dynamical models. By separating
reaction time and concentration we can greatly simplify the model,
circumventing typical problems of complex dynamical systems. The transfer
function transformation can be applied to mass-action kinetic models of signal
transduction. The paper shows that this approach yields significant insight,
while remaining an executable dynamical model for signal transduction. In
particular we can deconstruct the complex system into local transfer functions
between individual species. As an example, we examine modularity and signal
integration using a published model of striatal neural plasticity. The modules
that emerge correspond to a known biological distinction between
calcium-dependent and cAMP-dependent pathways. We also found that overall
interconnectedness depends on the magnitude of input, with high connectivity at
low input and less connectivity at moderate to high input. This general result,
which directly follows from the properties of individual transfer functions,
contradicts notions of ubiquitous complexity by showing input-dependent signal
transmission inactivation.


  We consider mutual information between release times and capture times for a
set of M identical quanta traveling independently from a source to a target.
The quanta are immediately captured upon arrival, first-passage times are
assumed independent and identically distributed and the quantum emission times
are constrained by a deadline. The primary application area is intended to be
inter/intracellular molecular signaling in biological systems whereby an
organelle, cell or group of cells must deliver some message (such as
transcription or developmental instructions) over distance with reasonable
certainty to another organelles, cells or group of cells. However, the model
can also be applied to communications systems wherein indistinguishable signals
have random transit latencies.


  The question of Jewish ancestry has been the subject of controversy for over
two centuries and has yet to be resolved. The "Rhineland Hypothesis" proposes
that Eastern European Jews emerged from a small group of German Jews who
migrated eastward and expanded rapidly. Alternatively, the "Khazarian
Hypothesis" suggests that Eastern European descended from Judean tribes who
joined the Khazars, an amalgam of Turkic clans that settled the Caucasus in the
early centuries CE and converted to Judaism in the 8th century. The Judaized
Empire was continuously reinforced with Mesopotamian and Greco-Roman Jews until
the 13th century. Following the collapse of their empire, the Judeo-Khazars
fled to Eastern Europe. The rise of European Jewry is therefore explained by
the contribution of the Judeo-Khazars. Thus far, however, their contribution
has been estimated only empirically; the absence of genome-wide data from
Caucasus populations precluded testing the Khazarian Hypothesis. Recent
sequencing of modern Caucasus populations prompted us to revisit the Khazarian
Hypothesis and compare it with the Rhineland Hypothesis. We applied a wide
range of population genetic analyses - including principal component,
biogeographical origin, admixture, identity by descent, allele sharing
distance, and uniparental analyses - to compare these two hypotheses. Our
findings support the Khazarian Hypothesis and portray the European Jewish
genome as a mosaic of Caucasus, European, and Semitic ancestries, thereby
consolidating previous contradictory reports of Jewish ancestry.


  Species introductions to new habitats can cause a decline in the population
size of competing native species and consequently also in their genetic
diversity. We are interested in why these adverse effects are weak in some
cases whereas in others the native species declines to the point of extinction.
While the introduction rate and the growth rate of the introduced species in
the new environment clearly have a positive relationship with invasion success
and impact, the influence of competition is poorly understood. Here, we
investigate how the intensity of interspecific competition influences the
persistence time of a native species in the face of repeated and ongoing
introductions of the nonnative species. We analyze two stochastic models: a
model for the population dynamics of both species and a model that additionally
includes the population genetics of the native species at a locus involved in
its adaptation to a changing environment. Counterintuitively, both models
predict that the persistence time of the native species is lowest for an
intermediate intensity of competition. This phenomenon results from the
opposing effects of competition at different stages of the invasion process:
With increasing competition intensity more introduction events are needed until
a new species can establish, but increasing competition also speeds up the
exclusion of the native species by an established nonnative competitor. By
comparing the ecological and the eco-genetic model, we detect and quantify a
synergistic feedback between ecological and genetic effects.


  We propose a precise definition of a continuous time dynamical system made up
of interacting open subsystems. The interconnections of subsystems are coded by
directed graphs. We prove that the appropriate maps of graphs called graph
fibrations give rise to maps of dynamical systems. Consequently surjective
graph fibrations give rise to invariant subsystems and injective graph
fibrations give rise to projections of dynamical systems.


  Weights of edges and nodes on food webs which are available from the
empirical data hide much information about energy flows and biomass
distributions in ecosystem. We define a set of variables related to weights for
each species $i$, including the throughflow $T_i$, the total biomass $X_i$, and
the dissipated flow $D_i$ (output to the environment) to uncover the following
common patterns in 19 empirical weighted food webs: (1) DGBD distributions
(Discrete version of a Generalized Beta Distribution), a kind of deformed
Zipf's law, of energy flow and storage biomass; (2) The allometric scaling law
$T_i\propto X_i^{\alpha}$, which can be viewed as the counterpart of the
Kleiber's 3/4 law at the population level; (3) The dissipation law $D_i\propto
T_i^{\beta}$; and (4) The gravity law, including univariate version
$f_{ij}\propto (T_iT_j)^{\gamma}$ and bivariate approvement $f_{ij}\propto
T_i^{\gamma_1}T_j^{\gamma_2}$. These patterns are very common and significant
in all collected webs, as a result, some remarkable regularities are hidden in
weights.


  Residue-residue interactions that fold a protein into a unique
three-dimensional structure and make it play a specific function impose
structural and functional constraints on each residue site. Selective
constraints on residue sites are recorded in amino acid orders in homologous
sequences and also in the evolutionary trace of amino acid substitutions. A
challenge is to extract direct dependences between residue sites by removing
indirect dependences through other residues within a protein or even through
other molecules. Recent attempts of disentangling direct from indirect
dependences of amino acid types between residue positions in multiple sequence
alignments have revealed that the strength of inferred residue pair couplings
is an excellent predictor of residue-residue proximity in folded structures.
Here, we report an alternative attempt of inferring co-evolving site pairs from
concurrent and compensatory substitutions between sites in each branch of a
phylogenetic tree. First, branch lengths of a phylogenetic tree inferred by the
neighbor-joining method are optimized as well as other parameters by maximizing
a likelihood of the tree in a mechanistic codon substitution model. Mean
changes of quantities, which are characteristic of concurrent and compensatory
substitutions, accompanied by substitutions at each site in each branch of the
tree are estimated with the likelihood of each substitution. Partial
correlation coefficients of the characteristic changes along branches between
sites are calculated and used to rank co-evolving site pairs. Accuracy of
contact prediction based on the present co-evolution score is comparable to
that achieved by a maximum entropy model of protein sequences for 15 protein
families taken from the Pfam release 26.0. Besides, this excellent accuracy
indicates that compensatory substitutions are significant in protein evolution.


  We use a stochastic Markovian dynamics approach to describe the spreading of
vector-transmitted diseases, like dengue, and the threshold of the disease. The
coexistence space is composed by two structures representing the human and
mosquito populations. The human population follows a
susceptible-infected-recovered (SIR) type dynamics and the mosquito population
follows a susceptible-infected-susceptible (SIS) type dynamics. The human
infection is caused by infected mosquitoes and vice-versa so that the SIS and
SIR dynamics are interconnected. We develop a truncation scheme to solve the
evolution equations from which we get the threshold of the disease and the
reproductive ratio. The threshold of the disease is also obtained by performing
numerical simulations. We found that for certain values of the infection rates
the spreading of the disease is impossible whatever is the death rate of
infected mosquito.


  Despite the significant progress made in recent years, the computation of the
complete set of elementary flux modes of large or even genome-scale metabolic
networks is still impossible. We introduce a novel approach to speed up the
calculation of elementary flux modes by including transcriptional regulatory
information into the analysis of metabolic network. Taking into account gene
regulation dramatically reduces the solution space and allows the presented
algorithm to constantly eliminate biologically infeasible modes at an early
stage of the computation procedure. Thereby, the computational costs, such as
runtime, memory usage and disk space are considerably reduced. Consequently,
using the presented mode elimination algorithm pushes the size of metabolic
networks that can be studied by elementary flux modes to new limits.


  Motivation: Annotations are a key feature of many biological databases, used
to convey our knowledge of a sequence to the reader. Ideally, annotations are
curated manually, however manual curation is costly, time consuming and
requires expert knowledge and training. Given these issues and the exponential
increase of data, many databases implement automated annotation pipelines in an
attempt to avoid un-annotated entries. Both manual and automated annotations
vary in quality between databases and annotators, making assessment of
annotation reliability problematic for users. The community lacks a generic
measure for determining annotation quality and correctness, which we look at
addressing within this article. Specifically we investigate word reuse within
bulk textual annotations and relate this to Zipf's Principle of Least Effort.
We use UniProt Knowledge Base (UniProtKB) as a case study to demonstrate this
approach since it allows us to compare annotation change, both over time and
between automated and manually curated annotations.
  Results: By applying power-law distributions to word reuse in annotation, we
show clear trends in UniProtKB over time, which are consistent with existing
studies of quality on free text English. Further, we show a clear distinction
between manual and automated analysis and investigate cohorts of protein
records as they mature. These results suggest that this approach holds distinct
promise as a mechanism for judging annotation quality.
  Availability: Source code is available at the authors website:
http://homepages.cs.ncl.ac.uk/m.j.bell1/annotation.
  Contact: phillip.lord@newcastle.ac.uk


  Hydrogen bonds are a common feature in protein folding and aggregation. Due
to their chemical peculiarities in terms of strength and directionality, a
particular attention must be paid to the definition of the hydrogen bond
potential itself. This global target has been tackled through a computational
approach based on a minimalist description of the protein and the proper design
of algorithms, mainly using Monte Carlo and Kinetic Monte Carlo methods. We
have designed a hydrogen bond potential, see J. Chem. Phys. 132, 235102 (2010).
We have been performed a complete study of sequenceless peptide systems under
different conditions, such as temperature and concentration. To carry out full
protein studies, we need additional potentials to describe tertiary
interactions. We have discussed two different points of view. The first one is
the combination of the hydrogen bond potential with a structure-based one. We
have evaluated the implications of a proper definition of hydrogen bonds in the
thermodynamic and dynamic aspects of protein folding. See Biophys. J. 101,
1474-1482 (2011). We have undertaken a second strategy, combining a generic
hydrophobic model with the hydrogen bond one. The two main factors of folding
and aggregation are merged, then, to create a simple but complete potential.
Thanks to it, we have re-studied peptide aggregation including sequence.
Besides, we have simulated complete proteins with different folded shapes. We
have analyzed the competition between folding and aggregation, and how sequence
and hydrogen bonds can influence the interplay between them. See J. Chem. Phys.
136, 215103 (2012).


  Comparisons of DNA sequences between Neandertals and present-day humans have
shown that Neandertals share more genetic variants with non-Africans than with
Africans. This could be due to interbreeding between Neandertals and modern
humans when the two groups met subsequent to the emergence of modern humans
outside Africa. However, it could also be due to population structure that
antedates the origin of Neandertal ancestors in Africa. We measure the extent
of linkage disequilibrium (LD) in the genomes of present-day Europeans and find
that the last gene flow from Neandertals (or their relatives) into Europeans
likely occurred 37,000-86,000 years before the present (BP), and most likely
47,000-65,000 years ago. This supports the recent interbreeding hypothesis, and
suggests that interbreeding may have occurred when modern humans carrying Upper
Paleolithic technologies encountered Neandertals as they expanded out of
Africa.


  Cellular biology exists embedded in a world dominated by random dynamics and
chance. Many vital molecules and pieces of cellular machinery diffuse within
cells, moving along random trajectories as they collide with the other
biomolecular inhabitants of the cell. Cellular components may block each
other's progress, be produced or degraded at random times, and become unevenly
separated as cells grow and divide. Cellular behaviour, including important
features of stem cells, tumours and infectious bacteria, is profoundly
influenced by the chaos which is the environment within the cell walls. Here we
will look at some important causes and effects of randomness in cellular
biology, and some ways in which researchers, helped by the vast amounts of data
that are now flowing in, have made progress in describing the randomness of
nature.


  We introduce and analyze a purely competitive dynamics for the evolution of
an infinite population subject to a 3-strategy game. We argue that this
dynamics represents a characterization of how certain systems, both natural and
artificial, are governed. In each period, the population is randomly sorted
into pairs, which engage in a once-off play of the game; the probability that a
member propagates its type to its offspring is proportional only to its payoff
within the pair. We show that if a type is dominant (obtains higher payoffs in
games with both other types), its 'pure' population state, comprising only
members of that type, is globally attracting. If there is no dominant type,
there is an unstable 'mixed' fixed point; the population state eventually
oscillates between the three near-pure states. We then allow for mutations,
where offspring have a non-zero probability of randomly changing their type. In
this case, the existence of a dominant type renders a point near its pure state
globally attracting. If no dominant type exists, a supercritical Hopf
bifurcation occurs at the unique mixed fixed point, and above a critical
(typically low) mutation rate, this fixed point becomes globally attracting:
the implication is that even very low mutation rates can stabilize a system
that would, in the absence of mutations, be unstable.


  Since the foundations of Population Genetics the notion of genetic
equilibrium (in close analogy to Classical Mechanics) has been associated to
the Hardy-Weinberg (HW) Principle and the identification of equilibrium is
currently assumed by stating that the HW axioms are valid if appropriate values
of Chi-Square (p<0.05) are observed in experiments. Here we show by numerical
experiments with the genetic system of one locus/two alleles that considering
large ensembles of populations the Chi-Square test is not decisive and may lead
to false negatives in random mating populations and false positives in
nonrandom mating populations. As a result we confirm the logical statement that
statistical tests can not be used to deduce if the genetic population is under
the HW conditions. Furthermore, we show that under the HW conditions
populations of any finite size evolve in time according to what can be
identified as neutral dynamics to which the very notion of equilibrium is
unattainable for any practical purpose. Therefore, under the HW conditions
equilibrium properties are not observable. We also show that by relaxing the
condition of random mating the dynamics acquires all the characteristics of
asymptotic stable equilibrium. As a consequence our results show that the
question of equilibrium in genetic systems should be approached in close
analogy to non-equilibrium statistical physics and its observability should be
focused on dynamical quantities like the typical decay properties of the
allelic auto correlation function in time. In this perspective one should
abandon the classical notion of genetic equilibrium and its relation to the HW
proportions and open investigations in the direction of searching for unifying
general principles of population genetic transformations capable to take in
consideration these systems in their full complexity.


  Intrinsically disordered proteins (IDPs) do not possess well-defined
three-dimensional structures in solution under physiological conditions. We
develop all-atom, united-atom, and coarse-grained Langevin dynamics simulations
for the IDP alpha-synuclein that include geometric, attractive hydrophobic, and
screened electrostatic interactions and are calibrated to the inter-residue
separations measured in recent smFRET experiments. We find that alpha-synuclein
is disordered with conformational statistics that are intermediate between
random walk and collapsed globule behavior. An advantage of calibrated
molecular simulations over constraint methods is that physical forces act on
all residues, not only on residue pairs that are monitored experimentally, and
these simulations can be used to study oligomerization and aggregation of
multiple alpha-synuclein proteins that may precede amyloid formation.


  Zero Determinant (ZD) strategies are a new class of probabilistic and
conditional strategies that are able to unilaterally set the expected payoff of
an opponent in iterated plays of the Prisoner's Dilemma irrespective of the
opponent's strategy, or else to set the ratio between a ZD player's and their
opponent's expected payoff. Here we show that while ZD strategies are weakly
dominant, they are not evolutionarily stable and will instead evolve into less
coercive strategies. We show that ZD strategies with an informational advantage
over other players that allows them to recognize other ZD strategies can be
evolutionarily stable (and able to exploit other players). However, such an
advantage is bound to be short-lived as opposing strategies evolve to
counteract the recognition.


  Consistently predicting biopolymer structure at atomic resolution from
sequence alone remains a difficult problem, even for small sub-segments of
large proteins. Such loop prediction challenges, which arise frequently in
comparative modeling and protein design, can become intractable as loop lengths
exceed 10 residues and if surrounding side-chain conformations are erased. This
article introduces a modeling strategy based on a 'stepwise ansatz', recently
developed for RNA modeling, which posits that any realistic all-atom molecular
conformation can be built up by residue-by-residue stepwise enumeration. When
harnessed to a dynamic-programming-like recursion in the Rosetta framework, the
resulting stepwise assembly (SWA) protocol enables enumerative sampling of a 12
residue loop at a significant but achievable cost of thousands of CPU-hours. In
a previously established benchmark, SWA recovers crystallographic conformations
with sub-Angstrom accuracy for 19 of 20 loops, compared to 14 of 20 by KIC
modeling with a comparable expenditure of computational power. Furthermore, SWA
gives high accuracy results on an additional set of 15 loops highlighted in the
biological literature for their irregularity or unusual length. Successes
include cis-Pro touch turns, loops that pass through tunnels of other
side-chains, and loops of lengths up to 24 residues. Remaining problem cases
are traced to inaccuracies in the Rosetta all-atom energy function. In five
additional blind tests, SWA achieves sub-Angstrom accuracy models, including
the first such success in a protein/RNA binding interface, the YbxF/kink-turn
interaction in the fourth RNA-puzzle competition. These results establish
all-atom enumeration as a systematic approach to protein structure that can
leverage high performance computing and physically realistic energy functions
to more consistently achieve atomic resolution.


  The quantum network model with real variables is usually used to describe the
excitation energy transfer (EET) in the Fenna-Matthews-Olson(FMO) complexes. In
this paper we add the quantum phase factors to the hopping terms and find that
the quantum phase factors play an important role in the EET. The quantum phase
factors allow us to consider the space structure of the pigments. It is found
that phase coherence within the complexes would allow quantum interference to
affect the dynamics of the EET. There exist some optimal phase regions where
the transfer efficiency takes its maxima, which indicates that when the
pigments are optimally spaced, the exciton can pass through the FMO with
perfect efficiency. Moreover, the optimal phase regions almost do not change
with the environments. In addition, we find that the phase factors are useful
in the EET just in the case of multiple-pathway. Therefore, we demonstrate
that, the quantum phases may bring the other two factors, the optimal space of
the pigments and multiple-pathway, together to contribute the EET in
photosynthetic complexes with perfect efficiency.


  Information has an entropic character which can be analyzed within the
Statistical Theory in molecular systems. R. Landauer and C.H. Bennett showed
that a logical copy can be carried out in the limit of no dissipation if the
computation is performed sufficiently slowly. Structural and recent
single-molecule assays have provided dynamic details of polymerase machinery
with insight into information processing. We introduce a rigorous
characterization of Shannon Information in biomolecular systems and apply it to
DNA replication in the limit of no dissipation. Specifically, we devise an
equilibrium pathway in DNA replication to determine the entropy generated in
copying the information from a DNA template in the absence of friction. Both
the initial state, the free nucleotides randomly distributed in certain
concentrations, and the final state, a polymerized strand, are mesoscopic
equilibrium states for the nucleotide distribution. We use empirical stacking
free energies to calculate the probabilities of incorporation of the
nucleotides. The copied strand is, to first order of approximation, a state of
independent and non-indentically distributed random variables for which the
nucleotide that is incorporated by the polymerase at each step is dictated by
the template strand, and to second order of approximation, a state of
non-uniformly distributed random variables with nearest-neighbor interactions
for which the recognition of secondary structure by the polymerase in the
resultant double-stranded polymer determines the entropy of the replicated
strand. Two incorporation mechanisms arise naturally and their biological
meanings are explained. It is known that replication occurs far from
equilibrium and therefore the Shannon entropy here derived represents an upper
bound for replication to take place. Likewise, this entropy sets a universal
lower bound for the copying fidelity in replication.


  For many traits, including susceptibility to common diseases in humans,
causal loci uncovered by genetic mapping studies explain only a minority of the
heritable contribution to trait variation. Multiple explanations for this
"missing heritability" have been proposed. Here we use a large cross between
two yeast strains to accurately estimate different sources of heritable
variation for 46 quantitative traits and to detect underlying loci with high
statistical power. We find that the detected loci explain nearly the entire
additive contribution to heritable variation for the traits studied. We also
show that the contribution to heritability of gene-gene interactions varies
among traits, from near zero to 50%. Detected two-locus interactions explain
only a minority of this contribution. These results substantially advance our
understanding of the missing heritability problem and have important
implications for future studies of complex and quantitative traits.


  The quaternary structure of Filamin A (FLNa) 16-23 was recently shown to
exhibit multiple domain-domain interactions that lead to a propeller-like
construction. Here we present single molecule force spectroscopy experiments to
show a wide variety of mechanical responses of this molecule and compare it
with its linear counterpart FLNa 1-8. The compact structure of FLNa 16-23 leads
to a broad distribution of rupture forces and end-to-end lengths in the
force-extension mode and multiple unraveling timescales in the force-clamp
mode. Moreover, a subset of force-extension trajectories reveals a mechanical
hierarchy in which the rupture of domain-domain interactions at high forces
(200 pN) liberates the unfolding of individual domains at low forces (100 pN).
This mechanism may also explain the order of magnitude difference in the rates
of the biexponential fits to the distribution of unfolding dwell times under
force-clamp. Overall, FLNa 16-23 under a force of 100 pN is more compliant than
the linear FLNa 1-8. Since a physiological role of FLNa is to crosslink actin
filaments, this range of responses allows it to accommodate a broad spectrum of
forces exerted by the cell and its environment.


  We investigate the length distribution of self-assembled, long and stiff
polymers at thermal equilibrium. Our analysis is based on calculating the
partition functions of stiff polymers of variable lengths in the elastic
regime. Our conclusion is that the length distribution of this self-assembled
system follows closely the exponential distribution, except at the short length
limit. We then discuss the implications of our results on the experimentally
observed length distributions in amyloid fibrils.


  In this work we develop a stochastic model of acute HIV infection, based on
the well-known standard model, that allows us to simulate the complex mutation
pathways of HIV escape from multiple CTL responses. Under this model, we
describe two computational inference methods. In one, we use a Bayesian
approach to construct posteriors for the parameters of our model. In the
second, we use hypothesis testing to determine the fit of the model to data.
The methods are applied to two CHAVI datasets, demonstrating the importance of
accounting for the interaction of multiple mutant variants and
multi-directional selection in analysing HIV dynamics under CTL response.


  This paper proposes an original theory of aging of multicellular organisms.
The cells of multicellular organisms, in contrast to unicellular organisms, are
burdened with a two- part genome: housekeeping and specialized (multicellular),
responsible for ontogenesis and terminal differentiation. The two parts of the
genome compete for limited adaptive resources thereby interfering with the
ability of the house-keeping part of the genome to adequately perform
reparative and adaptive functions in post mitotic cells. The necessity to
complete the ontgenesis program, leads to increased activity of the
multicellular components of the genome. As a result, the allocation of cellular
resources to specialized genome con-tinuously increases with time. This leads
to a deficit of reparative and adaptive capacity in post mitotic cells.
Suggestions for future research focus on identifying groups of genes
responsible for regulation of growth rate of specialized genome and suppressing
ability of the cell division. A better understanding of the relationship
between the two parts of the genome will not only help us to manipulate
ontogenesis and aging, but will also improve our understanding of cancer
development and ontogenesis.


  Modeling molecules as undirected graphs and chemical reactions as graph
rewriting operations is a natural and convenient approach tom odeling
chemistry. Graph grammar rules are most naturally employed to model elementary
reactions like merging, splitting, and isomerisation of molecules. It is often
convenient, in particular in the analysis of larger systems, to summarize
several subsequent reactions into a single composite chemical reaction. We use
a generic approach for composing graph grammar rules to define a chemically
useful rule compositions. We iteratively apply these rule compositions to
elementary transformations in order to automatically infer complex
transformation patterns. This is useful for instance to understand the net
effect of complex catalytic cycles such as the Formose reaction. The
automatically inferred graph grammar rule is a generic representative that also
covers the overall reaction pattern of the Formose cycle, namely two carbonyl
groups that can react with a bound glycolaldehyde to a second glycolaldehyde.
Rule composition also can be used to study polymerization reactions as well as
more complicated iterative reaction schemes. Terpenes and the polyketides, for
instance, form two naturally occurring classes of compounds of utmost
pharmaceutical interest that can be understood as "generalized polymers"
consisting of five-carbon (isoprene) and two-carbon units, respectively.


  We study the population genetics of Evolution in the important special case
of weak selection, in which all fitness values are assumed to be close to one
another. We show that in this regime natural selection is tantamount to the
multiplicative updates game dynamics in a coordination game between genes.
Importantly, the utility maximized in this game, as well as the amount by which
each allele is boosted, is precisely the allele's mixability, or average
fitness, a quantity recently proposed in [1] as a novel concept that is crucial
in understanding natural selection under sex, thus providing a rigorous
demonstration of that insight. We also prove that the equilibria in two-person
coordination games can have large supports, and thus genetic diversity does not
suffer much at equilibrium. Establishing large supports involves answering
through a novel technique the following question: what is the probability that
for a random square matrix A both systems Ax = 1 and A^T y = 1 have positive
solutions? Both the question and the technique may be of broader interest.
  [1] A. Livnat, C. Papadimitriou, J. Dushoff, and M.W. Feldman. A mixability
theory for the role of sex in evolution. Proceedings of the National Academy of
Sciences, 105(50):19803-19808, 2008.


  The genetic diversity of a species is shaped by its recent evolutionary
history and can be used to infer demographic events or selective sweeps. Most
inference methods are based on the null hypothesis that natural selection is a
weak or infrequent evolutionary force. However, many species, particularly
pathogens, are under continuous pressure to adapt in response to changing
environments. A statistical framework for inference from diversity data of such
populations is currently lacking. Toward this goal, we explore the properties
of genealogies in a model of continual adaptation in asexual populations. We
show that lineages trace back to a small pool of highly fit ancestors, in which
almost simultaneous coalescence of more than two lineages frequently occurs.
While such multiple mergers are unlikely under the neutral coalescent, they
create a unique genetic footprint in adapting populations. The site frequency
spectrum of derived neutral alleles, for example, is non-monotonic and has a
peak at high frequencies, whereas Tajima's D becomes more and more negative
with increasing sample size. Since multiple merger coalescents emerge in many
models of rapid adaptation, we argue that they should be considered as a
null-model for adapting populations.


  An early burst of speciation followed by a subsequent slowdown in the rate of
diversification is commonly inferred from molecular phylogenies. This pattern
is consistent with some verbal theory of ecological opportunity and adaptive
radiations. One often-overlooked source of bias in these studies is that of
sampling at the level of whole clades, as researchers tend to choose large,
speciose clades to study. In this paper, we investigate the performance of
common methods across the distribution of clade sizes that can be generated by
a constant-rate birth-death process. Clades which are larger than expected for
a given constant-rate branching process tend to show a pattern of an early
burst even when both speciation and extinction rates are constant through time.
All methods evaluated were susceptible to detecting this false signature when
extinction was low. Under moderate extinction, both the gamma-statistic and
diversity-dependent models did not detect such a slowdown but only because the
signature of a slowdown was masked by subsequent extinction. Some models which
estimate time-varying speciation rates are able to detect early bursts under
higher extinction rates, but are extremely prone to sampling bias. We suggest
that examining clades in isolation may result in spurious inferences that rates
of diversification have changed through time.


  Positive selection distorts the structure of genealogies and hence alters
patterns of genetic variation within a population. Most analyses of these
distortions focus on the signatures of hitchhiking due to hard or soft
selective sweeps at a single genetic locus. However, in linked regions of
rapidly adapting genomes, multiple beneficial mutations at different loci can
segregate simultaneously within the population, an effect known as clonal
interference. This leads to a subtle interplay between hitchhiking and
interference effects, which leads to a unique signature of rapid adaptation on
genetic variation both at the selected sites and at linked neutral loci. Here,
we introduce an effective coalescent theory (a "fitness-class coalescent") that
describes how positive selection at many perfectly linked sites alters the
structure of genealogies. We use this theory to calculate several simple
statistics describing genetic variation within a rapidly adapting population,
and to implement efficient backwards-time coalescent simulations which can be
used to predict how clonal interference alters the expected patterns of
molecular evolution.


  Motivation: The comparison of diverse genomic datasets is fundamental to
understanding genome biology. Researchers must explore many large datasets of
genome intervals (e.g., genes, sequence alignments) to place their experimental
results in a broader context and to make new discoveries. Relationships between
genomic datasets are typically measured by identifying intervals that
intersect: that is, they overlap and thus share a common genome interval. Given
the continued advances in DNA sequencing technologies, efficient methods for
measuring statistically significant relationships between many sets of genomic
features is crucial for future discovery.
  Results: We introduce the Binary Interval Search (BITS) algorithm, a novel
and scalable approach to interval set intersection. We demonstrate that BITS
outperforms existing methods at counting interval intersections. Moreover, we
show that BITS is intrinsically suited to parallel computing architectures such
as Graphics Processing Units (GPUs) by illustrating its utility for efficient
Monte-Carlo simulations measuring the significance of relationships between
sets of genomic intervals.


  Background: Recently developed techniques to study the spread of infectious
diseases through networks make assumptions that the initial proportion infected
is infinitesimal and the population behavior is static throughout the epidemic.
The models do not apply if the initial proportion is large (and fail whenever
R_0<1), and cannot measure the impact of an intervention.
  Methods: In this paper we adapt "edge-based compartmental models" to
situations having finite-sized initial conditions.
  Results: The resulting models remain simple and accurately capture the effect
of the initial conditions. It is possible to generalize the model to networks
whose partnerships change in time.
  Conclusions: The resulting models can be applied to a range of important
contexts. The models can be used to choose between different interventions that
affect the disease or the population structure.


  Our wellbeing depends as much on our personal success, as it does on the
success of our society. The realization of this fact makes cooperation a very
much needed trait. Experiments have shown that rewards can elevate our
readiness to cooperate, but since giving a reward inevitably entails paying a
cost for it, the emergence and stability of such behavior remain elusive. Here
we show that allowing for the act of rewarding to self-organize in dependence
on the success of cooperation creates several evolutionary advantages that
instill new ways through which collaborative efforts are promoted. Ranging from
indirect territorial battle to the spontaneous emergence and destruction of
coexistence, phase diagrams and the underlying spatial patterns reveal
fascinatingly reach social dynamics that explains why this costly behavior has
evolved and persevered. Comparisons with adaptive punishment, however, uncover
an Achilles heel of adaptive rewarding that is due to over-aggression, which in
turn hinders optimal utilization of network reciprocity. This may explain why,
despite of its success, rewarding is not as firmly weaved into our societal
organization as punishment.


  Here we propose the approach to design of the new type of hybrids of
oligonucleotides with fluorescein-functionalized single-walled carbon
nanotubes. The approach is based on stacking interactions of functionalized
nanotubes with pyrene residues in conjugates of oligonucleotides. The amino-
and fluorescein-modified single-walled carbon nanotubes were obtained, and
their physico-chemical properties were investigated. The effect of carbon
nanotubes functionalization type on the efficacy of sorption of pyrene
conjugates of oligonucleotides was examined. Proposed non-covalent hybrids of
fluorescein-labeled carbon nanotubes with oligonucleotides may be used for
intracellular transport of functional nucleic acids.


  Polyploidy is an important speciation mechanism, particularly in land plants.
Allopolyploid species are formed after hybridization between otherwise
intersterile parental species. Recent theoretical progress has led to
successful implementation of species tree models that take population genetic
parameters into account. However, these models have not included allopolyploid
hybridization and the special problems imposed when species trees of
allopolyploids are inferred. Here, two new models for the statistical inference
of the evolutionary history of allopolyploids are evaluated using simulations
and demonstrated on two empirical data sets. It is assumed that there has been
a single hybridization event between two diploid species resulting in a genomic
allotetraploid. The evolutionary history can be represented as a network or as
a multiply labeled tree, in which some pairs of tips are labeled with the same
species. In one of the models (AlloppMUL), the multiply labeled tree is
inferred directly. This is the simplest model and the most widely applicable,
since fewer assumptions are made. The second model (AlloppNET) incorporates the
hybridization event explicitly which means that fewer parameters need to be
estimated. Both models are implemented in the BEAST framework. Simulations show
that both models are useful and that AlloppNET is more accurate if the
assumptions it is based on are valid. The models are demonstrated on previously
analyzed data from the genus Pachycladon (Brassicaceae) and from the genus
Silene (Caryophyllaceae).


  Identification and quantification of condition-specific transcripts using
RNA-Seq is vital in transcriptomics research. While initial efforts using
mathematical or statistical modeling of read counts or per-base exonic signal
have been successful, they may suffer from model overfitting since not all the
reference transcripts in a database are expressed under a specific biological
condition. Standard shrinkage approaches, such as Lasso, shrink all the
transcript abundances to zero in a non-discriminative manner. Thus it does not
necessarily yield the set of condition-specific transcripts. Informed shrinkage
approaches, using the observed exonic coverage signal, are thus desirable.
Motivated by ubiquitous uncovered exonic regions in RNA-Seq data, termed as
"naked exons", we propose a new computational approach that first filters out
the reference transcripts not supported by splicing and paired-end reads, then
followed by fitting a new mathematical model of per-base exonic coverage signal
and the underlying transcripts structure. We introduce a tuning parameter to
penalize the specific regions of the selected transcripts that were not
supported by the naked exons. Our approach compares favorably with the selected
competing methods in terms of both time complexity and accuracy using simulated
and real-world data. Our method is implemented in SAMMate, a GUI software suite
freely available from http://sammate.sourceforge.net


  MicroRNAs (miRNAs) are an emerging class of biomarkers that are frequently
deregulated in cancer cells and have shown a great promise for cancer
classification and prognosis. In this work, we developed a three-way
electrochemical sensor for detection and quantitation of ultra-low amounts of
miRNAs in a wide dynamic range of measured concentrations. The sensor
facilitates three detection modes based on hybridization (H-SENS), p19 protein
binding (P-SENS) and protein displacement (D-SENS). The combined HPD sensor
(HPD-SENS) identifies as low as 5 aM or 90 molecules of miRNA per 30\muL of
sample without PCR amplification, and can be operated within the dynamic range
from 10 aM to 1 \muM. The HPD senor is made on a commercially-available gold
nanoparticles-modified electrode and suitable for analysing multiple miRNAs on
a single electrode. This three-way sensor exhibits high selectivity and
specificity and was successfully used for sequential analysis of miR-32 and
miR-122 on one electrode. Furthermore, the P-SENS was employed for profiling of
three endogenous microRNAs (miR-21, miR-32 and miR-122) in human blood serum.


  Force-clamp spectroscopy reveals the unfolding and disulfide bond rupture
times of single protein molecules as a function of the stretching force, point
mutations and solvent conditions. The statistics of these times reveal whether
the protein domains are independent of one another, the mechanical hierarchy in
the polyprotein chain, and the functional form of the probability distribution
from which they originate. It is therefore important to use robust statistical
tests to decipher the correct theoretical model underlying the process. Here we
develop multiple techniques to compare the well-established experimental data
set on ubiquitin with existing theoretical models as a case study. We show that
robustness against filtering, agreement with a maximum likelihood function that
takes into account experimental artifacts, the Kuiper statistic test and
alignment with synthetic data all identify the Weibull or stretched exponential
distribution as the best fitting model. Our results are inconsistent with
recently proposed models of Gaussian disorder in the energy landscape or noise
in the applied force as explanations for the observed non-exponential kinetics.
Since the physical model in the fit affects the characteristic unfolding time,
these results have important implications on our understanding of the
biological function of proteins.


  This paper presents a novel framework for the modeling of biological
networks. It makes use of recent tools analyzing the robust satisfaction of
properties of (hybrid) dynamical systems. The main challenge of this approach
as applied to biological systems is to get access to the relevant parameter
sets despite gaps in the available knowledge. An initial estimate of useful
parameters was sought by formalizing the known behavior of the biological
network in the STL logic using the tool Breach. Then, once a set of parameter
values consistent with known biological properties was found, we tried to
locally expand it into the largest possible valid region. We applied this
methodology in an effort to model and better understand the complex network
regulating iron homeostasis in mammalian cells. This system plays an important
role in many biological functions, including erythropoiesis, resistance against
infections, and proliferation of cancer cells.


  The evolutionary significance of hybridization and subsequent introgression
has long been appreciated, but evaluation of the genome-wide effects of these
phenomena has only recently become possible. Crop-wild study systems represent
ideal opportunities to examine evolution through hybridization. For example,
maize and the conspecific wild teosinte Zea mays ssp. mexicana, (hereafter,
mexicana) are known to hybridize in the fields of highland Mexico. Despite
widespread evidence of gene flow, maize and mexicana maintain distinct
morphologies and have done so in sympatry for thousands of years. Neither the
genomic extent nor the evolutionary importance of introgression between these
taxa is understood. In this study we assessed patterns of genome-wide
introgression based on 39,029 single nucleotide polymorphisms genotyped in 189
individuals from nine sympatric maize-mexicana populations and reference
allopatric populations. While portions of the maize and mexicana genomes were
particularly resistant to introgression (notably near known
cross-incompatibility and domestication loci), we detected widespread evidence
for introgression in both directions of gene flow. Through further
characterization of these regions and preliminary growth chamber experiments,
we found evidence suggestive of the incorporation of adaptive mexicana alleles
into maize during its expansion to the highlands of central Mexico. In
contrast, very little evidence was found for adaptive introgression from maize
to mexicana. The methods we have applied here can be replicated widely, and
such analyses have the potential to greatly informing our understanding of
evolution through introgressive hybridization. Crop species, due to their
exceptional genomic resources and frequent histories of spread into sympatry
with relatives, should be particularly influential in these studies.


  A mesoscopic model which allows us to identify and quantify the strength of
binding sites in DNA sequences is proposed. The model is based on the
Peyrard-Bishop-Dauxois model for the DNA chain coupled to a Brownian particle
which explores the sequence interacting more importantly with open base pairs
of the DNA chain. We apply the model to promoter sequences of different
organisms. The free energy landscape obtained for these promoters shows a
complex structure that is strongly connected to their biological behavior. The
analysis method used is able to quantify free energy differences of sites
within genome sequences.


  We perform an analysis of a recent spatial version of the classical
Lotka-Volterra model, where a finite scale controls individuals' interaction.
We study the behavior of the predator-prey dynamics in physical spaces higher
than one, showing how spatial patterns can emerge for some values of the
interaction range and of the diffusion parameter.


  Affymetrix GeneChip technology is used to detect gene expression levels in
samples of cells under different conditions. In this project, we analyzed the
gene expression profiling data for mouse induced pluripotent stem cell (iPSCs)
(Takahashi, 2006) on Affymetrix Mouse 430 2.0 GeneChip. Three biological
conditions were present: p53KO, microRNA mir34aKO, and wild type, each with
three biological replicates. The first part was devoted to identifying
differentially expressed genes from around 45,000 of them, and looking into
their biological meanings by pathway analysis. The second part dealt with
repetitive elements represented in the pool of mRNAs. We identified repetitive
elements that show a significant difference between two biological conditions.
Both the comparison of p53KO versus WT and mir34aKO versus WT were done.
However, the emphasis was on the former. Laboratory validation with qPCR
confirmed our findings. This work was done under the Overseas Research
Fellowship (ORF) Scheme 2012 for Science Students by the Faculty of Science,
The University of Hong Kong. Many thanks are due to the University for the
fellowship, and to Professors Terry Speed and Lin He and Drs Chao-po Lin and
Anne Biton of the University of California at Berkeley for their supervision
and generous support.


  Constructing gene interaction networks (GINs) from high-throughput gene
expression data is an important and challenging problem in systems biology.
Existing algorithms produce networks that either have undirected and unweighted
edges, or else are constrained to contain no cycles, both of which are
biologically unrealistic. In the present paper we propose a new algorithm,
based on a concept from probability theory known as the phi-mixing coefficient,
that produces networks whose edges are weighted and directed, and are permitted
to contain cycles. Because there is no "ground truth" for genome-wide networks
on a human scale, we analyzed the outcomes of several experiments on lung
cancer, and matched the predictions from the inferred networks with
experimental results. Specifically, we inferred three networks (NSCLC,
Neuro-endocrine NSCLC plus SCLC, and normal) from the gene expression
measurements of 157 lung cancer and 59 normal cell lines, compared with the
outcomes of siRNA screening of 19,000+ genes on 11 NSCLC cell lines, and
analyzed data from a ChIP-Seq experiment to determine putative downstream
targets of the lineage specific oncogenic transcription factor ASCL1. The
inferred networks displayed a scale-free or power law behavior between the
degree of a node and the number of nodes with that degree. There was a strong
correlation between the degree of a gene in the inferred NSCLC network and its
essentiality for the survival of the cells. The inferred downstream
neighborhood genes of ASCL1 in the SCLC network were significantly enriched by
ChIP-Seq determined putative target genes, while no such enrichment was found
in the inferred NSCLC network.


  Whether or not to change strategy depends not only on the personal success of
each individual, but also on the success of others. Using this as motivation,
we study the evolution of cooperation in games that describe social dilemmas,
where the propensity to adopt a different strategy depends both on individual
fitness as well as on the strategies of neighbors. Regardless of whether the
evolutionary process is governed by pairwise or group interactions, we show
that plugging into the "wisdom of groups" strongly promotes cooperative
behavior. The more the wider knowledge is taken into account the more the
evolution of defectors is impaired. We explain this by revealing a dynamically
decelerated invasion process, by means of which interfaces separating different
domains remain smooth and defectors therefore become unable to efficiently
invade cooperators. This in turn invigorates spatial reciprocity and
establishes decentralized decision making as very beneficial for resolving
social dilemmas.


  Gene annotation databases (compendiums maintained by the scientific community
that describe the biological functions performed by individual genes) are
commonly used to evaluate the functional properties of experimentally derived
gene sets. Overlap statistics, such as Fisher's Exact Test (FET), are often
employed to assess these associations, but don't account for non-uniformity in
the number of genes annotated to individual functions or the number of
functions associated with individual genes. We find FET is strongly biased
toward over-estimating overlap significance if a gene set has an unusually high
number of annotations. To correct for these biases, we develop Annotation
Enrichment Analysis (AEA), which properly accounts for the non-uniformity of
annotations. We show that AEA is able to identify biologically meaningful
functional enrichments that are obscured by numerous false-positive enrichment
scores in FET, and we therefore suggest it be used to more accurately assess
the biological properties of gene sets.


  In this short note we review deterministic simulation of biochemical
pathways, i.e. networks of biochemical reactions obeying the law of mass
action. It is meant as a basis for the MATLAB code, written by the author,
which permits easy input and simulation of general biochemical networks. This
work was carried out for the European Project `CardioWorkBench'.


  Many common diseases are highly polygenic, modulated by a large number
genetic factors with small effects on susceptibility to disease. These small
effects are difficult to map reliably in genetic association studies. To
address this problem, researchers have developed methods that aggregate
information over sets of related genes, such as biological pathways, to
identify gene sets that are enriched for genetic variants associated with
disease. However, these methods fail to answer a key question: which genes and
genetic variants are associated with disease risk? We develop a method based on
sparse multiple regression that simultaneously identifies enriched pathways,
and prioritizes the variants within these pathways, to locate additional
variants associated with disease susceptibility. A central feature of our
approach is an estimate of the strength of enrichment, which yields a coherent
way to prioritize variants in enriched pathways. We illustrate the benefits of
our approach in a genome-wide association study of Crohn's disease with
~440,000 genetic variants genotyped for ~4700 study subjects. We obtain strong
support for enrichment of IL-12, IL-23 and other cytokine signaling pathways.
Furthermore, prioritizing variants in these enriched pathways yields support
for additional disease-association variants, all of which have been
independently reported in other case-control studies for Crohn's disease.


  In research labs, there is often a need to customise software at every step
in a given bioinformatics workflow, but traditionally it has been difficult to
obtain both a high degree of customisability and good performance.
Performance-sensitive tools are often highly monolithic, which can make
research difficult. We present a novel set of software development principles
and a bioinformatics framework, Friedrich, which is currently in early
development. Friedrich applications support both early stage experimentation
and late stage batch processing, since they simultaneously allow for good
performance and a high degree of flexibility and customisability. These
benefits are obtained in large part by basing Friedrich on the multiparadigm
programming language Scala. We present a case study in the form of a basic
genome assembler and its extension with new functionality. Our architecture has
the potential to greatly increase the overall productivity of software
developers and researchers in bioinformatics.


  The ABO histo-blood group, the critical determinant of transfusion
incompatibility, was the first genetic polymorphism discovered in humans.
Remarkably, ABO antigens are also polymorphic in many other primates, with the
same two amino acid changes responsible for A and B specificity in all species
sequenced to date. Whether this recurrence of A and B antigens is the result of
an ancient polymorphism maintained across species or due to numerous, more
recent instances of convergent evolution has been debated for decades, with a
current consensus in support of convergent evolution. We show instead that
genetic variation data in humans and gibbons as well as in Old World Monkeys
are inconsistent with a model of convergent evolution and support the
hypothesis of an ancient, multi-allelic polymorphism of which some alleles are
shared by descent among species. These results demonstrate that the ABO
polymorphism is a trans-species polymorphism among distantly related species
and has remained under balancing selection for tens of millions of years, to
date, the only such example in Hominoids and Old World Monkeys outside of the
Major Histocompatibility Complex.


  A major challenge in biology is explaining how novel characters originate,
however, the molecular mechanisms that underlie the emergence of evolutionary
innovations are unclear. Here we show that while gene expression in the uterus
evolves at a slow and relatively constant rate, it has been punctuated by
periods of rapid change associated with the recruitment of thousands of genes
into uterine expression during the evolution of pregnancy in mammals. We found
that numerous genes and signaling pathways essential for the establishment of
pregnancy and maternal-fetal communication evolved uterine expression in
mammals. Remarkably the majority of genes recruited into endometrial expression
have cis-regulatory elements derived from lineage-specific transposons,
suggesting that that bursts of transposition facilitate adaptation and
speciation through genomic and regulatory reorganization.


  We analyze the gene expression data of Zebrafish under the combined framework
of complex networks and random matrix theory. The nearest neighbor spacing
distribution of the corresponding matrix spectra follows random matrix
predictions of Gaussian orthogonal statistics. Based on the eigenvector
analysis we can divide the spectra into two parts, first part for which the
eigenvector localization properties match with the random matrix theory
predictions, and the second part for which they show deviation from the theory
and hence are useful to understand the system dependent properties. Spectra
with the localized eigenvectors can be characterized into three groups based on
the eigenvalues. We explore the position of localized nodes from these
different categories. Using an overlap measure, we find that the top
contributing nodes in the different groups carry distinguished structural
features. Furthermore, the top contributing nodes of the different localized
eigenvectors corresponding to the lower eigenvalue regime form different
densely connected structure well separated from each other. Preliminary
biological interpretation of the genes, associated with the top contributing
nodes in the localized eigenvectors, suggests that the genes corresponding to
same vector share common features.


  It is believed that the quantum coherence itself cannot explain the very high
excitation energy transfer (EET) efficiency in the Fenna-Matthews-Olson (FMO)
complex. In this paper, we show that this is not the case if the inter-site
couplings take complex values. By phenomenologically introducing phases into
the inter-site couplings, we obtain the EET efficiency as high as 0.8972 in
contrast to 0.6781 with real inter-site couplings. Dependence of the excitation
energy transfer efficiency on the initial states is elaborated. Effects of
fluctuations in the site energies and inter-site couplings are also examined.


  (ABRIDGED) We report the genome sequencing of 139 wild-derived strains of D.
melanogaster, representing 22 population samples from the sub-Saharan ancestral
range of this species, along with one European population. Most genomes were
sequenced above 25X depth from haploid embryos. Results indicated a pervasive
influence of non-African admixture in many African populations, motivating the
development and application of a novel admixture detection method. Admixture
proportions varied among populations, with greater admixture in urban
locations. Admixture levels also varied across the genome, with localized peaks
and valleys suggestive of a non-neutral introgression process. Genomes from the
same location differed starkly in ancestry, suggesting that isolation
mechanisms may exist within African populations. After removing putatively
admixed genomic segments, the greatest genetic diversity was observed in
southern Africa (e.g. Zambia), while diversity in other populations was largely
consistent with a geographic expansion from this potentially ancestral region.
The European population showed different levels of diversity reduction on each
chromosome arm, and some African populations displayed chromosome arm-specific
diversity reductions. Inversions in the European sample were associated with
strong elevations in diversity across chromosome arms. Genomic scans were
conducted to identify loci that may represent targets of positive selection. A
disproportionate number of candidate selective sweep regions were located near
genes with varied roles in gene regulation. Outliers for Europe-Africa FST were
found to be enriched in genomic regions of locally elevated cosmopolitan
admixture, possibly reflecting a role for some of these loci in driving the
introgression of non-African alleles into African populations.