  The Caribbean basin is home to some of the most complex interactions in
recent history among previously diverged human populations. Here, by making use
of genome-wide SNP array data, we characterize ancestral components of
Caribbean populations on a sub-continental level and unveil fine-scale patterns
of population structure distinguishing insular from mainland Caribbean
populations as well as from other Hispanic/Latino groups. We provide genetic
evidence for an inland South American origin of the Native American component
in island populations and for extensive pre-Columbian gene flow across the
Caribbean basin. The Caribbean-derived European component shows significant
differentiation from parental Iberian populations, presumably as a result of
founder effects during the colonization of the New World. Based on demographic
models, we reconstruct the complex population history of the Caribbean since
the onset of continental admixture. We find that insular populations are best
modeled as mixtures absorbing two pulses of African migrants, coinciding with
early and maximum activity stages of the transatlantic slave trade. These two
pulses appear to have originated in different regions within West Africa,
imprinting two distinguishable signatures in present day Afro-Caribbean genomes
and shedding light on the genetic impact of the dynamics occurring during the
slave trade in the Caribbean.


  What is a population? This review considers how a population may be defined
in terms of understanding the structure of the underlying genetics of the
individuals involved. The main approach is to consider statistically
identifiable groups of randomly mating individuals, which is well defined in
theory for any type of (sexual) organism. We discuss generative models using
drift, admixture and spatial structure, and the ancestral recombination graph.
These are contrasted with statistical models for inference, principle component
analysis and other `non-parametric' methods. The relationships between these
approaches are explored with both simulated and real-data examples. The
state-of-the-art practical software tools are discussed and contrasted. We
conclude that populations are a useful theoretical construct that can be well
defined in theory and often approximately exist in practice.


  Sequence alignment data is often ordered by coordinate (id of the reference
sequence plus position on the sequence where the fragment was mapped) when
stored in BAM files, as this simplifies the extraction of variants between the
mapped data and the reference or of variants within the mapped data. In this
order paired reads are usually separated in the file, which complicates some
other applications like duplicate marking or conversion to the FastQ format
which require to access the full information of the pairs. In this paper we
introduce biobambam, an API for efficient BAM file reading supporting the
efficient collation of alignments by read name without performing a complete
resorting of the input file and some tools based on this API performing tasks
like marking duplicate reads and conversion to the FastQ format. In comparison
with previous approaches to problems involving the collation of alignments by
read name like the BAM to FastQ or duplication marking utilities in the Picard
suite the approach of biobambam can often perform an equivalent task more
efficiently in terms of the required main memory and run-time.


  In sexual populations, selection operates neither on the whole genome, which
is repeatedly taken apart and reassembled by recombination, nor on individual
alleles that are tightly linked to the chromosomal neighborhood. The resulting
interference between linked alleles reduces the efficiency of selection and
distorts patterns of genetic diversity. Inference of evolutionary history from
diversity shaped by linked selection requires an understanding of these
patterns. Here, we present a simple but powerful scaling analysis identifying
the unit of selection as the genomic "linkage block" with a characteristic
length determined in a self-consistent manner by the condition that the rate of
recombination within the block is comparable to the fitness differences between
different alleles of the block. We find that an asexual model with the strength
of selection tuned to that of the linkage block provides an excellent
description of genetic diversity and the site frequency spectra when compared
to computer simulations. This linkage block approximation is accurate for the
entire spectrum of strength of selection and is particularly powerful in
scenarios with many weakly selected loci. The latter limit allows us to
characterize coalescence, genetic diversity, and the speed of adaptation in the
infinitesimal model of quantitative genetics.


  Pervasive natural selection can strongly influence observed patterns of
genetic variation, but these effects remain poorly understood when multiple
selected variants segregate in nearby regions of the genome. Classical
population genetics fails to account for interference between linked mutations,
which grows increasingly severe as the density of selected polymorphisms
increases. Here, we describe a simple limit that emerges when interference is
common, in which the fitness effects of individual mutations play a relatively
minor role. Instead, molecular evolution is determined by the variance in
fitness within the population, defined over an effectively asexual segment of
the genome (a ``linkage block''). We exploit this insensitivity in a new
``coarse-grained'' coalescent framework, which approximates the effects of many
weakly selected mutations with a smaller number of strongly selected mutations
with the same variance in fitness. This approximation generates accurate and
efficient predictions for the genetic diversity that cannot be summarized by a
simple reduction in effective population size. However, these results suggest a
fundamental limit on our ability to resolve individual selection pressures from
contemporary sequence data alone, since a wide range of parameters yield nearly
identical patterns of sequence variability.


  We consider birth-and-death stochastic evolution of genotypes with different
lengths. The genotypes might mutate that provides a stochastic changing of
lengthes by a free diffusion law. The birth and death rates are length
dependent which corresponds to a selection effect. We study an asymptotic
behavior of a density for an infinite collection of genotypes. The cases of
space homogeneous and space heterogeneous densities are considered.


  Hierarchical autocorrelation in the error term of linear models arises when
sampling units are related to each other according to a tree. The residual
covariance is parametrized using the tree-distance between sampling units. When
observations are modeled using an Ornstein-Uhlenbeck (OU) process along the
tree, the autocorrelation between two tips decreases exponentially with their
tree distance. These models are most often applied in evolutionary biology,
when tips represent biological species and the OU process parameters represent
the strength and direction of natural selection. For these models, we show that
the mean is not microergodic: no estimator can ever be consistent for this
parameter and provide a lower bound for the variance of its MLE. For covariance
parameters, we give a general sufficient condition ensuring microergodicity.
This condition suggests that some parameters may not be estimated at the same
rate as others. We show that, indeed, maximum likelihood estimators of the
autocorrelation parameter converge at a slower rate than that of generally
microergodic parameters. We showed this theoretically in a symmetric tree
asymptotic framework and through simulations on a large real tree comprising
4507 mammal species.


  Finding optimal evolutionary trees from sequence data is typically an
intractable problem, and there is usually no way of knowing how close to
optimal the best tree from some search truly is. The problem would seem to be
particularly acute when we have many taxa and when that data has high levels of
homoplasy, in which the individual characters require many changes to fit on
the best tree. However, a recent mathematical result has provided a precise
tool to generate a short number of high-homoplasy characters for any given
tree, so that this tree is provably the optimal tree under the maximum
parsimony criterion. This provides, for the first time, a rigorous way to test
tree search algorithms on homoplasy-rich data, where we know in advance what
the `best' tree is. In this short note we consider just one search program
(TNT) but show that it is able to locate the globally optimal tree correctly
for 32,768 taxa, even though the characters in the dataset requires, on
average, 1148 state-changes each to fit on this tree, and the number of
characters is only 57.


  To determine the 3D conformation of proteins is a necessity to understand
their functions or interactions with other molecules. It is commonly admitted
that, when proteins fold from their primary linear structures to their final 3D
conformations, they tend to choose the ones that minimize their free energy. To
find the 3D conformation of a protein knowing its amino acid sequence,
bioinformaticians use various models of different resolutions and artificial
intelligence tools, as the protein folding prediction problem is a NP complete
one. More precisely, to determine the backbone structure of the protein using
the low resolution models (2D HP square and 3D HP cubic), by finding the
conformation that minimize free energy, is intractable exactly. Both the proof
of NP-completeness and the 2D prediction consider that acceptable conformations
have to satisfy a self-avoiding walk (SAW) requirement, as two different amino
acids cannot occupy a same position in the lattice. It is shown in this
document that the SAW requirement considered when proving NP-completeness is
different from the SAW requirement used in various prediction programs, and
that they are different from the real biological requirement. Indeed, the proof
of NP completeness and the predictions in silico consider conformations that
are not possible in practice. Consequences of this fact are investigated in
this research work.


  We study kinetic model of Nuclear Receptor Binding to Promoter Regions. This
model is written as a system of ordinary differential equations. Model
reduction techniques have been used to simplify chemical kinetics.In this case
study, the technique of Pseudo-first order approximation is applied to simplify
the reaction rates. CellDesigner has been used to draw the structures of
chemical reactions of Nuclear Receptor Binding to Promoter Regions. After model
reduction, the general analytical solution for reduced model is given and the
number of species and reactions are reduced from 9 species and 6 reactions to 6
species and 5 reactions.


  Self-avoiding walks (SAW) are the source of very difficult problems in
probabilities and enumerative combinatorics. They are also of great interest as
they are, for instance, the basis of protein structure prediction in
bioinformatics. Authors of this article have previously shown that, depending
on the prediction algorithm, the sets of obtained conformations differ: all the
self-avoiding walks can be reached using stretching-based algorithms whereas
only the folded SAWs can be attained with methods that iteratively fold the
straight line. A first study of (un)folded self-avoiding walks is presented in
this article. The contribution is majorly a survey of what is currently known
about these sets. In particular we provide clear definitions of various subsets
of self-avoiding walks related to pivot moves (folded or unfoldable SAWs, etc.)
and the first results we have obtained, theoretically or computationally, on
these sets. A list of open questions is provided too, and the consequences on
the protein structure prediction problem is finally investigated.


  Transcriptome assembly from RNA-Seq reads is an active area of bioinformatics
research. The ever-declining cost and the increasing depth of RNA-Seq have
provided unprecedented opportunities to better identify expressed transcripts.
However, the nonlinear transcript structures and the ultra-high throughput of
RNA-Seq reads pose significant algorithmic and computational challenges to the
existing transcriptome assembly approaches, either reference-guided or de novo.
While reference-guided approaches offer good sensitivity, they rely on
alignment results of the splice-aware aligners and are thus unsuitable for
species with incomplete reference genomes. In contrast, de novo approaches do
not depend on the reference genome but face a computational daunting task
derived from the complexity of the graph built for the whole transcriptome. In
response to these challenges, we present a hybrid approach to exploit an
incomplete reference genome without relying on splice-aware aligners. We have
designed a split-and-align procedure to efficiently localize the reads to
individual genomic loci, which is followed by an accurate de novo assembly to
assemble reads falling into each locus. Using extensive simulation data, we
demonstrate a high accuracy and precision in transcriptome reconstruction by
comparing to selected transcriptome assembly tools. Our method is implemented
in assemblySAM, a GUI software freely available at
http://sammate.sourceforge.net.


  Cheap high-throughput DNA sequencing may soon become routine not only for
human genomes but also for practically anything requiring the identification of
living organisms from their DNA: tracking of infectious agents, control of food
products, bioreactors, or environmental samples.
  We propose a novel general approach to the analysis of sequencing data in
which the reference genome does not have to be specified. Using a distributed
architecture we are able to query a remote server for hints about what the
reference might be, transferring a relatively small amount of data, and the
hints can be used for more computationally-demanding work.
  Our system consists of a server with known reference DNA indexed, and a
client with raw sequencing reads. The client sends a sample of unidentified
reads, and in return receives a list of matching references known to the
server. Sequences for the references can be retrieved and used for exhaustive
computation on the reads, such as alignment.
  To demonstrate this approach we have implemented a web server, indexing tens
of thousands of publicly available genomes and genomic regions from various
organisms and returning lists of matching hits from query sequencing reads. We
have also implemented two clients, one of them running in a web browser, in
order to demonstrate that gigabytes of raw sequencing reads of unknown origin
could be identified without the need to transfer a very large volume of data,
and on modestly powered computing devices.
  A web access is available at http://tapir.cbs.dtu.dk. The source code for a
python command-line client, a server, and supplementary data is available at
http://bit.ly/1aURxkc.


  We investigate the effect of spatial range expansions on the evolution of
fitness when beneficial and deleterious mutations co-segregate. We perform
individual-based simulations of a uniform linear habitat and complement them
with analytical approximations for the evolution of mean fitness at the edge of
the expansion. We find that deleterious mutations accumulate steadily on the
wave front during range expansions, thus creating an expansion load. Reduced
fitness due to the expansion load is not restricted to the wave front but
occurs over a large proportion of newly colonized habitats. The expansion load
can persist and represent a major fraction of the total mutation load thousands
of generations after the expansion. Our results extend qualitatively and
quantitatively to two-dimensional expansions. The phenomenon of expansion load
may explain growing evidence that populations that have recently expanded,
including humans, show an excess of deleterious mutations. To test the
predictions of our model, we analyze patterns of neutral and non-neutral
genetic diversity in humans and find an excellent fit between theory and data.


  A coarse-grained computational model is used to investigate the effect of a
fluctuating fluid membrane on the dynamics of patchy-particle assembly into
virus capsid-like cores. Results from simulations for a broad range of
parameters are presented, showing the effect of varying interaction strength,
membrane stiffness and membrane viscosity. Furthermore, the effect of
hydrodynamic interactions is investigated. Attraction to a membrane may promote
assembly, including for sub-unit interaction strengths for which it does not
occur in the bulk, and may also decrease single-core assembly time. The
membrane budding rate is strongly increased by hydrodynamic interactions. The
membrane deformation rate is important in determining the finite-time yield.
Higher rates may decrease the entropic penalty for assembly and help guide
sub-units towards each other but may also block partial cores from being
completed. For increasing sub-unit interaction strength, three regimes with
different effects of the membrane are identified.


  Background: Intra-tumour heterogeneity (ITH) is the result of ongoing
evolutionary change within each cancer. The expansion of genetically distinct
sub-clonal populations may explain the emergence of drug resistance and if so
would have prognostic and predictive utility. However, methods for objectively
quantifying ITH have been missing and are particularly difficult to establish
in cancers where predominant copy number variation prevents accurate
phylogenetic reconstruction owing to horizontal dependencies caused by long and
cascading genomic rearrangements.
  Results: To address these challenges we present MEDICC, a method for
phylogenetic reconstruction and ITH quantification based on a Minimum Event
Distance for Intra-tumour Copynumber Comparisons. Using a transducer-based
pairwise comparison function we determine optimal phasing of major and minor
alleles, as well as evolutionary distances between samples, and are able to
reconstruct ancestral genomes. Rigorous simulations and an extensive clinical
study show the power of our method, which outperforms state-of-the-art
competitors in reconstruction accuracy and additionally allows unbiased
numerical quantification of ITH.
  Conclusions: Accurate quantification and evolutionary inference are essential
to understand the functional consequences of ITH. The MEDICC algorithms are
independent of the experimental techniques used and are applicable to both
next-generation sequencing and array CGH data.


  In this work we focused on the kinetics of a one-dimensional Ising system
(1DIS) with constant nearest-neighbor interaction (NNI). The exact solution of
both thermodynamics and kinetics of this system under quasi-chemical
approximation (QCA) had been shown in the literature, and the equilibrium
solution was exact. In this work, it was discussed why QCA applied the best in
the case of 1DIS with constant NNI. Furthermore, extension had been made to
discuss that due to this special reason, perhaps the kinetics of the system
under QCA is the correct steady-state kinetics. Inspired by this observation,
the activity and activity coefficients of the system was studied closely to
re-examine the form of the equation of motion under QCA. A novel concept---the
instantaneous activities and the corresponding instantaneous activity
coefficients---was introduced, and in terms of these quantities the kinetics
seemed to be much simpler and physically more meaningful. The chevron plot of
this system was also discussed and new way of looking at the \emph{rollover} of
chevron plots was presented.


  Network inference approaches are widely used to shed light on regulatory
interplay between molecular players such as genes and proteins. Biochemical
processes underlying networks of interest (e.g. gene regulatory or protein
signalling networks) are generally nonlinear. In many settings, knowledge is
available concerning relevant chemical kinetics. However, existing network
inference methods for continuous, steady-state data are typically rooted in
statistical formulations, which do not exploit chemical kinetics to guide
inference. Herein, we present an approach to network inference for steady-state
data that is rooted in non-linear descriptions of biochemical mechanism. We use
equilibrium analysis of chemical kinetics to obtain functional forms that are
in turn used to infer networks using steady-state data. The approach we propose
is directly applicable to conventional steady-state gene expression or
proteomic data and does not require knowledge of either network topology or any
kinetic parameters. We illustrate the approach in the context of protein
phosphorylation networks, using data simulated from a recent mechanistic model
and proteomic data from cancer cell lines. In the former, the true network is
known and used for assessment, whereas in the latter, results are compared
against known biochemistry. We find that the proposed methodology is more
effective at estimating network topology than methods based on linear models.


  A fitness landscape is a mapping from the space of genetic sequences, which
is modeled here as a binary hypercube of dimension $L$, to the real numbers. We
consider random models of fitness landscapes, where fitness values are assigned
according to some probabilistic rule, and study the statistical properties of
pathways to the global fitness maximum along which fitness increases
monotonically. Such paths are important for evolution because they are the only
ones that are accessible to an adapting population when mutations occur at a
low rate. The focus of this work is on the block model introduced by A.S.
Perelson and C.A. Macken [Proc. Natl. Acad. Sci. USA 92:9657 (1995)] where the
genome is decomposed into disjoint sets of loci (`modules') that contribute
independently to fitness, and fitness values within blocks are assigned at
random. We show that the number of accessible paths can be written as a product
of the path numbers within the blocks, which provides a detailed analytic
description of the path statistics. The block model can be viewed as a special
case of Kauffman's NK-model, and we compare the analytic results to simulations
of the NK-model with different genetic architectures. We find that the mean
number of accessible paths in the different versions of the model are quite
similar, but the distribution of the path number is qualitatively different in
the block model due to its multiplicative structure. A similar statement
applies to the number of local fitness maxima in the NK-models, which has been
studied extensively in previous works. The overall evolutionary accessibility
of the landscape, as quantified by the probability to find at least one
accessible path to the global maximum, is dramatically lowered by the modular
structure.


  Mutational escape from vaccine induced immune responses has thwarted the
development of a successful vaccine against AIDS, whose causative agent is HIV,
a highly mutable virus. Knowing the virus' fitness as a function of its
proteomic sequence can enable rational design of potent vaccines, as this
information can focus vaccine induced immune responses to target mutational
vulnerabilities of the virus. Spin models have been proposed as a means to
infer intrinsic fitness landscapes of HIV proteins from patient-derived viral
protein sequences. These sequences are the product of non-equilibrium viral
evolution driven by patient-specific immune responses, and are subject to
phylogenetic constraints. How can such sequence data allow inference of
intrinsic fitness landscapes? We combined computer simulations and variational
theory \'{a} la Feynman to show that, in most circumstances, spin models
inferred from patient-derived viral sequences reflect the correct rank order of
the fitness of mutant viral strains. Our findings are relevant for diverse
viruses.


  Malaria is a life threatening disease, entirely preventable and treatable,
provided the currently recommended interventions are properly implemented.
These interventions include vector control through the use of
insecticide-treated nets (ITNs). However, ITN possession does not necessarily
translate into use. Human behavior change interventions, including information,
education, communication (IEC) campaigns and post-distribution hang-up
campaigns are strongly recommended. In this paper we consider a recent
mathematical model for the effects of ITNs on the transmission dynamics of
malaria infection, which takes into account the human behavior. We introduce in
this model a supervision control, representing IEC campaigns for improving the
ITN usage. We propose and solve an optimal control problem where the aim is to
minimize the number of infectious humans while keeping the cost low. Numerical
results are provided, which show the effectiveness of the optimal control
interventions.


  The indigenous oak species (Quercus spp.) of the Upper Rhine Valley have
migrated to their current distribution range in the area after the transition
to the Holocene interglacial. Since post-glacial recolonization, they have been
subjected to ecological changes and human impact. By using chloroplast
microsatellite markers (cpSSRs), we provide detailed phylogeographic
information and we address the contribution of natural and human-related
factors to the current pattern of chloroplast DNA (cpDNA) variation. 626
individual trees from 86 oak stands including all three indigenous oak species
of the region were sampled. In order to verify the refugial origin, reference
samples from refugial areas and DNA samples from previous studies with known
cpDNA haplotypes (chlorotypes) were used. Chlorotypes belonging to three
different maternal lineages, corresponding to the three main glacial refugia,
were found in the area. These were spatially structured and highly introgressed
among species, reflecting past hybridization which involved all three
indigenous oak species. Site condition heterogeneity was found among groups of
populations which differed in terms of cpDNA variation. This suggests that
different biogeographic subregions within the Upper Rhine Valley were colonized
during separate post-glacial migration waves. Genetic variation was higher in
Quercus robur than in Quercus petraea, which is probably due to more efficient
seed dispersal and the more pronounced pioneer character of the former species.
Finally, stands of Q. robur established in the last 70 years were significantly
more diverse, which can be explained by the improved transportation ability of
seeds and seedlings for artificial regeneration of stands during this period.


  Genetic analysis was carried out in order to provide insights into
differentiation among populations of two interfertile oak species, Quercus
petraea and Quercus robur. Gene flow between the two species, local adaptations
and speciation processes in general, may leave differential molecular
signatures across the genome. Three interspecific pairs of natural populations
from three ecologically different regions, one in central Europe (SW Germany)
and two in the Balkan Peninsula (Greece and Bulgaria) were sampled. Grouping of
highly informative SSR loci was made according to the component of variation
they express - interspecific or provenance specific. Species and provenance
discriminant loci were characterized based on FSTs. Locus specific FSTs were
tested for deviation from the neutral expectation both within and between
species. Data were then treated separately in a Bayesian analysis of genetic
structure. By using three species discriminant loci, high membership
probability to inferred species groups was achieved. On the other hand,
analysis of genetic structure based on five provenance discriminant loci was
correlated with geographic region and revealed shared genetic variation between
neighbouring Q. petraea and Q. robur. Small sets of highly variable nuclear
SSRs were sufficient to discriminate, either between species or between
provenances. Thus, an effective tool is provided for molecular identification
of both species and provenances. Furthermore, data suggest that a combination
of gene flow and natural selection forms these diversity patterns. Species
discriminant loci might represent genome regions affected by directional
selection, which maintains species identity. Provenance specific loci might
represent genome regions with high interspecific gene flow and common adaptive
patterns to local environmental factors.


  Gene trees record the combination of gene level events, such as duplication,
transfer and loss, and species level events, such as speciation and extinction.
Gene tree-species tree reconciliation methods model these processes by drawing
gene trees into the species tree using a series of gene and species level
events. The reconstruction of gene trees based on sequence alone almost always
involves choosing between statistically equivalent or weakly distinguishable
relationships that could be much better resolved based on a putative species
tree. To exploit this potential for accurate reconstruction of gene trees the
space of reconciled gene trees must be explored according to a joint model of
sequence evolution and gene tree-species tree reconciliation.
  Here we present amalgamated likelihood estimation (ALE), a probabilistic
approach to exhaustively explore all reconciled gene trees that can be
amalgamated as a combination of clades observed in a sample of trees. We
implement ALE in the context of a reconciliation model, which allows for the
duplication, transfer and loss of genes. We use ALE to efficiently approximate
the sum of the joint likelihood over amalgamations and to find the reconciled
gene tree that maximizes the joint likelihood.
  We demonstrate using simulations that gene trees reconstructed using the
joint likelihood are substantially more accurate than those reconstructed using
sequence alone. Using realistic topologies, branch lengths and alignment sizes,
we demonstrate that ALE produces more accurate gene trees even if the model of
sequence evolution is greatly simplified. Finally, examining 1099 gene families
from 36 cyanobacterial genomes we find that joint likelihood-based inference
results in a striking reduction in apparent phylogenetic discord, with 24%, 59%
and 46% percent reductions in the mean numbers of duplications, transfers and
losses.


  Here we present mendelFix, a Perl script for checking Mendelian errors in
genome-wide SNP data of trio designs. The program takes 12-recoded PLINK PED
and MAP files as input to calculate a series of summary statistics for
Mendelian errors, sets missing offspring genotypes that present Mendelian
inconsistencies, and implements a simplistic procedure to infer missing
genotypes using parent information. The program can be easily incorporated in
any pipeline for family-based SNP data analysis, and is distributed as free
software under the GNU General Public License.


  Self-avoiding random walks were performed on protein residue networks.
Compared with protein residue networks with randomized links, the probability
of a walk being successful is lower and the length of successful walks shorter
in (non-randomized) protein residue networks. Fewer successful walks and
shorter successful walks point to higher communication specificity between
protein residues, a conceivably favourable attribute for proteins to have. The
use of random walks instead of shortest paths also produced lower node
centrality, lower edge betweeness and lower edge load for (non-randomized)
protein residue networks than in their respective randomized counterparts. The
implications of these properties for protein residue networks are discussed in
terms of communication congestion and network vulnerability. The randomized
protein residue networks have lower network clustering than the
(non-randomized) protein residue networks. Hence, our findings also shed light
on a hitherto neglected aspect: the importance of high network clustering in
protein residue networks. High clustering increases navigability of a network
for local search and the combination of a local search process on a highly
clustered small-world network topology such as protein residue networks reduces
communication congestion and network vulnerability.


  Small Angle Scattering (SAS) of X-rays or neutrons is an experimental
technique that provides valuable structural information for biological
macromolecules under physiological conditions and with no limitation on the
molecular size. In order to refine molecular structure against experimental SAS
data, ab initio prediction of the scattering profile must be recomputed
hundreds of thousands of times, which involves the computation of the sinc
kernel over all pairs of atoms in a molecule. The quadratic computational
complexity of predicting the SAS profile limits the size of the molecules and
and has been a major impediment for integration of SAS data into structure
refinement protocols. In order to significantly speed up prediction of the SAS
profile we present a general purpose graphical processing unit (GPU) algorithm,
written in OpenCL, for the summation of the sinc kernel (Debye summation) over
all pairs of atoms. This program is an order of magnitude faster than a
parallel CPU algorithm, and faster than an FMM-like approximation method for
certain input domains. We show that our algorithm is currently the fastest
method for performing SAS computation for small and medium size molecules
(around 50000 atoms or less). This algorithm is critical for quick and accurate
SAS profile computation of elongated structures, such as DNA, RNA, and sparsely
spaced pseudo-atom molecules.


  This is an introduction to the special issue titled "Collective behavior and
evolutionary games" that is in the making at Chaos, Solitons & Fractals. The
term collective behavior covers many different phenomena in nature and society.
From bird flocks and fish swarms to social movements and herding effects, it is
the lack of a central planner that makes the spontaneous emergence of sometimes
beautifully ordered and seemingly meticulously designed behavior all the more
sensational and intriguing. The goal of the special issue is to attract
submissions that identify unifying principles that describe the essential
aspects of collective behavior, and which thus allow for a better
interpretation and foster the understanding of the complexity arising in such
systems. As the title of the special issue suggests, the later may come from
the realm of evolutionary games, but this is certainly not a necessity, neither
for this special issue, and certainly not in general. Interdisciplinary work on
all aspects of collective behavior, regardless of background and motivation,
and including synchronization and human cognition, is very welcome.


  Epigenetic modifications to histones may promote either activation or
repression of the transcription of nearby genes. Recent experimental studies
show that the promoters of many lineage-control genes in stem cells have
"bivalent domains" in which the nucleosomes contain both active (H3K4me3) and
repressive (H3K27me3) marks. It is generally agreed that bivalent domains play
an important role in stem cell differentiation, but the underlying mechanisms
remain unclear. Here we formulate a mathematical model to investigate the
dynamic properties of histone modification patterns. We then illustrate that
our modeling framework can be used to capture key features of experimentally
observed combinatorial chromatin states.


  An important step in a phospholipid membrane pore formation by melittin
antimicrobial peptide is a reorientation of the peptide from a surface into a
transmembrane conformation. In this work we perform umbrella sampling
simulations to calculate the potential of mean force (PMF) for the
reorientation of melittin from a surface-bound state to a transmembrane state
and provide a molecular level insight into understanding peptide and lipid
properties that influence the existence of the free energy barrier. The PMFs
were calculated for a peptide to lipid (P/L) ratio of 1/128 and 4/128. We
observe that the free energy barrier is reduced when the P/L ratio increased.
In addition, we study the cooperative effect; specifically we investigate if
the barrier is smaller for a second melittin reorientation, given that another
neighboring melittin was already in the transmembrane state. We observe that
indeed the barrier of the PMF curve is reduced in this case, thus confirming
the presence of a cooperative effect.


  Populations of replicating entities frequently experience sudden or cyclical
changes in environment. We explore the implications of this phenomenon via a
environmental switching parameter in several common evolutionary dynamics
models including the replicator dynamic for linear symmetric and asymmetric
landscapes, the Moran process, and incentive dynamics. We give a simple
relationship between the probability of environmental switching, the relative
fitness gain, and the effect on long term behavior in terms of fixation
probabilities and long term outcomes for deterministic dynamics. We also
discuss cases where the dynamic changes, for instance a population evolving
under a replicator dynamic switching to a best-reply dynamic and vice-versa,
giving Lyapunov stability results.


  We define the incentive process, a natural generalization of the Moran
process incorporating evolutionary updating mechanisms corresponding to
well-known evolutionary dynamics, such as the logit, projection, and best-reply
dynamics. Fixation probabilities and internal stable states are given for a
variety of incentives, including new closed-forms, as well as results relating
fixation probabilities for members of two one-parameter families of incentive
processes. We show that the behaviors of the incentive process can deviate
significantly from the analogous properties of deterministic evolutionary
dynamics in some ways but are similar in others. For example, while the
fixation probabilities change, their ratio remains constant.


  In this short note, a correction is made to the recently proposed solution
[1] to a 1D biased diffusion model for linear DNA translocation and a new
analysis will be given to the data in [1]. It was pointed out [2] by us
recently that this 1D linear translocation model is equivalent to the one that
was considered by Schrodinger [3] for the Enrenhaft-Millikan measurements [4,5]
on electron charge. Here we apply Schrodinger's first-passage-time distribution
formula to the data set in [1]. It is found that Schrodinger's formula can be
used to describe the time distribution of DNA translocation in solid-state
nanopores. These fittings yield two useful parameters: drift velocity of DNA
translocation and diffusion constant of DNA inside the nanopore. The results
suggest two regimes of DNA translocation: (I) at low voltages, there are clear
deviations from Smoluchowski's linear law of electrophoresis [6] which we
attribute to the entropic barrier effects; (II) at high voltages, the
translocation velocity is a linear function of the applied electric field. In
regime II, the apparent diffusion constant exhibits a quadratic dependence on
applied electric field, suggesting a mechanism of Taylor dispersion effect
likely due the electro-osmotic flow field in the nanopore channel. This
analysis yields a dispersion-free diffusion constant value for the segment of
DNA inside the nanopore which is in agreement with Stokes-Einstein theory
quantitatively. The implication of Schrodinger's formula for DNA sequencing is
discussed.


  In numerous systems in biophysics and related fields, scientists measure
(with very smart methods) individual molecules (e.g. biopolymers (proteins,
DNA, RNA, etc), nano - crystals, ion channels), aiming at finding a model from
the data. But the noise is not solved accurately in not so few cases and this
may lead to misleading models. Here, we solve the noise. We consider two state
photon trajectories from any on off kinetic scheme (KS): the process emitting
photons with a rate {\gamma}on when it is in the on state, and emitting with a
rate {\gamma}off when it is in the off state. We develop a filter that removes
the noise resulting in clean data also in cases where binning fails. The filter
is a numerical algorithm with various new statistical treatments. It is based
on a new general likelihood function developed here, with observable dependent
form. The filter can solve the noise, in the detectable region of the rate
space: that is, we also find a region where the data is "too" noisy.
Consistency tests will find the region's type from the data. If the data is
ruled "too noisy", binning obviously fails, and one should apply simpler
methods on the raw data and realizing that the extracted information is
partial. We show that not applying the filter while cleaning results in
erroneous rates. This filter (with minor adjustments) can solve the noise in
any discrete state trajectories, yet extensions are needed in "tackling" the
noise from other data, e.g. continuous data and FRET data.
  The filter developed here is complementary with our previous projects in this
field, where we have solved clean two state data with the development of
reduced dimensions forms (RDFs): only the combined procedures enabling building
the most accurate model from noisy trajectories from single molecules


  Hybridization and introgression between Quercus alnifolia Poech and Q.
coccifera L. is studied by analyzing morphological traits, nuclear and
chloroplast DNA markers. The study site is a mixed stand on Troodos Mountains
(Cyprus) and the analyzed material includes both adult trees and progenies of
specific mother trees. Multivariate analysis of morphological traits shows that
the two species can be well distinguished using simple leaf morphometric
parameters. A lower genetic diversity in Q. alnifolia than in Q. coccifera and
a high interspecific differentiation between the two species are supported by
an analysis of nuclear and chloroplast microsatellites. The intermediacy of the
four designated hybrids is verified by both leaf morphometric and genetic data.
Analysis of progeny arrays provides evidence that interspecific crossings are
rare. This finding is further supported by limited introgression of chloroplast
genomes. Reproductive barriers (e.g. asynchronous phenology, post-zygotic
incompatibilities) might account for this result. A directionality of
interspecific gene flow is indicated by a genetic assignment analysis of
effective pollen clouds with Q. alnifolia acting as pollen donor. Differences
in flowering phenology and species distribution in the stand may have
influenced the direction of gene flow and the genetic differentiation among
effective pollen clouds of different mother trees within species.


  The evergreen Quercus alnifolia and Q. coccifera form the only interfertile
pair of oak species growing in Cyprus. Hybridization between the two species
has already been observed and studied morphologically. However, little evidence
exists about the extent of genetic introgression. In the present study, we
aimed to study the effects of introgressive hybridization mutually on both
chloroplast and nuclear genomes. We sampled both pure and mixed populations of
Q. alnifolia and Q. coccifera from several locations across their distribution
area in Cyprus. We analyzed the genetic variation within and between species by
conducting Analysis of Molecular Variance (AMOVA) based on nuclear
microsatellites. Population genetic structure and levels of admixture were
studied by means of a Bayesian analysis (STRUCTURE simulation analysis).
Chloroplast DNA microsatellites were used for a spatial analysis of genetic
barriers. The main part of the nuclear genetic variation was explained by
partition into species groups. High interspecific differentiation and low
admixture of nuclear genomes, both in pure and mixed populations, support
limited genetic introgression between Q. alnifolia and Q. coccifera in Cyprus.
On the contrary, chloroplast DNA haplotypes were shared between the species and
were locally structured suggesting cytoplasmic introgression. Occasional
hybridization events followed by backcrossings with both parental species might
lead to this pattern of genetic differentiation.


  We define and study a class of (random) Boolean constraint satisfaction
problems representing minimal feasibility constraints for networks of chemical
reactions. The constraints we consider encode, respectively, for hard
mass-balance conditions (where the consumption and production fluxes of each
chemical species are matched) and for soft mass-balance conditions (where a net
production of compounds is in principle allowed). We solve these constraint
satisfaction problems under the Bethe approximation and derive the
corresponding Belief Propagation equations, that involve 8 different messages.
The statistical properties of ensembles of random problems are studied via the
population dynamics methods. By varying a chemical potential attached to the
activity of reactions, we find first order transitions and strong hysteresis,
suggesting a non-trivial structure in the space of feasible solutions.


  Cell-to-cell variations in protein abundance in clonal cell populations are
ubiquitous in living systems. Since protein composition determines responses in
individual cells, it stands to reason that the variations themselves are
subject to selective pressures. But the functional role of these cell-to-cell
differences is not well understood. One way to tackle questions regarding
relationships between form and function is to perturb the form (e.g., change
the protein abundances) and observe the resulting changes in some function.
Here we take on the form-function relationship from the inverse perspective,
asking instead what specific constraints on cell-to-cell variations in protein
abundance are imposed by a given functional phenotype. We develop a maximum
entropy (MaxEnt) based approach to posing questions of this type, and
illustrate the method by application to the well characterized chemotactic
response in Escherichia coli (E. coli). We find that full determination of
observed cell-to-cell variations in protein abundances is not inherent in
chemotaxis itself, but in fact appears to be jointly imposed by the chemotaxis
program in conjunction with other factors, e.g., the protein synthesis
machinery and/or additional non-chemotactic cell functions such as cell
metabolism. These results illustrate the power of MaxEnt as a tool for the
investigation of relationships between biological form and function.


  For many taxa, the current high rates of extinction are likely to result in a
significant loss of biodiversity. The evolutionary heritage of biodiversity is
frequently quantified by a measure called phylogenetic diversity (PD). We
predict the loss of PD under a wide class of phylogenetic tree models, where
speciation rates and extinction rates may be time-dependent, and assuming
independent random species extinctions at the present. We study the loss of PD
when $K$ contemporary species are selected uniformly at random from the $N$
extant species as the surviving taxa, while the remaining $N-K$ become extinct.
We consider two models of species sampling, the so-called field of bullets
model, where each species independently survives the extinction event at the
present with probability $p$, and a model for which the number of surviving
species is fixed.
  We provide explicit formulae for the expected remaining PD in both models,
conditional on $N=n$, conditional on $K=k$, or conditional on both events. When
$N=n$ is fixed, we show the convergence to an explicit deterministic limit of
the ratio of new to initial PD, as $n\to\infty$, both under the field of
bullets model, and when $K=k_n$ is fixed and depends on $n$ in such a way that
$k_n/n$ converges to $p$. We also prove the convergence of this ratio as
$T\to\infty$ in the supercritical, time-homogeneous case, where $N$
simultaneously goes to $\infty$, thereby strengthening previous results of
Mooers et al. (2012).


  We live within an increasingly technological, information-laden environment
for the first time in human evolution.This subjects us, and will continue to
subject us in an accelerating fashion, to an unremitting exposure to meaningful
information that requires action. Directly dependent upon this new environment
are novel evolutionary pressures, which can modify existing resource allocation
mechanisms and may eventually favor the survival of somatic cells,particularly
neurons, at the expense of germ line cells. Here it is argued that persistent,
structured information-sharing in both virtual and real domains, leads to
increased biological complexity and functionality, which reflects upon human
survival characteristics. Certain immortalisation mechanisms currently employed
by germ cells may thus need to be downgraded in order to enable somatic cells
to manage these new energy demands placed by our modern environment. Relevant
concepts from a variety of disciplines such as the evolution of complex
adaptive systems, information theory, digital hyper-connectivity, and cell
immortalisation will be reviewed.


  We reconsider the Moran model in continuous time with population size $N$,
two allelic types, and selection. We introduce a new particle representation,
which we call the labelled Moran model, and which has the same distribution of
type frequencies as the original Moran model, provided the initial values are
chosen appropriately. In the new model, individuals are labelled $1,2, \dots,
N$; neutral resampling events may take place between arbitrary labels, whereas
selective events only occur in the direction of increasing labels. With the
help of elementary methods only, we not only recover fixation probabilities,
but also obtain detailed insight into the number and nature of the selective
events that play a role in the fixation process forward in time.


  The notion of energy landscapes provides conceptual tools for understanding
the complexities of protein folding and function. Energy Landscape Theory
indicates that it is much easier to find sequences that satisfy the "Principle
of Minimal Frustration" when the folded structure is symmetric (Wolynes, P. G.
Symmetry and the Energy Landscapes of Biomolecules. Proc. Natl. Acad. Sci.
U.S.A. 1996, 93, 14249-14255). Similarly, repeats and structural mosaics may be
fundamentally related to landscapes with multiple embedded funnels. Here we
present analytical tools to detect and compare structural repetitions in
protein molecules. By an exhaustive analysis of the distribution of structural
repeats using a robust metric we define those portions of a protein molecule
that best describe the overall structure as a tessellation of basic units. The
patterns produced by such tessellations provide intuitive representations of
the repeating regions and their association towards higher order arrangements.
We find that some protein architectures can be described as nearly periodic,
while in others clear separations between repetitions exist. Since the method
is independent of amino acid sequence information we can identify structural
units that can be encoded by a variety of distinct amino acid sequences.


  Proteins experience a wide variety of conformational dynamics that can be
crucial for facilitating their diverse functions. How is the intrinsic
flexibility required for these motions encoded in their three-dimensional
structures? Here, the overall flexibility of a protein is demonstrated to be
tightly coupled to the total amount of surface area buried within its fold. A
simple proxy for this, the relative solvent accessible surface area (Arel),
therefore shows excellent agreement with independent measures of global protein
flexibility derived from various experimental and computational methods.
Application of Arel on a large scale demonstrates its utility by revealing
unique sequence and structural properties associated with intrinsic
flexibility. In particular, flexibility as measured by Arel shows little
correspondence with intrinsic disorder, but instead tends to be associated with
multiple domains and increased {\alpha}- helical structure. Furthermore, the
apparent flexibility of monomeric proteins is found to be useful for
identifying quaternary structure errors in published crystal structures. There
is also a strong tendency for the crystal structures of more flexible proteins
to be solved to lower resolutions. Finally, local solvent accessibility is
shown to be a primary determinant of local residue flexibility. Overall this
work provides both fundamental mechanistic insight into the origin of protein
flexibility and a simple, practical method for predicting flexibility from
protein structures.


  Simple RNA viruses efficiently encapsulate their genome into a nano-sized
protein shell: the capsid. Spontaneous co-assembly of the genome and the capsid
proteins is driven predominantly by electrostatic interactions between the
negatively charged RNA and the positively charged inner capsid wall. Using
field theoretic formulation we show that the inherently branched RNA secondary
structure allows viruses to maximize the amount of encapsulated genome and make
assembly more efficient, allowing viral RNAs to out-compete cellular RNAs
during replication in infected host cells.


  In this paper, we numerically investigate the robustness of cooperation
clusters in prisoner's dilemma played on scale-free networks, where the network
topologies change by continuous removal and addition of nodes. Each removal and
addition can be either random or intentional. We therefore have four different
strategies in changing network topology: random removal and random addition
(RR), random removal and preferential addition (RP), targeted removal and
random addition (TR), and targeted removal and preferential addition (TP). We
find that cooperation clusters are most fragile against TR, while they are most
robust against RP, even for large values of the temptation coefficient for
defection. The effect of the degree mixing pattern of the network is not the
primary factor for the robustness of cooperation under continuous change in
network topology, which is quite different from the cases observed in static
networks. Cooperation clusters become more robust as the number of links of
hubs occupied by cooperators increase. Our results might infer the fact that a
huge variety of individuals is needed for maintaining global cooperation in
social networks in the real world where each node representing an individual is
constantly removed and added.


  1. In recent years many studies have investigated the relationships between
different aspects of the immune system and ecology in various organisms. Yet,
it remains unclear why individuals differ in their ability to mount an immune
response against various antigens (often referred to as 'immunocompetence').
Different kinds of trade-offs may be involved and costs of mounting the immune
response often lead to condition dependent effects. 2. We investigated how
variation in condition, morphology and genetic variables influenced the amount
of antibodies produced against two novel antigens in a migrating bird, the
great snipe (Gallinago media). 3. We found no evidence for condition dependence
of the antibody response and no effect of MHC genetics. There was, however, a
weak negative correlation between body size and the amount of antibody
production, which may indicate a trade-off between growth and immune response
in this species.


  The aim of our study is to describe the dynamics of ant battles, with
reference to laboratory experiments, by means of a chemical stochastic model.
We focus on ants behavior as an interesting topic in order to predict the
ecological evolution of invasive species and their spreading. In our work we
want to describe the interactions between two groups of different ant species
with different war strategies. Our model considers the single ant individuals
and fighting groups in a way similar to atoms and molecules, respectively,
considering that ant fighting groups remain stable for a relative long time.
Starting from a system of differential non-linear equations (DE), derived from
the chemical reactions, we obtain a mean field description of the system. The
DE approach is valid when the number of individuals of each species is large in
the considered unit, while we consider battles of at most 10 vs. 10
individuals, due to the difficulties in following the individual behavior in a
large assembly. Therefore, we also adapt a Gillespie algorithm to reproduce the
fluctuations around the mean field. The DE scheme is exploited to characterize
the stochastic model. The set of parameters of chemical equations, obtained
using a minimization algorithm, are used by the Gillespie algorithm to generate
the stochastic trajectories. We then fit the stochastic paths with the DE, in
order to analyze the variability of the parameters and their variance. Finally,
we estimate the goodness of the applied methodology and we confirm that the
stochastic approach must be considered for a correct description of the ant
fighting dynamics. With respect to other war models, our chemical one considers
all phases of the battle and not only casualties. Thus, we can count on more
experimental data, but we also have more parameters to fit. In any case, our
model allows a much more detailed description of the fights.


  Candida albicans is responsible for a number of life-threatening infections
and causes considerable morbidity and mortality in immunocompromised patients.
Previous studies of C. albicans pathogenesis have suggested several steps must
occur before virulent infection, including early adhesion, invasion, and late
tissue damage. However, the mechanism that triggers C. albicans transformation
from yeast to hyphae form during infection has yet to be fully elucidated. This
study used a systems biology approach to investigate C. albicans infection in
zebrafish. The surviving fish were sampled at different post-infection time
points to obtain time-lapsed, genome-wide transcriptomic data from both
organisms, which were accompanied with in sync histological analyses. Principal
component analysis (PCA) was used to analyze the dynamic gene expression
profiles of significant variations in both C. albicans and zebrafish. The
results categorized C. albicans infection into three progressing phases:
adhesion, invasion, and damage. Such findings were highly supported by the
corresponding histological analysis. Furthermore, the dynamic interspecies
transcript profiling revealed that C. albicans activated its filamentous
formation during invasion and the iron scavenging functions during the damage
phases, whereas zebrafish ceased its iron homeostasis function following
massive hemorrhage during the later stages of infection. This was followed by
massive hemorrhaging toward the end stage of infection. Most of the immune
related genes were expressed as the infection progressed from invasion to the
damage phase. Such global, inter-species evidence of virulence-immune and iron
competition dynamics during C. albicans infection could be crucial in
understanding control fungal pathogenesis.


  There is a widening recognition that cancer cells are products of complex
developmental processes. Carcinogenesis and metastasis formation are
increasingly described as systems-level, network phenomena. Here we propose
that malignant transformation is a two-phase process, where an initial increase
of system plasticity is followed by a decrease of plasticity at late stages of
carcinogenesis as a model of cellular learning. We describe the hallmarks of
increased system plasticity of early, tumor initiating cells, such as increased
noise, entropy, conformational and phenotypic plasticity, physical
deformability, cell heterogeneity and network rearrangements. Finally, we argue
that the large structural changes of molecular networks during cancer
development necessitate a rather different targeting strategy in early and late
phase of carcinogenesis. Plastic networks of early phase cancer development
need a central hit, while rigid networks of late stage primary tumors or
established metastases should be attacked by the network influence strategy,
such as by edgetic, multi-target, or allo-network drugs. Cancer stem cells need
special diagnosis and targeting, since their dormant and rapidly proliferating
forms may have more rigid, or more plastic networks, respectively. The
extremely high ability to change their rigidity/plasticity may be a key
differentiating hallmark of cancer stem cells. The application of early
stage-optimized anti-cancer drugs to late-stage patients may be a reason of
many failures in anti-cancer therapies. Our hypotheses presented here underlie
the need for patient-specific multi-target therapies applying the correct ratio
of central hits and network influences -- in an optimized sequence.


  The Susceptible-Infected-Susceptible (SIS) model is one of the simplest
memoryless system for describing information/epidemic spreading phenomena with
competing creation and spontaneous annihilation reactions. The effect of
quenched disorder on the dynamical behavior has recently been compared to
quenched mean-field (QMF) approximations in scale-free networks. QMF can take
into account topological heterogeneity and clustering effects of the activity
in the steady state by spectral decomposition analysis of the adjacency matrix.
Therefore, it can provide predictions on possible rare-region effects, thus on
the occurrence of slow dynamics. I compare QMF results of SIS with simulations
on various large dimensional graphs. In particular, I show that for Erd\H
os-R\'enyi graphs this method predicts correctly the epidemic threshold and the
rare-region effects. Griffiths Phases emerge if the graph is fragmented or if
we apply strong, exponentially suppressing weighting scheme on the edges. The
latter model describes the connection time distributions in the face-to-face
experiments. In case of generalized Barab\'asi-Albert type of networks with
aging connections strong rare-region effects and numerical evidence for
Griffiths Phase dynamics are shown.


  The intricate pattern of chemical modifications on DNA and histones, the
"histone code", is considered to be a key gene regulation factor. Multivalency
is seen by many as an essential instrument to transmit the "encoded"
information to the transcription machinery via multi-domain effector proteins
and chromatin-associated complexes. However, as examples of multivalent histone
engagement accumulate, an apparent contradiction is emerging. The isolated
effector domains are notably weak binders, thus it is often asserted that the
entropic cost of orienting multiple domains can be "prepaid" by a rigid tether.
Meanwhile, evidence suggests that the tethers are largely disordered and offer
little rigidity. Here we consider a mechanism to "prepay" the entropic costs of
orienting the domains for binding, not through rigidity of the tether but
through the careful spacing of the modifications on chromatin. An all-atom
molecular dynamics study of the most fully characterized multivalent chromatin
effector conforms to the conditions for an optimal free-energy payout, as
predicted by the model discussed here.


  The reconstruction of phylogenetic trees based on viral genetic sequence data
sequentially sampled from an epidemic provides estimates of the past
transmission dynamics, by fitting epidemiological models to these trees. To our
knowledge, none of the epidemiological models currently used in phylogenetics
can account for recovery rates and sampling rates dependent on the time elapsed
since transmission.
  Here we introduce an epidemiological model where infectives leave the
epidemic, either by recovery or sampling, after some random time which may
follow an arbitrary distribution.
  We derive an expression for the likelihood of the phylogenetic tree of
sampled infectives under our general epidemiological model. The analytic
concept developed in this paper will facilitate inference of past
epidemiological dynamics and provide an analytical framework for performing
very efficient simulations of phylogenetic trees under our model. The main idea
of our analytic study is that the non-Markovian epidemiological model giving
rise to phylogenetic trees growing vertically as time goes by, can be
represented by a Markovian "coalescent point process" growing horizontally by
the sequential addition of pairs of coalescence and sampling times.
  As examples, we discuss two special cases of our general model, namely an
application to influenza and an application to HIV. Though phrased in
epidemiological terms, our framework can also be used for instance to fit
macroevolutionary models to phylogenies of extant and extinct species,
accounting for general species lifetime distributions.


  Ecosystems can undergo sudden shifts to undesirable states, but recent
studies with simple single species ecosystems have demonstrated that advance
warning can be provided by the slowing down of population dynamics near a
tipping point. However, it is not clear how this effect of critical slowing
down will manifest in ecosystems with strong interactions between their
components. Here we probe the dynamics of an experimental producer parasite
ecosystem as it approaches a catastrophic collapse. Surprisingly, the producer
population grows in size as the environment deteriorates, highlighting that
population size can be a misleading measure of ecosystem stability. By
analyzing the oscillatory producer parasite dynamics for over ~100 generations
in multiple environmental conditions, we found that the collective ecosystem
dynamics slows down as the tipping point is approached. Analysis of the coupled
dynamics of interacting populations may therefore be necessary to provide
advance warning of collapse in complex communities.


  We introduce Sashimi plots, a quantitative multi-sample visualization of mRNA
sequencing reads aligned to gene annotations. Sashimi plots are made using
alignments (stored in the SAM/BAM format) and gene model annotations (in GFF
format), which can be custom-made by the user or obtained from databases such
as Ensembl or UCSC. We describe two implementations of Sashimi plots: (1) a
stand-alone command line implementation aimed at making customizable
publication quality figures, and (2) an implementation built into the
Integrated Genome Viewer (IGV) browser, which enables rapid and dynamic
creation of Sashimi plots for any genomic region of interest, suitable for
exploratory analysis of alternatively spliced regions of the transcriptome.
Isoform expression estimates outputted by the MISO program can be optionally
plotted along with Sashimi plots. Sashimi plots can be used to quickly screen
differentially spliced exons along genomic regions of interest and can be used
in publication quality figures. The Sashimi plot software and documentation is
available from: http://genes.mit.edu/burgelab/miso/docs/sashimi.html


  We investigate the properties of a Wright-Fisher diffusion process started
from frequency x at time 0 and conditioned to be at frequency y at time T. Such
a process is called a bridge. Bridges arise naturally in the analysis of
selection acting on standing variation and in the inference of selection from
allele frequency time series. We establish a number of results about the
distribution of neutral Wright-Fisher bridges and develop a novel rejection
sampling scheme for bridges under selection that we use to study their
behavior.


  Blue Mesa Reservoir is the largest body of water in Colorado and is located
on the western slope of the Rocky Mountains at an elevation of 7520 feet. Blue
Mesa Reservoir contains recreationally important populations of lake trout
(Salvelinus namaycush), kokanee salmon (Oncorhynchus nerka), rainbow trout
(Onocorhynchus mykiss), and brown trout (Salmo trutta). A management challenge
in recent years has been the overpopulation of lake trout, which has led to a
steep decline in abundance of kokanee salmon (a preferred food source) followed
by a steep decline in body condition of the lake trout. Relative weight, Wr, is
commonly used to assess body condition and prey availability of fish in a
variety of ecosystems, but its seasonal variability has not been well
documented in cold water species. Relative weight is 100 times the actual
weight of a fish divided by expected weight (or standard weight) of that
species at the measured total length. The present study uses multi-year,
multi-season length and weight data provided by the Colorado Department of
Parks and Wildlife to quantify seasonal variations in mean relative weight in
these four species. It is found that lake trout consistently lose significant
body condition between the spring and the fall of a given year, with a mean
loss in relative weight of 11.7 +/- 2.4. In contrast, kokanee salmon tend to
gain significant body condition between spring and fall with a mean gain in
relative weight of 5.2 +/- 2.9. Rainbow trout tend to lose body condition
between summer and fall with a mean loss of relative weight of 7.4 +/- 1.7.
Brown trout tend to gain body condition between spring and fall with a mean
gain of 5.6 +/- 1.8. These findings can be interpreted in terms of the feeding
habits, forage availability, and metabolic activity of each species.


  This analysis of the variations of brown tiger prawn (Penaeus esculentus)
catch in the Moreton Bay multispecies trawl fishery estimated catchability
using a delay difference model. It integrated several factors responsible for
variations in catchability: targeting of fishing effort, increasing fishing
power and changing availability. An analysis of covariance was used to define
fishing events targeted at brown tiger prawns. A general linear model estimated
inter-annual variations of fishing power. Temperature induced changes in prawn
behaviour played an important role in the dynamic of this fishery. Maximum
likelihood estimates of targeted catchability ($3.92 \pm 0.40 \ 10^{-4}$
boat-days$^{-1}$) were twice as large as non-targeted catchability ($1.91 \pm
0.24 \ 10^{-4}$ boat-days$^{-1}$). The causes of recent decline in fishing
effort in this fishery were discussed.


  High-throughput sequencing is now regularly used for studies of the
transcriptome (RNA-seq), particularly for comparisons among experimental
conditions. For the time being, a limited number of biological replicates are
typically considered in such experiments, leading to low detection power for
differential expression. As their cost continues to decrease, it is likely that
additional follow-up studies will be conducted to re-address the same
biological question. We demonstrate how p-value combination techniques
previously used for microarray meta-analyses can be used for the differential
analysis of RNA-seq data from multiple related studies. These techniques are
compared to a negative binomial generalized linear model (GLM) including a
fixed study effect on simulated data and real data on human melanoma cell
lines. The GLM with fixed study effect performed well for low inter-study
variation and small numbers of studies, but was outperformed by the
meta-analysis methods for moderate to large inter-study variability and larger
numbers of studies. To conclude, the p-value combination techniques illustrated
here are a valuable tool to perform differential meta-analyses of RNA-seq data
by appropriately accounting for biological and technical variability within
studies as well as additional study-specific effects. An R package metaRNASeq
is available on the R Forge.


  The design of photo-active collagen systems is presented as a basis for
establishing biomimetic materials with varied network architecture and
programmable macroscopic properties. Following in-house isolation of type I
collagen, reaction with vinyl-bearing compounds of varied backbone rigidity,
i.e. 4-vinylbenzyl chloride (4VBC) and glycidyl methacrylate (GMA), was carried
out. TNBS colorimetric assay, 1H-NMR and ATR-FTIR confirmed covalent and
tunable functionalization of collagen lysines. Depending on the type and extent
of functionalization, controlled stability and thermal denaturation of triple
helices were observed via circular dichroism (CD), whereby the hydrogen-bonding
capability of introduced moieties was shown to play a major role. Full gel
formation was observed following photo-activation of functionalized collagen
solutions. The presence of a covalent network only slightly affected collagen
triple helix conformation (as observed by WAXS and ATR-FTIR), confirming the
structural organization of functionalized collagen precursors. Photo-activated
hydrogels demonstrated an increased denaturation temperature (DSC) with respect
to native collagen, suggesting that the formation of the covalent network
successfully stabilized collagen triple helices. Moreover, biocompatibility and
mechanical competence of obtained hydrogels were successfully demonstrated
under physiologically-relevant conditions. These results demonstrate that this
novel synthetic approach enabled the formation of biocompatible collagen
systems with defined network architecture and programmable macroscopic
properties, which can only partially be obtained with current synthetic
methods.


  As a consequence of the rugged landscape of RNA molecules their folding is
described by the kinetic partitioning mechanism according to which only a small
fraction ($\phi_F$) reaches the folded state while the remaining fraction of
molecules is kinetically trapped in misfolded intermediates. The transition
from the misfolded states to the native state can far exceed biologically
relevant time. Thus, RNA folding in vivo is often aided by protein cofactors,
called RNA chaperones, that can rescue RNAs from a multitude of misfolded
structures. We consider two models, based on chemical kinetics and chemical
master equation, for describing assisted folding. In the passive model,
applicable for class I substrates, transient interactions of misfolded
structures with RNA chaperones alone are sufficient to destabilize the
misfolded structures, thus entropically lowering the barrier to folding. For
this mechanism to be efficient the intermediate ribonucleoprotein (RNP) complex
between collapsed RNA and protein cofactor should have optimal stability. We
also introduce an active model (suitable for stringent substrates with small
$\phi_F$), which accounts for the recent experimental findings on the action of
CYT-19 on the group I intron ribozyme, showing that RNA chaperones does not
discriminate between the misfolded and the native states. In the active model,
the RNA chaperone system utilizes chemical energy of ATP hydrolysis to
repeatedly bind and release misfolded and folded RNAs, resulting in substantial
increase of yield of the native state. The theory outlined here shows, in
accord with experiments, that in the steady state the native state does not
form with unit probability.


  There is great scientific and popular interest in understanding the genetic
history of populations in the Americas. We wish to understand when different
regions of the continent were inhabited, where settlers came from, and how
current inhabitants relate genetically to earlier populations. Recent studies
unraveled parts of the genetic history of the continent using genotyping arrays
and uniparental markers. The 1000 Genomes Project provides a unique opportunity
for improving our understanding of population genetic history by providing over
a hundred sequenced low coverage genomes and exomes from Colombian (CLM),
Mexican-American (MXL), and Puerto Rican (PUR) populations. Here, we explore
the genomic contributions of African, European, and Native American ancestry to
these populations. Estimated Native American ancestry is 48% in MXL, 25% in
CLM, and 13% in PUR. Native American ancestry in PUR is most closely related to
populations surrounding the Orinoco River basin, confirming the Southern
America ancestry of the Ta\'ino people of the Caribbean. We present new methods
to estimate the allele frequencies in the Native American fraction of the
populations, and model their distribution using a demographic model for three
ancestral Native American populations. These ancestral populations likely split
in close succession: the most likely scenario, based on a peopling of the
Americas 16 thousand years ago (kya), supports that the MXL Ancestors split
12.2kya, with a subsequent split of the ancestors to CLM and PUR 11.7kya. The
model also features effective populations of 62,000 in Mexico, 8,700 in
Colombia, and 1,900 in Puerto Rico. Modeling Identity-by-descent and ancestry
tract length, we show that post-contact populations differ markedly in their
effective sizes and migration patterns, with Puerto Rico showing the smallest
effective size and the earlier migration from Europe.


  We consider voltage-driving DNA translocation through a nanopore in the
present study. By assuming the DNA is coaxial with the cylindrical nanopore, a
hydrodynamic model for determining effective force on a single DNA molecule in
a nanopore was presented, in which density functional theory (DFT) combined
with the continuum Navier-Stokes (NS) equations is utilized to investigate
electro-osmotic flow and the viscous drag force acting on the DNA inside a
nanopore. Surface charge on the walls of the nanopore is also taken into
account in our model. The consistence between our calculation and the previous
experimental measurement indicates that the present theoretical model is an
effective tool to predict the hydrodynamic resistance on DNA. Results show that
charge inversion, which cannot be obtained by the Poisson-Boltzmann (PB) model,
will reduce electro-osmotic velocity, or even lead to flow reversal for higher
salt concentration. This is helpful to raise the effective force profoundly in
the overscreening region.


  How mutations accumulate in genomes is the central question of molecular
evolution theories, however our understanding of this process is far from
complete. Drake's rule is a notoriously universal property of genomes from
microbes to mammals - the number of (functional) mutations per-genome
per-generation is approximately constant within a phylum, despite orders of
magnitude differences in genome sizes and diverse populations properties. So
far there is no concise explanation for this phenomenon. A formal model of
storage of genetic information suggests that a genome of any species operates
near its maximum informational storage capacity, and the mutation rate is near
its upper limit, providing a simple explanation for the rule with minimal
assumptions.


  We investigate a type of bistability where noise not only causes transitions
between stable states, but also constructs the states themselves. We focus on
the experimentally well-studied system of ants choosing between two food
sources to illustrate the essential points, but the ideas are more general. The
mean time for switching between the two bistable states of the system is
calculated. This suggests a procedure for estimating, in a real system, the
critical population size above which bistability ceases to occur.


  Various subsets of self-avoiding walks naturally appear when investigating
existing methods designed to predict the 3D conformation of a protein of
interest. Two such subsets, namely the folded and the unfoldable self-avoiding
walks, are studied computationally in this article. We show that these two sets
are equal and correspond to the whole $n$-step self-avoiding walks for
$n\leqslant 14$, but that they are different for numerous $n \geqslant 108$,
which are common protein lengths. Concrete counterexamples are provided and the
computational methods used to discover them are completely detailed. A tool for
studying these subsets of walks related to both pivot moves and proteins
conformations is finally presented.


  Using theoretical arguments and extensive Monte Carlo (MC) simulations of a
coarse-grained three-dimensional off-lattice model of a \beta-hairpin, we
demonstrate that the equilibrium critical force, $F_c$, needed to unfold the
biopolymer increases non-linearly with increasing volume fraction occupied by
the spherical macromolecular crowding agent. Both scaling arguments and MC
simulations show that the critical force increases as $F_c \approx
\varphi_c^{\alpha}$. The exponent $\alpha$ is linked to the Flory exponent
relating the size of the unfolded state of the biopolymer and the number of
amino acids. The predicted power law dependence is confirmed in simulations of
the dependence of the isothermal extensibility and the fraction of native
contacts on $\varphi_c$. We also show using MC simulations that $F_c$ is
linearly dependent on the average osmotic pressure ($\mathrm{P}$) exerted by
the crowding agents on the \beta-hairpin. The highly significant linear
correlation coefficient of 0.99657 between $F_c$ and $\mathrm{P}$ makes it
straightforward to predict the dependence of the critical force on the density
of crowders. Our predictions are amenable to experimental verification using
Laser Optical Tweezers.


  Protein structure prediction is one of the most important problems in
computational biology. The most successful computational approach, also called
template-based modeling, identifies templates with solved crystal structures
for the query proteins and constructs three dimensional models based on
sequence/structure alignments. Although substantial effort has been made to
improve protein sequence alignment, the accuracy of alignments between
distantly related proteins is still unsatisfactory. In this thesis, I will
introduce a number of statistical machine learning methods to build accurate
alignments between a protein sequence and its template structures, especially
for proteins having only distantly related templates. For a protein with only
one good template, we develop a regression-tree based Conditional Random Fields
(CRF) model for pairwise protein sequence/structure alignment. By learning a
nonlinear threading scoring function, we are able to leverage the correlation
among different sequence and structural features. We also introduce an
information-theoretic measure to guide the learning algorithm to better exploit
the structural features for low-homology proteins with little evolutionary
information in their sequence profile. For a protein with multiple good
templates, we design a probabilistic consistency approach to thread the protein
to all templates simultaneously. By minimizing the discordance between the
pairwise alignments of the protein and templates, we are able to construct a
multiple sequence/structure alignment, which leads to better structure
predictions than any single-template based prediction.


  We study the propagation of tension caused by an external force along a long
polymeric molecule in two different settings, namely along a free polymer in 3d
space being pulled from one end, and along a pre-stretched circular polymer,
confined in a narrow circular tube. We show that in both cases, the tension
propagation is governed by a diffusion equation, and in particular, the tension
front propagates as $t^{1/2}$ along the contour of the chain. The results are
confirmed numerically, and by molecular dynamics simulations in the case of the
3d polymer. We also compare our results with the previously suggested ones for
the translocation setting, and discuss why tension propagation is slower in
that case.


  Motivated by DNA electrophoresis near a nanopore, we consider the flow field
around an "elongated jet", a long thin source which injects momentum into a
liquid. This solution qualitatively describes the electro-osmotic flow around a
long rigid polymer, where due to electrohydrodynamic coupling, the solvent
receives momentum from the electric field. Based on the qualitative behavior of
the elongated jet solution, we develop a coarse-grained scheme which reproduces
the known theoretical results regarding the electrophoretic behavior of a long
rigid polymer and a polymer coil in a uniform field, which we then exploit to
analyze the electrophoresis of a polymer coil in the non-uniform field near a
nanopore.


  Based on our formulation of the DNA electrophoresis near a pore [P.
Rowghanian and A. Y. Grosberg, Phys. Rev. E 87, 042723 (2013)], we address the
electrophoretic DNA capture into a nanopore as a steady-state process of
particle absorption to a sink placed on top of an energy barrier. Reproducing
the previously observed diffusion-limited and barrier-limited regimes as two
different limits of the particle absorption process and matching the data, our
model suggests a slower growth of the capture rate with the DNA length for very
large DNA molecules than the previous model, motivating more experiments beyond
the current range of electric field and DNA length. At moderately weak electric
fields, our model predicts a different effect, stating that the DNA length
dependence of the capture rate first disappears as the field is reduced and
eventually reverses to a decreasing trend with $N$.


  Longstanding mechanistic questions about the role of protecting osmolyte
trimethylamine N- oxide (TMAO) which favors protein folding and the denaturing
osmolyte urea are addressed by studying their effects on the folding of
uncharged polymer chains. Using atomistic molecular dynamics simulations, we
show that 1-M TMAO and 7-M urea solutions act dramatically differently on these
model polymer chains. Their behaviors are sensitive to the strength of the
attractive dispersion interactions of the chain with its environment: when
these dispersion interactions are high enough, TMAO suppresses the formation of
extended conformations of the hydrophobic polymer as compared to water, while
urea promotes formation of extended conformations. Similar trends are observed
experimentally on real protein systems. Quite surprisingly, we find that both
protecting and denaturing osmolytes strongly interact with the polymer,
seemingly in contrast with existing explanations of the osmolyte effect on
proteins. We show that what really matters for a protective osmolyte is its
effective depletion as the polymer conformation changes, which leads to a
negative change in the preferential binding coefficient. For TMAO, there is a
much more favorable free energy of insertion of a single osmolyte near
collapsed conformations of the polymer than near extended conformations. By
contrast, urea is preferentially stabilized next to the extended conformation
and thus has a denaturing effect.


  We study the population genetics of two neutral alleles under reversible
mutation in the \Lambda-processes, a population model that features a skewed
offspring distribution. We describe the shape of the equilibrium allele
frequency distribution as a function of the model parameters. We show that the
mutation rates can be uniquely identified from the equilibrium distribution,
but that the form of the offspring distribution itself cannot be uniquely
identified. We also introduce an infinite-sites version of the \Lambda-process,
and we use it to study how reproductive skew influences standing genetic
diversity in a population. We derive asymptotic formulae for the expected
number of segregating sizes as a function of sample size. We find that the
Wright-Fisher model minimizes the equilibrium genetic diversity, for a given
mutation rate and variance effective population size, compared to all other
\Lambda-processes.


  For societies to produce or safeguard public goods, costly voluntary
contributions are often required. From the perspective of each individual,
however, it is advantageous not to volunteer such contributions, in the hope
that other individuals will carry the associated costs. This conflict can be
modeled as a volunteer's dilemma. To encourage rational individuals to make
voluntary contributions, a government or other social organizations can offer
rewards, to be shared among the volunteers. Here we apply such shared rewarding
to the generalized volunteer's dilemma, in which a threshold number of
volunteers is required for producing the public good. By means of theoretical
and numerical analyses, we show that without shared rewarding only two
evolutionary outcomes are possible: full defection or coexistence of volunteers
and non-volunteers. We show that already small rewards destabilize full
defection, stabilizing small fractions of volunteers instead. Furthermore, at
these intermediate reward levels, we find a hysteresis effect such that
increasing or decreasing group sizes can trigger different social outcomes. In
particular, when group size is increased, the fraction of volunteers first
increases gradually before jumping up abruptly; when group size is then
decreased again, the fraction of volunteers not only remains high, but even
continues to increase. As the shared reward is increased beyond a critical
level, the bi-stability underlying this hysteresis effect vanishes altogether,
and only a single social outcome remains, corresponding to the stable
coexistence of volunteers and non-volunteers. We find that this critical level
of shared rewarding is relatively small compared to the total cost of
contributing to the public good. These results show that the introduction of
shared rewarding is remarkably effective in overcoming defection traps in the
generalized volunteer's dilemma.


  The significant role of epigenetic mechanisms within natural systems has
become increasingly clear. This paper uses a recently presented abstract,
tunable Boolean genetic regulatory network model to explore aspects of
epigenetics. It is shown how dynamically controlling transcription via a DNA
methylation-inspired mechanism can be selected for by simulated evolution under
various single and multiple cell scenarios. Further, it is shown that the
effects of such control can be inherited without detriment to fitness.


  An active colloid is a suspension of particles that transduce free energy
from their environment and use the energy to engage in intrinsically
non-equilibrium activities such as growth, replication and self-propelled
motility. An obvious example of active colloids is a suspension of bacteria
such as Escherichia coli, their physical dimensions being almost invariably in
the colloidal range. Synthetic self-propelled particles have also become
available recently, such as two-faced, or Janus, particles propelled by
differential chemical reactions on their surfaces driving a self-phoretic
motion. In these lectures, I give a pedagogical introduction to the physics of
single-particle and collective properties of active colloids, focussing on self
propulsion. I will compare and contrast phenomena in suspensions of `swimmers'
with the behaviour of suspensions of passive particles, where only Brownian
motion (discovered by Robert Brown in granules from the pollen of the wild
flower {\it Clarkia pulchella}) is relevant. I will pay particular attention to
issues that pertain to performing experiments using these active particle
suspensions, such as how to characterise the suspension's swimming speed
distribution, and include an appendix to guide physicists wanting to start
culturing motile bacteria.


  Vast landscapes extending from southern Angola, Namibia, and South Africa
exhibit localised barren patches of vegetation called fairy circles. They
consist of isolated or randomly distributed circular areas devoid of any
vegetation. They are surrounded by fringes of tall grasses and are embedded in
a sparse grassland on sandy soils. When the aridity increases, the size of the
fairy circles increases and can reach diameters up to 10m. Although several
hypotheses have been proposed, the origin of this phenomenon remains unsolved.
We show that a simple non-local model of plant ecology based on realistic
assumptions provides a quantitative explanation of this phenomenon. Fairy
circles result from strong interaction between interfaces connecting two
uniform covers: the uniform grassland and the bare states, and their
stabilisation is attributed to the Lorentzian-like non-local coupling that
models the competition between plants. The cause of their formation in thus
rather inherent to the vegetation dynamics. Our analysis explain how a circular
shape and fringes are formed, and how the aridity level influences the size of
the fairy circles. In agreement with field observation, these theoretical
findings, provide a link between a strong non-local interaction of plants and
the formation of stable fairy circles. Finally, we show that the proposed
mechanism is model-independent: indeed, it also applies to the
reaction-diffusion type of model that emphasises the influence of water
transport on the vegetation dynamics.


  Progress in Whole Genome Sequencing (WGS) will soon allow a large number of
individuals to have their genome fully sequenced. This lays the foundations to
improve modern healthcare, enabling a new era of personalized medicine where
diagnosis and treatment is tailored to the patient's genetic makeup. It also
allows individuals motivated by personal curiosity to have access to their
genetic information, and use it, e.g., to trace their ancestry. However, the
very same progress also amplifies a number of ethical and privacy concerns,
that stem from the unprecedented sensitivity of genomic information and that
are not well studied. This paper presents an exploratory ethnographic study of
users' perception of privacy and ethical issues with WGS, as well as their
attitude toward different WGS programs. We report on a series of
semi-structured interviews, involving 16 participants, and analyze the results
both quantitatively and qualitatively. Our analysis shows that users exhibit
common trust concerns and fear of discrimination, and demand to retain strict
control over their genetic information. Finally, we highlight the need for
further research in the area and follow-up studies that build on our initial
findings.


  Many foraging animals find food using composite random search strategies,
which consist of intensive and extensive search modes. Models of composite
search can generate predictions about how optimal foragers should behave in
each search mode, and how they should determine when to switch between search
modes. Most of these models assume that foragers use resource encounters to
decide when to switch between search modes. Empirical observations indicate
that a variety of organisms use non-directional sensory cues to identify areas
that warrant intensive search. These cues are not precise enough to allow a
forager to directly orient itself to a resource, but can be used as a criterion
to determine the appropriate search mode. As a potential example, a forager
might use olfactory information, which could help it determine if an area is
worth searching carefully. We developed a model of composite search based on
non-directional sensory cues. With simulations, we compared the search
efficiencies of composite foragers that use resource encounters as their
mode-switching criterion with those that use non-directional sensory cues. We
identified optimal search patterns and mode-switching criteria on a variety of
resource distributions, characterized by different levels of resource
aggregation and density. On all resource distributions, foraging strategies
based on the non-directional sensory criterion were more efficient than those
based on the resource encounter criterion. Strategies based on the
non-directional sensory criterion were also more robust to changes in resource
distribution. Our results suggest that current assumptions about the role of
resource encounters in models of optimal composite search should be
re-examined. The search strategies predicted by our model can help bridge the
gap between random search theory and traditional patch-use foraging theory.


  Conservation of 16 nuclear microsatellite loci, originally developed for
Quercus macrocarpa (section Albae), Q. petraea, Q. robur (section Robur) and Q.
myrsinifolia, (subgenus Cyclobalanopsis) was tested in a Q. infectoria ssp.
veneris population from Cyprus. All loci could be amplified successfully and
displayed allele size and diversity patterns that match those of oak species
belonging to the section Robur. At least in one case, limited amplification and
high levels of homozygosity support the occurrence of 'null alleles', caused by
a possible mutation in the highly conserved primer areas, thus hindering PCR.
The sampled population exhibited high levels of diversity despite the very
limited distribution of this species in Cyprus and extended population
fragmentation. Allele sizes of Q. infectoria at locus QpZAG9 partially match
those of Q. alnifolia and Q. coccifera from neighboring populations. However,
sequencing showed homoplasy, excluding a case of interspecific introgression
with the latter, phylogenetically remote species. Q. infectoria ssp. veneris
sequences at this locus were concordant to those of other species of section
Robur, while sequences of Quercus alnifolia and Quercus coccifera were almost
identical to Q. cerris.


  Indian free ranging dogs live in a carbohydrate rich environment as
scavengers in and around human settlements. They rarely hunt and consequently
do not encounter rich sources of protein. Instead they have adapted to a diet
of primarily carbohydrates. As descendants of the exclusively carnivorous
wolves, they are subjected to the evolutionary load of a physiological demand
for proteins. To meet their protein needs they resort to a thumb rule, if it
smells like meat, eat it. Pups face high competition from group and non group
members and are in a phase of rapid growth with high protein demands. Following
the thumb rule, then they can acquire more protein at the cost of increased
competition and reduced supplementary non protein nutrition. However, if the
mother supplements their diet with protein rich regurgitates and milk, then the
pups can benefit by being generalists. Using a choice test in the field we show
that while adults have a clear preference for meat, pups have no such
preference, and they even eat degraded protein eagerly. Thus the thumb rule
used by adult dogs for efficient scavenging is not innate, and needs to be
learned. The thumb rule might be acquired by cultural transmission, through
exposure to meat in the regurgitate of the mother, or while accompanying her on
foraging trips.