  The innate immune system in insects consists of a conserved core signaling
network and rapidly diversifying effector and recognition components, often
containing a high proportion of taxonomically-restricted genes. In the absence
of functional annotation, genes encoding immune system proteins can thus be
difficult to identify, as homology-based approaches generally cannot detect
lineage-specific genes. Here, we use RNA-seq to compare the uninfected and
infection-induced transcriptome in the parasitoid wasp Nasonia vitripennis to
identify genes regulated by infection. We identify 183 genes significantly
up-regulated by infection and 61 genes significantly down-regulated by
infection. We also produce a new homology-based immune catalog in N.
vitripennis, and show that most infection-induced genes are not assigned an
immune function from homology alone, suggesting the potential for substantial
novel immune components in less-well-studied systems. Finally, we show that a
high proportion of these novel induced genes are taxonomically-restricted,
highlighting the rapid evolution of immune gene content. The combination of
functional annotation using RNA-seq and homology-based annotation provides a
robust method to characterize the innate immune response across a wide variety
of insects, and reveals significant novel features of the Nasonia immune
response.


  Gene regulatory networks (GRN) govern phenotypic adaptations and reflect the
trade-offs between physiological responses and evolutionary adaptation that act
at different time scales. To identify patterns of molecular function and
genetic diversity in GRNs, we studied the drought response of the common
sunflower, Helianthus annuus, and how the underlying GRN is related to its
evolution. We examined the responses of 32,423 expressed sequences to drought
and to abscisic acid and selected 145 co-expressed transcripts. We
characterized their regulatory relationships in nine kinetic studies based on
different hormones. From this, we inferred a GRN by meta-analyses of a Gaussian
graphical model and a random forest algorithm and studied the genetic
differentiation among populations (FST) at nodes. We identified two main hubs
in the network that transport nitrate in guard cells. This suggests that
nitrate transport is a critical aspect of sunflower physiological response to
drought. We observed that differentiation of the network genes in elite
sunflower cultivars is correlated with their position and connectivity. This
systems biology approach combined molecular data at different time scales and
identified important physiological processes. At the evolutionary level, we
propose that network topology could influence responses to human selection and
possibly adaptation to dry environments.


  Inference of individual admixture coefficients, which is important for
population genetic and association studies, is commonly performed using
compute-intensive likelihood algorithms. With the availability of large
population genomic data sets, fast versions of likelihood algorithms have
attracted considerable attention. Reducing the computational burden of
estimation algorithms remains, however, a major challenge. Here, we present a
fast and efficient method for estimating individual admixture coefficients
based on sparse non-negative matrix factorization algorithms. We implemented
our method in the computer program sNMF, and applied it to human and plant
genomic data sets. The performances of sNMF were then compared to the
likelihood algorithm implemented in the computer program ADMIXTURE. Without
loss of accuracy, sNMF computed estimates of admixture coefficients within
run-times approximately 10 to 30 times faster than those of ADMIXTURE.


  The possibility of complicated dynamic behaviour driven by non-linear
feedbacks in dynamical systems has revolutionized science in the latter part of
the last century. Yet despite examples of complicated frequency dynamics, the
possibility of long-term evolutionary chaos is rarely considered. The concept
of "survival of the fittest" is central to much evolutionary thinking and
embodies a perspective of evolution as a directional optimization process
exhibiting simple, predictable dynamics. This perspective is adequate for
simple scenarios, when frequency-independent selection acts on scalar
phenotypes. However, in most organisms many phenotypic properties combine in
complicated ways to determine ecological interactions, and hence
frequency-dependent selection. Therefore, it is natural to consider models for
the evolutionary dynamics generated by frequency-dependent selection acting
simultaneously on many different phenotypes. Here we show that complicated,
chaotic dynamics of long-term evolutionary trajectories in phenotype space is
very common in a large class of such models when the dimension of phenotype
space is large, and when there are epistatic interactions between the
phenotypic components. Our results suggest that the perspective of evolution as
a process with simple, predictable dynamics covers only a small fragment of
long-term evolution. Our analysis may also be the first systematic study of the
occurrence of chaos in multidimensional and generally dissipative systems as a
function of the dimensionality of phase space.


  Cell fate decisions in multicellular organisms are precisely coordinated,
leading to highly reproducible macroscopic outcomes of developmental processes.
The origins of this reproducibility can be found at the molecular level during
the earliest stages of development when spatial patterns of morphogen
(form-generating) molecules emerge reproducibly. However, the initial
conditions for these early stages are determined by the female during
oogenesis, and it is unknown whether reproducibility is passed on to the zygote
or whether it is reacquired by the zygote. Here we examine the earliest
reproducible pattern in the Drosophila embryo, the Bicoid protein gradient.
Using a unique combination of absolute molecule counting techniques, we show
that it is generated from a highly controlled source of mRNA molecules that is
reproducible from embryo to embryo to within ~8%. This occurs in a perfectly
linear feed-forward process: changes in the female's gene dosage lead to
proportional changes in the mRNA and protein counts in the embryo. In this
setup, noise is kept low in the transition from one molecular species to
another, allowing the female to precisely deposit the same absolute number of
mRNA molecules in each embryo and therefore confer reproducibility to the
Bicoid pattern. Our results indicate that the reproducibility of the
morphological structures that emerge in the embryo originates during oogenesis
when all initial patterning signals are controlled with precision similar to
what we observe for the Bicoid pattern.


  A quantum-based magnetic compass sensor, mediated through radical pair
reactions, has been suggested to underlie the sensory ability of migrating
birds to receive directional information from the geomagnetic field. Here we
extend the currently available models by considering the effects of slow
fluctuations in the nuclear spin environment on the directional signal. We
quantitatively evaluate the robustness of signals under fluctuations on a
timescale longer than the lifetime of a radical pair, utilizing two models of
radical pairs. Our results suggest design principles for building a
radical-pair based compass sensor that is both robust and highly directional
sensitive.


  The dimensions of neuronal synapses suggest that optical super-resolution
imaging methods are necessary for thorough investigation of protein
distributions and interactions. Nanoscopic evaluation of neuronal samples has
presented practical hurdles, but advancing methods are making synaptic protein
topology and quantification measurements feasible. This work explores the
application of Photoactivated Localization Microscopy (PALM) pointilistic
super-resolution imaging for investigation of the membrane bound sodium pump,
the Na,K-ATPase, in matured neurons.
  Two isoforms of the sodium pump (ATP1a1 and ATP1a3) were studied in cultured
neurons using the PALM-compatible fluorescent proteins PAGFP and mEos.
Nanoscopic imaging reveals a compartmentalized distribution of sodium pumps in
dendritic spines. Several nanoclusters of pumps are typically found in the
spine head and relatively few in the spine neck. The density of sodium pumps
was estimated from a quantification of detected single molecules to 600-800
pump copies/um2 in the spine heads.
  We have utilized PALM for dissection of nanoscale localization in mature
cultured neurons and demonstrated similar topology and quantification estimates
with PAGFP and mEos. PALM topology assessments of the sodium pump appeared
similar to previous STED studies, though quantification estimates varied
notably, implying that sample labeling strategies and choice of nanoscopic
imaging method are critical factors for correct quantification.


  The history of infections and epidemics holds famous examples where
understanding, containing and ultimately treating an outbreak began with
understanding its mode of spread. Influenza, HIV and most computer viruses,
spread person to person, device to device, through contact networks; Cholera,
Cancer, and seasonal allergies, on the other hand, do not. In this paper we
study two fundamental questions of detection: first, given a snapshot view of a
(perhaps vanishingly small) fraction of those infected, under what conditions
is an epidemic spreading via contact (e.g., Influenza), distinguishable from a
"random illness" operating independently of any contact network (e.g., seasonal
allergies); second, if we do have an epidemic, under what conditions is it
possible to determine which network of interactions is the main cause of the
spread -- the causative network -- without any knowledge of the epidemic, other
than the identity of a minuscule subsample of infected nodes?
  The core, therefore, of this paper, is to obtain an understanding of the
diagnostic power of network information. We derive sufficient conditions
networks must satisfy for these problems to be identifiable, and produce
efficient, highly scalable algorithms that solve these problems. We show that
the identifiability condition we give is fairly mild, and in particular, is
satisfied by two common graph topologies: the grid, and the Erdos-Renyi graphs.


  Human populations have experienced dramatic growth since the Neolithic
revolution. Recent studies that sequenced a very large number of individuals
observed an extreme excess of rare variants, and provided clear evidence of
recent rapid growth in effective population size, though estimates have varied
greatly among studies. All these studies were based on protein-coding genes, in
which variants are also impacted by natural selection. In this study, we
introduce targeted sequencing data for studying recent human history with
minimal confounding by natural selection. We sequenced loci very far from genes
that meet a wide array of additional criteria such that mutations in these loci
are putatively neutral. As population structure also skews allele frequencies,
we sequenced a sample of relatively homogeneous ancestry by first analyzing the
population structure of 9,716 European Americans. We employed very high
coverage sequencing to reliably call rare variants, and fit an extensive array
of models of recent European demographic history to the site frequency
spectrum. The best-fit model estimates ~3.4% growth per generation during the
last ~140 generations, resulting in a population size increase of two orders of
magnitude. This model fits the data very well, largely due to our observation
that assumptions of more ancient demography can impact estimates of recent
growth. This observation and results also shed light on the discrepancy in
demographic estimates among recent studies.


  What can relics of the past tell us about the thoughts and beliefs of the
people who invented and used them? Recent collaborations at the frontier of
archaeology, anthropology, and cognitive science are culminating in speculative
but nevertheless increasingly sophisticated efforts to unravel how modern human
cognition came about. By considering objects within their archaeological
context, we have begun to piece together something of the way of life of people
who inhabited particular locales, which in turn reflects their underlying
thought processes.


  Host mobility plays a fundamental role in the spatial spread of infectious
diseases. Previous theoretical works based on the integration of network theory
into the metapopulation framework have shown that the heterogeneities that
characterize real mobility networks favor the propagation of epidemics.
Nevertheless, the studies conducted so far assumed the mobility process to be
either Markovian or non-Markovian with a fixed traveling time scale. Available
statistics however show that the time spent by travelers at destination is
characterized by wide fluctuations, ranging between a single day up to several
months. Such varying length of stay crucially affects the chance and duration
of mixing events among hosts and may therefore have a strong impact on the
spread of an emerging disease. Here, we present an analytical and computational
study of epidemic processes on a complex subpopulation network where travelers
have memory of their origin and spend a heterogeneously distributed time
interval at their destination. Through analytical calculations and numerical
simulations we show that the heterogeneity of the length of stay alters the
expression of the threshold between local outbreak and global invasion, and,
moreover, it changes the epidemic behavior of the system in case of a global
outbreak. Additionally, our theoretical framework allows us to study the effect
of changes in the traveling behavior in response to the infection, by
considering a scenario in which sick individuals do not leave their home
location. Finally, we compare the results of our non-Markovian framework with
those obtained with a classic Markovian approach and find relevant differences
between the two, in the estimate of the epidemic invasion potential, as well as
of the timing and the pattern of its spatial spread.


  Recent experiments reveal both passive subdiffusion of various nanoparticles
and anomalous active transport of such particles by molecular motors in the
molecularly crowded environment of living biological cells. Passive and active
microrheology reveals that the origin of this anomalous dynamics is due to the
viscoelasticity of the intracellular fluid. How do molecular motors perform in
such a highly viscous, dissipative environment? Can we explain the observed
co-existence of the anomalous transport of relatively large particles of 100 to
500 nm in size by kinesin motors with the normal transport of smaller particles
by the same molecular motors? What is the efficiency of molecular motors in the
anomalous transport regime? Here we answer these seemingly conflicting
questions and consistently explain experimental findings in a generalization of
the well-known continuous diffusion model for molecular motors with two
conformational states in which viscoelastic effects are included.


  In this study, we forecast the population of the Philippines using a discrete
age-structured compartmental model. We estimated the future population
structure of the Philippines if the government imposes an n-child policy on top
of the declining current birth and death rate trend of the country.


  The chemical concrete machine is a graph rewriting system which uses only
local moves (rewrites), seen as chemical reactions involving molecules which
are graphs made up by 4 trivalent nodes. It is Turing complete, therefore it
might be used as a model of computation in algorithmic chemistry.


  We describe the Microbial Community Reconstruction ({\bf MCR}) Problem, which
is fundamental for microbiome analysis. In this problem, the goal is to
reconstruct the identity and frequency of species comprising a microbial
community, using short sequence reads from Massively Parallel Sequencing (MPS)
data obtained for specified genomic regions. We formulate the problem
mathematically as a convex optimization problem and provide sufficient
conditions for identifiability, namely the ability to reconstruct species
identity and frequency correctly when the data size (number of reads) grows to
infinity. We discuss different metrics for assessing the quality of the
reconstructed solution, including a novel phylogenetically-aware metric based
on the Mahalanobis distance, and give upper-bounds on the reconstruction error
for a finite number of reads under different metrics. We propose a scalable
divide-and-conquer algorithm for the problem using convex optimization, which
enables us to handle large problems (with $\sim10^6$ species). We show using
numerical simulations that for realistic scenarios, where the microbial
communities are sparse, our algorithm gives solutions with high accuracy, both
in terms of obtaining accurate frequency, and in terms of species phylogenetic
resolution.


  Recent studies uncovered important core/periphery network structures
characterizing complex sets of cooperative and competitive interactions between
network nodes, be they proteins, cells, species or humans. Better
characterization of the structure, dynamics and function of core/periphery
networks is a key step of our understanding cellular functions, species
adaptation, social and market changes. Here we summarize the current knowledge
of the structure and dynamics of "traditional" core/periphery networks,
rich-clubs, nested, bow-tie and onion networks. Comparing core/periphery
structures with network modules, we discriminate between global and local
cores. The core/periphery network organization lies in the middle of several
extreme properties, such as random/condensed structures, clique/star
configurations, network symmetry/asymmetry, network
assortativity/disassortativity, as well as network hierarchy/anti-hierarchy.
These properties of high complexity together with the large degeneracy of core
pathways ensuring cooperation and providing multiple options of network flow
re-channelling greatly contribute to the high robustness of complex systems.
Core processes enable a coordinated response to various stimuli, decrease
noise, and evolve slowly. The integrative function of network cores is an
important step in the development of a large variety of complex organisms and
organizations. In addition to these important features and several decades of
research interest, studies on core/periphery networks still have a number of
unexplored areas.


  Following transmission, HIV-1 evolves into a diverse population, and next
generation sequencing enables us to detect variants occurring at low
frequencies. Studying viral evolution at the level of whole genomes was
hitherto not possible because next generation sequencing delivers relatively
short reads. We here provide a proof of principle that whole HIV-1 genomes can
be reliably reconstructed from short reads, and use this to study the selection
of immune escape mutations at the level of whole genome haplotypes. Using
realistically simulated HIV-1 populations, we demonstrate that reconstruction
of complete genome haplotypes is feasible with high fidelity. We do not
reconstruct all genetically distinct genomes, but each reconstructed haplotype
represents one or more of the quasispecies in the HIV-1 population. We then
reconstruct 30 whole genome haplotypes from published short sequence reads
sampled longitudinally from a single HIV-1 infected patient. We confirm the
reliability of the reconstruction by validating our predicted haplotype genes
with single genome amplification sequences, and by comparing haplotype
frequencies with observed epitope escape frequencies. Phylogenetic analysis
shows that the HIV-1 population undergoes selection driven evolution, with
successive replacement of the viral population by novel dominant strains. We
demonstrate that immune escape mutants evolve in a dependent manner with
various mutations hitchhiking along with others. As a consequence of this
clonal interference, selection coefficients have to be estimated for complete
haplotypes and not for individual immune escapes.


  The vast majority of strategies aimed at controlling contagion processes on
networks considers the connectivity pattern of the system as either quenched or
annealed. However, in the real world many networks are highly dynamical and
evolve in time concurrently to the contagion process. Here, we derive an
analytical framework for the study of control strategies specifically devised
for time-varying networks. We consider the removal/immunization of individual
nodes according the their activity in the network and develop a block variable
mean-field approach that allows the derivation of the equations describing the
evolution of the contagion process concurrently to the network dynamic. We
derive the critical immunization threshold and assess the effectiveness of the
control strategies. Finally, we validate the theoretical picture by simulating
numerically the information spreading process and control strategies in both
synthetic networks and a large-scale, real-world mobile telephone call dataset


  Drosophila Y chromosomes are composed entirely of silent heterochromatin,
while male X chromosomes have highly accessible chromatin and are
hypertranscribed due to dosage compensation. Here, we dissect the molecular
mechanisms and functional pressures driving heterochromatin formation and
dosage compensation of the recently formed neo-sex chromosomes of Drosophila
miranda. We show that the onset of heterochromatin formation on the neo-Y is
triggered by an accumulation of repetitive DNA. The neo-X has evolved partial
dosage compensation and we find that diverse mutational paths have been
utilized to establish several dozen novel binding consensus motifs for the
dosage compensation complex on the neo-X, including simple point mutations at
pre-binding sites, insertion and deletion mutations, microsatellite expansions,
or tandem amplification of weak binding sites. Spreading of these silencing or
activating chromatin modifications to adjacent regions results in massive
mis-expression of neo-sex linked genes, and little correspondence between
functionality of genes and their silencing on the neo-Y or dosage compensation
on the neo-X. Intriguingly, the genomic regions being targeted by the dosage
compensation complex on the neo-X and those becoming heterochromatic on the
neo-Y show little overlap, possibly reflecting different propensities along the
ancestral chromosome to adopt active or repressive chromatin configurations.
Our findings have broad implications for current models of sex chromosome
evolution, and demonstrate how mechanistic constraints can limit evolutionary
adaptations. Our study also highlights how evolution can follow predictable
genetic trajectories, by repeatedly acquiring the same 21-bp consensus motif
for recruitment of the dosage compensation complex, yet utilizing a diverse
array of random mutational changes to attain the same phenotypic outcome.


  Chromatin is a highly compact and dynamic nuclear structure that consists of
DNA and associated proteins. The main organizational unit is the nucleosome,
which consists of a histone octamer with DNA wrapped around it. Histone
proteins are implicated in the regulation of eukaryote genes and they carry
numerous reversible post-translational modifications that control DNA-protein
interactions and the recruitment of chromatin binding proteins.
Heterochromatin, the transcriptionally inactive part of the genome, is densely
packed and contains histone H3 that is methylated at Lys 9 (H3K9me). The
propagation of H3K9me in nucleosomes along the DNA in chromatin is antagonizing
by methylation of H3 Lysine 4 (H3K4me) and acetylations of several lysines,
which is related to euchromatin and active genes. We show that the related
histone modifications form antagonized domains on a coarse scale. These histone
marks are assumed to be initiated within distinct nucleation sites in the DNA
and to propagate bi-directionally. We propose a simple computer model that
simulates the distribution of heterochromatin in human chromosomes. The
simulations are in agreement with previously reported experimental observations
from two different human cell lines. We reproduced different types of barriers
between heterochromatin and euchromatin providing a unified model for their
function. The effect of changes in the nucleation site distribution and of
propagation rates were studied. The former occurs mainly with the aim of
(de-)activation of single genes or gene groups and the latter has the power of
controlling the transcriptional programs of entire chromosomes. Generally, the
regulatory program of gene transcription is controlled by the distribution of
nucleation sites along the DNA string.


  The analysis of the structure of chemical reaction networks is crucial for a
better understanding of chemical processes. Such networks are well described as
hypergraphs. However, due to the available methods, analyses regarding network
properties are typically made on standard graphs derived from the full
hypergraph description, e.g.\ on the so-called species and reaction graphs.
However, a reconstruction of the underlying hypergraph from these graphs is not
necessarily unique. In this paper, we address the problem of reconstructing a
hypergraph from its species and reaction graph and show NP-completeness of the
problem in its Boolean formulation. Furthermore we study the problem
empirically on random and real world instances in order to investigate its
computational limits in practice.


  Human mobility is a key component of large-scale spatial-transmission models
of infectious diseases. Correctly modeling and quantifying human mobility is
critical for improving epidemic control policies, but may be hindered by
incomplete data in some regions of the world. Here we explore the opportunity
of using proxy data or models for individual mobility to describe commuting
movements and predict the diffusion of infectious disease. We consider three
European countries and the corresponding commuting networks at different
resolution scales obtained from official census surveys, from proxy data for
human mobility extracted from mobile phone call records, and from the radiation
model calibrated with census data. Metapopulation models defined on the three
countries and integrating the different mobility layers are compared in terms
of epidemic observables. We show that commuting networks from mobile phone data
well capture the empirical commuting patterns, accounting for more than 87% of
the total fluxes. The distributions of commuting fluxes per link from both
sources of data - mobile phones and census - are similar and highly correlated,
however a systematic overestimation of commuting traffic in the mobile phone
data is observed. This leads to epidemics that spread faster than on census
commuting networks, however preserving the order of infection of newly infected
locations. Match in the epidemic invasion pattern is sensitive to initial
conditions: the radiation model shows higher accuracy with respect to mobile
phone data when the seed is central in the network, while the mobile phone
proxy performs better for epidemics seeded in peripheral locations. Results
suggest that different proxies can be used to approximate commuting patterns
across different resolution scales in spatial epidemic simulations, in light of
the desired accuracy in the epidemic outcome under study.


  All dynamical systems of biological interest--be they food webs, regulation
of genes, or contacts between healthy and infectious individuals--have complex
network structure. Wigner's semicircular law and Girko's circular law describe
the eigenvalues of systems whose structure is a fully connected network.
However, these laws fail for systems with complex network structure. Here we
show that in these cases the eigenvalues are described by superellipses. We
also develop a new method to analytically estimate the dominant eigenvalue of
complex networks.


  Gene-regulatory enhancers have been identified by many lines of evidence,
including evolutionary conservation, regulatory protein binding, chromatin
modifications, and DNA sequence motifs. To integrate these different
approaches, we developed EnhancerFinder, a novel method for predicting
developmental enhancers and their tissue specificity. EnhancerFinder uses a
two-step multiple-kernel learning approach to integrate DNA sequence motifs,
evolutionary patterns, and thousands of diverse functional genomics datasets
from a variety of cell types and developmental stages. We trained
EnhancerFinder on hundreds of experimentally verified human developmental
enhancers from the VISTA Enhancer Browser, in contrast to histone mark or
sequence-based enhancer definitions commonly used. We comprehensively evaluated
EnhancerFinder, and found that our integrative approach improves enhancer
prediction accuracy over previous approaches that consider a single type of
data. Our evaluation highlights the importance of considering information from
many tissues when predicting specific types of enhancers. We find that VISTA
enhancers active in embryonic heart are easier to predict than enhancers active
in several other tissues due to their uniquely high GC content. We applied
EnhancerFinder to the entire human genome and predicted 84,301 developmental
enhancers and their tissue specificity. These predictions provide specific
functional annotations for large amounts of human non-coding DNA, and are
significantly enriched near genes with annotated roles in their predicted
tissues and hits from genome-wide association studies. We demonstrate the
utility of our enhancer predictions by identifying and validating a novel
cranial nerve enhancer in the ZEB2 locus. Our genome-wide developmental
enhancer predictions will be freely available as a UCSC Genome Browser track.


  Horseshoe crabs are marine arthropods with a fossil record extending back
approximately 450 million years. They exhibit remarkable morphological
stability over their long evolutionary history, retaining a number of ancestral
arthropod traits, and are often cited as examples of "living fossils." As
arthropods, they belong to the Ecdysozoa}, an ancient super-phylum whose
sequenced genomes (including insects and nematodes) have thus far shown more
divergence from the ancestral pattern of eumetazoan genome organization than
cnidarians, deuterostomes, and lophotrochozoans. However, much of ecdysozoan
diversity remains unrepresented in comparative genomic analyses. Here we use a
new strategy of combined de novo assembly and genetic mapping to examine the
chromosome-scale genome organization of the Atlantic horseshoe crab Limulus
polyphemus. We constructed a genetic linkage map of this 2.7 Gbp genome by
sequencing the nuclear DNA of 34 wild-collected, full-sibling embryos and their
parents at a mean redundancy of 1.1x per sample. The map includes 84,307
sequence markers and 5,775 candidate conserved protein coding genes. Comparison
to other metazoan genomes shows that the L. polyphemus genome preserves
ancestral bilaterian linkage groups, and that a common ancestor of modern
horseshoe crabs underwent one or more ancient whole genome duplications (WGDs)
~ 300 MYA, followed by extensive chromosome fusion.


  A model of the mechanisms underlying exploratory behaviour, based on
empirical research and refined using a computer simulation, is presented. The
behaviour of killifish from two lakes, one with killifish predators and one
without, was compared in the laboratory. Plotting average activity in a novel
environment versus time resulted in an inverted-U-shaped curve for both groups;
however, the curve for killifish from the lake without predators was (1)
steeper, (2) reached a peak value earlier, (S) reached a higher peak value, and
(4) subsumed less area than the curve for killifish from the lake with
predators. We hypothesize that the shape of the exploration curve reflects a
competition between motivational subsystems that excite and inhibit exploratory
behaviour in a way that is tuned to match the affordance probabilities of the
animal's environment. A computer implementation of this model produced curves
which differed along the same four dimensions as differentiate the two
killifish curves. All four differences were reproduced in the model by tuning a
single parameter: the time-dependent component of the decay-rate of the
exploration-inhibiting subsystem.


  The heat shock organizing protein (Hop) is important in modulating the
activity and co-interaction of two chaperones: heat shock protein 70 and 90
(Hsp70 and Hsp90). Recent research suggested that Plasmodium falciparum Hop
(PfHop), PfHsp70 and PfHsp90 form a complex in the trophozoite infective stage.
However, there has been little computational research on the malarial Hop
protein in complex with other malarial Hsps. Using in silico characterization
of the protein, this work showed that individual domains of Hop are evolving at
different rates within the protein. Differences between human Hop (HsHop) and
PfHop were identified by motif analysis. Homology modeling of PfHop and HsHop
in complex with their own cytosolic Hsp90 and Hsp70 C-terminal peptide partners
indicated excellent conservation of the Hop concave TPR sites bound to the
C-terminal motifs of partner proteins. Further, we analyzed additional binding
sites between Hop and Hsp90, and showed, for the first time, that they are
distinctly less conserved between human and malaria parasite. These sites are
located on the convex surface of Hop TPR2, and involved in interactions with
the Hsp90 middle domain. Since the convex sites are less conserved than the
concave sites, it makes their potential for malarial inhibitor design extremely
attractive (as opposed to the concave sites which have been the focus of
previous efforts).


  We report a theoretical study of DNA flexibility and quantitatively predict
the ring closure probability as a function of DNA contour length. Recent
experimental studies show that the flexibility of short DNA fragments (as
compared to the persistence length of DNA l_P~150 base pairs) cannot be
described by the traditional worm-like chain (WLC) model, e.g., the observed
ring closure probability is much higher than predicted. To explain these
observations, DNA flexibility is investigated with explicit considerations of a
new length scale l_D~10 base pairs, over which DNA local bend angles are
correlated. In this correlated worm-like chain (C-WLC) model, a finite length
correction term is analytically derived and the persistence length is found to
be contour length dependent. While our model reduces to the traditional
worm-like chain model when treating long DNA at length scales much larger than
l_P, it predicts that DNA becomes much more flexible at shorter sizes, which
helps explain recent cyclization measurements of short DNA fragments around 100
base pairs.


  We formulate a mathematical model of competition for resources between
representatives of different age groups. A nonlinear kinetic
integral-differential equation of the age aggression describes the process of
redistribution of resources. It is shown that the equation of the age
aggression has a stationary solution, in the absence of age-dependency in the
interaction of different age groups. A numerical simulation of the evolution of
resources for different initial distributions has done. It is shown the
instability of the system and the existence of regimes leading to a
concentration of resources in the certain age groups.


  Molecular motors are involved in key transport processes inside actin-based
cellular protrusions. The motors carry cargo proteins to the protrusion tip
which participate in regulating the actin polymerization, and play a key role
in facilitating the growth and formation of such protrusions. It is observed
that the motors accumulate at the tips of cellular protrusions, and in addition
form aggregates that are found to drift towards the protrusion base at the rate
of actin treadmilling. We present a one-dimensional driven lattice model, where
motors become inactive after delivering their cargo at the tip, or by loosing
their cargo to a cargo-less neighbor. The results suggest that the experimental
observations may be explained by the formation of traffic jams that form at the
tip. The model is solved using a novel application of mean-field and shock
analysis. We find a new class of shocks that undergo intermittent collapses,
and on average do not obey the Rankine-Hugoniot relation.


  We introduce a new rotationally invariant viewing angle classification method
for identifying, among a large number of Cryo-EM projection images, similar
views without prior knowledge of the molecule. Our rotationally invariant
features are based on the bispectrum. Each image is denoised and compressed
using steerable principal component analysis (PCA) such that rotating an image
is equivalent to phase shifting the expansion coefficients. Thus we are able to
extend the theory of bispectrum of 1D periodic signals to 2D images. The
randomized PCA algorithm is then used to efficiently reduce the dimensionality
of the bispectrum coefficients, enabling fast computation of the similarity
between any pair of images. The nearest neighbors provide an initial
classification of similar viewing angles. In this way, rotational alignment is
only performed for images with their nearest neighbors. The initial nearest
neighbor classification and alignment are further improved by a new
classification method called vector diffusion maps. Our pipeline for viewing
angle classification and alignment is experimentally shown to be faster and
more accurate than reference-free alignment with rotationally invariant K-means
clustering, MSA/MRA 2D classification, and their modern approximations.


  Over the last years, analyses performed on a stochastic model of catalytic
reaction networks have provided some indications about the reasons why wet-lab
experiments hardly ever comply with the phase transition typically predicted by
theoretical models with regard to the emergence of collectively
self-replicating sets of molecule (also defined as autocatalytic sets, ACSs), a
phenomenon that is often observed in nature and that is supposed to have played
a major role in the emergence of the primitive forms of life. The model at
issue has allowed to reveal that the emerging ACSs are characterized by a
general dynamical fragility, which might explain the difficulty to observe them
in lab experiments. In this work, the main results of the various analyses are
reviewed, with particular regard to the factors able to affect the generic
properties of catalytic reactions network, for what concerns, not only the
probability of ACSs to be observed, but also the overall activity of the
system, in terms of production of new species, reactions and matter.


  The bottom-up construction of synthetic cells is one of the most intriguing
and interesting research arenas in synthetic biology. Synthetic cells are built
by encapsulating biomolecules inside lipid vesicles (liposomes), allowing the
synthesis of one or more functional proteins. Thanks to the in situ synthesized
proteins, synthetic cells become able to perform several biomolecular
functions, which can be exploited for a large variety of applications. This
paves the way to several advanced uses of synthetic cells in basic science and
biotechnology, thanks to their versatility, modularity, biocompatibility, and
programmability. In the previous WIVACE (2012) we presented the
state-of-the-art of semi-synthetic minimal cell (SSMC) technology and
introduced, for the first time, the idea of chemical communication between
synthetic cells and natural cells. The development of a proper synthetic
communication protocol should be seen as a tool for the nascent field of
bio/chemical-based Information and Communication Technologies (bio-chem-ICTs)
and ultimately aimed at building soft-wet-micro-robots. In this contribution
(WIVACE, 2013) we present a blueprint for realizing this project, and show some
preliminary experimental results. We firstly discuss how our research goal
(based on the natural capabilities of biological systems to manipulate chemical
signals) finds a proper place in the current scientific and technological
contexts. Then, we shortly comment on the experimental approaches from the
viewpoints of (i) synthetic cell construction, and (ii) bioengineering of
microorganisms, providing up-to-date results from our laboratory. Finally, we
shortly discuss how autopoiesis can be used as a theoretical framework for
defining synthetic minimal life, minimal cognition, and as bridge between
synthetic biology and artificial intelligence.


  Epigenetic Tracking (ET) is an Artificial Embryology system which allows for
the evolution and development of large complex structures built from artificial
cells. In terms of the number of cells, the complexity of the bodies generated
with ET is comparable with the complexity of biological organisms. We have
previously used ET to simulate the growth of multicellular bodies with
arbitrary 3-dimensional shapes which perform computation using the paradigm of
"metabolic computing". In this paper we investigate the memory capacity of such
computational structures and analyse the trade-off between shape and
computation. We now plan to build on these foundations to create a
biologically-inspired model in which the encoding of the phenotype is efficient
(in terms of the compactness of the genome) and evolvable in tasks involving
non-trivial computation, robust to damage and capable of self-maintenance and
self-repair.


  Reactions forming a pathway can be rewritten by making explicit the different
molecular components involved in them. A molecular component represents a
biological entity (e.g. a protein) in all its states (free, bound, degraded,
etc.). In this paper we show the application of a component identification
algorithm to a number of real-world models to experimentally validate the
approach. Components identification allows subpathways to be computed to better
understand the pathway functioning.


  In this work we introduce some preliminary analyses on the role of a
semi-permeable membrane in the dynamics of a stochastic model of catalytic
reaction sets (CRSs) of molecules. The results of the simulations performed on
ensembles of randomly generated reaction schemes highlight remarkable
differences between this very simple protocell description model and the
classical case of the continuous stirred-tank reactor (CSTR). In particular, in
the CSTR case, distinct simulations with the same reaction scheme reach the
same dynamical equilibrium, whereas, in the protocell case, simulations with
identical reaction schemes can reach very different dynamical states, despite
starting from the same initial conditions.


  The Spatial Processes package enables an explicit definition of a spatial
environment on top of the normal dynamic modeling SBML capabilities. The
possibility of an explicit representation of spatial dynamics increases the
representation power of SBML. In this work we used those new SBML features to
define an extensive model of colonic crypts composed of the main cellular types
(from stem cells to fully differentiated cells), alongside their spatial
dynamics.


  An approach that combines Self-Organizing maps, hierarchical clustering and
network components is presented, aimed at comparing protein conformational
ensembles obtained from multiple Molecular Dynamic simulations. As a first
result the original ensembles can be summarized by using only the
representative conformations of the clusters obtained. In addition the network
components analysis allows to discover and interpret the dynamic behavior of
the conformations won by each neuron. The results showed the ability of this
approach to efficiently derive a functional interpretation of the protein
dynamics described by the original conformational ensemble, highlighting its
potential as a support for protein engineering.


  An extensive rewiring of cell metabolism supports enhanced proliferation in
cancer cells. We propose a systems level approach to describe this phenomenon
based on Flux Balance Analysis (FBA). The approach does not explicit a cell
biomass formation reaction to be maximized, but takes into account an ensemble
of alternative flux distributions that match the cancer metabolic rewiring
(CMR) phenotype description. The underlying concept is that the analysis the
common/distinguishing properties of the ensemble can provide indications on how
CMR is achieved and sustained and thus on how it can be controlled.


  In this paper we motivate and develop a model of gene expression for the
purpose of studying the interaction between synthetic gene circuits and the
chassis cell within which they are in- serted. This model focuses on the
translational aspect of gene expression as this is where the literature
suggests the crucial interaction between gene expression and shared resources
lies.


  We consider a model for substrate-depletion oscillations in genetic systems,
based on a stochastic differential equation with a slowly evolving external
signal. We show the existence of critical transitions in the system. We apply
two methods to numerically test the synthetic time series generated by the
system for early indicators of critical transitions: a detrended fluctuation
analysis method, and a novel method based on topological data analysis
(persistence diagrams).


  This paper will introduce a theory of emergent animal social complexity using
various results from computational models and empirical results. These results
will be organized into a vertical model of social complexity. This will support
the perspective that social complexity is in essence an emergent phenomenon
while helping to answer two interrelated questions. The first of these involves
how behavior is integrated at units of analysis larger than the individual
organism. The second involves placing aggregate social events into the context
of processes occurring within individual organisms over time (e.g. genomic and
physiological processes). By using a complex systems perspective, five
principles of social complexity can be identified. These principles suggest
that lower-level mechanisms give rise to high-level mechanisms, ultimately
resulting in metastable networks of social relations. These network structures
then constrain lower-level phenomena ranging from transient, collective social
groups to physiological regulatory mechanisms within individual organisms. In
conclusion, the broader implications and drawbacks of applying the theory to a
diversity of natural populations will be discussed.


  Inferences about the evolution of continuous traits based on reconstruction
of ancestral states has often been considered more error-prone than analysis of
independent contrasts. Here we show that both methods in fact yield identical
estimators for the correlation coefficient and regression gradient of
correlated traits, indicating that reconstructed ancestral states are a valid
source of information about correlated evolution. We show that the independent
contrast associated with a pair of sibling nodes on a phylogenetic tree can be
expressed in terms of the maximum likelihood ancestral state function at those
nodes and their common parent. This expression gives rise to novel formulae for
independent contrasts for any model of evolution admitting of a local
likelihood function. We thus derive new formulae for independent contrasts
applicable to traits evolving under directional drift, and use simulated data
to show that these directional contrasts provide better estimates of
evolutionary model parameters than standard independent contrasts, when traits
in fact evolve with a directional tendency.


  The interpretation of count data originating from the current generation of
DNA sequencing platforms requires special attention. In particular, the
per-sample library sizes often vary by orders of magnitude from the same
sequencing run, and the counts are overdispersed relative to a simple Poisson
model These challenges can be addressed using an appropriate mixture model that
simultaneously accounts for library size differences and biological
variability. This approach is already well-characterized and implemented for
RNA-Seq data in R packages such as edgeR and DESeq.
  We use statistical theory, extensive simulations, and empirical data to show
that variance stabilizing normalization using a mixture model like the negative
binomial is appropriate for microbiome count data. In simulations detecting
differential abundance, normalization procedures based on a Gamma-Poisson
mixture model provided systematic improvement in performance over crude
proportions or rarefied counts -- both of which led to a high rate of false
positives. In simulations evaluating clustering accuracy, we found that the
rarefying procedure discarded samples that were nevertheless accurately
clustered by alternative methods, and that the choice of minimum library size
threshold was critical in some settings, but with an optimum that is unknown in
practice. Techniques that use variance stabilizing transformations by modeling
microbiome count data with a mixture distribution, such as those implemented in
edgeR and DESeq, substantially improved upon techniques that attempt to
normalize by rarefying or crude proportions. Based on these results and
well-established statistical theory, we advocate that investigators avoid
rarefying altogether. We have provided microbiome-specific extensions to these
tools in the R package, phyloseq.


  Chenopodium album seedling emergence studies were conducted at nine European
and two North American locations comparing local populations with a common
population from Denmark. It is hypothesized that C. album seedling recruitment
timing and magnitude have adapted to environmental and cropping system
practices of a locality. Limitations in the habitat (filter 1) were reflected
in local C. album population recruitment season length. Generally, the duration
of seedling recruitment of both populations (local; DEN-COM) increased with
decreasing latitude, north-to-south. In general, compared to the local
population, DEN-COM recruitment at locations north of Denmark was longer and
south of Denmark was shorter, and ended sooner. Generally, the local cropping
system disturbances (CSD) period increased with decreasing latitude. The total
duration of the CSD period was over twice as long in the south as that in the
north. Recruitment at each locality possessed seasonal structure (time, number)
consisting of 2-4 discrete seasonal cohorts. This may be an adaptive means by
which C. album searches for, and exploits, recruitment opportunity just prior
to, and after, predictable disturbances. The control of C. album seedling
emergence is contained in the heteroblastic traits of its locally adapted
seeds, and is stimulated by a complex interaction of light, heat, water,
nitrate and oxygen signals inherent in the local environment. Our observations
of complex recruitment patterns occurring at critical cropping times is strong
evidence that C. album possesses a flexible and sensitive germination
regulation system adaptable to opportunity in many different Eurasian and North
American agricultural habitats.


  In this paper we unify, simplify, and extend previous work on the
evolutionary dynamics of symmetric $N$-player matrix games with two pure
strategies. In such games, gains from switching strategies depend, in general,
on how many other individuals in the group play a given strategy. As a
consequence, the gain function determining the gradient of selection can be a
polynomial of degree $N-1$. In order to deal with the intricacy of the
resulting evolutionary dynamics, we make use of the theory of polynomials in
Bernstein form. This theory implies a tight link between the sign pattern of
the gains from switching on the one hand and the number and stability
properties of the rest points of the replicator dynamics on the other hand.
While this relationship is a general one, it is most informative if gains from
switching have at most two sign changes, as it is the case for most
multi-player matrix games considered in the literature. We demonstrate that
previous results for public goods games are easily recovered and extended using
this observation. Further examples illustrate how focusing on the sign pattern
of the gains from switching obviates the need for a more involved analysis.


  The advent of molecular biology has led to the identification of definitive
causative factors for a number of diseases, most of which are monogenic. Causes
for most common diseases across the population, however, seem elusive and
cannot be pinpointed to a limited number of genes or genetic pathways. This
realization has led to the idea of personalized medicine and treating each case
individually. Nevertheless, since each common disease appears to have the same
endpoint and phenotypic features in all diagnosed individuals, the search for a
unifying cause will still continue. Given that multivariate scientific data is
of a correlative nature and causation is always inferred, a simple
formalization of the general structure of cause and correlation is presented
herein. Furthermore, the context in which a causal structure could take shape,
termed the 'pericause', is proposed as a tractable and uninvestigated concept
which could theoretically play a crucial role in determining the effects of a
cause.


  The phenomenon of preadaptation, or exaptation (wherein a trait that
originally evolved to solve one problem is co-opted to solve a new problem)
presents a formidable challenge to efforts to describe biological phenomena
using a classical (Kolmogorovian) mathematical framework. We develop a quantum
framework for exaptation with examples from both biological and cultural
evolution. The state of a trait is written as a linear superposition of a set
of basis states, or possible forms the trait could evolve into, in a complex
Hilbert space. These basis states are represented by mutually orthogonal unit
vectors, each weighted by an amplitude term. The choice of possible forms
(basis states) depends on the adaptive function of interest (e.g., ability to
metabolize lactose or thermoregulate), which plays the role of the observable.
Observables are represented by self-adjoint operators on the Hilbert space. The
possible forms (basis states) corresponding to this adaptive function
(observable) are called eigenstates. The framework incorporates key features of
exaptation: potentiality, contextuality, nonseparability, and emergence of new
features. However, since it requires that one enumerate all possible contexts,
its predictive value is limited, consistent with the assertion that there
exists no biological equivalent to "laws of motion" by which we can predict the
evolution of the biosphere.


  The Lenski experiment is a long term daily reproduction of Escherichia coli,
that has evidenced phenotypic and genetic evolutions along the years. Some
mathematical models, that could be usefull in understanding the results of that
experiment, are reviewed here: stochastic and deterministic growth, mutation
appearance and fixation, competition of species.


  This limited review is intended as an introduction to the fast growing
subject of mathematical modelling of cell metabolism and its biochemical
pathways, and more precisely on pathways linked to apoptosis of cancerous
cells. Some basic mathematical models of chemical kinetics, with emphasis on
stochastic models, are presented.


  The application of conventional phylogenetic techniques for inferring
cultural history is problematic due to differences in the nature of information
transmission in biological and cultural realms. In culture, units of
transmission are not just measurable attributes, but communicable concepts.
Therefore, relatedness amongst cultural elements often resides at the
conceptual level not captured by traditional phylogenetic methods. This paper
takes a cognitively inspired approach to analyzing material cultural history.
We show that combining data for physical attributes of cultural artifacts with
conceptual information can uncover cultural influences among different
ethnolinguistic groups, and reveal new patterns of cultural ancestry. Using the
Baltic psaltery, a musical instrument with a well-documented ethnographic and
archaeological record, we recovered a previously unacknowledged pattern of
historical relationship that is more congruent with geographical distribution
and temporal data than is obtained with other approaches.


  An intrinsically disordered protein (IDP) lacks a stable three-dimensional
structure, while it folds into a specific structure when it binds to a target
molecule. In some IDP-target complexes, not all target binding surfaces are
exposed on the outside, and intermediate states are observed in their binding
processes. We consider that stepwise target recognition via intermediate states
is a characteristic of IDP binding to targets with "hidden" binding sites. To
investigate IDP binding to hidden target binding sites, we constructed an IDP
lattice model based on the HP model. In our model, the IDP is modeled as a
chain and the target is modeled as a highly coarse-grained object. We
introduced motion and internal interactions to the target to hide its binding
sites. In the case of unhidden binding sites, a two-state transition between
the free states and a bound state is observed, and we consider that this
represents coupled folding and binding. Introducing hidden binding sites, we
found an intermediate bound state in which the IDP forms various structures to
temporarily stabilize the complex. The intermediate state provides a scaffold
for the IDP to access the hidden binding site. We call this process multiform
binding. We conclude that structural flexibility of IDPs enables them to access
hidden binding sites, and this is a functional advantage of IDPs.


  The increased availability of time series genetic variation data from
experimental evolution studies and ancient DNA samples has created new
opportunities to identify genomic regions under selective pressure and to
estimate their associated fitness parameters. However, it is a challenging
problem to compute the likelihood of nonneutral models for the population
allele frequency dynamics, given the observed temporal DNA data. Here, we
develop a novel spectral algorithm to analytically and efficiently integrate
over all possible frequency trajectories between consecutive time points. This
advance circumvents the limitations of existing methods which require
fine-tuning the discretization of the population allele frequency space when
numerically approximating requisite integrals. Furthermore, our method is
flexible enough to handle general diploid models of selection where the
heterozygote and homozygote fitness parameters can take any values, while
previous methods focused on only a few restricted models of selection. We
demonstrate the utility of our method on simulated data and also apply it to
analyze ancient DNA data from genetic loci associated with coat coloration in
horses. In contrast to previous studies, our exploration of the full fitness
parameter space reveals that a heterozygote advantage form of balancing
selection may have been acting on these loci.


  We cast the problem of protein transfer free energy within the formalism of
density functional theory (DFT), treating the protein as a source of external
potential that acts upon the solvent. Solvent excluded volume,
solvent-accessible surface area, and temperature-dependence of the transfer
free energy all emerge naturally within this formalism, and may be compared
with simplified "back of the envelope" models, which are also developed here.
Depletion contributions to osmolyte induced stability range from 5-$10\kB T$
for typical protein lengths. The general DFT transfer theory developed here may
be simplified to reproduce a Langmuir isotherm condensation mechanism on the
protein surface in the limits of short-ranged interactions, and dilute solute.
Extending the equation of state to higher solute densities results in
non-monotonic behavior of the free energy driving protein or polymer collapse.
Effective interaction potentials between protein backbone or sidechains and
TMAO are obtained, assuming a simple backbone/sidechain 2-bead model for the
protein with an effective 6-12 potential with the osmolyte. The transfer free
energy $\delta g$ shows significant entropy: $d(\delta g)/dT \approx 20 \kB$
for a 100 residue protein. The application of DFT to effective solvent forces
for use in implicit-solvent molecular dynamics is also developed. The simplest
DFT expressions for implicit-solvent forces contain both depletion interactions
and an "impeded-solvation" repulsive force at larger distances.


  The newly appeared idea of Logical Stochastic Resonance (LSR) is verified in
synthetic gene networks induced by non-Gaussian noise. We realize the switching
between two kinds of logic gates under optimal moderate noise intensity by
varying two different tunable parameters in a single gene network. Furthermore,
in order to obtain more logic operations, thus providing additional information
processing capacity, we obtain in a two-dimensional toggle switch model two
complementary logic gates and realize the transformation between two logic
gates via the methods of changing different parameters. These simulated results
contribute to improve the computational power and functionality of the
networks.


  Predator-prey coevolution is commonly thought to result in reciprocal arms
races that produce increasingly extreme and complex traits. However, such
directional change is not inevitable. Here, we provide evidence for a
previously undemonstrated dynamic that we call "strategy chases," wherein
populations explore strategies with similar levels of complexity, but differing
behaviorally. Indeed, in populations of evolving digital organisms, as prey
evolved more effective predator-avoidance strategies, they explored a wider
range of behavioral strategies in addition to exhibiting increased levels of
behavioral complexity. Furthermore, coevolved prey became more adept in
foraging, evidently through coopting components of explored sense-and-flee
avoidance strategies into sense-and-retrieve foraging strategies. Specifically,
we demonstrate that coevolution induced non-escalating exploration of
behavioral space, corresponding with significant evolutionary advancements,
including increasingly intelligent behavioral strategies.


  Mapping between sequence and structure is currently an open problem in
structural biology. Despite many experimental and computational efforts it is
not clear yet how the structure is encoded in the sequence. Answering this
question may pave the way for predicting a protein fold given its sequence.
  My doctoral studies have focused on a particular phenomenon relevant to the
protein sequence-structure relationship. It has been observed that many
proteins having apparently dissimilar sequences share the same native fold. The
phenomenon of mapping many divergent sequences into a single fold raises the
question of which positions along the sequence are important for the
conservation of fold and function in dissimilar sequences. Detecting those
positions, and classifying them according to role can help understand which
elements in a sequence are important for maintenance of structure and/or
function. In the course of my doctoral research I have attempted to discover
and characterize those positions.


  We study the effects of the shear force on the rupture mechanism on a double
stranded DNA. Motivated by recent experiments, we perform the atomistic
simulations with explicit solvent to obtain the distributions of extension in
hydrogen and covalent bonds below the rupture force. We obtain a significant
difference between the atomistic simulations and the existing results in the
iterature based on the coarse-grained models (theory and simulations). We
discuss the possible reasons and improve the coarse-grained model by
incorporating the consequences of semi-microscopic details of the nucleotides
in its description. The distributions obtained by the modified model
(simulations and theoretical) are qualitatively similar to the one obtained
using atomistic simulations.


  Nitrogen isotopic hetero- and homogeneity of D-{\alpha}-alanine and
L-{\alpha}-alanine were investigated in microbial processes in the domain
Bacteria and in chemical processes in symmetric organic synthesis. D-alanine is
an enantiomer that is physiologically essential for microbial growth and
metabolic maintenance. The nitrogen isotopic difference {\Delta}15ND-L (defined
as {\delta}15ND-Ala - {\delta}15NL-Ala) in peptidoglycan amino acids in
bacteria such as the representative gram-positive phyla Firmicutes and
Actinobacteria (Enterococcus faecalis, Staphylococcus aureus, Staphylococcus
staphylolyticus, Lactobacillus acidophilus, Bacillus subtilis, Micrococcus
luteus, and Streptomyces sp.) tended to be 15N-depleted in D-alanine
({\Delta}15ND-L < -2.0 permil). These results suggest that the composition of
isotopically heterogeneous components in these bacteria is primarily controlled
by enzymatic pathways prior to formation of the bacterial cell wall. In
contrast, the {\Delta}15ND-L of racemic alanine in the chemical pathway during
the nucleophilic substitution reaction (SN1 type) between 2-bromopropionic acid
and ammonia identified fully homogeneous components for each enantiomer. The
novel enantiomer-specific isotopic analysis (ESIA) method is useful in
determining the origins of chirality in biogenic and abiogenic processes and is
applicable to enantiomer studies.
  -- Keywords: D-alanine, microbial process, chemical process, nitrogen
isotopic composition, enantiomer-specific isotope analysis (ESIA)


  Proteins are biomolecules of life. They fold into a great variety of
three-dimensional (3D) shapes. Underlying these folding patterns are many
recurrent structural fragments or building blocks (analogous to `LEGO bricks').
This paper reports an innovative statistical inference approach to discover a
comprehensive dictionary of protein structural building blocks from a large
corpus of experimentally determined protein structures. Our approach is built
on the Bayesian and information-theoretic criterion of minimum message length.
To the best of our knowledge, this work is the first systematic and rigorous
treatment of a very important data mining problem that arises in the
cross-disciplinary area of structural bioinformatics. The quality of the
dictionary we find is demonstrated by its explanatory power -- any protein
within the corpus of known 3D structures can be dissected into successive
regions assigned to fragments from this dictionary. This induces a novel
one-dimensional representation of three-dimensional protein folding patterns,
suitable for application of the rich repertoire of character-string processing
algorithms, for rapid identification of folding patterns of newly-determined
structures. This paper presents the details of the methodology used to infer
the dictionary of building blocks, and is supported by illustrative examples to
demonstrate its effectiveness and utility.


  Before a cell divides into two daughter cells, chromosomes are replicated
resulting in two sister chromosomes embracing each other. Each sister
chromosome is bound to a separate proteinous structure, called kinetochore
(kt), that captures the tip of a filamentous protein, called microtubule (MT).
Two oppositely oriented MTs pull the two kts attached to two sister chromosomes
thereby pulling the two sisters away from each other. Here we theoretically
study an even simpler system, namely an isolated kt coupled to a single MT;
this system mimics an {\it in-vitro} experiment where a single kt-MT attachment
is reconstituted using purified extracts from budding yeast. Our models not
only account for the experimentally observed "catch-bond-like" behavior of the
kt-MT coupling, but also make new predictions on the probability distribution
of the lifetimes of the attachments. In principle, our new predictions can be
tested by analyzing the data collected in the {\it in-vitro} experiments
provided the experiment is repeated sufficiently large number of times. Our
theory provides a deep insight into the effects of (a) size, (b) energetics,
and (c) stochastic kinetics of the kt-MT coupling on the distribution of the
lifetimes of these attachments.


  Variable individual response to epidemics may be found within many contexts
in the study of infectious diseases (e.g., age structure or contact networks).
There are situations where the variability, in terms of epidemiological
parameter, cannot be neatly packaged along with other demographics of the
population like spatial location or life stage. Transport equations are a novel
method for handling this variability via a distributed parameter; where
particular parameter values are possessed by various proportions of the
population. Several authors (e.g., Kareva, Novozhilov, and Katriel) have
studied such systems in a closed population setting (no births/immigrations or
deaths/emigrations), but have cited restrictions to employing such methods when
entry and removal of individuals is added to the population. This paper
details, in the context of a simple susceptible-infectious-recovered (SIR)
epidemic, how the method works in the closed population setting and gives
conditions for initial, transient, and asymptotic results to be equivalent with
the nondistributed case. Additionally, I show how the method may be applied to
various forms of open SIR systems. Transport equations are used to transform an
infinite dimensional system for the open population case into a finite
dimensional system which is, at the very least, able to be numerically studied,
a model with direct inheritance of the distributed parameter is shown to be
qualitatively identical to the nondistributed case, and finally a model where
disease results in sterilization is fully analyzed.


  Over the past two decades, a series of works have aimed at studying the
problem of genome assembly: the process of reconstructing a genome from
sequence reads. An early formulation of the genome assembly problem showed that
genome reconstruction is NP-hard when framed as finding the shortest sequence
that contains all observed reads. Although this original formulation is very
simplistic and does not allow for mate-pair information, subsequent
formulations have also proven to be NP-hard, and/or may not be guaranteed to
return a correct assembly.
  In this paper, we provide an alternate perspective on the genome assembly
problem by showing genome assembly is easy when provided with sufficient
mate-pair information. Moreover, we quantify the number of mate-pair libraries
necessary and sufficient for accurate genome assembly, in terms of the length
of the longest repetitive region within a genome. In our analysis, we consider
an idealized sequencing model where each mate-pair library generates pairs of
error free reads with a fixed and known insert size at each position in the
genome.
  Even in this idealized model, we show that accurate genome reconstruction
cannot be guaranteed in the worst case unless at least roughly R/2L mate-pair
libraries are produced, where R is the length of the longest repetitive region
in the genome and L is the length of each read. On the other hand, if (R/L)+1
mate-pair libraries are provided, then a simple algorithm can be used to find a
correct genome assembly easily in polynomial time. Although (R/L)+1 mate-pair
libraries can be too much to produce in practice, the previous bounds only hold
in the worst case. In our last result, we show that if additional conditions
hold on a genome, a correct assembly can be guaranteed with only O(log (R/L))
mate-pair libraries.


  Simonton (2006) makes the unwarranted assumption that nonmonotonicity
supports a Darwinian view of creativity. Darwin's theory of natural selection
was motivated by a paradox that has no equivalent in creative thought: the
paradox of how change accumulates when acquired traits are not inherited. To
describe a process of cumulative change in which acquired traits are retained
is outside of the scope of the theory of natural selection. Even the early
evolution of life itself (prior to genetically mediated template replication)
cannot be described by natural selection. Specifically, natural selection
cannot describe change of state that involves horizontal (Lamarckian) exchange,
or occurs through interaction with an incompletely specified context. It cannot
describe change wherein variants are evaluated sequentially, and wherein this
evaluation can itself change the state space and/or fitness function, because
no two variants are ever evaluated according to the same selection criterion.
Concerns are also raised as to the methodology used in Simonton's study.


  Theoretical calculations of solvation contribution to hole energy in a
polynucleotide chain give very low hole mobility values at zero temperature,
\mu < 10^{-3} cm^2/(V s). We calculated hole mobility at physiological
temperature for the Poly G/Poly C DNA duplex, which gave substantially larger
mobility values. Mobility over the temperature range 20-400 K was calculated.
Taking stacking interaction into account substantially increased hole mobility.


  We investigate the dynamics of cancer initiation in a mathematical model with
one driver mutation and several passenger mutations. Our analysis is based on a
multi type branching process: We model individual cells which can either divide
or undergo apoptosis. In case of a cell division, the two daughter cells can
mutate, which potentially confers a change in fitness to the cell. In contrast
to previous models, the change in fitness induced by the driver mutation
depends on the genetic context of the cell, in our case on the number of
passenger mutations. The passenger mutations themselves have no or only a very
small impact on the cell's fitness. While our model is not designed as a
specific model for a particular cancer, the underlying idea is motivated by
clinical and experimental observations in Burkitt Lymphoma. In this tumor, the
hallmark mutation leads to deregulation of the MYC oncogene which increases the
rate of apoptosis, but also the proliferation rate of cells. This increase in
the rate of apoptosis hence needs to be overcome by mutations affecting
apoptotic pathways, naturally leading to an epistatic fitness landscape. This
model shows a very interesting dynamical behavior which is distinct from the
dynamics of cancer initiation in the absence of epistasis. Since the driver
mutation is deleterious to a cell with only a few passenger mutations, there is
a period of stasis in the number of cells until a clone of cells with enough
passenger mutations emerges. Only when the driver mutation occurs in one of
those cells, the cell population starts to grow rapidly.


  The Spanish government recently announced an official fast-track path to
citizenship for any individual who is Jewish and whose ancestors were expelled
from Spain during the inquisition-related dislocation of Spanish Jews in 1492.
It would seem that this policy targets a small subset of the global Jewish
population, i.e., restricted to individuals who retain cultural practices
associated with ancestral origins in Spain. However, the central contribution
of this manuscript is to demonstrate how and why the policy is far more likely
to apply to a very large fraction (i.e., the vast majority) of Jews. This claim
is supported using a series of genealogical models that include transmissable
"identities" and preferential intra-group mating. Model analysis reveals that
even when intra-group mating is strong and even if only a small subset of a
present-day population retains cultural practices typically associated with
that of an ancestral group, it is highly likely that nearly all members of that
population have direct geneaological links to that ancestral group, given
sufficient number of generations have elapsed. The basis for this conclusion is
that not having a link to an ancestral group must be a property of all of an
individual's ancestors, the probability of which declines (nearly)
superexponentially with each successive generation. These findings highlight
unexpected incongruities induced by genealogical dynamics between present-day
and ancestral identities.


  Many real networks exhibit a layered structure in which links in each layer
reflect the function of nodes on different environments. These multiple types
of links are usually represented by a multiplex network in which each layer has
a different topology. In real-world networks, however, not all nodes are
present on every layer. To generate a more realistic scenario, we use a
generalized multiplex network and assume that only a fraction $q$ of the nodes
are shared by the layers. We develop a theoretical framework for a branching
process to describe the spread of an epidemic on these partially overlapped
multiplex networks. This allows us to obtain the fraction of infected
individuals as a function of the effective probability that the disease will be
transmitted $T$. We also theoretically determine the dependence of the epidemic
threshold on the fraction $q > 0$ of shared nodes in a system composed of two
layers. We find that in the limit of $q \to 0$ the threshold is dominated by
the layer with the smaller isolated threshold. Although a system of two
completely isolated networks is nearly indistinguishable from a system of two
networks that share just a few nodes, we find that the presence of these few
shared nodes causes the epidemic threshold of the isolated network with the
lower propagating capacity to change discontinuously and to acquire the
threshold of the other network.


  Dynamic Monte Carlo simulations are used to study coupled transport
(co-transport) through sub-nanometer-diameter pores. In this classic
Hodgkin-Keynes mechanism, an ion species uses the large flux of an abundant ion
species to move against its concentration gradient. The efficiency of
co-transport is examined for various pore parameters so that synthetic
nanopores can be engineered to maximize this effect. In general, the pore must
be narrow enough that ions cannot pass each other and the charge of the pore
large enough to attract many ions so that they exchange momentum. Co-transport
efficiency increases as pore length increases, but even very short pores
exhibit co-transport, in contradiction to the usual perception that long pores
are necessary. The parameter ranges where co-transport occurs is consistent
with current and near-future synthetic nanopore geometry parameters, suggesting
that co-transport of ions may be a new application of nanopores.


  Purpose and meaning are necessary concepts for understanding mind and
culture, but appear to be absent from the physical world and are not part of
the explanatory framework of the natural sciences. Understanding how meaning
(in the broad sense of the term) could arise from a physical world has proven
to be a tough problem. The basic scheme of Darwinian evolution produces
adaptations that only represent apparent ("as if") goals and meaning. Here I
use evolutionary models to show that a slight, evolvable extension of the basic
scheme is sufficient to produce genuine goals. The extension, targeted
modulation of mutation rate, is known to be generally present in biological
cells, and gives rise to two phenomena that are absent from the non-living
world: intrinsic meaning and the ability to initiate goal-directed chains of
causation (active causation). The extended scheme accomplishes this by
utilizing randomness modulated by a feedback loop that is itself regulated by
evolutionary pressure. The mechanism can be extended to behavioural variability
as well, and thus shows how freedom of behaviour is possible. A further
extension to communication suggests that the active exchange of intrinsic
meaning between organisms may be the origin of consciousness, which in
combination with active causation can provide a physical basis for the
phenomenon of free will.


  Biological reaction networks often contain what might be called 'hub
molecules', which are involved in many reactions. For example, ATP is commonly
consumed and produced. When reaction networks contain molecules like ATP, they
are difficult to efficiently simulate, because every time such a molecule is
consumed or produced, the propensities of numerous reactions need to be
updated. In order to increase the speed of simulations, we developed 'Lazy
Updating', which postpones some propensity updates until some aspect of the
state of the system changes by more than a defined threshold. Lazy Updating
works with several existing stochastic simulation algorithms, including
Gillespie's direct method and the Next Reaction Method. We tested Lazy Updating
on two example models, and for the larger model it increased the speed of
simulations over eight-fold while maintaining a high level of accuracy. These
increases in speed will be larger for models with more widely connected hub
molecules. Thus Lazy Updating can contribute towards making models with a
limited computing time budget more realistic by including previously neglected
hub molecules.


  Solutes added to solutions often dramatically impact molecular processes
ranging from the suspension or precipitation of colloids to biomolecular
associations and protein folding. Here we revisit the origins of the effective
attractive interactions that emerge between and within macromolecules immersed
in solutions containing cosolutes that are preferentially excluded from the
macromolecular interfaces. Until recently, these depletion forces were
considered to be entropic in nature, resulting primarily from the tendency to
increase the space available to the cosolute. However, recent experimental
evidence indicates the existence of additional, energetically-dominated
mechanisms. In this review we follow the emerging characteristics of these
different mechanisms. By compiling a set of available thermodynamic data for
processes ranging from protein folding to protein-protein interactions, we show
that excluded cosolutes can act through two distinct mechanisms that correlate
to a large extent with their molecular properties. For many polymers at low to
moderate concentrations the steric interactions and molecular crowding effects
dominate, and the mechanism is entropic. To contrast, for many small excluded
solutes, such as naturally occurring osmolytes, the mechanism is dominated by
favorable enthalpy, whereas the entropic contribution is typically unfavorable.
We review the available models for these thermodynamic mechanisms, and comment
on the need for new models that would be able to explain the full range of
observed depletion forces.


  Ecology studies biodiversity in its variety and complexity. It describes how
species distribute and perform in response to environmental changes. Ecological
processes and structures are highly complex and adaptive. In order to quantify
emerging ecological patterns and investigate their hidden mechanisms, we need
to rely on the simplicity of mathematical language. This becomes especially
apparent when dealing with scaling patterns in ecology. Indeed, nearly all of
ecological patterns are scale dependent. Such scale dependence hampers our
predictive power and creates problems in our inference. This challenge calls
for a clear and fundamental understanding of how and why ecological patterns
change across scales. As Simon Levin stated in his MacArthur Award lecture, the
problem of relating phenomena across scales is the central problem in ecology
and other natural sciences. It has become clear that there is currently a drive
in ecology and complexity science to develop new quantitative approaches that
are suitable for analysing and forecasting patterns of ecological systems. Here
I provide a road map for future works on synthesizing the scaling patterns in
ecology, aiming (i) to collect and sort a diverse array of ecological patterns,
(ii) to present the dominant parametric forms of how these patterns change
across spatial and temporal scales, (iii) to detect the processes and
mechanisms using mathematical models, and finally (iv) to probe the physical
meaning of these scaling patterns. This road map is divided into three parts
and covers three main concepts of scale in ecology: heterogeneity, hierarchy
and size. Using scale as a thread, this road map and its following works weave
the kaleidoscope of ecological scaling patterns into a cohesive whole.


  Most of major algorithms for phylogenetic tree reconstruction assume that
sequences in the analyzed set either do not have any offspring, or that parent
sequences can maximally mutate into just two descendants. The graph resulting
from such assumptions forms therefore a binary tree, with all the nodes labeled
as leaves. However, these constraints are unduly restrictive as there are
numerous data sets with multiple offspring of the same ancestors. Here we
propose a solution to analyze and visualize such sets in a more intuitive
manner. The method reconstructs phylogenetic tree by assigning the sequences
with offspring as internal nodes, and the sequences without offspring as leaf
nodes. In the resulting tree there is no constraint for the number of adjacent
nodes, which means that the solution tree needs not to be a binary graph only.
The subsequent derivation of evolutionary pathways, and pair-wise mutations,
are then an algorithmically straightforward, with edge's length corresponding
directly to the number of mutations. Other tree reconstruction algorithms can
be extended in the proposed manner, to also give unbiased topologies.


  Motivation: Microarray results accumulated in public repositories are widely
re-used in meta-analytical studies and secondary databases. The quality of the
data obtained with this technology varies from experiment to experiment and
efficient method for quality assessment is neces-sary to ensure their
reliability. Results: The lack of a good benchmark has hampered evaluation of
existing methods for quality control. In this study we propose a new
inde-pendent quality metric that is based on evolutionary conservation of
expression profiles. We show, using 11 large organ-specific datasets, that
IQRray, a new quality metrics developed by us, exhibits the highest correlation
with this reference metric, among 14 metrics tested. IQRray outperforms other
methods in identification of poor quality arrays in dataset composed of arrays
from many independent experiments. In con-trast, the performance of methods
designed for detecting outliers in a single experiment like NUSE and RLE was
low because of the inability of these method to detect datasets containing only
low quality arrays, and the fact that the scores cannot be directly compared
between ex-periments. Availability: The R implementation of IQRray is available
at: ftp://lausanne.isb-sib.ch/pub/databases/Bgee/general/IQRray.R


  We show that the signal-processing paradigm known as compressed sensing (CS)
is applicable to genome-wide association studies (GWAS) and genomic selection
(GS). The aim of GWAS is to isolate trait-associated loci, whereas GS attempts
to predict the phenotypic values of new individuals on the basis of training
data. CS addresses a problem common to both endeavors, namely that the number
of genotyped markers often greatly exceeds the sample size. We show using CS
methods and theory that all loci of nonzero effect can be identified (selected)
using an efficient algorithm, provided that they are sufficiently few in number
(sparse) relative to sample size. For heritability h2 = 1, there is a sharp
phase transition to complete selection as the sample size is increased. For
heritability values less than one, complete selection can still occur although
the transition is smoothed. The transition boundary is only weakly dependent on
the total number of genotyped markers. The crossing of a transition boundary
provides an objective means to determine when true effects are being recovered;
we discuss practical methods for detecting the boundary. For h2 = 0.5, we find
that a sample size that is thirty times the number of nonzero loci is
sufficient for good recovery.


  Detecting communities, densely connected groups may contribute to unravel the
underlying relationships among the units present in diverse biological networks
(e.g., interactome, coexpression networks, ecological networks, etc.). We
recently showed that communities can be very precisely characterized by
maximizing Surprise, a global network parameter. Here we present SurpriseMe, a
tool that integrates the outputs of seven of the best algorithms available to
estimate the maximum Surprise value. SurpriseMe also generates distance
matrices that allow to visualize the relationships among the solutions
generated by the algorithms. We show that the communities present in small and
medium-sized networks, with up to 10.000 nodes, can be easily characterized: on
standard PC computers, these analyses take less than an hour. Also, four of the
algorithms may quite rapidly analyze networks with up to 100.000 nodes, given
enough memory resources. Because of its performance and simplicity, SurpriseMe
is a reference tool for community structure characterization.


  In the presence of ATP, molecular motors generate active force dipoles that
drive suspensions of protein filaments far from thermodynamic equilibrium,
leading to exotic dynamics and pattern formation. Microscopic modelling can
help to quantify the relationship between individual motors plus filaments to
organisation and dynamics on molecular and supra-molecular length scales. Here
we present results of extensive numerical simulations of active gels where the
motors and filaments are confined between two infinite parallel plates. Thermal
fluctuations and excluded-volume interactions between filaments are included. A
systematic variation of rates for motor motion, attachment and detachment,
including a differential detachment rate from filament ends, reveals a range of
non-equilibrium behaviour. Strong motor binding produces structured filament
aggregates that we refer to as asters, bundles or layers, whose stability
depends on motor speed and differential end-detachment. The gross features of
the dependence of the observed structures on the motor rate and the filament
concentration can be captured by a simple one-filament model. Loosely bound
aggregates exhibit super-diffusive mass transport, where filament translocation
scales with lag time with non-unique exponents that depend on motor kinetics.
An empirical data collapse of filament speed as a function of motor speed and
end-detachment is found, suggesting a dimensional reduction of the relevant
parameter space. We conclude by discussing the perspectives of microscopic
modelling in the field of active gels.


  Protein function depends on both protein structure and amino acid (aa)
sequence. Here we show that modular features of both structure and function can
be quantified from the aa sequences alone for the small (40,42 aa)
plaque-forming amyloid beta fragments. Some edge and center features of the
fragments are predicted. Contrasting results from the second order
hydropathicity scale based on evolutionary optimization (self-organized
criticality) and the first order scale based on complete protein (water-air)
unfolding show that fragmentation has mixed first- and second-order character.


  Control of complex processes is a major goal of network analyses. Most
approaches to control nonlinearly coupled systems require the network topology
and/or network dynamics. Unfortunately, neither the full set of participating
nodes nor the network topology is known for many important systems. On the
other hand, system responses to perturbations are often easily measured. We
show how the collection of such responses (a response surface) can be used for
network control. Analyses of model systems show that response surfaces are
smooth and hence can be approximated using low order polynomials. Importantly,
these approximations are largely insensitive to stochastic fluctuations in data
or measurement errors. They can be used to compute how a small set of nodes
need to be altered in order to direct the network close to a pre-specified
target state. These ideas, illustrated on a nonlinear electrical circuit, can
prove useful in many contexts including in reprogramming cellular states.


  Time-calibrated species phylogenies are critical for addressing a wide range
of questions in evolutionary biology, such as those that elucidate historical
biogeography or uncover patterns of coevolution and diversification. Because
molecular sequence data are not informative on absolute time, external data,
most commonly fossil age estimates, are required to calibrate estimates of
species divergence dates. For Bayesian divergence-time methods, the common
practice for calibration using fossil information involves placing arbitrarily
chosen parametric distributions on internal nodes, often disregarding most of
the information in the fossil record. We introduce the 'fossilized birth-death'
(FBD) process, a model for calibrating divergence-time estimates in a Bayesian
framework, explicitly acknowledging that extant species and fossils are part of
the same macroevolutionary process. Under this model, absolute node age
estimates are calibrated by a single diversification model and arbitrary
calibration densities are not necessary. Moreover, the FBD model allows for
inclusion of all available fossils. We performed analyses of simulated data and
show that node-age estimation under the FBD model results in robust and
accurate estimates of species divergence times with realistic measures of
statistical uncertainty, overcoming major limitations of standard divergence
time estimation methods. We then used this model to estimate the speciation
times for a dataset composed of all living bears, indicating that the genus
Ursus diversified in the late Miocene to mid Pliocene.


  We study the dependence of the largest component in regular networks on the
clustering coefficient, showing that its size changes smoothly without
undergoing a phase transition. We explain this behaviour via an analytical
approach based on the network structure, and provide an exact equation
describing the numerical results. Our work indicates that intrinsic structural
properties always allow the spread of epidemics on regular networks.


  Enzymes within biochemical pathways are often colocalized, yet the
consequences of specific spatial enzyme arrangements remain poorly understood.
We study the impact of enzyme arrangement on reaction efficiency within a
reaction-diffusion model. The optimal arrangement transitions from a cluster to
a distributed profile as a single parameter, which controls the probability of
reaction versus diffusive loss of pathway intermediates, is varied. We
introduce the concept of enzyme exposure to explain how this transition arises
from the stochastic nature of molecular reactions and diffusion.


  The spatial organization of enzymes often plays a crucial role in the
functionality and efficiency of enzymatic pathways. To fully understand the
design and operation of enzymatic pathways, it is therefore crucial to
understand how the relative arrangement of enzymes affects pathway function.
Here we investigate the effect of enzyme localization on the flux of a minimal
two-enzyme pathway within a reaction-diffusion model. We consider different
reaction kinetics, spatial dimensions, and loss mechanisms for intermediate
substrate molecules. Our systematic analysis of the different regimes of this
model reveals both universal features and distinct characteristics in the
phenomenology of these different systems. In particular, the distribution of
the second pathway enzyme that maximizes the reaction flux undergoes a generic
transition from co-localization with the first enzyme when the catalytic
efficiency of the second enzyme is low, to an extended profile when the
catalytic efficiency is high. However, the critical transition point and the
shape of the extended optimal profile is significantly affected by specific
features of the model. We explain the behavior of these different systems in
terms of the underlying stochastic reaction and diffusion processes of single
substrate molecules.


  We are entering the era of ubiquitous genetic information for research,
clinical care, and personal curiosity. Sharing these datasets is vital for
rapid progress in understanding the genetic basis of human diseases. However,
one growing concern is the ability to protect the genetic privacy of the data
originators. Here, we technically map threats to genetic privacy and discuss
potential mitigation strategies for privacy-preserving dissemination of genetic
data.


  forqs is a forward-in-time simulation of recombination, quantitative traits,
and selection. It was designed to investigate haplotype patterns resulting from
scenarios where substantial evolutionary change has taken place in a small
number of generations due to recombination and/or selection on polygenic
quantitative traits. forqs is implemented as a command- line C++ program.
Source code and binary executables for Linux, OSX, and Windows are freely
available under a permissive BSD license.


  The detection of genomic regions unusually rich in a given pattern is an
important undertaking in the analysis of next generation sequencing data.
Recent studies of chromosomal translocations in activated B lymphocytes have
identified regions that are frequently translocated to c-myc oncogene. A
quantitative method for the identification of translocation hotspots was
crucial to this study. Here we improve this analysis by using a simple
probabilistic model and the framework provided by scan statistics to define the
number and location of translocation breakpoint hotspots. A key feature of our
method is that it provides a global chromosome-wide significance level to
clustering, as opposed to previous methods based on local criteria. Whilst
being motivated by a specific application, the detection of unusual clusters is
a widespread problem in bioinformatics. We expect our method to be useful in
the analysis of data from other experimental approaches such as of ChIP-seq and
4C-seq.