  Assessing the statistical significance of an observed 2x2 contingency table
can easily be accomplished using Fisher's exact test (FET). However, if the
cell entries are continuous or represent values inferred from a continuous
parametric model, then FET cannot be applied. Such tables arise frequently in
areas of biostatistical research including population genetics and evolutionary
genomics, where cell entries are estimated by computational methods and result
in cell entries drawn from the non-negative real line R+. Simply rounding cell
entries to conform to the assumptions of FET is an ill-suited approach that we
show creates problems related to both type-I and type-II errors. Pearson's
chi^2 test for independence, while technically applicable, is not often
effective for these tables, as the test has several limiting assumptions that
make application of this method inadvisable in many common instances
(particularly with small cell entries). Here we develop a novel method for
tables with continuous entries, which we term continuous Fisher's Exact Test
(cFET). Through simulations, we show that cFET has a close-to-uniform
distribution of p-values under the null hypothesis of independence, and more
power when applied to tables where the null hypothesis is false (compared to
FET applied to rounded cell entries). We apply cFET to an example from
comparative genomics to confirm an overall increased evolutionary rate among
primates compared to rodents, and identify several genes that show particularly
elevated evolutionary rates in primates. Some of these genes exhibit signatures
of continued positive selection along the human lineage since our divergence
with chimpanzee 5-7 million years ago, as well as ongoing selection in modern
humans.


  Polymerases select nucleotides before incorporating them for chemical
synthesis during gene replication or transcription. How the selection proceeds
stepwise efficiently to achieve sufficiently high fidelity and speed is
essential for polymerase function. We examined step-by-step selections that
have conformational transition rates tuned one at time in the polymerase
elongation cycle, with a controlled differentiation free energy at each
checkpoint. The elongation is sustained at non-equilibrium steady state with
constant free energy input and heat dissipation. It is found that error
reduction capability does not improve for selection checkpoints down the
reaction path. Hence, it is essential to select early to achieve an efficient
fidelity control. In particular, for two consecutive selections that reject the
wrong substrate back and inhibit it forward from a same kinetic state, the same
error rates are obtained at the same free energy differentiation. The initial
screening is indispensible for maintaining the elongation speed high, as the
wrong nucleotides can be removed quickly and replaced by the right nucleotides
at the entry. Overall, the elongation error rate can be repeatedly reduced
through multiple selection checkpoints. The study provides a theoretical
framework to conduct further detailed researches, and assists engineering and
redesign of related enzymes.


  Methylation of CpG dinucleotides is a prevalent epigenetic modification that
is required for proper development in vertebrates, and changes in CpG
methylation are essential to cellular differentiation. Genome-wide DNA
methylation assays have become increasingly common, and recently distinct
stages across differentiating cellular lineages have been assayed. How- ever,
current methods for modeling methylation dynamics do not account for the
dependency structure between precursor and dependent cell types. We developed a
continuous-time Markov chain approach, based on the observation that changes in
methylation state over tissue differentiation can be modeled similarly to DNA
nucleotide changes over evolutionary time. This model explicitly takes
precursor to descendant relationships into account and enables inference of CpG
methylation dynamics. To illustrate our method, we analyzed a high-resolution
methylation map of the differentiation of mouse stem cells into several blood
cell types. Our model can successfully infer unobserved CpG methylation states
from observations at the same sites in related cell types (90% correct), and
this approach more accurately reconstructs missing data than imputation based
on neighboring CpGs (84% correct). Additionally, the single CpG resolution of
our methylation dynamics estimates enabled us to show that DNA sequence context
of CpG sites is informative about methylation dynamics across tissue
differentiation. Finally, we identified genomic regions with clusters of highly
dynamic CpGs and present a likely functional example. Our work establishes a
framework for inference and modeling that is well-suited to DNA methylation
data, and our success suggests that other methods for analyzing DNA nucleotide
substitutions will also translate to the modeling of epigenetic phenomena.


  Myosin II isoforms with varying mechanochemistry and filament size interact
with filamentous actin (F-actin) networks to generate contractile forces in
cells. How their properties control force production in environments with
varying stiffness is poorly understood. Here, we incorporated literature values
for properties of myosin II isoforms into a cross-bridge model. Similar actin
gliding speeds and force-velocity curves expected from previous experiments
were observed. Motor force output on an elastic load was regulated by two
timescales--that of their attachment to F-actin, which varied sharply with the
ensemble size, motor duty ratio, and external load, and that of force build up,
which scaled with ensemble stall force, gliding speed, and load stiffness.
While such regulation did not require force-dependent kinetics, the myosin
catch bond produced positive feedback between attachment time and force to
trigger switch-like transitions from short attachments and small forces to high
force-generating runs at threshold parameter values. Parameters representing
skeletal muscle myosin, non-muscle myosin IIB, and non-muscle myosin IIA
revealed distinct regimes of behavior respectively: (1) large assemblies of
fast, low-duty ratio motors rapidly build stable forces over a large range of
environmental stiffness, (2) ensembles of slow, high-duty ratio motors serve as
high-affinity cross-links with force build-up times that exceed physiological
timescales, and (3) small assemblies of low-duty ratio motors operating at
intermediate speeds may respond sharply to changes in mechanical context--at
low forces or stiffness, they serve as low affinity cross-links but they can
transition to effective force production via the positive feedback mechanism
described above. These results reveal how myosin isoform properties may be
tuned to produce force and respond to mechanical cues in their environment.


  The present work presents a density-functional microscopic model of soft
biological tissue. The model was based on a prototype molecular structure from
experimentally resolved collagen peptide residues and water clusters and has
the objective to capture some well-known experimental features of soft tissues.
It was obtained the optimized geometry, binding and coupling energies and
dipole moments. The results concerning the stability of the confined water
clusters, the water-water and water-collagen interactions within the CLBM
framework were successfully correlated to some important trends observed
experimentally in inflammatory tissues.


  Inside cells, various cargos are transported by teams of molecular motors.
Intriguingly, the motors involved generally have opposite pulling directions,
and the resulting cargo dynamics is a biased stochastic motion. It is an open
question how the cell can control this bias. Here we develop a model which
takes explicitly into account the elastic coupling of the cargo with each
motor. We show that bias can be simply controlled or even reversed in a
counterintuitive manner via a change in the external force exerted on the cargo
or a variation of the ATP binding rate to motors. Furthermore, the
superdiffusive behavior found at short time scales indicates the emergence of
motor cooperation induced by cargo-mediated coupling.


  Background: Biological effects of extra-low-frequency (ELF) magnetic fields
(MFs) have lacked a credible mechanism of interaction between MFs and living
material. Objectives: To examine the effect of ELF-MFs on cancer cells.
Methods: Five cancer cell lines were exposed to ELF-MFs within the range of
0.025-5 microTesla, and the cells were examined for karyotype changes after 6
days. Results: All cancer cells lines lost chromosomes from MF exposure, with a
mostly flat dose-response. Constant MF exposures for three weeks allow a rising
return to the baseline, unperturbed karyotypes. From this point, small MF
increases or decreases are again capable of inducing karyotype contractions
(KCs). Our data suggest that the KCs are caused by MF interference with
mitochondria's adenosine triphosphate synthase (ATPS), compensated by the
action of adenosine monophosphate-activated protein kinase (AMPK). The effects
of MFs are similar to those of the ATPS inhibitor, oligomycin. They are
amplified by metformin, an AMPK stimulator, and attenuated by resistin, an AMPK
inhibitor. Over environmental MFs, KCs of various cancer cell lines show
exceptionally wide and flat dose-responses, except for those of erythroleukemia
cells, which display a progressive rise from 0.025 to 0.4 microTesla.
Conclusions: The biological effects of MFs are connected to an alteration in
the structure of water that impedes the flux of protons in ATPS channels. These
results may be environment ally important, in view of the central roles played
in human physiology by ATPS and AMPK, particularly in their links to diabetes,
cancer and longevity.


  Inter-individual differences in many behaviors are partly due to genetic
differences, but the identification of the genes and variants that influence
behavior remains challenging. Here, we studied an F2 intercross of two outbred
lines of rats selected for tame and aggressive behavior towards humans for more
than 64 generations. By using a mapping approach that is able to identify
genetic loci segregating within the lines, we identified four times more loci
influencing tameness and aggression than by an approach that assumes fixation
of causative alleles, suggesting that many causative loci were not driven to
fixation by the selection. We used RNA sequencing in 150 F2 animals to identify
hundreds of loci that influence brain gene expression. Several of these loci
colocalize with tameness loci and may reflect the same genetic variants.
Through analyses of correlations between allele effects on behavior and gene
expression, differential expression between the tame and aggressive rat
selection lines, and correlations between gene expression and tameness in F2
animals, we identify the genes Gltscr2, Lgi4, Zfp40 and Slc17a7 as candidate
contributors to the strikingly different behavior of the tame and aggressive
animals.


  Social dilemmas are central to human society. Depletion of natural resources,
climate protection, security of energy supply, and workplace collaborations are
all examples of social dilemmas. Since cooperative behaviour in a social
dilemma is individually costly, Nash equilibrium predicts that humans should
not cooperate. Yet experimental studies show that people do cooperate even in
anonymous one-shot interactions. In spite of the large number of participants
in many modern social dilemmas, little is known about the effect of group size
on cooperation. Does larger group size favour or prevent cooperation? We
address this problem both experimentally and theoretically. Experimentally, we
find that there is no general answer: it depends on the strategic situation.
Specifically, we find that larger groups are more cooperative in the Public
Goods game, but less cooperative in the N-person Prisoner's dilemma.
Theoretically, we show that this behaviour is not consistent with either the
Fehr & Schmidt model or (a one-parameter version of) the Charness & Rabin
model, but it is consistent with the cooperative equilibrium model introduced
by the second author.


  The active sites of enzymes consist of residues necessary for catalysis, and
structurally important non-catalytic residues that together maintain the
architecture and function of the active site. Examples of evolutionary
interactions between active site residue have been difficult to define and
experimentally validate due to a general intolerance to substitution. Here,
using computational methods to predict co-evolving residues, we identify a
network of positions consisting of two catalytic metal-binding residues and two
adjacent non- catalytic residues in LAGLIDADG family homing endonucleases
(LHEs). Distinct combinations of residues in the network map to clades of LHE
sub-families, with a striking distribution of the metal-binding Asp (D) and Glu
(E) residues. Mutation of these four positions in two LHEs, I-LtrI and I-OnuI,
indicate that combinations of residues tolerated are specific to each enzyme.
Kinetic analyses under single-turnover conditions indicated a substantial
defect in either the catalytic rate (kcat*) or the Michaelis constant (KM*) for
I- LtrI variants with networks predicted to be suboptimal. Enzymatic activity
could be restored to near wild-type levels with a single compensatory mutation
that recreated a network of optimally predicted residues. Our results
demonstrate that LHE activity is constrained by an evolutionary barrier of
residues with strong context-dependent effects. Creation of optimal co-evolving
active site networks is therefore an important consideration in engineering of
LHEs and other enzymes.


  We consider the Moran process in a graph called the "star" and obtain the
asymptotic expression for the fixation probability of a single mutant when the
size of the graph is large. The expression obtained corrects the previously
known expression announced in reference [E Lieberman, C Hauert, and MA Nowak.
Evolutionary dynamics on graphs. Nature, 433(7023):312-316, 2005] and further
studied in [M. Broom and J. Rychtar. An analysis of the fixation probability of
a mutant on special classes of non-directed graphs. Proc. R. Soc. A-Math. Phys.
Eng. Sci., 464(2098):2609-2627, 2008]. We also show that the star graph is an
accelerator of evolution, if the graph is large enough.


  Discovering causal genetic variants from large genetic association studies
poses many difficult challenges. Assessing which genetic markers are involved
in determining trait status is a computationally demanding task, especially in
the presence of gene-gene interactions. A non-parametric Bayesian approach in
the form of a Bayesian neural network is proposed for use in analyzing genetic
association studies. Demonstrations on synthetic and real data reveal they are
able to efficiently and accurately determine which variants are involved in
determining case-control status. Using graphics processing units (GPUs) the
time needed to build these models is decreased by several orders of magnitude.
In comparison with commonly used approaches for detecting genetic interactions,
Bayesian neural networks perform very well across a broad spectrum of possible
genetic relationships while having the computational efficiency needed to
handle large datasets.


  The fitness contribution of an allele at one genetic site may depend on
alleles at other sites, a phenomenon known as epistasis. Epistasis can
profoundly influence the process of evolution in populations under selection,
and can shape the course of protein evolution across divergent species. Whereas
epistasis between adaptive substitutions has been the subject of extensive
study, relatively little is known about epistasis under purifying selection.
Here we use mechanistic models of thermodynamic stability in a ligand-binding
protein to explore the structure of epistatic interactions between
substitutions that fix in protein sequences under purifying selection. We find
that the selection coefficients of mutations that are nearly-neutral when they
fix are highly contingent on the presence of preceding mutations. Conversely,
mutations that are nearly-neutral when they fix are subsequently entrenched due
to epistasis with later substitutions. Our evolutionary model includes
insertions and deletions, as well as point mutations, and so it allows us to
quantify epistasis within each of these classes of mutations, and also to study
the evolution of protein length. We find that protein length remains largely
constant over time, because indels are more deleterious than point mutations.
Our results imply that, even under purifying selection, protein sequence
evolution is highly contingent on history and so it cannot be predicted by the
phenotypic effects of mutations assayed in the wild-type sequence.


  Several neurological disorders are associated with the aggregation of
aberrant proteins, often localized in intracellular organelles such as the
endoplasmic reticulum. Here we study protein aggregation kinetics by mean-field
reactions and three dimensional Monte carlo simulations of diffusion-limited
aggregation of linear polymers in a confined space, representing the
endoplasmic reticulum. By tuning the rates of protein production and
degradation, we show that the system undergoes a non-equilibrium phase
transition from a physiological phase with little or no polymer accumulation to
a pathological phase characterized by persistent polymerization. A combination
of external factors accumulating during the lifetime of a patient can thus
slightly modify the phase transition control parameters, tipping the balance
from a long symptomless lag phase to an accelerated pathological development.
The model can be successfully used to interpret experimental data on
amyloid-\b{eta} clearance from the central nervous system.


  Determining phenotype from genetic data is a fundamental challenge. Influenza
A viruses undergo rapid antigenic drift and identification of emerging
antigenic variants is critical to the vaccine selection process. Using former
seasonal influenza A(H1N1) viruses, hemagglutinin sequence and corresponding
antigenic data were analyzed in combination with 3-D structural information. We
attributed variation in hemagglutination inhibition to individual amino acid
substitutions and quantified their antigenic impact, validating a subset
experimentally using reverse genetics. Substitutions identified as low-impact
were shown to be a critical component of influenza antigenic evolution and by
including these, as well as the high-impact substitutions often focused on, the
accuracy of predicting antigenic phenotypes of emerging viruses from genotype
was doubled. The ability to quantify the phenotypic impact of specific amino
acid substitutions should help refine techniques that predict the fitness and
evolutionary success of variant viruses, leading to stronger theoretical
foundations for selection of candidate vaccine viruses.


  Dynamic extinction colonisation models (also called contact processes) are
widely studied in epidemiology and in metapopulation theory. Contacts are
usually assumed to be possible only through a network of connected patches.
This network accounts for a spatial landscape or a social organisation of
interactions. Thanks to social network literature, heterogeneous networks of
contacts can be considered. A major issue is to assess the influence of the
network in the dynamic model. Most work with this common purpose uses
deterministic models or an approximation of a stochastic
Extinction-Colonisation model (sEC) which are relevant only for large networks.
When working with a limited size network, the induced stochasticity is
essential and has to be taken into account in the conclusions. Here, a rigorous
framework is proposed for limited size networks and the limitations of the
deterministic approximation are exhibited. This framework allows exact
computations when the number of patches is small. Otherwise, simulations are
used and enhanced by adapted simulation techniques when necessary. A
sensitivity analysis was conducted to compare four main topologies of networks
in contrasting settings to determine the role of the network. A challenging
case was studied in this context: seed exchange of crop species in the R\'eseau
Semences Paysannes (RSP), an emergent French farmers' organisation. A
stochastic Extinction-Colonisation model was used to characterize the
consequences of substantial changes in terms of RSP's social organisation on
the ability of the system to maintain crop varieties.


  Structural bioinformatics and van der Waals density functional theory are
combined to investigate the mechanochemical impact of a major class of
histone-DNA interactions, namely the formation of salt bridges between arginine
residues in histones and phosphate groups on the DNA backbone. Principal
component analysis reveals that the configurational fluctuations of the
sugar-phosphate backbone display sequence-specific variability, and clustering
of nucleosomal crystal structures identifies two major salt bridge
configurations: a monodentate form in which the arginine end-group guanidinium
only forms one hydrogen bond with the phosphate, and a bidentate form in which
it forms two. Density functional theory calculations highlight that the
combination of sequence, denticity and salt bridge positioning enable the
histones to tunably activate specific backbone deformations via mechanochemical
stress. The results suggest that selection for specific placements of van der
Waals contacts, with high-precision control of the spatial distribution of
intermolecular forces, may serve as an underlying evolutionary design principle
for the structure and function of nucleosomes, a conjecture that is
corroborated by previous experimental studies.


  Motivation: Genome-wide association studies (GWAS) have identified many loci
implicated in disease susceptibility. Integration of GWAS summary statistics (p
values) and functional genomic datasets should help to elucidate mechanisms.
Results: We describe the extension of a previously described non-parametric
method to test whether GWAS signals are enriched in functionally defined loci
to a situation where only GWAS p values are available. The approach is
implemented in VSEAMS, a freely available software pipeline. We use VSEAMS to
integrate functional gene sets defined via transcription factor knock down
experiments with GWAS results for type 1 diabetes and find variant set
enrichment in gene sets associated with IKZF3, BATF and ESRRA. IKZF3 lies in a
known T1D susceptibility region, whilst BATF and ESRRA overlap other immune
disease susceptibility regions, validating our approach and suggesting novel
avenues of research for type 1 diabetes. Availability and implementation:
VSEAMS is available for download http://github.com/ollyburren/vseams.


  Within the domain of molecular communications, researchers mimic the
techniques in nature to come up with alternative communication methods for
collaborating nanomachines. This work investigates the channel transfer
function for molecular communication via diffusion. In nature,
information-carrying molecules are generally absorbed by the target node via
receptors. Using the concentration function, without considering the absorption
process, as the channel transfer function implicitly assumes that the receiver
node does not affect the system. In this letter, we propose a solid analytical
formulation and analyze the signal metrics (attenuation and propagation delay)
for molecular communication via diffusion channel with an absorbing receiver in
a 3-D environment. The proposed model and the formulation match well with the
simulations without any normalization.


  Conformational barcodes tag functional sites of proteins, and are decoded by
interacting molecules transmitting the incoming signal. Conformational barcodes
are modified by all co-occurring allosteric events induced by
post-translational modifications, pathogen, drug binding, etc. We argue that
fuzziness (plasticity) of conformational barcodes may be increased by
disordered protein structures, by integrative plasticity of
multi-phosphorylation events, by increased intracellular water content
(decreased molecular crowding) and by increased action of molecular chaperones.
This leads to increased plasticity of signaling and cellular networks.
Increased plasticity is both substantiated by and inducing an increased noise
level. Using the versatile network dynamics tool, Turbine
(www.turbine.linkgroup.hu), here we show that the 10% noise level expected in
cellular systems shifts a cancer-related signaling network of human cells from
its proliferative attractors to its largest, apoptotic attractor representing
their health-preserving response in the carcinogen containing and tumor
suppressor deficient environment modeled in our study. Thus, fuzzy
conformational barcodes may not only make the cellular system more plastic, and
therefore more adaptable, but may also stabilize the complex system allowing
better access to its largest attractor.


  Territorial behaviour is widespread in the animal kingdom, with creatures
seeking to gain parts of space for their exclusive use. It arises through a
complicated interplay of many different behavioural features. Extracting and
quantifying the processes that give rise to territorial patterns requires both
mathematical models of movement and interaction mechanisms, together with
statistical techniques for rigorously extracting parameters from data. Here, we
give a brisk, pedagogical overview of the techniques so far developed to tackle
the problem of territory formation. We give some examples of what has already
been achieved using these techniques, together with pointers as to where we
believe the future lies in this area of study. This progress is a single
example of a major aim for 21st century science: to construct quantitatively
predictive theory for ecological systems.


  During organogenesis, developmental programs governed by Gene Regulatory
Networks (GRN) define the functionality, size and shape of the different
constituents of living organisms. Robustness, thus, is an essential
characteristic that GRNs need to fulfill in order to maintain viability and
reproducibility in a species. In the present work we analyze the robustness of
the patterning for the ocellar complex formation in Drosophila melanogaster
fly. We have systematically pruned the GRN that drives the development of this
visual system to obtain the minimum pathway able to satisfy this pattern. We
found that the mechanism underlying the patterning obeys to the dynamics of a
3-nodes network motif with a double negative feedback loop fed by a
morphogenetic gradient that triggers the inhibition in a French flag problem
fashion. A Boolean modeling of the GRN confirms robustness in the patterning
mechanism showing the same result for different network complexity levels.
Interestingly, the network provides a steady state solution in the interocellar
part of the patterning and an oscillatory regime in the ocelli. This
theoretical result predicts that the ocellar pattern may underlie oscillatory
dynamics in its genetic regulation.


  Cooperation is ubiquitous ranging from multicellular organisms to human
societies. Population structures indicating individuals' limited interaction
ranges are crucial to understand this issue. But it is still at large to what
extend multiple interactions involving nonlinearity in payoff play a role on
cooperation in structured populations. Here we show a rule, which determines
the emergence and stabilization of cooperation, under multiple discounted,
linear, and synergistic interactions. The rule is validated by simulations in
homogenous and heterogenous structured populations. We find that the more
neighbors there are the harder for cooperation to evolve for multiple
interactions with linearity and discounting. For synergistic scenario, however,
distinct from its pairwise counterpart, moderate number of neighbors can be the
worst, indicating that synergistic interactions work with strangers but not
with neighbors. Our results suggest that the combination of different factors
which promotes cooperation alone can be worse than that with every single
factor.


  A major theme in constraint-based modeling is unifying experimental data,
such as biochemical information about the reactions that can occur in a system
or the composition and localization of enzyme complexes, with highthroughput
data including expression data, metabolomics, or DNA sequencing. The desired
result is to increase predictive capability resulting in improved understanding
of metabolism. The approach typically employed when only gene (or protein)
intensities are available is the creation of tissue-specific models, which
reduces the available reactions in an organism model, and does not provide an
objective function for the estimation of fluxes, which is an important
limitation in many modeling applications. We develop a method, flux assignment
with LAD (least absolute deviation) convex objectives and normalization
(FALCON), that employs metabolic network reconstructions along with expression
data to estimate fluxes. In order to use such a method, accurate measures of
enzyme complex abundance are needed, so we first present a new algorithm that
addresses quantification of complex abundance. Our extensions to prior
techniques include the capability to work with large models and significantly
improved run-time performance even for smaller models, an improved analysis of
enzyme complex formation logic, the ability to handle very large enzyme complex
rules that may incorporate multiple isoforms, and depending on the model
constraints, either maintained or significantly improved correlation with
experimentally measured fluxes. FALCON has been implemented in MATLAB and ATS,
and can be downloaded from: https://github.com/bbarker/FALCON. ATS is not
required to compile the software, as intermediate C source code is available,
and binaries are provided for Linux x86-64 systems. FALCON requires use of the
COBRA Toolbox, also implemented in MATLAB.


  Plant-pollinator associations are often seen as purely mutualistic, while in
reality they can be more complex. Indeed they may also display a diverse array
of antagonistic interactions, such as competition and victim--exploiter
interactions. In some cases mutualistic and antagonistic interactions are
carried-out by the same species but at different life-stages. As a consequence,
population structure affects the balance of inter-specific associations, a
topic that is receiving increased attention. In this paper, we developed a
model that captures the basic features of the interaction between a flowering
plant and an insect with a larval stage that feeds on the plant's vegetative
tissues (e.g. leaves) and an adult pollinator stage. Our model is able to
display a rich set of dynamics, the most remarkable of which involves
victim--exploiter oscillations that allow plants to attain abundances above
their carrying capacities, and the periodic alternation between states
dominated by mutualism or antagonism. Our study indicates that changes in the
insect's life cycle can modify the balance between mutualism and antagonism,
causing important qualitative changes in the interaction dynamics. These
changes in the life cycle could be caused by a variety of external drivers,
such as temperature, plant nutrients, pesticides and changes in the diet of
adult pollinators.
  Abstract Keywords: mutualism, pollination, herbivory, insects,
stage-structure, oscillations


  Mathematical models play an increasingly important role in the interpretation
of biological experiments. Studies often present a model that generates the
observations, connecting hypothesized process to an observed pattern. Such
generative models confirm the plausibility of an explanation and make testable
hypotheses for further experiments. However, studies rarely consider the broad
family of alternative models that match the same observed pattern. The
symmetries that define the broad class of matching models are in fact the only
aspects of information truly revealed by observed pattern. Commonly observed
patterns derive from simple underlying symmetries. This article illustrates the
problem by showing the symmetry associated with the observed rate of increase
in fitness in a constant environment. That underlying symmetry reveals how each
particular generative model defines a single example within the broad class of
matching models. Further progress on the relation between pattern and process
requires deeper consideration of the underlying symmetries.


  The efficient recognition of pathogens by the adaptive immune system relies
on the diversity of receptors displayed at the surface of immune cells. T-cell
receptor diversity results from an initial random DNA editing process, called
VDJ recombination, followed by functional selection of cells according to the
interaction of their surface receptors with self and foreign antigenic
peptides. To quantify the effect of selection on the highly variable elements
of the receptor, we apply a probabilistic maximum likelihood approach to the
analysis of high-throughput sequence data from the $\beta$-chain of human
T-cell receptors. We quantify selection factors for V and J gene choice, and
for the length and amino-acid composition of the variable region. Our approach
is necessary to disentangle the effects of selection from biases inherent in
the recombination process. Inferred selection factors differ little between
donors, or between naive and memory repertoires. The number of sequences shared
between donors is well-predicted by the model, indicating a purely stochastic
origin of such "public" sequences. We find a significant correlation between
biases induced by VDJ recombination and our inferred selection factors,
together with a reduction of diversity during selection. Both effects suggest
that natural selection acting on the recombination process has anticipated the
selection pressures experienced during somatic evolution.


  To support comparative genomics, population genetics, and medical genetics,
we propose that a reference genome should come with a scheme for mapping each
base in any DNA string to a position in that reference genome. We refer to a
collection of one or more reference genomes and a scheme for mapping to their
positions as a reference structure. Here we describe the desirable properties
of reference structures and give examples. To account for natural genetic
variation, we consider the more general case in which a reference genome is
represented by a graph rather than a set of phased chromosomes; the latter is
treated as a special case.


  Elucidation of possible pathways between folded (native) and unfolded states
of a protein is a challenging task, as the intermediates are often hard to
detect. Here we alter the solvent environment in a controlled manner by
choosing two different co-solvents of water, urea and dimethyl sulphoxide
(DMSO), and study unfolding of four different proteins to understand the
respective sequence of melting by computer simulation methods. We indeed find
interesting differences in the sequence of melting of alpha-helices and
beta-sheets in these two solvents. For example, at 8M urea solution, beta-sheet
parts of a protein is found to unfold preferentially, followed by the unfolding
of alpha helices. In contrast, 8M DMSO solution unfolds alpha helices first,
followed by the separation of beta-sheets for majority of proteins. Sequence of
unfolding events in four different alpha/beta proteins and also in chicken
villin head piece (HP-36) both in urea and DMSO solution demonstrate that the
unfolding pathways are determined jointly by relative exposure of polar and
non-polar residues of a protein and the mode of molecular action of a solvent
on that protein.


  Metabolic fluxes in cells are governed by physical, physiological, and
economic principles. Here I assume an optimal allocation of enzyme resources
and postulate a general principle for metabolism: each enzyme must convert less
valuable into more valuable metabolites to justify its own cost. The "values",
called economic potentials, describe the individual contributions of
metabolites to cell fitness. Local value production implies that the cost of an
enzyme must be balanced by a benefit, given by the economic potential
difference the catalysed reaction multiplied by the flux. Flux profiles that
satisfy this principle - i.e. for which consistent potentials can be found -
are called economical. Economical fluxes must lead from lower to higher
economic potentials, so certain flux cycles are incompatible with any choice of
economic potentials and can be excluded. To obtain economical flux profiles,
non-beneficial local patterns, called futile motifs, can be systematically
removed from a given flux distribution. The principle of local value production
resembles thermodynamic principles and complements them in models. Here I
describe a modelling framework called Value Balance Analysis (VBA) that uses
the two principles and yields the same solution as enzyme cost minimisation (in
kinetic models) and flux cost minimisation (in FBA). Given an economical flux
distribution, kinetic models in enzyme-optimal states and with these fluxes can
be constructed systematically. VBA justifies the principle of minimal fluxes
and the exclusion of futile cycles, predicts enzymes that could be plausible
targets for regulation, provides criteria for the usage of enzymes and
pathways, and explains the choice between high-yield and low-yield flux modes.


  Motivated by several recent experimental observations that vitamin-D could
interact with antigen presenting cells (APCs) and T-lymphocyte cells (T-cells)
to promote and to regulate different stages of immune response, we developed a
coarse grained kinetic model in an attempt to quantify the role of vitamin-D in
immunomodulatory responses. Our kinetic model, developed using the ideas of
chemical network theory, leads to a system of nine coupled equations that we
solve both by direct and by stochastic (Gillespie) methods. Both the analyses
consistently provide detail information on the dependence of immune response to
the variation of critical rate parameters. We find that although vitamin-D
plays a negligible role in the initial immune response, it exerts a profound
influence in the long term, especially in helping the system to achieve a new,
stable steady state. The study explores the role of vitamin-D in preserving an
observed bistability in the phase diagram (spanned by system parameters) of
immune regulation, thus allowing the response to tolerate a wide range of
pathogenic stimulation which could help in resisting autoimmune diseases. We
also study how vitamin-D affects the time dependent population of dendritic
cells that connect between innate and adaptive immune responses. Variations in
dose dependent response in anti-inflammatory and pro-inflammatory T-cell
populations to vitamin-D correlate well with recent experimental results. Our
kinetic model allows for an estimation of the range of optimum level of
vitamin-D required for smooth functioning of the immune system and for control
of both hyper-regulation and inflammation. Most importantly, the present study
reveals that an overdose or toxic level of vitamin-D or any steroid analogue
could give rise to too large a tolerant response, leading to an inefficacy in
adaptive immune function.


  The relatively low concentration of pathogen indicators, such Salmonella, in
composting sometimes causes a problem with detection when using the
conventional techniques. The presence of viable but non-culturable (VBNC)
organisms is also a potential problem with Salmonella detection when using
conventional techniques. In this study, the molecular approach for organism
recognition, known as Polymerase Chain Reaction (PCR), was used for
characterisation the Salmonella spp. used as inoculums in composting. This
study also provides a better understand about selecting the suitable primer for
Salmonella spp. The specificity of the probes and primers at the species level
were verified by performing NCBI-BLAST2 (Basic Local Alignment Search Tool).
Primers that target the invA gene for Salmonella were selected which produce
fragment lengths around 285bp (shown in Figure 1). The Salmonella spp.
solutions were tested and proved contained the sequences of invA gene by using
several steps of molecular techniques before used it as inoculums in
composting. The laboratory scale batch composting reactors were used to examine
the inactivation of Salmonella spp. The inoculate solution of Salmonella was
prepared by culturing Salmonella enteritidis (ATCC13076) in tryptone broth
solution for 24 hours before adding it directly into the compost material. The
reduction rate of Salmonella spp. was enumerated by conventional method
accordingly to the standard compost guidelines (Figure 2). The laboratory scale
study showed that after composting for 8 days the numbers of Salmonella spp.
were below the limits in the UK compost standard (PAS 100) which requires the
compost to be free of Salmonella spp.


  The metabolic state of a cell, comprising fluxes, metabolite concentrations
and enzyme levels, is shaped by a compromise between metabolic benefit and
enzyme cost. This hypothesis and its consequences can be studied by
computational models and using a theory of metabolic value. In optimal
metabolic states, any increase of an enzyme level must improve the metabolic
performance to justify its own cost, so each active enzyme must contribute to
the cell's benefit by producing valuable products. This principle of value
production leads to variation rules that relate metabolic fluxes and reaction
elasticities to enzyme costs. Metabolic value theory provides a language to
describe this. It postulates a balance of local values, which I derive here
from concepts of metabolic control theory. Economic state variables, called
economic potentials and loads, describe how metabolites, reactions, and enzymes
contribute to metabolic performance. Economic potentials describe the indirect
value of metabolite production, while economic loads describe the indirect
value of metabolite concentrations. These economic variables, and others, are
linked by local balance equations. These laws for optimal metabolic states
define conditions for metabolic fluxes that hold for a wide range of rate laws.
To produce metabolic value, fluxes run from lower to higher economic
potentials, must be free of futile cycles, and satisfy a principle of minimal
weighted fluxes. Given an economical flux mode, one can systematically
construct kinetic models in which all enzymes have positive effects on
metabolic performance.


  Modeling gene regulatory networks (GRNs) is an important topic in systems
biology. Although there has been much work focusing on various specific
systems, the generic behavior of GRNs with continuous variables is still
elusive. In particular, it is not clear typically how attractors partition
among the three types of orbits: steady state, periodic and chaotic, and how
the dynamical properties change with network's topological characteristics. In
this work, we first investigated these questions in random GRNs with different
network sizes, connectivity, fraction of inhibitory links and transcription
regulation rules. Then we searched for the core motifs that govern the dynamic
behavior of large GRNs. We show that the stability of a random GRN is typically
governed by a few embedding motifs of small sizes, and therefore can in general
be understood in the context of these short motifs. Our results provide
insights for the study and design of genetic networks.


  Recent methodological advances are enabling better examination of speciation
and extinction processes and patterns. A major open question is the origin of
large discrepancies in species number between groups of the same age. Existing
frameworks to model this diversity either focus on changes between lineages,
neglecting global effects such as mass extinctions, or focus on changes over
time which would affect all lineages. Yet it seems probable that both lineages
differences and mass extinctions affect the same groups. Here we used
simulations to test the performance of two widely used methods, under complex
scenarios. We report good performances, although with a tendency to
over-predict events when increasing the complexity of the scenario. Overall, we
find that lineage shifts are better detected than mass extinctions. This work
has significance for assessing the methods currently used for estimating
changes in diversification using phylogenies and developing new tests.


  In the applications of Boolean networks to modeling biological systems, an
important computational problem is the detection of the fixed points of these
networks. This is an NP-complete problem in general. There have been various
attempts to develop algorithms to address the computation need for large size
Boolean networks. The existing methods are usually based on known algorithms
and thus limited to the situations where these known algorithms can apply. In
this paper, we propose a novel approach to this problem. We show that any
system of Boolean equations is equivalent to one Boolean equation, and thus it
is possible to divide the polynomial equation system which defines the fixed
points of a Boolean network into subsystems that can be solved easily. After
solving these subsystems and thus reducing the number of states, we can combine
the solutions to obtain all fixed points of the given network. This approach
does not depend on other algorithms and it is straightforward and easy to
implement. We show that our method can handle large size Boolean networks, and
demonstrate its effectiveness by using MAPLE to compute the fixed points of
Boolean networks with hundreds of nodes and thousands of interactions.


  To model biological systems using networks, it is desirable to allow more
than two levels of expression for the nodes and to allow the introduction of
parameters. Various modeling and simulation methods addressing these needs
using Boolean models, both synchronous and asynchronous, have been proposed in
the literature. However, analytical study of these more general Boolean
networks models is lagging. This paper aims to develop a concise theory for
these different Boolean logic based modeling methods. Boolean models for
networks where each node can have more than two levels of expression and
Boolean models with parameters are defined algebraically with examples
provided. Certain classes of random asynchronous Boolean networks and
deterministic moduli asynchronous Boolean networks are investigated in detail
using the setting introduced in this paper. The derived theorems provide a
clear picture for the attractor structures of these asynchronous Boolean
networks.


  Multisite phosphorylation plays an important role in intracellular signaling.
There has been much recent work aimed at understanding the dynamics of such
systems when the phosphorylation/dephosphorylation mechanism is distributive,
that is, when the binding of a substrate and an enzyme molecule results in
addition or removal of a single phosphate group and repeated binding therefore
is required for multisite phosphorylation. In particular, such systems admit
bistability. Here we analyze a different class of multisite systems, in which
the binding of a substrate and an enzyme molecule results in addition or
removal of phosphate groups at all phosphorylation sites. That is, we consider
systems in which the mechanism is processive, rather than distributive. We show
that in contrast with distributive systems, processive systems modeled with
mass-action kinetics do not admit bistability and, moreover, exhibit rigid
dynamics: each invariant set contains a unique equilibrium, which is a global
attractor. Additionally, we obtain a monomial parametrization of the steady
states. Our proofs rely on a technique of Johnston for using "translated"
networks to study systems with "toric steady states", recently given sign
conditions for injectivity of polynomial maps, and a result from monotone
systems theory due to Angeli and Sontag.


  We study a discrete host-parasitoid system where the host population follows
the classical Ricker functional form and is also subject to Allee effects. We
determine basins of attraction of the local attractors of the single population
model when the host intrinsic growth rate is not large. In this situation,
existence and local stability of the interior steady states for the
host-parasitoid interaction are completely analyzed. If the host's intrinsic
growth rate is large, then the interaction may support multiple interior steady
states. Linear stability of these steady states is provided.


  Over the last years the vertex enumeration problem of polyhedra has seen a
revival in the study of metabolic networks, which increased the demand for
efficient vertex enumeration algorithms for high-dimensional polyhedra given by
inequalities. It is a famous and long standing open question in polyhedral
theory and computational geometry whether the vertices of a polytope (bounded
polyhedron), described by a set of linear constraints, can be enumerated in
total polynomial time. In this paper we apply the concept of
branch-decomposition to the vertex enumeration problem of polyhedra $P = \{x :
Ax = b, x \geq 0\}$. For this purpose, we introduce the concept of $k$-module
and show how it relates to the separators of the linear matroid generated by
the columns of $A$. We then use this to present a total polynomial time
algorithm for polytopes $P$ for which the branch-width of the linear matroid
generated by $A$ is bounded by a constant $k$.


  The concept of positional information is central to our understanding of how
cells in a multicellular structure determine their developmental fates.
Nevertheless, positional information has neither been defined mathematically
nor quantified in a principled way. Here we provide an information-theoretic
definition in the context of developmental gene expression patterns and examine
which features of expression patterns increase or decrease positional
information. We connect positional information with the concept of positional
error and develop tools to directly measure information and error from
experimental data. We illustrate our framework for the case of gap gene
expression patterns in the early Drosophila embryo and show how information
that is distributed among only four genes is sufficient to determine
developmental fates with single cell resolution. Our approach can be
generalized to a variety of different model systems; procedures and examples
are discussed in detail.


  Cultural transmission models are coming to the fore in explaining increases
in the Paleolithic toolkit richness and diversity. During the later
Paleolithic, technologies increase not only in terms of diversity but also in
their complexity and interdependence. As Mesoudi and O'Brien (2008) have shown,
selection broadly favors social learning of information that is hierarchical
and structured, and multiple studies have demonstrated that teaching within a
social learning environment can increase fitness. We believe that teaching also
provides the scaffolding for transmission of more complex cultural traits.
Here, we introduce an extension of the Axelrod (1997} model of cultural
differentiation in which traits have prerequisite relationships, and where
social learning is dependent upon the ordering of those prerequisites. We
examine the resulting structure of cultural repertoires as learning
environments range from largely unstructured imitation, to structured teaching
of necessary prerequisites, and we find that in combination with individual
learning and innovation, high probabilities of teaching prerequisites leads to
richer cultural repertoires. Our results point to ways in which we can build
more comprehensive explanations of the archaeological record of the Paleolithic
as well as other cases of technological change.


  Assembly of protein complexes like virus shells, the centriole, the nuclear
pore complex or the actin cytoskeleton is strongly determined by their spatial
structure. Moreover it is becoming increasingly clear that the reversible
nature of protein assembly is also an essential element for their biological
function. Here we introduce a computational approach for the Brownian dynamics
of patchy particles with anisotropic assemblies and fully reversible reactions.
Different particles stochastically associate and dissociate with microscopic
reaction rates depending on their relative spatial positions. The translational
and rotational diffusive properties of all protein complexes are evaluated
on-the-fly. Because we focus on reversible assembly, we introduce a scheme
which ensures detailed balance for patchy particles. We then show how the
macroscopic rates follow from the microscopic ones. As an instructive example,
we study the assembly of a pentameric ring structure, for which we find
excellent agreement between simulation results and a macroscopic kinetic
description without any adjustable parameters. This demonstrates that our
approach correctly accounts for both the diffusive and reactive processes
involved in protein assembly.


  An increasing number of studies have demonstrated correlations between
climate trends and body size change of organisms. In many cases, climate might
be expected to influence body size by altering thermoregulation, energetics or
food availability. However, observed body size changes can result from
ecological processes (e.g., growth, selection, population dynamics), yet may
also be due to imperfect observation. We used two extensive datasets to
evaluate alternative hypotheses for recently reported changes in the observed
body size of plethodontid salamanders. We found that mean adult body size of
salamanders can be highly sensitive to survey conditions, particularly
rainfall, such that smaller individuals are more likely to be sampled under dry
conditions. This systematic bias in the detection of individuals across a range
in body size would result in a signature of body size reduction in relation to
reported climate trends when it is simply observation error. We also identify
considerable variability in body size distributions among years that shows
correspondence with rainfall. This suggests that annual variation in growth or
shifting population age structure could also result in a correlation between
climate and mean body size. Finally, our study demonstrates that measures of
mean adult body size can be highly variable among surveys and that large sample
sizes may be required to make reliable inferences. Identifying the effects of
climate change is critical for future research in ecology and conservation.
Researchers should be aware that observed body size changes in some organisms
may be a result of true ecological processes or systematic bias due to
non-random sampling of populations. Ultimately, the credibility of ecological
research related to climate change depends on researchers demonstrating a
thorough consideration of alternative hypotheses for observed changes in
species.


  When there is conspicuous under-exploitation of a limited resource, it is
worth asking what mechanisms allow presumably valuable resources to be left
unused? Evolutionary biologists have generated a wide range of hypotheses to
explain this, ranging from interdemic group selection to selfishly prudent
individual restraint. We consider a situation in which, in spite of high
intraspecific competition, individuals leave most of a key resource
unexploited. The parasitic wasp that does this finds virtually all host egg
clusters in a landscape, but parasitizes only about a third of the eggs in
each, and then leaves a deterrent mark around the cluster. We first test, and
reject, a series of system-specific simple constraints that might limit full
host exploitation, such as asynchronous maturation of host eggs. We then
consider classical hypotheses for the evolution of restraint. Prudent predation
and bet-hedging fail as explanations because the wasp lives as a large
well-mixed population. Also, we find no individual benefits to the parasitoid
of developing in a sparsely parasitized host nest. However an optimal foraging
model, including empirically measured costs of superparasitism and
hyperparasitism, can explain through individual selection, both the
consistently low rate of parasitism and marking.


  As increasingly complex hypothesis-testing scenarios are considered in many
scientific fields, analytic derivation of null distributions is often out of
reach. To the rescue comes Monte Carlo testing, which may appear deceptively
simple: as long as you can sample test statistics under the null hypothesis,
the $p$-value is just the proportion of sampled test statistics that exceed the
observed test statistic. Sampling test statistics is often simple once you have
a Monte Carlo null model for your data, and defining some form of randomization
procedure is also, in many cases, relatively straightforward. However, there
may be several possible choices of a randomization null model for the data and
no clear-cut criteria for choosing among them. Obviously, different null models
may lead to very different $p$-values, and a very low $p$-value may thus occur
due to the inadequacy of the chosen null model. It is preferable to use
assumptions about the underlying random data generation process to guide
selection of a null model. In many cases, we may order the null models by
increasing preservation of the data characteristics, and we argue in this paper
that this ordering in most cases gives increasing $p$-values, that is, lower
significance. We denote this as the null complexity principle. The principle
gives a better understanding of the different null models and may guide in the
choice between the different models.


  A major part of membrane function is conducted by proteins, both integral and
peripheral. Peripheral membrane proteins temporarily adhere to biological
membranes, either to the lipid bilayer or to integral membrane proteins with
non-covalent interactions. The aim of this study was to construct and analyze
the interactions of the human plasma membrane peripheral proteins (peripherome
hereinafter). For this purpose, we collected a dataset of peripheral proteins
of the human plasma membrane. We also collected a dataset of experimentally
verified interactions for these proteins. The interaction network created from
this dataset has been visualized using Cytoscape. We grouped the proteins based
on their subcellular location and clustered them using the MCL algorithm in
order to detect functional modules. Moreover, functional and graph theory based
analyses have been performed to assess biological features of the network.
Interaction data with drug molecules show that ~10% of peripheral membrane
proteins are targets for approved drugs, suggesting their potential
implications in disease. In conclusion, we reveal novel features and properties
regarding the protein-protein interaction network created by peripheral
proteins of the human plasma membrane.


  Enhanced dynamical fluctuations of RNAs, facilitated by a network of water
molecules with strong interactions with RNA, are suspected to be critical in
their ability to respond to a variety of cellular signals. Using atomically
detailed molecular dynamics simulations at various temperatures of purine
(adenine)- and preQ$_1$ sensing riboswitch aptamers, we show that water
molecules in the vicinity of RNAs undergo complex dynamics depending on the
local structures of the RNAs. The overall lifetimes of hydrogen bonds (HBs) of
surface bound waters are more than at least 1-2 orders of magnitude longer than
bulk water. Slow hydration dynamics, revealed in non-Arrhenius behavior of the
relaxation time, arises from high activation barriers to break water hydrogen
bonds with a nucleotide and by reduced diffusion of water. The relaxation
kinetics at specific locations in the two RNAs show a broad spectrum of time
scales reminiscent of glass-like behavior, suggesting that the hydration
dynamics is highly heterogeneous. Both RNAs undergo dynamic transition at $T =
T_D \gtrsim 200$ K as assessed by the mean square fluctuation of hydrogen atoms
$\langle x^2\rangle$, which undergoes an abrupt harmonic-to-anharmonic
transition at $T_D$. The near universal value of $T_D$ found for these RNAs and
previously for tRNA is strongly correlated with changes in hydration dynamics
as $T$ is altered. Hierarchical dynamics of waters associated with the RNA
surface, revealed in the motions of distinct classes of water with
well-separated time scales, reflects the heterogeneous local environment on the
molecular surface of RNA. At low temperatures slow water dynamics predominates
over structural transitions. Our study demonstrates that the complex interplay
of dynamics between water and local environment in the RNA structures could be
a key determinant of the functional activities of RNA.


  Cooperation is a widespread natural phenomenon yet current evolutionary
thinking is dominated by the paradigm of selfish competition. Recent advanced
in many fronts of Biology and Non-linear Physics are helping to bring
cooperation to its proper place. In this contribution, the most important
controversies and open research avenues in the field of social evolution are
reviewed. It is argued that a novel theory of social evolution must integrate
the concepts of the science of Complex Systems with those of the Darwinian
tradition. Current gene-centric approaches should be reviewed and com-
plemented with evidence from multilevel phenomena (group selection), the
constrains given by the non-linear nature of biological dynamical systems and
the emergent nature of dissipative phenomena.


  The seasonality of respiratory diseases (common cold, influenza, etc.) is a
well-known phenomenon studied from ancient times. The development of predictive
models is still not only an actual unsolved problem of mathematical
epidemiology but also is very important for the safety of public health. Here
we show that SIRS (Susceptible-Infected-Recovered-Susceptible) model accurately
enough reproduces real curves of flu activity. It contains variable reaction
rate, which is a function of mean daily temperature. The proposed alternation
of variables represents SIRS equations as the second-order ODE with an outer
excitation. It reveals an origin of such predictive efficiency and explains
analytically the 1:1 dynamical resonance, which is known as a crucial property
of epidemic behavior. Our work opens the perspectives for the development of
instant short-time prediction of a normal level of flu activity based on the
weather forecast, and allow to estimate a current epidemic level more
precisely. The latter fact is based on the explicit difference between the
expected weather-based activity and instant anomalies.


  Analyzing qualitative behaviors of biochemical reactions using its associated
network structure has proven useful in diverse branches of biology. As an
extension of our previous work, we introduce a graph-based framework to
calculate steady state solutions of biochemical reaction networks with
synthesis and degradation. Our approach is based on a labeled directed graph
$G$ and the associated system of linear non-homogeneous differential equations
with first order degradation and zeroth order synthesis. We also present a
theorem which provides necessary and sufficient conditions for the dynamics to
engender a unique stable steady state.
  Although the dynamics are linear, one can apply this framework to nonlinear
systems by encoding nonlinearity into the edge labels. We answer open question
from our previous work concerning the non-positiveness of the elements in the
inverse of a perturbed Laplacian matrix. Moreover, we provide a graph
theoretical framework for the computation of the inverse of a such matrix. This
also completes our previous framework and makes it purely graph theoretical.
Lately, we demonstrate the utility of this framework by applying it to a
mathematical model of insulin secretion through ion channels and glucose
metabolism in pancreatic $\beta$-cells.


  Diversity measurement underpins the study of biological systems, but measures
used vary across disciplines. Despite their common use and broad utility, no
unified framework has emerged for measuring, comparing and partitioning
diversity. The introduction of information theory into diversity measurement
has laid the foundations, but the framework is incomplete without the ability
to partition diversity, which is central to fundamental questions across the
life sciences: How do we prioritise communities for conservation? How do we
identify reservoirs and sources of pathogenic organisms? How do we measure
ecological disturbance arising from climate change?
  The lack of a common framework means that diversity measures from different
fields have conflicting fundamental properties, allowing conclusions reached to
depend on the measure chosen. This conflict is unnecessary and unhelpful. A
mathematically consistent framework would transform disparate fields by
delivering scientific insights in a common language. It would also allow the
transfer of theoretical and practical developments between fields.
  We meet this need, providing a versatile unified framework for partitioning
biological diversity. It encompasses any kind of similarity between
individuals, from functional to genetic, allowing comparisons between
qualitatively different kinds of diversity. Where existing partitioning
measures aggregate information across the whole population, our approach
permits the direct comparison of subcommunities, allowing us to pinpoint
distinct, diverse or representative subcommunities and investigate population
substructure. The framework is provided as a ready-to-use R package to easily
test our approach.


  Host-to-host variability with respect to interactions between microorganisms
and multicellular hosts are commonly observed in infection and in homeostasis.
However, the majority of mechanistic models used in analyzing
host-microorganism relationships, as well as most of the ecological theories
proposed to explain co-evolution of host and microbes, are based on averages
across a host population. By assuming that observed variations are random and
independent, these models overlook the role of inter-host differences. Here we
analyze mechanisms underlying host-to-host variations, using the
well-characterized experimental infection model of polymicrobial otitis media
(OM) in chinchillas, in combination with population dynamic models and a
Maximum Entropy (MaxEnt) based inference scheme. We find that the nature of the
interactions among bacterial species critically regulates host-to-host
variations of these interactions. Surprisingly, seemingly unrelated phenomena,
such as the efficiency of individual bacterial species in utilizing nutrients
for growth and the microbe-specific host immune response, can become
interdependent in a host population. The latter finding suggests a potential
mechanism that could lead to selection of specific strains of bacterial species
during the coevolution of the host immune response and the bacterial species.


  We consider a population living in a patchy environment that varies
stochastically in space and time. The population is composed of two morphs
(that is, individuals of the same species with different genotypes). In terms
of survival and reproductive success, the associated phenotypes differ only in
their habitat selection strategies. We compute invasion rates corresponding to
the rates at which the abundance of an initially rare morph increases in the
presence of the other morph established at equilibrium. If both morphs have
positive invasion rates when rare, then there is an equilibrium distribution
such that the two morphs coexist; that is, there is a protected polymorphism
for habitat selection. Alternatively, if one morph has a negative invasion rate
when rare, then it is asymptotically displaced by the other morph under all
initial conditions where both morphs are present. We refine the
characterization of an evolutionary stable strategy for habitat selection from
[Schreiber, 2012] in a mathematically rigorous manner. We provide a necessary
and sufficient condition for the existence of an ESS that uses all patches and
determine when using a single patch is an ESS. We also provide an explicit
formula for the ESS when there are two habitat types. We show that adding
environmental stochasticity results in an ESS that, when compared to the ESS
for the corresponding model without stochasticity, spends less time in patches
with larger carrying capacities and possibly makes use of sink patches, thereby
practicing a spatial form of bet hedging.


  Recent literatures reported blue-green emission from amyloid fibril as
exclusive signature of fibril formation. This unusual visible luminescence is
regularly used to monitor fibril growth. Blue-green emission has also been
observed in crystalline protein and in solution. However, the origin of this
emission is not known exactly. Our spectroscopic study of serum proteins
reveals that the blue-green emission is a property of protein monomer.
Evidences suggest that semiconductor-like band structure of proteins with the
optical band-gap in the visible region is possibly the origin of this
phenomenon. We show here that the band structure of proteins is primarily the
result of electron delocalization through the peptide chain, rather than
through the hydrogen bond network in secondary structure.


  Unsupervised word segmentation methods were applied to analyze the protein
sequence. Protein sequences, such as 'MTMDKSELVQKA...', were used as input to
these methods. Segmented 'protein word' sequences, such as 'MTM DKSE LVQKA',
were then obtained. We compare the 'protein words' produced by unsupervised
segmentation and the protein secondary structure segmentation. An interesting
finding is that the unsupervised word segmentation is more efficient than
secondary structure segmentation in expressing information. Our experiment also
suggests there may be some 'protein ruins' in current noncoding regions.


  Street demonstrations occur across the world. In Rio de Janeiro, June/July
2013, they reach beyond one million people. A wrathful reader of \textit{O
Globo}, leading newspaper in the same city, published a letter \cite{OGlobo}
where many social questions are stated and answered Yes or No. These million
people of street demonstrations share opinion consensus about a similar set of
social issues. But they did not reach this consensus within such a huge
numbered meetings. Earlier, they have met in diverse small groups where some of
them could be convinced to change mind by other few fellows. Suddenly, a
macroscopic consensus emerges. Many other big manifestations are widespread all
over the world in recent times, and are supposed to remain in the future. The
interesting questions are: 1) How a binary-option opinion distributed among
some population evolves in time, through local changes occurred within
small-group meetings? and 2) Is there some natural selection rule acting upon?
Here, we address these questions through an agent-based model.


  We consider the problem of estimating the evolutionary history of a set of
species (phylogeny or species tree) from several genes. It is known that the
evolutionary history of individual genes (gene trees) might be topologically
distinct from each other and from the underlying species tree, possibly
confounding phylogenetic analysis. A further complication in practice is that
one has to estimate gene trees from molecular sequences of finite length. We
provide the first full data-requirement analysis of a species tree
reconstruction method that takes into account estimation errors at the gene
level. Under that criterion, we also devise a novel reconstruction algorithm
that provably improves over all previous methods in a regime of interest.


  Cooperation in one-shot anonymous interactions is a widely documented aspect
of human behaviour. Here we shed light on the motivations behind this behaviour
by experimentally exploring cooperation in a one-shot continuous-strategy
Prisoner's Dilemma (i.e. one-shot two-player Public Goods Game). We examine the
distribution of cooperation amounts, and how that distribution varies based on
the benefit-to-cost ratio of cooperation (b/c). Interestingly, we find a
trimodal distribution at all b/c values investigated. Increasing b/c decreases
the fraction of participants engaging in zero cooperation and increases the
fraction engaging in maximal cooperation, suggesting a role for efficiency
concerns. However, a substantial fraction of participants consistently engage
in 50% cooperation regardless of b/c. The presence of these persistent 50%
cooperators is surprising, and not easily explained by standard models of
social preferences. We present evidence that this behaviour is a result of
social preferences guided by simple decision heuristics, rather than the
rational examination of payoffs assumed by most social preference models. We
also find a strong correlation between play in the Prisoner's Dilemma and in a
subsequent Dictator Game, confirming previous findings suggesting a common
prosocial motivation underlying altruism and cooperation.


  In this review we summarize our recent efforts in trying to understand the
role of heterogeneity in cancer progression by using neural networks to
characterise different aspects of the mapping from a cancer cells genotype and
environment to its phenotype. Our central premise is that cancer is an evolving
system subject to mutation and selection, and the primary conduit for these
processes to occur is the cancer cell whose behaviour is regulated on multiple
biological scales. The selection pressure is mainly driven by the
microenvironment that the tumour is growing in and this acts directly upon the
cell phenotype. In turn, the phenotype is driven by the intracellular pathways
that are regulated by the genotype. Integrating all of these processes is a
massive undertaking and requires bridging many biological scales (i.e.
genotype, pathway, phenotype and environment) that we will only scratch the
surface of in this review. We will focus on models that use neural networks as
a means of connecting these different biological scales, since they allow us to
easily create heterogeneity for selection to act upon and importantly this
heterogeneity can be implemented at different biological scales. More
specifically, we consider three different neural networks that bridge different
aspects of these scales and the dialogue with the micro-environment, (i) the
impact of the micro-environment on evolutionary dynamics, (ii) the mapping from
genotype to phenotype under drug-induced perturbations and (iii) pathway
activity in both normal and cancer cells under different micro-environmental
conditions.


  We describe a new generation of algorithms capable of mapping the structure
and conformations of macromolecules and their complexes from large ensembles of
heterogeneous snapshots, and demonstrate the feasibility of determining both
discrete and continuous macromolecular conformational spectra. These algorithms
naturally incorporate conformational heterogeneity without resort to sorting
and classification, or prior knowledge of the type of heterogeneity present.
They are applicable to single-particle diffraction and image datasets produced
by X-ray lasers and cryo-electron microscopy, respectively, and particularly
suitable for systems not easily amenable to purification or crystallization.


  Background: A common modelling approach in public health and epidemiology
divides the population under study into compartments containing persons that
share the same status. Here we consider a three-state model with the
compartments: A, B and Dead. States A and B may be the states of any
dichotomous variable, for example, Healthy and Ill, respectively. The
transitions between the states are described by change rates (or synonymously:
densities), which depend on calendar time and on age. So far, a rigorous
mathematical calculation of the prevalence of property B has been difficult,
which has limited the use of the model in epidemiology and public health.
  Methods: We develop an equation that simplifies the use of the three-state
model. To demonstrate the broad applicability and the validity of the equation,
it is applied to simulated data and real world data from different
health-related topics.
  Results: The three-state model is governed by a partial differential equation
(PDE) that links the prevalence with the change rates between the states. The
validity of the PDE has been shown in two simulation studies, one about a
hypothetical chronic disease and one about dementia. In two further
applications, the equation may provide insights into smoking behaviour of males
in Germany and the knowledge about the ovulatory cycle in Egyptian women.
  Conclusions: We have found a simple equation that links the prevalence of a
dichotomous variable with the transmission rates in the three-state model. The
equation has a broad applicability in epidemiology and public health. Examples
are the estimation of incidence rates from cross-sectional surveys, the
prediction of the future prevalence of chronic diseases, and planning of
interventions against risky behaviour (e.g., smoking).


  Transcriptional regulation of gene expression is one of the main processes
that affect cell diversification from a single set of genes. Regulatory
proteins often interact with DNA regions located distally from the
transcription start sites (TSS) of the genes. We developed a computational
method that combines open chromatin and gene expression information for a large
number of cell types to identify these distal regulatory elements. Our method
builds correlation graphs for publicly available DNase-seq and exon array
datasets with matching samples and uses graph-based methods to filter findings
supported by multiple datasets and remove indirect interactions. The resulting
set of interactions was validated with both anecdotal information of known
long-range interactions and unbiased experimental data deduced from Hi-C and
CAGE experiments. Our results provide a novel set of high-confidence candidate
open chromatin regions involved in gene regulation, often located several Mb
away from the TSS of their target gene.


  Several recent works have shown that protein structure can predict
site-specific evolutionary sequence variation. In particular, sites that are
buried and/or have many contacts with other sites in a structure have been
shown to evolve more slowly, on average, than surface sites with few contacts.
Here, we present a comprehensive study of the extent to which numerous
structural properties can predict sequence variation. The quantities we
considered include buriedness (as measured by relative solvent accessibility),
packing density (as measured by contact number), structural flexibility (as
measured by B factors, root-mean-square fluctuations, and variation in dihedral
angles), and variability in designed structures. We obtained structural
flexibility measures both from molecular dynamics simulations performed on 9
non-homologous viral protein structures and from variation in homologous
variants of those proteins, where available. We obtained measures of
variability in designed structures from flexible-backbone design in the Rosetta
software. We found that most of the structural properties correlate with site
variation in the majority of structures, though the correlations are generally
weak (correlation coefficients of 0.1 to 0.4). Moreover, we found that
buriedness and packing density were better predictors of evolutionary variation
than was structural flexibility. Finally, variability in designed structures
was a weaker predictor of evolutionary variability than was buriedness or
packing density, but it was comparable in its predictive power to the best
structural flexibility measures. We conclude that simple measures of buriedness
and packing density are better predictors of evolutionary variation than are
more complicated predictors obtained from dynamic simulations, ensembles of
homologous structures, or computational protein design.


  Eukaryotic cell development has been optimized by natural selection to obey
maximal intracellular flux of messenger proteins. This, in turn, implies
maximum Fisher information on angular position about a target nuclear pore
complex (NPR). The cell is simply modeled as spherical, with cell membrane (CM)
diameter 10 micron and concentric nuclear membrane (NM) diameter 6 micron. The
NM contains about 3000 nuclear pore complexes (NPCs). Development requires
messenger ligands to travel from the CM-NPC-DNA target binding sites. Ligands
acquire negative charge by phosphorylation, passing through the cytoplasm over
Newtonian trajectories toward positively charged NPCs (utilizing positive
nuclear localization sequences). The CM-NPC channel obeys maximized mean
protein flux F and Fisher information I at the NPC, with first-order delta I =
0 and approximate 2nd-order delta I = 0 stability to environmental
perturbations. Many of its predictions are confirmed, including the dominance
of protein pathways of from 1-4 proteins, a 4nm size for the EGFR protein and
the approximate flux value F =10^16 proteins/m2-s. After entering the nucleus,
each protein ultimately delivers its ligand information to a DNA target site
with maximum probability, i.e. maximum Kullback-Liebler entropy HKL. In a
smoothness limit HKL approaches IDNA/2, so that the total CM-NPC-DNA channel
obeys maximum Fisher I. Thus maximum information approaches non-equilibrium,
one condition for life.


  Bow-tie or hourglass structure is a common architectural feature found in
biological and technological networks. A bow-tie in a multi-layered structure
occurs when intermediate layers have much fewer components than the input and
output layers. Examples include metabolism where a handful of building blocks
mediate between multiple input nutrients and multiple output biomass
components, and signaling networks where information from numerous receptor
types passes through a small set of signaling pathways to regulate multiple
output genes. Little is known, however, about how bow-tie architectures evolve.
Here, we address the evolution of bow-tie architectures using simulations of
multi-layered systems evolving to fulfill a given input-output goal. We find
that bow-ties spontaneously evolve when two conditions are met: (i) the
evolutionary goal is rank deficient, where the rank corresponds to the minimal
number of input features on which the outputs depend, and (ii) The effects of
mutations on interaction intensities between components are described by
product rule - namely the mutated element is multiplied by a random number.
Product-rule mutations are more biologically realistic than the commonly used
sum-rule mutations that add a random number to the mutated element. These
conditions robustly lead to bow-tie structures. The minimal width of the
intermediate network layers (the waist or knot of the bow-tie) equals the rank
of the evolutionary goal. These findings can help explain the presence of
bow-ties in diverse biological systems, and can also be relevant for machine
learning applications that employ multi-layered networks.


  Transposable elements can be categorised into DNA and RNA elements based on
their mechanism of transposition. Tyrosine recombinase elements (YREs) are
relatively rare and poorly understood, despite sharing characteristics with
both DNA and RNA elements. Previously, the Nematoda have been reported to have
a substantially different diversity of YREs compared to other animal phyla: the
Dirs1-like YRE retrotransposon was encountered in most animal phyla but not in
Nematoda, and a unique Pat1-like YRE retrotransposon has only been recorded
from Nematoda. We explored the diversity of YREs in Nematoda by sampling
broadly across the phylum and including 34 genomes representing the three
classes within Nematoda. We developed a method to isolate and classify YREs
based on both feature organization and phylogenetic relationships in an open
and reproducible workflow. We also ensured that our phylogenetic approach to
YRE classification identified truncated and degenerate elements, informatively
increasing the number of elements sampled. We identified Dirs1-like elements
(thought to be absent from Nematoda) in the nematode classes Enoplia and
Dorylaimia indicating that nematode model species do not adequately represent
the diversity of transposable elements in the phylum. Nematode Pat1-like
elements were found to be a derived form of another pat1-like element that is
present more widely in animals. Several sequence features used widely for the
classification of YREs were found to be homoplasious, highlighting the need for
a phylogenetically-based classification scheme. Nematode model species do not
represent the diversity of transposable elements in the phylum.


  Modeling of DNA-protein interactions is a complex process involving many
important time and length scales. This can be facilitated through the use of
coarse-grained models which reduce the number of degrees of freedom and allow
efficient exploration of binding configurations. It is known that the local
structure of DNA can significantly affect its protein-binding properties (i.e.
intrinsic curvature in DNA-histone complexes). In a step towards comprehensive
DNA-protein modeling, we expand the 3SPN.2 coarse-grained model to include
intrinsic shape, and validate the refined model against experimental data
including melting temperature, local flexibility, persistence length, and minor
groove width profile.


  Why are large, complex ecosystems stable? Both theory and simulations of
current models predict the onset of instability with growing size and
complexity, so for decades it has been conjectured that ecosystems must have
some unidentified structural property exempting them from this outcome. We show
that 'trophic coherence' -- a hitherto ignored feature of food webs which
current structural models fail to reproduce -- is a better statistical
predictor of linear stability than size or complexity. Furthermore, we prove
that a maximally coherent network with constant interaction strengths will
always be linearly stable. We also propose a simple model which, by correctly
capturing the trophic coherence of food webs, accurately reproduces their
stability and other basic structural features. Most remarkably, our model shows
that stability can increase with size and complexity. This suggests a key to
May's Paradox, and a range of opportunities and concerns for biodiversity
conservation.


  Introgressions from Neanderthals and Denisovans were detected in modern
humans. Introgressions from other archaic hominins were also implicated,
however, identification of which poses a great technical challenge. Here, we
introduced an approach in identifying introgressions from all possible archaic
hominins in Eurasian genomes, without referring to archaic hominin sequences.
We focused on mutations emerged in archaic hominins after their divergence from
modern humans (denoted as archaic-specific mutations), and identified
introgressive segments which showed significant enrichment of archaic-specific
mutations over the rest of the genome. Furthermore, boundaries of
introgressions were identified using a dynamic programming approach to
partition whole genome into segments which contained different levels of
archaic-specific mutations. We found that detected introgressions shared more
archaic-specific mutations with Altai Neanderthal than they shared with
Denisovan, and 60.3% of archaic hominin introgressions were from Neanderthals.
Furthermore, we detected more introgressions from two unknown archaic hominins
whom diverged with modern humans approximately 859 and 3,464 thousand years
ago. The latter unknown archaic hominin contributed to the genomes of the
common ancestors of modern humans and Neanderthals. In total, archaic hominin
introgressions comprised 2.4% of Eurasian genomes. Above results suggested a
complex admixture history among hominins. The proposed approach could also
facilitate admixture research across species.


  As part of a program to develop compounds with potential to treat cocaine
abuse, eleven (+/-)-threo-methylphenidate (TMP; Ritalin) derivatives were
synthesized and tested in rat striatal tissue preparations for inhibitory
potency against [3H]WIN 35,428 binding (WIN) to the dopamine (DA) transporter,
[3H]citalopram binding (CIT) to the serotonin transporter, and [3H]DA uptake.
The ester function was replaced by other functional groups in all of the
compounds; some also contained substituents on the phenyl ring and/or the
piperidine nitrogen. Potencies against WIN, measured as IC50, ranged from 27 nM
to 7,000 nM, compared to an IC50 of 83 nM for TMP itself. Potency against
[3H]DA uptake was approximately two-fold less than that against WIN, but
generally exhibited the same rank order. With one exception, the compounds were
significantly less potent against CIT than WIN. The one exception, which has a
rigid planar conformation at the altered ester position, is unique in that it
also is much less potent against [3H]DA uptake relative to WIN, compared to the
other derivatives. The three compounds with dichloro groups on the phenyl ring
did not exhibit positive cooperativity, as has been observed with several
previously synthesized halogenated TMP derivatives. Taken together, these
compounds should help to further our understanding of the stimulant binding
sites on both the dopamine and serotonin transporters.


  We introduce a framework within evolutionary game theory for studying the
distinction between objective and subjective rationality and apply it to the
evolution of cooperation on 3-regular random graphs. In our simulations, agents
evolve misrepresentations of objective reality that help them cooperate and
maintain higher social welfare in the Prisoner's dilemma. These agents act
rationally on their subjective representations of the world, but irrationally
from the perspective of an external observer. We model misrepresentations as
subjective perceptions of payoffs and quasi-magical thinking as an inferential
bias, finding that the former is more conducive to cooperation. This highlights
the importance of internal representations, not just observed behavior, in
evolutionary thought. Our results provide support for the interface theory of
perception and suggest that the individual's interface can serve not only the
individual's aims, but also society as a whole, offering insight into social
phenomena such as religion.


  Batch effects, artificial sources of variation due to experimental design,
are a widespread phenomenon in high throughput data. Therefore, mechanisms for
detection of batch effects are needed requiring comparison of multiple samples.
We apply hierarchical clustering (HC) on DNA k-mer counts of multiple RNA-seq
derived Fastq files. Ideally, HC generated trees reflect experimental treatment
groups and thus may indicate experimental effects, but clustering of
preparation groups indicates the presence of batch effects. In order to provide
a simple applicable tool we implemented sequential analysis of Fastq reads with
low memory usage in an R package (seqTools) available on Bioconductor. DNA
k-mer counts were analysed on 61 Fastq files containing RNA-seq data from two
cell types (dermal fibroblasts and Jurkat cells) sequenced on 8 different
Illumina Flowcells. Results: Pairwise comparison of all Flowcells with
hierarchical clustering revealed strong Flowcell based tree separation in 6 (21
%) and detectable Flowcell based clustering in 17 (60.7 %) of 28 Flowcell
comparisons. In our samples, batch effects were also present in reads mapped to
the human genome. Filtering reads for high quality (Phred >30) did not remove
the batch effects. Conclusions: Hierarchical clustering of DNA k-mer counts
provides a quality criterion and an unspecific diagnostic tool for RNA-seq
experiments.


  Many structures inside the cell such as nucleosomes and protein-mediated DNA
loops contain sharply bent double-stranded (ds) DNA. Therefore, the energetics
of strong dsDNA bending constitutes an essential part of cellular
thermodynamics. Although the thermomechanical behavior of long dsDNA is well
described by the worm-like chain (WLC) model, the length limit of such elastic
behavior remains controversial. To investigate the energetics of strong dsDNA
bending, we measured the opening rate of small dsDNA loops with contour lengths
of 40-200 bp using Fluorescence Resonance Energy Transfer (FRET). From the
measured relationship of loop stability to loop size, we observed a transition
between two separate bending regimes at a critical loop size below 100 bp.
Above this loop size, the loop lifetime decreased with decreasing loop size in
a manner consistent with an elastic bending stress. Below the critical loop
size, however, the loop lifetime became less sensitive to loop size, indicative
of softening of the double helix. The critical loop size was measured to be ~60
bp with sodium only and ~100 bp with 5 mM magnesium, which suggests that
magnesium facilitates the softening transition. We show that our results are in
quantitative agreement with the kinkable worm-like chain model. Furthermore,
the model parameters constrained by our data can reproduce previously measured
J factors between 50 and 200 bp. Our work provides powerful means to study
dsDNA bending in the strong bending regime.


  Tandem duplications are an essential source of genetic novelty, and their
variation in natural populations is expected to influence adaptive walks. Here,
we describe evolutionary impacts of recently-derived, segregating tandem
duplications in Drosophila yakuba and Drosophila simulans. We observe an excess
of duplicated genes involved in defense against pathogens, insecticide
resistance, chorion development, cuticular peptides, and lipases or
endopeptidases associated with the accessory glands, suggesting that
duplications function in Red Queen dynamics and rapid evolution. We document
evidence of widespread selection on the D. simulans X, suggesting adaptation
through duplication is common on the X. Despite the evidence for positive
selection, duplicates display an excess of low frequency variants consistent
with largely detrimental impacts, limiting the variation that can effectively
facilitate adaptation. Although we observe hundreds of gene duplications, we
show that segregating variation is insufficient to provide duplicate copies of
the entire genome, and the number of duplications in the population spans
13.4\% of major chromosome arms in D. yakuba and 9.7\% in D. simulans. Whole
gene duplication rates are low at $1.17\times10^{-9}$ per gene per generation
in D. yakuba and $6.03\times10^{-10}$ per gene per generation in D. simulans,
suggesting long wait times for new mutations on the order of thousands of years
for the establishment of sweeps. Hence, in cases where adaption depends on
individual tandem duplications, evolution will be severely limited by mutation.
We observe low levels of parallel recruitment of the same duplicated gene in
different species, suggesting that the span of standing variation will define
evolutionary outcomes in spite of convergence across gene ontologies consistent
with rapidly evolving phenotypes.} }


  The dynamics of stochastic reaction networks within cells are inevitably
modulated by factors considered extrinsic to the network such as for instance
the fluctuations in ribsome copy numbers for a gene regulatory network. While
several recent studies demonstrate the importance of accounting for such
extrinsic components, the resulting models are typically hard to analyze. In
this work we develop a general mathematical framework that allows to uncouple
the network from its dynamic environment by incorporating only the
environment's effect onto the network into a new model. More technically, we
show how such fluctuating extrinsic components (e.g., chemical species) can be
marginalized in order to obtain this decoupled model. We derive its
corresponding process- and master equations and show how stochastic simulations
can be performed. Using several case studies, we demonstrate the significance
of the approach. For instance, we exemplarily formulate and solve a marginal
master equation describing the protein translation and degradation in a
fluctuating environment.


  Maximum likelihood estimation (MLE) methods are widely used for evolutionary
tree. As evolutionary tree is not a smooth parameter, the consistency of its
MLE has been a topic of debate. It has been noted without proof that the
classical proof of consistency by Wald holds for the MLE of evolutionary tree.
Other proofs of consistency under various models were also proposed. Here we
will discuss some shortcomings in some of these proofs and comment on the
applicability of Wald's proof.


  RNA-sequencing has revolutionized biomedical research and, in particular, our
ability to study gene alternative splicing. The problem has important
implications for human health, as alternative splicing may be involved in
malfunctions at the cellular level and multiple diseases. However, the
high-dimensional nature of the data and the existence of experimental biases
pose serious data analysis challenges. We find that the standard data summaries
used to study alternative splicing are severely limited, as they ignore a
substantial amount of valuable information. Current data analysis methods are
based on such summaries and are hence suboptimal. Further, they have limited
flexibility in accounting for technical biases. We propose novel data summaries
and a Bayesian modeling framework that overcome these limitations and determine
biases in a nonparametric, highly flexible manner. These summaries adapt
naturally to the rapid improvements in sequencing technology. We provide
efficient point estimates and uncertainty assessments. The approach allows to
study alternative splicing patterns for individual samples and can also be the
basis for downstream analyses. We found a severalfold improvement in estimation
mean square error compared popular approaches in simulations, and substantially
higher consistency between replicates in experimental data. Our findings
indicate the need for adjusting the routine summarization and analysis of
alternative splicing RNA-seq studies. We provide a software implementation in
the R package casper.