  A DNA polymerase (DNAP) replicates a template DNA strand. It also exploits
the template as the track for its own motor-like mechanical movement. In the
polymerase mode it elongates the nascent DNA by one nucleotide in each step.
But, whenever it commits an error by misincorporating an incorrect nucleotide,
it can switch to an exonuclease mode. In the latter mode it excises the wrong
nucleotide before switching back to its polymerase mode. We develop a
stochastic kinetic model of DNA replication that mimics an {\it in-vitro}
experiment where a single-stranded DNA, subjected to a mechanical tension $F$,
is converted to a double-stranded DNA by a single DNAP. The $F$-dependence of
the average rate of replication, which depends on the rates of both polymerase
and exonuclease activities of the DNAP, is in good qualitative agreement with
the corresponding experimental results. We introduce 9 novel distinct {\it
conditional dwell times} of a DNAP. Using the methods of first-passage times,
we also derive the exact analytical expressions for the probability
distributions of these conditional dwell times. The predicted $F$-dependence of
these distributions are, in principle, accessible to single-molecule
experiments.


  Starting from fitness correlation functions, we calculate exact expressions
for the amplitude spectra of fitness landscapes as defined by P.F. Stadler [J.
Math. Chem. 20, 1 (1996)] for common landscape models, including Kauffman's
NK-model, rough Mount Fuji landscapes and general linear superpositions of such
landscapes. We further show that correlations decaying exponentially with the
Hamming distance yield exponentially decaying spectra similar to those reported
recently for a model of molecular signal transduction. Finally, we compare our
results for the model systems to the spectra of various experimentally measured
fitness landscapes. We claim that our analytical results should be helpful when
trying to interpret empirical data and guide the search for improved fitness
landscape models.


  Revealing the clonal composition of a single tumor is essential for
identifying cell subpopulations with metastatic potential in primary tumors or
with resistance to therapies in metastatic tumors. Sequencing technologies
provide an overview of an aggregate of numerous cells, rather than
subclonal-specific quantification of aberrations such as single nucleotide
variants (SNVs). Computational approaches to de-mix a single collective signal
from the mixed cell population of a tumor sample into its individual components
are currently not available. Herein we propose a framework for deconvolving
data from a single genome-wide experiment to infer the composition, abundance
and evolutionary paths of the underlying cell subpopulations of a tumor. The
method is based on the plausible biological assumption that tumor progression
is an evolutionary process where each individual aberration event stems from a
unique subclone and is present in all its descendants subclones. We have
developed an efficient algorithm (TrAp) for solving this mixture problem. In
silico analyses show that TrAp correctly deconvolves mixed subpopulations when
the number of subpopulations and the measurement errors are moderate. We
demonstrate the applicability of the method using tumor karyotypes and somatic
hypermutation datasets. We applied TrAp to SNV frequency profile from Exome-Seq
experiment of a renal cell carcinoma tumor sample and compared the mutational
profile of the inferred subpopulations to the mutational profiles of twenty
single cells of the same tumor. Despite the large experimental noise, specific
co-occurring mutations found in clones inferred by TrAp are also present in
some of these single cells. Finally, we deconvolve Exome-Seq data from three
distinct metastases from different body compartments of one melanoma patient
and exhibit the evolutionary relationships of their subpopulations.


  Interactions among living organisms, from bacteria colonies to human
societies, are inherently more complex than interactions among particles and
nonliving matter. Group interactions are a particularly important and
widespread class, representative of which is the public goods game. In
addition, methods of statistical physics have proven valuable for studying
pattern formation, equilibrium selection, and self-organisation in evolutionary
games. Here we review recent advances in the study of evolutionary dynamics of
group interactions on structured populations, including lattices, complex
networks and coevolutionary models. We also compare these results with those
obtained on well-mixed populations. The review particularly highlights that the
study of the dynamics of group interactions, like several other important
equilibrium and non-equilibrium dynamical processes in biological, economical
and social sciences, benefits from the synergy between statistical physics,
network science and evolutionary game theory.


  Tumor development is characterized by a compromised balance between cell life
and death decision mechanisms, which are tighly regulated in normal cells.
Understanding this process provides insights for developing new treatments for
fighting with cancer. We present a study of a mathematical model describing
cellular choice between survival and two alternative cell death modalities:
apoptosis and necrosis. The model is implemented in discrete modeling formalism
and allows to predict probabilities of having a particular cellular phenotype
in response to engagement of cell death receptors. Using an original parameter
sensitivity analysis developed for discrete dynamic systems, we determine the
critical parameters affecting cellular fate decision variables that appear to
be critical in the cellular fate decision and discuss how they are exploited by
existing cancer therapies.


  The conformational change of biological macromolecule is investigated from
the point of quantum transition. A quantum theory on protein folding is
proposed. Compared with other dynamical variables such as mobile electrons,
chemical bonds and stretching-bending vibrations the molecular torsion has the
lowest energy and can be looked as the slow variable of the system.
Simultaneously, from the multi-minima property of torsion potential the local
conformational states are well defined. Following the idea that the slow
variables slave the fast ones and using the nonadiabaticity operator method we
deduce the Hamiltonian describing conformational change. It is proved that the
influence of fast variables on the macromolecule can fully be taken into
account through a phase transformation of slow variable wave function. Starting
from the conformation- transition Hamiltonian the nonradiative matrix element
is calculated in two important cases: A, only electrons are fast variables and
the electronic state does not change in the transition process; B, fast
variables are not limited to electrons but the perturbation approximation can
be used. Then, the general formulas for protein folding rate are deduced. The
analytical form of the formula is utilized to study the temperature dependence
of protein folding rate and the curious non-Arrhenius temperature relation is
interpreted. The decoherence time of quantum torsion state is estimated and the
quantum coherence degree of torsional angles in the protein folding is studied
by using temperature dependence data. The proposed folding rate formula gives a
unifying approach for the study of a large class problems of biological
conformational change.


  We present a multi-scale model to study the attachment of spherical particles
with a rigid core, coated with binding ligands and in equilibrium with the
surrounding, quiescent fluid medium. This class of fluid-immersed adhesion is
widespread in many natural and engineering settings. Our theory highlights how
the micro-scale binding kinetics of these ligands, as well as the attractive /
repulsive surface potential in an ionic medium effects the eventual macro-scale
size distribution of the particle aggregates (flocs). The results suggest that
the presence of elastic ligands on the particle surface allow large floc
aggregates by inducing efficient inter-floc collisions (i.e., a large, non-zero
collision factor). Strong electrolytic composition of the surrounding fluid
favors large floc formation as well.


  Systems Biology has emerged in the last years as a new holistic approach
based on the global understanding of cells instead of only being focused on
their individual parts (genes or proteins), to better understand the complexity
of human cells. Since the Systems Biology still does not provide the most
accurate answers to our questions due to the complexity of cells and the
limited quality of available information to perform a good gene/protein map
analysis, we have created simpler models to ensure easier analysis of the map
that represents the human cell. Therefore, a virtual organism has been designed
according to the main physiological rules for humans in order to replicate the
human organism and its vital functions. This toy model was constructed by
defining the topology of its genes/proteins and the biological functions
associated to it. There are several examples of these toy models that emulate
natural processes to perform analysis of the virtual life in order to design
the best strategy to understand real life. The strategy applied in this study
combines topological and functional analysis integrating the knowledge about
the relative position of a node among the others in the map with the
conclusions generated by mathematical models that reproduce functional data of
the virtual organism. Our results demonstrate that the combination of both
strategies allows better understanding of our virtual organism even with the
lower input of information needed and therefore it can be a potential tool to
better understand the real life.


  Evolutionary graph theory studies the evolutionary dynamics of populations
structured on graphs. A central problem is determining the probability that a
small number of mutants overtake a population. Currently, Monte Carlo
simulations are used for estimating such fixation probabilities on general
directed graphs, since no good analytical methods exist. In this paper, we
introduce a novel deterministic framework for computing fixation probabilities
for strongly connected, directed, weighted evolutionary graphs under neutral
drift. We show how this framework can also be used to calculate the expected
number of mutants at a given time step (even if we relax the assumption that
the graph is strongly connected), how it can extend to other related models
(e.g. voter model), how our framework can provide non-trivial bounds for
fixation probability in the case of an advantageous mutant, and how it can be
used to find a non-trivial lower bound on the mean time to fixation. We provide
various experimental results determining fixation probabilities and expected
number of mutants on different graphs. Among these, we show that our method
consistently outperforms Monte Carlo simulations in speed by several orders of
magnitude. Finally we show how our approach can provide insight into synaptic
competition in neurology.


  Breakage-Fusion-Bridge (BFB) is a mechanism of genomic instability
characterized by the joining and subsequent tearing apart of sister chromatids.
When this process is repeated during multiple rounds of cell division, it leads
to patterns of copy number increases of chromosomal segments as well as
fold-back inversions where duplicated segments are arranged head-to-head. These
structural variations can then drive tumorigenesis.
  BFB can be observed in progress using cytogenetic techniques, but generally
BFB must be inferred from data like microarrays or sequencing collected after
BFB has ceased. Making correct inferences from this data is not
straightforward, particularly given the complexity of some cancer genomes and
BFB's ability to generate a wide range of rearrangement patterns.
  Here we present algorithms to aid the interpretation of evidence for BFB. We
first pose the BFB count vector problem: given a chromosome segmentation and
segment copy numbers, decide whether BFB can yield a chromosome with the given
segment counts. We present the first linear-time algorithm for the problem,
improving a previous exponential-time algorithm. We then combine this algorithm
with fold-back inversions to develop tests for BFB. We show that, contingent on
assumptions about cancer genome evolution, count vectors and fold-back
inversions are sufficient evidence for detecting BFB. We apply the presented
techniques to paired-end sequencing data from pancreatic tumors and confirm a
previous finding of BFB as well as identify a new chromosomal region likely
rearranged by BFB cycles, demonstrating the practicality of our approach.


  Two blind source separation methods (Independent Component Analysis and
Non-negative Matrix Factorization), developed initially for signal processing
in engineering, found recently a number of applications in analysis of
large-scale data in molecular biology. In this short review, we present the
common idea behind these methods, describe ways of implementing and applying
them and point out to the advantages compared to more traditional statistical
approaches. We focus more specifically on the analysis of gene expression in
cancer. The review is finalized by listing available software implementations
for the methods described.


  We demonstrate simultaneous measurements of DNA translocation into glass
nanopores using ionic current detection and fluorescent imaging. We verify the
correspondence between the passage of a single DNA molecule through the
nanopore and the accompanying characteristic ionic current blockage. By
tracking the motion of individual DNA molecules in the nanocapillary
perpendicular to the optical axis and using a model, we can extract an
effective mobility constant for DNA in our geometry under high electric fields.


  The premise of genetic analysis is that a causal link exists between
phenotypic and allelic variation. Yet it has long been documented that mutant
phenotypes are not a simple result of a single DNA lesion, but rather are due
to interactions of the focal allele with other genes and the environment.
Although an experimentally rigorous approach, focusing on individual mutations
and isogenic control strains, has facilitated amazing progress within genetics
and related fields, a glimpse back suggests that a vast complexity has been
omitted from our current understanding of allelic effects. Armed with
traditional genetic analyses and the foundational knowledge they have provided,
we argue that the time and tools are ripe to return to the under-explored
aspects of gene function and embrace the context-dependent nature of genetic
effects. We assert that a broad understanding of genetic effects and the
evolutionary dynamics of alleles requires identifying how mutational outcomes
depend upon the wild-type genetic background. Furthermore, we discuss how best
to exploit genetic background effects to broaden genetic research programs.


  The evolutionary spread of cheater strategies can destabilize populations
engaging in cooperative behaviors, thus demonstrating that evolutionary changes
can have profound implications for populations dynamics. At the same time, the
relative fitness of cooperative traits often depends upon population density,
thus leading to the potential for bidirectional coupling between population
density and the evolution of a cooperative trait. Despite the potential
importance of these eco evolutionary feedback loops, they have not yet been
demonstrated experimentally in social species and their ecological implications
are poorly understood. Here, we demonstrate the presence of a strong feedback
loop between population dynamics and the evolutionary dynamics of a social
microbial gene, SUC2, in laboratory yeast populations whose cooperative growth
is mediated by the SUC2 gene. We directly visualize eco evolutionary
trajectories of hundreds of populations over fifty generations, allowing us to
characterize the phase space describing the interplay of evolution and ecology
in this system. Small populations collapse despite continual evolution towards
increased cooperative allele frequencies; large populations with a sufficient
number of cooperators spiral to a stable state of coexistence between
cooperator and cheater strategies. The presence of cheaters does not
significantly affect the equilibrium population density, but it does reduce the
resilience of the population as well as its ability to adapt to a rapidly
deteriorating environment. Our results demonstrate the potential ecological
importance of coupling between evolutionary dynamics and the population
dynamics of cooperatively growing organisms, particularly in microbes.


  This article shows the asymptotics of distribution and moments of the size
$X_n$ of the minimal clade of a randomly chosen individual in a
Bolthausen-Sznitman $n$-coalescent for $n\to\infty$. The Bolthausen-Sznitman
$n$-coalescent is a Markov process taking states in the set of partitions of
$\left\{1,\ldots,n\right\}$, where $1,\ldots,n$ are referred to as individuals.
The minimal clade of an individual is the equivalence class the individual is
in at the time of the first coalescence event this individual participates
in.\\ The main tool used is the connection of the Bolthausen-Sznitman
$n$-coalescent with random recursive trees introduced by Goldschmidt and Martin
(see \cite{goldschmidtmartin}). This connection shows that $X_n-1$ is
distributed as the number $M_n$ of all individuals not in the equivalence class
of individual 1 shortly before the time of the last coalescence event. Both
functionals are distributed like the size $RT_{n-1}$ of an uniformly chosen
table in a standard Chinese restaurant process with $n-1$ customers.We give
exact formulae for these distributions.\\ Using the asymptotics of $M_n$ shown
by Goldschmidt and Martin in \cite{goldschmidtmartin}, we see $(\log
n)^{-1}\log X_n$ converges in distribution to the uniform distribution on [0,1]
for $n\to\infty$.\\ We provide the complimentary information that $\frac{\log
n}{n^k}E(X_n^k)\to \frac{1}{k}$ for $n\to\infty$, which is also true for $M_n$
and $RT_n$.


  The annual occurrence of many infectious diseases remains a constant burden
to public health systems. The seasonal patterns in respiratory disease
incidence observed in temperate regions have been attributed to the impact of
environmental conditions on pathogen survival. A model describing the
transmission of an infectious disease by means of a pathogenic state capable of
surviving in an environmental reservoir outside of its host organism is
presented in this paper. The ratio of pathogen lifespan to the duration of the
infectious disease state is found to be a critical parameter in determining
disease dynamics. The introduction of a seasonally forced pathogen inactivation
rate identifies a time delay between peak pathogen survival and peak disease
incidence. The delay is dependent on specific disease parameters and, for
influenza, decreases with increasing reproduction number. The observed seasonal
oscillations are found to have a period identical to that of the seasonally
forced inactivation rate and which is independent of the duration of infection
acquired immunity


  The evolutionary pressures that underlie the large-scale functional
organization of the genome are not well understood in eukaryotes. Recent
evidence suggests that functionally similar genes may colocalize (cluster) in
the eukaryotic genome, suggesting the role of chromatin-level gene regulation
in shaping the physical distribution of coordinated genes. However, few of the
bioinformatic tools currently available allow for a systematic study of gene
colocalization across several, evolutionarily distant species. Kerfuffle is a
web tool designed to help discover, visualize, and quantify the physical
organization of genomes by identifying significant gene colocalization and
conservation across the assembled genomes of available species (currently up to
47, from humans to worms). Kerfuffle only requires the user to specify a list
of human genes and the names of other species of interest. Without further
input from the user, the software queries the e!Ensembl BioMart server to
obtain positional information and discovers homology relations in all genes and
species specified. Using this information, Kerfuffle performs a multi-species
clustering analysis, presents downloadable lists of clustered genes, performs
Monte Carlo statistical significance calculations, estimates how conserved gene
clusters are across species, plots histograms and interactive graphs, allows
users to save their queries, and generates a downloadable visualization of the
clusters using the Circos software. These analyses may be used to further
explore the functional roles of gene clusters by interrogating the enriched
molecular pathways associated with each cluster.


  It is established that prion protein is the sole causative agent in a number
of diseases in humans and animals. However, the nature of conformational
changes that the normal cellular form PrPC undergoes in the conversion process
to a self-replicating state is still not fully understood. The ordered
C-terminus of PrPC proteins has three helices (H1, H2, and H3). Here, we use
the Statistical Coupling Analysis (SCA) to infer co-variations at various
locations using a family of evolutionarily related sequences, and the response
of mouse and human PrPCs to mechanical force to decipher the initiation sites
for transition from PrPC to an aggregation prone PrP* state. The sequence-based
SCA predicts that the clustered residues in non-mammals are localized in the
stable core (near H1) of PrPC whereas in mammalian PrPC they are localized in
the frustrated helices H2 and H3 where most of the pathogenic mutations are
found. Force-extension curves and free energy profiles as a function of
extension of mouse and human PrPC in the absence of disulfide (SS) bond between
residues Cys179 and Cys214, generated by applying mechanical force to the ends
of the molecule, show a sequence of unfolding events starting first with
rupture of H2 and H3. This is followed by disruption of structure in two
strands. Helix H1, stabilized by three salt-bridges, resists substantial force
before unfolding. Force extension profiles and the dynamics of rupture of
tertiary contacts also show that even in the presence of SS bond the
instabilities in most of H3 and parts of H2 still determine the propensity to
form the PrP* state. In mouse PrPC with SS bond there are about ten residues
that retain their order even at high forces.


  To understand sensory processing in neuronal populations, it is necessary to
deliver stimuli to the sensory organs of animals and record evoked population
activity downstream. However, the pathways from sensory input to synaptic
currents in cells that are several synapses removed from sensory organs are
complex. Intrinsic noise and uncontrolled modulatory input from other brain
regions can interfere with the delivery of well-defined stimuli. Here we
investigate the ability of channelrhodopsins to deliver precise time-varying
currents to neurons at any point along the sensory-motor pathway. To do this,
we first deduce the amplitude response function of channelrhodopsin-2 (ChR2)
using a three state Markov model of channel kinetics. With biophysically
realistic parameters, this function supports a relatively broad signal passband
and contains a resonance. We confirm the validity of our predicted amplitude
response function using time-varying optical stimulation of excitatory neurons
that express either wild type ChR2 or the ChR2(H134R) mutant. Together, our
results indicate that ChR2-derived optogenetic tools are useful for delivering
repeatable, time-varying currents to genetically-specified populations over a
physiologically-relevant frequency band.


  SLiM is an efficient forward population genetic simulation designed for
studying the effects of linkage and selection on a chromosome-wide scale. The
program can incorporate complex scenarios of demography and population
substructure, various models for selection and dominance of new mutations,
arbitrary gene and chromosomal structure, and user-defined recombination maps.


  This is part of a series of weekly influenza forecasts made during the
2012-2013 influenza season. Here we present results of forecasts initiated
following assimilation of observations for Week 1 (i.e. the forecast begins
January 6, 2013) for municipalities in the United States. These forecasts were
performed on January 11, 2013. Results from forecasts initiated the six
previous weeks (Weeks 47-52) are also presented. The accuracy of these
predictions will not be known for certain until the conclusion of the current
influenza season; however, at the moment a number of the forecasted peaks
appear to be inaccurate. This inaccuracy may be due to the virulence of
influenza this season, which appears to be sending more influenza-infected
persons to seek medical attention and inflates ILI levels (and possibly the
proportion testing influenza positive) relative to years with milder flu
strains. New forecasts that adjust, or scale, for this difference and match the
two focus cities that appear to have already peaked are identified. These new
forecasts will be used, in addition to the previously scaled forms, to make
influenza predictions for the remainder of the season.


  The Rock-Paper-Scissors (RPS) game is a widely used model system in game
theory. Evolutionary game theory predicts the existence of persistent cycles in
the evolutionary trajectories of the RPS game, but experimental evidence has
remained to be rather weak. In this work we performed laboratory experiments on
the RPS game and analyzed the social-state evolutionary trajectories of twelve
populations of N=6 players. We found strong evidence supporting the existence
of persistent cycles. The mean cycling frequency was measured to be $0.029 \pm
0.009$ period per experimental round. Our experimental observations can be
quantitatively explained by a simple non-equilibrium model, namely the
discrete-time logit dynamical process with a noise parameter. Our work
therefore favors the evolutionary game theory over the classical game theory
for describing the dynamical behavior of the RPS game.


  The TNF initiated processes are simulated using our integrative computational
model with new quantitative proteomics data obtained on HeLa cell line. In
spite of fact, that the model development is based on limited experimental
information, with missing data estimated by simulations and parametrization
using indirect experiments on different cell lines, the actual model with HeLa
experimental data provides pleasant agreement between the computationally
predicted behavior and experiments. Since no significant departures from
experimentally determined functionality of HeLa cells has been observed, this
can be taken as encouraging sign, that substantial features of the cell line
behavior are reproduced correctly and indicate new possibilities for cell-line
specific computational experiments.


  Synonymous sites are generally assumed to be subject to weak selective
constraint. For this reason, they are often neglected as a possible source of
important functional variation. We use site frequency spectra from deep
population sequencing data to show that, contrary to this expectation, 22% of
four-fold synonymous (4D) sites in D. melanogaster evolve under very strong
selective constraint while few, if any, appear to be under weak constraint.
Linking polymorphism with divergence data, we further find that the fraction of
synonymous sites exposed to strong purifying selection is higher for those
positions that show slower evolution on the Drosophila phylogeny. The function
underlying the inferred strong constraint appears to be separate from splicing
enhancers, nucleosome positioning, and the translational optimization
generating canonical codon bias. The fraction of synonymous sites under strong
constraint within a gene correlates well with gene expression, particularly in
the mid-late embryo, pupae, and adult developmental stages. Genes enriched in
strongly constrained synonymous sites tend to be particularly functionally
important and are often involved in key developmental pathways. Given that the
observed widespread constraint acting on synonymous sites is likely not limited
to Drosophila, the role of synonymous sites in genetic disease and adaptation
should be reevaluated.


  The evolution of animals involved acquisition of an emergent gene repertoire
for gastrulation. Whether loss of genes also co-evolved with this developmental
reprogramming has not yet been addressed. Here, we identify twenty-four genetic
functions that are retained in fungi and choanoflagellates but undetectable in
animals. These lost genes encode: (i) sixteen distinct biosynthetic functions;
(ii) the two ancestral eukaryotic ClpB disaggregases, Hsp78 and Hsp104, which
function in the mitochondria and cytosol, respectively; and (iii) six other
assorted functions. We present computational and experimental data that are
consistent with a joint function for the differentially localized ClpB
disaggregases, and with the possibility of a shared client/chaperone
relationship between the mitochondrial Fe/S homoaconitase encoded by the lost
LYS4 gene and the two ClpBs. Our analyses lead to the hypothesis that the
evolution of gastrulation-based multicellularity in animals led to efficient
extraction of nutrients from dietary sources, loss of natural selection for
maintenance of energetically expensive biosynthetic pathways, and subsequent
loss of their attendant ClpB chaperones.


  Wang-Landau sampling has been applied to investigate the thermodynamics and
structural properties of a lattice hydrophobic-polar heteropolymer (the HP
protein model) interacting with an attractive substrate. For simplicity, we
consider a short HP sequence consisting of only 36 monomers interacting with a
substrate which attracts all monomers in the sequence. The conformational
"phase transitions" have been identified by a canonical analysis of the
specific heat and suitable structural observables. Three major "transitions",
namely, adsorption, hydrophobic core formation and "flattening" of adsorbed
structures, are observed. Depending on the surface attractive strength relative
to the intra-protein attraction among the H monomers, these processes take
place in different sequences upon cooling.


  The thermodynamic behavior and structural properties of hydrophobic-polar
(HP) lattice proteins interacting with attractive surfaces are studied by means
of Wang-Landau sampling. Three benchmark HP sequences (48mer, 67mer, and
103mer) are considered with different types of surfaces, each of which attract
either all monomers, only hydrophobic (H) monomers, or only polar (P) monomers,
respectively. The diversity of folding behavior in dependence of surface
strength is discussed. Analyzing the combined patterns of various structural
observables, such as, e.g., the derivatives of the numbers of surface contacts,
together with the specific heat, we are able to identify generic categories of
folding and transition hierarchies. We also infer a connection between these
transition categories and the relative surface strengths, i.e., the ratio of
the surface attractive strength to the interchain attraction among H monomers.
The validity of our proposed classification scheme is reinforced by the
analysis of additional benchmark sequences. We thus believe that the folding
hierarchies and identification scheme are generic for HP proteins interacting
with attractive surfaces, regardless of chain length, sequence, or surface
attraction.


  We describe a class of "bare bones" models of homopolymers which undergo
coil-globule collapse and proteins which fold into their native states in free
space or into denatured states when captured by an attractive substrate as the
temperature is lowered. We then show how, with the use of a properly chosen
trial move set, Wang-Landau Monte Carlo sampling can be used to study the rough
free energy landscape and ground (native) states of these intriguingly simple
systems and thus elucidate their thermodynamic complexity.


  We introduce the anti-profile Support Vector Machine (apSVM) as a novel
algorithm to address the anomaly classification problem, an extension of
anomaly detection where the goal is to distinguish data samples from a number
of anomalous and heterogeneous classes based on their pattern of deviation from
a normal stable class. We show that under heterogeneity assumptions defined
here that the apSVM can be solved as the dual of a standard SVM with an
indirect kernel that measures similarity of anomalous samples through
similarity to the stable normal class. We characterize this indirect kernel as
the inner product in a Reproducing Kernel Hilbert Space between representers
that are projected to the subspace spanned by the representers of the normal
samples. We show by simulation and application to cancer genomics datasets that
the anti-profile SVM produces classifiers that are more accurate and stable
than the standard SVM in the anomaly classification setting.


  Fast and cheaper next generation sequencing technologies will generate
unprecedentedly massive and highly-dimensional genomic and epigenomic variation
data. In the near future, a routine part of medical record will include the
sequenced genomes. A fundamental question is how to efficiently extract genomic
and epigenomic variants of clinical utility which will provide information for
optimal wellness and interference strategies. Traditional paradigm for
identifying variants of clinical validity is to test association of the
variants. However, significantly associated genetic variants may or may not be
usefulness for diagnosis and prognosis of diseases. Alternative to association
studies for finding genetic variants of predictive utility is to systematically
search variants that contain sufficient information for phenotype prediction.
To achieve this, we introduce concepts of sufficient dimension reduction and
coordinate hypothesis which project the original high dimensional data to very
low dimensional space while preserving all information on response phenotypes.
We then formulate clinically significant genetic variant discovery problem into
sparse SDR problem and develop algorithms that can select significant genetic
variants from up to or even ten millions of predictors with the aid of dividing
SDR for whole genome into a number of subSDR problems defined for genomic
regions. The sparse SDR is in turn formulated as sparse optimal scoring
problem, but with penalty which can remove row vectors from the basis matrix.
To speed up computation, we develop the modified alternating direction method
for multipliers to solve the sparse optimal scoring problem which can easily be
implemented in parallel. To illustrate its application, the proposed method is
applied to simulation data and the NHLBI's Exome Sequencing Project dataset


  James et al. (2012) presented simulations that apparently falsify the
analytical result by Bastolla et al. (2009), who showed that nested mutualistic
interactions decrease interspecific competition and increase biodiversity in
model ecosystems. This contradiction, however, mainly stems from the incorrect
application of formulas derived for fully connected networks to empirical,
sparse networks.


  Antimicrobial peptides are a class of small, usually positively charged
amphiphilic peptides that are used by the innate immune system to combat
bacterial infection in multicellular eukaryotes. Antimicrobial peptides are
known for their broad-spectrum antimicrobial activity and thus can be used as a
basis for a development of new antibiotics against multidrug-resistant
bacteria. The most challengeous task on the way to a therapeutic use of
antimicrobial peptides is a rational design of new peptides with enhanced
activity and reduced toxicity. Here we report a molecular dynamics and circular
dichroism study of a novel synthetic antimicrobial peptide D51. This peptide
was earlier designed by Loose et al. using a linguistic model of natural
antimicrobial peptides. Molecular dynamics simulation of the peptide folding in
explicit solvent shows fast formation of two antiparallel beta strands
connected by a beta-turn that is confirmed by circular dichroism measurements.
Obtained from simulation amphipatic conformation of the peptide is analysed and
possible mechanism of its interaction with bacterial membranes together with
ways to enhance its antibacterial activity are suggested.


  Although RNAs play many cellular functions little is known about the dynamics
and thermodynamics of these molecules. In principle, all-atom molecular
dynamics simulations can investigate these issues, but with current computer
facilities, these simulations have been limited to small RNAs and to short
times.
  HiRe-RNA, a recently proposed high-resolution coarse-grained for RNA that
captures many geometric details such as base pairing and stacking, is able to
fold RNA molecules to near-native structures in a short computational time. So
far it had been applied to simple hairpins, and here we present its application
to duplexes of a couple dozen nucleotides and show how with our model and with
Replica Exchange Molecular Dynamics (REMD) we can easily predict the correct
double helix from a completely random configuration and study the dissociation
curve. To show the versatility of our model, we present an application to a
double stranded DNA molecule as well.
  A reconstruction algorithm allows us to obtain full atom structures from the
coarse-grained model. Through atomistic Molecular Dynamics (MD) we can compare
the dynamics starting from a representative structure of a low temperature
replica or from the experimental structure, and show how the two are
statistically identical, highlighting the validity of a coarse-grained approach
for structured RNAs and DNAs.


  Background: Significance analysis plays a major role in identifying and
ranking genes, transcription factor binding sites, DNA methylation regions, and
other high-throughput features for association with disease. We propose a new
approach, called gene set bagging, for measuring the stability of ranking
procedures using predefined gene sets. Gene set bagging involves resampling the
original high-throughput data, performing gene-set analysis on the resampled
data, and confirming that biological categories replicate. This procedure can
be thought of as bootstrapping gene-set analysis and can be used to determine
which are the most reproducible gene sets. Results: Here we apply this approach
to two common genomics applications: gene expression and DNA methylation. Even
with state-of-the-art statistical ranking procedures, significant categories in
a gene set enrichment analysis may be unstable when subjected to resampling.
Conclusions: We demonstrate that gene lists are not necessarily stable, and
therefore additional steps like gene set bagging can improve biological
inference of gene set analysis.


  Molecular phenotypes are important links between genomic information and
organismic functions, fitness, and evolution. Complex phenotypes, which are
also called quantitative traits, often depend on multiple genomic loci. Their
evolution builds on genome evolution in a complicated way, which involves
selection, genetic drift, mutations and recombination. Here we develop a
coarse-grained evolutionary statistics for phenotypes, which decouples from
details of the underlying genotypes. We derive approximate evolution equations
for the distribution of phenotype values within and across populations. This
dynamics covers evolutionary processes at high and low recombination rates,
that is, it applies to sexual and asexual populations. In a fitness landscape
with a single optimal phenotype value, the phenotypic diversity within
populations and the divergence between populations reach evolutionary
equilibria, which describe stabilizing selection. We compute the equilibrium
distributions of both quantities analytically and we show that the ratio of
mean divergence and diversity depends on the strength of selection in a
universal way: it is largely independent of the phenotype's genomic encoding
and of the recombination rate. This establishes a new method for the inference
of selection on molecular phenotypes beyond the genome level. We discuss the
implications of our findings for the predictability of evolutionary processes.


  We study the first passage time for a polymer, that we call the narrow
encounter time (NETP), to reach a small target located on the surface of a
microdomain. The polymer is modeled as a Freely Joint Chain (beads connected by
springs with a resting non zero length) and we use Brownian simulations to
study two cases: when (i) any of the monomer or (ii) only one can be absorbed
at the target window. Interestingly, we find that {in the first case} the NETP
is an increasing function of the polymer length until a critical length, after
which it decreases. Moreover, in the long polymer regime, we identified an
exponential scaling law for the NETP as a function of the polymer length. {In
the second case, the position of the absorbed monomer along the polymer chain
strongly influences the NETP}. Our analysis can be applied to estimate the mean
first time of a DNA fragment to a small target in the chromatin structure or
for mRNA to find a small target.


  Next-generation sequencing technologies provide a revolutionary tool for
generating gene expression data. Starting with a fixed RNA sample, they
construct a library of millions of differentially abundant short sequence tags
or "reads", which constitute a fundamentally discrete measure of the level of
gene expression. A common limitation in experiments using these technologies is
the low number or even absence of biological replicates, which complicates the
statistical analysis of digital gene expression data. Analysis of this type of
data has often been based on modified tests originally devised for analysing
microarrays; both these and even de novo methods for the analysis of RNA-seq
data are plagued by the common problem of low replication. We propose a novel,
non-parametric Bayesian approach for the analysis of digital gene expression
data. We begin with a hierarchical model for modelling over-dispersed count
data and a blocked Gibbs sampling algorithm for inferring the posterior
distribution of model parameters conditional on these counts. The algorithm
compensates for the problem of low numbers of biological replicates by
clustering together genes with tag counts that are likely sampled from a common
distribution and using this augmented sample for estimating the parameters of
this distribution. The number of clusters is not decided a priori, but it is
inferred along with the remaining model parameters. We demonstrate the ability
of this approach to model biological data with high fidelity by applying the
algorithm on a public dataset obtained from cancerous and non-cancerous neural
tissues.


  The evolutionary process has been modelled in many ways using both stochastic
and deterministic models. We develop an algebraic model of evolution in a
population of asexually reproducing organisms in which we represent a
stochastic walk in phenotype space, constrained to the edges of an underlying
graph representing the genotype, with a time-homogeneous Markov Chain. We show
its equivalence to a more standard, explicit stochastic model and show the
algebraic model's superiority in computational efficiency. Because of this
increase in efficiency, we offer the ability to simulate the evolution of much
larger populations in more realistic genotype spaces. Further, we show how the
algebraic properties of the Markov Chain model can give insight into the
evolutionary process and allow for analysis using familiar linear algebraic
methods.


  Phylogeographic studies frequently reveal multiple morphologically-cryptic
lineages within species. What is yet unclear is whether such lineages represent
nascent species or evolutionary ephemera. To address this question, we compare
five contact zones, each of which occurs between eco-morphologically cryptic
lineages of rainforest skinks from the rainforests of the Australian Wet
Tropics. Although the contacts likely formed concurrently in response to
Holocene expansion from glacial refugia, we estimate that the divergence times
(t) of the lineage-pairs range from 3.1 to 11.5 Myr. Multilocus analyses of the
contact zones yielded estimates of reproductive isolation that are tightly
correlated with divergence time and, for longer-diverged lineages (t > 5 Myr),
substantial. These results show that phylogeographic splits of increasing depth
can represent stages along the speciation continuum, even in the absence of
overt change in ecologically relevant morphology.


  Biomedical advances have led to a relaxation of natural selection in the
human population in developed countries. In the absence of strong purifying
selection, spontaneous and frequently deleterious mutations tend to accumulate
in the human genome and gradually increase the genetic load; that is, the
frequency of potentially lethal genes in the gene pool. Gradual increase in
incidence of many hereditary disorders suggests deleterious impact of the
genetic load on human well being. Many complex diseases such as type 1 and 2
diabetes, autism, bipolar disorder, allergies, Alzheimer disease, and some
cancers show significantly higher concordance in monozygotic (MZ) twins than in
dizygotic twins (DZ) or parent-child pairs, suggesting their etiology is
strongly influenced by genetics. Preventing these diseases based on genetic
data alone is frequently impossible. We hypothesize that the incidence of
complex diseases could be significantly reduced in the future through a
strategy based on time-separated twinning. This strategy involves the
collection and fertilization of human oocytes followed by several rounds of
artificial twinning. If preimplantation genetic screening (PGS) reports no
aneuploidy or known Mendelian disorders, one of the MZ siblings would be
implanted and the remaining embryos cryoconserved. Once the health of the adult
MZ sibling(s) is established, subsequent parenthood with the cryoconserved
twins could substantially lower the incidence of hereditary disorders with
Mendelian or complex etiology. The proposed method of artificial twinning has
the potential to alleviate suffering and reduce the negative social impact
induced by dysgenic effects associated with known and unknown genetic factors.
Time-separated twinning has the capacity to prevent further accumulation of the
genetic load and to provide source of isogenic embryonic stem cells for future
regenerative therapies.


  We present a mathematical model, and the corresponding mathematical analysis,
that justifies and quantifies the use of principal component analysis of
biallelic genetic marker data for a set of individuals to detect the number of
subpopulations represented in the data. We indicate that the power of the
technique relies more on the number of individuals genotyped than on the number
of markers.


  The time required to transcribe genes with long primary transcripts may limit
their ability to be expressed in cells with short mitotic cycles, a phenomenon
termed intron delay. As such short cycles are a hallmark of the earliest stages
of insect development, we used Drosophila developmental timecourse expression
data to test whether intron delay affects gene expression genome-wide, and to
determine its consequences for the evolution of gene structure. We find that
long zygotically expressed, but not maternally deposited, genes show
substantial delay in expression relative to their shorter counterparts and that
this delay persists over a substantial portion of the ~24 hours of
embryogenesis. Patterns of RNA-seq coverage from the 5' and 3' ends of
transcripts show that this delay is consistent with their inability to
terminate transcription, but not with transcriptional initiation-based
regulatory control. Highly expressed zygotic genes are subject to purifying
selection to maintain compact transcribed regions, allowing conservation of
embryonic expression patterns across the Drosophila phylogeny. We propose that
intron delay is an underappreciated physical mechanism affecting both patterns
of expression as well as gene structure of many genes across Drosophila.


  Side chain flexibility is an important factor in ligand binding. In order to
determine the extent to which side chain flexibility is involved in ligand
binding, a knowledge-based approach was taken. A database composed of examples
of protein structures in the presence or absence of a given ligand is used to
analyze which side chains undergo side chain conformational changes. Such an
analysis has determined that up to 40% of binding site do not present side
chain conformational changes. A total of three residues undergoing side chain
conformational changes encompass approximately 85% of the binding sites
studied. When analyzing the propensities of different amino acids to undergo
side chain conformational changes we find that there are considerable
differences between different amino acids. A support vector machine learning
approach was used to create a classifier system utilizing information about the
solvent accessible area as well as flexibility scale value of each specific
side chain to be predicted together with its neighboring side chains. An
accuracy level of 70% is reached using this approach. The fact that a small
number of residues undergo side chain conformational changes in the majority of
binding sites makes it feasible to introduce side chain flexibility in docking
simulations. An algorithm has been developed for introducing side chain
flexibility utilizing a hybrid genetic-algorithm/exhaustive-search procedure
and a surface complementarity based scoring function. This approach is
implemented in the software tool FlexAID. FlexAID utilizes a rotamer library to
create alternative conformations for a list of residues that are exhaustively
searched during the docking simulation. The performance of FlexAID in rigid
local as well as global simulations falls in the 70-80% range for both local
and global simulations.


  The cell cycle is a tightly controlled process, yet its underlying genetic
network shows marked differences across species. Which of the associated
structural features follow solely from the ability to impose the appropriate
gene expression patterns? We tackle this question in silico by examining the
ensemble of all regulatory networks which satisfy the constraint of producing a
given sequence of gene expressions. We focus on three cell cycle profiles
coming from baker's yeast, fission yeast and mammals. First, we show that the
networks in each of the ensembles use just a few interactions that are
repeatedly reused as building blocks. Second, we find an enrichment in network
motifs that is similar in the two yeast cell cycle systems investigated. These
motifs do not have autonomous functions, but nevertheless they reveal a
regulatory logic for cell cycling based on a feed-forward cascade of activating
interactions.


  Phylogenetic tree comparison metrics are an important tool in the study of
evolution, and hence the definition of such metrics is an interesting problem
in phylogenetics. In a paper in Taxon fifty years ago, Sokal and Rohlf proposed
to measure quantitatively the difference between a pair of phylogenetic trees
by first encoding them by means of their half-matrices of cophenetic values,
and then comparing these matrices. This idea has been used several times since
then to define dissimilarity measures between phylogenetic trees but, to our
knowledge, no proper metric on weighted phylogenetic trees with nested taxa
based on this idea has been formally defined and studied yet. Actually, the
cophenetic values of pairs of different taxa alone are not enough to single out
phylogenetic trees with weighted arcs or nested taxa. In this paper we define a
family of cophenetic metrics that compare phylogenetic trees on a same set of
taxa by encoding them by means of their vectors of cophenetic values of pairs
of taxa and depths of single taxa, and then computing the $L^p$ norm of the
difference of the corresponding vectors. Then, we study, either analytically or
numerically, some of their basic properties: neighbors, diameter, distribution,
and their rank correlation with each other and with other metrics.


  For a nonnegative matrix A and real diagonal matrix D, two known inequalities
on the spectral radius, r(A^2 D^2) >= r(AD)^2 and r(A) r(A D^2) >= r(AD)^2,
leave open the question of what determines the order of r(A^2 D^2) with respect
to r(A) r(A D^2). This is a special case of a broad class of problems that
arise repeatedly in ecological and evolutionary dynamics. Here, sufficient
conditions are found on A that determine orders in either direction. For a
diagonally symmetrizable nonnegative matrix A with all positive eigenvalues and
nonnegative D, r(A^2 D) <= r(A) r(AD). The reverse holds if all of the
eigenvalues of A are negative besides the Perron root. This is a particular
case of the more general result that r(A [(1-m) r(B) I + m B] D) is monotonic
in m when all non-Perron eigenvalues of A have the same sign --- decreasing for
positive signs and increasing for negative signs, for symmetrizable nonnegative
A and B that commute. Commuting matrices include the Kronecker products A, B in
{\otimes_{i=1}^L M_i^{t_i}}, t_i in {0, 1, 2, ...}, which comprise a class of
application for these results. This machinery allows analysis of the sign of
d/d m_j r({\otimes_{i=1}^L[(1-m_i) r(A_i) I_i + m_i A_i]}D). The eigenvalue
sign conditions also provide lower or upper bounds to the harmonic mean of the
expected sojourn times of Markov chains. These inequalities appear in the
asymptotic growth rates of viral quasispecies, models for the evolution of
dispersal in random environments, and the evolution of site-specific mutation
rates over the entire genome.


  Probabilistic models over strings have played a key role in developing
methods allowing indels to be treated as phylogenetically informative events.
There is an extensive literature on using automata and transducers on
phylogenies to do inference on these probabilistic models, in which an
important theoretical question in the field is the complexity of computing the
normalization of a class of string-valued graphical models. This question has
been investigated using tools from combinatorics, dynamic programming, and
graph theory, and has practical applications in Bayesian phylogenetics. In this
work, we revisit this theoretical question from a different point of view,
based on linear algebra. The main contribution is a new proof of a known result
on the complexity of inference on TKF91, a well-known probabilistic model over
strings. Our proof uses a different approach based on classical linear algebra
results, and is in some cases easier to extend to other models. The proving
method also has consequences on the implementation and complexity of inference
algorithms.


  Regulatory relationships of 686 intronic miRNA and 784 intergenic miRNAs with
mRNAs of 51 intronic miRNA coding genes were established. Interaction features
of studied miRNAs with 5'UTR, CDS and 3'UTR of mRNA of each gene were revealed.
Functional regions of mRNA were shown to be significantly heterogenous
according to the number of binding sites of miRNA and to the location density
of these sites.


  The cophenetic metrics $d_{\varphi,p}$, for $p\in {0}\cup[1,\infty[$, are a
recent addition to the kit of available distances for the comparison of
phylogenetic trees. Based on a fifty years old idea of Sokal and Rohlf, these
metrics compare phylogenetic trees on a same set of taxa by encoding them by
means of their vectors of cophenetic values of pairs of taxa and depths of
single taxa, and then computing the $L^p$ norm of the difference of the
corresponding vectors. In this paper we compute the expected value of the
square of $d_{\varphi,2}$ on the space of fully resolved rooted phylogenetic
trees with $n$ leaves, under the Yule and the uniform probability
distributions.


  We have recently used genetic programming to automatically generate an
improved version of Langmead's DNA read alignment tool Bowtie2 Sect.5.3
RN/12/09. We find it runs more than four times faster than the Bioinformatics
sequencing tool (BWA) currently used with short next generation paired end DNA
sequences by the Cancer Institute, takes less memory and yet finds similar
matches in the human genome.


  Withdrawn by arXiv administrators due to content entirely plagiarized from
other authors (not in arXiv).


  High-throughput sequencing of RNA transcripts (RNA-seq) has become the method
of choice for detection of differential expression (DE). Concurrent with the
growing popularity of this technology there has been a significant research
effort devoted towards understanding the statistical properties of this data
and the development of analysis methods. We report on a comprehensive
evaluation of the commonly used DE methods using the SEQC benchmark data set.
We evaluate a number of key features including: assessment of normalization,
accuracy of DE detection, modeling of genes expressed in only one condition,
and the impact of sequencing depth and number of replications on identifying DE
genes. We find significant differences among the methods with no single method
consistently outperforming the others. Furthermore, the performance of
array-based approach is comparable to methods customized for RNA-seq data.
Perhaps most importantly, our results demonstrate that increasing the number of
replicate samples provides significantly more detection power than increased
sequencing depth.


  Three steps aid in the analysis of selection. First, describe phenotypes by
their component causes. Components include genes, maternal effects, symbionts,
and any other predictors of phenotype that are of interest. Second, describe
fitness by its component causes, such as an individual's phenotype, its
neighbors' phenotypes, resource availability, and so on. Third, put the
predictors of phenotype and fitness into an exact equation for evolutionary
change, providing a complete expression of selection and other evolutionary
processes. The complete expression separates the distinct causal roles of the
various hypothesized components of phenotypes and fitness. Traditionally, those
components are given by the covariance, variance, and regression terms of
evolutionary models. I show how to interpret those statistical expressions with
respect to information theory. The resulting interpretation allows one to read
the fundamental equations of selection and evolution as sentences that express
how various causes lead to the accumulation of information by selection and the
decay of information by other evolutionary processes. The interpretation in
terms of information leads to a deeper understanding of selection and
heritability, and a clearer sense of how to formulate causal hypotheses about
evolutionary process. Kin selection appears as a particular type of causal
analysis that partitions social effects into meaningful components.


  Background - The process of generating raw genome sequence data continues to
become cheaper, faster, and more accurate. However, assembly of such data into
high-quality, finished genome sequences remains challenging. Many genome
assembly tools are available, but they differ greatly in terms of their
performance (speed, scalability, hardware requirements, acceptance of newer
read technologies) and in their final output (composition of assembled
sequence). More importantly, it remains largely unclear how to best assess the
quality of assembled genome sequences. The Assemblathon competitions are
intended to assess current state-of-the-art methods in genome assembly. Results
- In Assemblathon 2, we provided a variety of sequence data to be assembled for
three vertebrate species (a bird, a fish, and snake). This resulted in a total
of 43 submitted assemblies from 21 participating teams. We evaluated these
assemblies using a combination of optical map data, Fosmid sequences, and
several statistical methods. From over 100 different metrics, we chose ten key
measures by which to assess the overall quality of the assemblies. Conclusions
- Many current genome assemblers produced useful assemblies, containing a
significant representation of their genes, regulatory sequences, and overall
genome structure. However, the high degree of variability between the entries
suggests that there is still much room for improvement in the field of genome
assembly and that approaches which work well in assembling the genome of one
species may not necessarily work well for another.


  Extremely low frequency electromagnetic fields aren't considered as a real
carcinogenic agent despite the fact that some studies have showed impairment of
the DNA integrity in different cells lines. The aim of this study was
evaluation of the late effects of a 100 Hz and 5.6 mT electromagnetic field,
applied continuously or discontinuously, on the DNA integrity of Vero cells
assessed by alkaline Comet assay and by cell cycle analysis. Normal Vero cells
were exposed to extremely low frequency electromagnetic fields (100 Hz, 5.6 mT)
for 45 minutes. The Comet assay and cell cycle analysis were performed 48 hours
after the treatment. Exposed samples presented an increase of the number of
cells with high damaged DNA as compared with non-exposed cells. Quantitative
evaluation of the comet assay showed a significantly ($<$0.001) increase of the
tail lengths, of the quantity of DNA in tail and of Olive tail moments,
respectively. The analysis of the registered comet indices showed that an
extremely low frequency electromagnetic field of 100 Hz and 5.6 mT had a
genotoxic impact on Vero cells. Cell cycle analysis showed an increase of the
frequency of the cells in S phase, proving the occurrence of single strand
breaks. The most probable mechanism of induction of the registered effects is
the production of different types of reactive oxygen species.


  We investigate the relationship between communication and search efficiency
in a biological context by proposing a model of Brownian searchers with
long-range pairwise interaction. After a general study of the properties of the
model, we show an application to the particular case of acoustic communication
among Mongolian gazelle, for which data are available, searching for good
habitat areas. Using Monte Carlo simulations and density equations, our results
point out that the search is optimal (i.e. the mean first hitting time among
searchers is minimum) at intermediate scales of communication, showing that
both an excess and a lack of information may worsen it.


  Under dilute in vitro conditions transcription factors rapidly locate their
target sequence on DNA by using the facilitated diffusion mechanism. However,
whether this strategy of alternating between three-dimensional bulk diffusion
and one-dimensional sliding along the DNA contour is still beneficial in the
crowded interior of cells is highly disputed. Here we use a simple model for
the bacterial genome inside the cell and present a semi-analytical model for
the in vivo target search of transcription factors within the facilitated
diffusion framework. Without having to resort to extensive simulations we
determine the mean search time of a lac repressor in a living E. coli cell by
including parameters deduced from experimental measurements. The results agree
very well with experimental findings, and thus the facilitated diffusion
picture emerges as a quantitative approach to gene regulation in living
bacteria cells. Furthermore we see that the search time is not very sensitive
to the parameters characterizing the DNA configuration and that the cell seems
to operate very close to optimal conditions for target localization. Local
searches as implied by the colocalization mechanism are only found to mildly
accelerate the mean search time within our model.


  Integration of data from genome-wide single nucleotide polymorphism (SNP)
association studies of different traits should allow researchers to disentangle
the genetics of potentially related traits within individually associated
regions. Formal statistical colocalisation testing of individual regions, which
requires selection of a set of SNPs summarizing the association in a region. We
show that the SNP selection method greatly affects type 1 error rates, with
published studies having used methods expected to result in substantially
inflated type 1 error rates. We show that either avoiding variable selection
and instead testing the most informative principal components or integrating
over variable selection using Bayesian model averaging can lead to correct
control of type 1 error rates. Application to data from Graves' disease and
Hashimoto's thyroiditis reveals a common genetic signature across seven regions
shared between the diseases, and indicates that in five of six regions
associated with Graves' disease and not Hashimoto's thyroiditis, this more
likely reflects genuine absence of association with the latter rather than lack
of power. Our examination, by simulation, of the performance of colocalisation
tests and associated software will foster more widespread adoption of formal
colocalisation testing. Given the increasing availability of large expression
and genetic association data sets from disease-relevant tissue and purified
cell populations, coupled with identification of regulatory sequences by
projects such as ENCODE, colocalisation analysis has the potential to reveal
both shared genetic signatures of related traits and causal disease genes and
tissues.


  A popular line of research in evolutionary biology is the use of
time-calibrated phylogenies for the inference of diversification processes.
This requires computing the likelihood of a given ultrametric tree as the
reconstructed tree produced by a given model of diversification. Etienne &
Rosindell (2012) proposed a lineage-based model of diversification, called
protracted speciation, where species remain incipient during a random duration
before turning good species, and showed that this can explain the slowdown in
lineage accumulation observed in real phylogenies. However, they were unable to
provide a general likelihood formula. Here, we present a likelihood formula for
protracted speciation models, where rates at which species turn good or become
extinct can depend both on their age and on time. Our only restrictive
assumption is that speciation rate does not depend on species status.
  Our likelihood formula utilizes a new technique, based on the contour of the
phylogenetic tree and first developed in Lambert (2010). We consider the
reconstructed trees spanned by all extant species, by all good extant species,
or by all representative species, which are either good extant species or
incipient species representative of some good extinct species. Specifically, we
prove that each of these trees is a coalescent point process, that is, a
planar, ultrametric tree where the coalescence times between two consecutive
tips are independent, identically distributed random variables. We characterize
the common distribution of these coalescence times in some, biologically
meaningful, special cases for which the likelihood reduces to an elegant
analytical formula or becomes numerically tractable.


  We examine the two-dimensional extension of the model of Kessler and Sander
of competition between two species identical except for dispersion rates. In
this class of models, the spatial inhomogeneity of reproduction rates gives
rise to an implicit cost of dispersal, due to the tendency to leave favorable
locations. Then, as in the Hamilton-May model with its explicit dispersal cost,
the tradeoff between dispersal case and the beneficial role of dispersal in
limiting fluctuations, leads to an advantage of one dispersal rate over
another, and the eventual extinction of the disadvantaged species. In two
dimensions we find that while the competition leads to the elimination of one
species at high and low population density, at intermediate densities the two
species can coexist essentially indefinitely. This is a new phenomenon not
present in either the one-dimensional form of the Kessler-Sander model nor in
the totally connected Hamilton-May model, and points to the importance of
geometry in the question of dispersal.


  The evolution of antibiotic resistance among bacteria threatens our continued
ability to treat infectious diseases. The need for sustainable strategies to
cure bacterial infections has never been greater. So far, all attempts to
restore susceptibility after resistance has arisen have been unsuccessful,
including restrictions on prescribing [1] and antibiotic cycling [2,3]. Part of
the problem may be that those efforts have implemented different classes of
unrelated antibiotics, and relied on removal of resistance by random loss of
resistance genes from bacterial populations (drift). Here, we show that
alternating structurally similar antibiotics can restore susceptibility to
antibiotics after resistance has evolved. We found that the resistance
phenotypes conferred by variant alleles of the resistance gene encoding the TEM
{\beta}-lactamase (blaTEM) varied greatly among 15 different {\beta}-lactam
antibiotics. We captured those differences by characterizing complete adaptive
landscapes for the resistance alleles blaTEM-50 and blaTEM-85, each of which
differs from its ancestor blaTEM-1 by four mutations. We identified pathways
through those landscapes where selection for increased resistance moved in a
repeating cycle among a limited set of alleles as antibiotics were alternated.
Our results showed that susceptibility to antibiotics can be sustainably
renewed by cycling structurally similar antibiotics. We anticipate that these
results may provide a conceptual framework for managing antibiotic resistance.
This approach may also guide sustainable cycling of the drugs used to treat
malaria and HIV.


  We characterize the class of symmetric two-player games in which tit-for-tat
cannot be beaten even by very sophisticated opponents in a repeated game. It
turns out to be the class of exact potential games. More generally, there is a
class of simple imitation rules that includes tit-for-tat but also
imitate-the-best and imitate-if-better. Every decision rule in this class is
essentially unbeatable in exact potential games. Our results apply to many
interesting games including all symmetric 2x2 games, and standard examples of
Cournot duopoly, price competition, public goods games, common pool resource
games, and minimum effort coordination games.


  Before cell division, two identical copies of chromosomes are pulled apart by
microtubule (MT) filaments that approach the chromosomes from the opposite
poles a mitotic spindle. Connection between the MTs and the chromosomes are
mediated by a molecular complex called kinetochore. An externally applied
tension can lead to detachment of the MTs from the kinetochore; the mean
lifetime of such an attachment is essentially a mean first-passage time. In
their in-vitro pioneering single-kinetochore experiments, Akiyoshi et al.
(Nature 468, 576 (2010)), observed that the mean lifetimes of reconstituted
MT-kinetochore attachments vary non-monotonically with increasing tension. The
counter-intuitive stabilization of the attachments by small load forces was
interpreted in terms of a catch-bond-like mechanism based on a phenomenological
2-state kinetic model. Here we develop the first detailed microscopic model for
studying the dependence of the lifetime of the MT-kinetochore attachment on (a)
the structure, (b) energetics, and (c) kinetics of the coupling. The
catch-bond-like mechanism emerges naturally from this model. Moreover,
in-silico experiments on this model reveal further interesting phenomena,
arising from the subtle effects of competing sub-processes, which are likely to
motivate new experiments in this emerging area of single-particle biophysics.


  We present a maximum entropy framework to separate intrinsic and extrinsic
contributions to noisy gene expression solely from the profile of expression.
We express the experimentally accessible probability distribution of the copy
number of the gene product (mRNA or protein) by accounting for possible
variations in extrinsic factors. The distribution of extrinsic factors is
estimated using the maximum entropy principle. Our results show that extrinsic
factors qualitatively and quantitatively affect the probability distribution of
the gene product. We work out, in detail, the transcription of mRNA from a
constitutively expressed promoter in {\it E. coli}. We suggest that the
variation in extrinsic factors may account for the observed {\it wider than
Poisson} distribution of mRNA copy numbers. We successfully test our framework
on a numerical simulation of a simple gene expression scheme that accounts for
the variation in extrinsic factors. We also make falsifiable predictions, some
of which are tested on previous experiments in {\it E. coli} while others need
verification. Application of the current framework to more complex situations
is also discussed.


  We present schema redescription as a methodology to characterize canalization
in automata networks used to model biochemical regulation and signalling. In
our formulation, canalization becomes synonymous with redundancy present in the
logic of automata. This results in straightforward measures to quantify
canalization in an automaton (micro-level), which is in turn integrated into a
highly scalable framework to characterize the collective dynamics of
large-scale automata networks (macro-level). This way, our approach provides a
method to link micro- to macro-level dynamics -- a crux of complexity. Several
new results ensue from this methodology: uncovering of dynamical modularity
(modules in the dynamics rather than in the structure of networks),
identification of minimal conditions and critical nodes to control the
convergence to attractors, simulation of dynamical behaviour from incomplete
information about initial conditions, and measures of macro-level canalization
and robustness to perturbations. We exemplify our methodology with a well-known
model of the intra- and inter cellular genetic regulation of body segmentation
in Drosophila melanogaster. We use this model to show that our analysis does
not contradict any previous findings. But we also obtain new knowledge about
its behaviour: a better understanding of the size of its wild-type attractor
basin (larger than previously thought), the identification of novel minimal
conditions and critical nodes that control wild-type behaviour, and the
resilience of these to stochastic interventions. Our methodology is applicable
to any complex network that can be modelled using automata, but we focus on
biochemical regulation and signalling, towards a better understanding of the
(decentralized) control that orchestrates cellular activity -- with the
ultimate goal of explaining how do cells and tissues 'compute'.


  We propose a simple model with two infective classes in order to model the
cholera epidemic in Haiti. We include the impact of environmental events
(rainfall, temperature and tidal range) on the epidemic in the Artibonite and
Ouest regions by introducing terms in the transmission rate that vary with
environmental conditions. We fit the model on weekly data from the beginning of
the epidemic until December 2013, including the vaccination programs that were
recently undertaken in the Ouest and Artibonite regions. We then modified these
projections excluding vaccination to assess the programs' effectiveness. Using
real-time daily rainfall, we found lag times between precipitation events and
new cases that range from 3.4 to 8.4 weeks in Artibonite and 5.1 to 7.4 in
Ouest. In addition, it appears that, in the Ouest region, tidal influences play
a significant role in the dynamics of the disease. Intervention efforts of all
types have reduced case numbers in both regions; however, persistent outbreaks
continue. In Ouest, where the population at risk seems particularly besieged
and the overall population is larger, vaccination efforts seem to be taking
hold more slowly than in Artibonite, where a smaller core population was
vaccinated. The models including the vaccination programs predicted that a year
and six months later, the mean number of cases in Artibonite would be reduced
by about two thousand cases, and in Ouest by twenty four hundred cases below
that predicted by the models without vaccination. We also found that
vaccination is best when done in the early spring, and as early as possible in
the epidemic. Comparing vaccination between the first spring and the second,
there is a drop of about 40% in the case reduction due to the vaccine and about
10% per year after that.


  Giant viruses contain large genomes, encode many proteins atypical for
viruses, replicate in large viral factories, and tend to infect protists. The
giant virus replication factories can in turn be infected by so called
virophages, which are smaller viruses that negatively impact giant virus
replication. An example are Mimiviruses that infect the protist Acanthamoeba
and that are themselves infected by the virophage Sputnik. This paper examines
the evolutionary dynamics of this system, using mathematical models. While the
models suggest that the virophage population will evolve to increasing degrees
of giant virus inhibition, it further suggests that this renders the virophage
population prone to extinction due to dynamic instabilities over wide parameter
ranges. Implications and conditions required to avoid extinction are discussed.
Another interesting result is that virophage presence can fundamentally alter
the evolutionary course of the giant virus. While the giant virus is predicted
to evolve towards increasing its basic reproductive ratio in the absence of the
virophage, the opposite is true its presence. Therefore, virophages can not
only benefit the host population directly by inhibiting the giant viruses, but
also indirectly by causing giant viruses to evolve towards weaker phenotypes.
Experimental tests for this model are suggested.


  Practical application of genomic-based risk stratification to clinical
diagnosis is appealing yet performance varies widely depending on the disease
and genomic risk score (GRS) method. Celiac disease (CD), a common
immune-mediated illness, is strongly genetically determined and requires
specific HLA haplotypes. HLA testing can exclude diagnosis but has low
specificity, providing little information suitable for clinical risk
stratification. Using six European CD cohorts, we provide a proof-of-concept
that statistical learning approaches which simultaneously model all SNPs can
generate robust and highly accurate predictive models based on genome-wide SNP
profiles. The high predictive capacity replicated both in cross-validation
within each cohort (AUC of 0.87-0.89) and in independent replication across
cohorts (AUC of 0.86-0.9), despite differences in ethnicity. The models
explained 30-35% of disease variance and up to $\sim43\%$ of heritability. The
GRS's utility was assessed in different clinically relevant settings.
Comparable to HLA typing, the GRS can be used to identify individuals without
CD with $\geq99.6\%$ negative predictive value however, unlike HLA typing,
patients can also be stratified into categories of higher-risk for CD who would
benefit from more invasive and costly definitive testing. The GRS is flexible
and its performance can be adapted to the clinical situation by adjusting the
threshold cut-off. Despite explaining a minority of disease heritability, our
findings indicate a predictive GRS provides clinically relevant information to
improve upon current diagnostic pathways for CD, and support further studies
evaluating the clinical utility of this approach in CD and other complex
diseases.


  Indirect reciprocity is one of the main mechanisms to explain the emergence
and sustainment of altruism in societies. The standard approach to indirect
reciprocity are reputation models. These are games in which players base their
decisions on their opponent's reputation gained in past interactions with other
players (moral assessment). The combination of actions and moral assessment
leads to a large diversity of strategies, thus determining the stability of any
of them against invasions by all the others is a difficult task. We use a
variant of a previously introduced reputation-based model that let us
systematically analyze all these invasions and determine which ones are
successful. Accordingly we are able to identify the third-order strategies
(those which, apart from the action, judge considering both the reputation of
the donor and that of the recipient) that are evolutionarily stable. Our
results reveal that if a strategy resists the invasion of any other one sharing
its same moral assessment, it can resist the invasion of any other strategy.
However, if actions are not always witnessed, cheaters (i.e., individuals with
a probability of defecting regardless of the opponent's reputation) have a
chance to defeat the stable strategies for some choices of the probabilities of
cheating and of being witnessed. Remarkably, by analyzing this issue with
adaptive dynamics we find that whether a honest population resists the invasion
of cheaters is determined by a Hamilton-like rule---with the probability that
the cheat is discovered playing the role of the relatedness parameter.


  Mirtrons are a special type of pre-miRNA which originate from intronic
regions and are spliced directly from the transcript instead of being processed
by Drosha. The splicing mechanism is better understood for the processing of
mRNA for which was established that there is a characteristic CG content around
splice sites. Here we analyse the CG-content ratio of pre-miRNAs and mirtrons
and compare them with their genomic neighbourhood in an attempt to establish
key properties which are easy to evaluate and to understand their biogenesis.
We propose a simple log-ratio of the CG-content comparing the precursor
sequence and is flanking region. We discovered that Caenorhabditis elegans and
Drosophila melanogaster mirtrons, so far without exception, have smaller
CG-content than their genomic neighbourhood. This is markedly different from
usual pre-miRNAs which mostly have larger CG-content when compared to their
genomic neighbourhood. We also analysed some mammalian and primate mirtrons
which, in contrast the invertebrate mirtrons, have higher CG-content ratio.


  Many models have been developed to study the role of branching actin networks
in motility. One important component of those models is the distribution of
filament orientations relative to the cell membrane. Two mean-field models
previously proposed are generalized and analyzed. In particular, we find that
both models uniquely select for a dominant orientation pattern. In the linear
case, the pattern is the eigenfunction associated with the principal
eigenvalue. In the nonlinear case, we show there exists a unique equilibrium
and that the equilibrium is locally stable. Approximate techniques are then
used to provide evidence for global stability.


  In this paper we are investigating the long time behaviour of the solution of
a mutation competition model of Lotka-Volterra's type. Our main motivation
comes from the analysis of the Lotka-Volterra's competition system with
mutation which simulates the demo-genetic dynamics of diverse virus in their
host : $$ \frac{dv_{i}(t)}{dt}=v_i\[r_i-\frac{1}{K}\Psi_i(v)\]+\sum_{j=1}^{N}
\mu_{ij}(v_j-v_i). $$ In a first part we analyse the case where the competition
terms $\Psi_i$ are independent of the virus type $i$. In this situation and
under some rather general assumptions on the functions $\Psi_i$, the
coefficients $r_i$ and the mutation matrix $\mu_{ij}$ we prove the existence of
a unique positive globally stable stationary solution i.e. the solution
attracts the trajectory initiated from any nonnegative initial datum. Moreover
the unique steady state $\bar v$ is strictly positive in the sense that $\bar
v_i>0$ for all $i$. These results are in sharp contrast with the behaviour of
Lotka-Volterra without mutation term where it is known that multiple non
negative stationary solutions exist and an exclusion principle occurs (i.e For
all $i\neq i_0, \bar v_{i}=0$ and $\bar v_{i_0}>0$). Then we explore a typical
example that has been proposed to explain some experimental data. For such
particular models we characterise the speed of convergence to the equilibrium.
In a second part, under some additional assumption, we prove the existence of a
positive steady state for the full system and we analyse the long term
dynamics. The proofs mainly rely on the construction of a relative entropy
which plays the role of a Lyapunov functional.


  Although species longevity is subject to a diverse range of selective forces,
the mortality curves of a wide variety of organisms are rather similar. We
argue that aging and its universal characteristics may have evolved by means of
a gradual increase in the systemic interdependence between a large collection
of biochemical or mechanical components. Modeling the organism as a dependency
network which we create using a constructive evolutionary process, we age it by
allowing nodes to be broken or repaired according to a probabilistic algorithm
that accounts for random failures/repairs and dependencies. Our simulations
show that the network slowly accumulates damage and then catastrophically
collapses. We use our simulations to fit experimental data for the time
dependent mortality rates of a variety of multicellular organisms and even
complex machines such as automobiles. Our study suggests that aging is an
emergent finite-size effect in networks with dynamical dependencies and that
the qualitative and quantitative features of aging are not sensitively
dependent on the details of system structure.


  In experiments, the free energy of transferring the peptide group from water
to an osmolyte solution is obtained using the transfer free energy of (Gly)_n
with the added assumption that a constant incremental change in free energy
with n implies that each additional unit makes an independent contribution to
the free energy. Here we test this assumption and uncover its limitations.
Together with results for cyclic-diglycine, we show that, in principle, it is
not possible to obtain a peptide group transfer free energy that is independent
of the model system. We calculate the hydration free energy of
acetyl-(Gly)_n-methyl amide (n=1..7) peptides modeled in the extended
conformation in water and osmolyte solutions and find that the hydration free
energy is linear in n, suggestive of independent, additive group-contributions.
To probe the observed linearity further, we study the hydration of the solute
bereft of water molecules in the first hydration shell. This conditioned solute
arises naturally in the theoretical formulation and helps us focus on hydration
effects uncluttered by the complexities of short-range solute-water
interactions. We subdivide the conditioned solute into n+1 peptide groups and a
methyl end group. The binding energy of each of these groups with the solvent
is Gaussian distributed, but the near neighbor binding energies are themselves
correlated: the (i,i+1) correlation is the strongest and tends to lower the
free energy over the independent group case. We show that the observed
linearity can be explained by the similarity of near neighbor correlations.
Implications for group additive transfer free energy models are indicated.


  Historical sciences like evolutionary biology reconstruct past events by
using the traces that the past has bequeathed to the present. The Markov Chain
Convergence Theorem and the Data Processing Inequality describe how the mutual
information between present and past is affected by how much time there is in
between. These two results are very general; they pertain to any process, not
just to the biological processes that occur in evolution. To study the
specifically biological question of how the present state of a lineage provides
information about its evolutionary past, we use a Moran process framework and
consider how the kind of evolutionary process (drift, and selection of various
kinds) at work in a lineage affects the epistemological relation of present to
past.


  The genome of bacterial species is much more flexible than that of
eukaryotes. Moreover, the distributed genome hypothesis for bacteria states
that the total number of genes present in a bacterial population is greater
than the genome of every single individual. The pangenome, i.e. the set of all
genes of a bacterial species (or a sample), comprises the core genes which are
present in all living individuals, and accessory genes, which are carried only
by some individuals. In order to use accessory genes for adaptation to
environmental forces, genes can be transferred horizontally between
individuals. Here, we extend the infinitely many genes model from Baumdicker,
Hess and Pfaffelhuber (2010) for horizontal gene transfer. We take a
genealogical view and give a construction -- called the Ancestral Gene Transfer
Graph -- of the joint genealogy of all genes in the pangenome. As application,
we compute moments of several statistics (e.g. the number of differences
between two individuals and the gene frequency spectrum) under the infinitely
many genes model with horizontal gene transfer.


  Fungicide mixtures produced by the agrochemical industry often contain
low-risk fungicides, to which fungal pathogens are fully sensitive, together
with high-risk fungicides known to be prone to fungicide resistance. Can these
mixtures provide adequate disease control while minimizing the risk for the
development of resistance? We present a population dynamics model to address
this question. We found that the fitness cost of resistance is a crucial
parameter to determine the outcome of competition between the sensitive and
resistant pathogen strains and to assess the usefulness of a mixture. If
fitness costs are absent, then the use of the high-risk fungicide in a mixture
selects for resistance and the fungicide eventually becomes nonfunctional. If
there is a cost of resistance, then an optimal ratio of fungicides in the
mixture can be found, at which selection for resistance is expected to vanish
and the level of disease control can be optimized.


  We consider energy relaxation of the long-range spin glass model with sparse
couplings, the so-called dilute Sherrington-Kirkpatrick (SK) model, starting
from a random initial state. We consider the effect that modularity of the
coupling matrix has on this relaxation dynamics. In the absence of finite size
effects, the relaxation dynamics appears independent of modularity. For finite
sizes, a more modular system reaches a less favorable energy at long times. For
small sizes, a more modular system also has a less favorable energy at short
times. For large sizes, modularity appears to lead to slightly more favorable
energies at intermediate times. We discuss these results in the context of
evolutionary theory, where horizontal gene transfer, absent in the Glauber
equilibration dynamics of the SK model studied here, endows modular organisms
with larger response functions at short times.


  We evaluate analytically and numerically the size of the frozen core and
various scaling laws for critical Boolean networks that have a power-law in-
and/or out-degree distribution. To this purpose, we generalize an efficient
method that has previously been used for conventional random Boolean networks
and for networks with power-law in-degree distributions. With this
generalization, we can also deal with power-law out-degree distributions. When
the power-law exponent is between 2 and 3, the second moment of the
distribution diverges with network size, and the scaling exponent of the
nonfrozen nodes depends on the degree distribution exponent. Furthermore, the
exponent depends also on the dependence of the cutoff of the degree
distribution on the system size. Altogether, we obtain an impressive number of
different scaling laws depending on the type of cutoff as well as on the
exponents of the in- and out-degree distributions. We confirm our scaling
arguments and analytical considerations by numerical investigations.


  We investigate numerically and analytically the formation of the frozen core
in critical random Boolean networks with biased functions. We demonstrate that
a previously used efficient algorithm for obtaining the frozen core, which
starts from the nodes with constant functions, fails when the number of inputs
per node exceeds 4. We present computer simulation data for the process of
formation of the frozen core and its robustness, and we show that several
important features of the data can be derived by using a mean-field
calculation.


  Swarming is an essential part of honeybee behaviour, wherein thousands of
bees cling onto each other to form a dense cluster that may be exposed to the
environment for several days. This cluster has the ability to maintain its core
temperature actively without a central controller, and raises the question of
how this is achieved. We suggest that the swarm cluster is akin to an active
porous structure whose functional requirement is to adjust to outside
conditions by varying its porosity to control its core temperature. Using a
continuum model that takes the form of a set of advection-diffusion equations
for heat transfer in a mobile porous medium, we show that the equalization of
an effective "behavioural pressure", which propagates information about the
ambient temperature through variations in density, leads to effective
thermoregulation. Our model extends and generalizes previous models by focusing
the question of mechanism on the form and role of the behavioural pressure, and
allows us to explain the vertical asymmetry of the cluster (as a consequence of
buoyancy driven flows), the ability of the cluster to overpack at low ambient
temperatures without breaking up at high ambient temperatures, and the relative
insensitivity to large variations in the ambient temperature. Finally, our
theory makes testable hypotheses for how the cluster bee density should respond
to externally imposed temperature inhomogeneities, and suggests strategies for
biomimetic thermoregulation.


  Following recent shootings in the USA, a debate has erupted, one side
favoring stricter gun control, the other promoting protection through more
weapons. We provide a scientific foundation to inform this debate, based on
mathematical, epidemiological models that quantify the dependence of
firearm-related death rates of people on gun policies. We assume a shooter
attacking a single individual or a crowd. Two strategies can minimize deaths in
the model, depending on parameters: either a ban of private firearms
possession, or a policy allowing the general population to carry guns. In
particular, the outcome depends on the fraction of offenders that illegally
possess a gun, on the degree of protection provided by gun ownership, and on
the fraction of the population who take up their right to own a gun and carry
it with them when attacked, parameters that can be estimated from statistical
data. With the measured parameters, the model suggests that if the gun law is
enforced at a level similar to that in the United Kingdom, gun-related deaths
are minimized if private possession of firearms is banned. If such a policy is
not practical or possible due to constitutional or cultural constraints, the
model and parameter estimation indicate that a partial reduction in firearm
availability can lead to a reduction in gun-induced death rates, even if they
are not minimized. Most importantly, our analysis identifies the crucial
parameters that determine which policy reduces the death rates, providing
guidance for future statistical studies that will be necessary for more refined
quantitative predictions.


  Development is a complex phenomenon where the forming phenotype interacts
with genetic and environmental inputs. Teeth are an important model for
developmental studies and their development has been thoroughly investigated.
However, despite of a large literature on the genetics of dental development,
no studies have focused yet on the environmental influences on dental
morphology. Here, we aim to test whether and to what extent the environment
plays a role in producing morphological variation in human teeth. We selected a
sample of modern human skulls and we used dental enamel hypoplasia as an
environmental stress marker in order to identify two groups with different
stress level, referred to as SG ("stressed" group) and NSG ("non-stressed"
group). We collected data on the occurrence and the relative development of 15
morphological traits (5 for each molar) on upper molars using a standard
methodology commonly used in dental anthropology (ASU-DAS system) and then we
compared the frequencies of the traits in the two groups. Overall, the results
suggest that: (i) stressors like malnutrition and/or systemic diseases have a
significant effect on upper molar morphology; (ii) stress generates a
developmental response which increases the morphological variability of the SG;
(iii) the increase of variability is directional, since individuals belonging
to the SG have increased cusps dimensions and number. These results are
consistent with the expectations of the morphodynamic model of dental
development (Jernvall, 1995).


  We develop a model of amoeboid cell motility based on active gel theory.
Modeling the motile apparatus of a eukaryotic cell as a confined layer of
finite length of poroelastic active gel permeated by a solvent, we first show
that, due to active stress and gel turnover, an initially static and
homogeneous layer can undergo a contractile-type instability to a polarized
moving state in which the rear is enriched in gel polymer. This agrees
qualitatively with motile cells containing an actomyosin-rich uropod at their
rear. We find that the gel layer settles into a steadily moving, inhomogeneous
state at long times, sustained by a balance between contractility and filament
turnover. In addition, our model predicts an optimal value of the
gel-susbstrate adhesion leading to maximum layer speed, in agreement with cell
motility assays. The model may be relevant to motility of cells translocating
in complex, confining environments that can be mimicked experimentally by cell
migration through microchannels.