  Epigenetic Tracking is a mathematical model of biological cells, originally
conceived to study embryonic development. Computer simulations proved the
capacity of the model to generate complex 3-dimensional cellular structures,
and the potential to reproduce the complexity typical of living beings. The
most distinctive feature of this model is the presence in the body of a
homogeneous distribution of stem cells, which are dinamically and continuously
created during development from non-stem cells and reside in niches. Embryonic
stem cells orchestrate early developmental events, adult stem cells direct late
developmental and regeneration events, ageing stem cells cause ageing and
cancer stem cells are responsible for cancer growth. The conceptual backbone
provided by Epigenetic Tracking brings together a wide range of biological
phenomena: for this reason, we think it can be proposed as a general model for
multicellular biology. Despite, or perhaps due to its theoretical origin, the
model allowed us to make predictions relevant to very diverse fields of
biology, such as transposable elements, and cancer-related patterns of gene
mutations. This paper contains a summary of the model and its implications.


  In order to associate complex traits with genetic polymorphisms, genome-wide
association studies process huge datasets involving tens of thousands of
individuals genotyped for millions of polymorphisms. When handling these
datasets, which exceed the main memory of contemporary computers, one faces two
distinct challenges: 1) Millions of polymorphisms come at the cost of hundreds
of Gigabytes of genotype data, which can only be kept in secondary storage; 2)
the relatedness of the test population is represented by a covariance matrix,
which, for large populations, can only fit in the combined main memory of a
distributed architecture. In this paper, we present solutions for both
challenges: The genotype data is streamed from and to secondary storage using a
double buffering technique, while the covariance matrix is kept across the main
memory of a distributed memory system. We show that these methods sustain
high-performance and allow the analysis of enormous dataset


  Genetically identical cells in the same population can take on phenotypically
variable states, leading to differentiated responses to external signals, such
as nutrients and drug-induced stress. Many models and experiments have focused
on a description based on discrete phenotypic states. Here we consider the
effects of selection acting on a single trait, which we explicitly link to the
variable number of proteins expressed by a gene. Considering different
regulatory models for the gene under selection, we calculate the steady-state
distribution of expression levels and show how the population adapts its
expression to enhance its fitness. We quantitatively relate the overall fitness
of the population to the heritability of expression levels, and their diversity
within the population. We show how selection can increase or decrease the
variability in the population, alter the stability of bimodal states, and
impact the switching rates between metastable attractors.


  Bacteria and their viruses ("bacteriophages") coexist in natural environments
forming complex infection networks. Recent empirical findings suggest that
phage-bacteria infection networks often possess a nested structure such that
there is a hierarchical relationship among who can infect whom. Here we
consider how nested infection networks may affect phage and bacteria dynamics
using a multi-type Lotka-Volterra framework with cross-infection. Analysis of
similar models have, in the past, assumed simpler interaction structures as a
first step towards tractability. We solve the proposed model, finding trade-off
conditions on the life-history traits of both bacteria and viruses that allow
coexistence in communities with nested infection networks. First, we find that
bacterial growth rate should decrease with increasing defense against
infection. Second, we find that the efficiency of viral infection should
decrease with host range. Next, we establish a relationship between relative
densities and the curvature of life history trade-offs. We compare and contrast
the current findings to the "Kill-the-Winner" model of multi-species
phage-bacteria communities. Finally, we discuss a suite of testable hypotheses
stemming from the current model concerning relationships between infection
range, life history traits and coexistence in complex phage-bacteria
communities.


  A phase space theory for population dynamics in Ecology is presented. This
theory applies for a certain class of dynamical systems, that will be called
M-systems, for which a conserved quantity, the M-function, can be defined in
phase space. This M-function is the generator of time displacements and
contains all the dynamical information of the system. In this sense the
M-function plays the role of the hamiltonian function for mechanical systems.
In analogy with Hamilton theory we derive equations of motion as derivatives
over the resource function in phase space. A M-bracket is defined which allows
one to perform a geometrical approach in analogy to Poisson bracket of
hamiltonian systems. We show that the equations of motion can be derived from a
variational principle over a functional J of the trajectories. This functional
plays for M-systems the same role than the action S for hamiltonian systems.
Finally, three important systems in population dynamics, namely,
Lotka-Volterra, self-feeding and logistic evolution, are shown to be M-systems.


  We investigate, using the Hierarchy method, the entanglement and the
excitation transfer efficiency of the Fenna-Matthews-Olson complex under two
different local modifications: the suppression of transitions between
particular sites and localized changes to the protein environment. We find that
inhibiting the connection between the site-5 and site-6, or disconnecting
site-5 from the complex completely, leads to an dramatic enhancement of the
entanglement between site-6 and site-7. Similarly, the transfer efficiency
actually increases if site-5 is disconnected from the complex entirely. We
further show that if site-5 and site-7 are conjointly removed, the efficiency
falls. This suggests that while not contributing to the transport efficiency in
a normal complex, site-5 introduces a redundant transport route in case of
damage to site-7. Our results suggest an overall robustness of excitation
energy transfer in the FMO complex under mutations, local defects, and other
abnormal situations.


  The character diagnosis of Tenontosaurus tilletti has been revised and
redefined into a more robust and quantifiable state. Significant emphasis is
placed on constructing phylogenetic definition in such a method, as it prevents
occlusion of true character states by alleviating potential individual
interpretational bias. Previous placement within the Iguanodontia is refuted
based on the lack of character affinity with the defining synapomorphies of the
clade. The clade Hypsilophodontidae (=Hypsilophodontia), along with
Iguanodontia, however is deemed to be in critical need of refinement to account
for recent discoveries and re-classifications of certain euornithopods. Several
of the synapomorphies are out-dated and deemed redundant in favour of a more
quantifiable approach. Re- definition of these clades is critical if the
current state of basal euornithopodan relationships is to be resolved.
Phylogenetic studies must be approached from a multidisciplinary perspective;
integration of tectonostratigraphical, ontogenetic, palaeoecological, and
biomechanical data with sets of well-defined primary homologies are essential
in increasing phylogenetic resolution and generating stratigraphically feasible
ancestor-descendant relationships. Material attributed to Tenontosaurus
tilletti is in need of strict re-analysis; the significant quantity of
specimens attributed to this species is potentially the result of poor
stratigraphic constraints and the vast spatiotemporal span occupied. Future
revision of this material is expected to reveal temporal variations on the
species -level inherently linked to environmental evolution, as well as
possibly provide clues to sexual dimorphism in contemporaneous, yet
morphologically distinct tenontosaurs.


  In finite-size population models, one can derive Fokker-Planck equations to
describe the fluctuations of the species numbers about the deterministic
behaviour. In the steady state of populations comprising two or more species,
it is permissible for a probability current to flow. In such a case, the system
does not relax to equilibrium but instead reaches a non-equilibrium steady
state. In a two-species model, these currents form cycles (e.g., ellipses) in
probability space. We investigate the conditions under which such currents are
solely responsible for macroscopically-observable cycles in population
abundances. We find that this can be achieved when the deterministic limit
yields a circular neutrally-stable manifold. We further discuss the efficacy of
one-dimensional approximation to the diffusion on the manifold, and obtain
estimates for the macroscopically-observable current around this manifold by
appealing to Kramers' escape rate theory.


  Genetic linkage may result in the expression of multiple products from a
polycistronic transcript, under the control of a single promoter. In animals,
protein-coding polycistronic transcripts are rare. However, microRNAs are
frequently clustered in the genomes of animals, and these clusters are often
transcribed as a single unit. The evolution of microRNA clusters has been the
subject of much speculation, and a selective advantage of clusters of
functionally related microRNAs is often proposed. However, the origin of
microRNA clusters has not been so far explored. Here we study the evolution of
microRNA clusters in Drosophila melanogaster. We observed that the majority of
microRNA clusters arose by the de novo formation of new microRNA-like hairpins
in existing microRNA transcripts. Some clusters also emerged by tandem
duplication of a single microRNA. Comparative genomics show that these clusters
are unlikely to split or undergo rearrangements. We did not find any instances
of clusters appearing by rearrangement of pre-existing microRNA genes. We
propose a model for microRNA cluster evolution in which selection over one of
the microRNAs in the cluster interferes with the evolution of the other linked
microRNAs. Our analysis suggests that the study of microRNAs and small RNAs
must consider linkage associations.


  We establish the fundamental limits of DNA shotgun sequencing under noisy
reads. We show a surprising result: for the i.i.d. DNA model, noisy reads are
as good as noiseless reads, provided that the noise level is below a certain
threshold which can be surprisingly high. As an example, for a uniformly
distributed DNA sequence and a symmetric substitution noisy read channel, the
threshold is as high as 19%.


  Sex in higher diploids carries a two-fold cost of males that should reduce
its fitness relative to cloning and result in its extinction. Instead, sex is
widespread and it is clonal species that face early obsolescence. One possible
reason is that sex is an adaptation to resist ubiquitous parasites, which
evolve rapidly and potentially antagonistically. We use a heuristic approach to
model mutation-selection in finite populations where a parasitic haploid mounts
a negative frequency-dependent attack on a diploid host. The host evolves
reflexively to reduce parasitic load. Both host and parasite populations
generate novel alleles by mutation and have access to large allele spaces. Sex
outcompetes cloning by two overlapping mechanisms. First, sexual diploids adopt
advantageous homozygous mutations more rapidly than clonal diploids under
conditions of lag load. This rate advantage can offset the lesser fecundity of
sex. Second, a relative advantage to sex emerges under host mutation rates that
are fast enough to retain fitness in a rapidly mutating parasite environment
and increase host polymorphism. Clonal polymorphic populations
disproportionately experience interference with selection at high mutation
rates, both between and within loci. This slows clonal population adaptation to
a changing parasite environment and reduces clonal population fitness relative
to sex. The interference increases markedly with the number of loci under
independent selection. Rates of parasite mutation exist that not only allow sex
to survive despite the two-fold cost of males but which enable sexual and
clonal populations to have equal fitness and co-exist. Since all higher
organisms carry parasitic loads, the model is of general applicability.


  The phenomenon of apoptosis observed in cell cultures consists in the fact
that unless cells permanently receive a "Stay alive" signal from their
neighbors, they are bound to die. A natural question is whether manifestations
of this apoptosis paradigm can also be observed in other organizations of
living organisms. In this paper we report results from a two-year long campaign
of experiments on three species of ants and one species of (tephritid) fruit
flies. In these experiments individuals were separated from their colony and
kept in isolation either alone or in groups of 10 individuals. The overall
conclusion is that "singles" have a shorter life expectancy than individuals in
the groups of 10. This observation holds for ants as well as for fruit flies.
The paper also provides compelling evidence of a similar effect in married
versus unmarried (i.e. single, widowed or divorced) people. A natural question
concerns the dynamic of the transition between the two regimes. Observation
suggests an abrupt (rather than smooth) transition and this conclusion seems to
hold for ants, fruit flies and humans as well. We call it a shock transition.
In addition, for red fire ants Solenopsis invicta, it was observed that
individuals in groups of 10 that also comprise one queen, die much faster than
those in similar groups without queens. The paper also examines the
corresponding survivorship curves from the perspective of the standard
classification into 3 types. The survivorship curves of ants (whether single or
in groups of 10) are found to be of type II whereas those of the fruit fly
Bactrocera dorsalis are rather of type III. In this connection it is recalled
that the survivorship curve of the fruit fly Drosophila melanogaster is of type
I, i.e. of same type as for humans.


  In a population practicing consanguineous marriage, rare recessive lethal
alleles (RRLA) have higher chances of affecting phenotypes. As inbreeding
causes more homozygosity and subsequently more deaths, the loss of individuals
with RRLA decreases the frequency of these alleles. Although this phenomenon is
well studied in general, here some hitherto unstudied cases are presented. An
analytical formula for the RRLA frequency is presented for infinite monoecious
population practicing several different types of inbreeding. In finite diecious
populations, it is found that more severe inbreeding leads to quicker RRLA
losses, making the upcoming generations healthier. A population of size 10,000
practicing 30% half-sib marriages loses more than 95% of its RRLA in 100
generations; a population practicing 30% cousin marriages loses about 75% of
its RRLA. Our findings also suggest that given enough resources to grow, a
small inbred population will be able to rebound while losing the RRLA.


  We have created a genome-scale network reconstruction of Escherichia coli
metabolism. Existing reconstructions were improved in terms of annotation
standards, to facilitate their subsequent use in dynamic modelling. The
resultant network is available from EcoliNet (http://ecoli.sf.net/).


  Herein we provide a closed form perturbative solution to a general $M$-node
network SIS model using the transport rates between nodes as a perturbation
parameter. We separate the dynamics into a short-time regime and a
medium/long-time regime. We solve the short-time dynamics of the system and
provide a limit before which our explicit, analytical result of the first-order
perturbation for the medium/long-time regime is to be employed. These stitched
calculations provide an approximation to the full temporal dynamics for rather
general initial conditions.
  To further corroborate our results, we solve the mean-field equations
numerically for an infectious SIS outbreak in New Zealand (NZ, \emph{Aotearoa})
recomposed into 23 subpopulations where the virus is spread to different
subpopulations via (documented) air traffic data, and the country is
internationally quarantined. We demonstrate that our analytical predictions
compare well to the numerical solution.


  We have created a genome-scale network reconstruction of chinese hamster
ovary (CHO) cell metabolism. Existing reconstructions were improved in terms of
annotation standards, to facilitate their subsequent use in dynamic modelling.
The resultant network is available from ChoNet (http://cho.sf.net/).


  We investigate the consequences of adopting the criteria used by the state of
California, as described by Myers et al. (2011), for conducting familial
searches. We carried out a simulation study of randomly generated profiles of
related and unrelated individuals with 13-locus CODIS genotypes and YFiler
Y-chromosome haplotypes, on which the Myers protocol for relative
identification was carried out. For Y-chromosome sharing first degree
relatives, the Myers protocol has a high probability (80 - 99%) of identifying
their relationship. For unrelated individuals, there is a low probability that
an unrelated person in the database will be identified as a first-degree
relative. For more distant Y-haplotype sharing relatives (half-siblings, first
cousins, half-first cousins or second cousins) there is a substantial
probability that the more distant relative will be incorrectly identified as a
first-degree relative. For example, there is a 3 - 18% probability that a first
cousin will be identified as a full sibling, with the probability depending on
the population background. Although the California familial search policy is
likely to identify a first degree relative if his profile is in the database,
and it poses little risk of falsely identifying an unrelated individual in a
database as a first-degree relative, there is a substantial risk of falsely
identifying a more distant Y-haplotype sharing relative in the database as a
first-degree relative, with the consequence that their immediate family may
become the target for further investigation. This risk falls disproportionately
on those ethnic groups that are currently overrepresented in state and federal
databases.


  Ecological networks are often composed of different sub-communities (often
referred to as modules). Identifying such modules has the potential to develop
a better understanding of the assembly of ecological communities and to
investigate functional overlap or specialisation. The most informative form of
networks are quantitative or weighted networks. Here we introduce an algorithm
to identify modules in quantitative bipartite (or two-mode) networks. It is
based on the hierarchical random graphs concept of Clauset et al. (2008 Nature
453: 98-101) and is extended to include quantitative information and adapted to
work with bipartite graphs. We define the algorithm, which we call QuaBiMo,
sketch its performance on simulated data and illustrate its potential
usefulness with a case study.


  Y haplogroup analyses are an important component of genealogical
reconstruction, population genetic analyses, medical genetics and forensics.
These fields are increasingly moving towards use of low-coverage, high
throughput sequencing. However, there is as yet no software available for using
sequence data to assign Y haplogroup groups probabilistically, such that the
posterior probability of assignment fully reflects the information present in
the data, and borrows information across all samples sequenced from a
population. YHap addresses this problem.


  Discrete dynamic models are a powerful tool for the understanding and
modeling of large biological networks. Although a lot of progress has been made
in developing analysis tools for these models, there is still a need to find
approaches that can directly relate the network structure to its dynamics. Of
special interest is identifying the stable patterns of activity, i.e., the
attractors of the system. This is a problem for large networks, because the
state space of the system increases exponentially with network size. In this
work we present a novel network reduction approach that is based on finding
network motifs that stabilize in a fixed state. Notably, we use a topological
criterion to identify these motifs. Specifically, we find certain types of
strongly connected components in a suitably expanded representation of the
network. To test our method we apply it to a dynamic network model for a type
of cytotoxic T cell cancer and to an ensemble of random Boolean networks of
size up to 200. Our results show that our method goes beyond reducing the
network and in most cases can actually predict the dynamical repertoire of the
nodes (fixed states or oscillations) in the attractors of the system.


  Animals that scavenge in and around human localities need to utilize a broad
range of resources. Preference for any one kind of food, under such
circumstances, might be inefficient. Indian free-ranging dogs, Canis lupus
familiaris are scavengers that are heavily dependent on humans for sustaining
their omnivorous diet. The current study suggests that because of evolutionary
load, these dogs, which are descendants of the decidedly carnivorous gray wolf,
still retain a preference for meat though they live on carbohydrate-rich
resources. The plasticity in their diet probably fosters efficient scavenging
in a competitive environment, while a thumb rule for preferentially acquiring
specific nutrients enables them to sequester proteins from the
carbohydrate-rich environment.


  We present a nonparametric Bayesian method for disease subtype discovery in
multi-dimensional cancer data. Our method can simultaneously analyse a wide
range of data types, allowing for both agreement and disagreement between their
underlying clustering structure. It includes feature selection and infers the
most likely number of disease subtypes, given the data.
  We apply the method to 277 glioblastoma samples from The Cancer Genome Atlas,
for which there are gene expression, copy number variation, methylation and
microRNA data. We identify 8 distinct consensus subtypes and study their
prognostic value for death, new tumour events, progression and recurrence. The
consensus subtypes are prognostic of tumour recurrence (log-rank p-value of
$3.6 \times 10^{-4}$ after correction for multiple hypothesis tests). This is
driven principally by the methylation data (log-rank p-value of $2.0 \times
10^{-3}$) but the effect is strengthened by the other 3 data types,
demonstrating the value of integrating multiple data types.
  Of particular note is a subtype of 47 patients characterised by very low
levels of methylation. This subtype has very low rates of tumour recurrence and
no new events in 10 years of follow up. We also identify a small gene
expression subtype of 6 patients that shows particularly poor survival
outcomes. Additionally, we note a consensus subtype that showly a highly
distinctive data signature and suggest that it is therefore a biologically
distinct subtype of glioblastoma.
  The code is available from https://sites.google.com/site/multipledatafusion/


  Influenza viruses undergo continual antigenic evolution allowing mutant
viruses to evade host immunity acquired to previous virus strains. Antigenic
phenotype is often assessed through pairwise measurement of cross-reactivity
between influenza strains using the hemagglutination inhibition (HI) assay.
Here, we extend previous approaches to antigenic cartography, and
simultaneously characterize antigenic and genetic evolution by modeling the
diffusion of antigenic phenotype over a shared virus phylogeny. Using HI data
from influenza lineages A/H3N2, A/H1N1, B/Victoria and B/Yamagata, we determine
patterns of antigenic drift across viral lineages, showing that A/H3N2 evolves
faster and in a more punctuated fashion than other influenza lineages. We also
show that year-to-year antigenic drift appears to drive incidence patterns
within each influenza lineage. This work makes possible substantial future
advances in investigating the dynamics of influenza and other
antigenically-variable pathogens by providing a model that intimately combines
molecular and antigenic evolution.


  We develop a path-based approach to continuous-time random walks on networks
with arbitrarily weighted edges. We describe an efficient numerical algorithm
for calculating statistical properties of the stochastic path ensemble. After
demonstrating our approach on two reaction rate problems, we present a
biophysical model that describes how proteins evolve new functions while
maintaining thermodynamic stability. We use our methodology to characterize
dynamics of evolutionary adaptation, reproducing several key features observed
in directed evolution experiments. We find that proteins generally fall into
two qualitatively different regimes of adaptation depending on their binding
and folding energetics.


  Identifiability of evolutionary tree models has been a recent topic of
discussion and some models have been shown to be non-identifiable. A
coalescent-based rooted population tree model, originally proposed by Nielsen
et al. 1998 [2], has been used by many authors in the last few years and is a
simple tool to accurately model the changes in allele frequencies in the tree.
However, the identifiability of this model has never been proven. Here we prove
this model to be identifiable by showing that the model parameters can be
expressed as functions of the probability distributions of subsamples. This a
step toward proving the consistency of the maximum likelihood estimator of the
population tree based on this model.


  The concept of a nested canalizing Boolean function has been studied over the
course of the last decade in the context of understanding the regulatory logic
of molecular interaction networks, such as gene regulatory networks. Such
functions appear preferentially in published models of such networks. Recently,
this concept has been generalized to include multi-state functions, and a
recursive formula has been derived for their number, as a function of the
number of variables. This paper carries out a detailed analysis of the class of
nested canalizing functions over an arbitrary finite field. Furthermore, the
paper generalizes the concept further, and derives a closed formula for the
number of such generalized functions. The paper also derives a closed formula
for the number of equivalence classes under permutation of variables. This is
motivated by the fact that two nested canalizing functions that differ by a
permutation of the variables share many important properties with each other.
The paper contributes to the effort of identifying a class of functions over
finite fields that are of interest in biology and also have interesting
mathematical properties.


  The Joint United Nations Programme on HIV and AIDS (UNAIDS) has set a new
90-90-90 global target for the coverage of anti-retroviral therapy (ART) to be
reached by 2020. This would mean that 90% of all people infected with HIV know
their status, 90% of them are on ART and 90% of them will have full viral load
suppression. Here we first estimate the case reproduction number, R0, for
countries in sub-Saharan Africa and for India using data on the rate at which
the prevalence of HIV increased at the start of the epidemic and the life
expectancy of people living with HIV who are not on ART. R0 determines the
magnitude of the control problem, that is to say, the extent to which
transmission must be reduced to eliminate HIV. We show that in sub-Saharan
Africa the median value of R0 is 4.6 and in all but five countries R0 is less
than 6.3. If the 90-90-90 target is reached, 73% of all those living with HIV
will have full viral load suppression. If this is maintained it should
guarantee elimination in 70% of all countries in sub-Saharan Africa and will
reduce R0 to less than 2 in the remaining 12 countries, making elimination easy
to achieve by increasing the availability of other high impact methods of
prevention.


  The factors that influence genetic architecture shape the structure of the
fitness landscape, and therefore play a large role in the evolutionary
dynamics. Here the NK model is used to investigate how epistasis and pleiotropy
-- key components of genetic architecture -- affect the structure of the
fitness landscape, and how they affect the ability of evolving populations to
adapt despite the difficulty of crossing valleys present in rugged landscapes.
Populations are seen to make use of epistatic interactions and pleiotropy to
attain higher fitness, and are not inhibited by the fact that valleys have to
be crossed to reach peaks of higher fitness.


  Background: The short reads output by first- and second-generation DNA
sequencing instruments cannot completely reconstruct microbial chromosomes.
Therefore, most genomes have been left unfinished due to the significant
resources required to manually close gaps in draft assemblies.
Third-generation, single-molecule sequencing addresses this problem by greatly
increasing sequencing read length, which simplifies the assembly problem.
  Results: To measure the benefit of single-molecule sequencing on microbial
genome assembly, we sequenced and assembled the genomes of six bacteria and
analyzed the repeat complexity of 2,267 complete bacteria and archaea. Our
results indicate that the majority of known bacterial and archaeal genomes can
be assembled without gaps, at finished-grade quality, using a single PacBio RS
sequencing library. These single-library assemblies are also more accurate than
typical short-read assemblies and hybrid assemblies of short and long reads.
  Conclusions: Automated assembly of long, single-molecule sequencing data
reduces the cost of microbial finishing to $1,000 for most genomes, and future
advances in this technology are expected to drive the cost lower. This is
expected to increase the number of completed genomes, improve the quality of
microbial genome databases, and enable high-fidelity, population-scale studies
of pan-genomes and chromosomal organization.


  We investigate a special part of the human immune system, namely the
activation of T-Cells, using stochastic tools, especially sharp large deviation
results. T-Cells have to distinguish reliably between foreign and self peptides
which are both presented to them by antigen presenting cells. Our work is based
on a model studied by Zint, Baake, and den Hollander, and originally proposed
by van den Berg, Rand, and Burroughs. We are able to dispense with some
restrictive distribution assumptions that were used previously, i.e. we
establish a higher robustness of the model. A central issue is the analysis of
two new perspectives to the scenario (two different quenched systems) in
detail. This means that we do not only analyse the total probability of a
T-Cell activation (the annealed case) but also consider the probability of an
activation of one certain T-Cell type and the probability of a T-Cell
activation by a certain antigen presenting cell (the quenched cases). Finally,
we see analytically that the probability of T-Cell activation increases with
the number of presented foreign peptides in all three cases.


  Cooperation played a significant role in the self-organization and evolution
of living organisms. Both network topology and the initial position of
cooperators heavily affect the cooperation of social dilemma games. We
developed a novel simulation program package, called 'NetworGame', which is
able to simulate any type of social dilemma games on any model, or real world
networks with any assignment of initial cooperation or defection strategies to
network nodes. The ability of initially defecting single nodes to break overall
cooperation was called as 'game centrality'. The efficiency of this measure was
verified on well-known social networks, and was extended to 'protein games',
i.e. the simulation of cooperation between proteins, or their amino acids. Hubs
and in particular, party hubs of yeast protein-protein interaction networks had
a large influence to convert the cooperation of other nodes to defection.
Simulations on methionyl-tRNA synthetase protein structure network indicated an
increased influence of nodes belonging to intra-protein signaling pathways on
breaking cooperation. The efficiency of single, initially defecting nodes to
convert the cooperation of other nodes to defection in social dilemma games may
be an important measure to predict the importance of nodes in the integration
and regulation of complex systems. Game centrality may help to design more
efficient interventions to cellular networks (in forms of drugs), to ecosystems
and social networks. The NetworGame algorithm is downloadable from here:
www.NetworGame.linkgroup.hu


  Propulsion by growing actin networks is a universal mechanism used in many
different biological systems. Although the core molecular machinery for actin
network growth is well preserved in most cases, the geometry of the propelled
obstacle can vary considerably. In recent years, filament orientation
distribution has emerged as an important observable characterizing the
structure and dynamical state of the growing network. Here we derive several
continuum equations for the orientation distribution of filaments growing
behind stiff obstacles of various shapes and validate the predicted steady
state orientation patterns by stochastic computer simulations based on discrete
filaments. We use an ordinary differential equation approach to demonstrate
that for flat obstacles of finite size, two fundamentally different orientation
patterns peaked at either +35/-35 or +70/0/-70 degrees exhibit mutually
exclusive stability, in agreement with earlier results for flat obstacles of
very large lateral extension. We calculate and validate phase diagrams as a
function of model parameters and show how this approach can be extended to
obstacles with piecewise straight contours. For curved obstacles, we arrive at
a partial differential equation in the continuum limit, which again is in good
agreement with the computer simulations. In all cases, we can identify the same
two fundamentally different orientation patterns, but only within an
appropriate reference frame, which is adjusted to the local orientation of the
obstacle contour. Our results suggest that two fundamentally different network
architectures compete with each other in growing actin networks, irrespective
of obstacle geometry, and clarify how simulated and electron tomography data
have to be analyzed for non-flat obstacle geometries.


  DNA microarrays are a relatively new technology that can simultaneously
measure the expression level of thousands of genes. They have become an
important tool for a wide variety of biological experiments. One of the most
common goals of DNA microarray experiments is to identify genes associated with
biological processes of interest. Conventional statistical tests often produce
poor results when applied to microarray data due to small sample sizes, noisy
data, and correlation among the expression levels of the genes. Thus, novel
statistical methods are needed to identify significant genes in DNA microarray
experiments. This article discusses the challenges inherent in DNA microarray
analysis and describes a series of statistical techniques that can be used to
overcome these challenges. The problem of multiple hypothesis testing and its
relation to microarray studies is also considered, along with several possible
solutions.


  We report the first realization of a biomolecular AND gate function with
double-sigmoid response (sigmoid in both inputs). Two enzyme biomarker inputs
activate the gate output signal which can then be used as indicating liver
injury, but only when both of these inputs have elevated pathophysiological
concentrations, effectively corresponding to logic-1 of the binary gate
functioning. At lower, normal physiological concentrations, defined as logic-0
inputs, the liver-injury output levels are not obtained. High-quality gate
functioning in handling of various sources of noise, on time scales of
relevance to potential applications is enabled by utilizing "filtering"
effected by a simple added biocatalytic process. The resulting gate response is
sigmoid in both inputs when proper system parameters are chosen, and the gate
properties are theoretically analyzed within a model devised to evaluate its
noise-handling properties.


  Stochastic modeling and simulation provide powerful predictive methods for
the intrinsic understanding of fundamental mechanisms in complex biochemical
networks. Typically, such mathematical models involve networks of coupled jump
stochastic processes with a large number of parameters that need to be suitably
calibrated against experimental data. In this direction, the parameter
sensitivity analysis of reaction networks is an essential mathematical and
computational tool, yielding information regarding the robustness and the
identifiability of model parameters. However, existing sensitivity analysis
approaches such as variants of the finite difference method can have an
overwhelming computational cost in models with a high-dimensional parameter
space. We develop a sensitivity analysis methodology suitable for complex
stochastic reaction networks with a large number of parameters. The proposed
approach is based on Information Theory methods and relies on the
quantification of information loss due to parameter perturbations between
time-series distributions. For this reason, we need to work on path-space,
i.e., the set consisting of all stochastic trajectories, hence the proposed
approach is referred to as "pathwise". The pathwise sensitivity analysis method
is realized by employing the rigorously-derived Relative Entropy Rate (RER),
which is directly computable from the propensity functions. A key aspect of the
method is that an associated pathwise Fisher Information Matrix (FIM) is
defined, which in turn constitutes a gradient-free approach to quantifying
parameter sensitivities. The structure of the FIM turns out to be
block-diagonal, revealing hidden parameter dependencies and sensitivities in
reaction networks.


  Chemoreceptors McpB and McpC in Salmonella enterica have been reported to
promote chemotaxis in LB motility-plate assays. Of the chemicals tested as
potential effectors of these receptors, the only response was towards
L-cysteine and its oxidized form, L-cystine. Although enhanced radial migration
in plates suggested positive chemotaxis to both amino acids, capillary assays
failed to show an attractant response to either, in cells expressing only these
two chemoreceptors. In vivo fluorescence resonance energy transfer (FRET)
measurements of kinase activity revealed that in wild-type bacteria, cysteine
and cystine are chemoeffectors of opposing sign, the reduced form being a
chemoattractant and the oxidized form a repellent. The attractant response to
cysteine was mediated primarily by Tsr, as reported earlier for E. coli. The
repellent response to cystine was mediated by McpB / C. Adaptive recovery upon
cystine exposure required the methyl-transferase/-esterase pair, CheR / CheB,
but restoration of kinase activity was never complete (i.e. imperfect
adaptation). We provide a plausible explanation for the attractant-like
responses to both cystine and cysteine motility plates, and speculate that the
opposing signs of response to this redox pair might afford Salmonella a
mechanism to gauge and avoid oxidative environments.


  We have considered as a theoretical possibility for the development of a
nanobiochip the operation principle of which is based on measuring conductance
in single-stranded and double-stranded DNA. Calculations have demonstrated that
in the majority of cases the conductance of double-stranded nucleotides
considerably exceeds that of single-stranded ones. The results obtained are in
agreement with recent experiments on measuring the oligonucleotide conductance.
It has been shown that an electronic biochip containing 11 nucleotide pairs
will recognize approximately 97% sequences. It has also been demonstrated that
the percentage of identifiable sequences will grow with the sequence length.


  We report the phenomenon of vibrational resonance in a single species and a
two species models of groundwater-dependent plant ecosystems with a biharmonic
oscillation (with two widely different frequencies \omega and \Omega, \Omega >>
\omega) of the water table depth. In these two systems, the response amplitude
of the species biomass shows multiple resonances with different mechanisms. The
resonance occurs at both low- and high-frequencies of the biharmonic force. In
the single species bistable system, the resonance occurs at discrete values of
the amplitude g of the high-frequency component of the water table.
Furthermore, the best synchronization of biomass and its carrying capacity with
the biharmonic force occurs at the resonance. In the two species excitable and
time-delay model, the response amplitude (Q) profile shows several plateau
regions of resonance, where the period of evolution of the species biomass
remains the same and the value of Q is inversely proportional to it. The
response amplitude is highly sensitive to the time-delay parameter \tau and
shows two distinct sequences of resonance intervals with a decreasing amplitude
with \tau.


  We present a thermodynamically robust coarse-grained model to simulate
folding of RNA in monovalent salt solutions. The model includes stacking,
hydrogen bond and electrostatic interactions as fundamental components in
describing the stability of RNA structures. The stacking interactions are
parametrized using a set of nucleotide-specific parameters, which were
calibrated against the thermodynamic measurements for single-base stacks and
base-pair stacks. All hydrogen bonds are assumed to have the same strength,
regardless of their context in the RNA structure. The ionic buffer is modeled
implicitly, using the concept of counterion condensation and the Debye-H\"uckel
theory. The three adjustable parameters in the model were determined by fitting
the experimental data for two RNA hairpins and a pseudoknot. A single set of
parameters provides good agreement with thermodynamic data for the three RNA
molecules over a wide range of temperatures and salt concentrations. In the
process of calibrating the model, we establish the extent of counterion
condensation onto the single-stranded RNA backbone. The reduced backbone charge
is independent of the ionic strength and is 60% of the RNA bare charge at 37
degrees Celsius. Our model can be used to predict the folding thermodynamics
for any RNA molecule in the presence of monovalent ions.


  We review the effects of macromolecular crowding on the folding of RNA by
considering the simplest scenario when excluded volume interactions between
crowding particles and RNA dominate. Using human telomerase enzyme as an
example, we discuss how crowding can alter the equilibrium between pseudoknot
and hairpin states of the same RNA molecule - a key aspect of crowder-RNA
interactions. We summarize data showing that the crowding effect is significant
only if the size of the spherical crowding particle is smaller than the radius
of gyration of the RNA in the absence of crowding particles. The implication
for function of the wild type and mutants of human telomerase is outlined by
using a relationship between enzyme activity and its conformational
equilibrium. In addition, we discuss the interplay between macromolecular
crowding and ionic strength of the RNA buffer. Finally, we briefly review
recent experiments which illustrate the connection between excluded volume due
to macromolecular crowding and the thermodynamics of RNA folding.


  Stimulation of LDH initial reaction rates by timed pre-irradiation of
crystalline sodium pyruvate and lithium lactate is reported for enzymes
isolated from rabbit muscle, pig heart, human erythrocytes and chicken heart.
The phenomenon investigated is referred to as the Comorosan effect. For the
mammalian source enzymes, the pyruvate irradiation stimulations occurred at
irradiation times of 5 and 35 sec. and the lactate irradiation times at 15 and
45 seconds. In contrast, for the chicken heart enzyme, the pyruvate irradiation
stimulations occurred at 15 and 35 sec., while those for lactate occurred at 5
and 20 sec. Thus, a shift in stimulatory irradiation times is found on going
from the mammalian enzymes to the avian enzyme. A similar shift between
mammalian and yeast enzymes has been established by Comorosan and co-workers.
For the chicken heart LDH, the separation between successive irradiation times
is different for the forward and reverse reactions. This is the first reported
incidence of the separation not being the same.


  Infection dynamics have been studied extensively on complex networks,
yielding insight into the effects of heterogeneity in contact patterns on
disease spread. Somewhat separately, metapopulations have provided a paradigm
for modeling systems with spatially extended and "patchy" organization. In this
paper we expand on the use of multitype networks for combining these paradigms,
such that simple contagion models can include complexity in the agent
interactions and multiscale structure. We first present a generalization of the
Volz-Miller mean-field approximation for Susceptible-Infected-Recovered (SIR)
dynamics on multitype networks. We then use this technique to study the special
case of epidemic fronts propagating on a one-dimensional lattice of
interconnected networks - representing a simple chain of coupled population
centers - as a necessary first step in understanding how macro-scale disease
spread depends on micro-scale topology. Using the formalism of front
propagation into unstable states, we derive the effective transport
coefficients of the linear spreading: asymptotic speed, characteristic
wavelength, and diffusion coefficient for the leading edge of the pulled
fronts, and analyze their dependence on the underlying graph structure. We also
derive the epidemic threshold for the system and study the front profile for
various network configurations. To our knowledge, this is the first such
application of front propagation concepts to random network models.


  Under the effect of strong genetic drift, it is highly probable to observe
gene fixation or gene loss in a population, shown by infinite peaks on a
coherently constructed potential energy landscape. It is then important to ask
what such singular peaks imply, with or without the effects of other biological
factors. We studied the stochastic escape time from the infinite potential
peaks in the Wright-Fisher model, where the typical two-scale diffusion
dynamics was observed via computer simulations. We numerically found the
average escape time for all the bi-stable cases and analytically approximated
the results under weak mutations and selections by calculating the mean first
passage time (MFPT) in singular potential peak. Our results showed that
Kramers' classical escape formula can be extended to the models with
non-Gaussian probability distributions, overcoming constraints in previous
methods. The constructed landscape provides a global and coherent description
for system's evolutionary dynamics, allowing new biological results to be
generated.


  RNA molecules follow a succession of enzyme-mediated processing steps from
transcription until maturation. The participating enzymes, for example the
spliceosome for mRNAs and Drosha and Dicer for microRNAs, are also produced in
the cell and their copy-numbers fluctuate over time. Enzyme copy-number changes
affect the processing rate of the substrate molecules; high enzyme numbers
increase the processing probability, low enzyme numbers decrease it. We study
different RNA processing cascades where enzyme copy-numbers are either fixed or
fluctuate. We find that for fixed enzyme-copy numbers the substrates at
steady-state are Poisson-distributed, and the whole RNA cascade dynamics can be
understood as a single birth-death process of the mature RNA product. In this
case, solely fluctuations in the timing of RNA processing lead to variation in
the number of RNA molecules. However, we show analytically and numerically that
when enzyme copy-numbers fluctuate, the strength of RNA fluctuations increases
linearly with the RNA transcription rate. This linear effect becomes stronger
as the speed of enzyme dynamics decreases relative to the speed of RNA
dynamics. Interestingly, we find that under certain conditions, the RNA cascade
can reduce the strength of fluctuations in the expression level of the mature
RNA product. Finally, by investigating the effects of processing polymorphisms
we show that it is possible for the effects of transcriptional polymorphisms to
be enhanced, reduced, or even reversed. Our results provide a framework to
understand the dynamics of RNA processing.


  Background Oryzias latipes (Medaka) has been established as a vertebrate
genetic model for over a century, and has recently been rediscovered outside
its native Japan. The power of new sequencing methods now makes it possible to
reinvigorate Medaka genetics, in particular by establishing a near-isogenic
panel derived from a single wild population. Results Here we characterise the
genomes of wild Medaka catches obtained from a single Southern Japanese
population in Kiyosu as a precursor for the establishment of a near isogenic
panel of wild lines. The population is free of significant detrimental
population structure, and has advantageous linkage disequilibrium properties
suitable for establishment of the proposed panel. Analysis of morphometric
traits in five representative inbred strains suggests phenotypic mapping will
be feasible in the panel. In addition high throughput genome sequencing of
these Medaka strains confirms their evolutionary relationships on lines of
geographic separation and provides further evidence that there has been little
significant interbreeding between the Southern and Northern Medaka population
since the Southern/Northern population split. The sequence data suggest that
the Southern Japanese Medaka existed as a larger older population which went
through a relatively recent bottleneck around 10,000 years ago. In addition we
detect patterns of recent positive selection in the Southern population.
Conclusions These data indicate that the genetic structure of the Kiyosu Medaka
samples are suitable for the establishment of a vertebrate near isogenic panel
and therefore inbreeding of 200 lines based on this population has commenced.
Progress of this project can be tracked at
http://www.ebi.ac.uk/birney-srv/medaka-ref-panel


  This paper answers an open question of Chen, Doty, and Soloveichik [1], who
showed that a function f:N^k --> N^l is deterministically computable by a
stochastic chemical reaction network (CRN) if and only if the graph of f is a
semilinear subset of N^{k+l}. That construction crucially used "leaders": the
ability to start in an initial configuration with constant but non-zero counts
of species other than the k species X_1,...,X_k representing the input to the
function f. The authors asked whether deterministic CRNs without a leader
retain the same power.
  We answer this question affirmatively, showing that every semilinear function
is deterministically computable by a CRN whose initial configuration contains
only the input species X_1,...,X_k, and zero counts of every other species. We
show that this CRN completes in expected time O(n), where n is the total number
of input molecules. This time bound is slower than the O(log^5 n) achieved in
[1], but faster than the O(n log n) achieved by the direct construction of [1]
(Theorem 4.1 in the latest online version of [1]), since the fast construction
of that paper (Theorem 4.4) relied heavily on the use of a fast, error-prone
CRN that computes arbitrary computable functions, and which crucially uses a
leader.


  Background: Coevolution within a protein family is often predicted using
statistics that measure the degree of covariation between positions in the
protein sequence. Mutual Information is a measure of dependence between two
random variables that has been used extensively to predict intra-protein
coevolution.
  Results: Here we provide an algorithm for the efficient calculation of Mutual
Information within a protein family. The algorithm uses linked lists which are
directly accessed by a pointer array. The linked list allows efficient storage
of sparse count data caused by protein conservation. The direct access array of
pointers prevents the linked list from being traversed each time it is
modified.
  Conclusions: This algorithm is implemented in the software MIpToolset, but
could also be easily implemented in other Mutual Information based standalone
software or web servers. The current implementation in the MIpToolset has been
critical in large-scale protein family analysis and real-time coevolution
calculations during alignment editing and curation.
  The MIpToolset is available at: https://sourceforge.net/projects/miptoolset/


  Background: Computational analysis of next-generation sequencing data is
outpaced by data generation in many cases. In one such case, paired-end reads
can be produced from the Illumina sequencing method faster than they can be
overlapped by downstream analysis. The advantages in read length and accuracy
provided by overlapping paired-end reads demonstrates the necessity for
software to efficiently solve this problem.
  Results: XORRO is an extremely efficient paired-end read overlapping program.
XORRO can overlap millions of short paired-end reads in a few minutes. It uses
64-bit registers with a two bit alphabet to represent sequences and does
comparisons using low-level logical operations like XOR, AND, bitshifting and
popcount.
  Conclusions: As of the writing of this manuscript, XORRO provides the fastest
solution to the paired-end read overlap problem. XORRO is available for
download at: sourceforge.net/projects/xorro-overlap/


  Canids display a wide diversity of social systems, from solitary to pairs to
packs, and hence they have been extensively used as model systems to understand
social dynamics in natural systems. Among canids, the dog can show various
levels of social organization due to the influence of humans on their lives.
Though the dog is known as man's best friend and has been studied extensively
as a pet, studies on the natural history, ecology and behaviour of dogs in a
natural habitat are rare. Here we report results of an extensive
population-level study conducted through one-time censuses in urban India to
understand the ecoethology of free-ranging dogs. We built a model to test if
the observed groups could have been formed through random associations while
foraging. Our modeling results suggest that the dogs, like all efficient
scavengers, tend to forage singly but also form random uncorrelated groups. A
closer inspection of the group compositions however reveals that the foraging
associations are non-random events. The tendency of adults to associate with
the opposite sex in the mating season and of juveniles to stay close to adults
in the non-mating season drives the population towards aggregation, in spite of
the apparently random nature of the group size distribution. Hence we conclude
that to be or not to be social is a matter of choice for the free-ranging dogs,
and not a matter of chance.


  Neurons are thought of as the building blocks of excitable brain tissue.
However, at the single neuron level, the neuronal membrane, the dendritic arbor
and the axonal projections can also be considered an extended active medium.
Active dendritic branchlets enable the propagation of dendritic spikes, whose
computational functions, despite several proposals, remain an open question.
Here we propose a concrete function to the active channels in large dendritic
trees. By using a probabilistic cellular automaton approach, we model the
input-output response of large active dendritic arbors subjected to complex
spatio-temporal inputs and exhibiting non-stereotyped dendritic spikes. We find
that, if dendritic spikes have a non-deterministic duration, the dendritic
arbor can undergo a continuous phase transition from a quiescent to an active
state, thereby exhibiting spontaneous and self-sustained localized activity as
suggested by experiments. Analogously to the critical brain hypothesis, which
states that neuronal networks self-organize near a phase transition to take
advantage of specific properties of the critical state, here we propose that
neurons with large dendritic arbors optimize their capacity to distinguish
incoming stimuli at the critical state. We suggest that "computation at the
edge of a phase transition" is more compatible with the view that dendritic
arbors perform an analog rather than a digital dendritic computation.


  Epidemic models currently play a central role in our attempts to understand
and control infectious diseases. Here, we derive a model for the diffusion
limit of stochastic susceptible-infectious-removed (SIR) epidemic dynamics on a
heterogeneous network. Using this, we consider analytically the early
asymptotic exponential growth phase of such epidemics, showing how the higher
order moments of the network degree distribution enter into the stochastic
behaviour of the epidemic. We find that the first three moments of the network
degree distribution are needed to specify the variance in disease prevalence
fully, meaning that the skewness of the degree distribution affects the
variance of the prevalence of infection. We compare these asymptotic results to
simulation and find a close agreement for city-sized populations.


  Summary: We present an extension of CUSHAW2 for fast and accurate alignments
of SOLiD color-space short-reads. Our extension introduces a double-seeding
approach to improve mapping sensitivity, by combining maximal exact match seeds
and variable-length seeds derived from local alignments. We have compared the
performance of CUSHAW2 to SHRiMP2 and BFAST by aligning both simulated and real
color-space mate-paired reads to the human genome. The results show that
CUSHAW2 achieves comparable or better alignment quality compared to SHRiMP2 and
BFAST at an order-of-magnitude faster speed and significantly smaller peak
resident memory size. Availability: CUSHAW2 and all simulated datasets are
available at http://cushaw2.sourceforge.net. Contact: liuy@uni-mainz.de;
bertil.schmidt@uni-mainz.de


  The expected level of linkage disequilibrium (LD) in a finite ideal
population at equilibrium is of relevance for many applications in population
and quantitative genetics. Several recursion formulae have been proposed during
the last decades, whose derivations mostly contain heuristic parts and
therefore remain mathematically questionable. We propose a more justifiable
approach, including an alternative recursion formula for the expected LD. Since
the exact formula depends on the distribution of allele frequencies in a very
complicated manner, we suggest an approximate solution and analyze its validity
extensively in a simulation study. Compared to the widely used formula of Sved,
the proposed formula performs better for all parameter constellations
considered. We then analyze the expected LD at equilibrium using the theory on
discrete-time Markov chains based on the linear recursion formula, with
equilibrium being defined as the steady-state of the chain, which finally leads
to a formula for the effective population size N_e. An additional analysis
considers the effect of non-exactness of a recursion formula on the
steady-state, demonstrating that the resulting error in expected LD can be
substantial. In an application to the HapMap data of two human populations we
illustrate the dependency of the N_e-estimate on the distribution of minor
allele frequencies (MAFs), showing that the estimated N_e can vary by up to 30%
when a uniform instead of a skewed distribution of MAFs is taken as a basis to
select SNPs for the analyses. Our analyses provide new insights into the
mathematical complexity of the problem studied.


  Modern DNA sequencing technologies enable geneticists to rapidly identify
genetic variation among many human genomes. However, isolating the minority of
variants underlying disease remains an important, yet formidable challenge for
medical genetics. We have developed GEMINI (GEnome MINIng), a flexible software
package for exploring all forms of human genetic variation. Unlike existing
tools, GEMINI integrates genetic variation with a diverse and flexible set of
genome annotations (e.g., dbSNP, ENCODE, UCSC, ClinVar, KEGG) into a unified
database to facilitate interpretation and data exploration. Whereas other
methods provide an inflexible set of variant filters or variant prioritization
methods, GEMINI allows researchers to compose complex queries based on sample
genotypes, inheritance patterns, and both pre-installed and custom genome
annotations. GEMINI also provides methods for ad hoc queries and data
exploration, a simple programming interface for custom analyses that leverage
the underlying database, and both command line and graphical tools for common
analyses. We demonstrate the utility of GEMINI for exploring variation in
personal genomes and family based genetic studies, and illustrate its ability
to scale to studies involving thousands of human samples. GEMINI is designed
for reproducibility and flexibility and our goal is to will provide researchers
with a standard framework for medical genomics.


  Punishment may deter antisocial behavior. Yet to punish is costly, and the
costs often do not offset the gains that are due to elevated levels of
cooperation. However, the effectiveness of punishment depends not only on how
costly it is, but also on the circumstances defining the social dilemma. Using
the snowdrift game as the basis, we have conducted a series of economic
experiments to determine whether severe punishment is more effective than mild
punishment. We have observed that severe punishment is not necessarily more
effective, even if the cost of punishment is identical in both cases. The
benefits of severe punishment become evident only under extremely adverse
conditions, when to cooperate is highly improbable in the absence of sanctions.
If cooperation is likely, mild punishment is not less effective and leads to
higher average payoffs, and is thus the much preferred alternative. Presented
results suggest that the positive effects of punishment stem not only from
imposed fines, but may also have a psychological background. Small fines can do
wonders in motivating us to chose cooperation over defection, but without the
paralyzing effect that may be brought about by large fines. The later should be
utilized only when absolutely necessary.


  A modeling formalism is proposed for the description and study of living and
life-like systems. It provides an abstract conceptual model framework for real
life and evolution of biological organisms. It is proposed, that this model
formalism provides a novel system view and immediately applicable conceptual
tools for understanding complex life and evolution. The modeling principle is
very generic suggesting that it can be directly applied also to the study of
engineered and artificial systems.


  In this study extremely halophilic and moderately thermophilic microorganisms
from a hypersaline microbial mat were screened for their ability to produce
antibacterial, antidiatom, antialgal and quorum sensing (QS) inhibitory
compounds. Five bacterial strains belonging to the genera Marinobacter and
Halomonas and one archaeal strain belonging to the genus Haloterrigena were
isolated from a microbial mat. The strains were able to grow at a maximum
salinity of 22-25% and a maximum temperature of 45-60{\deg}C. Hexanes,
dichloromethane and butanol extracts from the strains inhibited the growth of
at least one out of nine human pathogens. Only butanol extracts of supernatants
of Halomonas sp. SK-1 inhibited growth of the microalga Dunaliella salina. Most
extracts from isolates inhibited QS of the acyl homoserine lactone producer and
reporter Chromobacterium violaceum CV017. Purification of QS inhibitory
dichloromethane extracts of Marinobacter sp. SK-3 resulted in isolation of four
related diketopiperazines (DKPs): cyclo(L-Pro-L-Phe), cyclo(L-Pro-L-Leu),
cyclo(L-Pro-L-isoLeu) and cyclo(L-Pro-D-Phe). QS inhibitory properties of these
DKPs were tested using C. violaceum CV017 and Escherichia coli-based QS
reporters (pSB401 and pSB1075) deficient in AHL production. Cyclo(L-Pro-L-Phe)
and cyclo(L-Pro-L-isoLeu) inhibited QS dependent production of violacein by C.
violaceum CV017. Cyclo(L-Pro-L-Phe), cyclo(L-Pro-L-Leu), and
cyclo(L-Pro-L-isoLeu) reduced QS dependent luminescence of the reporter E. coli
pSB401 induced by 3-oxo-C6-HSL. Our study demonstrated the ability of
halophilic and moderately thermophilic strains from a hypersaline microbial mat
produce biotechnologically-relevant compounds that could be used as antifouling
agents


  Reaction networks are systems in which the populations of a finite number of
species evolve through predefined interactions. Such networks are found as
modeling tools in many biological disciplines such as biochemistry, ecology,
epidemiology, immunology, systems biology and synthetic biology. It is now
well-established that, for small population sizes, stochastic models for
biochemical reaction networks are necessary to capture randomness in the
interactions. The tools for analyzing such models, however, still lag far
behind their deterministic counterparts. In this paper, we bridge this gap by
developing a constructive framework for examining the long-term behavior and
stability properties of the reaction dynamics in a stochastic setting. In
particular, we address the problems of determining ergodicity of the reaction
dynamics, which is analogous to having a globally attracting fixed point for
deterministic dynamics. We also examine when the statistical moments of the
underlying process remain bounded with time and when they converge to their
steady state values. The framework we develop relies on a blend of ideas from
probability theory, linear algebra and optimization theory. We demonstrate that
the stability properties of a wide class of biological networks can be assessed
from our sufficient theoretical conditions that can be recast as efficient and
scalable linear programs, well-known for their tractability. It is notably
shown that the computational complexity is often linear in the number of
species. We illustrate the validity, the efficiency and the wide applicability
of our results on several reaction networks arising in biochemistry, systems
biology, epidemiology and ecology. The biological implications of the results
as well as an example of a non-ergodic biological network are also discussed.


  Expansion of the neocortex is a hallmark of human evolution. However, it
remains an open question what adaptive mechanisms facilitated its expansion.
Here we show, using gyrencephaly index (GI) and other physiological and
life-history data for 102 mammalian species, that gyrencephaly is an ancestral
mammalian trait. We provide evidence that the evolution of a highly folded
neocortex, as observed in humans, requires the traversal of a threshold of 10^9
neurons, and that species above and below the threshold exhibit a bimodal
distribution of physiological and life-history traits, establishing two
phenotypic groups. We identify, using discrete mathematical models,
proliferative divisions of progenitors in the basal compartment of the
developing neocortex as evolutionarily necessary and sufficient for generating
a fourteen-fold increase in daily prenatal neuron production and thus traversal
of the neuronal threshold. We demonstrate that length of neurogenic period,
rather than any novel progenitor-type, is sufficient to distinguish cortical
neuron number between species within the same phenotypic group.


  Charge transfer is considered in systems composed of a donor, an acceptor and
bridge sites of (AT) nucleotide pairs. For a bridge consisting of 180 (AT)
pairs, three cases are dealt with: a uniform case, when all the nucleotides in
each strand are identical; an ordered case, when nucleotides in each DNA strand
are arranged in an orderly fashion; a chaotic case, when (AT) and (TA) pairs
are arranged randomly. It is shown that in all the cases a charge transfer from
a donor to an acceptor can take place. All other factors being equal, the
transfer is the most efficient in the uniform case, the ordered and chaotic
cases are less and the least efficient, accordingly. The results obtained are
in agreement with experimental data on long-range charge transfer in DNA.


  Numerous experiments on charge transfer in DNA yield a contradictory picture
of the transfer: on the one hand they suggest that it is a very slow process
and the charge is almost completely localized on one Watson-Crick pair, but on
the other hand they demonstrate that the charge can travel a very large
distance. To explain this contradiction we propose that superfast charge
transitions are possible between base pairs on individual DNA fragments
resulting in the establishment of a quasi-equilibrium charge distribution
during the time less than that of charge solvation. In other words, we
hypothesize these states irrespective of the nature of a mechanism responsible
for their establishment, whether it be a hopping mechanism, or a band
mechanism, or superexchange, or polaron transport, etc., leaving aside the
debates of which one is more advantageous. We discuss qualitative differences
between the charge transfer in a dry DNA and that in a solution. In a solution,
of great importance is the charge solvation which decreases the transfer rate
$10^7 - 10^8$ times as compared with a dry DNA. We consider the conditions
under which the superfast charge transfer in a DNA leading to quasi-equilibrium
distributions of polarons in a duplex is possible. Comparison of calculated
quasi-equilibrium distributions with the experiment testifies to the
possibility of superfast tunnel transitions of a hole in a DNA duplex in a
solution.


  In this paper, we attempt to provide mathematical models of Mendelian and
Non-Mendelian inheritances of the bisexual population system having Fisher's
{\textbf{1:1}} principle. In our model, we always assume that distributions of
the same phenotype of female and male populations are equal. We study the
evolution of a Mendelian trait. As an application of a non-Mendelian
inheritance, we construct a quadratic stochastic operator that describes
transmission of {\textbf{ABO}} and Rh blood groups.


  One of the outstanding challenges in comparative genomics is to interpret the
evolutionary importance of regulatory variation between species. Rigorous
molecular evolution-based methods to infer evidence for natural selection from
expression data are at a premium in the field, and to date, phylogenetic
approaches have not been well-suited to address the question in the small sets
of taxa profiled in standard surveys of gene expression. We have developed a
strategy to infer evolutionary histories from expression profiles by analyzing
suites of genes of common function. In a manner conceptually similar to
molecular evolution models in which the evolutionary rates of DNA sequence at
multiple loci follow a gamma distribution, we modeled expression of the genes
of an \emph{a priori}-defined pathway with rates drawn from an inverse gamma
distribution. We then developed a fitting strategy to infer the parameters of
this distribution from expression measurements, and to identify gene groups
whose expression patterns were consistent with evolutionary constraint or rapid
evolution in particular species. Simulations confirmed the power and accuracy
of our inference method. As an experimental testbed for our approach, we
generated and analyzed transcriptional profiles of four \emph{Saccharomyces}
yeasts. The results revealed pathways with signatures of constrained and
accelerated regulatory evolution in individual yeasts and across the phylogeny,
highlighting the prevalence of pathway-level expression change during the
divergence of yeast species. We anticipate that our pathway-based phylogenetic
approach will be of broad utility in the search to understand the evolutionary
relevance of regulatory change.


  The two-regime method (TRM) has been recently developed for optimizing
stochastic reaction-diffusion simulations. It is a multiscale (hybrid)
algorithm which uses stochastic reaction-diffusion models with different levels
of detail in different parts of the computational domain. The coupling
condition on the interface between different modelling regimes of the TRM was
previously derived for one-dimensional models. In this paper, the TRM is
generalized to higher dimensional reaction-diffusion systems. Coupling Brownian
dynamics models with compartment-based models on regular (square)
two-dimensional lattices is studied in detail. In this case, the interface
between different modelling regimes contain either flat parts or right-angled
corners. Both cases are studied in the paper. For flat interfaces, it is shown
that the one-dimensional theory can be used along the line perpendicular to the
TRM interface. In the direction tangential to the interface, two choices of the
TRM parameters are presented. Their applicability depends on the compartment
size and the time step used in the molecular-based regime. The two-dimensional
generalization of the TRM is also discussed in the case of corners.


  The new dynamical game theoretic model of sex ratio evolution emphasizes the
role of males as passive carriers of sex ratio genes. This shows inconsistency
between population genetic models of sex ratio evolution and classical
strategic models. In this work a novel technique of change of coordinates will
be applied to the new model. This will reveal new aspects of the modelled
phenomenon which cannot be shown or proven in the original formulation. The
underlying goal is to describe the dynamics of selection of particular genes in
the entire population, instead of in the same sex subpopulation, as in the
previous paper and earlier population genetics approaches. This allows for
analytical derivation of the unbiased strategic model from the model with
rigorous non-simplified genetics. In effect, an alternative system of
replicator equations is derived. It contains two subsystems: the first
describes changes in gene frequencies (this is an alternative unbiased
formalization of the Fisher-Dusing argument), whereas the second describes
changes in the sex ratios in subpopulations of carriers of genes for each
strategy. An intriguing analytical result of this work is that fitness of a
gene depends on the current sex ratio in the subpopulation of its carriers, not
on the encoded individual strategy. Thus, the argument of the gene fitness
function is not constant but is determined by the trajectory of the sex ratio
among carriers of that gene. This aspect of the modelled phenomenon cannot be
revealed by the static analysis. Dynamics of the sex ratio among gene carriers
is driven by a dynamic "tug of war" between female carriers expressing the
encoded strategic trait value and random partners of male carriers expressing
the average population strategy (a primary sex ratio). This mechanism can be
called "double level selection". Therefore, gene interest perspective leads to
multi-level selection.


  Popular sequence alignment tools such as BWA convert a reference genome to an
indexing data structure based on the Burrows-Wheeler Transform (BWT), from
which matches to individual query sequences can be rapidly determined. However
the utility of also indexing the query sequences themselves remains relatively
unexplored.
  Here we show that an all-against-all comparison of two sequence collections
can be computed from the BWT of each collection with the BWTs held entirely in
external memory, i.e. on disk and not in RAM. As an application of this
technique, we show that BWTs of transcriptomic and genomic reads can be
compared to obtain reference-free predictions of splice junctions that have
high overlap with results from more standard reference-based methods.
  Code to construct and compare the BWT of large genomic data sets is available
at http://beetl.github.com/BEETL/ as part of the BEETL library.


  Understanding large molecular networks consisting of entities such as genes,
proteins or RNAs that interact in complex ways to drive the cellular machinery
has been an active focus of systems biology. Computational approaches have
played a key role in systems biology by complementing theoretical and
experimental approaches. Here we roadmap some key contributions of
computational methods developed over the last decade in the reconstruction of
biological pathways. We position these contributions in a 'systems biology
perspective' to reemphasize their roles in unraveling cellular mechanisms and
to understand 'systems biology diseases' including cancer.


  A model based on a thermodynamic approach is proposed for predicting the
dynamics of communicable epidemics in a city, when the epidemic is governed by
controlling efforts of multiple scales so that an entropy is associated with
the system. All the epidemic details are factored into a single parameter that
is determined by maximizing the rate of entropy production. Despite the
simplicity of the final model, it predicts the number of hospitalized cases
with a reasonable accuracy, using the data of SARS of the year 2003, once the
inflexion point characterizing the effect of multiple controlling efforts is
known. This model is supposed to be of potential usefulness since epidemics
such as avian influenza like H7H9 in China this year have the risk to become
communicable among human beings.


  Coevolution of two species is typically thought to favour the evolution of
faster evolutionary rates helping a species keep ahead in the Red Queen race,
where `it takes all the running you can do to stay where you are'. In contrast,
if species are in a mutualistic relationship, it was proposed that the Red King
effect may act, where it can be beneficial to evolve slower than the
mutualistic species. The Red King hypothesis proposes that the species which
evolves slower can gain a larger share of the benefits. However, the
interactions between the two species may involve multiple individuals. To
analyse such a situation, we resort to evolutionary multiplayer games. Even in
situations where evolving slower is beneficial in a two-player setting, faster
evolution may be favoured in a multiplayer setting. The underlying features of
multiplayer games can be crucial for the distribution of benefits. They also
suggest a link between the evolution of the rate of evolution and group size.


  Genome assembly tools based on the de Bruijn graph framework rely on a
parameter k, which represents a trade-off between several competing effects
that are difficult to quantify. There is currently a lack of tools that would
automatically estimate the best k to use and/or quickly generate histograms of
k-mer abundances that would allow the user to make an informed decision.
  We develop a fast and accurate sampling method that constructs approximate
abundance histograms with a several orders of magnitude performance improvement
over traditional methods. We then present a fast heuristic that uses the
generated abundance histograms for putative k values to estimate the best
possible value of k. We test the effectiveness of our tool using diverse
sequencing datasets and find that its choice of k leads to some of the best
assemblies.
  Our tool KmerGenie is freely available at: http://kmergenie.bx.psu.edu/


  We report a realization of an associative memory signal/information
processing system based on simple enzyme-catalyzed biochemical reactions.
Optically detected chemical output is always obtained in response to the
triggering input, but the system can also "learn" by association, to later
respond to the second input if it is initially applied in combination with the
triggering input as the "training" step. This second chemical input is not
self-reinforcing in the present system, which therefore can later "unlearn" to
react to the second input if it is applied several times on its own. Such
processing steps realized with (bio)chemical kinetics promise applications of
bio-inspired/memory-involving components in "networked" (concatenated)
biomolecular processes for multi-signal sensing and complex information
processing.


  Evolutionary events such as incomplete lineage sorting and lateral gene
transfer constitute major problems for inferring species trees from gene trees,
as they can sometimes lead to gene trees which conflict with the underlying
species tree. One particularly simple and efficient way to infer species trees
from gene trees under such conditions is to combine three-taxon analyses for
several genes using a majority vote approach. For incomplete lineage sorting
this method is known to be statistically consistent, however, in the case of
lateral gene transfer it is known that a zone of inconsistency does exist for a
specific four-taxon tree topology. In this paper we analyze all remaining
four-taxon topologies and show that no other inconsistencies exist.


  The dynamic behaviors of microRNA and mRNA under external stress are studied
with biological experiments and mathematics models. In this study, we developed
a mathematic model to describe the biological phenomenon and for the first time
reported that, as responses to external stress, the expression levels of
microRNA and mRNA sustained oscillation. And the period of the oscillation is
much shorter than several reported transcriptional regulation negative feedback
loop.


  Ductal carcinoma {\em in situ} (DCIS) lesions are non-invasive tumours of the
breast which are thought to precede most invasive breast cancers (IBC). As
individual DCIS lesions are initiated, grow and invade (i.e. become IBC) the
size distribution of the DCIS lesions present in a given human population will
evolve. We derive a differential equation governing this evolution and show,
for given assumptions about growth and invasion, that there is a unique
distribution which does not vary with time. Further, we show that any initial
distribution converges to this stationary distribution exponentially quickly.
It is therefore reasonable to assume that the stationary distribution is equal
to the true DCIS size distribution, at least for human populations which are
relatively stable with respect to the determinants of breast cancer. Based on
this assumption and the size data of 110 DCIS lesions detected in a
mammographic screening program between 1993 and 2000, we produce maximum
likelihood estimates for certain growth and invasion parameters. Assuming that
DCIS size is proportional to a positive power $p$ of the time since tumour
initiation we estimate $p$ to be 0.50 with a 95% confidence interval of $(0.35,
0.71)$. Therefore we estimate that DCIS lesions follow a square-root growth law
and hence that they grow rapidly when small and relatively slowly when large.
Our approach and results should be useful for other mathematical studies of
cancer, especially those investigating biological mechanisms of invasion.


  The Yule process generates a class of binary trees which is fundamental to
population genetic models and other applications in evolutionary biology. In
this paper, we introduce a family of sub-classes of ranked trees, called
Omega-trees, which are characterized by imbalance of internal nodes. The degree
of imbalance is defined by an integer 0 <= w. For caterpillars, the extreme
case of unbalanced trees, w = 0. Under models of neutral evolution, for
instance the Yule model, trees with small w are unlikely to occur by chance.
Indeed, imbalance can be a signature of permanent selection pressure, such as
observable in the genealogies of certain pathogens. From a mathematical point
of view it is interesting to observe that the space of Omega-trees maintains
several statistical invariants although it is drastically reduced in size
compared to the space of unconstrained Yule trees. Using generating functions,
we study here some basic combinatorial properties of Omega-trees. We focus on
the distribution of the number of subtrees with two leaves. We show that
expectation and variance of this distribution match those for unconstrained
trees already for very small values of w.


  So far, more than 82,000 protein structures have been reported in the Protein
Data Bank, but the driving force and structures that allow for protein
functions have not been elucidated at the atomic level for even one protein. We
have been able to clarify that the inter-subunit hydrophobic interaction
driving the electrostatic opening of the pore in aquaporin 0 (AQP0). Aquaporins
are membrane channels for water and small non-ionic solutes found in animals,
plants, and microbes. The structures of aquaporins have high homology and
consist of homotetramers, each monomer of which has one pore for a water
channel. Each pore has two narrow portions: one is the narrowest constriction
region consisting of aromatic residues and an arginine (ar/R), and another is
two asparagine-proline-alanine (NPA) homolog portions. Here we show that an
inter-subunit hydrophobic interaction in AQP0 drives a stick portion consisting
of four amino acids toward the pore and the tip of the stick portion,
consisting of a nitrogen atom, opens the pore: that movement is the swing
mechanism (this http URL). The energetics and conformational change of amino
acids participating in the swing mechanism confirm this view. The swing
mechanism in which inter-subunit hydrophobic interactions in the tetramer drive
the on-off switching of the pore explains why aquaporins consist of tetramers.
Here, we report that experimental and molecular dynamics findings using various
mutants support this view of the swing mechanism. The finding that mutants of
amino acids in AQP2 corresponding to the stick of the swing mechanism cause
severe recessive nephrogenic diabetes insipidus (NDI) demonstrates the critical
role of the swing mechanism for the aquaporin function. We report first that
the inter-subunit hydrophobic interaction in aquaporin 0 drives the
electrostatic opening of the aquaporin pore at the atomic level.


  Recent advances in computational power and simulation programs finally
delivered the first examples of reversible folding for small proteins with an
all-atom description. But having at hand the atomistic details of the process
did not lead to a straightforward interpretation of the mechanism. For the case
of the Fip35 WW-domain where multiple long trajectories of 100 {\mu}s are
available from D. E. Shaw Research, different interpretations emerged. Some of
those are in clear contradiction with each other while others are in
qualitative agreement. Here, we present a network-based analysis of the same
data by looking at the local fluctuations of conventional order parameters for
folding. We found that folding occurs through two major pathways, one almost
four times more populated than the other. Each pathway involves the formation
of an intermediate with one of the two hairpins in a native configuration. The
quantitative agreement of our results with a state-of-the-art reaction
coordinate optimization procedure as well as qualitative agreement with other
Markov-state-models and different simulation schemes provides strong evidence
for a multiple folding pathways scenario with the presence of intermediates.


  Specific interactions between cations and proteins have a strong impact on
peptide and protein structure. We here shed light on the nature of the
underlying interactions, especially regarding the effects on the polyamide
backbone structure. To do so, we compare the conformational ensembles of model
peptides in isolation and in the presence of either Li+ or Na+ cations by
state-of-the-art density-functional theory (including van der Waals effects)
and gas-phase infrared spectroscopy. These monovalent cations have a drastic
effect on the local backbone conformation of turn-forming peptides, by
disruption of the H bonding networks and the resulting severe distortion of the
backbone conformations. In fact, Li+ and Na+ can even have different
conformational effects on the same peptide. We also assess the predictive power
of current approximate density functionals for peptide-cation systems and
compare to results from established protein force fields as well as to
high-level quantum chemistry (CCSD(T)).


  The development of novel high-throughput sequencing (HTS) methods for RNA
(RNA-Seq) has provided a very powerful mean to study splicing under multiple
conditions at unprecedented depth. However, the complexity of the information
to be analyzed has turned this into a challenging task. In the last few years,
a plethora of tools have been developed, allowing researchers to process
RNA-Seq data to study the expression of isoforms and splicing events, and their
relative changes under different conditions. We provide an overview of the
methods available to study splicing from short RNA-Seq data. We group the
methods according to the different questions they address: 1) Assignment of the
sequencing reads to their likely gene of origin. This is addressed by methods
that map reads to the genome and/or to the available gene annotations. 2)
Recovering the sequence of splicing events and isoforms. This is addressed by
transcript reconstruction and de novo assembly methods. 3) Quantification of
events and isoforms. Either after reconstructing transcripts or using an
annotation, many methods estimate the expression level or the relative usage of
isoforms and/or events. 4) Providing an isoform or event view of differential
splicing or expression. These include methods that compare relative
event/isoform abundance or isoform expression across two or more conditions. 5)
Visualizing splicing regulation. Various tools facilitate the visualization of
the RNA-Seq data in the context of alternative splicing. In this review, we do
not describe the specific mathematical models behind each method. Our aim is
rather to provide an overview that could serve as an entry point for users who
need to decide on a suitable tool for a specific analysis. We also attempt to
propose a classification of the tools according to the operations they do, to
facilitate the comparison and choice of methods.


  The dynamics leading to extinction or coexistence of competing species is of
great interest in ecology and related fields. Recently a model of intra- and
interspecific competition between two species was proposed by Gabel et al.
[Phys. Rev. E 87 (2013) 010101], in which the scarcer species (i.e., with
smaller stationary population size) can be more resistant to extinction when it
holds a competitive advantage; the latter study considered populations without
spatial variation. Here we verify this phenomenon in populations distributed in
space. We extend the model of Gabel et al. to a d-dimensional lattice, and
study its population dynamics both analytically and numerically. Survival of
the scarcer in space is verified for situations in which the more competitive
species is closer to the threshold for extinction than is the less competitive
species, when considered in isolation. The conditions for survival of the
scarcer species, as obtained applying renormalization group analysis and Monte
Carlo simulation, differ in detail from those found in the spatially
homogeneous case. Simulations highlight the speed of invasion waves in
determining the survival times of the competing species.


  We introduce an atomistic approach to the dissipative quantum dynamics of
charged or neutral excitations propagating through macromolecular systems.
Using the Feynman-Vernon path integral formalism, we analytically trace out
from the density matrix the atomic coordinates and the heat bath degrees of
freedom. This way we obtain an effective field theory which describes the
real-time evolution of the quantum excitation and is fully consistent with the
fluctuation-dissipation relation. The main advantage of the field-theoretic
approach is that it allows to avoid using the Keldysh contour formulation. This
simplification makes it straightforward to derive Feynman diagrams to
analytically compute the effects of the interaction of the propagating quantum
excitation with the heat bath and with the molecular atomic vibrations. For
illustration purposes, we apply this formalism to investigate the loss of
quantum coherence of holes propagating through a poly(3-alkylthiophene) polymer


  Sequence alignment is one of the oldest and the most famous problems in
bioinformatics. Even after 45 years, for one reason or another, this problem is
still actual; current solutions are trade-offs between execution time, memory
consumption and accuracy. We purpose SW#, a new CUDA GPU enabled and memory
efficient implementation of dynamic programming algorithms for local alignment.
In this implementation indels are treated using the affine gap model. Although
there are other GPU implementations of the Smith-Waterman algorithm, SW# is the
only publicly available implementation that can produce sequence alignments on
genome-wide scale. For long sequences, our implementation is at least a few
hundred times faster than a CPU version of the same algorithm.


  Mendez et al. recently report the identification of a Y chromosome lineage
from an African American that is an outgroup to all other known Y haplotypes,
and report a time to most recent common ancestor, TMRCA, for human Y lineages
that is substantially longer than any previous estimate. The identification of
a novel Y haplotype is always exciting, and this haplotype, in particular, is
unique in its basal position on the Y haplotype tree. However, at 338 (237-581)
thousand years ago, kya, the extremely ancient TMRCA reported by Mendez et al.
is inconsistent with the known human fossil record (which estimate the age of
anatomically modern humans at 195 +- 5 kya), with estimates from mtDNA (176.6
+- 11.3 kya, and 204.9 (116.8-295.7) kya) and with population genetic theory.
The inflated TMRCA can quite easily be attributed to the extremely low Y
chromosome mutation rate used by the authors.


  We study a stochastic community model able to interpolate from a neutral
regime to a niche partitioned regime upon varying a single parameter tuning the
intensity of niche stabilization, namely the difference between intraspecific
and interspecific competition. By means of a self-consistent approach, we
obtain an analytical expression for the species abundance distribution, in
excellent agreement with stochastic simulations of the model. In the neutral
limit, the Fisher log-series is recovered, while upon increasing the
stabilization strength the species abundance distribution develops a maximum
for species at intermediate abundances, corresponding to the emergence of a
carrying capacity. Numerical studies of species extinction-time distribution
show that niche-stabilization strongly affects also the dynamical properties of
the system by increasing the average species lifetimes, while suppressing their
fluctuations. The results are discussed in view of the niche-neutral debate and
of their potential relevance to field data.


  The dispersal stages of organisms with sessile adults must be able to select
habitats with suitable conditions for establishment and survival, and must also
be able to reach those locations. For marine planktonic larvae, movement due to
currents is often orders of magnitude greater than movement due to swimming
behaviour, so transport is largely passive.
  Current patterns are determined by the interaction of geography with tidal
forces, modified by meteorological conditions. These linkages impose an
area-specific focus to connectivity studies. Yet, how geographical features and
meteorological forcing combine to produce specific current patterns and
resultant connectivity among populations remains unclear.
  In this study, we followed tracks of generic particles driven by modelled
hydrodynamic currents to investigate how connectivity between evenly spaced
habitat sites varies in relation to coastal topography. Larvae released from
regions of high current velocity, open coastline and low local habitat
availability travelled furthest but were less likely to disperse successfully.
Extensive natal habitat in the vicinity of a site generally had a positive
impact on the number of larvae arriving, as did low current velocities.
However, relationships between numbers of arriving larvae and local
geographical indices were complex, particularly at longer larval durations.
  Local geography alone explained up to 50% of the variance in success of
larvae released and around 10% of the variation in the number of larvae
arriving at each site, but coastline properties fall short of predicting
dispersal measures for particular locations. The study shows that
meteorological variation and broad scale current patterns interact strongly
with local geography to determine connectivity in coastal areas.


  We propose a simple model for genetic adaptation to a changing environment,
describing a fitness landscape characterized by two maxima. One is associated
with "specialist" individuals that are adapted to the environment; this maximum
moves over time as the environment changes. The other maximum is static, and
represents "generalist" individuals not affected by environmental changes. The
rest of the landscape is occupied by "maladapted" individuals. Our analysis
considers the evolution of these three subpopulations. Our main result is that,
in presence of a sufficiently stable environmental feature, as in the case of
an unchanging aspect of a physical habitat, specialists can dominate the
population. By contrast, rapidly changing environmental features, such as
language or cultural habits, are a moving target for the genes; here,
generalists dominate, because the best evolutionary strategy is to adopt
neutral alleles not specialized for any specific environment. The model we
propose is based on simple assumptions about evolutionary dynamics and
describes all possible scenarios in a non-trivial phase diagram. The approach
provides a general framework to address such fundamental issues as the Baldwin
effect, the biological basis for language, or the ecological consequences of a
rapid climate change.