  Molecular phylogeny has focused mainly on improving models for the
reconstruction of gene trees based on sequence alignments. Yet, most
phylogeneticists seek to reveal the history of species. Although the histories
of genes and species are tightly linked, they are seldom identical, because
genes duplicate, are lost or horizontally transferred, and because alleles can
co-exist in populations for periods that may span several speciation events.
Building models describing the relationship between gene and species trees can
thus improve the reconstruction of gene trees when a species tree is known, and
vice-versa. Several approaches have been proposed to solve the problem in one
direction or the other, but in general neither gene trees nor species trees are
known. Only a few studies have attempted to jointly infer gene trees and
species trees. In this article we review the various models that have been used
to describe the relationship between gene trees and species trees. These models
account for gene duplication and loss, transfer or incomplete lineage sorting.
Some of them consider several types of events together, but none exists
currently that considers the full repertoire of processes that generate gene
trees along the species tree. Simulations as well as empirical studies on
genomic data show that combining gene tree-species tree models with models of
sequence evolution improves gene tree reconstruction. In turn, these better
gene trees provide a better basis for studying genome evolution or
reconstructing ancestral chromosomes and ancestral gene sequences. We predict
that gene tree-species tree methods that can deal with genomic data sets will
be instrumental to advancing our understanding of genomic evolution.


  Interacting RNA complexes are studied via bicellular maps using a filtration
via their topological genus. Our main result is a new bijection for RNA-RNA
interaction structures and linear time uniform sampling algorithm for RNA
complexes of fixed topological genus. The bijection allows to either reduce the
topological genus of a bicellular map directly, or to lose connectivity by
decomposing the complex into a pair of single stranded RNA structures. Our main
result is proved bijectively. It provides an explicit algorithm of how to
rewire the corresponding complexes and an unambiguous decomposition grammar.
Using the concept of genus induction, we construct bicellular maps of fixed
topological genus $g$ uniformly in linear time. We present various statistics
on these topological RNA complexes and compare our findings with biological
complexes. Furthermore we show how to construct loop-energy based complexes
using our decomposition grammar.


  The phylogenetic tree space, introduced by Billera, Holmes, and Vogtmann, is
a cone over a simplicial complex. In this short article, we construct this
complex from local gluings of classical polytopes, the associahedron and the
permutohedron. Its homotopy is also reinterpreted and calculated based on
polytope data.


  The molecular responses of macrophages to copper-based nanoparticles have
been investigated via a combination of proteomic and biochemical approaches,
using the RAW264.7 cell line as a model. Both metallic copper and copper oxide
nanoparticles have been tested, with copper ion and zirconium oxide
nanoparticles used as controls. Proteomic analysis highlighted changes in
proteins implicated in oxidative stress responses (superoxide dismutases and
peroxiredoxins), glutathione biosynthesis, the actomyosin cytoskeleton, and
mitochondrial proteins (especially oxidative phosphorylation complex subunits).
Validation studies employing functional analyses showed that the increases in
glutathione biosynthesis and in mitochondrial complexes observed in the
proteomic screen were critical to cell survival upon stress with copper-based
nanoparticles; pharmacological inhibition of these two pathways enhanced cell
vulnerability to copper-based nanoparticles, but not to copper ions.
Furthermore, functional analyses using primary macrophages derived from bone
marrow showed a decrease in reduced glutathione levels, a decrease in the
mitochondrial transmembrane potential, and inhibition of phagocytosis and of
lipopolysaccharide-induced nitric oxide production. However, only a fraction of
these effects could be obtained with copper ions. In conclusion, this study
showed that macrophage functions are significantly altered by copper-based
nanoparticles. Also highlighted are the cellular pathways modulated by cells
for survival and the exemplified cross-toxicities that can occur between
copper-based nanoparticles and pharmacological agents.


  We report a study of a system which involves an enzymatic cascade realizing
an AND logic gate, with an added photochemical processing of the output
allowing to make the gate's response sigmoid in both inputs. New functional
forms are developed for quantifying the kinetics of such systems, specifically
designed to model their response in terms of signal and information processing.
These theoretical expressions are tested for the studied system, which also
allows us to consider aspects of biochemical information processing such as
noise transmission properties and control of timing of the chemical and
physical steps.


  We demonstrate a vast expansion of the theory of evolutionary stability to
finite populations with mutation, connecting the theory of the stationary
distribution of the Moran process with the Lyapunov theory of evolutionary
stability. We define the notion of stationary stability for the Moran process
with mutation and generalizations, as well as a generalized notion of
evolutionary stability that includes mutation called an incentive stable state
(ISS) candidate. For sufficiently large populations, extrema of the stationary
distribution are ISS candidates and we give a family of Lyapunov quantities
that are locally minimized at the stationary extrema and at ISS candidates. In
various examples, including for the Moran and Wright-Fisher processes, we show
that the local maxima of the stationary distribution capture the
traditionally-defined evolutionarily stable states. The classical stability
theory of the replicator dynamic is recovered in the large population limit.
Finally we include descriptions of possible extensions to populations of
variable size and populations evolving on graphs.


  Analysis of genomic segments shared identical-by-descent (IBD) between
individuals is fundamental to many genetic applications, from demographic
inference to estimating the heritability of diseases, but IBD detection
accuracy in non-simulated data is largely unknown. In principle, it can be
evaluated using known pedigrees, as IBD segments are by definition inherited
without recombination down a family tree. We extracted 25,432 genotyped
European individuals containing 2,952 father-mother-child trios from the
23andMe, Inc. dataset. We then used GERMLINE, a widely used IBD detection
method, to detect IBD segments within this cohort. Exploiting known familial
relationships, we identified a false positive rate over 67% for 2-4 centiMorgan
(cM) segments, in sharp contrast with accuracies reported in simulated data at
these sizes. Nearly all false positives arose from the allowance of haplotype
switch errors when detecting IBD, a necessity for retrieving long (> 6 cM)
segments in the presence of imperfect phasing. We introduce HaploScore, a
novel, computationally efficient metric that scores IBD segments proportional
to the number of switch errors they contain. Applying HaploScore filtering to
the IBD data at a precision of 0.8 produced a 13-fold increase in recall when
compared to length-based filtering. We replicate the false IBD findings and
demonstrate the generalizability of HaploScore to alternative data sources
using an independent cohort of 555 European individuals from the 1000 Genomes
project. HaploScore can improve the accuracy of segments reported by any IBD
detection method, provided that estimates of the genotyping error rate and
switch error rate are available.


  Nitrogen (N) removal from estuaries is driven in part by sedimentary
microbial processes. The processes of denitrification and anaerobic ammonium
oxidation (anammox) remove N from estuaries by producing di-nitrogen gas, and
each can be coupled to N recycling pathways such as nitrification and
dissimilatory nitrate reduction to ammonium (DNRA). Environmental conditions in
estuaries influence sedimentary N cycling processes; therefore, seawater
intrusion may affect the coupling of N cycling processes in the freshwater
portions of estuaries. This study investigated the potential effects of
seawater intrusion on these process couplings through a comparative modeling
approach. We applied environ analysis, a form of ecological network analysis,
to two N cycling mass-balance network models constructed at freshwater
(oligohaline) and saltwater (polyhaline) sites in the Cape Fear River Estuary,
North Carolina. We used a space-for-time substitution to predict the effects of
seawater intrusion on the sedimentary N cycle. Further, we conducted an
uncertainty analysis using linear inverse modeling to evaluate the effects of
parameterization uncertainty on model results. Nitrification coupled to both
denitrification and anammox was 2.5 times greater in the oligohaline model,
while DNRA coupled to anammox was 2.7 times greater in the polyhaline model.
However, the total amount of N2 gas produced relative to the nitrogen inputs to
each network was 4.7% and 4.6% at the oligohaline and polyhaline sites,
respectively. These findings suggest that changes in water chemistry from
seawater intrusion may favor direct over coupled nitrogen removal, but may not
substantially change the N removal capacity of the sedimentary microbial
processes.


  We introduce a pipeline to easily generate collections of web accessible UCSC
genome browsers interrelated by an alignment. Using the alignment, all
annotations and the alignment itself can be efficiently viewed with reference
to any genome in the collection, symmetrically. A new, intelligently scaled
alignment display makes it simple to view all changes between the genomes at
all levels of resolution, from substitutions to complex structural
rearrangements, including duplications.


  One of the most challenging and long-standing problems in computational
biology is the prediction of three-dimensional protein structure from amino
acid sequence. A promising approach to infer spatial proximity between residues
is the study of evolutionary covariance from multiple sequence alignments,
especially in light of recent algorithmic improvements and the fast growing
size of sequence databases. In this paper, we present a simple, fast and
accurate algorithm for the prediction of residue-residue contacts based on
regularized least squares. The basic assumption is that spatially proximal
residues in a protein coevolve to maintain the physicochemical complementarity
of the amino acids involved in the contact. Our regularized inversion of the
sample covariance matrix allows the computation of partial correlations between
pairs of residues, thereby removing the effect of spurious transitive
correlations. The method also accounts for low number of observations by means
of a regularization parameter that depends on the effective number of sequences
in the alignment. When tested on a set of protein families from Pfam, we found
the RLS algorithm to have performance comparable to state-of-the-art methods
for contact prediction, while at the same time being faster and conceptually
simpler.


  A Flory theory is constructed for a long polymer ring in a melt of unknotted
and non-concatenated rings. The theory assumes that the ring forms an effective
annealed branched object and computes its primitive path. It is shown that the
primitive path follows self-avoiding statistics and is characterized by the
corresponding Flory exponent of a polymer with excluded volume. Based on that,
it is shown that rings in the melt are compact objects with overall size
proportional to their length raised to the 1/3 power. Furthermore, the contact
probability exponent is estimated, albeit by a poorly controlled approximation,
with the result consistent with both numerical and experimental data.


  The appearance of a novel coronavirus named Middle East (ME) Respiratory
Syndrome Coronavirus (MERS-CoV) has raised global public health concerns
regarding the current situation and its future evolution. Here we propose an
integrative maximum likelihood analysis of both cluster data in the ME region
and importations in Europe to assess transmission scenario and incidence of
sporadic infections. Our approach is based on a spatial-transmission model
integrating mobility data worldwide and allows for variations in the
zoonotic/environmental transmission and underascertainment. Maximum likelihood
estimates for the ME region indicate the occurrence of a subcritical epidemic
(R=0.50, 95% confidence interval (CI) 0.30-0.77) associated with a 0.28 (95% CI
0.12-0.85) daily rate of sporadic introductions. Infections in the region
appear to be mainly dominated by zoonotic/environmental transmissions, with
possible underascertainment (95% CI of estimated to observed sporadic cases in
the range 1.03-7.32). No time evolution of the situation emerges. Analyses of
flight passenger data from the region indicate areas at high risk of
importation. While dismissing an immediate threat for global health security,
this analysis provides a baseline scenario for future reference and updates,
suggests reinforced surveillance to limit underascertainment, and calls for
increased alertness in high-risk areas worldwide.


  Despite the general acknowledgment of the role of niche and fitness
differences in community dynamics, species abundance has been coined as a
relevant feature not just regarding niche perspectives, but also according to
neutral perspectives. Here we explore a minimum probabilistic stochastic model
to evaluate the role of populations relative initial and total abundances on
species chances to outcompete each other and their persistence in time (i.e.,
unstable coexistence). The present results show that taking into account the
stochasticity in demographic properties and conservation of individuals in
closed communities (zero-sum assumption), population initial abundance can
strongly influence species chances to outcompete each other, and also influence
the period of coexistence of these species in a particular time interval.
Systems carrying capacity can have an important role in species coexistence by
exacerbating fitness inequalities and affecting the size of the period of
coexistence. This study also shows that populations initial abundances can act
as an equalizing mechanism, reducing fitness inequalities, which can favor
species coexistence and even make less fitted species to be more likely to
outcompete better fitted species, and thus to dominate ecological communities
in the absence of niche mechanisms.


  Myoglobin (Mb) binds diatomic ligands, like O$_2$, CO, and NO, in a cavity
that is only transiently accessible. Crystallography and molecular simulations
show that the ligands can migrate through an extensive network of transiently
connected cavities, but disagree on the locations and occupancy of internal
hydration sites. Here, we use water $^2$H and $^{17}$O magnetic relaxation
dispersion (MRD) to characterize the internal water molecules in Mb under
physiological conditions. We find that equine carbonmonoxy Mb contains 4.5
$\pm$ 1.0 ordered internal water molecules with a mean survival time of 5.6
$\pm$ 0.5 $\mu$s at 25 $^\circ$C. The likely location of these water molecules
are the four polar hydration sites, including one of the xenon-binding
cavities, that are fully occupied in all high-resolution crystal structures of
equine Mb. The finding that water escapes from these sites, located 17 -- 31
{\AA} apart in the protein, on the same $\mu$s time scale suggests a global
exchange mechanism. We propose that this mechanism involves transient
penetration of the protein by H-bonded water chains. Such a mechanism could
play a functional role by eliminating trapped ligands. In addition, the MRD
results indicate that two or three of the 11 histidine residues of equine Mb
undergo intramolecular hydrogen exchange on a $\mu$s time scale.


  Network ecologists investigate the structure, function, and evolution of
ecological systems using network models and analyses. For example, network
techniques have been used to study community interactions (i.e., food-webs,
mutualisms), gene flow across landscapes, and the sociality of individuals in
populations. The work presented here uses a bibliographic and network approach
to (1) document the rise of Network Ecology, (2) identify the diversity of
topics addressed in the field, and (3) map the structure of scientific
collaboration among contributing scientists. Our aim is to provide a broad
overview of this emergent field that highlights its diversity and to provide a
foundation for future advances. To do this, we searched the ISI Web of Science
database for ecology publications between 1900 and 2012 using the search terms
for research areas of Environmental Sciences & Ecology and Evolutionary Biology
and the topic tag ecology. From these records we identified the Network Ecology
publications using the topic terms network, graph theory, and web while
controlling for the usage of misleading phrases. The resulting corpus entailed
29,513 publications between 1936 and 2012. We document the rapid rise in
Network Ecology publications per year reaching a magnitude of over 5% of the
ecological publications in 2012. Drawing topical information from the
publication record content (titles, abstracts, keywords) and collaboration
information from author listing, our analysis highlights the diversity and
clustering of topics addressed within Network Ecology and reveals the highly
collaborative approach of scientists publishing in this field. We conclude that
Network Ecology is a large and rapidly growing area of ecology, and we expect
continued growth of this research field.


  Analysis of demographic catastrophes shows that, with the exception of
perhaps only two critical events, they were too weak to influence the growth of
human population. These results reinforce the conclusion that the concept of
the Epoch of Malthusian Stagnation, the alleged first stage of growth claimed
by the Demographic Transition Theory, is not supported by empirical evidence.
They show that even if we assume that Malthusian positive checks are capable of
suppressing the growth of population their impact was too weak to create the
Epoch of Malthusian Stagnation.


  We present a novel model that describes the within-host evolutionary dynamics
of parasites undergoing antigenic variation. The approach uses a multi-type
branching process with two types of entities defined according to their
relationship with the immune system: clans of resistant parasitic cells (i.e.
groups of cells sharing the same antigen not yet recognized by the immune
system) that may become sensitive, and individual sensitive cells that can
acquire a new resistance thus giving rise to the emergence of a new clan. The
simplicity of the model allows analytical treatment to determine the
subcritical and supercritical regimes in the space of parameters. By
incorporating a density-dependent mechanism the model is able to capture
additional relevant features observed in experimental data, such as the
characteristic parasitemia waves. In summary our approach provides a new
general framework to address the dynamics of antigenic variation which can be
easily adapted to cope with broader and more complex situations.


  Analyses and visualizations by the ISSCOR method of influenza virus
hemagglutinin genes of different A-subtypes revealed some rather striking
temporal relationships between groups of individual gene subsets. Based on
these findings we consider application of the ISSCOR-PCA method for analyses of
large sets of homologous genes to be a worthwhile addition to a toolbox of
genomics - allowing for a rapid diagnostics of trends, and ultimately even
aiding an early warning of newly emerging epidemiological threats.


  Phylogenetic analyses based on small to moderately sized sets of sequential
data lead to overestimating mutation rates in influenza hemagglutinin (HA) by
at least an order of magnitude. Two major underlying reasons are: the
incomplete lineage sorting, and a possible absence in the analyzed sequences
set some of key missing ancestors. Additionally, during neighbor joining tree
reconstruction each mutation is considered equally important, regardless of its
nature. Here we have implemented a heuristic method optimizing site dependent
factors weighting differently 1st, 2nd, and 3rd codon position mutations,
allowing to extricate incorrectly attributed sub-clades. The least squares
regression analysis of distribution of frequencies for all mutations observed
on a partially disentangled tree for a large set of unique 3243 HA sequences,
along all nucleotide positions, was performed for all mutations as well as for
non-equivalent amino acid mutations: in both cases demonstrating almost flat
gradients, with a very slight downward slope towards the 3'-end positions. The
mean mutation rates per sequence per year were 3.83*10^-4 for the all
mutations, and 9.64*10^-5 for the non-equivalent ones.


  Gross chromosomal rearrangements have the potential to be evolutionarily
advantageous to an adapting organism. The generation of a hybrid species
increases opportunity for recombination by bringing together two homologous
genomes. We sought to define the location of genomic rearrangements in three
strains of Saccharomyces pastorianus, a natural lager-brewing yeast hybrid of
Saccharomyces cerevisiae and Saccharomyces eubayanus, using whole genome
shotgun sequencing. Each strain of S. pastorianus has lost species-specific
portions of its genome and has undergone extensive recombination, producing
chimeric chromosomes. We predicted 30 breakpoints that we confirmed at the
single nucleotide level by designing species-specific primers that flank each
breakpoint, and then sequencing the PCR product. These rearrangements are the
result of recombination between areas of homology between the two subgenomes,
rather than repetitive elements such as transposons or tRNAs. Interestingly,
28/30 S. cerevisiae- S. eubayanus recombination breakpoints are located within
genic regions, generating chimeric genes. Furthermore we show evidence for the
reuse of two breakpoints, located in HSP82 and KEM1, in strains of proposed
independent origin.


  Animals foraging alone are hypothesized to optimize the encounter rates with
resources through L\'evy walks. However, the issue of how the interactions
between multiple foragers influence their search efficiency is still not
completely understood. To address this, we consider a model to study the
optimal strategy for a group of foragers searching for targets distributed
heterogeneously. In our model foragers move on a square lattice containing
immobile but regenerative targets. At any instant a forager is able to detect
only those targets that happen to be in the same site. However, we allow the
foragers to have information about the state of other foragers. A forager who
has not detected any target walks towards the nearest location, where another
forager has detected a target, with probability $\exp{\left(-\alpha d\right)}$,
where $d$ is the distance and $\alpha$ is a parameter. The model reveals that
neither overcrowding ($\alpha\to 0$) nor independent searching
($\alpha\to\infty$) is beneficial for the group. For patchy distribution of
targets the efficiency is maximum for intermediate values of $\alpha$. Also, in
the limit $\alpha\to 0$, the length of the walks can become scale-free.


  Theoretical modeling of computer virus/worm epidemic dynamics is an important
problem that has attracted many studies. However, most existing models are
adapted from biological epidemic ones. Although biological epidemic models can
certainly be adapted to capture some computer virus spreading scenarios
(especially when the so-called homogeneity assumption holds), the problem of
computer virus spreading is not well understood because it has many important
perspectives that are not necessarily accommodated in the biological epidemic
models. In this paper we initiate the study of such a perspective, namely that
of adaptive defense against epidemic spreading in arbitrary networks. More
specifically, we investigate a non-homogeneous
Susceptible-Infectious-Susceptible (SIS) model where the model parameters may
vary with respect to time. In particular, we focus on two scenarios we call
semi-adaptive defense and fully-adaptive} defense, which accommodate implicit
and explicit dependency relationships between the model parameters,
respectively. In the semi-adaptive defense scenario, the model's input
parameters are given; the defense is semi-adaptive because the adjustment is
implicitly dependent upon the outcome of virus spreading. For this scenario, we
present a set of sufficient conditions (some are more general or succinct than
others) under which the virus spreading will die out; such sufficient
conditions are also known as epidemic thresholds in the literature. In the
fully-adaptive defense scenario, some input parameters are not known (i.e., the
aforementioned sufficient conditions are not applicable) but the defender can
observe the outcome of virus spreading. For this scenario, we present adaptive
control strategies under which the virus spreading will die out or will be
contained to a desired level.


  Fluctuating environments pose tremendous challenges to bacterial populations.
It is observed in numerous bacterial species that individual cells can
stochastically switch among multiple phenotypes for the population to survive
in rapidly changing environments. This kind of phenotypic heterogeneity with
stochastic phenotype switching is generally understood to be an adaptive
bet-hedging strategy. Mathematical models are essential to gain a deeper
insight into the principle behind bet-hedging and the pattern behind
experimental data. Traditional deterministic models cannot provide a correct
description of stochastic phenotype switching and bet-hedging, and traditional
Markov chain models at the cellular level fail to explain their underlying
molecular mechanisms. In this paper, we propose a nonlinear stochastic model of
multistable bacterial systems at the molecular level. It turns out that our
model not only provides a clear description of stochastic phenotype switching
and bet-hedging within isogenic bacterial populations, but also provides a
deeper insight into the analysis of multidimensional experimental data.
Moreover, we use some deep mathematical theories to show that our stochastic
model and traditional Markov chain models are essentially consistent and
reflect the dynamic behavior of the bacterial system at two different time
scales. In addition, we provide a quantitative characterization of the critical
state of multistable bacterial systems and develop an effective data-driven
method to identify the critical state without resorting to specific
mathematical models.


  Tracking the dynamics of fluorescent nanoparticles during embryonic
development allows insights into the physical state of the embryo and,
potentially, molecular processes governing developmental mechanisms. In this
work, we investigate the motion of individual fluorescent nanodiamonds
micro-injected into Drosophila melanogaster embryos prior to cellularisation.
Fluorescence correlation spectroscopy and wide-field imaging techniques are
applied to individual fluorescent nanodiamonds in blastoderm cells during stage
5 of development to a depth of ~40 \mu m. The majority of nanodiamonds in the
blastoderm cells during cellularisation exhibit free diffusion with an average
diffusion coefficient of (6 $\pm$ 3) x 10$^{-3}$ \mu m$^2$/s, (mean $\pm$ SD).
Driven motion in the blastoderm cells was also observed with an average
velocity of 0.13 $\pm$ 0.10 \mu m/s (mean $\pm$ SD) \mu m/s and an average
applied force of 0.07 $\pm$ 0.05 pN (mean $\pm$ SD). Nanodiamonds in the
periplasm between the nuclei and yolk were also found to undergo free diffusion
with a significantly larger diffusion coefficient of (63 $\pm$ 35) x10$^{-3}$
\mu m$^2$/s (mean $\pm$ SD). Driven motion in this region exhibited similar
average velocities and applied forces compared to the blastoderm cells
indicating the transport dynamics in the two cytoplasmic regions are analogous.


  Cancer is a disease involving many genes, consequently it has been difficult
to design anticancer drugs that are efficacious over a broad range of cancers.
The robustness of cellular responses to gene knockout and the need to reduce
undesirable side effects also contribute to the problem of effective
anti-cancer drug design. To promote the successful selection of drug targets,
each potential target should be subjected to a systems biology scrutiny to
locate effective and specific targets while minimizing undesirable side
effects. The gene p53 is considered a good candidate for such a target because
it has been implicated in 50% of all cancers and is considered to encode a hub
protein that is highly connected to other proteins. Using P53 as a test case,
this paper explores the capacity of the systems biology tool, CellDesigner, to
aid in the selection of anticancer drug targets and to serve as a teaching
resource for human resource development.


  The accumulation of adaptive mutations is essential for survival in novel
environments. However, in clonal populations with a high mutational supply, the
power of natural selection is expected to be limited. This is due to clonal
interference - the competition of clones carrying different beneficial
mutations - which leads to the loss of many small effect mutations and fixation
of large effect ones. If interference is abundant, then mechanisms for
horizontal transfer of genes, which allow the immediate combination of
beneficial alleles in a single background, are expected to evolve. However, the
relevance of interference in natural complex environments, such as the gut, is
poorly known. To address this issue, we studied the invasion of beneficial
mutations responsible for Escherichia coli's adaptation to the mouse gut and
demonstrate the pervasiveness of clonal interference. The observed dynamics of
change in frequency of beneficial mutations are consistent with soft sweeps,
where a similar adaptive mutation arises repeatedly on different haplotypes
without reaching fixation. The genetic basis of the adaptive mutations revealed
a striking parallelism in independently evolving populations. This was mainly
characterized by the insertion of transposable elements in both coding and
regulatory regions of a few genes. Interestingly in most populations, we
observed a complete phenotypic sweep without loss of genetic variation. The
intense clonal interference during adaptation to the gut environment, here
demonstrated, may be important for our understanding of the levels of strain
diversity of E. coli inhabiting the human gut microbiota and of its
recombination rate.


  Biological network alignment aims to identify similar regions between
networks of different species. Existing methods compute node "similarities" to
rapidly identify from possible alignments the "high-scoring" alignments with
respect to the overall node similarity. However, the accuracy of the alignments
is then evaluated with some other measure that is different than the node
similarity used to construct the alignments. Typically, one measures the amount
of conserved edges. Thus, the existing methods align similar nodes between
networks hoping to conserve many edges (after the alignment is constructed!).
  Instead, we introduce MAGNA to directly "optimize" edge conservation while
the alignment is constructed. MAGNA uses a genetic algorithm and our novel
function for crossover of two "parent" alignments into a superior "child"
alignment to simulate a "population" of alignments that "evolves" over time;
the "fittest" alignments survive and proceed to the next "generation", until
the alignment accuracy cannot be optimized further. While we optimize our new
and superior measure of the amount of conserved edges, MAGNA can optimize any
alignment accuracy measure. In systematic evaluations against existing
state-of-the-art methods (IsoRank, MI-GRAAL, and GHOST), MAGNA improves
alignment accuracy of all methods.


  The reference annotations made for a genome sequence provide the framework
for all subsequent analyses of the genome. Correct annotation is particularly
important when interpreting the results of RNA-seq experiments where short
sequence reads are mapped against the genome and assigned to genes according to
the annotation. Inconsistencies in annotations between the reference and the
experimental system can lead to incorrect interpretation of the effect on RNA
expression of an experimental treatment or mutation in the system under study.
Until recently, the genome-wide annotation of 3-prime untranslated regions
received less attention than coding regions and the delineation of intron/exon
boundaries. In this paper, data produced for samples in Human, Chicken and A.
thaliana by the novel single-molecule, strand-specific, Direct RNA Sequencing
technology from Helicos Biosciences which locates 3-prime polyadenylation sites
to within +/- 2 nt, were combined with archival EST and RNA-Seq data. Nine
examples are illustrated where this combination of data allowed: (1) gene and
3-prime UTR re-annotation (including extension of one 3-prime UTR by 5.9 kb);
(2) disentangling of gene expression in complex regions; (3) clearer
interpretation of small RNA expression and (4) identification of novel genes.
While the specific examples displayed here may become obsolete as genome
sequences and their annotations are refined, the principles laid out in this
paper will be of general use both to those annotating genomes and those seeking
to interpret existing publically available annotations in the context of their
own experimental data


  Due to recent advances in synthetic biology and artificial life, the origin
of life is currently a hot topic of research. We review the literature and
argue that the two traditionally competing "replicator-first" and
"metabolism-first" approaches are merging into one integrated theory of
individuation and evolution. We contribute to the maturation of this more
inclusive approach by highlighting some problematic assumptions that still lead
to an impoverished conception of the phenomenon of life. In particular, we
argue that the new consensus has so far failed to consider the relevance of
intermediate timescales. We propose that an adequate theory of life must
account for the fact that all living beings are situated in at least four
distinct timescales, which are typically associated with metabolism, motility,
development, and evolution. On this view, self-movement, adaptive behavior and
morphological changes could have already been present at the origin of life. In
order to illustrate this possibility we analyze a minimal model of life-like
phenomena, namely of precarious, individuated, dissipative structures that can
be found in simple reaction-diffusion systems. Based on our analysis we suggest
that processes in intermediate timescales could have already been operative in
prebiotic systems. They may have facilitated and constrained changes occurring
in the faster- and slower-paced timescales of chemical self-individuation and
evolution by natural selection, respectively.


  Nucleic acids, like DNA and RNA, are molecules that are present in any life
form. Their most notable function is to encode biological information. Why then
would a physicist be interested in these molecules? As we will see, DNA is an
interesting molecular tool for physicists to test and explore physical laws and
theories, like the ergodic theorem, the theory of elasticity and information
theory. DNA also has unique material properties, which attract material
scientists, nanotechnologists and engineers. Among interesting developments in
this field are DNA-based hybrid materials and DNA origami.


  Mutational robustness is the extent to which an organism has evolved to
withstand the effects of deleterious mutations. We explored the extent of
mutational robustness in the budding yeast by genome wide dosage suppressor
analysis of 53 conditional lethal mutations in cell division cycle and RNA
synthesis related genes, revealing 660 suppressor interactions of which 642 are
novel. This collection has several distinctive features, including high
co-occurrence of mutant-suppressor pairs within protein modules, highly
correlated functions between the pairs, and higher diversity of functions among
the co-suppressors than previously observed. Dosage suppression of essential
genes encoding RNA polymerase subunits and chromosome cohesion complex suggest
a surprising degree of functional plasticity of macromolecular complexes and
the existence od degenerate pathways for circumventing potentially lethal
mutations. The utility of dosage-suppressor networks is illustrated by the
discovery of a novel connection between chromosome cohesion-condensation
pathways involving homologous recombination, and Huntington's disease.


  Neuronal membrane potentials fluctuate stochastically due to conductance
changes caused by random transitions between the open and close states of ion
channels. Although it has previously been shown that channel noise can
nontrivially affect neuronal dynamics, it is unknown whether ion-channel noise
is strong enough to act as a noise source for hypothesised noise-enhanced
information processing in real neuronal systems, i.e. 'stochastic
facilitation.' Here, we demonstrate that biophysical models of channel noise
can give rise to two kinds of recently discovered stochastic facilitation
effects in a Hodgkin-Huxley-like model of auditory brainstem neurons. The
first, known as slope-based stochastic resonance (SBSR), enables phasic neurons
to emit action potentials that can encode the slope of inputs that vary slowly
relative to key time-constants in the model. The second, known as inverse
stochastic resonance (ISR), occurs in tonically firing neurons when small
levels of noise inhibit tonic firing and replace it with burst-like dynamics.
Consistent with previous work, we conclude that channel noise can provide
significant variability in firing dynamics, even for large numbers of channels.
Moreover, our results show that possible associated computational benefits may
occur due to channel noise in neurons of the auditory brainstem. This holds
whether the firing dynamics in the model are phasic (SBSR can occur due to
channel noise) or tonic (ISR can occur due to channel noise).


  Cooperation is ubiquitous in nature, but explaining its existence remains a
central interdisciplinary challenge. Cooperation is most difficult to explain
in the Prisoner's Dilemma game, where cooperators always lose in direct
competition with defectors despite increasing mean fitness. Here we demonstrate
how spatial population expansion, a widespread natural phenomenon, promotes the
evolution of cooperation. We engineer an experimental Prisoner's Dilemma game
in the budding yeast Saccharomyces cerevisiae to show that, despite losing to
defectors in nonexpanding conditions, cooperators increase in frequency in
spatially expanding populations. Fluorescently labeled colonies show genetic
demixing of cooperators and defectors, followed by increase in cooperator
frequency as cooperator sectors overtake neighboring defector sectors. Together
with lattice-based spatial simulations, our results suggest that spatial
population expansion drives the evolution of cooperation by (1) increasing
positive genetic assortment at population frontiers and (2) selecting for
phenotypes maximizing local deme productivity. Spatial expansion thus creates a
selective force whereby cooperator-enriched demes overtake neighboring
defector-enriched demes in a "survival of the fastest". We conclude that colony
growth alone can promote cooperation and prevent defection in microbes. Our
results extend to other species with spatially restricted dispersal undergoing
range expansion, including pathogens, invasive species, and humans.


  The inositol trisphosphate receptor (IPR) is a crucial ion channel that
regulates the Ca$^{2+}$ influx from the endoplasmic reticulum (ER) to the
cytoplasm. A thorough study of the IPR channel contributes to a better
understanding of calcium oscillations and waves. It has long been observed that
the IPR channel is a typical biological system which performs adaptation.
However, recent advances on the physical essence of adaptation show that
adaptation systems with a negative feedback mechanism, such as the IPR channel,
must break detailed balance and always operate out of equilibrium with energy
dissipation. Almost all previous IPR models are equilibrium models assuming
detailed balance and thus violate the physical essence of adaptation. In this
article, we constructed a nonequilibrium allosteric model of single IPR
channels based on the patch-clamp experimental data obtained from the IPR in
the outer membranes of isolated nuclei of the \emph{Xenopus} oocyte. It turns
out that our model reproduces the patch-clamp experimental data reasonably well
and produces both the correct steady-state and dynamic properties of the
channel. Particularly, our model successfully describes the complicated bimodal
[Ca$^{2+}$] dependence of the mean open duration at high [IP$_3$], a
steady-state behavior which fails to be correctly described in previous IPR
models. Finally, we used the patch-clamp experimental data to validate that the
IPR channel indeed breaks detailed balance and thus is a nonequilibrium system
which consumes energy.


  With the assumption that the effect is a mathematical function of the cause
in a causal relationship, FunChisq, a chi-square test defined on a
non-parametric representation of interactions, infers network topology
considering both interaction directionality and nonlinearity. Here we show that
both experimental and in silico biological network data suggest the importance
of directionality as evidence for causality. Counter-intuitively, patterns in
those interactions effectively revealed by FunChisq enlist an experimental
design principle essential to network inference -- perturbations to a
biological system shall make it transits between linear and nonlinear working
zones, instead of operating only in a linear working zone.


  DNA copy number aberrations (CNAs) can be found in the majority of cancer
genomes, and are crucial for understanding the potential mechanisms underlying
tumor initiation and progression. Since the first release in 2001, the
Progenetix project (http://www.progenetix.org) has provided a reference
resource dedicated to provide the most comprehensive collection of genome-wide
CNA profiles. Reflecting the application of comparative genomic hybridization
(CGH) techniques to tens of thousands of cancer genomes, over the past 12 years
our data curation efforts have resulted in a more than 60-fold increase in the
number of cancer samples presented through Progenetix. In addition, new data
exploration tools and visualization options have been added. In particular, the
gene specific CNA frequency analysis should facilitate the assignment of cancer
genes to related cancer types. Additionally, the new user file processing
interface allows users to take advantage of the online tools, including various
data representation options for proprietary data pre-publication. In this
update article, we report recent improvements of the database in terms of
content, user interface and online tools.


  One of the main challenges for biomedical research lies in the
computer-assisted integrative study of large and increasingly complex
combinations of data in order to understand molecular mechanisms. The
preservation of the materials and methods of such computational experiments
with clear annotations is essential for understanding an experiment, and this
is increasingly recognized in the bioinformatics community. Our assumption is
that offering means of digital, structured aggregation and annotation of the
objects of an experiment will provide necessary meta-data for a scientist to
understand and recreate the results of an experiment. To support this we
explored a model for the semantic description of a workflow-centric Research
Object (RO), where an RO is defined as a resource that aggregates other
resources, e.g., datasets, software, spreadsheets, text, etc. We applied this
model to a case study where we analysed human metabolite variation by
workflows.


  Current theories about mechanisms promoting species co-existence in diverse
communities assume that species only interact ecologically. Species are treated
as discrete evolutionary entities, even though abundant empirical evidence
indicates that patterns of gene flow such as selfing and hybridization
frequently occur in plant and animal groups. Here, we allow mating behavior to
respond to local species composition and abundance in a data-driven
meta-community model of species co-existence. While individuals primarily
out-cross, they also maintain some small capacity for selfing and
hybridization. Mating choice is treated as a 'fuzzy' behavior, determined by an
interaction between intrinsic properties affecting mate choice and extrinsic
properties of the local community, primarily the density and availability of
sympatric inter-fertile species. When mate choice is strongly limited, even low
survivorship of selfed offspring (<10%) can prevent extinction of rare species.
With increasing mate choice, low hybridization success rates (~20%) maintain
community level diversity for substantially longer periods of time. Given the
low species densities and high diversity of tropical tree communities, the
evolutionary costs of competition among sympatric congeneric species are
negligible because direct competition is infrequent. In diverse communities,
many species are chronically rare and thus vulnerable to stochastic extinction,
which occurs rapidly if species are completely reproductively isolated. By
incorporating fuzzy mating behavior into models of species co-existence, a more
realistic understanding of the extinction process can be developed. Fuzzy
mating strategies, potentially an important mechanism for rare species to
escape extinction and gain local adaptations, should be incorporated into
forest management strategies.


  A malaria model is formulated which includes the enhanced attractiveness of
infectious humans to mosquitoes, as result of host manipulation by malaria
parasite, and the human behavior, represented by insecticide-treated bed nets
usage. The occurrence of a backward bifurcation at R_0=1 is shown to be
possible, which implies that multiple endemic equilibria co-exist with a stable
disease-free equilibrium when the basic reproduction number is less than unity.
This phenomenon is found to be caused by disease-induced human mortality. The
global asymptotic stability of the endemic equilibrium for R_0>1 is proved, by
using the geometric approach to global stability. Therefore, the disease
becomes endemic for $R_0>1$ regardless of the number of initial cases in both
the human and vector populations. Finally, the impact of vector's host
preferences and bed-net usage behavior on system dynamics is investigated.


  Understanding the influence of an environment on the evolution of its
resident population is a major challenge in evolutionary biology. Great
progress has been made in homogeneous population structures while heterogeneous
structures have received relatively less attention. Here we present a
structured population model where different individuals are best suited to
different regions of their environment. The underlying structure is a graph:
individuals occupy vertices, which are connected by edges. If an individual is
suited for their vertex, they receive an increase in fecundity. This framework
allows attention to be restricted to the spatial arrangement of suitable
habitat. We prove some basic properties of this model and find some
counter-intuitive results. Notably, 1) the arrangement of suitable sites is as
important as their proportion, and, 2) decreasing the proportion of suitable
sites may result in a decrease in the fixation time of an allele.


  Knowledge of the exact distribution of meiotic crossovers (COs) and gene
conversions (GCs) is essential for understanding many aspects of population
genetics and evolution, from haplotype structure and long-distance genetic
linkage to the generation of new allelic variants of genes. To this end, we
resequenced the four products of 13 meiotic tetrads along with 10 doubled
haploids derived from Arabidopsis thaliana hybrids. GC detection through short
reads has previously been confounded by genomic rearrangements. Rigid filtering
for misaligned reads allowed GC identification at high accuracy and revealed an
~80-kb transposition, which undergoes copy-number changes mediated by meiotic
recombination. Non-crossover associated GCs were extremely rare most likely due
to their short average length of ~25-50 bp, which is significantly shorter than
the length of CO associated GCs. Overall, recombination preferentially targeted
non-methylated nucleosome-free regions at gene promoters, which showed
significant enrichment of two sequence motifs.


  The phenomenon of oral tolerance refers to a local and systemic state of
tolerance, induced in the gut associated lymphoid tissues, after its exposure
to innocuous antigens, such as food proteins. While recent findings shed light
in the cellular and molecular basis of oral tolerance, the network of
interactions between the components mediating oral tolerance has not been
investigated yet. Our work brings a complex systems theory approach, aiming to
identify the contribution of each element in an oral tolerance network. We also
propose a model that allows dynamical plus topological quantifying which must
encompass functional responses as the local host involved on the oral
tolerance. To keep track of reality of our model, we test knockout (KO) of
immunological components (i. e. silencing a vertex) and see how it diverges
when the system is topologically health. The results from these simulated KO's
are then compared to real molecular knock-outs. To infer from these processing
we apply a new implementation of a random walk algorithm for directed graphs,
which ultimately generate statistical quantities provided by the dynamical
behaviour of the simulated KO's. It was observed that the a specifics KO caused
the greatest impact on network standard flux. In a brief analysis, the results
obtained correspond to biological data. Our model addresses both topological
proprieties and dynamical relations. The construction of a qualitative dynamic
model for oral tolerance could reflect empirical observations, through the
standard flux results and relative error based on individual knockout.


  The size of mature oocytes is similar across mammalian species, yet the size
of ovarian follicles increases with species size, with some ovarian follicles
reaching diameters more than 1000-fold the size of the enclosed oocyte. Here we
show that the different follicular sizes can be explained with diffusion-based
limitations on the thickness of the hormone-secreting granulosa layer. By
analysing published data on human follicular growth and granulosa cell
expansion during follicular maturation we find that the 4-fold increase of the
antral follicle diameter is entirely driven by an increase in the follicular
fluid volume, while the thickness of the surrounding granulosa layer remains
constant at about 45+/-10 mkm. Based on the measured kinetic constants, the
model reveals that the observed fall in the gonadotropin concentration from
peripheral blood circulation to the follicular antrum is a result of
sequestration in the granulosa. The model further shows that as a result of
sequestration, an increased granulosa thickness cannot substantially increase
estradiol production but rather deprives the oocyte from gonadotropins. Larger
animals (with a larger blood volume) require more estradiol as produced by the
ovaries to downregulate FSH-secretion in the pituitary. Larger follicle
diameters result in larger follicle surface areas for constant granulosa layer
thickness. The reported increase in follicular surface area in larger species
indeed correlates linearly both with species mass and with the predicted
increase in estradiol output. In summary, we propose a structural role for the
antrum in that it determines the volume of the granulosa layer and thus the
level of estrogen production.


  Many database resources, such as Reactome, collect manually annotated
reactions, interactions, and pathways from peer-reviewed publications. The
interactors (e.g., a protein), interactions, and pathways in these data
resources are often represented as instances in using BioPAX, a standard
pathway data exchange format. However, these interactions are better
represented as classes (or universals) since they always occur given
appropriate conditions. This study aims to represent various human interaction
pathways and networks as classes via a formal ontology aligned with the Basic
Formal Ontology (BFO). Towards this goal, the Human Interaction Network
Ontology (HINO) was generated by extending the BFO-aligned Interaction Network
Ontology (INO). All human pathways and associated processes and interactors
listed in Reactome and represented in BioPAX were first converted to ontology
classes by aligning them under INO. Related terms and associated relations and
hierarchies from external ontologies (e.g., CHEBI and GO) were also retrieved
and imported into HINO. HINO ontology terms were resolved in the linked
ontology data server Ontobee. The RDF triples stored in the RDF triple store
are queryable through a SPARQL program. Such an ontology system supports
advanced pathway data integration and applications.


  Metagenomics is an emerging field of molecular biology concerned with
analyzing the genomes of environmental samples comprising many different
diverse organisms. Given the nature of metagenomic data, one usually has to
sequence the genomic material of all organisms in a batch, leading to a mix of
reads coming from different DNA sequences. In deep high-throughput sequencing
experiments, the volume of the raw reads is extremely high, frequently
exceeding 600 Gb. With an ever increasing demand for storing such reads for
future studies, the issue of efficient metagenomic compression becomes of
paramount importance. We present the first known approach to metagenome read
compression, termed MCUIUC (Metagenomic Compression at UIUC). The gist of the
proposed algorithm is to perform classification of reads based on unique
organism identifiers, followed by reference-based alignment of reads for
individually identified organisms, and metagenomic assembly of unclassified
reads. Once assembly and classification are completed, lossless reference based
compression is performed via positional encoding. We evaluate the performance
of the algorithm on moderate sized synthetic metagenomic samples involving 15
randomly selected organisms and describe future directions for improving the
proposed compression method.


  The use of an extra SDS separation in a different buffer system provide a
technique for deconvoluting 2D gel spots made of several proteins (Colignon et
al. Proteomics, 2013, 13, 2077-2082). This technique keeps the quantitative
analysis of the protein amounts and combines it with a strongly improved
identification process by mass spectrometry, removing identification
ambiguities in most cases. In some favorable cases, posttranslational variants
can be separated by this procedure. This versatile and easy to use technique is
anticipated to be a very valuable addition to the toolbox used in 2D gel-based
proteomics.


  While prefractionation has previously been shown to improve results in MS
analysis, a novel combination provides an additional dimension of separation:
protein fractionation by SDS-PAGE followed by IEF of tryptic peptides before
separation by RP-LC [Atanassov and Urlaub, Proteomics 2013, 13, 2947-2955].
This three-step separation procedure prior to MS/MS substantially increases
proteome coverage and represents a further step toward a more comprehensive
analysis of complex proteomes.


  We introduce a parallel algorithmic architecture for metagenomic sequence
assembly, termed MetaPar, which allows for significant reductions in assembly
time and consequently enables the processing of large genomic datasets on
computers with low memory usage. The gist of the approach is to iteratively
perform read (re)classification based on phylogenetic marker genes and
assembler outputs generated from random subsets of metagenomic reads. Once a
sufficiently accurate classification within genera is performed, de novo
metagenomic assemblers (such as Velvet or IDBA-UD) or reference based
assemblers may be used for contig construction. We analyze the performance of
MetaPar on synthetic data consisting of 15 randomly chosen species from the
NCBI database through the effective gap and effective coverage metrics.


  Growth of human population shows no signs of stagnation. The only small
disturbance is identified as being probably associated with the coinciding
impacts of five demographic catastrophes. The concept of the Epoch of
Malthusian Stagnation is convincingly contradicted by empirical evidence.


  Given a finite set $X$, a collection $\mathcal{T}$ of rooted phylogenetic
trees on $X$ and an integer $k$, the Hybridization Number problem asks if there
exists a phylogenetic network on $X$ that displays all trees from $\mathcal{T}$
and has reticulation number at most $k$. We show two kernelization algorithms
for Hybridization Number, with kernel sizes $4k(5k)^t$ and $20k^2(\Delta^+-1)$
respectively, with $t$ the number of input trees and $\Delta^+$ their maximum
outdegree. Experiments on simulated data demonstrate the practical relevance of
these kernelization algorithms. In addition, we present an $n^{f(k)}t$-time
algorithm, with $n=|X|$ and $f$ some computable function of $k$.


  The dynamics of economies and infectious disease are inexorably linked:
economic well-being influences health (sanitation, nutrition, treatment
capacity, etc.) and health influences economic well-being (labor productivity
lost to sickness and disease). Often societies are locked into "poverty traps"
of poor health and poor economy. Here, using a simplified coupled
disease-economic model with endogenous capital growth we demonstrate the
formation of poverty traps, as well as ways to escape them. We suggest two
possible mechanisms of escape both motivated by empirical data: one, through an
influx of capital (development aid), and another through changing the
percentage of GDP spent on healthcare. We find that a large influx of capital
is successful in escaping the poverty trap, but increasing health spending
alone is not. Our results demonstrate that escape from a poverty trap may be
possible, and carry important policy implications in the world-wide
distribution of aid and within-country healthcare spending.


  The present study reports on the bio-organic composition of a deep-sea
venting hydrothermal system originating from arc volcanism; the origin of the
particulates in hydrothermal fluids from the Suiyo Seamount in the southern
Izu-Bonin (Ogasawara) Arc is discussed with regard to amino compounds. Chimney
samples on deep-sea hydrothermal systems and core samples at Suiyo Seamount
were determined for amino acids, and occasionally amino sugars. Two types of
chimney samples were obtained from active hydrothermal systems by submersible
vehicles: one was natural chimney (NC) on a hydrothermal natural vent; the
other was artificial chimneys (AC), mainly formed by the growth and deposition
of sulfide-rich particulate components in a Kuwabara-type in situ incubator (KI
incubator). Total hydrolyzed amino acids (THAA) and hydrolyzed hexosamines (HA)
in AC ranged from 10.7 nmol/g to 64.0 nmol/g and from 0 nmol/g to 8.1 nmol/g,
respectively, while THAA in hydrothermally altered core samples ranged from
26.0 nmol/g to 107.4 nmol/g-rock. The THAA was thus of the same order of
magnitude in AC and core samples. The representative biochemical indicator
ratios of {\beta}-alanine/aspartic acid and {\gamma}-aminobutyric acid/glutamic
acid were low. This is consistent with a large microbial population and labile
subterranean biogenic organic compounds. A micro energy dispersive X-ray
spectrometer ({\mu}EDX) was used to show that the formation of massive pyrite-
and chalcopyrite-ores, the major matrixes of the natural chimney structures,
were processed during hydrothermal activity. These data suggest that labile
particulate organic matter (POM) and/or dissolved organic matter (DOM) from
hydrothermal fluid were continuously trapped to form concentrated organic
matter in a sulfide-rich matrix.


  Due to the availability of large molecular data-sets, covariance models are
increasingly used to describe the structure of genetic variation as an
alternative to more heavily parametrised biological models. We focus here on a
class of parametric covariance models that received sustained attention lately
and show that the conditions under which they are valid mathematical models
have been overlooked so far. We provide rigorous results for the construction
of valid covariance models in this family. We also outline how to construct
alternative covariance models for the analysis of geographical variation that
are both mathematically well behaved and easily implementable.


  Online genealogy datasets contain extensive information about millions of
people and their past and present family connections. This vast amount of data
can assist in identifying various patterns in human population. In this study,
we present methods and algorithms which can assist in identifying variations in
lifespan distributions of human population in the past centuries, in detecting
social and genetic features which correlate with human lifespan, and in
constructing predictive models of human lifespan based on various features
which can easily be extracted from genealogy datasets.
  We have evaluated the presented methods and algorithms on a large online
genealogy dataset with over a million profiles and over 9 million connections,
all of which were collected from the WikiTree website. Our findings indicate
that significant but small positive correlations exist between the parents'
lifespan and their children's lifespan. Additionally, we found slightly higher
and significant correlations between the lifespans of spouses. We also
discovered a very small positive and significant correlation between longevity
and reproductive success in males, and a small and significant negative
correlation between longevity and reproductive success in females. Moreover,
our machine learning algorithms presented better than random classification
results in predicting which people who outlive the age of 50 will also outlive
the age of 80.
  We believe that this study will be the first of many studies which utilize
the wealth of data on human populations, existing in online genealogy datasets,
to better understand factors which influence human lifespan. Understanding
these factors can assist scientists in providing solutions for successful
aging.


  Neutral models aspire to explain biodiversity patterns in ecosystems where
species difference can be neglected, as it might occur at a specific trophic
level, and perfect symmetry is assumed between species. Voter-like models
capture the essential ingredients of the neutral hypothesis and represent a
paradigm for other disciplines like social studies and chemical reactions. In a
system where each individual can interact with all the other members of the
community, the typical time to reach an absorbing state with a single species
scales linearly with the community size. Here we show, by using a rigorous
approach within a large deviation principle and confirming previous approximate
and numerical results, that in a heterogeneous voter model the typical time to
reach an absorbing state scales exponentially with the system size, suggestive
of an asymptotic active phase.


  Acute Respiratory Distress Syndrome (ARDS) is a very severe syndrome leading
to respiratory failure and subsequent mortality. Sepsis is one of the leading
causes of ARDS. Thus, extracellular bacteria play an important role in the
pathophysiology of ARDS. Overactivated neutrophils are the major effector cells
in ARDS. Thus, extracellular bacteria triggered TH17-like innate immunity with
neutrophil activation might accounts for the etiology of ARDS. Here, microarray
analysis was employed to describe TH17-like innate immunity-related cytokine
including TGF-{\beta} and IL-6 up-regulation in whole blood of ARDS patients.
It was found that the innate TH17-related TLR1,2,4,5,8, HSP70, G-CSF, GM-CSF,
complements, defensin, PMN chemokines, cathepsins, Fc receptors, NCFs, FOS,
JunB, CEBPs, NFkB, and leukotriene B4 are all up-regulated. TGF-{\beta}
secreting Treg cells play important roles in lung fibrosis. Up-regulation of
Treg associated STAT5B and TGF-{\beta} with down-regulation of MHC genes, TCR
genes, and co-stimulation molecule CD86 are noted. Key TH17 transcription
factors, STAT3 and ROR{\alpha}, are down-regulated. Thus, the full adaptive
TH17 helper CD4 T cells may not be successfully triggered. Many fibrosis
promoting genes are also up-regulated including MMP8, MMP9, FGF13, TIMP1,
TIMP2, PLOD1, P4HB, P4HA1, PDGFC, HMMR, HS2ST1, CHSY1, and CSGALNACT. Failure
to induce successful adaptive immunity could also attribute to ARDS
pathogenesis. Thus, ARDS is actually a TH17-like and Treg immune disorder.


  We compare the likelihood of different socially relevant features to allow
the evolutionary emergence and maintenance of cooperation in a generalized
variant of the iterated Prisoners Dilemma game. Results show that the average
costs/benefit balance of cooperation is the primary constraint for its
establishment and maintenance. Behavior increasing inclusive fitness such as
assortation, homophily, kin-selection and tagging of individuals, is second in
importance. Networks characteristics were the least important in favoring the
establishment and maintenance of cooperation, despite being the most popular in
recent research on the subject. Results suggest that inclusive fitness theory
with its expansions to include assortative and economic considerations is more
general, powerful and relevant in analyzing social phenomena than kin selection
theory with its emphasis on genetic relatedness. Merging economics with
evolutionary theory will be necessary to reveal more about the nature of social
dynamics.


  We explore the probabilistic structure of DNA in a number of bacterial
genomes and conclude that a form of Markovianness is present at the boundaries
between coding and non-coding regions, that is, the sequence of START and STOP
codons annotated for the bacterial genome. This sequence is shown to satisfy a
conditional independence property which allows its governing Markov chain to be
uniquely identified from the abundances of START and STOP codons. Furthermore,
the annotated sequence is shown to comply with Chargaff's second parity rule at
the codon level.


  Molecular communication emerges as a promising communication paradigm for
nanotechnology. However, solid mathematical foundations for
information-theoretic analysis of molecular communication have not yet been
built. In particular, no one has ever proven that the channel coding theorem
applies for molecular communication, and no relationship between information
rate capacity (maximum mutual information) and code rate capacity (supremum
achievable code rate) has been established. In this paper, we focus on a major
subclass of molecular communication - the diffusion-based molecular
communication. We provide solid mathematical foundations for information theory
in diffusion-based molecular communication by creating a general
diffusion-based molecular channel model in measure-theoretic form and prove its
channel coding theorems. Various equivalence relationships between statistical
and operational definitions of channel capacity are also established, including
the most classic information rate capacity and code rate capacity. As
byproducts, we have shown that the diffusion-based molecular channel is with
"asymptotically decreasing input memory and anticipation" and "d-continuous".
Other properties of diffusion-based molecular channel such as stationarity or
ergodicity are also proven.


  A phenomenon that strongly influences the demography of small introduced
populations and thereby potentially their genetic diversity is the Allee
effect, a reduction in population growth rates at small population sizes. We
take a stochastic modeling approach to investigate levels of genetic diversity
in populations that successfully overcame a strong demographic Allee effect, a
scenario in which populations smaller than a certain critical size are expected
to decline. Our results indicate that compared to successful populations
without Allee effect, successful Allee-effect populations tend to 1) derive
from larger founder population sizes and thus have a higher initial amount of
genetic variation, 2) spend fewer generations at small population sizes where
genetic drift is particularly strong, and 3) spend more time around the
critical population size and thus experience more drift there. Altogether, the
Allee effect can either increase or decrease genetic diversity, depending on
the average founder population size. In the case of multiple introduction
events, there is an additional increase in diversity because Allee-effect
populations tend to derive from a larger number of introduction events than
other populations. Finally, we show that given genetic data from sufficiently
many populations, we can statistically infer the critical population size.


  It is now well established that the cell is a highly crowded environment.
Yet, the effects of crowding on the dynamics of signaling pathways, gene
regulation networks and metabolic networks are still largely unknown. Crowding
can alter both molecular diffusion and the equilibria of biomolecular
reactions. In this review, we first discuss how diffusion can affect
biochemical networks. Diffusion of transcription factors can increase noise in
gene expression, while diffusion of proteins between intracellular compartments
or between cells can reduce concentration fluctuations. In push-pull networks
diffusion can impede information transmission, while in multi-site protein
modification networks diffusion can qualitatively change the macroscopic
response of the system, such as the loss or emergence of bistability. Moreover,
diffusion can directly change the metabolic flux. We describe how crowding
affects diffusion, and thus how all these phenomena are influenced by crowding.
Yet, a potentially more important effect of crowding on biochemical networks is
mediated via the shift in the equilibria of bimolecular reactions, and we
provide computational evidence that supports this idea. Finally, we discuss how
the effects of crowding can be incorporated in models of biochemical networks.


  The evolution of ants is marked by remarkable adaptations that allowed the
development of very complex social systems. To identify how ant-specific
adaptations are associated with patterns of molecular evolution, we searched
for signs of positive selection on amino-acid changes in proteins. We
identified 24 functional categories of genes which were enriched for positively
selected genes in the ant lineage. We also reanalyzed genome-wide datasets in
bees and flies with the same methodology, to check whether positive selection
was specific to ants or also present in other insects. Notably, genes
implicated in immunity were enriched for positively selected genes in the three
lineages, ruling out the hypothesis that the evolution of hygienic behaviors in
social insects caused a major relaxation of selective pressure on immune genes.
Our scan also indicated that genes implicated in neurogenesis and olfaction
started to undergo increased positive selection before the evolution of
sociality in Hymenoptera. Finally, the comparison between these three lineages
allowed us to pinpoint molecular evolution patterns that were specific to the
ant lineage. In particular, there was ant-specific recurrent positive selection
on genes with mitochondrial functions, suggesting that mitochondrial activity
was improved during the evolution of this lineage. This might have been an
important step toward the evolution of extreme lifespan that is a hallmark of
ants.


  Currently, there are two major theories for the pathogenesis of sepsis:
hyperimmune and hypoimmune. Hyperimmune theory suggests that cytokine storm
causes the symptoms of sepsis. On the contrary, hypoimmune theory suggests that
immunosuppression causes the manifestations of sepsis. By using microarray
study, this study implies that hyperactivity of TH17-like innate immunity and
failure of adaptive immunity are noted in sepsis patients. I find out that
innate immunity related genes are significantly up-regulated including CD14,
TLR1,2,4,5,8, HSP70, CEBP proteins, AP1(JUNB, FOSL2), TGF-{\beta}, IL-6,
TGF-{\alpha}, CSF2 receptor, TNFRSF1A, S100A binding proteins, CCR2, formyl
peptide receptor2, amyloid proteins, pentraxin, defensins, CLEC5A, whole
complement machinery, CPD, NCF, MMP, neutrophil elastase, caspases, IgG and IgA
Fc receptors(CD64, CD32), ALOX5, PTGS, LTB4R, LTA4H, and ICAM1. Majority of
adaptive immunity genes are down-regulated including MHC related genes, TCR
genes, granzymes/perforin, CD40, CD8, CD3, TCR signaling, BCR signaling, T & B
cell specific transcription factors, NK killer receptors, and TH17 helper
specific transcription factors(STAT3, RORA, REL). In addition, Treg related
genes are up-regulated including TGF{\beta}, IL-15, STAT5B, SMAD2/4, CD36, and
thrombospondin. Thus, both hyperimmune and hypoimmune play important roles in
the pathophysiology of sepsis.


  In mathematical studies of the dynamics of multi-strain diseases caused by
antigenically diverse pathogens, there is a substantial interest in analytical
insights. Using the example of a generic model of multi-strain diseases with
cross-immunity between strains, we show that a significant understanding of the
stability of steady states and possible dynamical behaviours can be achieved
when the symmetry of interactions between strains is taken into account.
Techniques of equivariant bifurcation theory allow one to identify the type of
possible symmetry-breaking Hopf bifurcation, as well as to classify different
periodic solutions in terms of their spatial and temporal symmetries. The
approach is also illustrated on other models of multi-strain diseases, where
the same methodology provides a systematic understanding of bifurcation
scenarios and periodic behaviours. The results of the analysis are quite
generic, and have wider implications for understanding the dynamics of a large
class of models of multi-strain diseases.


  Lipid monolayers and bilayers have been used as experimental models for the
investigation of membrane thermal transitions. The main transition takes place
near ambient temperatures for several lipids and reflects the order-disorder
transition of lipid hydrocarbonic chains, which is accompanied by a small
density gap. Equivalence between the transitions in the two systems has been
argued by several authors. The two-state statistical model adopted by numerous
authors for different properties of the membrane, such as permeability,
diffusion, mixture or insertion of cholesterol or protein, is inadequate for
the description of charged membranes, since it lacks a proper description of
surface density. We propose a lattice solution model which adds interactions
with water molecules to lipid-lipid interactions and obtain its thermal
properties under a mean-field approach. Density variations, although
concomitant with chain order variations, are independent of the latter. The
model presents both chain order and gas-liquid transitions, and extends the
range of applicability of previous models, yielding Langmuir isotherms in the
full range of pressures and areas.


  Matrix population models have proved popular and useful for studying
stage-structured populations in quantitative ecology. The goal of this paper is
to develop and analyze a general matrix population model that is applicable to
the population dynamics of the invasive round goby fish that incorporates both
stage structure---larvae, juveniles, and adults---and dispersion. Specifically,
we address the issue of whether or not dispersion can amplify the net
reproductive number. To this end, we will first review the mathematics of
matrix population models with dispersion, particularly those with an Usher
demography matrix. Techniques for computing the net reproductive number, like
the graph reduction method of de-Camino-Beck and Lewis, will be discussed. A
common theme will be the usefulness of submatrices of relevant matrices
obtained by the expunging of rows and columns corresponding to non-newborns.
Finally, examples will be provided, including the calculation of the net
reproductive number for multiple regions using the graph reduction method of
de-Camino-Beck and Lewis, examples where dispersion results in a total
population flourishing when the populations would otherwise go extinct, and an
application to the invasive round goby fish.


  Annotations of gene structures and regulatory elements can inform genome-wide
association studies (GWAS). However, choosing the relevant annotations for
interpreting an association study of a given trait remains challenging. We
describe a statistical model that uses association statistics computed across
the genome to identify classes of genomic element that are enriched or depleted
for loci that influence a trait. The model naturally incorporates multiple
types of annotations. We applied the model to GWAS of 18 human traits,
including red blood cell traits, platelet traits, glucose levels, lipid levels,
height, BMI, and Crohn's disease. For each trait, we evaluated the relevance of
450 different genomic annotations, including protein-coding genes, enhancers,
and DNase-I hypersensitive sites in over a hundred tissues and cell lines. We
show that the fraction of phenotype-associated SNPs that influence protein
sequence ranges from around 2% (for platelet volume) up to around 20% (for LDL
cholesterol); that repressed chromatin is significantly depleted for SNPs
associated with several traits; and that cell type-specific DNase-I
hypersensitive sites are enriched for SNPs associated with several traits (for
example, the spleen in platelet volume). Finally, by re-weighting each GWAS
using information from functional genomics, we increase the number of loci with
high-confidence associations by around 5%.


  A quasispecies is the stationary state of a set of interrelated genotypes
that evolve according to the usual principles of selection and mutation.
Quasispecies studies have invariably concentrated on the possibility of errors
during genotype replication and their role in promoting either the survival or
the demise of the quasispecies. In a previous work [V. C. Barbosa, R.
Donangelo, and S. R. Souza, J. Theor. Biol. 312, 114 (2012)], we introduced a
network model of quasispecies dynamics, based on a single probability parameter
($p$) and capable of addressing several plausibility issues of previous models.
Here we extend that model by pairing its network with another one aimed at
modeling the dynamics of the immune system when confronted with the
quasispecies. The new network is based on the idiotypic-network model of
immunity and, together with the previous one, constitutes a network model of
interacting genotypes and idiotypes. The resulting model requires further
parameters and as a consequence leads to a vast phase space. We have focused on
a particular niche in which it is possible to observe the trade-offs involved
in the quasispecies' survival or destruction. Within this niche, we give
simulation results that highlight some key preconditions for quasispecies
survival. These include a minimum initial abundance of genotypes relative to
that of the idiotypes and a minimum value of $p$. The latter, in particular, is
to be contrasted with the stand-alone quasispecies network of our previous
work, in which arbitrarily low values of $p$ constitute a guarantee of
quasispecies survival.


  Hysteresis and bet-hedging (random choice of phenotypes) are two different
observations typically linked with multiplicity of phenotypes in biological
systems. Hysteresis can be viewed as form of the system's persistent memory of
past environmental conditions, while bet-hedging is a diversification strategy
not necessarily associated with any memory. It has been shown that bet-hedging
can increase population growth when phenotype adjusts its switching probability
in response to environmental inputs. Although memory and hysteresis have been
used to generate environment dependent phenotype switching probabilities, their
exact connection to bet-hedging have remained unclear. In this work, using a
simple model that takes into account phenotype switching as well as lag phase
in the population growth occurring after the phenotype switching, it is shown
that memory and hysteresis emerge naturally and are firmly linked to
bet-hedging as organisms attempt to optimize their population growth rate. The
optimal ``magnitude'' of hysteresis is explained to be associated with
stochastic resonance where the characteristic time between subsequent phenotype
switching events is linked to the lag phase delay. Furthermore, hysteretic
switching strategy is shown not to confer any additional population growth
advantage if the environment varies periodically in a deterministic fashion.
This suggests that, while bet-hedging may evolve under some conditions even in
deterministic environments, memory and hysteresis is probably the result of
environmental uncertainty in the presence of lag phase in the switching of
phenotypes.


  Here we introduce a general class of multiple calibration birth-death tree
priors for use in Bayesian phylogenetic inference. All tree priors in this
class separate ancestral node heights into a set of "calibrated nodes" and
"uncalibrated nodes" such that the marginal distribution of the calibrated
nodes is user-specified whereas the density ratio of the birth-death prior is
retained for trees with equal values for the calibrated nodes. We describe two
formulations, one in which the calibration information informs the prior on
ranked tree topologies, through the (conditional) prior, and the other which
factorizes the prior on divergence times and ranked topologies, thus allowing
uniform, or any arbitrary prior distribution on ranked topologies. While the
first of these formulations has some attractive properties the algorithm we
present for computing its prior density is computationally intensive. On the
other hand, the second formulation is always computationally efficient. We
demonstrate the utility of the new class of multiple-calibration tree priors
using both small simulations and a real-world analysis and compare the results
to existing schemes. The two new calibrated tree priors described in this paper
offer greater flexibility and control of prior specification in calibrated
time-tree inference and divergence time dating, and will remove the need for
indirect approaches to the assessment of the combined effect of calibration
densities and tree process priors in Bayesian phylogenetic inference.


  Unstable angina is common clinical manifestation of atherosclerosis. However,
the detailed pathogenesis of unstable angina is still not known. Here, I
propose that unstable angina is a mixed TH17 and TH1 immune disorder. By using
microarray analysis, I find out that TH1 and TH17 related cytokine, cytokine
receptor, chemokines, complement, immune-related transcription factors,
anti-bacterial genes, Toll-like receptors, and heat shock proteins are all
up-regulated in peripheral leukocytes of unstable angina. In addition,
H-ATPase, glycolytic genes, platelet and RBC related genes are also
up-regulated in peripheral leukocytes of during unstable angina. It also
implies that atherosclerosis is a mixed TH17 and TH1 autoimmune disease. If we
know the etiology of unstable angina as well as atherosclerosis better, we can
have better methods to control and prevent this detrimental illness.


  We review pro and contra of the hypothesis that generic polymer properties of
topological constraints are behind many aspects of chromatin folding in
eukaryotic cells. For that purpose, we review, first, recent theoretical and
computational findings in polymer physics related to concentrated,
topologically-simple (unknotted and unlinked) chains or a system of chains.
Second, we review recent experimental discoveries related to genome folding.
Understanding in these fields is far from complete, but we show how looking at
them in parallel sheds new light on both.


  A strong demographic Allee effect in which the expected population growth
rate is negative below a certain critical population size can cause high
extinction probabilities in small introduced populations. However, many species
are repeatedly introduced to the same location and eventually one population
may overcome the Allee effect by chance. With the help of stochastic models, we
investigate how much genetic diversity such successful populations harbour on
average and how this depends on offspring-number variation, an important source
of stochastic variability in population size. We find that with increasing
variability, the Allee effect increasingly promotes genetic diversity in
successful populations. Successful Allee-effect populations with highly
variable population dynamics escape rapidly from the region of small population
sizes and do not linger around the critical population size. Therefore, they
are exposed to relatively little genetic drift. We show that here---unlike in
classical population genetics models---the role of offspring-number variation
cannot be accounted for by an effective-population-size correction. Thus, our
results highlight the importance of detailed biological knowledge, in this case
on the probability distribution of family sizes, when predicting the
evolutionary potential of newly founded populations or when using genetic data
to reconstruct their demographic history.


  Nearly all living organisms, including some bacterial species, exhibit
biological processes with a period of about 24 h called circadian (from the
Latin circa, about, dies, day) rhythms. These rhythms allow living organisms to
anticipate the daily alternation of light and darkness. Experiments carried out
in cyanobacteria have shown the adaptive value of circadian clocks. In these
experiments a wild type cyanobacterial strain (with a 24 h circadian rhythm)
and a mutant strain (with a longer or shorter period) grow in competition. In
different experiments the external light dark cycle was changed in order to
match the circadian period of the different strains, revealing that the strain
whose circadian period matches the light-dark has a larger fitness. As a
consequence the initial population of one strain grows while the other decays.
  These experiments were made under fixed light and dark intervals. However, in
Nature this relationship changes according to the season. Therefore, seasonal
changes in light could affect the results of the competition. Using a
theoretical model we analyze how modulation of light can change the survival of
the different cyanobacterial strains. Our results show that there is a clear
shift in the competition due to the modulation of light, which could be
verified experimentally.


  Background: The reconstruction of the phylogenetic tree topology of four taxa
is, still nowadays, one of the main challenges in phylogenetics. Its
difficulties lie in considering not too restrictive evolutionary models, and
correctly dealing with the long-branch attraction problem. The correct
reconstruction of 4-taxon trees is crucial for making quartet-based methods
work and being able to recover large phylogenies.
  Results: In this paper we consider an expectation-maximization method for
maximizing the likelihood of (time nonhomogeneous) evolutionary Markov models
on trees. We study its success on reconstructing 4-taxon topologies and its
performance as input method in quartet-based phylogenetic reconstruction
methods such as QFIT and QuartetSuite. Our results show that the method
proposed here outperforms neighbor-joining and the usual (time-homogeneous
continuous-time) maximum likelihood methods on 4-leaved trees with
among-lineage instantaneous rate heterogeneity, and perform similarly to usual
continuous-time maximum-likelihood when data satisfies the assumptions of both
methods.
  Conclusions: The method presented in this paper is well suited for
reconstructing the topology of any number of taxa via quartet-based methods and
is highly accurate, specially regarding largely divergent trees and time
nonhomogeneous data.


  We report on a genome-wide scan for introgression in the house mouse (Mus
musculus domesticus) involving the Algerian mouse (Mus spretus), using samples
from the ranges of sympatry and allopatry in Africa and Europe. Our analysis
reveals wide variability in introgression signatures along the genomes, as well
as across the samples. We find that fewer than half of the autosomes in each
genome harbor all detectable introgression, while the X chromosome has none.
Further, European mice carry more M. spretus alleles than the sympatric African
ones. Using the length distribution and sharing patterns of introgressed
genomic tracts across the samples, we infer, first, that at least three
distinct hybridization events involving M. spretus have occurred, one of which
is ancient, and the other two are recent (one presumably due to warfarin
rodenticide selection). Second, several of the inferred introgressed tracts
contain genes that are likely to confer adaptive advantage. Third, introgressed
tracts might contain driver genes that determine the evolutionary fate of those
tracts. Further, functional analysis revealed introgressed genes that are
essential to fitness, including the Vkorc1 gene, which is implicated in
rodenticide resistance, and olfactory receptor genes. Our findings highlight
the extent and role of introgression in nature, and call for careful analysis
and interpretation of house mouse data in evolutionary and genetic studies.


  Recon 2 is a highly curated reconstruction of the human metabolic network.
Whilst the network is state of the art, it has shortcomings, including the
presence of unbalanced reactions involving generic metabolites. By replacing
these generic molecules with each of their specific instances, we can ensure
full elemental balancing, in turn allowing constraint-based analyses to be
performed. The resultant model, called Recon 2.1, is an order of magnitude
larger than the original.


  Full likelihood inference under Kingman's coalescent is a computationally
challenging problem to which importance sampling (IS) and the product of
approximate conditionals (PAC) method have been applied successfully. Both
methods can be expressed in terms of families of intractable conditional
sampling distributions (CSDs), and rely on principled approximations for
accurate inference. Recently, more general $\Lambda$- and $\Xi$-coalescents
have been observed to provide better modelling fits to some genetic data sets.
We derive families of approximate CSDs for finite sites $\Lambda$- and
$\Xi$-coalescents, and use them to obtain "approximately optimal" IS and PAC
algorithms for $\Lambda$-coalescents, yielding substantial gains in efficiency
over existing methods.


  Zero-order ultrasensitivity (ZOU) is a long known and interesting phenomenon
in enzyme networks. Here, a substrate is reversibly modified by two
antagonistic enzymes (a "push-pull" system) and the fraction in modified state
undergoes a sharp switching from near-zero to near-unity at a critical value of
the ratio of the enzyme concentrations, under saturation conditions. ZOU and
its extensions have been studied for several decades now, ever since the
seminal paper of Goldbeter and Koshland (1981); however, a complete
probabilistic treatment, important for the study of fluctuations in finite
populations, is still lacking. In this paper, we study ZOU using a modular
approach, akin to the total quasi-steady state approximation (tQSSA). This
approach leads to a set of Fokker-Planck (drift-diffusion) equations for the
probability distributions of the intermediate enzyme-bound complexes, as well
as the modified/unmodified fractions of substrate molecules. We obtain explicit
expressions for various average fractions and their fluctuations in the linear
noise approximation (LNA). The emergence of a "critical point" for the
switching transition is rigorously established. New analytical results are
derived for the average and variance of the fractional substrate concentration
in various chemical states in the near-critical regime. For the total fraction
in the modified state, the variance is shown to be a maximum near the critical
point and decays algebraically away from it, similar to a second-order phase
transition. The new analytical results are compared with existing ones as well
as detailed numerical simulations using a Gillespie algorithm.


  Malaria is a serious global health problem which is especially devastating to
the developing world. Mosquitoes are the carriers of the parasite responsible
for the disease, and hence malaria control programs focus on controlling
mosquito populations. This is done primarily through the spraying of
insecticides, or through the use of insecticide treated bed nets. However,
usage of these insecticides exerts massive selection pressure on mosquitoes,
resulting in insecticide resistant mosquito breeds. Hence, developing
alternative strategies is crucial for sustainable malaria control. Here we
explore the usefulness of paratransgenesis, i.e., introducing genetically
engineered bacteria which secrete anti-plasmodium molecules, inside the
mosquito midgut. The bacteria enter a mosquito's midgut when it drinks from a
sugar bait, i.e., a sugar solution containing the bacterium. We formulate a
mathematical model for evaluating the number of such baits required for
preventing an outbreak. We study scenarios where vectors and hosts mix
homogeneously as well as heterogeneously. We perform a full stability analysis
and calculate the basic reproductive number for both the cases. Additionally,
for the heterogeneous mixing scenario, we propose a targeted bait distribution
strategy. The optimal bait allocation is calculated and is found to be
extremely efficient in terms of bait usage. Our analyses suggest that
paratransgenesis can prevent an outbreak, and hence it offers a viable and
sustainable path to malaria control.


  The conventional cancer stem cell (CSC) theory indicates a hierarchy of CSCs
and non-stem cancer cells (NSCCs), that is, CSCs can differentiate into NSCCs
but not vice versa. However, an alternative paradigm of CSC theory with
reversible cell plasticity among cancer cells has received much attention very
recently. Here we present a generalized multi-phenotypic cancer model by
integrating cell plasticity with the conventional hierarchical structure of
cancer cells. We prove that under very weak assumption, the nonlinear dynamics
of multi-phenotypic proportions in our model has only one stable steady state
and no stable limit cycle. This result theoretically explains the phenotypic
equilibrium phenomena reported in various cancer cell lines. Furthermore,
according to the transient analysis of our model, it is found that cancer cell
plasticity plays an essential role in maintaining the phenotypic diversity in
cancer especially during the transient dynamics. Two biological examples with
experimental data show that the phenotypic conversions from NCSSs to CSCs
greatly contribute to the transient growth of CSCs proportion shortly after the
drastic reduction of it. In particular, an interesting overshooting phenomenon
of CSCs proportion arises in three-phenotypic example. Our work may pave the
way for modeling and analyzing the multi-phenotypic cell population dynamics
with cell plasticity.


  When organisms adapt to spatially heterogeneous environments, selection may
drive divergence at multiple genes. If populations under divergent selection
also exchange migrants, we expect genetic differentiation to be high at
selected loci, relative to the baseline caused by migration and genetic drift.
Indeed, empirical studies have found peaks of putatively adaptive
differentiation. These are highly variable in length, some of them extending
over several hundreds of thousands of base pairs. How can such 'islands of
divergence' be explained? Physical linkage produces elevated levels of
differentiation at loci close to genes under selection. However, whether this
is enough to account for the observed patterns of divergence is not well
understood. Here, we investigate the fate of a locally beneficial mutation that
arises in linkage to an existing migration-selection polymorphism and derive
two important quantities: the probability that the mutation becomes
established, and the expected time to its extinction. We find that intermediate
levels of recombinations are sometimes favourable, and that physical linkage
can lead to strongly elevated invasion probabilities and extinction times. We
provide a rule of thumb for when this is the case. Moreover, we quantify the
long-term effect of polygenic local adaptation on linked neutral variation.