  The interrelationships of the fundamental biological processes natural
selection, mutation, and stochastic drift are quantified by the entropy rate of
Moran processes with mutation, measuring the long-run variation of a Markov
process. The entropy rate is shown to behave intuitively with respect to
evolutionary parameters such as monotonicity with respect to mutation
probability (for the neutral landscape), relative fitness, and strength of
selection. Strict upper bounds, depending only on the number of replicating
types, for the entropy rate are given and the neutral fitness landscape attains
the maximum in the large population limit. Various additional limits are
computed including small mutation, weak and strong selection, and large
population holding the other parameters constant, revealing the individual
contributions and dependences of each evolutionary parameter on the long-run
outcomes of the processes.


  Cells often exhibit different and stable phenotypes from the same DNA
sequence. Robustness and plasticity of such cellular states are controlled by
diverse transcriptional and epigenetic mechanisms, among them the modification
of biochemical marks on chromatin. Here, we develop a stochastic model that
describes the dynamics of epigenetic marks along a given DNA region. Through
mathematical analysis, we show the emergence of bistable and persistent
epigenetic states from the cooperative recruitment of modifying enzymes. We
also find that the dynamical system exhibits a critical point and displays, in
presence of asymmetries in recruitment, a bifurcation diagram with hysteresis.
These results have deep implications for our understanding of epigenetic
regulation. In particular, our study allows to reconcile within the same
formalism the robust maintenance of epigenetic identity observed in
differentiated cells, the epigenetic plasticity of pluripotent cells during
differentiation and the effects of epigenetic misregulation in diseases.
Moreover, it suggests a possible mechanism for developmental transitions where
the system is shifted close to the critical point to benefit from high
susceptibility to developmental cues.


  More than twenty coarse-grained (CG) DNA models have been developed for
simulating the behavior of this molecule under various conditions, including
those required for nanotechnology. However, none of these models reproduces the
DNA polymorphism associated with conformational changes in the ribose rings of
the DNA backbone. These changes make an essential contribution to the DNA local
deformability and provide the possibility of the transition of the DNA double
helix from the B-form to the A-form during interactions with biological
molecules. We propose a CG representation of the ribose conformational
flexibility. We substantiate the choice of the CG sites (6 per nucleotide)
needed for the "sugar" GC DNA model, and obtain the potentials of the CG
interactions between the sites by the "bottom-up" approach using the all-atom
AMBER force field. We show that the representation of the ribose flexibility
requires one non-harmonic and one three-particle potential, the forms of both
the potentials being different from the ones generally used. The model also
includes (i) explicit representation of ions (in an implicit solvent) and (ii)
sequence dependence. With these features, the sugar CG DNA model reproduces
(with the same parameters) both the B- and A- stable forms under corresponding
conditions and demonstrates both the A to B and the B to A phase transitions.


  We examine the modular structure of the metabolic network when combined with
the regulatory network representing direct regulation of enzymes by small
metabolites in E.coli. In order to identify the modular structure we introduce
clustering algorithm based on a novel vertex similarity measure for bipartite
graphs. We also apply a standard module identification method based on
simulated annealing. Both methods identify the same modular core each of them
with different resolution. We observe slight but still statistically
significant increase of modularity after regulatory interactions addition.
Enrichment of the metabolic network with the regulatory information leads to
identification of new functional modules, which cannot be detected in the
metabolic network only. Regulatory loops in the modules are the source of their
self-control, i.e. autonomy, and allow to make hypothesis about module
function. This study demonstrates that incorporation of regulatory information
is important component in defining functional units of the metabolic network.


  The capacity for aerial maneuvering shaped the evolution of flying animals.
Here we evaluate consequences of aviaian morphology for aerial performance
(1,2) by quantifying static stability and control effectiveness of physical
models (3) for numerous taxa sampled from within the lineage leading to birds
(Paraves, 4). Results of aerodynamic testing are mapped phylogenetically (5-9)
to examine how maneuvering characteristics correlate with tail shortening,
fore- and hindwing elaboration, and other morphological features (10). In the
evolution of the Avialae we observe shifts from static stability to inherently
unstable aerial planforms; control effectiveness also migrated from tails to
the forewings. These shifts suggest that some degree of aerodynamic control and
and capacity for maneuvering preceded the evolution of strong power stroke. The
timing of shifts also suggests some features normally considered in light of
development of a power stroke may play important roles in control.


  Chargaff's second parity rule for short oligonucleotides states that the
frequency of any short nucleotide sequence on a strand is approximately equal
to the frequency of its reverse complement on the same strand. Recent studies
have shown that, with the exception of organellar DNA, this parity rule
generally holds for double stranded DNA genomes and fails to hold for
single-stranded genomes. While Chargaff's first parity rule is fully explained
by the Watson-Crick pairing in the DNA double helix, a definitive explanation
for the second parity rule has not yet been determined. In this work, we
propose a model based on a hidden Markov process for approximating the
distributional structure of primitive DNA sequences. Then, we use the model to
provide another possible theoretical explanation for Chargaff's second parity
rule, and to predict novel distributional aspects of bacterial DNA sequences.


  Most sensory cells use surface receptors to detect environmental stimuli and
initiate downstream signaling. Cooperative interactions among sensory receptors
is known to play a crucial role in enhancing the sensitivity of biochemical
processes such as oxygen sensing by hemoglobin, but whether cooperativity
enhances the fidelity with which a system with multiple receptors can
accurately and quickly detect a signal is poorly understood. We model the
kinetics of small clusters of receptors in the presence of ligand, where the
receptors act independently or cooperatively. We show that the interaction
strength and how it is coupled to the dynamics influences the macroscopic
observables. Contrary to recent reports, our analysis shows that receptor
cooperativity can increase the signal-to-noise ratio, but this increase depends
on the underlying dynamics of the signaling receptor cluster.


  As noninvasive sampling techniques for animal populations have become more
popular, there has been increasing interest in the development of
capture-recapture models that can accommodate both imperfect detection and
misidentification of individuals (e.g., due to genotyping error). However,
current methods do not allow for individual variation in parameters, such as
detection or survival probability. Here we develop misidentification models for
capture-recapture data that can simultaneously account for temporal variation,
behavioral effects and individual heterogeneity in parameters. To facilitate
Bayesian inference using our approach, we extend standard probit regression
techniques to latent multinomial models where the dimension and zeros of the
response cannot be observed. We also present a novel Metropolis-Hastings within
Gibbs algorithm for fitting these models using Markov chain Monte Carlo. Using
closed population abundance models for illustration, we re-visit a DNA
capture-recapture population study of black bears in Michigan, USA and find
evidence of misidentification due to genotyping error, as well as temporal,
behavioral and individual variation in detection probability. We also estimate
a salamander population of known size from laboratory experiments evaluating
the effectiveness of a marking technique commonly used for amphibians and fish.
Our model was able to reliably estimate the size of this population and
provided evidence of individual heterogeneity in misidentification probability
that is attributable to variable mark quality. Our approach is more
computationally demanding than previously proposed methods, but it provides the
flexibility necessary for a much broader suite of models to be explored while
properly accounting for uncertainty introduced by misidentification and
imperfect detection. In the absence of misidentification, our probit
formulation also provides a convenient and efficient Gibbs sampler for Bayesian
analysis of traditional closed population capture-recapture data.


  We study a simple model of a forager as a walk that modifies a relaxing
substrate. Within it simplicity, this provides an insight on a number of
relevant and non-intuitive facts. Even without memory of the good places to
feed and no explicit cost of moving, we observe the emergence of a finite home
range. We characterize the walks and the use of resources in several
statistical ways, involving the behavior of the average used fraction of the
system, the length of the cycles followed by the walkers, and the frequency of
visits to plants. Preliminary results on population effects are explored by
means of a system of two non directly interacting animals. Properties of the
overlap of home ranges show the existence of a set of parameters that provides
the best utilization of the shared resource.


  Whole exome sequencing was performed on HLA-matched stem cell donors and
transplant recipients to measure sequence variation contributing to minor
histocompatibility antigen differences between the two. A large number of
nonsynonymous single nucleotide polymorphisms were identified in each of the
nine unique donor-recipient pairs tested. This variation was greater in
magnitude in unrelated donors as compared with matched related donors.
Knowledge of the magnitude of exome variation between stem cell transplant
recipients and donors may allow more accurate titration of immunosuppressive
therapy following stem cell transplantation.


  Complex networks are intrinsically modular. Resolving small modules is
particularly difficult when the network is densely connected; wide variation of
link weights invites additional complexities. In this article we present an
algorithm to detect community structure in densely connected weighted networks.
First, modularity of the network is calculated by erasing the links having
weights smaller than a cutoff $q.$ Then one takes all the disjoint components
obtained at $q=q_c,$ where the modularity is maximum, and modularize the
components individually using Newman Girvan's algorithm for weighted networks.
We show, taking microRNA (miRNA) co-target network of Homo sapiens as an
example, that this algorithm could reveal miRNA modules which are known to be
relevant in biological context.


  We survey and introduce concepts and tools located at the intersection of
information theory and network biology. We show that Shannon's information
entropy, compressibility and algorithmic complexity quantify different local
and global aspects of synthetic and biological data. We show examples such as
the emergence of giant components in Erdos-Renyi random graphs, and the
recovery of topological properties from numerical kinetic properties simulating
gene expression data. We provide exact theoretical calculations, numerical
approximations and error estimations of entropy, algorithmic probability and
Kolmogorov complexity for different types of graphs, characterizing their
variant and invariant properties. We introduce formal definitions of complexity
for both labeled and unlabeled graphs and prove that the Kolmogorov complexity
of a labeled graph is a good approximation of its unlabeled Kolmogorov
complexity and thus a robust definition of graph complexity.


  In a randomly-mating biparental population of size $N$ there are, with high
probability, individuals who are genealogical ancestors of every extant
individual within approximately $\log_2(N)$ generations into the past. We use
this result of J. Chang to prove a curious corollary under standard models of
recombination: there exist, with high probability, individuals within a
constant multiple of $ \log_2(N)$ generations into the past who are
simultaneously (i) genealogical ancestors of {\em each} of the individuals at
the present, and (ii) genetic ancestors to {\em none} of the individuals at the
present. Such ancestral individuals - ancestors of everyone today that left no
genetic trace -- represent `ghost' ancestors in a strong sense. In this short
note, we use simple analytical argument and simulations to estimate how many
such individuals exist in finite Wright-Fisher populations.


  fwdpp is a C++ library of routines intended to facilitate the development of
forward-time simulations under arbitrary mutation and fitness models. The
library design provides a combination of speed, low memory overhead, and
modeling flexibility not currently available from other forward simulation
tools. The library is particularly useful when the simulation of large
populations is required, as programs implemented using the library are much
more efficient that other available forward simulation programs.


  In a recent work [arXiv:1307.1382] it was pointed out that the link-weight
distribution of microRNA (miRNA) co-target network of a wide class of species
are universal up to scaling. The number cell types, widely accepted as a
measure of complexity, turns out to be proportional to these scale-factor. In
this article we discuss additional universal features of these networks and
show that, this universality splits if one considers distribution of number of
common targets of three or more number of miRNAs. These distributions for
different species can be collapsed onto two distinct set of universal
functions, revealing the fact that the species which appeared in early
evolution have different complexity measure compared to those appeared late.


  We describe a compression-based distance for genomic sequences. Instead of
using the usual conjoint information content, as in the classical Normalized
Compression Distance (NCD), it uses the conditional information content. To
compute this Normalized Conditional Compression Distance (NCCD), we need a
normal conditional compressor, that we built using a mixture of static and
dynamic finite-context models. Using this approach, we measured chromosomal
distances between Hominidae primates and also between Muroidea (rat and mouse),
observing several insights of evolution that so far have not been reported in
the literature.


  Large whole-genome sequencing projects have provided access to much of the
rare variation in human populations, which is highly informative about
population structure and recent demography. Here, we show how the age of rare
variants can be estimated from patterns of haplotype sharing and how these ages
can be related to historical relationships between populations. We investigate
the distribution of the age of variants occurring exactly twice (f2 variants)
in a worldwide sample sequenced by the 1000 Genomes Project, revealing enormous
variation across populations. The median age of haplotypes carrying f2 variants
is 50 to 160 generations across populations within Europe or Asia, and 170 to
320 generations within Africa. Haplotypes shared between continents are much
older with median ages for haplotypes shared between Europe and Asia ranging
from 320 to 670 generations. The distribution of the ages of f2 haplotypes is
informative about their demography, revealing recent bottlenecks, ancient
splits, and more modern connections between populations. We see the signature
of selection in the observation that functional variants are significantly
younger than nonfunctional variants of the same frequency. This approach is
relatively insensitive to mutation rate and complements other nonparametric
methods for demographic inference.


  The Internet has enabled the emergence of collective problem solving, also
known as crowdsourcing, as a viable option for solving complex tasks. However,
the openness of crowdsourcing presents a challenge because solutions obtained
by it can be sabotaged, stolen, and manipulated at a low cost for the attacker.
We extend a previously proposed crowdsourcing dilemma game to an iterated game
to address this question. We enumerate pure evolutionarily stable strategies
within the class of so-called reactive strategies, i.e., those depending on the
last action of the opponent. Among the 4096 possible reactive strategies, we
find 16 strategies each of which is stable in some parameter regions. Repeated
encounters of the players can improve social welfare when the damage inflicted
by an attack and the cost of attack are both small. Under the current
framework, repeated interactions do not really ameliorate the crowdsourcing
dilemma in a majority of the parameter space.


  In this work we describe a non-parametric disease model that links the
temporal change of the prevalence of an infectious disease to the incidence and
the recovery rates. The model is only based on the common epidemiological
measures incidence and recovery rate. As an application, the model is used to
calculate the prevalence of influenza in Germany for a hypothetical birth
cohort during 2001 and 2013.


  In the framework of a 1D cyclic competition model, the Rock-Paper-Scissor
model, where bacteria are allowed to mutate and move in space, we study the
formation of stochastic patterns, where all the bacteria species do coexist. We
modelled the problem using an individual-based setting and using the system
size van Kampen expansion to deal with the Master Equation, we have been able
to characterise the spatio-temporal patterns using the power spectrum of the
fluctuations. We proved that such patterns are robust against the intrinsic
noise and they can be found for parameters values beyond the ones fixed by the
deterministic approach. We complement such analytical results with numerical
simulations based on the Gillespie's algorithm.


  In this paper the problem of surface charge of the lipid membrane immersed in
the physiological solution is considered. It is shown that both side of the
bilayer phospholipid membrane surface are negatively charged. A self-consistent
model of the potential in solution is developed, and a stationary charge
density on the membrane surface is found. It is shown that the ions of the
surface charge are in a relatively deep (as compared to kBT) potential wells,
which are localized near the dipole heads of phospholipid membrane. It makes
impossible for ions to slip along the membrane surface. Simple experiments for
verifying the correctness of the considered model are proposed. A developed
approach can be used for estimations of the surface charges on the outer and
inner membrane of the cell.


  Finite-context modeling is a powerful tool for compressing and hence for
representing DNA sequences. We describe an algorithm to detect genomic
regularities, within a blind discovery strategy. The algorithm uses information
profiles built using suitable combinations of finite-context models. We used
the genome of the fission yeast Schizosaccharomyces pombe strain 972 h- for
illustration, unveilling locations of low information content, which are
usually associated with DNA regions of potential biological interest.


  Networks are widely used to represent interaction pattern among the
components in complex systems. Structures of real networks from differ- ent
domains may vary quite significantly. Since there is an interplay be- tween
network architecture and dynamics, structure plays an important role in
communication and information spreading on a network. Here we investigate the
underlying undirected topology of different biological networks which support
faster spreading of information and are better in communication. We analyze the
good expansion property by using the spectral gap and communicability between
nodes. Different epidemic models are also used to study the transmission of
information in terms of disease spreading through individuals (nodes) in those
networks. More- over, we explore the structural conformation and properties
which may be responsible for better communication. Among all biological
networks studied here, the undirected structure of neuronal networks not only
pos- sesses the small-world property but the same is expressed remarkably to a
higher degree than any randomly generated network which possesses the same
degree sequence. A relatively high percentage of nodes, in neuronal networks,
form a higher core in their structure. Our study shows that the underlying
undirected topology in neuronal networks is significantly qualitatively
different than the same from other biological networks and that they may have
evolved in such a way that they inherit a (undirected) structure which is
excellent and robust in communication.


  Orphans are genes restricted to a single phylogenetic lineage and emerge at
high rates. While this predicts an accumulation of genes, the gene number has
remained remarkably constant through evolution. This paradox has not yet been
resolved. Because orphan genes have been mainly analyzed over long evolutionary
time scales, orphan loss has remained unexplored. Here we study the patterns of
orphan turnover among close relatives in the Drosophila obscura group. We show
that orphans are not only emerging at a high rate, but that they are also
rapidly lost. Interestingly, recently emerged orphans are more likely to be
lost than older ones. Furthermore, highly expressed orphans with a strong
male-bias are more likely to be retained. Since both lost and retained orphans
show similar evolutionary signatures of functional conservation, we propose
that orphan loss is not driven by high rates of sequence evolution, but
reflects lineage specific functional requirements.


  The complex interplay between population movements in space and
non-homogeneous mixing patterns have so far hindered the fundamental
understanding of the conditions for spatial invasion through a general
theoretical framework. To address this issue, we present an analytical
modelling approach taking into account such interplay under general conditions
of mobility and interactions, in the simplifying assumption of two population
classes. We describe a spatially structured population with non-homogeneous
mixing and travel behaviour through a multi-host stochastic epidemic
metapopulation model. Different population partitions, mixing patterns and
mobility structures are considered, along with a specific application for the
study of the role of age partition in the early spread of the 2009 H1N1
pandemic influenza.
  We provide a complete mathematical formulation of the model and derive a
semi-analytical expression of the threshold condition for global invasion of an
emerging infectious disease in the metapopulation system. A rich solution space
is found that depends on the social partition of the population, the pattern of
contacts across groups and their relative social activity, the travel attitude
of each class, and the topological and traffic features of the mobility
network. Reducing the activity of the less social group and reducing the
cross-group mixing are predicted to be the most efficient strategies for
controlling the pandemic potential in the case the less active group
constitutes the majority of travellers. If instead traveling is dominated by
the more social class, our model predicts the existence of an optimal
across-groups mixing that maximises the pandemic potential of the disease,
whereas the impact of variations in the activity of each group is less
important.


  We calculate measures of hierarchy in gene and tissue networks of breast
cancer patients. We find that the likelihood of metastasis in the future is
correlated with increased values of network hierarchy for expression networks
of cancer-associated genes, due to correlated expression of cancer-specific
pathways. Conversely, future metastasis and quick relapse times are negatively
correlated with values of network hierarchy in the expression network of all
genes, due to dedifferentiation of gene pathways and circuits. These results
suggest that hierarchy of gene expression may be useful as an additional
biomarker for breast cancer prognosis.


  Bacteria and archaea have evolved an adaptive, heritable immune system that
recognizes and protects against viruses or plasmids. This system, known as the
CRISPR-Cas system, allows the host to recognize and incorporate short foreign
DNA or RNA sequences, called `spacers' into its CRISPR system. Spacers in the
CRISPR system provide a record of the history of bacteria and phage
coevolution. We use a physical model to study the dynamics of this coevolution
as it evolves stochastically over time. We focus on the impact of mutation and
recombination on bacteria and phage evolution and evasion. We discuss the
effect of different spacer deletion mechanisms on the coevolutionary dynamics.
We make predictions about bacteria and phage population growth, spacer
diversity within the CRISPR locus, and spacer protection against the phage
population.


  Motivation: Open-source bacterial genome assembly remains inaccessible to
many biologists due to its complexity. Few software solutions exist that are
capable of automating all steps in the process of de novo genome assembly from
Illumina data.
  Results: A5-miseq can produce high quality microbial genome assemblies from
as little as 20-fold sequence data coverage on a laptop computer without any
parameter tuning. A5-miseq does this by automating the process of adapter
trimming, quality filtering, error correction, contig and scaffold generation,
and detection of misassemblies. Unlike the original A5 pipeline, A5-miseq can
use long reads from the Illumina MiSeq, use read pairing information during
contig generation, and includes several improvements to read trimming. Together
these changes result in substantially improved assemblies that recover a more
complete set of reference genes than previous methods.
  Availability: A5-miseq is licensed under the GPL open source license. Source
code and precompiled binaries for Mac OS X 10.6+ and Linux 2.6.15+ are
available from http://sourceforge.net/projects/ngopt


  We derive and apply a partial differential equation for the moment generating
function of the Wright-Fisher model of population genetics.


  Population genetics theory has laid the foundations for genomics analyses
including the recent burst in genome scans for selection and statistical
inference of past demographic events in many prokaryote, animal and plant
species. Identifying SNPs under natural selection and underpinning species
adaptation relies on disentangling the respective contribution of random
processes (mutation, drift, migration) from that of selection on nucleotide
variability. Most theory and statistical tests have been developed using the
Kingman coalescent theory based on the Wright-Fisher population model. However,
these theoretical models rely on biological and life-history assumptions which
may be violated in many prokaryote, fungal, animal or plant species. Recent
theoretical developments of the so called multiple merger coalescent models are
reviewed here ({\Lambda}-coalescent, beta-coalescent, Bolthausen-Snitzman,
{\Xi}-coalescent). We explicit how these new models take into account various
pervasive ecological and biological characteristics, life history traits or
life cycles which were not accounted in previous theories such as 1) the skew
in offspring production typical of marine species, 2) fast adapting
microparasites (virus, bacteria and fungi) exhibiting large variation in
population sizes during epidemics, 3) the peculiar life cycles of fungi and
bacteria alternating sexual and asexual cycles, and 4) the high rates of
extinction-recolonization in spatially structured populations. We finally
discuss the relevance of multiple merger models for the detection of SNPs under
selection in these species, for population genomics of very large sample size
and advocate to potentially examine the conclusion of previous population
genetics studies.


  The de Bruijn graph plays an important role in bioinformatics, especially in
the context of de novo assembly. However, the representation of the de Bruijn
graph in memory is a computational bottleneck for many assemblers. Recent
papers proposed a navigational data structure approach in order to improve
memory usage. We prove several theoretical space lower bounds to show the
limitation of these types of approaches. We further design and implement a
general data structure (DBGFM) and demonstrate its use on a human whole-genome
dataset, achieving space usage of 1.5 GB and a 46% improvement over previous
approaches. As part of DBGFM, we develop the notion of frequency-based
minimizers and show how it can be used to enumerate all maximal simple paths of
the de Bruijn graph using only 43 MB of memory. Finally, we demonstrate that
our approach can be integrated into an existing assembler by modifying the
ABySS software to use DBGFM.


  The force induced unzipping transition of a double stranded DNA is considered
from a purely thermodynamic point of view. This analysis provides us with a set
of relations that can be used to test microscopic theories and experiments. The
thermodynamic approach is based on the hypothesis of impenetrability of the
force in the zipped state. The melting and the unzipping transitions are
considered in the same framework and compared with the existing statistical
model results. The analysis is then extended to a possible continuous unzipping
transition.


  The three-dimensional conformations of non-coding RNAs underpin their
biochemical functions but have largely eluded experimental characterization.
Here, we report that integrating a classic mutation/rescue strategy with
high-throughput chemical mapping enables rapid RNA structure inference with
unusually strong validation. We revisit a paradigmatic 16S rRNA domain for
which SHAPE (selective 2`-hydroxyl acylation with primer extension) suggested a
conformational change between apo- and holo-ribosome conformations.
Computational support estimates, data from alternative chemical probes, and
mutate-and-map (M2) experiments expose limitations of prior methodology and
instead give a near-crystallographic secondary structure. Systematic
interrogation of single base pairs via a high-throughput mutation/rescue
approach then permits incisive validation and refinement of the M2-based
secondary structure and further uncovers the functional conformation as an
excited state (25+/-5% population) accessible via a single-nucleotide register
shift. These results correct an erroneous SHAPE inference of a ribosomal
conformational change and suggest a general mutate-map-rescue approach for
dissecting RNA dynamic structure landscapes.


  We present a stochastic formalism for signal transduction processes in
bacterial two-component system. Using elementary mass action kinetics, the
proposed model takes care of signal transduction in terms of phosphotransfer
mechanism between the cognate partners of a two-component system, viz, the
sensor kinase and the response regulator. Based on the difference in
functionality of the sensor kinase, the noisy phosphotransfer mechanism has
been studied for monofunctional and bifunctional two component system using the
formalism of linear noise approximation. Steady state analysis of both models
quantifies different physically realizable quantities, e.g., variance,
coefficient of variation, mutual information. The resultant data reveals that
both systems reliably transfer information of extra-cellular environment under
low external stimulus and at high kinase and phosphatase regime. We extend our
analysis further by studying the role of two-component system in downstream
gene regulation.


  Model reduction is a central topic in systems biology and dynamical systems
theory, for reducing the complexity of detailed models, finding important
parameters, and developing multi-scale models for instance. While perturbation
theory is a standard mathematical tool to analyze the different time scales of
a dynamical system, and decompose the system accordingly, tropical methods
provide a simple algebraic framework to perform these analyses systematically
in polynomial systems. The crux of these tropicalization methods is in the
computation of tropical equilibrations. In this paper we show that
constraint-based methods, using reified constraints for expressing the
equilibration conditions, make it possible to numerically solve non-linear
tropical equilibration problems, out of reach of standard computation methods.
We illustrate this approach first with the reduction of simple biochemical
mechanisms such as the Michaelis-Menten and Goldbeter-Koshland models, and
second, with performance figures obtained on a large scale on the model
repository \texttt{biomodels.net}


  Immigration can rescue local populations from extinction, helping to
stabilise a metapopulation. Local population dynamics is important for
determining the strength of this rescue effect, but the mechanistic link
between local demographic parameters and the rescue effect at the
metapopulation level has received very little attention by modellers. We
develop an analytical framework that allows us to describe the emergence of the
rescue effect from interacting local stochastic dynamics. We show this
framework to be applicable to a wide range of spatial scales, providing a
powerful and convenient alternative to individual-based models for making
predictions concerning the fate of metapopulations. We show that the rescue
effect plays an important role in minimising the increase in local extinction
probability associated with high demographic stochasticity, but its role is
more limited in the case of high local environmental stochasticity of
recruitment or survival. While most models postulate the rescue effect, our
framework provides an explicit mechanistic link between local dynamics and the
emergence of the rescue effect, and more generally the stability of the whole
metapopulation.


  A number of methods have been developed to infer differential rates of
species diversification through time and among clades using time-calibrated
phylogenetic trees. However, we lack a general framework that can delineate and
quantify heterogeneous mixtures of dynamic processes within single phylogenies.
I developed a method that can identify arbitrary numbers of time-varying
diversification processes on phylogenies without specifying their locations in
advance. The method uses reversible-jump Markov Chain Monte Carlo to move
between model subspaces that vary in the number of distinct diversification
regimes. The model assumes that changes in evolutionary regimes occur across
the branches of phylogenetic trees under a compound Poisson process and
explicitly accounts for rate variation through time and among lineages. Using
simulated datasets, I demonstrate that the method can be used to quantify
complex mixtures of time-dependent, diversity-dependent, and constant-rate
diversification processes. I compared the performance of the method to the
MEDUSA model of rate variation among lineages. As an empirical example, I
analyzed the history of speciation and extinction during the radiation of
modern whales. The method described here will greatly facilitate the
exploration of macroevolutionary dynamics across large phylogenetic trees,
which may have been shaped by heterogeneous mixtures of distinct evolutionary
processes.


  We study an epidemic propagation between $M$ population centra. The novelty
of the model is in analyzing the migration of host (remaining in the same
centre) and guest (migrated to another centre) populations separately. Even in
the simplest case $M=2$, this modification is justified because it gives a more
realistic description of migration processes. This becomes evident in a purely
migration model with vanishing epidemic parameters. It is important to account
for a certain number of guest susceptible present in non-host cenrta because
these susceptible may be infected and return to the host node as infectives.
The flux of such infectives is not negligible and is comparable with the flux
of host infectives migrated to other centra, because the return rate of a guest
individual will, by nature, tend to be high. It is shown that taking account of
both fluxes of infectives noticeably increases the speed of epidemic spread in
a 1D lattice of identical SIR centra.


  Consider a graph in which each site is endowed with a value called
\emph{fitness}. A path in the graph is said to be "open" or "accessible" if the
fitness values along that path is strictly increasing. We say that there is
accessibility percolation between two sites when such a path between them
exists. Motivated by the so called House-of-Cards model from evolutionary
biology, we consider this question on the $L$-hypercube $\{0,1\}^L$ where the
fitness values are independent random variables. We show that, in the large $L$
limit, the probability that an accessible path exists from an arbitrary
starting point to the (random) fittest site is no more than $x^*_{1/2}=
1-\frac12\sinh^{-1}(2) =0.27818\ldots$ and we conjecture that this probability
does converge to $x^*_{1/2}$. More precisely, there is a phase transition on
the value of the fitness $x$ of the starting site: assuming that the fitnesses
are uniform in $[0,1]$, we show that, in the large $L$ limit, there is almost
surely no path to the fittest site if $x>x^*_{1/2}$ and we conjecture that
there are almost surely many paths if $x<x^*_{1/2}$.
  If one conditions on the fittest site to be on the opposite corner of the
starting site rather than being randomly chosen, the picture remains the same
but with the critical point being now $x^*_1= 1-\sinh^{-1}(1)= 0.11863\ldots$.
Along the way, we obtain a large $L$ estimation for the number of self-avoiding
paths joining two opposite corners of the $L$-hypercube.


  We demonstrate how nanojunctions between lysosomes and sarcoplasmic reticulum
(L-SR junctions) serve to couple lysosomal activation to regenerative,
ryanodine receptor-mediated cellular calcium (Ca2+) waves. In pulmonary artery
smooth muscle cells (PASMCs) nicotinic acid adenine dinucleotide phosphate
(NAADP) may trigger increases in cytoplasmic Ca2+ via L-SR junctions, in a
manner that requires initial Ca2+ release from lysosomes and subsequent
Ca2+-induced Ca2+ release (CICR) via ryanodine receptor (RyR) subtype 3 on the
SR membrane proximal to lysosomes. L-SR junction membrane separation has been
estimated to be <400 nm and thus beyond the resolution of light microscopy.
This study utilizes transmission electron microscopy to provide a thorough
ultrastructural characterization of the L-SR junctions in PASMCs. These
junctions are prominent features in these cells and we estimate that the
membrane separation and extension are about 15 nm and 300 nm, respectively. We
also develop a quantitative model of the L-SR junction using these
measurements, prior kinetic and specific Ca2+ signal information as input data.
Simulations of NAADP-dependent junctional Ca2+ transients show that the
magnitude of these signals can breach the threshold for CICR via RyR3. By
correlation analysis of live cell Ca2+ signals and simulated L-SR junctional
Ca2+ transients, we estimate that "trigger zones" with a 60-100 junctions are
required to confer a signal of similar magnitude. This is compatible with the
130 lysosomes/cell estimated from our ultrastructural observations. Most
importantly, our model shows that increasing the L-SR junctional width above 50
nm lowers the magnitude of junctional [Ca2+] such that there is a failure to
breach the threshold for CICR via RyR3. L-SR junctions are therefore a
pre-requisite for efficient Ca2+ signal coupling and may contribute to cellular
function in health and disease.


  We investigated the relationships between search efficiency, movement
strategy, and non-local communication in the biological context of animal
foraging. We considered situations where the members of a population of
foragers perform either Gaussian jumps or Levy flights, and show that the
search time is minimized when communication among individuals occurs at
intermediate ranges, independently of the type of movement. Additionally, while
Brownian strategies are more strongly influenced by the communication
mechanism, Levy flights still result in shorter overall search durations.


  A novel uracil-DNA degrading protein factor (termed UDE) was identified in
Drosophila melanogaster with no significant structural and functional homology
to other uracil-DNA binding or processing factors. Determination of the 3D
structure of UDE will be a true breakthrough in description of the molecular
mechanism of action of UDE catalysis, as well as in general uracil-recognition
and nuclease action. The revolutionary ESPRIT technology was applied to the
novel protein UDE to overcome problems in identifying soluble expressing
constructs given the absence of precise information on domain content and
arrangement. Nine specimen from the created numerous truncated constructs of
UDE were choosen to dechiper structural and functional relationships. VUVCD
with neural network was performed to define the secondary structure content and
location of UDE and its truncated variants. The quantitative analysis
demonstrated exclusive {\alpha}-helical content for the full-length protein,
which is preserved in the truncated constructs. Partition of {\alpha}-helical
boundles comparison with the truncated protein segments denoted new domain
boundaries which differ from the conserved motifs determined by sequence-based
alignment of UDE homologues in pupating insects. Here we demonstrate that
combination of ESPRIT and VUVCD with NN provided significant structural
description of UDE and resulted useful truncated constructs for further
detailed functional studies.


  We have used whole genome paired-end Illumina sequence data to identify
tandem duplications in 20 isofemale lines of D. yakuba, and 20 isofemale lines
of D. simulans and performed genome wide validation with PacBio long molecule
sequencing. We identify 1,415 tandem duplications that are segregating in D.
yakuba as well as 975 duplications in D. simulans, indicating greater variation
in D. yakuba. Additionally, we observe high rates of secondary deletions at
duplicated sites, with 8% of duplicated sites in D. simulans and 17% of sites
in D. yakuba modified with deletions. These secondary deletions are consistent
with the action of the large loop mismatch repair system acting to remove
polymorphic tandem duplication, resulting in rapid dynamics of gain and loss in
duplicated alleles and a richer substrate of genetic novelty than has been
previously reported. Most duplications are present in only single strains,
suggesting deleterious impacts are common. D. simulans shows larger numbers of
whole gene duplications in comparison to larger proportions of gene fragments
in D. yakuba. D. simulans displays an excess of high frequency variants on the
X chromosome, consistent with adaptive evolution through duplications on the D.
simulans X or demographic forces driving duplicates to high frequency. We
identify 78 chimeric genes in D. yakuba and 38 chimeric genes in D. simulans,
as well as 143 cases of recruited non-coding sequence in D. yakuba and 96 in D.
simulans, in agreement with rates of chimeric gene origination in D.
melanogaster. Together, these results suggest that tandem duplications often
result in complex variation beyond whole gene duplications that offers a rich
substrate of standing variation that is likely to contribute both to
detrimental phenotypes and disease, as well as to adaptive evolutionary change.


  The relationship between epidemiology, mathematical modeling and
computational tools allows to build and test theories on the development and
battling of a disease. This PhD thesis is motivated by the study of
epidemiological models applied to infectious diseases in an Optimal Control
perspective, giving particular relevance to Dengue. Dengue is a subtropical and
tropical disease transmitted by mosquitoes, that affects about 100 million
people per year and is considered by the World Health Organization a major
concern for public health. The mathematical models developed and tested in this
work, are based on ordinary differential equations that describe the dynamics
underlying the disease, including the interaction between humans and
mosquitoes. An analytical study is made related to equilibrium points, their
stability and basic reproduction number. The spreading of Dengue can be
attenuated through measures to control the transmission vector, such as the use
of specific insecticides and educational campaigns. Since the development of a
potential vaccine has been a recent global bet, models based on the simulation
of a hypothetical vaccination process in a population are proposed. Based on
Optimal Control theory, we have analyzed the optimal strategies for using these
controls, and respective impact on the reduction/eradication of the disease
during an outbreak in the population, considering a bioeconomic approach. The
formulated problems are numerically solved using direct and indirect methods.
The first discretize the problem turning it into a nonlinear optimization
problem. Indirect methods use the Pontryagin Maximum Principle as a necessary
condition to find the optimal curve for the respective control. In these two
strategies several numerical software packages are used.


  Advances in DNA sequencing technology have stimulated the development of
algorithms and tools for processing very large collections of short strings
(reads). Short-read alignment and assembly are among the most well-studied
problems. Many state-of-the-art aligners, at their core, have used the
Burrows-Wheeler transform (BWT) as a main-memory index of a reference genome
(typical example, NCBI human genome). Recently, BWT has also found its use in
string-graph assembly, for indexing the reads (i.e., raw data from DNA
sequencers). In a typical data set, the volume of reads is tens of times of the
sequenced genome and can be up to 100 Gigabases. Note that a reference genome
is relatively stable and computing the index is not a frequent task. For reads,
the index has to computed from scratch for each given input. The ability of
efficient BWT construction becomes a much bigger concern than before. In this
paper, we present a practical method called CX1 for constructing the BWT of
very large string collections. CX1 is the first tool that can take advantage of
the parallelism given by a graphics processing unit (GPU, a relative cheap
device providing a thousand or more primitive cores), as well as simultaneously
the parallelism from a multi-core CPU and more interestingly, from a cluster of
GPU-enabled nodes. Using CX1, the BWT of a short-read collection of up to 100
Gigabases can be constructed in less than 2 hours using a machine equipped with
a quad-core CPU and a GPU, or in about 43 minutes using a cluster with 4 such
machines (the speedup is almost linear after excluding the first 16 minutes for
loading the reads from the hard disk). The previously fastest tool BRC is
measured to take 12 hours to process 100 Gigabases on one machine; it is
non-trivial how BRC can be parallelized to take advantage a cluster of
machines, let alone GPUs.


  When long-lived, balancing selection can lead to trans-species polymorphisms
that are shared by two or more species identical by descent. In this case, the
gene genealogies at the selected sites cluster by allele instead of by species
and, because of linkage, nearby neutral sites also have unusual genealogies.
Although it is clear that this scenario should lead to discernible footprints
in genetic variation data, notably the presence of additional neutral
polymorphisms shared between species and the absence of fixed differences, the
effects remain poorly characterized. We focus on the case of a single site
under long-lived balancing selection and derive approximations for summaries of
the data that are sensitive to a trans-species polymorphism: the length of the
segment that carries most of the signals, the expected number of shared neutral
SNPs within the segment and the patterns of allelic associations among them.
Coalescent simulations of ancient balancing selection confirm the accuracy of
our approximations. We further show that for humans and chimpanzees, and more
generally for pairs of species with low genetic diversity levels, the patterns
of genetic variation on which we focus are highly unlikely to be generated by
neutral recurrent mutations, so these statistics are specific as well as
sensitive. We discuss the implications of our results for the design and
interpretation of genome scans for ancient balancing selection in apes and
other taxa.


  In addressing the origins of Darwinian evolution, recent experimental work
has been focussed on the discovery of simple physical effects which would
provide a relevant selective advantage to protocells competing with each other
for a limited supply of lipid. In particular, data coming from Szostak's lab
suggest that the transition from simple prebiotically plausible lipid membranes
to more complex and heterogeneous ones, closer to real biomembranes, may have
been driven by changes in the fluidity of the membrane and its affinity for the
available amphiphilic compound, which in turn would involve changes in vesicle
growth dynamics. Earlier work from the same group reported osmotically-driven
competition effects, whereby swelled vesicles grow at the expense of isotonic
ones. In this paper, we try to expand on these experimental studies by
providing a simple mathematical model of a population of competing vesicles,
studied at the level of lipid kinetics. In silico simulations of the model are
able to reproduce qualitatively and often quantitatively the experimentally
reported competition effects in both scenarios. We also develop a method for
numerically solving the equilibrium of a population of competing model
vesicles, which is quite general and applicable to different vesicle kinetics
schemes.


  The mechanosensitive channel of large conductance, which serves as a model
system for mechanosensitive channels, has previously been crystallized in the
closed form, but not in the open form. Ensemble measurements and
electrophysiological sieving experiments show that the open-diameter of the
channel pore is >25{\AA}, but the exact size and whether the conformational
change follows a helix-tilt or barrel-stave model are unclear. Here we report
measurements of the distance changes on liposome-reconstituted MscL
transmembrane {\alpha}-helices, using a "virtual sorting" single-molecule
fluorescence energy transfer. We observed directly that the channel opens via
the helix-tilt model and the open pore reaches 2.8 nm in diameter. In addition,
based on the measurements, we developed a molecular dynamics model of the
channel structure in the open state which confirms our direct observations.


  Unstable dynamics characterizes the evolution of most solid tumors. Because
of an increased failure of maintaining genome integrity, a cumulative increase
in the levels of gene mutation and loss is observed. Previous work suggests
that instability thresholds to cancer progression exist, defining phase
transition phenomena separating tumor-winning scenarios from tumor extinction
or coexistence phases. Here we present an integral equation approach to the
quasispecies dynamics of unstable cancer. The model exhibits two main phases,
characterized by either the success or failure of cancer tissue. Moreover, the
model predicts that tumor failure can be due to either a reduced selective
advantage over healthy cells or excessive instability. We also derive an
approximate, analytical solution that predicts the front speed of aggressive
tumor populations on the instability space.


  Background: Trace quantities of contaminating DNA are widespread in the
laboratory environment, but their presence has received little attention in the
context of high throughput sequencing. This issue is highlighted by recent
works that have rested controversial claims upon sequencing data that appear to
support the presence of unexpected exogenous species.
  Results: I used reads that preferentially aligned to alternate genomes to
infer the distribution of potential contaminant species in a set of independent
sequencing experiments. I confirmed that dilute samples are more exposed to
contaminating DNA, and, focusing on four single-cell sequencing experiments,
found that these contaminants appear to originate from a wide diversity of
clades. Although negative control libraries prepared from "blank" samples
recovered the highest-frequency contaminants, low-frequency contaminants, which
appeared to make heterogeneous contributions to samples prepared in parallel
within a single experiment, were not well controlled for. I used these results
to show that, despite heavy replication and plausible controls, contamination
can explain all of the observations used to support a recent claim that
complete genes pass from food to human blood.
  Conclusions: Contamination must be considered a potential source of signals
of exogenous species in sequencing data, even if these signals are replicated
in independent experiments, vary across conditions, or indicate a species which
seems a priori unlikely to contaminate. Negative control libraries processed in
parallel are essential to control for contaminant DNAs, but their limited
ability to recover low-frequency contaminants must be recognized.


  During Saccharomyces cerevisiae closed mitosis, parental identity is
sustained by the asymmetric segregation of ageing factors. Such asymmetry has
been hypothesized to occur via diffusion barriers, constraining protein lateral
exchange in cellular membranes. Diffusion barriers have been extensively
studied in the plasma membrane, but their identity and organization within the
nucleus remain unknown. Here, we propose how sphingolipid domains, protein
rings, and morphological changes of the nucleus may coordinate to restrict
protein exchange between nuclear lobes. Our spatial stochastic model is based
on several lines of experimental evidence and predicts that, while a
sphingolipid domain and a protein ring could constitute the barrier during
early anaphase; a sphingolipid domain spanning the bridge between lobes during
late anaphase would be entirely sufficient. Additionally, we explore the
structural organization of plausible diffusion barriers. Our work shows how
nuclear diffusion barriers in closed mitosis may be emergent properties of
simple nanoscale biophysical interactions.


  Inverse parameter estimation of process-based models is a long-standing
problem in many scientific disciplines. A key question for inverse parameter
estimation is how to define the metric that quantifies how well model
predictions fit to the data. This metric can be expressed by general cost or
objective functions, but statistical inversion methods require a particular
metric, the probability of observing the data given the model parameters, known
as the likelihood.
  For technical and computational reasons, likelihoods for process-based
stochastic models are usually based on general assumptions about variability in
the observed data, and not on the stochasticity generated by the model. Only in
recent years have new methods become available that allow the generation of
likelihoods directly from stochastic simulations. Previous applications of
these approximate Bayesian methods have concentrated on relatively simple
models. Here, we report on the application of a simulation-based likelihood
approximation for FORMIND, a parameter-rich individual-based model of tropical
forest dynamics.
  We show that approximate Bayesian inference, based on a parametric likelihood
approximation placed in a conventional Markov chain Monte Carlo (MCMC) sampler,
performs well in retrieving known parameter values from virtual inventory data
generated by the forest model. We analyze the results of the parameter
estimation, examine its sensitivity to the choice and aggregation of model
outputs and observed data (summary statistics), and demonstrate the application
of this method by fitting the FORMIND model to field data from an Ecuadorian
tropical forest. Finally, we discuss how this approach differs from approximate
Bayesian computation (ABC), another method commonly used to generate
simulation-based likelihood approximations.
  Our results demonstrate that simulation-based inference, [...]


  Zinc-finger nucleases (ZFNs) and transcription activator-like effector
nucleases (TALENs) comprise a powerful class of tools that are redefining the
boundaries of biological research. Although these technologies have begun to
enable targeted genome modifications, there remains a need for new technologies
that are affordable, scalable, and easy to engineer. In this paper, we propose
a new tool for genetic engineering, the pseudocomplementary peptide nucleic
acid nucleases (pcPNANs), which are composed of a pseudocomplementary PNA
(pcPNA) specific for a DNA target sequence, a FokI nuclease cleavage domain and
a nuclear localization signal. pcPNANs may induce targeted DNA double-strand
breaks that activate DNA damage response pathways and enable custom
alterations. Their cleavage-site is determined by simple Watson-Crick rule, and
thus pcPNANs for aimed cleavage of genomes can be straightforwardly designed
and synthesized without any selection procedure. Accordingly, the cleavage-site
and site-specificity are freely chosen by changing the sequences and the
lengths of pcPNA strands. We believe that the potentiality of pcPNAN as a new
tool for genetic engineering will be confirmed in the future.


  The asymmetric Hopfield model is used to simulate signaling dynamics in
gene/transcription factor networks. The model allows for a direct mapping of a
gene expression pattern into attractor states. We analyze different control
strategies aiming at disrupting attractor patterns using selective local fields
representing therapeutic interventions. The control strategies are based on the
identification of signaling $bottlenecks$, which are single nodes or strongly
connected clusters of nodes that have a large impact on the signaling. We
provide a theorem with bounds on the minimum number of nodes that guarantee
controllability of bottlenecks consisting of strongly connected components. The
control strategies are applied to the identification of sets of proteins that,
when inhibited, selectively disrupt the signaling of cancer cells while
preserving the signaling of normal cells. We use an experimentally validated
non-specific network and a specific B cell interactome reconstructed from gene
expression data to model cancer signaling in lung and B cells, respectively.
This model could help in the rational design of novel robust therapeutic
interventions based on our increasing knowledge of complex gene signaling
networks.


  Extortion strategies can dominate any opponent in an iterated prisoner's
dilemma game. But if players are able to adopt the strategies performing
better, extortion becomes widespread and evolutionary unstable. It may
sometimes act as a catalyst for the evolution of cooperation, and it can also
emerge in interactions between two populations, yet it is not the evolutionary
stable outcome. Here we revisit these results in the realm of spatial games. We
find that pairwise imitation and birth-death dynamics return known evolutionary
outcomes. Myopic best response strategy updating, on the other hand, reveals
new counterintuitive solutions. Defectors and extortioners coarsen
spontaneously, which allows cooperators to prevail even at prohibitively high
temptations to defect. Here extortion strategies play the role of a Trojan
horse. They may emerge among defectors by chance, and once they do, cooperators
become viable as well. These results are independent of the interaction
topology, and they highlight the importance of coarsening, checkerboard
ordering, and best response updating in evolutionary games.


  Next-generation sequencing technology enables the identification of thousands
of gene regulatory sequences in many cell types and organisms. We consider the
problem of testing if two such sequences differ in their number of binding site
motifs for a given transcription factor (TF) protein. Binding site motifs
impart regulatory function by providing TFs the opportunity to bind to genomic
elements and thereby affect the expression of nearby genes. Evolutionary
changes to such functional DNA are hypothesized to be major contributors to
phenotypic diversity within and between species; but despite the importance of
TF motifs for gene expression, no method exists to test for motif loss or gain.
Assuming that motif counts are Binomially distributed, and allowing for
dependencies between motif instances in evolutionarily related sequences, we
derive the probability mass function of the difference in motif counts between
two nucleotide sequences. We provide a method to numerically estimate this
distribution from genomic data and show through simulations that our estimator
is accurate. Finally, we introduce the R package {\tt motifDiverge} that
implements our methodology and illustrate its application to gene regulatory
enhancers identified by a mouse developmental time course experiment. While
this study was motivated by analysis of regulatory motifs, our results can be
applied to any problem involving two correlated Bernoulli trials.


  The paper discusses a connection between asymmetric reproduction -- that is
reproduction in a parent-child relationship where the parent does not mutate
during reproduction --, the fact that all non-viral lifeforms bear genes of
their reproduction machinery and how this could relate to evolutionary
mechanisms behind aging. In a highly simplified model of the evolution process
rules are derived under which aging is an important factor of the adaption in
the evolution process and what groups of life-forms necessarily have to age and
where exceptions from that rule are possible.


  The current capacity of computers makes it possible to perform simulations of
small systems with portable, explicit-solvent potentials achieving high degree
of accuracy. However, simplified models must be employed to exploit the
behaviour of large systems or to perform systematic scans of smaller systems.
While powerful algorithms are available to facilitate the sampling of the
conformational space, successful applications of such models are hindered by
the availability of simple enough potentials able to satisfactorily reproduce
known properties of the system. We develop an interatomic potential to account
for a number of properties of proteins in a computationally economic way. The
potential is defined within an all-atom, implicit solvent model by contact
functions between the different atom types. The associated numerical values can
be optimised by an iterative Monte Carlo scheme on any available experimental
data, provided that they are expressible as thermal averages of some
conformational properties. We test this model on three different proteins, for
which we also perform a scan of all possible point mutations with explicit
conformational sampling. The resulting models, optimised solely on a subset of
native distances, not only reproduce the native conformations within a few
Angstroms from the experimental ones, but show the cooperative transition
between native and denatured state and correctly predict the measured
free--energy changes associated with point mutations. Moreover, differently
from other structure-based models, our method leaves a residual degree of
frustration, which is known to be present in protein molecules.


  Biopolymer length regulation is a complex process that involves a large
number of subprocesses acting simultaneously across multiple spatial and
temporal scales. An illustrative example important for genomic stability is the
length regulation of telomeres---nucleo-protein structures at the ends of
linear chromosomes. Maintenance of telomeres is often facilitated by the enzyme
telomerase but, particularly in telomerase-free systems, the maintenance of
chromosomal termini depends on alternative lengthening of telomeres (ALT)
mechanisms mediated by recombination. Various linear and circular DNA
structures were identified to participate in ALT, however, dynamics of the
whole process is still poorly understood. We propose a chemical kinetics model
of ALT with kinetic rates systematically derived from the biophysics of DNA
diffusion and looping. The reaction system is reduced to a
coagulation-fragmentation system by quasi-steady state approximation. The
detailed treatment of kinetic rates yields explicit formulae for expected size
distributions of telomeres that demonstrate the key role played by the
J-factor, a quantitative measure of bending of polymers. The results are in
agreement with experimental data and point out interesting phenomena: an
appearance of very long telomeric circles if the total telomere density exceeds
a critical value (excess mass) and a nonlinear response of the telomere size
distributions to the amount of telomeric DNA in the system. The results can be
of general importance for understanding dynamics of telomeres in
telomerase-independent systems as this mode of telomere maintenance is similar
to the situation in tumor cells lacking telomerase activity. Furthermore, due
to its universality, the model may also serve as a prototype of an interaction
between linear and circular DNA structures in various settings.


  Jim Shapiro synthesizes a great many observations about the mechanisms of
evolution to reach the remarkable conclusion that large-scale modification,
exchange, and rearrangement of the genome are common and should be viewed as
fundamental features of life. In other words, the genome should be viewed not
as mostly read-only with a few rare mutations, but rather as a fully-fledged
read-write library of genetic functions under continuous revision. Revision of
the genome occurs during cellular replication, during multicellular
development, and during evolution of a population of individuals. DNA
formatting controls the timing and location of genetic rearrangements, gene
expression, and genetic repair. Each of these events is under the control of
precise cellular circuits. Shapiro reviews the toolbox of natural genetic
engineering that provides the functionalities necessary for efficient long-term
genome restructuring.


  Human associated microbial communities exert tremendous influence over human
health and disease. With modern metagenomic sequencing methods it is possible
to follow the relative abundance of microbes in a community over time. These
microbial communities exhibit rich ecological dynamics and an important goal of
microbial ecology is to infer the interactions between species from sequence
data. Any algorithm for inferring species interactions must overcome three
obstacles: 1) a correlation between the abundances of two species does not
imply that those species are interacting, 2) the sum constraint on the relative
abundances obtained from metagenomic studies makes it difficult to infer the
parameters in timeseries models, and 3) errors due to experimental uncertainty,
or mis-assignment of sequencing reads into operational taxonomic units, bias
inferences of species interactions. Here we introduce an approach, Learning
Interactions from MIcrobial Time Series (LIMITS), that overcomes these
obstacles. LIMITS uses sparse linear regression with boostrap aggregation to
infer a discrete-time Lotka-Volterra model for microbial dynamics. We tested
LIMITS on synthetic data and showed that it could reliably infer the topology
of the inter-species ecological interactions. We then used LIMITS to
characterize the species interactions in the gut microbiomes of two individuals
and found that the interaction networks varied significantly between
individuals. Furthermore, we found that the interaction networks of the two
individuals are dominated by distinct "keystone species", Bacteroides fragilis
and Bacteroided stercosis, that have a disproportionate influence on the
structure of the gut microbiome even though they are only found in moderate
abundance. Based on our results, we hypothesize that the abundances of certain
keystone species may be responsible for individuality in the human gut
microbiome.


  High-throughput DNA sequencers are becoming indispensable in our
understanding of diseases at molecular level, in marker-assisted selection in
agriculture and in microbial genetics research. These sequencing instruments
produce enormous amount of data (often terabytes of raw data in a month) that
requires efficient analysis, management and interpretation. The commonly used
sequencing instrument today produces billions of short reads (upto 150 bases)
from each run. The first step in the data analysis step is alignment of these
short reads to the reference genome of choice. There are different open source
algorithms available for sequence alignment to the reference genome. These
tools normally have a high computational overhead, both in terms of number of
processors and memory. Here, we propose a hybrid-computing environment called
MUSIC (Mapping USIng hybrid Computing) for one of the most popular open source
sequence alignment algorithm, BWA, using accelerators that show significant
improvement in speed over the serial code.


  Enormous volumes of short reads data from next-generation sequencing (NGS)
technologies have posed new challenges to the area of genomic sequence
comparison.
  The multiple sequence alignment approach is hardly applicable to NGS data due
to the challenging problem of short read assembly.
  Thus alignment-free methods need to be developed for the comparison of NGS
samples of short reads.
  Recently, new $k$-mer based distance measures such as {\it CVTree},
$d_{2}^{S}$, {\it co-phylog} have been proposed to address this problem.
  However, those distances depend considerably on the parameter $k$, and how to
choose the optimal $k$ is not trivial since it may depend on different aspects
of the sequence data.
  Hence, in this paper we consider an alternative parameter-free approach:
compression-based distance measures.
  These measures have shown impressive performance on long genome sequences in
previous studies, but they have not been tested on NGS short reads.
  In this study we perform extensive validation and show that the
compression-based distances are highly consistent with those distances obtained
from the $k$-mer based methods, from the alignment-based approach, and from
existing benchmarks in the literature.
  Moreover, as these measures are parameter-free, no optimization is required
and they still perform consistently well on multiple types of sequence data,
for different kinds of species and taxonomy levels.
  The compression-based distance measures are assembly-free, alignment-free,
parameter-free, and thus represent useful tools for the comparison of long
genome sequences and NGS samples of short reads.


  A systematic study of electron or hole transfer along DNA dimers, trimers and
polymers is presented with a tight-binding approach at the base-pair level,
using the relevant on-site energies of the base-pairs and the hopping
parameters between successive base-pairs. A system of $N$ coupled differential
equations is solved numerically with the eigenvalue method, allowing the
temporal and spatial evolution of electrons or holes along a $N$ base-pair DNA
segment to be determined. Useful physical quantities are defined and calculated
including the maximum transfer percentage $p$ and the pure maximum transfer
rate $\frac{p}{T}$ for cases where a period $T$ can be defined, as well as the
pure mean carrier transfer rate $k$ and the speed of charge transfer $u=kd$,
where $d = N \times$ 3.4 {\AA} is the charge transfer distance. The inverse
decay length $\beta$ used for the exponential fit $k = k_0 \exp(-\beta d)$ and
the exponent $\eta$ used for the power law fit $k = k_0' N^{-\eta}$ are
computed. The electron and hole transfer along polymers including
poly(dG)-poly(dC), poly(dA)-poly(dT), GCGCGC..., ATATAT... is studied, too.
$\beta$ falls in the range $\approx$ 0.2-2 {\AA}$^{-1}$, $k_0$ is usually
10$^{-2}$-10$^{-1}$ PHz although, generally, it falls in the wider range
10$^{-4}$-10 PHz. $\eta$ falls in the range $\approx$ 1.7-17, $k_0'$ is usually
$\approx 10^{-2}$-10$^{-1}$ PHz, although generally, it falls in the wider
range $\approx10^{-4}$-10$^3$ PHz. Finally, the results are compared with the
predictions of Wang et al. Phys. Rev. Lett. 93, (2004) 016401, as well as
experiments, including Murphy et al. Science 262, 1025 (1993); Arkin et al.
Science 273, 475 (1996); Giese et al. Angew. Chem. Int. Ed. 38, 996 (1999);
Giese et al. Nature 412, 318 (2001). This method allows to assess the extent at
which a specific DNA segment can serve as an efficient medium for charge
transfer.


  We show that in the generic situation where a biological network, e.g. a
protein interaction network, is in fact a subnetwork embedded in a larger
"bulk" network, the presence of the bulk causes not just extrinsic noise but
also memory effects. This means that the dynamics of the subnetwork will depend
not only on its present state, but also its past. We use projection techniques
to get explicit expressions for the memory functions that encode such memory
effects, for generic protein interaction networks involving binary and unary
reactions such as complex formation and phosphorylation, respectively.
Remarkably, in the limit of low intrinsic copy-number noise such expressions
can be obtained even for nonlinear dependences on the past. We illustrate the
method with examples from a protein interaction network around epidermal growth
factor receptor (EGFR), which is relevant to cancer signalling. These examples
demonstrate that inclusion of memory terms is not only important conceptually
but also leads to substantially higher quantitative accuracy in the predicted
subnetwork dynamics.


  Cancer cells are known to modify their micro-environment such that it can
sustain a larger population, or, in ecological terms, they construct a niche
which increases the carrying capacity of the population. It has however been
argued that niche construction, which benefits all cells in the tumour, would
be selected against since cheaters could reap the benefits without paying the
cost. We have investigated the impact of niche specificity on tumour evolution
using an individual based model of breast tumour growth, in which the carrying
capacity of each cell consists of two components: an intrinsic,
subclone-specific part and a contribution from all neighbouring cells. Analysis
of the model shows that the ability of a mutant to invade a resident population
depends strongly on the specificity. When specificity is low selection is
mostly on growth rate, while high specificity shifts selection towards
increased carrying capacity. Further, we show that the long-term evolution of
the system can be predicted using adaptive dynamics. By comparing the results
from a spatially structured vs.\ well-mixed population we show that spatial
structure restores selection for carrying capacity even at zero specificity,
which a poses solution to the niche construction dilemma. Lastly, we show that
an expanding population exhibits spatially variable selection pressure, where
cells at the leading edge exhibit higher growth rate and lower carrying
capacity than those at the centre of the tumour.


  The detection of somatic single nucleotide variants is a crucial component to
the characterization of the cancer genome. Mutation calling algorithms thus far
have focused on comparing the normal and tumor genomes from the same
individual. In recent years, it has become routine for projects like The Cancer
Genome Atlas (TCGA) to also sequence the tumor RNA. Here we present RADIA (RNA
and DNA Integrated Analysis), a method that combines the patient-matched normal
and tumor DNA with the tumor RNA to detect somatic mutations. The inclusion of
the RNA increases the power to detect somatic mutations, especially at low DNA
allelic frequencies. By integrating the DNA and RNA, we are able to rescue back
calls that would be missed by traditional mutation calling algorithms that only
examine the DNA.
  RADIA was developed for the identification of somatic mutations using both
DNA and RNA from the same individual. We demonstrate high sensitivity (84%) and
very high specificity (98% and 99%) in real data from endometrial carcinoma and
lung adenocarcinoma from TCGA. Mutations with both high DNA and RNA read
support have the highest validation rate of over 99%. We also introduce a
simulation package that spikes in artificial mutations to real data, rather
than simulating sequencing data from a reference genome. We evaluate
sensitivity on the simulation data and demonstrate our ability to rescue back
calls at low DNA allelic frequencies by including the RNA. Finally, we
highlight mutations in important cancer genes that were rescued back due to the
incorporation of the RNA.
  Software available at https://github.com/aradenbaugh/radia/


  The bacterial transcription factor LacI loops DNA by binding to two separate
locations on the DNA simultaneously. Despite being one of the best-studied
model systems for transcriptional regulation, the number and conformations of
loop structures accessible to LacI remain unclear, though the importance of
multiple co-existing loops has been implicated in interactions between LacI and
other cellular regulators of gene expression. To probe this issue, we have
developed a new analysis method for tethered particle motion, a versatile and
commonly-used in vitro single-molecule technique. Our method, vbTPM, performs
variational Bayesian inference in hidden Markov models. It learns the number of
distinct states (i.e., DNA-protein conformations) directly from tethered
particle motion data with better resolution than existing methods, while easily
correcting for common experimental artifacts. Studying short (roughly 100 bp)
LacI-mediated loops, we provide evidence for three distinct loop structures,
more than previously reported in single-molecule studies. Moreover, our results
confirm that changes in LacI conformation and DNA binding topology both
contribute to the repertoire of LacI-mediated loops formed in vitro, and
provide qualitatively new input for models of looping and transcriptional
regulation. We expect vbTPM to be broadly useful for probing complex
protein-nucleic acid interactions.


  We study a system of dynamical units, each of which shows excitable or
oscillatory behavior, depending on the choice of parameters. When we couple
these units with repressive bonds, we can control the duration of collective
oscillations for an intermediate period between collective fixed-point
behavior. The control mechanism works by monotonically increasing a single
bifurcation parameter. Both the onset and arrest of oscillations are due to
bifurcations. Depending on the coupling strength, the network topology and the
tuning speed, our numerical simulations reveal a rich dynamics
out-of-equilibrium with multiple inherent time scales, long transients towards
the stationary states and interesting transient patterns like self-organized
pacemakers. Zooming into the transition regime, we pursue the arrest of
oscillations along the arms of spirals. We point out possible relations to the
genetic network of segmentation clocks.


  Pasteur has demonstrated that the chiral components in a racemic mixture can
separate in homochiral crystals. But with a strong chiral discrimination the
chiral components in a concentrated mixture can also phase separate into
homochiral fluid domains, and the isomerization kinetics can then perform a
symmetry breaking into one thermodynamical stable homochiral system.
Glyceraldehyde has a sufficient chiral discrimination to perform such a
symmetry breaking. The requirement of a high concentration of the chiral
reactant(s) in an aqueous solution in order to perform and $\textit{maintain}$
homochirality; the appearance of phosphorylation of almost all carbohydrates in
the central machinery of life; the basic ideas that the biochemistry and the
glycolysis and gluconeogenesis contain the trace of the biochemical evolution,
all point in the direction of that homochirality was obtained just after- or at
a phosphorylation of the very first products of the formose reaction, at high
concentrations of the reactants in phosphate rich compartments in submarine
hydrothermal vents. A racemic solution of D,L-glyceraldehyde-3-phosphate could
be the template for obtaining homochiral D-glyceraldehyde-3-phosphate(aq) as
well as L-amino acids.


  An aspatial version for the famous Schelling's segregation model has recently
been proposed, which, called two-room model, is still in an agent-based format
like the original Schelling model. In the present study, we propose a new,
state equation version of the Schelling model. The new equation is based on the
two-room model and is derived in terms of a set of discrete maps. Fixed point
solutions for the new equation are found analytically and confirmed
numerically. Especially, we show that the extremely simple state equations can
reasonably reveal the essence of the Schelling dynamics: integration,
segregation and tipping. In addition to the fixed point solutions, periodic
solutions are identified and conditions of the limit cycles are derived
analytically.


  The minimal ecological requirements for formation of regular vegetation
patterns in semiarid systems have been recently questioned. Against the general
belief that a combination of facilitative and competitive interactions is
necessary, recent theoretical studies suggest that, under broad conditions,
nonlocal competition among plants alone may induce patterns. In this paper, we
review results along this line, presenting a series of models that yield
spatial patterns when finite-range competition is the only driving force. A
preliminary derivation of this type of model from a more detailed one that
considers water-biomass dynamics is also presented. Keywords: Vegetation
patterns, nonlocal interactions


  We derived the critical neighborhood demand in the Schelling's segregation
model by studying the conditions for which a chain reaction of migrations of
unsatisfied agents occurs. The essence of Schelling dynamics was approximated
in two simplified models: (1) a random walk model for the initial stage of the
migrations to illustrate the power-law behavior of chain reaction lengths under
critical conditions, and (2) a two-room model for the whole process to
represent a non-spatial version of segregation dynamics in the Schelling model.
Our theoretical results showed good agreements with numerical results obtained
from agent-based simulations.


  Circadian rhythms are acquired through evolution to increase the chances for
survival through synchronizing with the daylight cycle. Reliable
synchronization is realized through two trade-off properties: regularity to
keep time precisely, and entrainability to synchronize the internal time with
daylight. We found by using a phase model with multiple inputs that achieving
the maximal limit of regularity and entrainability entails many inherent
features of the circadian mechanism. At the molecular level, we demonstrate the
role sharing of two light inputs, phase advance and delay, as is well observed
in mammals. At the behavioral level, the optimal phase-response curve
inevitably contains a dead zone, a time during which light pulses neither
advance nor delay the clock. We reproduce the results of phase-controlling
experiments entrained by two types of periodic light pulses. Our results
indicate that circadian clocks are designed optimally for reliable clockwork
through evolution.


  We investigate the stochastic dynamics of entities which are confined to a
set of islands, between which they migrate. They are assumed to be one of two
types, and in addition to migration, they also reproduce and die. Systems which
fall into this class are common in biology and social science, occurring in
ecology, population genetics, epidemiology, biochemistry, linguistics, opinion
dynamics, and other areas. In all these cases the governing equations are
intractable, consisting as they do of multidimensional Fokker-Planck equations
or, equivalently, coupled nonlinear stochastic differential equations with
multiplicative noise. We develop a methodology which exploits a separation in
time scales between fast and slow variables to reduce these equations so that
they resemble those for a single island, which are amenable to analysis. The
technique is generally applicable, but we choose to discuss it in the context
of population genetics, in part because of the extra features that appear due
to selection. The idea behind the method is simple, its application systematic,
and the results in very good agreement with simulations of the full model for a
range of parameter values.


  Genotype-phenotype (GP) maps specify how the random mutations that change
genotypes generate variation by altering phenotypes, which, in turn, can
trigger selection. Many GP maps share the following general properties: 1) The
number of genotypes $N_G$ is much larger than the number of selectable
phenotypes; 2) Neutral exploration changes the variation that is accessible to
the population; 3) The distribution of phenotype frequencies $F_p=N_p/N_G$,
with $N_p$ the number of genotypes mapping onto phenotype $p$, is highly
biased: the majority of genotypes map to only a small minority of the
phenotypes. Here we explore how these properties affect the evolutionary
dynamics of haploid Wright-Fisher models that are coupled to a simplified and
general random GP map or to a more complex RNA sequence to secondary structure
map. For both maps the probability of a mutation leading to a phenotype $p$
scales to first order as $F_p$, although for the RNA map there are further
correlations as well. By using mean-field theory, supported by computer
simulations, we show that the discovery time $T_p$ of a phenotype $p$ similarly
scales to first order as $1/F_p$ for a wide range of population sizes and
mutation rates in both the monomorphic and polymorphic regimes. These
differences in the rate at which variation arises can vary over many orders of
magnitude. Phenotypic variation with a larger $F_p$ is therefore be much more
likely to arise than variation with a small $F_p$. We show, using the RNA
model, that frequent phenotypes (with larger $F_p$) can fix in a population
even when alternative, but less frequent, phenotypes with much higher fitness
are potentially accessible. In other words, if the fittest never `arrive' on
the timescales of evolutionary change, then they can't fix. We call this highly
non-ergodic effect the `arrival of the frequent'.


  The need of new diagnostic methods satisfying, as an early detection, a low
invasive procedure and a cost-efficient value, is orienting the technological
research toward the use of bio-integrated devices, in particular bio-sensors.
The set of know-why necessary to achieve this goal is wide, from biochemistry
to electronics and is summarized in an emerging branch of electronics, called
\textit{proteotronics}. Proteotronics is here here applied to state a
comparative analysis of the electrical responses coming from type-1 and type-2
opsins. In particular, the procedure is used as an early investigation of a
recently discovered family of opsins, the proteorhodopsins activated by blue
light, BPRs. The results reveal some interesting and unexpected similarities
between proteins of the two families, suggesting the global electrical response
are not strictly linked to the class identity.


  ATP-binding cassette (ABC) transporters are integral membrane proteins that
mediate the exchange of diverse substrates across membranes powered by ATP
hydrolysis. We report results of coarse-grained dynamical simulations performed
for the bacterial heme transporter HmuUV. Based on the nucleotide-free
structure, we have constructed a ligand-elastic-network description for this
protein and investigated ATP-induced conformational motions in structurally
resolved computer experiments. As we found, interactions with nucleotides
resulted in generic motions which are functional and robust. Upon binding of
ATP-mimicking ligands the structure changed from a conformation in which the
nucleotide-binding domains formed an open shape, to a conformation in which
they were found in tight contact and the transmembrane domains were rotated.
The heme channel was broadened in the ligand-bound complex and the gate to the
cytoplasm, which was closed in the nucleotide-free conformation, was rendered
open by a mechanism that involved tilting motions of essential transmembrane
helices. Based on our findings we propose that the HmuUV transporter behaves
like a `simple' mechanical device in which, induced by binding of ATP ligands,
linear motions of the nucleotide-binding domains are translated into rotational
motions and internal tilting dynamics of the transmembrane domains that control
gating inside the heme pathway.


  Within the field of phylogenetics there is great interest in distance
measures to quantify the dissimilarity of two trees. Here, based on an idea of
Bruen and Bryant, we propose and analyze a new distance measure: the Maximum
Parsimony (MP) distance. This is based on the difference of the parsimony
scores of a single character on both trees under consideration, and the goal is
to find the character which maximizes this difference. In this article we show
that this new distance is a metric and provides a lower bound to the well-known
Subtree Prune and Regraft (SPR) distance. We also show that to compute the MP
distance it is sufficient to consider only characters that are convex on one of
the trees, and prove several additional structural properties of the distance.
On the complexity side, we prove that calculating the MP distance is in general
NP-hard, and identify an interesting island of tractability in which the
distance can be calculated in polynomial time.


  We study a mosquito-borne epidemic model where the vector population is
distinct in aquatic and adult stages and a saturating effect of disease
transmission is assumed to ocurr when the number of infectious (humans and
mosquitoes) becomes large enough. Several techniques, including center manifold
analysis and sensitivity analysis, have been used to reveal relevant features
of the model dynamics. We determine the existence of stability-instability
thresholds and the individual role played in such thresholds by the model
parameters.


  Anomalous surface enhanced Raman scattering (SERS) peaks are identified for
liquid sample stored in polypropylene centrifuge tubes (PP tube) for months. We
observed the unexpected Raman peaks during experiments for Thiamine
Hydrochloride aqueous solution stored in PP tube for two months. In order to
identify the contaminants we have performed SERS experiments for de-ionized
water (DI water) stored in polypropylene centrifuge tube for two months and
compared them with fresh DI water sample. We have also carried out Ultra Violet
(UV) absorption spectra for both fresh and contaminated water. We believe that
the water is contaminated because of chemicals leaching from the PP tube. From
the GC-MS data the main contaminant was found to be Phthalic acid and its
derivatives. Further SERS and UV absorption experiment for Phthalic acid
correlates well with the anomalous peaks identified earlier. We qualitatively
confirmed the identification and quantitatively estimated the concentration of
suspect contaminants as between 1uM and 10uM with both SERS and UV absorption
spectroscopy. With UV absorption spectroscopy, we precisely estimate the
concentration as 2.1uM. We have shown that sample in PP tube can be
contaminated due to leaching chemicals upon long term storage and suggested
SERS and UV-absorption spectroscopy as two quick and simple techniques to
detect the contamination


  Genetic variation affecting gene regulation is a central driver of phenotypic
differences between individuals and can be used to uncover how biological
processes are organized in a cell. Although detecting cis-eQTLs is now routine,
trans-eQTLs have proven more challenging to find due to the modest variance
explained and the multiple tests burden of testing millions of SNPs for
association to thousands of transcripts. Here, we successfully map trans-eQTLs
with the complementary approach of looking for SNPs associated to the
expression of multiple genes simultaneously. We find 732 trans- eQTLs that
replicate across two continental populations; each trans-eQTL controls large
groups of target transcripts (regulons), which are part of interacting networks
controlled by transcription factors. We are thus able to uncover co-regulated
gene sets and begin describing the cell circuitry of gene regulation.


  Exploiting the power of the expectation operator and indicator (or Bernoulli)
random variables, we present the exact governing equations for both the SIR and
SIS epidemic models on \emph{networks}. Although SIR and SIS are basic epidemic
models, deductions from their exact stochastic equations \textbf{without}
making approximations (such as the common mean-field approximation) are scarce.
An exact analytic solution of the governing equations is highly unlikely to be
found (for any network) due to the appearing pair (and higher order)
correlations. Nevertheless, the maximum average fraction $y_{I}$ of infected
nodes in both SIS and SIR can be written as a quadratic form of the graph's
Laplacian. Only for regular graphs, the expression for the maximum of $y_{I}$
can be simplied to exhibit the explicit dependence on the spectral radius. From
our new Laplacian expression, we deduce a general \textbf{upper} bound for the
epidemic SIS threshold in any graph.


  We demonstrated fabrication of black silicon with slanted nanocone array on
both planar and 3D micro and meso scale structures produced by a
high-throughput lithography-free oblique-angle plasma etching process.
Nanocones with gradual change in height were created on the same piece of
silicon. The relation between the slanted angle of nanocones and incident angle
of directional plasma is experimentally investigated. In order to demonstrate
the monolithic integration of nanostructures on micro and meso scale non-planar
surfaces, nanocone forest is fabricated on non-planar silicon surfaces in
various morphologies such as silicon atomic force microscopy (AFM) tips and
pyramidal pits. By integrating nanocones on inverse silicon micro-pyramid array
devices, we further improved the surface enhanced Raman scattering (SERS)
enhancement property of this optimized commercial SERS substrate by several
folds even when using 66% less noble metal coating. We investigated the length
gradient dependence and asymmetric properties of SERS effects for slanted
nanocone with polarized excitation. This versatile and angle-controllable
nanocone fabrication and monolithic 3D nano-micro-meso integration method
provides new dimensions for production and optimization of SERS and other
nanophotonic sensors.


  For a beneficial allele which enters a large unstructured population and
eventually goes to fixation, it is known that the time to fixation is
approximately $2\log(\alpha)/\alpha$ for a large selection coefficient
$\alpha$. For a population that is distributed over finitely many colonies,
with migration between these colonies, we detect various regimes of the
migration rate $\mu$ for which the fixation times have different asymptotics as
$\alpha \to \infty$. If $\mu$ is of order $\alpha$, the allele fixes (as in the
spatially unstructured case) in time $\sim 2\log(\alpha)/\alpha$. If $\mu$ is
of order $\alpha^\gamma, 0\leq \gamma \leq 1$, the fixation time is $\sim (2 +
(1-\gamma)\Delta) \log(\alpha)/\alpha$, where $\Delta$ is the number of
migration steps that are needed to reach all other colonies starting from the
colony where the beneficial allele appeared. If $\mu = 1/\log(\alpha)$, the
fixation time is $\sim (2+S)\log(\alpha)/\alpha$, where $S$ is a random time in
a simple epidemic model.
  The main idea for our analysis is to combine a new moment dual for the
process conditioned to fixation with the time reversal in equilibrium of a
spatial version of Neuhauser and Krone's ancestral selection graph.