  Free-riding on a joint venture bears the risk of losing personal endowment as
the group may fail to reach the collective target due to insufficient
contributions. A collective-risk social dilemma emerges, which we here study in
the realm of the spatial public goods game with group-performance-dependent
risk levels. Instead of using an overall fixed value, we update the risk level
in each group based on the difference between the actual contributions and the
declared target. A single parameter interpolates between a step-like risk
function and virtual irrelevance of the group's performance in averting the
failure, thus bridging the two extremes constituting maximal and minimal
feedback. We show that stronger feedback between group performance and risk
level is in general more favorable for the successful evolution of public
cooperation, yet only if the collective target to be reached is moderate.
Paradoxically, if the goals are overambitious, intermediate feedback strengths
yield optimal conditions for cooperation. This can be explained by the
propagation of players that employ identical strategies but experience
different individual success while trying to cope with the collective-risk
social dilemma.


  A simple model for the force-dependent unwinding and rewinding rates of the
nucleosome inner turn is constructed and quantitatively compared to the results
of recent measurements [A. H. Mack et al., J. Mol. Biol. 423, 687 (2012)].
First, a coarse-grained model for the histone-DNA free energy landscape that
incorporates both an elastic free energy barrier and specific histone-DNA bonds
is developed. Next, a theoretical expression for the rate of transitions across
a piecewise linear free energy landscape with multiple minima and maxima is
presented. Then, the model free energy landscape, approximated as a piecewise
linear function, and the theoretical expression for the transition rates are
combined to construct a model for the force-dependent unwinding and re-winding
rates of the nucleosome inner turn. Least-mean-squares fitting of the model
rates to the rates observed in recent experiments rates demonstrates that this
model is able to well describe the force-dependent unwinding and rewinding
rates of the nucleosome inner turn, observed in the recent experiments, except
at the highest forces studied, where an additional ad hoc term is required to
describe the data, which may be interpreted as an indication of an alternate
high-force nucleosome disassembly pathway, that bypasses simple unwinding. The
good agreement between the measurements and the model at lower forces
demonstrates that both specific histone-DNA contacts and an elastic free energy
barrier play essential roles for nucleosome winding and unwinding, and
quantifies their relative contributions.


  High-throughput sequencing allows the detection and quantification of
frequencies of somatic single nucleotide variants (SNV) in heterogeneous tumor
cell populations. In some cases, the evolutionary history and population
frequency of the subclonal lineages of tumor cells present in the sample can be
reconstructed from these SNV frequency measurements. However, automated methods
to do this reconstruction are not available and the conditions under which
reconstruction is possible have not been described.
  We describe the conditions under which the evolutionary history can be
uniquely reconstructed from SNV frequencies from single or multiple samples
from the tumor population and we introduce a new statistical model, PhyloSub,
that infers the phylogeny and genotype of the major subclonal lineages
represented in the population of cancer cells. It uses a Bayesian nonparametric
prior over trees that groups SNVs into major subclonal lineages and
automatically estimates the number of lineages and their ancestry. We sample
from the joint posterior distribution over trees to identify evolutionary
histories and cell population frequencies that have the highest probability of
generating the observed SNV frequency data. When multiple phylogenies are
consistent with a given set of SNV frequencies, PhyloSub represents the
uncertainty in the tumor phylogeny using a partial order plot. Experiments on a
simulated dataset and two real datasets comprising tumor samples from acute
myeloid leukemia and chronic lymphocytic leukemia patients demonstrate that
PhyloSub can infer both linear (or chain) and branching lineages and its
inferences are in good agreement with ground truth, where it is available.


  MicroRNAs (miRNAs) are small RNA molecules composed of 19-22 nt, which play
important regulatory roles in post-transcriptional gene regulation by
inhibiting the translation of the mRNA into proteins or otherwise cleaving the
target mRNA. Inferring miRNA targets provides useful information for
understanding the roles of miRNA in biological processes that are potentially
involved in complex diseases. Statistical methodologies for point estimation,
such as the Least Absolute Shrinkage and Selection Operator (LASSO) algorithm,
have been proposed to identify the interactions of miRNA and mRNA based on
sequence and expression data. In this paper, we propose using the Bayesian
LASSO (BLASSO) and the non-negative Bayesian LASSO (nBLASSO) to analyse the
interactions between miRNA and mRNA using expression data. The proposed
Bayesian methods explore the posterior distributions for those parameters
required to model the miRNA-mRNA interactions. These approaches can be used to
observe the inferred effects of the miRNAs on the targets by plotting the
posterior distributions of those parameters. For comparison purposes, the Least
Squares Regression (LSR), Ridge Regression (RR), LASSO, non-negative LASSO
(nLASSO), and the proposed Bayesian approaches were applied to four public
datasets. We concluded that nLASSO and nBLASSO perform best in terms of
sensitivity and specificity. Compared to the point estimate algorithms, which
only provide single estimates for those parameters, the Bayesian methods are
more meaningful and provide credible intervals, which take into account the
uncertainty of the inferred interactions of the miRNA and mRNA. Furthermore,
Bayesian methods naturally provide statistical significance to select
convincing inferred interactions, while point estimate algorithms require a
manually chosen threshold, which is less meaningful, to choose the possible
interactions.


  The high throughput and cost-effectiveness afforded by short-read sequencing
technologies, in principle, enable researchers to perform 16S rRNA profiling of
complex microbial communities at unprecedented depth and resolution. Existing
Illumina sequencing protocols are, however, limited by the fraction of the 16S
rRNA gene that is interrogated and therefore limit the resolution and quality
of the profiling. To address this, we present the design of a novel protocol
for shotgun Illumina sequencing of the bacterial 16S rRNA gene, optimized to
capture more than 90% of sequences in the Greengenes database and with nearly
twice the resolution of existing protocols. Using several in silico and
experimental datasets, we demonstrate that despite the presence of multiple
variable and conserved regions, the resulting shotgun sequences can be used to
accurately quantify the diversity of complex microbial communities. The
reconstruction of a significant fraction of the 16S rRNA gene also enabled high
precision (>90%) in species-level identification thereby opening up potential
application of this approach for clinical microbial characterization.


  Elucidating the spatial dynamic and core constituents of the microbial
communities found in association with arthropod hosts is of crucial importance
for insects that may vector human or agricultural pathogens. The hematophagous
Cimex lectularius, known as the common bed bug, has made a recent resurgence in
North America, as well as worldwide, potentially owing to increased travel and
resistance to insecticides. A comprehensive survey of the bed bug microbiome
has not been performed to date, nor has an assessment of the spatial dynamics
of its microbiome. Here we present a survey of bed bug microbial communities by
amplifying the V4-V6 hypervariable region of the 16S rDNA gene region followed
by 454 Titanium sequencing using 31 individuals from eight natural populations
collected from residences in Cincinnati, OH. Across all samples, 97% of the
microbial community is made up of two dominant OTUs identified as the
{\alpha}-proteobacterium Wolbachia and an unnamed {\gamma}-proteobacterium from
the Enterobacteriaceae. Microbial communities varied among host populations for
measures of community diversity and exhibited significant population structure.
We also uncovered a strong negative correlation in the abundance of the two
dominant OTUs, suggesting they may fulfill similar roles as nutritional
mutualists. This broad survey represents the most comprehensive assessment, to
date, of the microbes that associate with bed bugs, and uncovers evidence for
potential antagonism between the two dominant members of the bed bug
microbiome.


  Biochemical signaling pathways can be insulated from impedance and
competition effects through enzymatic "futile cycles" which consume energy,
typically in the form of ATP. We hypothesize that better insulation necessarily
requires higher energy consumption, and provide evidence, through the
computational analysis of a simplified physical model, to support this
hypothesis.


  Graphical modelling has a long history in statistics as a tool for the
analysis of multivariate data, starting from Wright's path analysis and Gibbs'
applications to statistical physics at the beginning of the last century. In
its modern form, it was pioneered by Lauritzen and Wermuth and Pearl in the
1980s, and has since found applications in fields as diverse as bioinformatics,
customer satisfaction surveys and weather forecasts.
  Genetics and systems biology are unique among these fields in the dimension
of the data sets they study, which often contain several hundreds of variables
and only a few tens or hundreds of observations. This raises problems in both
computational complexity and the statistical significance of the resulting
networks, collectively known as the "curse of dimensionality". Furthermore, the
data themselves are difficult to model correctly due to the limited
understanding of the underlying mechanisms. In the following, we will
illustrate how such challenges affect practical graphical modelling and some
possible solutions.


  Chromosomal aberrations are very frequent in leukemias and several recurring
mutations capable of malignant transformation have been described. These
mutations usually occur in hematopoietic stem cells (HSC), transforming them
into leukemia stem cells. NUP98 gene translocations are an example of such
chromosomal aberrations; these translocations produce a fusion protein
containing the N-terminal portion of Nup98 and the C-terminal of a fusion
partner. Over 75% of Nup98 fusions can interact with chromatin, and lead to
changes in gene expression. Therefore, I hypothesize that nup98 fusions act as
rogue transcriptional regulators in the cell. Collecting previously published
gene expression data (microarray) from HSCs expressing Nup98 fusions, we can
generate data to corroborate this hypothesis. Several different fusions affect
the expression of similar genes; these are involved in a few biological
processes in the cell: embryonic development, immune system formation and
chromatin organization. Deregulated genes also present similar transcription
factor binding sites in their regulatory regions. These putative regulatory
transcription factors are highly interconnected through protein-protein
interactions and transcriptional regulation among themselves, and they have
important roles in cell cycle regulation, embryonic development, hematopoiesis,
apoptosis and chromatin modification.


  In this paper we present a discrete dynamical population modeling of invasive
species, with reference to the swamp crayfish Procambarus clarkii. Since this
species can cause environmental damage of various kinds, it is necessary to
evaluate its expected in not yet infested areas. A structured discrete model is
built, taking into account all biological information we were able to find,
including the environmental variability implemented by means of stochastic
parameters (coefficients of fertility, death, etc.). This model is based on a
structure with 7 age classes, i.e. a Leslie mathematical population modeling
type and it is calibrated with laboratory data provided by the Department of
Evolutionary Biology (DEB) of Florence (Italy). The model presents many
interesting aspects: the population has a high initial growth, then it
stabilizes similarly to the logistic growth, but then it exhibits oscillations
(a kind of limit-cycle attractor in the phase plane). The sensitivity analysis
shows a good resilience of the model and, for low values of reproductive female
fraction, the fluctuations may eventually lead to the extinction of the
species: this fact might be exploited as a controlling factor. Moreover, the
probability of extinction is valuated with an inverse Gaussian that indicates a
high resilience of the species, confirmed by experimental data and field
observation: this species has diffused in Italy since 1989 and it has shown a
natural tendency to grow. Finally, the spatial mobility is introduced in the
model, simulating the movement of the crayfishes in a virtual lake of
elliptical form by means of simple cinematic rules encouraging the movement
towards the banks of the catchment (as it happens in reality) while a random
walk is imposed when the banks are reached.


  F$_\mathrm{o}$F$_1$-ATP synthase is a factory for synthesizing ATP in
virtually all cells. Its core machinery is the subcomplex F$_1$-motor
(F$_1$-ATPase) and performs the reversible mechanochemical coupling. Isolated
F$_1$-motor hydrolyzes ATP, which is accompanied by unidirectional rotation of
its central $\gamma$-shaft. When a strong opposing torque is imposed, the
$\gamma$-shaft rotates in the opposite direction and drives the F$_1$-motor to
synthesize ATP. This mechanical-to-chemical free-energy transduction is the
final and central step of the multistep cellular ATP-synthetic pathway. Here,
we determined the amount of mechanical work exploited by the F$_1$-motor to
synthesize an ATP molecule during forced rotations using methodology combining
a nonequilibrium theory and single molecule measurements of responses to
external torque. We found that the internal dissipation of the motor is
negligible even during rotations far from a quasistatic process.


  Evolutionary dynamics and patterns of molecular evolution are strongly
influenced by selection on linked regions of the genome, but our quantitative
understanding of these effects remains incomplete. Recent work has focused on
predicting the distribution of fitness within an evolving population, and this
forms the basis for several methods that leverage the fitness distribution to
predict the patterns of genetic diversity when selection is strong. However, in
weakly selected populations random fluctuations due to genetic drift are more
severe, and neither the distribution of fitness nor the sequence diversity
within the population are well understood. Here, we briefly review the
motivations behind the fitness-distribution picture, and summarize the general
approaches that have been used to analyze this distribution in the
strong-selection regime. We then extend these approaches to the case of weak
selection, by outlining a perturbative treatment of selection at a large number
of linked sites. This allows us to quantify the stochastic behavior of the
fitness distribution and yields exact analytical predictions for the sequence
diversity and substitution rate in the limit that selection is weak.


  One of application that needs high performance computing resources is
molecular d ynamic. There is some software available that perform molecular
dynamic, one of these is a well known GROMACS. Our previous experiment
simulating molecular dynamics of Indonesian grown herbal compounds show
sufficient speed up on 32 n odes Cluster computing environment. In order to
obtain a reliable simulation, one usually needs to run the experiment on the
scale of hundred nodes. But this is expensive to develop and maintain. Since
the invention of Graphical Processing Units that is also useful for general
programming, many applications have been developed to run on this. This paper
reports our experiments that evaluate the performance of GROMACS that runs on
two different environment, Cluster computing resources and GPU based PCs. We
run the experiment on BRV-1 and REM2 compounds. Four different GPUs are
installed on the same type of PCs of quad cores; they are Gefore GTS 250, GTX
465, GTX 470 and Quadro 4000. We build a cluster of 16 nodes based on these
four quad cores PCs. The preliminary experiment shows that those run on GTX 470
is the best among the other type of GPUs and as well as the cluster computing
resource. A speed up around 11 and 12 is gained, while the cost of computer
with GPU is only about 25 percent that of Cluster we built.


  The level of supercoiling in the chromosome can affect gene expression. To
clarify the basis of supercoiling sensitivity, we analyzed the structural
features of nucleotide sequences in the vicinity of promoters for the genes
with expression enhanced and decreased in response to loss of chromosomal
supercoiling in E. coli. Fourier analysis of promoter sequences for
supercoiling-sensitive genes reveals the tendency in selection of sequences
with helical periodicities close to 10 nt for relaxation-induced genes and to
11 nt for relaxation-repressed genes. The helical periodicities in the subsets
of promoters recognized by RNA polymerase with different sigma factors were
also studied. A special procedure was developed for study of correlations
between the intensities of periodicities in promoter sequences and the
expression levels of corresponding genes. Significant correlations of
expression with the AT content and with AT periodicities about 10, 11, and 50
nt indicate their role in regulation of supercoiling-sensitive genes.


  Models of infectious diseases have been developed since the first half of the
twentieth century. Most models haven't considered the role that emotional
factors of the individual may play on the population's behavioral adaptation
during the spread of a pandemic disease. Considering that local interactions
among individuals generate patterns that -at a large scale- govern the action
of masses, we have studied the behavioral adaptation of a population induced by
the spread of an infectious disease. Therefore, we have developed a rule-based
model of a hypothetical zombie outbreak, written in Kappa language, and
simulated using Guillespie's stochastic approach. Our study addresses the
specificity and heterogeneity of the system at the individual level, a highly
desirable characteristic, mostly overlooked in classic epidemic models.
Together with the basic elements of a typical epidemiological model, our model
includes an individual representation of the disease progression and the
traveling of agents among cities being affected. It also introduces an
approximation to measure the effect of panic in the population as a function of
the individual situational awareness. In addition, the effect of two possible
countermeasures to overcome the zombie threat is considered: the availability
of medical treatment and the deployment of special armed forces. However, due
to the special characteristics of this hypothetical infectious disease, even
using exaggerated numbers of countermeasures, only a small percentage of the
population can be saved at the end of the simulations. As expected from a
rule-based model approach, the global dynamics of our model resulted primarily
governed by the mechanistic description of local interactions occurring at the
individual level. As a whole, people's situational awareness resulted essential
to modulate the inner dynamics of the system.


  Weak purifying selection, acting on many linked mutations, may play a major
role in shaping patterns of molecular evolution in natural populations. Yet
efforts to infer these effects from DNA sequence data are limited by our
incomplete understanding of weak selection on local genomic scales. Here, we
demonstrate a natural symmetry between weak and strong selection, in which the
effects of many weakly selected mutations on patterns of molecular evolution
are equivalent to a smaller number of more strongly selected mutations. By
introducing a coarse-grained "effective selection coefficient," we derive an
explicit mapping between weakly selected populations and their strongly
selected counterparts, which allows us to make accurate and efficient
predictions across the full range of selection strengths. This suggests that an
effective selection coefficient and effective mutation rate --- not an
effective population size --- is the most accurate summary of the effects of
selection over locally linked regions. Moreover, this correspondence places
fundamental limits on our ability to resolve the effects of weak selection from
contemporary sequence data alone.


  We consider multiple diseases spreading in a static Configuration Model
network. We make standard assumptions that infection transmits from neighbor to
neighbor at a disease-specific rate and infected individuals recover at a
disease-specific rate. Infection by one disease confers immediate and permanent
immunity to infection by any disease. Under these assumptions, we find a
simple, low-dimensional ordinary differential equations model which captures
the global dynamics of the infection. The dynamics depend strongly on initial
conditions. Although we motivate this article with infectious disease, the
model may be adapted to the spread of other infectious agents such as competing
political beliefs, rumors, or adoption of new technologies if these are
influenced by contacts. As an example, we demonstrate how to model an
infectious disease which can be prevented by a behavior change.


  Electrostatic interaction is the driving force for the encapsulation by virus
coat proteins of nanoparticles such as quantum dots, gold particles and
magnetic beads for, e.g., imaging and therapeutic purposes. In recent
experimental work, Daniel et al. [ACS Nano 4 (2010), 3853-3860] found the
encapsulation efficiency to sensitively depend on the interplay between the
surface charge density of negatively charged gold nanoparticles and the number
of positive charges on the RNA binding domains of the proteins. Surprisingly,
these experiments reveal that despite the highly cooperative nature of the
co-assembly at low pH, the efficiency of encapsulation is a gradual function of
their surface charge density. We present a simple all-or-nothing mass action
law combined with an electrostatic interaction model to explain the
experiments. We find quantitative agreement with experimental observations,
supporting the existence of a natural statistical charge distribution between
nanoparticles.


  An important task in quantitative biology is to understand the role of
stochasticity in biochemical regulation. Here, as an extension of our recent
work [Phys. Rev. Lett. 107, 148101 (2011)], we study how input fluctuations
affect the stochastic dynamics of a simple biological switch. In our model, the
on transition rate of the switch is directly regulated by a noisy input signal,
which is described as a nonnegative mean-reverting diffusion process. This
continuous process can be a good approximation of the discrete birth-death
process and is much more analytically tractable. Within this new setup, we
apply the Feynman-Kac theorem to investigate the statistical features of the
output switching dynamics. Consistent with our previous findings, the input
noise is found to effectively suppress the input-dependent transitions. We show
analytically that this effect becomes significant when the input signal
fluctuates greatly in amplitude and reverts slowly to its mean.


  Source-sink systems are metapopulations of habitat patches with different,
and possibly temporally varying, habitat qualities, which are commonly used in
ecology to study the fate of spatially extended natural populations. We propose
new techniques that allow to disentangle the respective contributions of
demography and dispersal to the dynamics and fate of a single species in a
source-sink metapopulation. Our approach is valid for a general class of
stochastic, individual-based, stepping-stone models, with density-independent
demography and dispersal, provided the metapopulation is finite or else enjoys
some transitivity property. We provide 1) a simple criterion of persistence, by
studying the motion of a single random disperser until it returns to its
initial position; 2) a joint characterization of the long-term growth rate and
of the asymptotic occupancy frequencies of the ancestral lineage of a random
survivor, by using large deviations theory. Both techniques yield formulae
decoupling demography and dispersal, and can be adapted to the case of periodic
or random environments, where habitat qualities are autocorrelated in space and
possibly in time. In this last case, we display examples of coupled
time-averaged sinks allowing survival, as was previously known in the absence
of demographic stochasticity for fully mixing (Jansen and Yoshimura, 1998) and
even partially mixing (Evans et al., 2012; Schreiber, 2010) metapopulations.


  Emergence of robustness in biological networks is a paramount feature of
evolving organisms, but a study of this property in vivo, for any level of
representation such as Genetic, Metabolic, or Neuronal Networks, is a very hard
challenge. In the case of Genetic Networks, mathematical models have been used
in this context to provide insights on their robustness, but even in relatively
simple formulations, such as Boolean Networks (BN), it might not be feasible to
compute some measures for large system sizes. We describe in this work a Monte
Carlo approach to calculate the size of the largest basin of attraction of a
BN, which is intrinsically associated with its robustness, that can be used
regardless the network size. We show the stability of our method through
finite-size analysis and validate it with a full search on small networks.


  Compartmentalization of self-replicating molecules (templates) in protocells
is a necessary step towards the evolution of modern cells. However, coexistence
between distinct template types inside a protocell can be achieved only if
there is a selective pressure favoring protocells with a mixed template
composition. Here we study analytically a group selection model for the
coexistence between two template types using the diffusion approximation of
population genetics. The model combines competition at the template and
protocell levels as well as genetic drift inside protocells. At the steady
state, we find a continuous phase transition separating the coexistence and
segregation regimes, with the order parameter vanishing linearly with the
distance to the critical point. In addition, we derive explicit analytical
expressions for the critical steady-state probability density of protocell
compositions.


  Over the past several years, DNA sequencing has emerged as one of the driving
forces in life-sciences, paving the way for affordable and accurate whole
genome sequencing. As genomes represent the entirety of an organism's
hereditary information, the availability of complete human genomes prompts a
wide range of revolutionary applications. The hope for improving modern
healthcare and better understanding the human genome propels many interesting
and challenging research frontiers. Unfortunately, however, the proliferation
of human genomes amplifies worrisome privacy concerns, since a genome
represents a treasure trove of highly personal and sensitive information. In
this article, we provide an overview of positive results and biomedical
advances in the field, and discuss privacy issues associated with human genomic
information. Finally, we survey available privacy-enhancing technologies and
list a number of open research challenges.


  Sensory rhodopsin II (SRII) is a seven helix protein that belongs to the
rhodopsin protein family. Light induced conformational changes govern SRII's
function. These changes are related to the photo cycle of the protein that is
comprised of various metastable states. After the completion of this cycle the
protein returns to its ground state. Mutational studies of key residues will
reveal the mechanism that underlies the function of the protein. The result
will allow us to determine key structures at various PHs involved in the photo
cycle of SRII and to understand the protein's function and mechanism.


  The goal of phylodynamics, an area on the intersection of phylogenetics and
population genetics, is to reconstruct population size dynamics from genetic
data. Recently, a series of nonparametric Bayesian methods have been proposed
for such demographic reconstructions. These methods rely on prior
specifications based on Gaussian processes and proceed by approximating the
posterior distribution of population size trajectories via Markov chain Monte
Carlo (MCMC) methods. In this paper, we adapt an integrated nested Laplace
approximation (INLA), a recently proposed approximate Bayesian inference for
latent Gaussian models, to the estimation of population size trajectories. We
show that when a genealogy of sampled individuals can be reliably estimated
from genetic data, INLA enjoys high accuracy and can replace MCMC entirely. We
demonstrate significant computational efficiency over the state-of-the-art MCMC
methods. We illustrate INLA-based population size inference using simulations
and genealogies of hepatitis C and human influenza viruses.


  One of the challenges of the proteomic analysis by 2D-gel is to visualize the
low abundance proteins, particularly those localized in organelles. An
additional problem with nuclear proteins lies in their strong interaction with
nuclear acids. Several experimental procedures have been tested to increase, in
the nuclear extract, the ratio of nuclear proteins compared to contaminant
proteins, and also to obtain reproducible conditions compatible with 2D-gel
electrophoresis. The NaCl procedure has been chosen. To test the interest of
this procedure, the nuclear protein expression profiles of macrophages and
dendritic cells have been compared with a proteomic approach by 2D-gel
electrophoresis. Delta 2D software and mass spectrometry analyses have allowed
pointing out some proteins of interest. We have chosen some of them, involved
in transcriptional regulation and/or chromatin structure for further
validations. The immunoblotting experiments have shown that most of observed
changes are due to post-translational modifications, thereby a exemplifying the
interest of the 2D gel approach. Finally, this approach allowed us to reach not
only high abundance nuclear proteins but also lower abundance proteins, such as
the HP1 proteins and reinforces the interest of using 2DE-gel in proteomics
because of its ability to visualize intact proteins with their modifications.


  Background: Much has been written about introduced rainbow trout
(Oncorhynchus mykiss) interbreeding and outcompeting cutthroat trout
(Oncorhynchus clarkii). However, the specific mechanisms by which rainbow trout
and their hybrids outcompete cutthroat trout have not been thoroughly explored,
and the published data is limited to lotic ecosystems. Materials and Methods:
Samples of rainbow trout and cutthroat trout were obtained from a lentic
ecosystem by angling. The total length and weight of each fish was measured and
the relative weight of each fish was computed (Anderson R.O., Neumann R.M.
1996. Length, Weight, and Associated Structural Indices, Pp. 447-481. In:
Murphy B.E. and Willis D.W. (eds.) Fisheries Techniques, second edition.
American Fisheries Society.), along with the mean and uncertainty in the mean
for each species. Data from an independent source (K.D. Carlander, 1969.
Handbook of Freshwater Fishery Biology, Volume One, Iowa University Press,
Ames.) was also used to generate mean weight-length curves, as well as 25th and
75th percentile curves for each species to allow further comparison. Results:
The mean relative weight of the rainbow trout was 72.5 (+/- 2.1); whereas, the
mean relative weight of the cutthroat trout was 101.0 (+/- 4.9). The rainbow
trout were thin; 80% weighed below the 25th percentile. The cutthroat trout
were plump; 86% weighed above the 75th percentile, and 29% were above the
heaviest recorded specimens at a given length in the Carlander (1969) data set.
Conclusion: This data casts doubt on the hypothesis that rainbow trout are
strong food competitors with cutthroat trout in lentic ecosystems. On the
contrary, in the lake under study, the cutthroat trout seem to be outcompeting
rainbow trout for the available food.


  Many growth processes lead to intriguing stochastic patterns and complex
fractal structures which exhibit local scale invariance properties. Such
structures can often be described effectively by space-time trajectories of
interacting particles, and their large scale behaviour depends on the overall
growth geometry. We establish an exact relation between statistical properties
of structures in uniformly expanding and fixed geometries, which preserves the
local scale invariance and is independent of other properties such as the
dimensionality. This relation generalizes standard conformal transformations as
the natural symmetry of self-affine growth processes. We illustrate our main
result numerically for various structures of coalescing L\'evy flights and
fractional Brownian motions, including also branching and finite particle
sizes. One of the main benefits of this new approach is a full, explicit
description of the asymptotic statistics in expanding domains, which are often
non-trivial and random due to amplification of initial fluctuations.


  Living cells often need to extract information from biochemical signals that
are noisy. We study how accurately cells can measure chemical concentrations
with signaling networks that are linear. For stationary signals of long
duration, they can reach, but not beat, the Berg-Purcell limit, which relies on
uniformly averaging in time the fluctuations in the input signal. For short
times or nonstationary signals, however, they can beat the Berg-Purcell limit,
by non-uniformly time-averaging the input. We derive the optimal weighting
function for time averaging and use it to provide the fundamental limit of
measuring chemical concentrations with linear signaling networks.


  The evolution of cooperation has been a perennial problem in evolutionary
biology because cooperation can be undermined by selfish cheaters who gain an
advantage in the short run, while compromising the long-term viability of the
population. Evolutionary game theory has shown that under certain conditions,
cooperation nonetheless evolves stably, for example if players have the
opportunity to punish cheaters that benefit from a public good yet refuse to
pay into the common pool. However, punishment has remained enigmatic because it
is costly, and difficult to maintain. On the other hand, cooperation emerges
naturally in the Public Goods game if the synergy of the public good (the
factor multiplying the public good investment) is sufficiently high. In terms
of this synergy parameter, the transition from defection to cooperation can be
viewed as a phase transition with the synergy as the critical parameter. We
show here that punishment reduces the critical value at which cooperation
occurs, but also creates the possibility of meta-stable phase transitions,
where populations can "tunnel" into the cooperating phase below the critical
value. At the same time, cooperating populations are unstable even above the
critical value, because a group of defectors that are large enough can
"nucleate" such a transition. We study the mean-field theoretical predictions
via agent-based simulations of finite populations using an evolutionary
approach where the decisions to cooperate or to punish are encoded genetically
in terms of evolvable probabilities. We recover the theoretical predictions and
demonstrate that the population shows hysteresis, as expected in systems that
exhibit super-heating and super-cooling. We conclude that punishment can
stabilize populations of cooperators below the critical point, but it is a
two-edged sword: it can also stabilize defectors above the critical point.


  As a widely used method in metabolic network studies, Monte-Carlo sampling in
the steady state flux space is known for its flexibility and convenience of
carrying out different purposes, simply by alternating constraints or objective
functions, or appending post processes. Recently the concept of a non-linear
constraint based on the second thermodynamic law, known as "Loop Law", is
challenging current sampling algorithms which will inevitably give rise to the
internal loops. A generalized method is proposed here to eradicate the
probability of the appearance of internal loops during sampling process. Based
on Artificial Centered Hit and Run (ACHR) method, each step of the new sampling
process will avoid entering "loop-forming" subspaces. This method has been
applied on the metabolic network of Helicobacter pylori with three different
objective functions: uniform sampling, optimizing biomass synthesis, optimizing
biomass synthesis efficiency over resources ingested. Comparison between
results from the new method and conventional ACHR method shows effective
elimination of loop fluxes without affecting non-loop fluxes.


  Phenotypic fluctuations and plasticity can generally affect the course of
evolution, a process known as the Baldwin effect. Several studies have recast
this effect and claimed that phenotypic plasticity acceler- ates evolutionary
rate (the Baldwin expediting effect); however, the validity of this claim is
still controversial. In this study, we investi- gate the evolutionary
population dynamics of a quantitative genetic model under a multi-peaked
fitness landscape, in order to evaluate the validity of the effect. We provide
analytical expressions for the evolutionary rate and average population
fitness. Our results indicate that under a multi-peaked fitness landscape,
phenotypic fluctuation always accelerates evolutionary rate, but it decreases
the average fit- ness. As an extreme case of the trade-off between the rate of
evolution and average fitness, phenotypic fluctuation is shown to accelerate
the error catastrophe, in which a population fails to sustain a high-fitness
peak. In the context of our findings, we discuss the role of phenotypic
plasticity in adaptive evolution.


  Basic body-part names (BBPNs) were defined as body-part names in Swadesh
basic 200 words. Non-Mayan cognates of Mayan (MY) BBPNs were extensively
searched for, by comparing with non-MY vocabulary, including ca.1300 basic
words of 82 AN languages listed by Tryon (1985), etc. Thus found cognates (CGs)
in non-MY are listed in Table 1, as classified by language groups to which most
similar cognates (MSCs) of MY BBPNs belong. CGs of MY are classified to 23
mutually unrelated CG-items, of which 17.5 CG-items have their MSCs in
Austronesian (AN), giving its closest similarity score (CSS), CSS(AN) = 17.5,
which consists of 10.33 MSCs in Formosan, 1.83 MSCs in Western
Malayo-Polynesian (W.MP), 0.33 in Central MP, 0.0 in SHWNG, and 5.0 in Oceanic
[i.e., CSS(FORM)= 10.33, CSS(W.MP) = 1.88, ..., CSS(OC)= 5.0]. These CSSs for
language (sub)groups are also listed in the underline portion of every section
of (Section1 - Section 6) in Table 1. Chi-squar test (degree of freedom = 1)
using [Eq 1] and [Eqs.2] revealed that MSCs of MY BBPNs are distributed in
Formosan in significantly higher frequency (P < 0.001) than in other subgroups
of AN, as well as than in non-AN languages. MY is thus concluded to have been
derived from Formosan of AN. Eskimo shows some BBPN similarities to FORM and
MY.


  Combining extensive molecular dynamics simulations of lipid bilayer systems
of varying chemical composition with single-trajectory analyses we
systematically elucidate the stochastic nature of the lipid motion. We observe
subdiffusion over more than four orders of magnitude in time, clearly
stretching into the sub-microsecond domain. The lipid motion delicately depends
on the lipid chemistry, the lipid phase, and especially on the presence of
cholesterol. We demonstrate that fractional Langevin equation motion
universally describes the lipid motion in all phases including the gel phase,
and in the presence of cholesterol. The results underline the relevance of
anomalous diffusion in lipid bilayers and the strong effects of the membrane
composition.


  The basic mechanics of evolution have been understood since Darwin. But
debate continues over whether macroevolutionary phenomena are driven primary by
the fitness structure of genotype space or by ecological interaction. In this
paper we propose a simple, abstract model capturing some key features of
fitness-landscape and ecological models of evolution. Our model describes
evolutionary dynamics in a high-dimensional, structured genotype space with a
significant role for interspecific interaction. We find some promising
qualitative similarity with the empirical facts about macroevolution, including
broadly distributed extinction sizes and realistic exploration of the genotype
space. The abstraction of our model permits numerous interpretations and
applications beyond macroevolution, including molecular evolution and
technological innovation.


  We combine incentive, adaptive, and time-scale dynamics to study
multipopulation dynamics on the simplex equipped with a large class of
Riemmanian metrics, simultaneously generalizing and extending many dynamics
commonly studied in dynamic game theory and evolutionary dynamics. Each
population has its own geometry, method of adaptation (incentive), and
time-scale (discrete, continuous, and others). Using an information-theoretic
measure of distance we give a widely-applicable Lyapunov result for the
dynamic. We include a wealth of examples leading up to and beyond the main
results.


  The adaptive evolution of large asexual populations is generally
characterized by competition between clones carrying different beneficial
mutations. This interference phenomenon slows down the adaptation speed and
makes the theoretical description of the dynamics more complex with respect to
the successional occurrence and fixation of beneficial mutations typical of
small populations. A simplified modeling framework considering multiple
beneficial mutations with equal and constant fitness advantage captures some of
the essential features of the actual complex dynamics, and some key predictions
from this model are verified in laboratory evolution experiments. However, in
these experiments the relative advantage of a beneficial mutation is generally
dependent on the genetic background. In particular, the general pattern is
that, as mutations in different loci accumulate, the relative advantage of new
mutations decreases, trend often referred to as "diminishing return" epistasis.
  In this paper, we propose a phenomenological model that generalizes the
fixed-advantage framework to include in a simple way this feature. To evaluate
the quantitative consequences of diminishing returns on the evolutionary
dynamics, we approach the model analytically as well as with direct
simulations. Finally, we show how the model parameters can be matched with data
from evolutionary experiments in order to infer the mean effect of epistasis
and derive order-of-magnitude estimates of the rate of beneficial mutations.
Applying this procedure to two experimental data sets gives values of the
beneficial mutation rate within the range of previous measurements.


  Living systems are typically characterized by irreversible processes. A
condition equivalent to the reversibility is the detailed balance, whose
absence is an obstacle for analytically solving ecological models. We revisit a
promising model with an elegant field-theoretic analytic solution and show that
the theoretical analysis is invalid because of an implicit assumption of
detailed balance. A signature of the difficulties is evident in the
inconsistencies appearing in the many-point correlation functions and in the
analytical formula for the species area relationship.


  Global mobility flow data are at the heart of spatial epidemiological models
used to predict infectious disease behavior but this wealth of data on human
mobility has been largely neglected by reconstructions of pathogen evolutionary
dynamics using viral genetic data. Although stochastic models of viral
evolution may potentially be informed by such data, a major challenge lies in
deciding which mobility processes are critical and to what extent they
contribute to shaping contemporaneous distributions of pathogen diversity.
Here, we develop a framework to integrate predictors of viral diffusion with
phylogeographic inference and estimate human influenza H3N2 migration history
while simultaneously testing and quantifying the factors that underly it. We
provide evidence for air travel governing the global dynamics of human
influenza whereas other processes act at a more local scale.


  We propose a model of parameter learning for signal transduction, where the
objective function is defined by signal transmission efficiency. We apply this
to learn kinetic rates as a form of evolutionary learning, and look for
parameters which satisfy the objective. This is a novel approach compared to
the usual technique of adjusting parameters only on the basis of experimental
data. The resulting model is self-organizing, i.e. perturbations in protein
concentrations or changes in extracellular signaling will automatically lead to
adaptation. We systematically perturb protein concentrations and observe the
response of the system. We find compensatory or co-regulation of protein
expression levels. In a novel experiment, we alter the distribution of
extracellular signaling, and observe adaptation based on optimizing signal
transmission. We also discuss the relationship between signaling with and
without transients. Signaling by transients may involve maximization of signal
transmission efficiency for the peak response, but a minimization in
steady-state responses. With an appropriate objective function, this can also
be achieved by concentration adjustment. Self-organizing systems may be
predictive of unwanted drug interference effects, since they aim to mimic
complex cellular adaptation in a unified way.


  The frequent outbreak of severe foodborne diseases warns of a potential
threat that the global trade networks could spread fatal pathogens. The global
trade network is a typical overlay network, which compounds multiple standalone
trade networks representing the transmission of a single product and connecting
the same set of countries and territories through their own set of trade
interactions. Although the epidemic dynamic implications of overlay networks
have been debated in recent studies, some general answers for the overlay of
multiple and diverse standalone networks remain elusive, especially the
relationship between the heterogeneity and diversity of a set of standalone
networks and the behavior of the overlay network. In this paper, we establish a
general analysis framework for multiple overlay networks based on diversity
theory. The framework could reveal the critical epidemic mechanisms beyond
overlay processes. Applying the framework to global trade networks, we found
that, although the distribution of connectivity of standalone trade networks
was highly heterogeneous, epidemic behavior on overlay networks is more
dependent on cooperation among standalone trade networks rather than on a few
high-connectivity networks as the general property of complex systems with
heterogeneous distribution. Moreover, the analysis of overlay trade networks
related to 7 real pathogens also suggested that epidemic behavior is not
controlled by high-connectivity goods but that the actual compound mode of
overlay trade networks plays a critical role in spreading pathogens. Finally,
we study the influence of cooperation mechanisms on the stability of overlay
networks and on the control of global epidemics. The framework provides a
general tool to study different problems on overlay networks.


  We perform a spatially resolved simulation study of an AND gate based on DNA
strand displacement using several lengths of the toehold and the adjacent
domains. DNA strands are modelled using a coarse-grained dynamic bonding model
{[}C. Svaneborg, Comp. Phys. Comm. 183, 1793 (2012){]}. We observe a complex
transition path from the initial state to the final state of the AND gate. This
path is strongly influenced by non-ideal effects due to transient bubbles
revealing undesired toeholds and thermal melting of whole strands. We have also
characterized the bound and unbound kinetics of single strands, and in
particular the kinetics of the total AND operation and the three distinct
distinct DNA transitions that it is based on. We observe a exponential kinetic
dependence on the toehold length of the competitive displacement operation, but
that the gate operation time is only weakly dependent on both the toehold and
adjacent domain length. Our gate displays excellent logical fidelity in three
input states, and quite poor fidelity in the fourth input state. This
illustrates how non-ideality can have very selective effects on fidelity.
Simulations and detailed analysis such as those presented here provide
molecular insights into strand displacement computation, that can be also be
expected in chemical implementations.


  Biological networks have two modes. The first mode is static: a network is a
passage on which something flows. The second mode is dynamic: a network is a
pattern constructed by gluing functions of entities constituting the network.
In this paper, first we discuss that these two modes can be associated with the
category theoretic duality (adjunction) and derive a natural network structure
(a path notion) for each mode by appealing to the category theoretic
universality. The path notion corresponding to the static mode is just the
usual directed path. The path notion for the dynamic mode is called lateral
path which is the alternating path considered on the set of arcs. Their general
functionalities in a network are transport and coherence, respectively. Second,
we introduce a betweenness centrality of arcs for each mode and see how the two
modes are embedded in various real biological network data. We find that there
is a trade-off relationship between the two centralities: if the value of one
is large then the value of the other is small. This can be seen as a kind of
division of labor in a network into transport on the network and coherence of
the network. Finally, we propose an optimization model of networks based on a
quality function involving intensities of the two modes in order to see how
networks with the above trade-off relationship can emerge through evolution. We
show that the trade-off relationship can be observed in the evolved networks
only when the dynamic mode is dominant in the quality function by numerical
simulations. We also show that the evolved networks have features qualitatively
similar to real biological networks by standard complex network analysis.


  We investigate the motion of active semiflexible filament with shape
kinematics and hydrodynamic interaction including. Three types of filament
motion are found: Translation, snaking and rotation. Change of flexibility will
induce instability of shape kinematics and further result in asymmetry of shape
kinematics respect to the motion of mass center, which are responsible to a
continuous-like transition from translation to snaking and a first-order-like
transition from snaking to rotation, respectively. Of particular interest, we
find that long-range hydrodynamic interaction is not necessary for filament
rotation, but can enhance remarkably the parameter region for its appearance.
This finding may provide an evidence that the experimentally found collective
rotation of active filaments is more likely to arise from the individual
property even without the long-range hydrodynamic interaction.


  In large asexual populations, multiple beneficial mutations arise in the
population, compete, interfere with each other, and accumulate on the same
genome, before any of them fix. The resulting dynamics, although studied by
many authors, is still not fully understood, fundamentally because the effects
of fluctuations due to the small numbers of the fittest individuals are large
even in enormous populations. In this paper, branching processes and various
asymptotic methods for analyzing the stochastic dynamics are further developed
and used to obtain information on fluctuations, time dependence, and the
distributions of sizes of subpopulations, jumps in the mean fitness, and other
properties. The focus is on the behavior of a broad class of models: those with
a distribution of selective advantages of available beneficial mutations that
falls off more rapidly than exponentially. For such distributions, many aspects
of the dynamics are universal - quantitatively so for extremely large
populations. On the most important time scale that controls coalescent
properties and fluctuations of the speed, the dynamics is reduced to a simple
stochastic model that couples the peak and the high-fitness "nose" of the
fitness distribution. Extensions to other models and distributions of available
mutations are discussed briefly.


  We developed NameMyGene, a web tool and a stand alone program to easily
generate putative family-based names for small RNA sequences so that
laboratories can easily organize, analyze, and observe patterns from, the
massive amount of data generated by next-generation sequencers. NameMyGene,
also applicable to other emerging methods such as RNA-Seq, and Chip-Seq, solely
uses the input small RNA sequence and does not require any additional data such
as other sequence data sets. The web server and software is freely available
(http://www.johnlab.org/NameMyGene) and is based on Java to ensure platform
independency.


  We study the combined impact of a colored environmental noise and demographic
noise on the extinction risk of a long-lived and well-mixed isolated stochastic
population which exhibits the Allee effect. The environmental noise modulates
the population birth and death rates. Assuming that the Allee effect is strong,
and the environmental noise is positively correlated and Gaussian, we derive a
Fokker-Planck equation for the joint probability distribution of the population
sizes and environmental fluctuations. In WKB approximation this equation
reduces to an effective two-dimensional Hamiltonian mechanics, where the most
likely path to extinction and the most likely environmental fluctuation are
encoded in an instanton-like trajectory in the phase space. The mean time to
extinction (MTE) is related to the mechanical action along this trajectory. We
obtain new analytic results for short-correlated, long-correlated and
relatively weak environmental noise. The population-size dependence of the MTE
changes from exponential for weak environmental noise to no dependence for
strong noise, implying a greatly increased extinction risk. The theory is
readily extendable to population switches between different metastable states,
and to stochastic population explosion, due to a combined action of demographic
and environmental noise.


  Developmental constraints have been postulated to limit the space of feasible
phenotypes and thus shape animal evolution. These constraints have been
suggested to be the strongest during either early or mid-embryogenesis, which
corresponds to the early conservation model or the hourglass model,
respectively. Conflicting results have been reported, but in recent studies of
animal transcriptomes the hourglass model has been favored. Studies usually
report descriptive statistics calculated for all genes over all developmental
time points. This introduces dependencies between the sets of compared genes,
and may lead to biased results. Here we overcome this problem using an
alternative modular analysis. We used the Iterative Signature Algorithm to
identify distinct modules of genes co-expressed specifically in consecutive
stages of zebrafish development. We then performed a detailed comparison of
several gene properties between modules, allowing for a less biased and more
powerful analysis. Notably, our analysis corroborated the hourglass pattern
only at the regulatory level, with sequences of regulatory regions being most
conserved for genes expressed in mid-development, but not at the level of gene
sequence, age or expression, in contrast to some previous studies. The early
conservation model was supported with gene duplication and birth that were the
most rare for genes expressed in early development. Finally, for all gene
properties we observed the least conservation for genes expressed in late
development or adult, consistent with both models. Overall, with the modular
approach, we showed that different levels of molecular evolution follow
different patterns of developmental constraints. Thus both models are valid,
but with respect to different genomic features.


  The gene networks that comprise the circadian clock modulate biological
function across a range of scales, from gene expression to performance and
adaptive behaviour. The clock functions by generating endogenous rhythms that
can be entrained to the external 24-h day?night cycle, enabling organisms to
optimally time biochemical processes relative to dawn and dusk. In recent
years, computational models based on differential equations have become useful
tools for dissecting and quantifying the complex regulatory relationships
underlying the clock's oscillatory dynamics. However, optimizing the large
parameter sets characteristic of these models places intense demands on both
computational and experimental resources, limiting the scope of in silico
studies. Here, we develop an approach based on Boolean logic that dramatically
reduces the parametrization, making the state and parameter spaces finite and
tractable. We introduce efficient methods for fitting Boolean models to
molecular data, successfully demonstrating their application to synthetic time
courses generated by a number of established clock models, as well as
experimental expression levels measured using luciferase imaging. Our results
indicate that despite their relative simplicity, logic models can (i) simulate
circadian oscillations with the correct, experimentally observed phase
relationships among genes and (ii) flexibly entrain to light stimuli,
reproducing the complex responses to variations in daylength generated by more
detailed differential equation formulations. Our work also demonstrates that
logic models have sufficient predictive power to identify optimal regulatory
structures from experimental data. By presenting the first Boolean models of
circadian circuits together with general techniques for their optimization, we
hope to establish a new framework for the systematic modelling of complex
clocks.


  We propose a new method for molecular dynamics and Monte Carlo simulations,
which is referred to as the replica-permutation method (RPM), to realize more
efficient sampling than the replica-exchange method (REM).In RPM not only
exchanges between two replicas but also permutations among more than two
replicas are performed. Furthermore, instead of the Metropolis algorithm, the
Suwa-Todo algorithm is employed for replica-permutation trials to minimize its
rejection ratio. We applied RPM to particles in a double-well potential energy,
Met-enkephalin in vacuum, and a C-peptide analog of ribonuclease A in explicit
water. For a comparison purposes, replica-exchange molecular dynamics
simulations were also performed. As a result, RPM sampled not only the
temperature space but also the conformational space more efficiently than REM
for all systems. From our simulations of C-peptide, we obtained the
\alpha-helix structure with salt-bridges between Gly2 and Arg10 which is known
in experiments. Calculating its free-energy landscape, the folding pathway was
revealed from an extended structure to the \alpha-helix structure with the
salt-bridges. We found that the folding pathway consists of the two steps: The
first step is the "salt-bridge formation step", and the second step is the
"\alpha-helix formation step".


  Menzerath's law, the tendency of Z, the mean size of the parts, to decrease
as X, the number of parts, increases is found in language, music and genomes.
Recently, it has been argued that the presence of the law in genomes is an
inevitable consequence of the fact that Z = Y/X, which would imply that Z
scales with X as Z ~ 1/X. That scaling is a very particular case of
Menzerath-Altmann law that has been rejected by means of a correlation test
between X and Y in genomes, being X the number of chromosomes of a species, Y
its genome size in bases and Z the mean chromosome size. Here we review the
statistical foundations of that test and consider three non-parametric tests
based upon different correlation metrics and one parametric test to evaluate if
Z ~ 1/X in genomes. The most powerful test is a new non-parametric based upon
the correlation ratio, which is able to reject Z ~ 1/X in nine out of eleven
taxonomic groups and detect a borderline group. Rather than a fact, Z ~ 1/X is
a baseline that real genomes do not meet. The view of Menzerath-Altmann law as
inevitable is seriously flawed.


  The evidence that Flowering Locus T (FT) protein and its paralog Twin Sister
of FT, act as the long distance floral stimulus, or at least that they are part
of it in diverse plant species, has attracted much attention in recent years.
Studies to understand the physiological and molecular apparatuses that
integrate spatial and temporal signals to regulate developmental transition in
plants have occupied countless scientists and have resulted in an unmanageably
large amount of research data. Analysis of these data has helped to identify
multiple systemic florigenic and antiflorigenic regulators. This study gives an
overview of the recent research on gene products, phytohormones and other
metabolites that have been demonstrated to have florigenic or antiflorigenic
functions in plants.


  A collective-risk social dilemma implies that personal endowments will be
lost if contributions to the common pool within a group are too small. Failure
to reach the collective target thus has dire consequences for all group
members, independently of their strategies. Wanting to move away from
unfavorable locations is therefore all but surprising. Inspired by these
observations, we here propose and study a collective-risk social dilemma where
players are allowed to move if the collective failure becomes too probable.
More precisely, this so-called risk-driven migration is launched depending on
the difference between the actual contributions and the declared target.
Mobility therefore becomes an inherent property that is utilized in an entirely
self-organizing manner. We show that under these assumptions cooperation is
promoted much more effectively than under the action of manually determined
migration rates. For the latter, we in fact identify parameter regions where
the evolution of cooperation is incredibly inhibited. Moreover, we find
unexpected spatial patterns where cooperators that do not form compact clusters
outperform those that do, and where defectors are able to utilize strikingly
different ways of invasion. The presented results support the recently revealed
importance of percolation for the successful evolution of public cooperation,
while at the same time revealing surprisingly simple ways of self-organization
towards socially desirable states.


  Pulsed-field-gradient nuclear magnetic resonance (PFG-NMR) is used to obtain
the true hydrodynamic size of complexes of peptides with sodium dodecyl sulfate
SDS micelles. The peptide used in this study is a 19-residue antimicrobial
peptide, GAD-2. Two smaller dipeptides, alanine-glycine (Ala-Gly) and
tyrosine-leucine (Tyr-Leu), are used for comparison. We use PFG-NMR to
simultaneously measure diffusion coefficients of both peptide and surfactant.
These two inputs, as a function of SDS concentration, are then fit to a simple
two species model that neglects hydrodynamic interactions between complexes.
From this we obtain the fraction of free SDS, and the hydrodynamic size of
complexes in a GAD-2--SDS system as a function of SDS concentration. These
results are compared to those for smaller dipeptides and for peptide-free
solutions. At low SDS concentrations ([SDS] $\leq$ 25 mM), the results
self-consistently point to a GAD-2--SDS complex of fixed hydrodynamic size R
=(5.5 $\pm$ 0.3) nm. At intermediate SDS concentrations (25 mM $<$ [SDS] $<$ 60
mM), the apparent size of a GAD-2--SDS complex shows almost a factor of two
increase without a significant change in surfactant-to-peptide ratio within a
complex, most likely implying an increase in the number of peptides in a
complex. For peptide-free solutions, the self-diffusion coefficients of SDS
with and without buffer are significantly different at low SDS concentrations
but merge above [SDS]=60 mM. We find that in order to obtain unambiguous
information about the hydrodynamic size of a peptide-surfactant complex from
diffusion measurements, experiments must be carried out at or below [SDS] = 25
mM.


  Protein interaction networks are a promising type of data for studying
complex biological systems. However, despite the rich information embedded in
these networks, they face important data quality challenges of noise and
incompleteness that adversely affect the results obtained from their analysis.
Here, we explore the use of the concept of common neighborhood similarity
(CNS), which is a form of local structure in networks, to address these issues.
Although several CNS measures have been proposed in the literature, an
understanding of their relative efficacies for the analysis of interaction
networks has been lacking. We follow the framework of graph transformation to
convert the given interaction network into a transformed network corresponding
to a variety of CNS measures evaluated. The effectiveness of each measure is
then estimated by comparing the quality of protein function predictions
obtained from its corresponding transformed network with those from the
original network. Using a large set of S. cerevisiae interactions, and a set of
136 GO terms, we find that several of the transformed networks produce more
accurate predictions than those obtained from the original network. In
particular, the $HC.cont$ measure proposed here performs particularly well for
this task. Further investigation reveals that the two major factors
contributing to this improvement are the abilities of CNS measures, especially
$HC.cont$, to prune out noisy edges and introduce new links between
functionally related proteins.


  The public goods game is one of the most famous models for studying the
evolution of cooperation in sizable groups. The multiplication factor in this
game can characterize the investment return from the public good, which may be
variable depending on the interactive environment in realistic situations.
Instead of using the same universal value, here we consider that the
multiplication factor in each group is updated based on the differences between
the local and global interactive environments in the spatial public goods game,
but meanwhile limited to within a certain range. We find that the adaptive and
bounded investment returns can significantly promote cooperation. In
particular, full cooperation can be achieved for high feedback strength when
appropriate limitation is set for the investment return. Also, we show that the
fraction of cooperators in the whole population can become larger if the lower
and upper limits of the multiplication factor are increased. Furthermore, in
comparison to the traditionally spatial public goods game where the
multiplication factor in each group is identical and fixed, we find that
cooperation can be better promoted if the multiplication factor is constrained
to adjust between one and the group size in our model. Our results highlight
the importance of the locally adaptive and bounded investment returns for the
emergence and dominance of cooperative behavior in structured populations.


  In the study of the evolution of cooperation, resource limitations are
usually assumed just to provide a finite population size. Recently, however,
agent-based models have pointed out that resource limitation may modify the
original structure of the interactions and allow for the survival of
unconditional cooperators in well-mixed populations. Here, we present
analytical simplified versions of two types of agent-based models recently
published: one in which the limiting resource constrains the ability of
reproduction of individuals but not their survival, and a second one where the
limiting resource is necessary for both reproduction and survival. One finds
that the analytical models display, with a few differences, the same
qualitative behavior of the more complex agent-based models. In addition, the
analytical models allow us to expand the study and identify the dimensionless
parameters governing the final fate of the system, such as coexistence of
cooperators and defectors, or dominance of defectors or of cooperators. We
provide a detailed analysis of the occurring phase transitions as these
parameters are varied.


  Riboswitches, structured elements in the untranslated regions of messenger
RNAs, regulate gene expression by binding specific metabolites. We introduce a
kinetic network model that describes the functions of riboswitches at the
systems level. Using experimental data for flavin mono nucleotide riboswitch as
a guide we show that efficient function, implying a large dynamic range without
compromising the requirement to suppress transcription, is determined by a
balance between the transcription speed, the folding and unfolding rates of the
aptamer, and the binding rates of the metabolite. We also investigated the
effect of negative feedback accounting for binding to metabolites, which are
themselves the products of genes that are being regulated. For a range of
transcription rates negative feedback suppresses gene expression by nearly 10
fold. Negative feedback speeds the gene expression response time, and
suppresses the change of steady state protein concentration by half relative to
that without feedback, when there is a modest spike in DNA concentration. A
dynamic phase diagram expressed in terms of transcription speed, folding rates,
and metabolite binding rates predicts different scenarios in
riboswitch-mediated transcription regulation.


  The spatial rock-paper-scissors ecosystem, where three species interact
cyclically, is a model example of how spatial structure can maintain
biodiversity. We here consider such a system for a broad range of interaction
rates. When one species grows very slowly, this species and its prey dominate
the system by self-organizing into a labyrinthine configuration in which the
third species propagates. The cluster size distributions of the two dominating
species have heavy tails and the configuration is stabilized through a complex,
spatial feedback loop. We introduce a new statistical measure that quantifies
the amount of clustering in the spatial system by comparison with its mean
field approximation. Hereby, we are able to quantitatively explain how the
labyrinthine configuration slows down the dynamics and stabilizes the system.


  The emergence and promotion of cooperation are two of the main issues in
evolutionary game theory, as cooperation is amenable to exploitation by
defectors, which take advantage of cooperative individuals at no cost, dooming
them to extinction. It has been recently shown that the existence of purely
destructive agents (termed jokers) acting on the common enterprises (public
goods games) can induce stable limit cycles among cooperation, defection, and
destruction when infinite populations are considered. These cycles allow for
time lapses in which cooperators represent a relevant fraction of the
population, providing a mechanism for the emergence of cooperative states in
nature and human societies. Here we study analytically and through agent-based
simulations the dynamics generated by jokers in finite populations for several
selection rules. Cycles appear in all cases studied, thus showing that the
joker dynamics generically yields a robust cyclic behavior not restricted to
infinite populations. We also compute the average time in which the population
consists mostly of just one strategy and compare the results with numerical
simulations.


  Although mechanical properties of DNA are well characterized at the kilo
base-pair range, a number of recent experiments have suggested that DNA is more
flexible at shorter length scales, which correspond to the regime that is
crucial for cellular processes such as DNA packaging and gene regulation. Here,
we perform a systematic study of the effective elastic properties of DNA at
different length scales by probing the conformation and fluctuations of DNA
from single base-pair level up to four helical turns, using trajectories from
atomistic simulation. We find evidence that supports cooperative softening of
the stretch modulus and identify the essential modes that give rise to this
effect. The bend correlation exhibits modulations that reflect the helical
periodicity, while it yields a reasonable value for the effective persistence
length, and the twist modulus undergoes a smooth crossover---from a relatively
smaller value at the single base-pair level to the bulk value---over half a
DNA-turn.


  Generalized linear mixed-effects models in the context of genome-wide
association studies (GWAS) represent a formidable computational challenge: the
solution of millions of correlated generalized least-squares problems, and the
processing of terabytes of data. We present high performance in-core and
out-of-core shared-memory algorithms for GWAS: By taking advantage of
domain-specific knowledge, exploiting multi-core parallelism, and handling data
efficiently, our algorithms attain unequalled performance. When compared to
GenABEL, one of the most widely used libraries for GWAS, on a 12-core processor
we obtain 50-fold speedups. As a consequence, our routines enable genome
studies of unprecedented size.


  A method is presented for measuring anisotropic THz response for small
crystals, Crystal Anisotropy Terahertz Microscopy (CATM). Sucrose CATM
measurements find the expected anisotropic phonon resonances. CATM measurements
of protein crystals find the expected broadband water absorption is suppressed
and strong orientation and hydration dependent resonant features.


  We introduce Modulated Orientation Sensitive Terahertz Spectroscopy (MOSTS)
that can both increase contrast between adjacent modes and remove the broad
band librational background to allow mode identification for large macro
molecules such as proteins.


  The detection of molecular signatures of selection is one of the major
concerns of modern population genetics. A widely used strategy in this context
is to compare samples from several populations, and to look for genomic regions
with outstanding genetic differentiation between these populations. Genetic
differentiation is generally based on allele frequency differences between
populations, which are measured by Fst or related statistics. Here we introduce
a new statistic, denoted hapFLK, which focuses instead on the differences of
haplotype frequencies between populations. In contrast to most existing
statistics, hapFLK accounts for the hierarchical structure of the sampled
populations. Using computer simulations, we show that each of these two
features - the use of haplotype information and of the hierarchical structure
of populations - significantly improves the detection power of selected loci,
and that combining them in the hapFLK statistic provides even greater power. We
also show that hapFLK is robust with respect to bottlenecks and migration and
improves over existing approaches in many situations. Finally, we apply hapFLK
to a set of six sheep breeds from Northern Europe, and identify seven regions
under selection, which include already reported regions but also several new
ones. We propose a method to help identifying the population(s) under selection
in a detected region, which reveals that in many of these regions selection
most likely occurred in more than one population. Furthermore, several of the
detected regions correspond to incomplete sweeps, where the favourable
haplotype is only at intermediate frequency in the population(s) under
selection.


  Global dynamical behaviors of the competitive Lotka-Volterra system even in
3-dimension are not fully understood. The Lyapunov function can provide us such
knowledge once it is constructed. In this paper, we construct explicitly the
Lyapunov function in three examples of the competitive Lotka-Volterra system
for the whole state space: (1) the general 2-dimensional case; (2) a
3-dimensional model; (3) the model of May-Leonard. The dynamics of these
examples include bistable case and cyclical behavior. The first two examples
are the generalized gradient system defined in the Appendixes, while the model
of May-Leonard is not. Our method is helpful to understand the limit cycle
problems in general 3-dimensional case.


  In many scientific and engineering applications, one has to solve not one but
a sequence of instances of the same problem. Often times, the problems in the
sequence are linked in a way that allows intermediate results to be reused. A
characteristic example for this class of applications is given by the
Genome-Wide Association Studies (GWAS), a widely spread tool in computational
biology. GWAS entails the solution of up to trillions ($10^{12}$) of correlated
generalized least-squares problems, posing a daunting challenge: the
performance of petaflops ($10^{15}$ floating-point operations) over terabytes
of data.
  In this paper, we design an algorithm for performing GWAS on multi-core
architectures. This is accomplished in three steps. First, we show how to
exploit the relation among successive problems, thus reducing the overall
computational complexity. Then, through an analysis of the required data
transfers, we identify how to eliminate any overhead due to input/output
operations. Finally, we study how to decompose computation into tasks to be
distributed among the available cores, to attain high performance and
scalability. With our algorithm, a GWAS that currently requires the use of a
supercomputer may now be performed in matter of hours on a single multi-core
node.
  The discussion centers around the methodology to develop the algorithm rather
than the specific application. We believe the paper contributes valuable
guidelines of general applicability for computational scientists on how to
develop and optimize numerical algorithms.


  A new statistical mechanics formulation of characterizing the structural
fluctuation of protein correlated with that of water is presented based on the
generalized Langevin equation and the 3D-RISM/RISM theory of molecular liquids.
The displacement vector of atom positions and their conjugated momentum, are
chosen for the dynamic variables for protein, while the density fields of atoms
and their momentum fields are chosen for water. Projection of other degrees of
freedom onto those dynamic variables using the standard projection operator
method produces essentially two equations which describe the time evolution of
fluctuation concerning the density field of solvent and the conformation of
protein around an equilibrium state, which are coupled with each other. The
equation concerning the protein dynamics is formally akin to that of the
coupled Langevin oscillators, and is a generalization of the latter, to atomic
level. The most intriguing feature of the new equation is that it contains the
variance-covariance matrix as the "Hessian" term describing the "force"
restoring an equilibrium conformation, which is the second moment of the
fluctuation of atom positions. The "Hessian" matrix is naturally identified as
the second derivative of the free energy surface around the equilibrium. A
method to evaluate the Hessian matrix based on the 3D-RISM/RISM theory is
proposed. Proposed also is an application of the present formulation to the
molecular recognition, in which the conformational fluctuation of protein
around its native state becomes an important factor as exemplified by so called
"induced fitting".


  Some mathematical models for the estimation of the effects of Cry1Ab and
Cry1F Bt-maize exposure in the biodiversity are discussed. Novel results about
these models are obtained and described in the note. The exact formula for the
proportion of the population that suffers mortality exposed to Cry1Ab pollen,
underlining its dependence on the margin from the Bt crop edge, is derived. In
addition, regarding Cry1F pollen effects, it is proposed a procedure, using a
probabilistic and statistical approach, that computes the width of the non
Bt-stripes used as mitigation measures. Finally, it has been derived a lower
bound, using probabilistic consideration, on the species sensitivity of
Lepidoptera.


  When small RNAs are loaded onto Argonaute proteins they can form the
RNA-induced silencing complexes (RISCs), which mediate RNA interference.
RISC-formation is dependent on a shared pool of Argonaute proteins and RISC
loading factors, and is thus susceptible to competition among small RNAs for
loading. We present a mathematical model that aims to understand how small RNA
competition for the PTR resources affects target gene repression. We discuss
that small RNA activity is limited by RISC-formation, RISC-degradation and the
availability of Argonautes. Together, these observations explain a number of
PTR saturation effects encountered experimentally. We show that different
competition conditions for RISC-loading result in different signatures of PTR
activity determined also by the amount of RISC-recycling taking place. In
particular, we find that the small RNAs less efficient at RISC-formation, using
fewer resources of the PTR pathway, can perform in the low RISC-recycling range
equally well as their more effective counterparts. Additionally, we predict a
novel signature of PTR in target expression levels. Under conditions of low
RISC-loading efficiency and high RISC-recycling, the variation in target levels
increases linearly with the target transcription rate. Furthermore, we show
that RISC-recycling determines the effect that Argonaute scarcity conditions
have on target expression variation. Our observations taken together offer a
framework of predictions which can be used in order to infer from experimental
data the particular characteristics of underlying PTR activity.


  We identify three dimeric proteins whose mechanostability is anisotropic and
should exceed 1 nN along some directions. They come with distinct mechanical
clamps: shear-based, involving a cystine slipknot, and due to dragging of a
cystine plug through a cystine ring. The latter two mechanisms are topological
in nature and the cystine plug mechanism has not yet been discussed but it
turns out to provide the largest resistance to stretching. Its possible
applications in elastomers are discussed.


  We construct a two dimensional Cellular Automata based model for the
description of pedestrian dynamics. Wide range of complicated pattern formation
phenomena in pedestrian dynamics are described in the model, e.g. lane
formation, jams in a counterflow and egress behavior. Mean-field solution of
the densely populated case and numerical solution of the sparsely populated
case are provided. This model has the potential to describe more flow
phenomena.


  Species diversity in ecosystems is often accompanied by the self-organisation
of the population into fascinating spatio-temporal patterns. Here, we consider
a two-dimensional three-species population model and study the spiralling
patterns arising from the combined effects of generic cyclic dominance,
mutation, pair-exchange and hopping of the individuals. The dynamics is
characterised by nonlinear mobility and a Hopf bifurcation around which the
system's phase diagram is inferred from the underlying complex Ginzburg-Landau
equation derived using a perturbative multiscale expansion. While the dynamics
is generally characterised by spiralling patterns, we show that spiral waves
are stable in only one of the four phases. Furthermore, we characterise a phase
where nonlinearity leads to the annihilation of spirals and to the spatially
uniform dominance of each species in turn. Away from the Hopf bifurcation, when
the coexistence fixed point is unstable, the spiralling patterns are also
affected by nonlinear diffusion.


  The cytoskeleton of eukaryotic cells provides mechanical support and governs
intracellular transport. These functions rely on the complex mechanical
properties of networks of semiflexible protein filaments. Recent theoretical
interest has focused on mesoscopic properties of such networks and especially
on the effect of local, non-affine bending deformations on mechanics. Here, we
study the impact of local network deformations on the scale-dependent mobility
of probe particles in entangled networks of semiflexible actin filaments by
high-bandwidth microrheology. We find that micron-sized particles in these
networks experience two opposing non-continuum elastic effects: entropic
depletion reduces the effective network rigidity, while local non-affine
deformations of the network substantially enhance the rigidity at low
frequencies. We show that a simple model of lateral bending of filaments
embedded in a viscoelastic background leads to a scaling regime for the
apparent elastic modulus G'(\omega) \sim \omega^{9/16}, closely matching the
experiments. These results provide quantitative evidence for how different a
semiflexible polymer network can feel for small objects, and they demonstrate
how non-affine bending deformations can be dominant for the mobility of
vesicles and organelles in the cell.


  In spatial games players typically alter their strategy by imitating the most
successful or one randomly selected neighbor. Since a single neighbor is taken
as reference, the information stemming from other neighbors is neglected, which
begets the consideration of alternative, possibly more realistic approaches.
Here we show that strategy changes inspired not only by the performance of
individual neighbors but rather by entire neighborhoods introduce a
qualitatively different evolutionary dynamics that is able to support the
stable existence of very small cooperative clusters. This leads to phase
diagrams that differ significantly from those obtained by means of pairwise
strategy updating. In particular, the survivability of cooperators is possible
even by high temptations to defect and over a much wider uncertainty range. We
support the simulation results by means of pair approximations and analysis of
spatial patterns, which jointly highlight the importance of local information
for the resolution of social dilemmas.


  In the classic view introduced by R. A. Fisher, a quantitative trait is
encoded by many loci with small, additive effects. Recent advances in QTL
mapping have begun to elucidate the genetic architectures underlying vast
numbers of phenotypes across diverse taxa, producing observations that
sometimes contrast with Fisher's blueprint. Despite these considerable
empirical efforts to map the genetic determinants of traits, it remains poorly
understood how the genetic architecture of a trait should evolve, or how it
depends on the selection pressures on the trait. Here we develop a simple,
population-genetic model for the evolution of genetic architectures. Our model
predicts that traits under moderate selection should be encoded by many loci
with highly variable effects, whereas traits under either weak or strong
selection should be encoded by relatively few loci. We compare these
theoretical predictions to qualitative trends in the genetics of human traits,
and to systematic data on the genetics of gene expression levels in yeast. Our
analysis provides an evolutionary explanation for broad empirical patterns in
the genetic basis of traits, and it introduces a single framework that unifies
the diversity of observed genetic architectures, ranging from Mendelian to
Fisherian.


  Population genomic studies have shown that genetic draft and background
selection can profoundly affect the genome-wide patterns of molecular
variation. We performed forward simulations under realistic gene-structure and
selection scenarios to investigate whether such linkage effects impinge on the
ability of the McDonald-Kreitman (MK) test to infer the rate of positive
selection (\alpha) from polymorphism and divergence data. We find that in the
presence of slightly deleterious mutations, MK estimates of \alpha\ severely
underestimate the true rate of adaptation even if all polymorphisms with
population frequencies under 50% are excluded. Furthermore, already under
intermediate rates of adaptation, genetic draft substantially distorts the site
frequency spectra at neutral and functional sites from the expectations under
mutation-selection-drift balance. MK-type approaches that first infer
demography from synonymous sites and then use the inferred demography to
correct the estimation of \alpha\ obtain almost the correct \alpha\ in our
simulations. However, these approaches typically infer a severe past population
expansion although there was no such expansion in the simulations, casting
doubt on the accuracy of methods that infer demography from synonymous
polymorphism data. We suggest a simple asymptotic extension of the MK test that
should yield accurate estimates of \alpha\ even in the presence of linkage
effects.


  We present a generic analytical scheme for the quantification of fluctuations
due to bifunctionality-induced signal transduction within the members of
bacterial two-component system. The proposed model takes into account
post-translational modifications in terms of elementary phosphotransfer
kinetics. Sources of fluctuations due to autophosphorylation, kinase and
phosphatase activity of the sensor kinase have been considered in the model via
Langevin equations, which are then solved within the framework of linear noise
approximation. The resultant analytical expression of phosphorylated response
regulators are then used to quantify the noise profile of biologically
motivated single and branched pathways. Enhancement and reduction of noise in
terms of extra phosphate outflux and influx, respectively, have been analyzed
for the branched system. Furthermore, role of fluctuations of the network
output in the regulation of a promoter with random activation/deactivation
dynamics has been analyzed.


  Magnetic hyperthermia is currently an EU-approved clinical therapy against
tumor cells that uses magnetic nanoparticles under a time varying magnetic
field (TVMF). The same basic principle seems promising against trypanosomatids
causing Chagas disease and sleeping sickness, since therapeutic drugs available
display severe side effects and drug-resistant strains. However, no
applications of this strategy against protozoan-induced diseases have been
reported so far. In the present study, Crithidia fasciculata, a widely used
model for therapeutic strategies against pathogenic trypanosomatids, was
targeted with Fe_{3}O_{4} magnetic nanoparticles (MNPs) in order to remotely
provoke cell death using TVMFs. The MNPs with average sizes of d approx. 30 nm
were synthesized using a precipitation of FeSO_{4}4 in basic medium. The MNPs
were added to Crithidia fasciculata choanomastigotes in exponential phase and
incubated overnight. The amount of uploaded MNPs per cell was determined by
magnetic measurements. Cell viability using the MTT colorimetric assay and flow
cytometry showed that the MNPs were incorporated by the cells with no
noticeable cell-toxicity effects. When a TVMF (f = 249 kHz, H = 13 kA/m) was
applied to MNP-bearing cells, massive cell death was induced via a
non-apoptotic mechanism. No effects were observed by applying a TVMF on control
(without loaded MNPs) cells. No macroscopic rise in temperature was observed in
the extracellular medium during the experiments. Scanning Electron Microscopy
showed morphological changes after TVMF experiments. These data indicate (as a
proof of principle) that intracellular hyperthermia is a suitable technology to
induce the specific death of protozoan parasites bearing MNPs. These findings
expand the possibilities for new therapeutic strategies that combat parasitic
infections.


  Neuronal growth cones are the most sensitive amongst eukaryotic cells in
responding to directional chemical cues. Although a dynamic microtubule
cytoskeleton has been shown to be essential for growth cone turning, the
precise nature of coupling of the spatial cue with microtubule polarization is
less understood. Here we present a computational model of microtubule
polarization in a turning neuronal growth cone (GC). We explore the limits of
directional cues in modifying the spatial polarization of microtubules by
testing the role of microtubule dynamics, gradients of regulators and
retrograde forces along filopodia. We analyze the steady state and transition
behavior of microtubules on being presented with a directional stimulus. The
model makes novel predictions about the minimal angular spread of the chemical
signal at the growth cone and the fastest polarization times. A regulatory
reaction-diffusion network based on the cyclic
phosphorylation-dephosphorylation of a regulator predicts that the receptor
signal magnitude can generate the maximal polarization of microtubules and not
feedback loops or amplifications in the network. Using both the
phenomenological and network models we have demonstrated some of the physical
limits within which the MT polarization system works in turning neuron.


  Long-range migrations and the resulting admixtures between populations have
been important forces shaping human genetic diversity. Most existing methods
for detecting and reconstructing historical admixture events are based on
allele frequency divergences or patterns of ancestry segments in chromosomes of
admixed individuals. An emerging new approach harnesses the exponential decay
of admixture-induced linkage disequilibrium (LD) as a function of genetic
distance. Here, we comprehensively develop LD-based inference into a versatile
tool for investigating admixture. We present a new weighted LD statistic that
can be used to infer mixture proportions as well as dates with fewer
constraints on reference populations than previous methods. We define an
LD-based three-population test for admixture and identify scenarios in which it
can detect admixture events that previous formal tests cannot. We further show
that we can uncover phylogenetic relationships among populations by comparing
weighted LD curves obtained using a suite of references. Finally, we describe
several improvements to the computation and fitting of weighted LD curves that
greatly increase the robustness and speed of the calculations. We implement all
of these advances in a software package, ALDER, which we validate in
simulations and apply to test for admixture among all populations from the
Human Genome Diversity Project (HGDP), highlighting insights into the admixture
history of Central African Pygmies, Sardinians, and Japanese.


  MicroRNAs (miRNAs) critically modulate stem cell properties like
pluripotency, but the fundamental mechanism remains largely unknown. This study
systematically analyzes multiple-omics data and builds a systems physical
network including genome-wide interactions between miRNAs and their targets to
reveal the systems mechanism of miRNA functions in mouse pluripotent stem
cells. Globally, miRNAs directly repress the pluripotent core factors during
differentiation state. Surprisingly, during pluripotent state, the top
important miRNAs do not directly regulate the pluripotent core factors as
thought, but they only directly target the pluripotent signal pathways and
directly repress developmental processes. Furthermore, at pluripotent state
miRNAs predominately repress DNA methyltransferases, the core enzymes for DNA
methylation. The decreasing methylation repressed by miRNAs in turn activates
the top miRNAs and pluripotent core factors, creating an active circuit system
to modulate pluripotency. MiRNAs vary their functions with different stem cell
states. While miRNAs directly repress pluripotent core factors to facilitate
the differentiation during cell differentiation, they also help stem cells to
maintain pluripotency by activating pluripotent cores through directly
repressing DNA methylation systems and primarily inhibiting development.


  The simplest possible beginning of abiogenesis has been a riddle from the
last century, which is most successfully solved by the Lipid World hypothesis.
However, origin of the next stages of evolution starting form lipids is still
in dark. We propose a 'Lipid-RNA World Scenario' based on the assumption that
modern stable lipid-RNA interactions are molecular fossils of an ancient stage
of evolution when RNA World originated from Lipid World. In accordance to the
faint young sun conditions, we present an 'ice-covered hydrothermal vent' model
of Hadean Ocean. Our hypothetical model suggests that faint young sun condition
probably provided susceptible physical conditions for an evolutionary route
from Lipid-World to Protein-RNA World, through an intermediate Lipid-RNA World.
Ancient ribozymes were 'protected' by lipids assuring their survival in
prebiotic ocean. The origin of natural selection ensures transition of
Lipid-RNA World to Protein-RNA World after the origin of ribosome. Assuming the
modern peptidyltransferase as the proto-ribosome structure, we have presented a
hypothetical translation mechanism: proto-ribosome randomly polymerized amino
acids being attached to the inner layer of a lipid-vesicle, using only physical
energies available from our Hadean Ocean model. In accordance to the strategy
of chemical evolution, we also have described the possible evolutionary
behavior of this proto-ribosome, which explains the contemporary
three-dimensional structure of 50S subunit and supports the predictions
regarding the ancient regions of it. It also explains the origin of
membrane-free 'minimal ribosome' in the time of LUCA.