  Axonemal dyneins are the molecular motors responsible for the beating of
cilia and flagella. These motors generate sliding forces between adjacent
microtubule doublets within the axoneme, the motile cytoskeletal structure
inside the flagellum. To create regular, oscillatory beating patterns, the
activities of the axonemal dyneins must be coordinated both spatially and
temporally. It is thought that coordination is mediated by stresses or strains
that build up within the moving axoneme, but it is not known which components
of stress or strain are involved, nor how they feed back on the dyneins. To
answer this question, we used isolated, reactivate axonemes of the unicellular
alga Chlamydomonas as a model system. We derived a theory for beat regulation
in a two-dimensional model of the axoneme. We then tested the theory by
measuring the beat waveforms of wild type axonemes, which have asymmetric
beats, and mutant axonemes, in which the beat is nearly symmetric, using
high-precision spatial and temporal imaging. We found that regulation by
sliding forces fails to account for the measured beat, due to the short lengths
of Chlamydomonas cilia. We found that regulation by normal forces (which tend
to separate adjacent doublets) cannot satisfactorily account for the symmetric
waveforms of the mbo2 mutants. This is due to the model's failure to produce
reciprocal inhibition across the axes of the symmetrically beating axonemes.
Finally, we show that regulation by curvature accords with the measurements.
Unexpectedly, we found that the phase of the curvature feedback indicates that
the dyneins are regulated by the dynamic (i.e. time-varying) component of
axonemal curvature, but not by the static one. We conclude that a high-pass
filtered curvature signal is a good candidate for the signal that feeds back to
coordinate motor activity in the axoneme.


  Although networks provide a powerful approach to study a large variety of
ecological systems, their formulation does not typically account for multiple
interaction types, interactions that vary in space and time, and interconnected
systems such as networks of networks. The emergent field of `multilayer
networks' provides a natural framework for extending analyses of ecological
systems to include such multiple layers of complexity, as it specifically
allows one to differentiate and model `intralayer' and `interlayer'
connectivity. The framework provides a set of concepts and tools that can be
adapted and applied to ecology, facilitating research on high-dimensional,
heterogeneous systems in nature. Here, we formally define ecological multilayer
networks based on a review of previous and related approaches, illustrate their
application and potential with analyses of existing data, and discuss
limitations, challenges, and future applications. The integration of multilayer
network theory into ecology offers largely untapped potential to further
address ecological complexity, to ultimately provide new theoretical and
empirical insights into the architecture and dynamics of ecological systems.


  Bovine spongiform encephalopathy, otherwise known as mad cow disease, can
spread when an individual cow consumes feed containing the infected tissues of
another individual, forming a one-species feedback loop. Such feedback is the
primary means of transmission for BSE during epidemic conditions. Following
outbreaks in the European Union and elsewhere, many governments enacted
legislation designed to limit the spread of such diseases via elimination or
reduction of one-species feedback loops in agricultural systems. However,
two-species feedback loops---those in which infectious material from
one-species is consumed by a secondary species whose tissue is then consumed by
the first species---were not universally prohibited and have not been studied
before. Here we present a basic ecological disease model which examines the
role feedback loops may play in the spread of BSE and related diseases. Our
model shows that there are critical thresholds between the infection's
expansion and decrease related to the lifespan of the hosts, the growth rate of
the prions, and the amount of prions circulating between hosts. The ecological
disease dynamics can be intrinsically oscillatory, having outbreaks as well as
refractory periods which can make it appear that the disease is under control
while it is still increasing. We show that non-susceptible species that have
been intentionally inserted into a feedback loop to stop the spread of disease
do not, strictly by themselves, guarantee its control, though they may give
that appearance by increasing the refractory period of an epidemic's
oscillations. We suggest ways in which age-related dynamics and cross-species
coupling should be considered in continuing evaluations aimed at maintaining a
safe food supply.


  We present a nonparametric framework to model a short sequence of probability
distributions that vary both due to underlying effects of sequential
progression and confounding noise. To distinguish between these two types of
variation and estimate the sequential-progression effects, our approach
leverages an assumption that these effects follow a persistent trend. This work
is motivated by the recent rise of single-cell RNA-sequencing experiments over
a brief time course, which aim to identify genes relevant to the progression of
a particular biological process across diverse cell populations. While
classical statistical tools focus on scalar-response regression or
order-agnostic differences between distributions, it is desirable in this
setting to consider both the full distributions as well as the structure
imposed by their ordering. We introduce a new regression model for ordinal
covariates where responses are univariate distributions and the underlying
relationship reflects consistent changes in the distributions over increasing
levels of the covariate. This concept is formalized as a "trend" in
distributions, which we define as an evolution that is linear under the
Wasserstein metric. Implemented via a fast alternating projections algorithm,
our method exhibits numerous strengths in simulations and analyses of
single-cell gene expression data.


  Originally developed to elucidate the mechanisms of natural selection in
bacteria, the Luria-Delbr\"uck model assumed that cells are intrinsically
capable of dividing an unlimited number of times. This assumption however, is
not true for human somatic cells which undergo replicative senescence.
Replicative senescence is thought to act as a mechanism to protect against
cancer and the escape from it is a rate-limiting step in cancer progression.
Here we introduce a Luria-Delbr\"uck model that explicitly takes into account
cellular replication limits in the wild type cell population and models the
emergence of mutants that escape replicative senescence. We present results on
the mean, variance, distribution, and asymptotic behavior of the mutant
population in terms of three classical formulations of the problem. More
broadly the paper introduces the concept of incorporating replicative limits as
part of the Luria-Delbr\"uck mutational framework. Guidelines to extend the
theory to include other types of mutations and possible applications to the
modeling of telomere crisis and fluctuation analysis are also discussed.


  The organizations of insect societies, such as division of labor, task
allocation, collective regulation, mass action responses, have been considered
as main reasons for the ecological success. In this article, we propose and
study a general modeling framework that includes the following three features:
(a) the average internal response threshold for each task (the internal
factor); (b) social network communications that could lead to task switching
(the environmental factor); and (c) dynamical changes of task demands (the
external factor). Since workers in many social insect species exhibit \emph{age
polyethism}, we also extend our model to incorporate \emph{age polyethism} in
which worker task preferences change with age. We apply our general modeling
framework to the cases of two task groups: the inside colony task versus the
outside colony task. Our analytical study of the models provides important
insights and predictions on the effects of colony size, social communication,
and age related task preferences on task allocation and division of labor in
the adaptive dynamical environment. Our study implies that the smaller size
colony invests its resource for the colony growth and allocates more workers in
the risky tasks such as foraging while the larger colony shifts more workers to
perform the safer tasks inside the colony. Social interactions among different
task groups play an important role in shaping task allocation depending on the
relative cost and demands of the tasks.


  We consider the correction of errors from nucleotide sequences produced by
next-generation targeted amplicon sequencing. The next-generation sequencing
(NGS) platforms can provide a great deal of sequencing data thanks to their
high throughput, but the associated error rates often tend to be high.
Denoising in high-throughput sequencing has thus become a crucial process for
boosting the reliability of downstream analyses. Our methodology, named
DUDE-Seq, is derived from a general setting of reconstructing finite-valued
source data corrupted by a discrete memoryless channel and effectively corrects
substitution and homopolymer indel errors, the two major types of sequencing
errors in most high-throughput targeted amplicon sequencing platforms. Our
experimental studies with real and simulated datasets suggest that the proposed
DUDE-Seq not only outperforms existing alternatives in terms of
error-correction capability and time efficiency, but also boosts the
reliability of downstream analyses. Further, the flexibility of DUDE-Seq
enables its robust application to different sequencing platforms and analysis
pipelines by simple updates of the noise model. DUDE-Seq is available at
http://data.snu.ac.kr/pub/dude-seq.


  Motivated by the need for fast and accurate classification of unlabeled
nucleotide sequences on a large scale, we developed NASCUP, a new
classification method that captures statistical structures of nucleotide
sequences by compact context-tree models and universal probability from
information theory. NASCUP achieved BLAST-like classification accuracy
consistently for several large-scale databases in orders-of-magnitude reduced
runtime, and was applied to other bioinformatics tasks such as outlier
detection and synthetic sequence generation.


  Multiply inverted balancer chromosomes that suppress exchange with their
homologs are an essential part of the genetic toolkit in Drosophila
melanogaster. Despite their widespread use, the organization of balancer
chromosomes has not been characterized at the molecular level, and the degree
of sequence variation among copies of any given balancer chromosome is unknown.
To map inversion breakpoints and study potential sequence diversity in the
descendants of a structurally identical balancer chromosome, we sequenced a
panel of laboratory stocks containing the most widely used X-chromosome
balancer, First Multiple 7 (FM7). We mapped the locations of FM7 breakpoints to
precise euchromatic coordinates and identified the flanking sequence of
breakpoints in heterochromatic regions. Analysis of SNP variation revealed
megabase-scale blocks of sequence divergence among currently used FM7 stocks.
We present evidence that this divergence arose by rare double crossover events
that replaced a female-sterile allele of the singed gene (sn[X2]) on FM7c with
wild type sequence from balanced chromosomes, and propose that many FM7c
chromosomes in the Bloomington Drosophila Stock Center have lost sn[X2] by this
mechanism. Finally, we characterize the original allele of the Bar gene (B[1])
that is carried on FM7 and validate the hypothesis that the origin and
subsequent reversion of the B1 duplication is mediated by unequal exchange. Our
results reject a simple non-recombining, clonal mode for the laboratory
evolution of balancer chromosomes and have implications for how balancer
chromosomes should be used in the design and interpretation of genetic
experiments in Drosophila.


  We investigate a partial differential equation model of a cancer cell
population, which is structured with respect to age and telomere length of
cells. We assume a continuous telomere length structure, which is applicable to
the clonal evolution model of cancer cell growth. This model has a non-standard
non-local boundary condition. We establish global existence of solutions and
study their qualitative behaviour. We study the effect of telomere restoration
on cancer cell dynamics. Our results indicate that without telomere
restoration, the cell population extinguishes. With telomere restoration,
exponential growth occurs in the linear model. We further characterise the
specific growth behaviour of the cell population for special cases. We also
study the effects of crowding induced mortality on the qualitative behaviour,
and the existence and stability of steady states of a nonlinear model
incorporating crowding effect. We present examples and extensive numerical
simulations, which illustrate the rich dynamic behaviour of the linear and
nonlinear models.


  We employ Monte Carlo simulations to study a stochastic Lotka-Volterra model
on a two-dimensional square lattice with periodic boundary conditions. If the
(local) prey carrying capacity is finite, there exists an extinction threshold
for the predator population that separates a stable active two-species
coexistence phase from an inactive state wherein only prey survive. Holding all
other rates fixed, we investigate the non-equilibrium relaxation of the
predator density in the vicinity of the critical predation rate. As expected,
we observe critical slowing-down, i.e., a power law dependence of the
relaxation time on the predation rate, and algebraic decay of the predator
density at the extinction critical point. The numerically determined critical
exponents are in accord with the established values of the directed percolation
universality class. Following a sudden predation rate change to its critical
value, one finds critical aging for the predator density autocorrelation
function that is also governed by universal scaling exponents. This aging
scaling signature of the active-to-absorbing state phase transition emerges at
significantly earlier times than the stationary critical power laws, and could
thus serve as an advanced indicator of the (predator) population's proximity to
its extinction threshold.


  Convolutional neural networks (CNN) have achieved state of the art
performance on both classification and segmentation tasks. Applying CNNs to
microscopy images is challenging due to the lack of datasets labeled at the
single cell level. We extend the application of CNNs to microscopy image
classification and segmentation using multiple instance learning (MIL). We
present the adaptive Noisy-AND MIL pooling function, a new MIL operator that is
robust to outliers. Combining CNNs with MIL enables training CNNs using full
resolution microscopy images with global labels. We base our approach on the
similarity between the aggregation function used in MIL and pooling layers used
in CNNs. We show that training MIL CNNs end-to-end outperforms several previous
methods on both mammalian and yeast microscopy images without requiring any
segmentation steps.


  We propose a mathematical framework for natural selection in finite
populations. Traditionally, many of the selection-based processes used to
describe cultural and genetic evolution (such as imitation and birth-death
models) have been studied on a case-by-case basis. Over time, these models have
grown in sophistication to include population structure, differing phenotypes,
and various forms of interaction asymmetry, among other features. Furthermore,
many processes inspired by natural selection, such as evolutionary algorithms
in computer science, possess characteristics that should fall within the realm
of a "selection process," but so far there is no overarching theory
encompassing these evolutionary processes. The framework of $\textit{stochastic
selection processes}$ we present here provides such a theory and consists of
three main components: a $\textit{population state space}$, an
$\textit{aggregate payoff function}$, and an $\textit{update rule}$. A
population state space is a generalization of the notion of population
structure, and it can include non-spatial information such as strategy-mutation
rates and phenotypes. An aggregate payoff function allows one to generically
talk about the fitness of traits without explicitly specifying a method of
payoff accounting or even the nature of the interactions that determine
payoff/fitness. An update rule is a fitness-based function that updates a
population based on its current state, and it includes as special cases the
classical update mechanisms (Moran, Wright-Fisher, etc.) as well as more
complicated mechanisms involving chromosomal crossover, mutation, and even
complex cultural syntheses of strategies of neighboring individuals. Our
framework covers models with variable population size as well as with
arbitrary, measurable trait spaces.


  Co-evolutionary adaptive mechanisms are not only ubiquitous in nature, but
also beneficial for the functioning of a variety of systems. We here consider
an adaptive network of oscillators with a stochastic, fitness-based, rule of
connectivity, and show that it self-organizes from fragmented and incoherent
states to connected and synchronized ones. The synchronization and percolation
are associated to abrupt transitions, and they are concurrently (and
significantly) enhanced as compared to the non-adaptive case. Finally we
provide evidence that only partial adaptation is sufficient to determine these
enhancements. Our study, therefore, indicates that inclusion of simple adaptive
mechanisms can efficiently describe some emergent features of networked
systems' collective behaviors, and suggests also self-organized ways to control
synchronization and percolation in natural and social systems.


  We consider a model of "isolation with an initial period of migration" (IIM),
where an ancestral population instantaneously split into two descendant
populations which exchanged migrants symmetrically at a constant rate for a
period of time but which are now completely isolated from each other. A method
of Maximum Likelihood estimation of the parameters of the model is implemented,
for data consisting of the number of nucleotide differences between two DNA
sequences at each of a large number of independent loci, using the explicit
analytical expressions for the likelihood obtained in Wilkinson-Herbots (2012).
The method is demonstrated on a large set of DNA sequence data from two species
of Drosophila, as well as on simulated data. The method is extremely fast,
returning parameter estimates in less than 1 minute for a data set consisting
of the numbers of differences between pairs of sequences from 10,000s of loci,
or in a small fraction of a second if all loci are trimmed to the same
estimated mutation rate. It is also illustrated how the maximized likelihood
can be used to quickly distinguish between competing models describing
alternative evolutionary scenarios, either by comparing AIC scores or by means
of likelihood ratio tests. The present implementation is for a simple version
of the model, but various extensions are possible and are briefly discussed.


  Background: Identification of biological specimens is a major requirement for
a range of applications. Reference-free methods analyse unprocessed sequencing
data without relying on prior knowledge, but generally do not scale to
arbitrarily large genomes and arbitrarily large phylogenetic distances.
Results: We present Cnidaria, a practical tool for clustering genomic and
transcriptomic data with no limitation on genome size or phylogenetic
distances. We successfully simultaneously clustered 169 genomic and
transcriptomic datasets from 4 kingdoms, achieving 100% identification accuracy
at supra-species level and 78% accuracy for species level. Discussion: CNIDARIA
allows for fast, resource-efficient comparison and identification of both raw
and assembled genome and transcriptome data. This can help answer both
fundamental (e.g. in phylogeny, ecological diversity analysis) and practical
questions (e.g. sequencing quality control, primer design).


  Genome sizes have evolved to vary widely, from 250 bases in viroids to 670
billion bases in amoeba. This remarkable variation in genome size is the
outcome of complex interactions between various evolutionary factors such as
point mutation rate, population size, insertions and deletions, and genome
editing mechanisms that may be specific to certain taxonomic lineages. While
comparative genomics analyses have uncovered some of the relationships between
these diverse evolutionary factors, we still do not understand what drives
genome size evolution. Specifically, it is not clear how primordial mutational
processes of base substitutions, insertions, and deletions influence genome
size evolution in asexual organisms. Here, we use digital evolution to
investigate genome size evolution by tracking genome edits and their fitness
effects in real time. In agreement with empirical data, we find that mutation
rate is inversely correlated with genome size in asexual populations. We show
that at low point mutation rate, insertions are significantly more beneficial
than deletions, driving genome expansion and acquisition of phenotypic
complexity. Conversely, high mutational load experienced at high mutation rates
inhibits genome growth, forcing the genomes to compress genetic information.
Our analyses suggest that the inverse relationship between mutation rate and
genome size is a result of the tradeoff between evolving phenotypic innovation
and limiting the mutational load.


  Coastal marine waters are ranked among the most important aquatic ecosystems
on earth in terms of ecological and economic significance. Since the Industrial
Revolution, human activities have drastically changed coastal marine
ecosystems. The development of rules and regulations to protect these
ecosystems against human activities requires availability of environmental
assessment standard. This necessitates the identification of the key parameters
that reflect condition of the coastal water ecosystem. Macrobenthic assemblages
are recognized to rapidly respond to changes in the quality of water or
habitat. Therefore, it would be useful to study the population of macrobenthos
and assess the influential factors on the growth of this species. This study is
categorized as multidisciplinary approach which contains two perspectives,
ecological and mathematical. In the ecological section, the effect of the water
quality parameters (e.g. pH, temperature, dissolved oxygen and salinity) and
the sediment characteristics on the macrobenthic abundance is studied. A total
of 432 samples were collected and analyzed from four touristic costal locations
(at various distances form the coast) of Penang National Park to investigate
the spatial change of macrobenthic assemblage. In a mathematical perspective,
this paper pursues a new algorithm based on the performance evaluation methods.
For this purpose, first, Data Envelopment Analysis (DEA) is employed to
evaluate a group of Decision Making Units (DMUs). Consequently, inputs of the
mentioned DMUs are considered as alternatives (or candidates), and using a
modified DEA model that is categorized as aggregating preference ranking
method, the influence of inputs in efficiency of DMUs is investigated.


  We propose a new geometric approach to describe the qualitative dynamics of
chemical reactions networks. By this method we identify metastable regimes,
defined as low dimensional regions of the phase space close to which the
dynamics is much slower compared to the rest of the phase space. Given the
network topology and the orders of magnitude of kinetic parameters, the number
of such metastable regimes is finite. The dynamics of the network can be
described as a sequence of jumps from one metastable regime to another. We show
that a geometrically computed connectivity graph restricts the set of possible
jumps. We also provide finite state machine (Markov chain) models for such
dynamic changes. Applied to signal transduction models, our approach unravels
dynamical and functional capacities of signaling pathways, as well as
parameters responsible for specificity of the pathway response. In particular,
for a model of TGF$\beta$ signalling, we find that the ratio of TGFBR1 to
TGFBR2 concentrations can be used to discriminate between metastable regimes.
Using expression data from the NCI60 panel of human tumor cell lines, we show
that aggressive and non-aggressive tumour cell lines function in different
metastable regimes and can be distinguished by measuring the relative
concentrations of receptors of the two types.


  The division of labor (DOL) and task allocation among groups of ants living
in a colony is thought to be highly efficient, and key to the robust survival
of a colony. A great deal of experimental and theoretical work has been done
toward gaining a clear understanding of the evolution of, and underlying
mechanisms of these phenomena. Much of this research has utilized mathematical
modeling. Here we continue this tradition by developing a mathematical model
for a particular aspect of task allocation, known as age-related repertoire
expansion, that has been observed in the minor workers of the ant species
\emph{Pheidole dentata}. In fact, we present a relatively broad mathematical
modeling framework based on the dynamics of the frequency with which members of
specific age groups carry out distinct tasks. We apply our modeling approach to
a specific task allocation scenario, and compare our theoretical results with
experimental data. It is observed that the model predicts perceived behavior,
and provides a possible explanation for the aforementioned experimental
results.


  Characterizing the link between small-scale chromatin structure and
large-scale chromosome folding during interphase is a prerequisite for
understanding transcription. Yet, this link remains poorly investigated. Here,
we introduce a simple biophysical model where interphase chromosomes are
described in terms of the folding of chromatin sequences composed of
alternating blocks of fibers with different thicknesses and flexibilities, and
we use it to study the influence of sequence disorder on chromosome behaviors
in space and time. By employing extensive computer simulations,we thus
demonstrate that chromosomes undergo noticeable conformational changes only on
length-scales smaller than $10^5$ basepairs and time-scales shorter than a few
seconds, and we suggest there might exist effective upper bounds to the
detection of chromosome reorganization in eukaryotes. We prove the relevance of
our framework by modeling recent experimental FISH data on murine chromosomes.


  We consider a multi-type Moran model (in continuous time) with selection and
type-dependent mutation. This paper is concerned with the evolution of
genealogical information forward in time. For this purpose we define and
analytically characterize a path-valued Markov process that contains in its
state at time $t$ the extended ancestral lines (adding genealogical distances)
of the population alive at time $t$. The main result is a representation for
the conditional distribution of the extended ancestral lines of a subpopulation
alive at a fixed time $T$ (present time) given the type information of the
subpopulation at time $T$ in terms of the distribution of the sample paths (up
to time $T$) of a special Markov process (different from the ancestral
selection graph) to which we refer as backward process. This representation
allows us both to prove that the extended ancestral lines converge in the limit
$T \to \infty$ if the type information converges in the limit $T \to \infty$
and to study the resulting limit of the extended ancestral lines by means of
the backward process. The limit theorem has two applications: First, we can
represent the stationary type distribution of the common ancestor type process
in terms of the equilibrium distribution of a functional of the backward
process, where in the two type case we recover the common ancestor process of
Fearnhead if we let the population size tend to infinity. Second, we obtain
that the conditioned genealogical distance of two individuals given the types
of the two individuals is distributed as a certain stopping time of a further
functional of the backward process which is a new approach towards a proof that
genealogical distances are stochastically smaller under selection.


  The function of membrane-embedded proteins such as ion channels depends
crucially on their conformation. We demonstrate how conformational changes in
asymmetric membrane proteins may be inferred from measurements of their
diffusion. Such proteins cause local deformations in the membrane, which induce
an extra hydrodynamic drag on the protein. Using membrane tension to control
the magnitude of the deformations and hence the drag, measurements of
diffusivity can be used to infer--- via an elastic model of the protein--- how
conformation is changed by tension. Motivated by recent experimental results
[Quemeneur et al., Proc. Natl. Acad. Sci. USA, 111 5083 (2014)] we focus on
KvAP, a voltage-gated potassium channel. The conformation of KvAP is found to
change considerably due to tension, with its `walls', where the protein meets
the membrane, undergoing significant angular strains. The torsional stiffness
is determined to be 26.8 kT at room temperature. This has implications for both
the structure and function of such proteins in the environment of a
tension-bearing membrane.


  We present a comparison between two different mathematical models used in the
description of the Ebola virus propagation currently occurring in West Africa.
In order to improve the prediction and the control of the propagation of the
virus, numerical simulations and optimal control of the two models for Ebola
are investigated. In particular, we study when the two models generate similar
results.


  Monitoring with implicated punishment is common in human societies to avert
freeriding on common goods. But is it effective in promoting public
cooperation? We show that the introduction of monitoring and implicated
punishment is indeed effective, as it transforms the public goods game to a
coordination game, thus rendering cooperation viable in infinite and finite
well-mixed populations. We also show that the addition of within-group
enforcement further promotes the evolution of public cooperation. However,
although the group size in this context has nonlinear effects on collective
action, an intermediate group size is least conductive to cooperative
behaviour. This contradicts recent field observations, where an intermediate
group size was declared optimal with the conjecture that group-size effects and
within-group enforcement are responsible. Our theoretical research thus
clarifies key aspects of monitoring with implicated punishment in human
societies, and additionally, it reveals fundamental group-size effects that
facilitate prosocial collective action.


  Data curation presents a challenge to all scientific disciplines to ensure
public availability and reproducibility of experimental data. Standards for
data preservation and exchange are central to addressing this challenge: the
Investigation-Study-Assay Tabular (ISA-Tab) project has developed a widely used
template for such standards in biological research. This paper describes the
application of ISA-Tab to protein titration data. Despite the importance of
titration experiments for understanding protein structure, stability, and
function and for testing computational approaches to protein electrostatics, no
such mechanism currently exists for sharing and preserving biomolecular
titration data. We have adapted the ISA-Tab template to provide a structured
means of supporting experimental structural chemistry data with a particular
emphasis on the calculation and measurement of pKa values. This activity has
been performed as part of the broader pKa Cooperative effort, leveraging data
that has been collected and curated by the Cooperative members. In this
article, we present the details of this specification and its application to a
broad range of pKa and electrostatics data obtained for multiple protein
systems. The resulting curated data is publicly available at
http://pkacoop.org.


  Bond graphs can be used to build thermodynamically-compliant hierarchical
models of biomolecular systems. As bond graphs have been widely used to model,
analyse and synthesise engineering systems, this paper suggests that they can
play the same role in the modelling, analysis and synthesis of biomolecular
systems. The particular structure of bond graphs arising from biomolecular
systems is established and used to elucidate the relation between
thermodynamically closed and open systems. Block diagram representations of the
dynamics implied by these bond graphs are used to reveal implicit feedback
structures and are linearised to allow the application of control-theoretical
methods.
  Two concepts of modularity are examined: computational modularity where
physical correctness is retained and behavioural modularity where module
behaviour (such as ultrasensitivity) is retained. As well as providing
computational modularity, bond graphs provide a natural formulation of
behavioural modularity and reveal the sources of retroactivity. A bond graph
approach to reducing retroactivity, and thus inter-module interaction, is shown
to require a power supply such as that provided by the ATP = ADP + Pi reaction.
  The MAPK cascade (Raf-MEK-ERK pathway) is used as an illustrative example.


  Rapidly assaying the diversity of a bacterial species present in a sample
obtained from a hospital patient or an evironmental source has become possible
after recent technological advances in DNA sequencing. For several applications
it is important to accurately identify the presence and estimate relative
abundances of the target organisms from short sequence reads obtained from a
sample. This task is particularly challenging when the set of interest includes
very closely related organisms, such as different strains of pathogenic
bacteria, which can vary considerably in terms of virulence, resistance and
spread. Using advanced Bayesian statistical modelling and computation
techniques we introduce a novel pipeline for bacterial identification that is
shown to outperform the currently leading pipeline for this purpose. Our
approach enables fast and accurate sequence-based identification of bacterial
strains while using only modest computational resources. Hence it provides a
useful tool for a wide spectrum of applications, including rapid clinical
diagnostics to distinguish among closely related strains causing nosocomial
infections. The software implementation is available at
https://github.com/PROBIC/BIB


  The hydrophobic effect stabilizes the native structure of proteins by
minimizing the unfavourable interactions between hydrophobic residues and water
through the formation of a hydrophobic core. Here we include the entropic and
enthalpic contributions of the hydrophobic effect explicitly in an implicit
solvent model. This allows us to capture two important effects: a length-scale
dependence and a temperature dependence for the solvation of a hydrophobic
particle. This consistent treatment of the hydrophobic effect explains cold
denaturation and heat capacity measurements of solvated proteins.


  Restriction site Associated DNA (RAD) tagging (also known as RAD-seq, etc.)
is an emerging method for analyzing an organism's genome without completely
sequencing it. This can be applied to a non-model organism without a reference
genome, though this creates the problem of how to begin data analysis on
unmapped and unannotated reads. Our program, Hot RAD, presents a
straightforward and easy-to-use method to take raw Illumina data that has been
RAD tagged and produce consensus contigs or sequence stacks using a distributed
framework, creating a basis on which to begin analyzing an organism's DNA. The
GUI (graphical user interface) element of our tool makes it easy for those not
familiar with the command line to take raw sequence files and produce usable
data in a timely manner.


  Recently, we proposed an state model (compartment model) to describe the
progression of a chronic disease with an pre-clinical (undiagnosed) state
before clinical diagnosis. It is an open question, if a sequence of
cross-sectional studies with mortality follow-up is sufficient to estimate the
true incidence rate of the disease, i.e. the incidence of the undiagnosed and
diagnosed disease. In this note, we construct a counterexample and show that
this cannot be achieved in general.


  In this projecting work we propose a mass spectrometric patch-clamp equipment
with the capillary performing both a local potential registration at the cell
membrane and the analyte suction simultaneously. This paper provides a current
literature analysis comparing the possibilities of the novel approach proposed
with the known methods, such as scanning patch-clamp, scanning ion conductance
microscopy, patch clamp based on scanning probe microscopy technology,
quantitative subcellular secondary ion mass spectrometry or "ion microscopy",
live single-cell mass spectrometry, in situ cell-by-cell imaging, single-cell
video-mass spectrometry, etc. We also consider the ways to improve the
informativeness of these methods and particularly emphasize the trend at the
increasing of the analysis complexity. We propose here the way to improve the
efficiency of the cell trapping to the capillary during MS-path-clamp, as well
as to provide laser surface ionization using laser trapping and tweezing of
cells with the laser beam transmitted through the capillary as a waveguide. It
is also possible to combine the above system with the microcolumn separation
system or capillary electrophoresis as an optional direction of further
development of the complex of analytical techniques emerging from the MS
variation of patch-clamp.


  Noise through its interaction with the nonlinearity of the living systems can
give rise to counter-intuitive phenomena. In this paper we shortly review noise
induced effects in different ecosystems, in which two populations compete for
the same resources. We also present new results on spatial patterns of two
populations, while modeling real distributions of anchovies and sardines. The
transient dynamics of these ecosystems are analyzed through generalized
Lotka-Volterra equations in the presence of multiplicative noise, which models
the interaction between the species and the environment. We find noise induced
phenomena such as quasi-deterministic oscillations, stochastic resonance, noise
delayed extinction, and noise induced pattern formation. In addition, our
theoretical results are validated with experimental findings. Specifically the
results, obtained by a coupled map lattice model, well reproduce the spatial
distributions of anchovies and sardines, observed in a marine ecosystem.
Moreover, the experimental dynamical behavior of two competing bacterial
populations in a meat product and the probability distribution at long times of
one of them are well reproduced by a stochastic microbial predictive model.


  We inquire how structure emerges during the process of protein folding. For
this we scrutinise col- lective many-atom motions during all-atom molecular
dynamics simulations. We introduce, develop and employ various topological
techniques, in combination with analytic tools that we deduce from the concept
of integrable models and structure of discrete nonlinear Schroedinger equation.
The example we consider is an alpha-helical subunit of the HIV envelope
glycoprotein gp41. The helical structure is stable when the subunit is part of
the biological oligomer. But in isolation the helix becomes unstable, and the
monomer starts deforming. We follow the process computationally. We interpret
the evolving structure both in terms of a backbone based Heisenberg spin chain
and in terms of a side chain based XY spin chain. We find that in both cases
the formation of protein super-secondary structure is akin the formation of a
topological Bloch domain wall along a spin chain. During the process we
identify three individual Bloch walls and we show that each of them can be
modelled with a very high precision in terms of a soliton solution to a
discrete nonlinear Schroedinger equation.


  We study the weak convergence of iterates of so-called centred kernel
quadratic stochastic operators. These iterations, in a population evolution
setting, describe the additive perturbation of the arithmetic mean of the
traits of an individual's parents and correspond to certain weighted sums of
independent random variables. We show that one can obtain weak convergence
results under rather mild assumptions on the kernel. Essentially it is
sufficient for the distribution of the perturbing random variable to have a
finite variance or have tails controlled by a power function. The advantage of
these conditions is that in many cases they are easily verifiable by an applied
user. Additionally, the representation by sums of random variables implies an
efficient simulation algorithm to obtain random variables approximately
following the law of the iterates of the quadratic stochastic operator, with
full control of the degree of approximation. Our results also indicate where
lies an intrinsic difficulty in the analysis of the behaviour of quadratic
stochastic operators.


  The subtree prune-and-regraft (SPR) distance metric is a fundamental way of
comparing evolutionary trees. It has wide-ranging applications, such as to
study lateral genetic transfer, viral recombination, and Markov chain Monte
Carlo phylogenetic inference. Although the rooted version of SPR distance can
be computed relatively efficiently between rooted trees using
fixed-parameter-tractable maximum agreement forest (MAF) algorithms, no MAF
formulation is known for the unrooted case. Correspondingly, previous
algorithms are unable to compute unrooted SPR distances larger than 7.
  In this paper, we substantially advance understanding of and computational
algorithms for the unrooted SPR distance. First we identify four properties of
optimal SPR paths, each of which suggests that no MAF formulation exists in the
unrooted case. Then we introduce the replug distance, a new lower bound on the
unrooted SPR distance that is amenable to MAF methods, and give an efficient
fixed-parameter algorithm for calculating it. Finally, we develop a
"progressive A*" search algorithm using multiple heuristics, including the TBR
and replug distances, to exactly compute the unrooted SPR distance. Our
algorithm is nearly two orders of magnitude faster than previous methods on
small trees, and allows computation of unrooted SPR distances as large as 14 on
trees with 50 leaves.


  Adhesion processes of biological membranes that enclose cells and cellular
organelles are essential for immune responses, tissue formation, and signaling.
These processes depend sensitively on the binding constant K2D of the
membrane-anchored receptor and ligand proteins that mediate adhesion, which is
difficult to measure in the 'two-dimensional' (2D) membrane environment of the
proteins. An important problem therefore is to relate K2D} to the binding
constant K3D} of soluble variants of the receptors and ligands that lack the
membrane anchors and are free to diffuse in three dimensions (3D). In this
article, we present a general theory for the binding constants K2D and K3D of
rather stiff proteins whose main degrees of freedom are translation and
rotation, along membranes and around anchor points 'in 2D', or unconstrained
'in 3D'. The theory generalizes previous results by describing how K2D depends
both on the average separation and thermal nanoscale roughness of the apposing
membranes, and on the length and anchoring flexibility of the receptors and
ligands. Our theoretical results for the ratio K2D/K3D of the binding constants
agree with detailed results from Monte Carlo simulations without any data
fitting, which indicates that the theory captures the essential features of the
'dimensionality reduction' due to membrane anchoring. In our Monte Carlo
simulations, we consider a novel coarse-grained model of biomembrane adhesion
in which the membranes are represented as discretized elastic surfaces, and the
receptors and ligands as anchored molecules that diffuse continuously along the
membranes and rotate at their anchor points.


  The adhesion of biological membranes is mediated by the binding of
membrane-anchored receptor and ligand proteins. Central questions are how the
binding kinetics of these proteins is affected by the membranes and by the
membrane anchoring of the proteins. In this article, we (i) present detailed
data for the binding of membrane-anchored proteins from coarse-grained
molecular dynamics simulations, and (ii) provide a theory that describes how
the binding kinetics depends on the average separation and thermal roughness of
the adhering membranes, and on the anchoring, lengths, and length variations of
the proteins. An important element of our theory is the tilt of bound
receptor-ligand complexes and transition-state complexes relative to the
membrane normals. This tilt results from an interplay of the anchoring energy
and rotational entropy of the complexes and facilitates the formation of
receptor-ligand bonds at membrane separations smaller than the preferred
separation for binding. In our simulations, we have considered both
lipid-anchored and transmembrane receptor and ligand proteins. We find that the
binding equilibrium constant and binding on-rate constant of lipid-anchored
proteins are considerably smaller than the binding constant and on-rate
constant of rigid transmembrane proteins with identical binding domains.


  Cilia and flagella are hairlike organelles that propel cells through fluid.
The active motion of the axoneme, the motile structure inside cilia and
flagella, is powered by molecular motors of the dynein family. These motors
generate forces and torques that slide and bend the microtubule doublets within
the axoneme. To create regular waveforms the activities of the dyneins must be
coordinated. It is thought that coordination is mediated by stresses due to
radial, transverse, or sliding deformations, that build up within the moving
axoneme. However, which particular component of the stress regulates the motors
to produce the observed flagellar waveforms remains an open question. To
address this question, we describe the axoneme as a three-dimensional bundle of
filaments and characterize its mechanics. We show that regulation of the motors
by radial and transverse stresses can lead to a coordinated flagellar motion
only in the presence of twist. By comparison, regulation by shear stress is
possible without twist. We calculate emergent beating patterns in twisted
axonemes resulting from regulation by transverse stresses. The waveforms are
similar to those observed in flagella of Chlamydomonas and sperm. Due to the
twist, the waveform has non-planar components, which result in swimming
trajectories such as twisted ribbons and helices that agree with observations.


  We present an extensive, systematic study of the Prisoner's Dilemma and
Snowdrift games on a square lattice under a synchronous, noiseless imitation
dynamics. We show that for both the occupancy of the network and the (random)
mobility of the agents there are intermediate values that may increase the
amount of cooperators in the system and new phases appear. We analytically
determine the transition lines between these phases and compare with the mean
field prediction and the observed behavior on a square lattice. We point out
which are the more relevant microscopic processes that entitle cooperators to
invade a population of defectors in the presence of mobility and discuss the
universality of these results.


  The protein-protein interaction (PPI) network is crucial for cellular
information processing and decision-making. With suitable inputs, PPI networks
drive the cells to diverse functional outcomes such as cell proliferation or
cell death. Here we characterize the structural controllability of a large
directed human PPI network comprised of 6,339 proteins and 34,813 interactions.
This allows us to classify proteins as "indispensable", "neutral" or
"dispensable", which correlates to increasing, no effect, or decreasing the
number of driver nodes in the network upon removal of that protein. We find
that 21% of the proteins in the PPI network are indispensable. Interestingly,
these indispensable proteins are the primary targets of disease-causing
mutations, human viruses, and drugs, suggesting that altering a network's
control property is critical for the transition between healthy and disease
states. Furthermore, analyzing copy number alterations data from 1,547 cancer
patients reveals that 56 genes that are frequently amplified or deleted in nine
different cancers are indispensable. Among the 56 genes, 46 of them have not
been previously associated with cancer. This suggests that controllability
analysis is very useful in identifying novel disease genes and potential drug
targets.


  The role of historical contingency in the origin of life is one of the great
unknowns in modern science. Only one example of life exists--one that proceeded
from a single self-replicating organism (or a set of replicating hyper-cycles)
to the vast complexity we see today in Earth's biosphere. We know that emergent
life has the potential to evolve great increases in complexity, but it is
unknown if evolvability is automatic given any self-replicating organism. At
the same time, it is difficult to test such questions in biochemical systems.
Laboratory studies with RNA replicators have had some success with exploring
the capacities of simple self-replicators, but these experiments are still
limited in both capabilities and scope. Here, we use the digital evolution
system Avida to explore the interaction between emergent replicators (rare
randomly-assembled self-replicators) and evolvability. We find that we can
classify fixed-length emergent replicators in Avida into two classes based on
functional analysis. One class is more evolvable in the sense of optimizing
their replication abilities. However, the other class is more evolvable in the
sense of acquiring evolutionary innovations. We tie this trade-off in
evolvability to the structure of the respective classes' replication machinery,
and speculate on the relevance of these results to biochemical replicators.


  Genetic Algorithms (GA) are a powerful set of tools for search and
optimization that mimic the process of natural selection, and have been used
successfully in a wide variety of problems, including evolving neural networks
to solve cognitive tasks. Despite their success, GAs sometimes fail to locate
the highest peaks of the fitness landscape, in particular if the landscape is
rugged and contains multiple peaks. Reaching distant and higher peaks is
difficult because valleys need to be crossed, in a process that (at least
temporarily) runs against the fitness maximization objective. Here we propose
and test a number of information-theoretic (as well as network-based) measures
that can be used in conjunction with a fitness maximization objective
(so-called ``neuro-correlates") to evolve neural controllers for two widely
different tasks: a behavioral task that requires information integration, and a
cognitive task that requires memory and logic. We find that judiciously chosen
neuro-correlates can significantly aid GAs to find the highest peaks.


  Validation of research methodology is critical in research design.
Correlation between experimental observables must be established before
undertaking extensive experiments or propose mechanisms. This article shows
that, observables in the popular calcium flux strength assay used in the
characterization of sweetener-sweet taste receptor (STR) interaction are
uncorrelated. In pursuit to find potential sweeteners and enhancers, calcium
flux generated via G-protein coupling for wildtype and mutant receptors
expressed on cell surface is measured to identify and localize sweetener
binding sites. Results are channeled for sweetener development with direct
impact on public health.
  We show that flux strength is independent of EC50 and sweet potency. Sweet
potency-EC50 relation is non-linear and anti-correlated. Single point mutants
affecting receptor efficiency, without significant shift in EC50 have been
published, indicating flux strength is independent of ligand binding. G-protein
coupling step is likely observed in the assay. Thus, years have been spent
generating uncorrelated data. Data from uncorrelated observables does not give
meaningful results. Still, majority of research in the field, uses change in
calcium flux strength to study the receptor. Methodology away from flux
strength monitor is required for sweetener development, reestablish binding
localization of sweeteners established by flux strength method. This article
serves to remind researchers to validate methodology before plunging into long
term projects. Ignoring validation test on methodology, have been a costly
mistake in the field. Concepts discussed here is applicable, whenever
observable in biological systems are many steps moved from the event of
interest.


  Despite the significant increase in computational power, molecular modeling
of protein structure using classical all-atom approaches remains inefficient,
at least for most of the protein targets in the focus of biomedical research.
Perhaps the most successful strategy to overcome the inefficiency problem is
multiscale modeling to merge all-atom and coarse-grained models. This chapter
describes a well-established CABS coarse-grained protein model. The CABS
(C-Alpha, C-Beta and Side chains) model assumes a 2-4 united-atom
representation of amino acids, knowledge-based force field (derived from the
statistical regularities seen in known protein sequences and structures) and
efficient Monte Carlo sampling schemes (MC dynamics, MC replica-exchange, and
combinations). A particular emphasis is given to the unique design of the CABS
force-field, which is largely defined using one-dimensional structural
properties of proteins, including protein secondary structure. This chapter
also presents CABS-based modeling methods, including multiscale tools for de
novo structure prediction, modeling of protein dynamics and prediction of
protein-peptide complexes. CABS-based tools are freely available at
http://biocomp.chem.uw.edu.pl/tools


  Bacterial chemotaxis systems are as diverse as the environments that bacteria
inhabit, but how much environmental variation can cells tolerate with a single
system? Diversification of a single chemotaxis system could serve as an
alternative, or even evolutionary stepping-stone, to switching between multiple
systems. We hypothesized that mutations in gene regulation could lead to
heritable control of chemotactic diversity. By simulating foraging and
colonization of E. coli using a single-cell chemotaxis model, we found that
different environments selected for different behaviors. The resulting
trade-offs show that populations facing diverse environments would ideally
diversify behaviors when time for navigation is limited. We show that
advantageous diversity can arise from changes in the distribution of protein
levels among individuals, which could occur through mutations in gene
regulation. We propose experiments to test our prediction that chemotactic
diversity in a clonal population could be a selectable trait that enables
adaptation to environmental variability.


  Metagenomics offers a way to analyze biotopes at the genomic level and to
reach functional and taxonomical conclusions. The bio-analyzes of large
metagenomic projects face critical limitations: complex metagenomes cannot be
assembled and the taxonomical or functional annotations are much smaller than
the real biological diversity. This motivated the development of de novo
metagenomic read comparison approaches to extract information contained in
metagenomic datasets.
  However, these new approaches do not scale up large metagenomic projects, or
generate an important number of large intermediate and result files. We
introduce Commet ("COmpare Multiple METagenomes"), a method that provides
similarity overview between all datasets of large metagenomic projects.
  Directly from non-assembled reads, all against all comparisons are performed
through an efficient indexing strategy. Then, results are stored as bit
vectors, a compressed representation of read files, that can be used to further
combine read subsets by common logical operations. Finally, Commet computes a
clusterization of metagenomic datasets, which is visualized by dendrogram and
heatmaps.
  Availability: http://github.com/pierrepeterlongo/commet


  Polymerases are protein enzymes that move along nucleic acid chains and
catalyze template-based polymerization reactions during gene transcription and
replication. The polymerases also substantially improve transcription or
replication fidelity through the non-equilibrium enzymatic cycles. We briefly
review computational efforts that have been made toward understanding
mechano-chemical coupling and fidelity control mechanisms of the polymerase
elongation. The polymerases are regarded as molecular information motors during
the elongation process. It requires a full spectrum of computational approaches
from multiple time and length scales to understand the full polymerase
functional cycle. We keep away from quantum mechanics based approaches to the
polymerase catalysis due to abundant former surveys, while address only
statistical physics modeling approach and all-atom molecular dynamics
simulation approach. We organize this review around our own modeling and
simulation practices on a single-subunit T7 RNA polymerase, and summarize
commensurate studies on structurally similar DNA polymerases. For multi-subunit
RNA polymerases that have been intensively studied in recent years, we leave
detailed discussions on the simulation achievements to other computational
chemical surveys, while only introduce very recently published representative
studies, including our own preliminary work on structure-based modeling on
yeast RNA polymerase II. In the end, we quickly go through kinetic modeling on
elongation pauses and backtracking activities. We emphasize the fluctuation and
control mechanisms of the polymerase actions, highlight the non-equilibrium
physical nature of the system, and try to bring some perspectives toward
understanding replication and transcription regulation from single molecular
details to a genome-wide scale.


  We perform a statistical study of the distances between successive
occurrencies of a given dinucleotide in the DNA sequence for a number of
organisms of different complexity. Our analysis highlights peculiar features of
the dinucleotide CG distribution in mammalian DNA, pointing towards a
connection with the role of such dinucleotide in DNA methylation. While the CG
distributions of mammals exhibit exponential tails with comparable parameters,
the picture for the other organisms studied (e.g., fish, insects, bacteria and
viruses) is more heterogeneous, possibly because in these organisms DNA
methylation has different functional roles. Our analysis suggests that the
distribution of the distances between dinucleotides CG provides useful insights
in characterizing and classifying organisms in terms of methylation
functionalities.


  The susceptible-infected-recovered (SIR) model has been used extensively to
model disease spread and other processes. Despite the widespread usage of this
ordinary differential equation (ODE) based model which represents the
mean-field approximation of the underlying stochastic SIR process on contact
networks, only few rigorous approaches exist and these use complex semigroup
and martingale techniques to prove that the expected fraction of the
susceptible and infected nodes of the stochastic SIR process on a complete
graph converges as the number of nodes increases to the solution of the
mean-field ODE model. Extending the elementary proof of convergence for the SIS
process introduced by Armbruster and Beck (2016) to the SIR process, we show
convergence using only a system of three ODEs, simple probabilistic
inequalities, and basic ODE theory. Our approach can also be generalized to
many other types of compartmental models (e.g.,
susceptible-infected-recovered-susceptible (SIRS)) which are linear ODEs with
the addition of quadratic terms for the number of new infections similar to the
SI term in the SIR model.


  It was proposed earlier that Pfcrmp (Plasmodium falciparum chloroquine
resistance marker protein) may be the chloroquine's target protein in nucleus.
In this communication, further evidence is presented to support the view that
Pfcrmp may play a key role in chloroquine antimalarial actions as well as
resistance development.


  MetaScope is a fast and accurate tool for analyzing (host-associated)
metagenome datasets. Sequence alignment of reads against the host genome (if
requested) and against microbial Genbank is performed using a new DNA aligner
called SASS. The output of SASS is processed so as to assign all microbial
reads to taxa and genes, using a new weighted version of the LCA algorithm.
MetaScope is the winner of the 2013 DTRA software challenge entitled "Identify
Organisms from a Stream of DNA Sequences".


  Time kinetics of interaction of natural graphite (GR) to colloidal graphene
(G) collagen (C) nanocomposites was studied at ambient conditions, and observed
that just one day at ambient conditions is enough to form colloidal graphene
directly from graphite using the protein collagen. Neither controlled
temperature and pressure ambiance nor sonication was needed for the same;
thereby rendering the process biomimetic. Detailed spectroscopy, X ray
diffraction, electron microscopy as well as fluorescence and luminescence
assisted characterization of the colloidal dispersions on day one and day seven
reveals graphene and collagen interaction and subsequent rearrangement to form
an open structure. Detailed confocal microscopy, in the liquid state, reveals
the initial attack at the zigzag edges of GR, the enhancement of auto
fluorescence and finally the opening up of graphitic stacks of GR to form near
transparent G. Atomic Force Microscopy studies prove the existence of both
collagen and graphene and the disruption of periodicity at the atomic level.
Thermal conductivity of the colloid shows a 17% enhancement for a volume
fraction of less than 0.00005 of G. Time variant increase in thermal
conductivity provides qualitative evidence for the transient exfoliation of GR
to G. The composite reveals interesting properties that could propel it as a
future material for advanced bio applications including therapeutics.


  In this article, we construct a generalization of the Blum-Fran\c{c}ois
Beta-splitting model for evolutionary trees, which was itself inspired by
Aldous' Beta-splitting model on cladograms. The novelty of our approach allows
for asymmetric shares of diversification rates (or diversification `potential')
between two sister species in an evolutionarily interpretable manner, as well
as the addition of extinction to the model in a natural way. We describe the
incremental evolutionary construction of a tree with n leaves by splitting or
freezing extant lineages through the Generating, Organizing and Deleting
processes. We then give the probability of any (binary rooted) tree under this
model with no extinction, at several resolutions: ranked planar trees giving
asymmetric roles to the first and second offspring species of a given species
and keeping track of the order of the speciation events occurring during the
creation of the tree, unranked planar trees, ranked non-planar trees and
finally (unranked non-planar) trees. We also describe a continuous-time
equivalent of the Generating, Organizing and Deleting processes where tree
topology and branch-lengths are jointly modeled and provide code in
SageMath/python for these algorithms.


  Biological organisms have evolved a wide range of immune mechanisms to defend
themselves against pathogens. Beyond molecular details, these mechanisms differ
in how protection is acquired, processed and passed on to subsequent
generations -- differences that may be essential to long-term survival. Here,
we introduce a mathematical framework to compare the long-term adaptation of
populations as a function of the pathogen dynamics that they experience and of
the immune strategy that they adopt. We find that the two key determinants of
an optimal immune strategy are the frequency and the characteristic timescale
of the pathogens. Depending on these two parameters, our framework identifies
distinct modes of immunity, including adaptive, innate, bet-hedging and
CRISPR-like immunities, which recapitulate the diversity of natural immune
systems.


  Horizontal transfer (HT) of heritable information or `traits' (carried by
genetic elements, endosymbionts, or culture) is widespread among living
organisms. Yet current ecological and evolutionary theory addressing HT is
limited. We present a modeling framework for the dynamics of two populations
that compete for resources and exchange horizontally (transfer) an otherwise
vertically inherited trait. Competition influences individual demographics,
affecting population size, which feeds back on the dynamics of transfer. We
capture this feedback with a stochastic individual-based model, from which we
derive a deterministic approximation for large populations. The interaction
between horizontal transfer and competition makes possible the stable (or
bi-stable) polymorphic maintenance of deleterious traits (including costly
plasmids). When transfer rates are of a general density-dependent form,
transfer stochasticity contributes strongly to population fluctuations. For an
initially rare trait, we describe the probabilistic dynamics of invasion and
fixation. Acceleration of fixation by HT is faster when competition is weak in
the resident population. Thus, HT can have a major impact on the distribution
of mutational effects that are fixed, and our model provides a basis for a
general theory of the influence of HT on eco-evolutionary dynamics and
adaptation.


  Directional growth caused by gravitropism and corresponding bending of plant
cells has been explored since 19th century, however, many aspects of mechanisms
underlying the perception of gravity at the molecular level are still not well
known. Perception of gravity in root and shoot gravitropisms is usually
attributed to gravisensitive cells, called statocytes, which exploit
sedimentation of macroscopic and heavy organelles, amyloplasts, to sense the
direction of gravity. Gravity stimulus is then transduced into distal
elongation zone, which is several mm far from statocytes, where it causes
stretching. It is suggested that gravity stimulus is conveyed by gradients in
auxin flux. We propose a theoretical model that may explain how concentration
gradients and/or stretching may indirectly affect the global orientation of
cortical microtubules, attached to the cell membrane and induce their dynamic
reorientation perpendicular to the gradients. In turn, oriented microtubules
arrays direct the growth and orientation of cellulose microfibrils, forming
part of the cell external skeleton and determine the shape of the cell.
Reorientation of microtubules is also observed in reaction to light in
phototropism and mechanical bending, thus suggesting universality of the
proposed mechanism.


  A fundamental problem in community ecology is to understand how ecological
processes such as selection, drift, and immigration give rise to observed
patterns in species composition and diversity. Here, we present a simple,
analytically tractable, presence-absence (PA) model for community assembly and
use it to ask how ecological traits such as the strength of competition, the
amount of diversity, and demographic and environmental stochasticity affect
species composition in a community. In the PA model, species are treated as
stochastic binary variables that can either be present or absent in a
community: species can immigrate into the community from a regional species
pool and can go extinct due to competition and stochasticity. Despite its
simplicity, the PA model reproduces the qualitative features of more
complicated models of community assembly. In agreement with recent work on
large, competitive Lotka-Volterra systems, the PA model exhibits distinct
ecological behaviors organized around a special ("critical") point
corresponding to Hubbell's neutral theory of biodiversity. These results
suggest that the concepts of ecological "phases" and phase diagrams can provide
a powerful framework for thinking about community ecology and that the PA model
captures the essential ecological dynamics of community assembly.


  Protein contacts contain important information for protein structure and
functional study, but contact prediction from sequence information remains very
challenging. Recently evolutionary coupling (EC) analysis, which predicts
contacts by detecting co-evolved residues (or columns) in a multiple sequence
alignment (MSA), has made good progress due to better statistical assessment
techniques and high-throughput sequencing. Existing EC analysis methods predict
only a single contact map for a given protein, which may have low accuracy
especially when the protein under prediction does not have a large number of
sequence homologs. Analogous to ab initio folding that usually predicts a few
possible 3D models for a given protein sequence, this paper presents a novel
structure learning method that can predict a set of diverse contact maps for a
given protein sequence, in which the best solution usually has much better
accuracy than the first one. Our experimental tests show that for many test
proteins, the best out of 5 solutions generated by our method has accuracy at
least 0.1 better than the first one when the top L/5 or L/10 (L is the sequence
length) predicted long-range contacts are evaluated, especially for protein
families with a small number of sequence homologs. Our best solutions also have
better quality than those generated by the two popular EC methods Evfold and
PSICOV.


  Rare long distance dispersal events are thought to have a disproportionate
impact on the spread of invasive species. Modelling using integrodifference
equations suggests that, when long distance contacts are represented by a
fat-tailed dispersal kernel, an accelerating wave of advance can ensue.
Invasions spreading in this manner could have particularly dramatic effects.
Recently, various authors have suggested that demographic stochasticity
disrupts wave acceleration. Integrodifference models have been widely used in
movement ecology, and as such a clearer understanding of stochastic effects is
needed. Here, we present a stochastic non-linear one-dimensional lattice model
in which demographic stochasticity and the dispersal regime can be
systematically varied. Extensive simulations show that stochasticity has a
profound effect on model behaviour, and usually breaks acceleration for
fat-tailed kernels. Exceptions are seen for some power law kernels, $K(l)
\propto |l|^{-\beta}$ with $\beta < 3$, for which acceleration persists despite
stochasticity. Such kernels lack a second moment and are important in
`accelerating' phenomena such as L\'{e}vy flights. Furthermore, for long-range
kernels the approach to the continuum limit behaviour as stochasticity is
reduced is generally slow. Given that real-world populations are finite,
stochastic models may give better predictive power when long-range dispersal is
important. Insights from mean-field models such as integrodifference equations
should be applied with caution in such circumstances.


  Gene regulatory networks present a wide variety of dynamical responses to
intrinsic and extrinsic perturbations. Arguably, one of the most important of
such coordinated responses is the one of amplification cascades, in which
activation of a few key-responsive transcription factors (termed master
regulators, MRs) lead to a large series of transcriptional activation events.
This is so since master regulators are transcription factors controlling the
expression of other transcription factor molecules and so on. MRs hold a
central position related to transcriptional dynamics and control of gene
regulatory networks and are often involved in complex feedback and feedforward
loops inducing non-trivial dynamics. Recent studies have pointed out to the
myocyte enhancing factor 2C (MEF2C, also known as MADS box transcription
enhancer factor 2, polypeptide C) as being one of such master regulators
involved in the pathogenesis of primary breast cancer. In this work, we perform
an integrative genomic analysis of the transcriptional regulation activity of
MEF2C and its target genes to evaluate to what extent are these molecules
inducing collective responses leading to gene expression deregulation and
carcinogenesis. We also analyzed a number of induced dynamic responses, in
particular those associated with transcriptional bursts, and nonlinear
cascading to evaluate the influence they may have in malignant phenotypes and
cancer.


  In this paper we review and extend our earlier recent work on thermostated
systems. A description of nano-biological systems by Markov chains in
coordinate space in the strongly overdamped limit is presented.
Characterization of the most probable path is given and a new formula for the
probability of this special path is provided from recursion formulae. The
deterministic limit is derived and the significance of Lagrange multipliers
introduced when constructing the most probable path is elucidated. The
characterization of the generation of path entropy by the most probable path is
given an equivalent interpretation relating to the rate of entropy production
by the most probable path. The paper concludes with an account of the
biological implications. Here we address why the origin of life and its
subsequent evolution took place, not the particular chemical details of how it
happened.


  Transformation of protein 3D structures from their all-atom representation
(AAR) to the double-centroid reduced representation (DCRR) is a prerequisite to
the implementation of both the tetrahedral three-dimensional search motif (3D
SM) method for predicting specific ligand binding sites (LBS) in proteins, and
the 3D interface search motif tetrahedral pair (3D ISMTP) method for predicting
binary protein-protein interaction (PPI) partners (Reyes, V.M., 2015a & c,
2015b, 2009a, b & c). Here we describe results demonstrating the efficacy of
the set of FORTRAN 77 and 90 source codes used in the transformation from AAR
to DCRR and the implementation of the 3D SM and 3D ISMTP methods. Precisely, we
show here the construction of the 3D SM for the biologically important ligands,
GTP and sialic acid, from a training set composed of experimentally solved
structures of proteins complexed with the pertinent ligand, and their
subsequent use in the screening for potential receptor proteins of the two
ligands. We also show here the construction of the 3D ISMTP for the binary
complexes, RAC:P67PHOX and KAP:phospho-CDK2, from a training set composed of
the experimentally solved complexes, and their subsequent use in the screening
for potential protomers of the two complexes. The 15 FORTRAN programs used in
the AAR to DCRR transformation and the implementation of the two said methods
are: get_bbn.f, get_sdc.f, res2cm_bbn.f, res2cm_sdc.f, nrst_ngbr.f,
find_Hbonds.f, find_VDWints.f, find_clusters.f90, find_trees.f90,
find_edgenodes.f90, match_nodes.f, fpBS.f90, Gen_Chain_Separ.f,
remove_H_atoms.f and resd_num_reduct.f. Two flowcharts - one showing how to
implement the tetrahedral 3D SM method to find LBSs in proteins, and another
how to implement the 3D ISMTP method to find binary PPI partners - are
presented in our two companion papers (Fig. 2, Reyes, V.M., 2015a, Fig. 1 & 2,
Reyes, V.M., 2015c).


  Population expansions trigger many biomedical and ecological transitions,
from tumor growth to invasions of non-native species. Although population
spreading often selects for more invasive phenotypes, we show that this outcome
is far from inevitable. In cooperative populations, mutations reducing
dispersal have a competitive advantage. Such mutations then steadily accumulate
at the expansion front bringing invasion to a halt. Our findings are a rare
example of evolution driving the population into an unfavorable state and could
lead to new strategies to combat unwelcome invaders. In addition, we obtain an
exact analytical expression for the fitness advantage of mutants with different
dispersal rates.


  Fisher developed his geometric model to support the micro-mutationalism
hypothesis which claims that small mutations are more likely to be beneficial
and therefore to contribute to evolution and adaptation. While others have
provided a general solution to the model using geometric approaches, we derive
an equivalent general solution using a probabilistic approach. Our approach to
Fisher's geometric model provides alternative intuition and interpretation of
the solution in terms of the model's parameters: for mutation to improve a
phenotype, its relative beneficial effect must be larger than the ratio of its
total effect and twice the difference between the current phenotype and the
optimal one. Our approach provides new insight into this classical model of
adaptive evolution.


  A fundamental result in metabolic pathway analysis states that every flux
mode can be decomposed into a sum of elementary modes. However, only a
decomposition without cancelations is biochemically meaningful, since a
reversible reaction cannot have different directions in the contributing
elementary modes. This essential requirement has been largely overlooked by the
metabolic pathway community.
  Indeed, every flux mode can be decomposed into elementary modes without
cancelations. The result is an immediate consequence of a theorem by
Rockafellar which states that every element of a linear subspace is a conformal
sum (a sum without cancelations) of elementary vectors (support-minimal
vectors). In this work, we extend the theorem, first to "subspace cones" and
then to general polyhedral cones and polyhedra. Thereby, we refine Minkowski's
and Carath\'eodory's theorems, two fundamental results in polyhedral geometry.
We note that, in general, elementary vectors need not be support-minimal, in
fact, they are conformally non-decomposable and form a unique minimal set of
conformal generators.
  Our treatment is mathematically rigorous, but suitable for systems
biologists, since we give self-contained proofs for our results and use
concepts motivated by metabolic pathway analysis. In particular, we study cones
defined by linear subspaces and nonnegativity conditions - like the flux cone -
and use them to analyze general polyhedral cones and polyhedra.
  Finally, we review applications of elementary vectors and conformal sums in
metabolic pathway analysis.


  Network analysis is a powerful way of modeling chromatin interactions.
Assortativity is a network property used in social sciences to identify factors
affecting how people establish social ties. We propose a new approach, using
chromatin assortativity to integrate the epigenomic landscape of a specific
cell type with its chromatin interaction network and thus investigate which
proteins or chromatin marks mediate genomic contacts. We use high-resolution
Promoter Capture Hi-C and Hi-Cap data as well as ChIA-PET data from mouse
embryonic stem cells to investigate promoter-centered chromatin interaction
networks and calculate the presence of specific epigenomic features in the
chromatin fragments constituting the nodes of the network. We estimate the
association of these features to the topology of four chromatin interaction
networks and identify features localized in connected areas of the network.
Polycomb Group proteins and associated histone marks are the features with the
highest chromatin assortativity in promoter-centred networks. We then ask which
features distinguish contacts amongst promoters from contacts between promoters
and other genomic elements. We observe higher chromatin assortativity of the
actively elongating form of RNA Polymerase 2 (RNAPII) compared to inactive
forms only in interactions between promoters and other elements. Contacts among
promoters, and between promoters and other elements have different
characteristic epigenomic features. We identify a possible role for the
elongating form of RNAPII in mediating interactions among promoters, enhancers
and transcribed gene bodies. Our approach facilitates the study of multiple
genome-wide epigenomic profiles, considering network topology and allowing the
comparison of chromatin interaction networks.


  Conventional kinesin walks by a hand-over-hand mechanism on the microtubule
(MT) by taking $\sim$ 8$nm$ discrete steps, and consumes one ATP molecule per
step. The time needed to complete a single step is on the order of twenty
microseconds. We show, using simulations of a coarse-grained model of the
complex containing the two motor heads, the MT, and the coiled coil that in
order to obtain quantitative agreement with experiments for the stepping
kinetics hydrodynamic interactions (HI) have to be included. In simulations
without hydrodynamic interactions spanning nearly twenty $\mu$s not a single
step was completed in hundred trajectories. In sharp contrast, nearly 14\% of
the steps reached the target binding site within 6$\mu$s when HI were included.
Somewhat surprisingly, there are qualitative differences in the diffusion
pathways in simulations with and without HI. The extent of movement of the
trailing head of kinesin on the MT during the diffusion stage of stepping is
considerably greater in simulations with HI than in those without HI. Our
results suggest that inclusion of HI is crucial in the accurate description of
motility of other motors as well.


  This paper is concerned with stochastic SIR and SEIR epidemic models on
random networks in which individuals may rewire away from infected neighbors at
some rate $\omega$ (and reconnect to non-infectious individuals with
probability $\alpha$ or else simply drop the edge if $\alpha=0$), so-called
preventive rewiring. The models are denoted SIR-$\omega$ and SEIR-$\omega$, and
we focus attention on the early stages of an outbreak, where we derive
expression for the basic reproduction number $R_0$ and the expected degree of
the infectious nodes $E(D_I)$ using two different approximation approaches. The
first approach approximates the early spread of an epidemic by a branching
process, whereas the second one uses pair approximation. The expressions are
compared with the corresponding empirical means obtained from stochastic
simulations of SIR-$\omega$ and SEIR-$\omega$ epidemics on Poisson and
scale-free networks. Without rewiring of exposed nodes, the two approaches
predict the same epidemic threshold and the same $E(D_I)$ for both types of
epidemics, the latter being very close to the mean degree obtained from
simulated epidemics over Poisson networks. Above the epidemic threshold,
pairwise models overestimate the value of $R_0$ computed from simulations,
which turns out to be very close to the one predicted by the branching process
approximation. When exposed individuals also rewire with $\alpha > 0$ (perhaps
unaware of being infected), the two approaches give different epidemic
thresholds, with the branching process approximation being more in agreement
with simulations.


  Microbial communities play important roles in the function and maintenance of
various biosystems, ranging from human body to the environment. Current methods
for analysis of microbial communities are typically based on taxonomic
phylogenetic alignment using 16S rRNA metagenomic or Whole Genome Sequencing
data. In typical characterizations of microbial communities, studies deal with
billions of micobial sequences, aligning them to a phylogenetic tree. We
introduce a new approach for the efficient analysis of microbial communities.
Our new reference-free analysis tech- nique is based on n-gram sequence
analysis of 16S rRNA data and reduces the processing data size dramatically (by
105 fold), without requiring taxonomic alignment. The proposed approach is
applied to characterize phenotypic microbial community differ- ences in
different settings. Specifically, we applied this approach in classification of
microbial com- munities across different body sites, characterization of oral
microbiomes associated with healthy and diseased individuals, and
classification of microbial communities longitudinally during the develop- ment
of infants. Different dimensionality reduction methods are introduced that
offer a more scalable analysis framework, while minimizing the loss in
classification accuracies. Among dimensionality re- duction techniques, we
propose a continuous vector representation for microbial communities, which can
widely be used for deep learning applications in microbial informatics.


  Ligand burial depth is an indicator of protein flexibility, as the extent of
receptor conformational change required to bind a ligand in general varies
directly with its depth of burial. In a companion paper (Reyes, V.M. 2015a), we
report on the Tangent Sphere (TS) and Cutting Plane (CP) methods --
complementary methods to quantify, independent of protein size, the degree of
ligand burial in a protein receptor. In this report, we present results that
demonstrate the effectiveness of a set of FORTRAN 77 and 90 source codes used
in the implementation of the two related procedures, as well as the precise
implementation of the procedures. Particularly, we show here that application
of the TS and CP methods on a theoretical model protein in the form of a
spherical grid of points accurately portrays the behavior of the TS and CP
indices, the predictive parameters obtained from the two methods. We also show
that results of the implementation of the TS and CP methods on six protein
receptors (Laskowski et al. 1996) are inagreement with their findings regarding
cavity sizes in these proteins. The six FORTRAN programs we present here are:
find_molec_centr.f, tangent_sphere.f, find_CP_coeffs.f, CPM_Neg_Side.f,
CPM_Pos_Side.f and CPM_Zero_Side.f. The first program calculates the x-, y- and
z-coordinates of the molecular geometric centroid of the protein (global
centroid, GC), the center of the TS. Its radius is the distance between the GC
and the local centroid (LC), the centroid of the bound ligand or a portion of
its binding site. The second program finds the number of protein atoms inside,
outside and on the TS. The third determines the four coefficients A, B, C and D
of the equation of the CP, Ax + By + Cz + D = 0. The CP is tangent to the TS at
GC. The fourth, fifth and sixth programs determine the number of protein atoms
lying on the negative side, positive side, and on the CP.


  We previously described the representation of protein 3D structures in
spherical coordinates (rho, phi, theta) and two of its applications: separation
of the outer layer (OL) from the inner core (IC) of proteins, and assessment of
protein surface protrusions and invaginations (Reyes, V.M., 2011& 2009). Here
we present results demonstrating the performance success of the FORTRAN 77 and
90 programs used in the implementation of the two said applications, and how to
implement both applications. In particular, we show here data that demonstrate
the success of our OL-IC separation procedure using a subset of the Laskowski
et al. (1996) dataset. Using a theoretical model protein in the form of a
scalene ellipsoid grid of points with and without an artificially constructed
protrusion or invagination, we also show results demonstrating that protrusions
and invaginations on the protein surface maybe predicted. The nine programs we
present here and their respective functions are: find_molec_centr.f: finds the
x-, y- and z-coordinates of the protein molecular geometric centroid,
cart2sphere_degrees.f90: converts PDB protein coordinates to spherical, with
phi and theta in degrees, cart2sphere_radians.f90: does the same thing as the
second program, but with phi and theta in radians, spher2cart_degrees.f90:
converts the coordinates from spherical back to PDB, where input phi and theta
are in degrees, spher2cart_radians.f90: does the same thing as the fourth
program, but with phi and theta in radians, find_rho_cutoff.f: determines the
rho cut-off for finding the boundary between OL and IC,
phi6_theta8_binning.f90: performs the binning of phi in six- and theta in
eight-degree increments, phi10_theta10_binning.f90: performs the binning of phi
and theta both in ten-degree increments, and bin_rho.f90: performs the binning
of rho values for plotting the frequency distribution of maximum rho values.


  We study the spreading of two mutually cooperative diseases on different
network topologies, and with two microscopic realizations, both of which are
stochastic versions of an SIR type model studied by us recently in mean field
approximation. There it had been found that cooperativity can lead to
first-order spreading/extinction transitions. However, due to the rapid mixing
implied by the mean field assumption, first order transitions required non-zero
initial densities of sick individuals. For the stochastic model studied here
the results depend strongly on the underlying network. First order transitions
are found when there are few short but many long loops: (i) No first order
transitions exist on trees and on 2-d lattices with local contacts (ii) They do
exist on Erdos-Renyi (ER) networks, on d-dimensional lattices with d >= 4, and
on 2-d lattices with sufficiently long-ranged contacts; (iii) On 3-d lattices
with local contacts the results depend on the microscopic details of the
implementation; (iv) While single infected seeds can always lead to infinite
epidemics on regular lattices, on ER networks one sometimes needs finite
initial densities of infected nodes; (v) In all cases the first order
transitions are actually "hybrid", i.e. they display also power law scaling
usually associated with second order transitions. On regular lattices, our
model can also be interpreted as the growth of an interface due to cooperative
attachment of two species of particles. Critically pinned interfaces in this
model seem to be in different universality classes than standard critically
pinned interfaces in models with forbidden overhangs. Finally, the detailed
results mentioned above hold only when both diseases propagate along the same
network of links. If they use different links, results can be rather different
in detail, but are similar overall.


  We consider the dynamics of a population of organisms containing two mutually
inhibitory gene regulatory networks, that can result in a bistable switch-like
behaviour. We completely characterize their local and global dynamics in the
absence of any noise, and then go on to consider the effects of either noise
coming from bursting (transcription or translation), or Gaussian noise in
molecular degradation rates when there is a dominant slow variable in the
system. We show analytically how the steady state distribution in the
population can range from a single unimodal distribution through a bimodal
distribution and give the explicit analytic form for the invariant stationary
density which is globally asymptotically stable. Rather remarkably, the
behaviour of the stationary density with respect to the parameters
characterizing the molecular behaviour of the bistable switch is qualitatively
identical in the presence of noise coming from bursting as well as in the
presence of Gaussian noise in the degradation rate. This implies that one
cannot distinguish between either the dominant source or nature of noise based
on the stationary molecular distribution in a population of cells. We finally
show that the switch model with bursting but two dominant slow genes has an
asymptotically stable stationary density.


  The classical game of rock-paper-scissors have inspired experiments and
spatial model systems that address robustness of biological diversity. In
particular the game nicely illustrates that cyclic interactions allow multiple
strategies to coexist for long time intervals. When formulated in terms of a
one-dimensional cellular automata, the spatial distribution of strategies
exhibits coarsening with algebraically growing domain size over time, while the
two-dimensional version allows domains to break and thereby opens for long-time
coexistence. We here consider a quasi-one-dimensional implementation of the
cyclic competition, and study the long-term dynamics as a function of rare
invasions between parallel linear ecosystems. We find that increasing the
complexity from two to three parallel subsystems allows a transition from
complete coarsening to an active steady state where the domain size stays
finite. We further find that this transition happens irrespective of whether
the update is done in parallel for all sites simultaneously, or done randomly
in sequential order. In both cases the active state is characterized by
localized bursts of dislocations, followed by longer periods of coarsening. In
the case of the parallel dynamics, we find that there is another phase
transition between the active steady state and the coarsening state within the
three-line system when the invasion rate between the subsystems is varied. We
identify the critical parameter for this transition, and show that the density
of active boundaries have critical exponents that are consistent with the
directed percolation universality class. On the other hand, numerical
simulations with the random sequential dynamics suggest that the system may
exhibit an active steady state as long as the invasion rate is finite.


  Protein secondary structure (SS) prediction is important for studying protein
structure and function. When only the sequence (profile) information is used as
input feature, currently the best predictors can obtain ~80% Q3 accuracy, which
has not been improved in the past decade. Here we present DeepCNF (Deep
Convolutional Neural Fields) for protein SS prediction. DeepCNF is a Deep
Learning extension of Conditional Neural Fields (CNF), which is an integration
of Conditional Random Fields (CRF) and shallow neural networks. DeepCNF can
model not only complex sequence-structure relationship by a deep hierarchical
architecture, but also interdependency between adjacent SS labels, so it is
much more powerful than CNF. Experimental results show that DeepCNF can obtain
~84% Q3 accuracy, ~85% SOV score, and ~72% Q8 accuracy, respectively, on the
CASP and CAMEO test proteins, greatly outperforming currently popular
predictors. As a general framework, DeepCNF can be used to predict other
protein structure properties such as contact number, disorder regions, and
solvent accessibility.


  We present an efficient and flexible method for computing likelihoods of
phenotypic traits on a phylogeny. The method does not resort to Monte-Carlo
computation but instead blends Felsenstein's discrete character pruning
algorithm with methods for numerical quadrature. It is not limited to Gaussian
models and adapts readily to model uncertainty in the observed trait values. We
demonstrate the framework by developing efficient algorithms for likelihood
calculation and ancestral state reconstruction under Wright's threshold model,
applying our methods to a dataset of trait data for extrafloral nectaries
(EFNs) across a phylogeny of 839 Labales species.


  A theoretical model was developed to investigate the rod phototransduction
process in the mouse. In particular, we explored the biochemical reactions of
several chemical components that contribute to the signaling process
into/around the membranous disks in the outer segments of the rod cells. We
constructed a rate equation model incorporating the molecular crowding effects
of rhodopsin according to experimental results, which may hinder the diffusion
of molecules on the disk mem- brane. The present model could effectively
reproduce and explain the mechanisms of the following phenomena observed in
experiments. First, the activations and relaxation of the wild-type mouse rod
cell progressed more slowly than those of mutant cells containing half the
amount of rhodopsin on the disk membrane. Second, the strong photoactivated
state of the cell was sustained for a longer period when the light stimuli were
strong. Finally, the lifetime of photoactivation exhibited a logarithmic
increase with increasing light strength given exposure to strong light stimuli.


  The widely used genetic pleiotropic analysis of multiple phenotypes are often
designed for examining the relationship between common variants and a few
phenotypes. They are not suited for both high dimensional phenotypes and high
dimensional genotype (next-generation sequencing) data. To overcome these
limitations, we develop sparse structural equation models (SEMs) as a general
framework for a new paradigm of genetic analysis of multiple phenotypes. To
incorporate both common and rare variants into the analysis, we extend the
traditional multivariate SEMs to sparse functional SEMs. To deal with high
dimensional phenotype and genotype data, we employ functional data analysis and
the alternative direction methods of multiplier (ADMM) techniques to reduce
data dimension and improve computational efficiency. Using large scale
simulations we showed that the proposed methods have higher power to detect
true causal genetic pleiotropic structure than other existing methods.
Simulations also demonstrate that the gene-based pleiotropic analysis has
higher power than the single variant-based pleiotropic analysis. The proposed
method is applied to exome sequence data from the NHLBI Exome Sequencing
Project (ESP) with 11 phenotypes, which identifies a network with 137 genes
connected to 11 phenotypes and 341 edges. Among them, 114 genes showed
pleiotropic genetic effects and 45 genes were reported to be associated with
phenotypes in the analysis or other cardiovascular disease (CVD) related
phenotypes in the literature.