  The origin and shape of metabolic scaling has been controversial since
Kleiber found that basal metabolic rate of animals seemed to vary as a power
law of their body mass with exponent 3/4, instead of 2/3, as a
surface-to-volume argument predicts. The universality of exponent 3/4 -claimed
in terms of the fractal properties of the nutrient network- has recently been
challenged according to empirical evidence that observed a wealth of robust
exponents deviating from 3/4. Here we present a conceptually simple
thermodynamic framework, where the dependence of metabolic rate with body mass
emerges from a trade-off between the energy dissipated as heat and the energy
efficiently used by the organism to maintain its metabolism. This balance tunes
the shape of an additive model from which different effective scalings can be
recovered as particular cases, thereby reconciling previously inconsistent
empirical evidence in mammals, birds, insects and even plants under a unified
framework. This model is biologically motivated, fits remarkably well the data,
and also explains additional features such as the relation between energy lost
as heat and mass, the role and influence of different climatic environments or
the difference found between endotherms and ectotherms.


  Phenomenological screening of small molecule libraries for anticancer
activity yields potentially interesting candidate molecules, with a bottleneck
in the determination of drug targets and the mechanism of anticancer action. A
novel approach to drug target deconvolution compares the abundance profiles of
proteins expressed in a panel of cells treated with different drugs, and
identifies proteins with cell-type independent and drug-specific regulation
that is exceptionally strong in relation to the other proteins. Mapping top
candidates on known protein networks reveals the mechanism of drug action,
while abundant proteins provide a signature of cellular death/survival
pathways. The above approach can significantly shorten drug target
identification, and thus facilitate the emergence of novel anticancer
treatments.


  This study analyzes rates and ways of data sharing regarding mitochondrial, Y
chromosomal and autosomal polymorphisms in a total of 162 papers on human
ancient DNA published between 1988 and 2013. For the most part, data are
available in such a way as to make their scrutiny and reuse possible. The
estimated sharing rate is not far from totality (97.6% +/- 2.1%) and
substantially higher than observed in other fields of genetic research
(Evolutionary, Medical and Forensic Genetics). A questionnaire-based survey
suggests that the authors awareness of the importance of openness and
transparency for scientific progress is a fundamental factor for the
achievement of such a high sharing rate. Most data were made available through
body text, but the use of primary databases increased with the application of
complete mitochondrial and next generation sequencing methods. Our study
highlights three important aspects. First, we provide evidence that researchers
motivations are as necessary as stakeholders policies and norms to achieve very
high sharing rates. Second, careful analyses of the ways in which data are made
available are an important first step to maximize data findability,
accessibility, useability and preservation. Third and finally, the case of
human ancient DNA studies demonstrates how Open Science can foster scientific
advancements, showing that openness and transparency can help build rigorous
and reliable scientific practices even in the presence of complex experimental
challenges.


  The biomolecules in and around a living cell -- proteins, nucleic acids,
lipids, carbohydrates -- continuously sample myriad conformational states that
are thermally accessible at physiological temperatures. Simultaneously, a given
biomolecule also samples (and is sampled by) a rapidly fluctuating local
environment comprised of other biopolymers, small molecules, water, ions, etc.
that diffuse to within a few nanometers, leading to inter-molecular contacts
that stitch together large supramolecular assemblies. Indeed, all biological
systems can be viewed as dynamic networks of molecular interactions. As a
complement to experimentation, molecular simulation offers a uniquely powerful
approach to analyze biomolecular structure, mechanism, and dynamics; this is
possible because the molecular contacts that define a complicated biomolecular
system are governed by the same physical principles (forces, energetics) that
characterize individual small molecules, and these simpler systems are
relatively well-understood. With modern algorithms and computing capabilities,
simulations are now an indispensable tool for examining biomolecular assemblies
in atomic detail, from the conformational motion in an individual protein to
the diffusional dynamics and inter-molecular collisions in the early stages of
formation of cellular-scale assemblies such as the ribosome. This text
introduces the physicochemical foundations of molecular simulations and
docking, largely from the perspective of biomolecular interactions.


  Cycles in population dynamics are widely found in nature. These cycles are
understood as emerging from the interaction between two or more coupled
species. Here, we argue that data regarding population dynamics are prone to
misinterpretation when sampling is conducted at a slow rate compared to the
population cycle period. This effect, known as aliasing, is well described in
other areas, such as signal processing and computer graphics. However, to the
best of our knowledge, aliasing has never been addressed in the population
dynamics context or in coupled oscillatory systems. To illustrate aliasing, the
Lotka-Volterra model oscillatory regime is numerically sampled, creating
prey-predator cycles. Inadequate sampling periods produce inversions in the
cause-effect relationship and an increase in cycle period, as reported in the
well-known hare-lynx paradox. More generally, slow acquisition rates may
distort data, producing deceptive patterns and eventually leading to data
misinterpretation.


  Viruses are nanoscale entities containing a nucleic acid genome encased in a
protein shell called a capsid, and in some cases surrounded by a lipid bilayer
membrane. This review summarizes the physics that govern the processes by which
capsids assembles within their host cells and in vitro. We describe the
thermodynamics and kinetics for assembly of protein subunits into icosahedral
capsid shells, and how these are modified in cases where the capsid assembles
around a nucleic acid or on a lipid bilayer. We present experimental and
theoretical techniques that have been used to characterize capsid assembly, and
we highlight aspects of virus assembly which are likely to receive significant
attention in the near future.


  The Yule-Harding-Kingman (YHK) model and the proportional to distinguishable
arrangements (PDA) model are two binary tree generating models that are widely
used in evolutionary biology. Understanding the distributions of clade sizes
under these two models provides valuable insights into macro-evolutionary
processes, and is important in hypothesis testing and Bayesian analyses in
phylogenetics. Here we show that these distributions are log-convex, which
implies that very large clades or very small clades are more likely to occur
under these two models. Moreover, we prove that there exists a critical value
$\kappa(n)$ for each $n\geqslant 4$ such that for a given clade with size $k$,
the probability that this clade is contained in a random tree with $n$ leaves
generated under the YHK model is higher than that under the PDA model if
$1<k<\kappa(n)$, and lower if $\kappa(n)<k<n$. Finally, we extend our results
to binary unrooted trees, and obtain similar results for the distributions of
clan sizes.


  We develop a new methodology for the efficient computation of epidemic final
size distributions for a broad class of Markovian models. We exploit a
particular representation of the stochastic epidemic process to derive a method
which is both computationally efficient and numerically stable. The algorithms
we present are also physically transparent and so allow us to extend this
method from the basic SIR model to a model with a phase-type infectious period
and another with waning immunity. The underlying theory is applicable to many
Markovian models where we wish to efficiently calculate hitting probabilities.


  Genomic subsequences conserved between closely related species such as human
and chimpanzee exhibit an exponential length distribution, in contrast to the
algebraic length distribution observed for sequences shared between distantly
related genomes. We find that the former exponential can be further decomposed
into an exponential component primarily composed of orthologous sequences, and
a truncated algebraic component primarily composed of paralogous sequences.


  We explore the interplay between tunneling process and chiral interactions in
the discrimination of chiral states for an ensemble of molecules in a
biological environment. Each molecule is described by an asymmetric double-well
potential and the environment is modeled as a bath of harmonic oscillators. We
carefully analyze different time-scales appearing in the resulting master
equation at both weak- and strong-coupling limits. The corresponding results
are accompanied by a set of coupled differential equations characterizing
optical activity of the molecules. We show that, at the weak-coupling limit,
chiral interactions prohibit the coherent racemization induced by decoherence
effects and thus preserve the initial chiral state. At the strong-coupling
limit, considering the memory effects of the environment, Markovian behavior is
observed at long times.


  Flavescence dor\'ee (FD) is a disease of grapevine transmitted by an insect
vector, $Scaphoideus$ $titanus$ Ball. At present, no prophylaxis exists, so
mandatory control procedures (e.g. removal of infected plants, and insecticidal
sprays to avoid transmission) are in place in Italy and other European
countries. We propose a model of the epidemiology of FD by taking into account
the different aspects involved into the transmission process (acquisition of
the disease, latency and expression of symptoms, recovery rate, removal and
replacement of infected plants, insecticidal treatments, and the effect of
hotbeds). The model was constructed as a system of first order nonlinear ODEs
in four compartment variables. We perform a bifurcation analysis of the
equilibria of the model using the severity of the hotbeds as the control
parameter. Depending on the non-dimensional grapevine density of the vineyard
we find either a single family of equilibria in which the health of the
vineyard gradually deteriorates for progressively more severe hotbeds, or
multiple equilibria that give rise to sudden transitions from a nearly healthy
vineyard to a severely deteriorated one when the severity of the hotbeds
crosses a critical value. These results suggest some lines of intervention for
limiting the spread of the disease.


  The chemical reaction network (CRN) is a widely used formalism to describe
macroscopic behavior of chemical systems. Available tools for CRN modelling and
simulation require local access, installation, and often involve local file
storage, which is susceptible to loss, lacks searchable structure, and does not
support concurrency. Furthermore, simulations are often single-threaded, and
user interfaces are non-trivial to use. Therefore there are significant hurdles
to conducting efficient and collaborative chemical research. In this paper, we
introduce a new enterprise chemistry simulation framework, COEL, which
addresses these issues. COEL is the first web-based framework of its kind. A
visually pleasing and intuitive user interface, simulations that run on a large
computational grid, reliable database storage, and transactional services make
COEL ideal for collaborative research and education. COEL's most prominent
features include ODE-based simulations of chemical reaction networks and
multicompartment reaction networks, with rich options for user interactions
with those networks. COEL provides DNA-strand displacement transformations and
visualization (and is to our knowledge the first CRN framework to do so), GA
optimization of rate constants, expression validation, an application-wide
plotting engine, and SBML/Octave/Matlab export. We also present an overview of
the underlying software and technologies employed and describe the main
architectural decisions driving our development. COEL is available at
http://coel-sim.org for selected research teams only. We plan to provide a part
of COEL's functionality to the general public in the near future.


  One need only compare the number of three-dimensional molecular illustrations
in the first (1990) and third (2004) editions of Voet & Voet's "Biochemistry"
in order to appreciate this field's profound communicative value in modern
biological sciences -- ranging from medicine, physiology, and cell biology, to
pharmaceutical chemistry and drug design, to structural and computational
biology. The clich\'e about a picture being worth a thousand words is quite
poignant here: The information 'content' of an effectively-constructed piece of
molecular graphics can be immense. Because biological function arises from
structure, it is difficult to overemphasize the utility of visualization and
graphics in molding our current understanding of the molecular nature of
biological systems. Nevertheless, creating effective molecular graphics is not
easy -- neither conceptually, nor in terms of effort required. The present
collection of Rules is meant as a guide for those embarking upon their first
molecular illustrations.


  More than 1000 complete or near-complete mitochondrial DNA (mtDNA) sequences
have been deposited in GenBank for eight common domestic animals (i.e. cattle,
dog, goat, horse, pig, sheep, yak and chicken) and their close wild ancestors
or relatives. Nevertheless, few efforts have been performed to evaluate the
sequence data quality, which heavily impact the original conclusion. Herein, we
conducted a phylogenetic survey of these complete or near-complete mtDNA
sequences based on mtDNA haplogroup trees for the eight animals. We show that,
errors due to artificial recombination, surplus of mutations, and phantom
mutations, do exist in 14.5% (194/1342) of mtDNA sequences and shall be treated
with wide caution. We propose some caveats for mtDNA studies of domestic
animals in the future.


  Dominance hierarchy among animals is widespread in various species and
believed to serve to regulate resource allocation within an animal group.
Unlike small groups, however, detection and quantification of linear hierarchy
in large groups of animals are a difficult task. Here, we analyse
aggression-based dominance hierarchies formed by worker ants in Diacamma sp. as
large directed networks. We show that the observed dominance networks are
perfect or approximate directed acyclic graphs, which are consistent with
perfect linear hierarchy. The observed networks are also sparse and random but
significantly different from networks generated through thinning of the perfect
linear tournament (i.e., all individuals are linearly ranked and dominance
relationship exists between every pair of individuals). These results pertain
to global structure of the networks, which contrasts with the previous studies
inspecting frequencies of different types of triads. In addition, the
distribution of the out-degree (i.e., number of workers that the focal worker
attacks), not in-degree (i.e., number of workers that attack the focal worker),
of each observed network is right-skewed. Those having excessively large
out-degrees are located near the top, but not the top, of the hierarchy. We
also discuss evolutionary implications of the discovered properties of
dominance networks.


  Many different types of cellular cargos are transported bidirectionally along
microtubules by teams of molecular motors. The motion of this cargo-motors
system has been experimentally characterized in vivo as processive with rather
persistent directionality. Different theoretical approaches have been suggested
in order to explore the origin of this kind of motion. An effective theoretical
approach, introduced by M\"uller et al., describes the cargo dynamics as a
tug-of-war between different kinds of motors. An alternative approach has been
suggested recently by Kunwar et al., who considered the coupling between motor
and cargo in more detail. Based on this framework we introduce a model
considering single motor positions which we propagate in continuous time.
Furthermore, we analyze the possible influence of the discrete time update
schemes used in previous publications on the system's dynamic.


  Ecological networks allow us to study the structure and function of
ecosystems and gain insights on species resilience/stability. The study of this
ecological networks is usually a snapshop focused in a limited specific range
of space and time, prevent us to perceive the real dynamics of ecological
processes. By definition, an alien species has some ecological strategies and
traits that permit it to compete better than the native species (e.g. absence
of predators, different bloom period, high grow rate, etc.). Plant-pollinator
networks provide valuable services to whole ecosystems and the introduction of
an alien species may have different effects on the native network (competitive
facilitation, native species extinction, etc.). While scientists acknowledge
the significance of network connectivity in driving ecosystem services, the
inclusion of temporary networks in ecological models is still in its infancy.
We propose to use existing data on seasonality to develop a simulation platform
that show inference between temporality of networks and invasions traits. Our
focus is only to pick up some simple model to show, that theoretically temporal
aspect play a role (different extinction patterns) to encourage ecologist to
get involved in temporal networks. Moreover, the derived simulations could be
further extended and adjust to other ecological questions.


  Target identification, one of the steps of drug discovery, aims at
identifying biomolecules whose function should be therapeutically altered in
order to cure the considered pathology. This work proposes an algorithm for in
silico target identification using Boolean network attractors. It assumes that
attractors of dynamical systems, such as Boolean networks, correspond to
phenotypes produced by the modeled biological system. Under this assumption,
and given a Boolean network modeling a pathophysiology, the algorithm
identifies target combinations able to remove attractors associated with
pathological phenotypes. It is tested on a Boolean model of the mammalian cell
cycle bearing a constitutive inactivation of the retinoblastoma protein, as
seen in cancers, and its applications are illustrated on a Boolean model of
Fanconi anemia. The results show that the algorithm returns target combinations
able to remove attractors associated with pathological phenotypes and then
succeeds in performing the proposed in silico target identification. However,
as with any in silico evidence, there is a bridge to cross between theory and
practice, thus requiring it to be used in combination with wet lab experiments.
Nevertheless, it is expected that the algorithm is of interest for target
identification, notably by exploiting the inexpensiveness and predictive power
of computational approaches to optimize the efficiency of costly wet lab
experiments.


  Thresholds occur in the dynamics of many biological communities. Here we
model a persistence type threshold which has been shown experimentally to exist
in hyperparasitised flukes in the cockle, a shellfish. Our model consists of a
periodically driven slow-fast host-parasite system of equations for a slow
flukes population and a fast Unikaryon hyperparasite population. The model
exhibits two branches of the critical curve crossing in a transcritical
bifurcation scenario. We discuss two thresholds due to immediate and delayed
exchange of stability effects; and we derive algebraic relationships for
parameters of the periodic solution in the limit of the infinite ratio of the
time scales. Flukes parasitise cockles and in turn are hyperparasitised by the
microsporidian Unikaryon legeri; the life cycle of flukes includes several life
stages and a number of different hosts. That is, the flukes-hyperparasite
system in a cockle is part of a larger estuarine ecosystem of species involving
parasites, shellfish and birds which prey on shellfish. A population dynamics
model which accounts for such multi-species interactions and includes the
fluke-hyperparasite model in a cockle as a subsystem is presented. We provide
evidence that the threshold effect we observed in the flukes-hyperparasite
subsystem remains apparent in the multi-species system. Assuming that flukes
damage cockles, and taking into account that the hyperparasite is detrimental
to flukes, it is natural to suggest that the hyperparasitism may support the
abundance of cockles and, thereby, the persistence of the ecosystem, including
shellfish and birds. We confirm the possibility of this scenario in our model
by removing the hyperparasite and demonstrating that this may result in a
substantial drop in cockle numbers. The result indicates a possible significant
role for the microparasite in this estuarine ecosystem.


  The determination of a 'folding core' can help to provide insight into the
structure, flexibility, mobility and dynamics, and hence, ultimately, function
of a protein - a central concern of structural biology. Changes in the folding
core upon ligand binding are of particular interest because they may be
relevant to drug-induced functional changes. Cyclophilin A is a
multi-functional ligand-binding protein and a significant drug target. It acts
principally as an enzyme during protein folding, but also as the primary
binding partner for the immunosuppressant drug cyclosporin A (CsA). Here, we
have used hydrogen-deuterium exchange (HDX) NMR spectroscopy to determine the
folding core of the CypA-CsA complex. We also use the rapid computational tool
of rigidity analysis, implemented in FIRST, to determine a theoretical folding
core of the complex. In addition we generate a theoretical folding core for the
unbound protein and compare this with previously published HDX data. The FIRST
method gives a good prediction of the HDX folding core, but we find that it is
not yet sufficiently sensitive to predict the effects of ligand binding on
CypA.


  Determining the folding core of a protein yields information about its
folding process and dynamics. The experimental procedures for identifying the
amino acids which make up the folding core include hydrogen-deuterium exchange
and $\Phi$-value analysis and can be expensive and time consuming. As such
there is a desire to improve upon existing methods for determining protein
folding cores theoretically. Here, we use a combined method of rigidity
analysis alongside coarse-grained simulations of protein motion in order to
improve folding core predictions for unbound CypA and for the CypA-CsA complex.
We find that the most specific prediction of folding cores in CypA and CypA-CsA
comes from the intersection of the results of static rigidity analysis,
implemented in the FIRST software suite, and simulations of the propensity for
flexible motion, using the FRODA tool.


  In subdivided populations, migration acts together with selection and genetic
drift and determines their evolution. Building up on a recently proposed
method, which hinges on the emergence of a time scale separation between local
and global dynamics, we study the fixation properties of subdivided populations
in the presence of balancing selection. The approximation implied by the method
is accurate when the effective selection strength is small and the number of
subpopulations is large. In particular, it predicts a phase transition between
species coexistence and biodiversity loss in the infinite-size limit and, in
finite populations, a nonmonotonic dependence of the mean fixation time on the
migration rate. In order to investigate the fixation properties of the
subdivided population for stronger selection, we introduce an effective coarser
description of the dynamics in terms of a voter model with intermediate states,
which highlights the basic mechanisms driving the evolutionary process.


  Linking networks of molecular interactions to cellular functions and
phenotypes is a key goal in systems biology. Here, we adapt concepts of spatial
statistics to assess the functional content of molecular networks. Based on the
guilt-by-association principle, our approach (called SANTA) quantifies the
strength of association between a gene set and a network, and functionally
annotates molecular networks like other enrichment methods annotate lists of
genes. As a general association measure, SANTA can (i) functionally annotate
experimentally derived networks using a collection of curated gene sets, and
(ii) annotate experimentally derived gene sets using a collection of curated
networks, as well as (iii) prioritize genes for follow-up analyses. We
exemplify the efficacy of SANTA in several case studies using the \emph{S.
cerevisiae} genetic interaction network and genome-wide RNAi screens in cancer
cell lines. Our theory, simulations and applications show that SANTA provides a
principled statistical way to quantify the association between molecular
networks and cellular functions and phenotypes. SANTA is available from
http://bioconductor.org/packages/release/bioc/html/SANTA.html.


  Motivated by the wide range of known self-replicating systems, some far from
genetics, we study a system composed by individuals having an internal dynamics
with many possible states that are partially stable, with varying mutation
rates. Individuals reproduce and die with a rate that is a property of each
state, not necessarily related to its stability, and the offspring is born on
the parent's state. The total population is limited by resources or space, as
for example in a chemostat or a Petri dish. Our aim is to show that mutation
rate and fitness become more correlated, \emph{even if they are completely
uncorrelated for an isolated individual}, underlining the fact that the
interaction induced by limitation of resources is by itself efficient for
generating collective effects.


  Modulation of the properties of AMPA receptors at the post-synaptic membrane
is one of the main suggested mechanisms behind synaptic plasticity in the
central nervous system of vertebrates. Electrophysiological recordings of
single channels stimulated with agonists showed that both recombinant and
native AMPA receptors visit multiple conductance states in an agonist
concentration dependent manner. We propose an allosteric model of the multiple
conductance states based on concerted conformational transitions of the four
subunits, as an iris diaphragm. Our model predicts that the thermodynamic
behaviour of the conductance states upon full and partial agonist stimulations
can be described with increased affinity of receptors as they progress to
higher conductance states. The model also predicts existence of AMPA receptors
in non-liganded conductive substates. However, spontaneous openings probability
decreases with increasing conductances. Finally, we predict that the large
conductance states are stabilized within the rise phase of a whole-cell EPSC in
glutamatergic hippocampal neurons. Our model provides a mechanistic link
between ligand concentration and conductance states that can explain
thermodynamic and kinetic features of AMPA receptor gating.


  It has been recently observed that the dynamical properties of mass action
systems arising from many models of biochemical reaction networks can be
derived by considering the corresponding properties of a related generalized
mass action system. The correspondence process known as network translation in
particular has been shown to be useful in characterizing a system's steady
states. In this paper, we further develop the theory of network translation
with particular focus on a subclass of translations known as improper
translations. For these translations, we derive conditions on the network
topology of the translated network which are sufficient to guarantee the
original and translated systems share the same steady states. We then present a
mixed-integer linear programming (MILP) algorithm capable of determining
whether a mass action system can be corresponded to a generalized system
through the process of network translation.


  Biological network alignment identifies topologically and functionally
conserved regions between networks of different species. It encompasses two
algorithmic steps: node cost function (NCF), which measures similarities
between nodes in different networks, and alignment strategy (AS), which uses
these similarities to rapidly identify high-scoring alignments. Different
methods use both different NCFs and different ASs. Thus, it is unclear whether
the superiority of a method comes from its NCF, its AS, or both. We already
showed on MI-GRAAL and IsoRankN that combining NCF of one method and AS of
another method can lead to a new superior method. Here, we evaluate MI-GRAAL
against newer GHOST to potentially further improve alignment quality. Also, we
approach several important questions that have not been asked systematically
thus far. First, we ask how much of the node similarity information in NCF
should come from sequence data compared to topology data. Existing methods
determine this more-less arbitrarily, which could affect the resulting
alignment(s). Second, when topology is used in NCF, we ask how large the size
of the neighborhoods of the compared nodes should be. Existing methods assume
that larger neighborhood sizes are better.
  We find that MI-GRAAL's NCF is superior to GHOST's NCF, while the performance
of the methods' ASs is data-dependent. Thus, the combination of MI-GRAAL's NCF
and GHOST's AS could be a new superior method for certain data. Also, which
amount of sequence information is used within NCF does not affect alignment
quality, while the inclusion of topological information is crucial. Finally,
larger neighborhood sizes are preferred, but often, it is the second largest
size that is superior, and using this size would decrease computational
complexity.
  Together, our results give several general recommendations for a fair
evaluation of network alignment methods.


  Network self-similarity or fractality are widely accepted as an important
topological property of metabolic networks; however, recent studies cast doubt
on the reality of self-similarity in the networks. Therefore, we perform a
comprehensive evaluation of metabolic network fractality using a box-covering
method with an earlier version and the latest version of metabolic networks,
and demonstrate that the latest metabolic networks are almost self-dissimilar,
while the earlier ones are fractal, as reported in a number of previous
studies. This result may be because the networks were randomized because of an
increase in network density due to database updates, suggesting that the
previously observed network fractality was due to a lack of available data on
metabolic reactions. This finding may not entirely discount the importance of
self-similarity of metabolic networks. Rather, it highlights the need for a
more suitable definition of network fractality and a more careful examination
of self-similarity of metabolic networks.


  Background: With the ever increasing use of computational models in the
biosciences, the need to share models and reproduce the results of published
studies efficiently and easily is becoming more important. To this end, various
standards have been proposed that can be used to describe models, simulations,
data or other essential information in a consistent fashion. These constitute
various separate components required to reproduce a given published scientific
result.
  Results: We describe the Open Modeling EXchange format (OMEX). Together with
the use of other standard formats from the Computational Modeling in Biology
Network (COMBINE), OMEX is the basis of the COMBINE Archive, a single file that
supports the exchange of all the information necessary for a modeling and
simulation experiment in biology. An OMEX file is a ZIP container that includes
a manifest file, listing the content of the archive, an optional metadata file
adding information about the archive and its content, and the files describing
the model. The content of a COMBINE Archive consists of files encoded in
COMBINE standards whenever possible, but may include additional files defined
by an Internet Media Type. Several tools that support the COMBINE Archive are
available, either as independent libraries or embedded in modeling software.
  Conclusions: The COMBINE Archive facilitates the reproduction of modeling and
simulation experiments in biology by embedding all the relevant information in
one file. Having all the information stored and exchanged at once also helps in
building activity logs and audit trails. We anticipate that the COMBINE Archive
will become a significant help for modellers, as the domain moves to larger,
more complex experiments such as multi-scale models of organs, digital
organisms, and bioengineering.


  The PyMOL molecular graphics program has been modified to introduce a new
'putty' cartoon representation, akin to the 'sausage'-style representation of
the MOLMOL molecular visualization (MolVis) software package. This document
outlines the development and implementation of the putty representation.


  To facilitate flexible and efficient structural bioinformatics analyses, new
functionality for three-dimensional structure processing and analysis has been
introduced into PyCogent -- a popular feature-rich framework for sequence-based
bioinformatics, but one which has lacked equally powerful tools for handling
stuctural/coordinate-based data. Extensible Python modules have been developed,
which provide object-oriented abstractions (based on a hierarchical
representation of macromolecules), efficient data structures (e.g. kD-trees),
fast implementations of common algorithms (e.g. surface-area calculations),
read/write support for Protein Data Bank-related file formats and wrappers for
external command-line applications (e.g. Stride). Integration of this code into
PyCogent is symbiotic, allowing sequence-based work to benefit from
structure-derived data and, reciprocally, enabling structural studies to
leverage PyCogent's versatile tools for phylogenetic and evolutionary analyses.


  The structure of a network dramatically affects the spreading phenomena
unfolding upon it. The contact distribution of the nodes has long been
recognized as the key ingredient in influencing the outbreak events. However,
limited knowledge is currently available on the role of the weight of the edges
on the persistence of a pathogen. At the same time, recent works showed a
strong influence of temporal network dynamics on disease spreading. In this
work we provide an analytical understanding, corroborated by numerical
simulations, about the conditions for infected stable state in weighted
networks. In particular, we reveal the role of heterogeneity of edge weights
and of the dynamic assignment of weights on the ties in the network in driving
the spread of the epidemic. In this context we show that when weights are
dynamically assigned to ties in the network an heterogeneous distribution is
able to hamper the diffusion of the disease, contrary to what happens when
weights are fixed in time.


  Regulating the stability of microtubule(MT)-kinetochore attachments is
fundamental to avoiding mitotic errors and ensure proper chromosome segregation
during cell division. While biochemical factors involved in this process have
been identified, its mechanics still needs to be better understood. Here we
introduce and simulate a mechanical model of MT-kinetochore interactions in
which the stability of the attachment is ruled by the geometrical conformations
of curling MT-protofilaments entangled in kinetochore fibrils. The model allows
us to reproduce with good accuracy in vitro experimental measurements of the
detachment times of yeast kinetochores from MTs under external pulling forces.
Numerical simulations suggest that geometrical features of MT-protofilaments
may play an important role in the switch between stable and unstable
attachments.


  We present a probabilistic T-cell model that includes negative selection and
takes contrasting models of tissue-restricted antigen (TRA) expression in the
thymus into account. We start from the basic model of van den Berg, Rand, and
Burroughs (2001) and include negative selection via individual-based T-cell
modelling, in which each T-cell is defined by its stimulation rates to all
relevant self antigens. We present a simulation approach based on partial
tilting of the stimulation rates recognized by a single T-cell. We investigate
the effects of negative selection for diverging modes of thymic antigen
presentation, namely arbitrary TRA presentation, and more or less strict
emulation of tissue-specific cell lines. We observe that negative selection
leads to truncation of the tail of the distribution of the stimulation rates
mature T-cells receive from self antigens, i.e., the self background is
reduced. This increases the activation probabilities of single T-cells in the
presence of non-self antigens.


  We study the global spatio-temporal patterns of influenza dynamics. This is
achieved by analysing and modelling weekly laboratory confirmed cases of
influenza A and B from 138 countries between January 2006 and May 2014. The
data were obtained from FluNet, the surveillance network compiled by the the
World Health Organization. We report a pattern of {\it skip-and-resurgence}
behavior between the years 2011 and 2013 for influenza H1N1/09, the strain
responsible for the 2009 pandemic, in Europe and Eastern Asia. In particular,
the expected H1N1/09 epidemic outbreak in 2011 failed to occur (or"skipped") in
many countries across the globe, although an outbreak occurred in the following
year. We also report a pattern of {\it well-synchronized} 2010 winter wave of
H1N1/09 in the Northern Hemisphere countries, and a pattern of replacement of
strain H1N1/77 by H1N1/09 between the 2009 and 2012 influenza seasons. Using
both a statistical and a mechanistic mathematical model, and through fitting
the data of 108 countries (108 countries in a statistical model and 10 large
populations with a mechanistic model), we discuss the mechanisms that are
likely to generate these events taking into account the role of multi-strain
dynamics. A basic understanding of these patterns has important public health
implications and scientific significance.


  When a population is engaged in successive prisoner's dilemmas, indirect
reciprocity through reputation fosters cooperation through the emergence of
moral and action rules. A simplified model has recently been proposed where
individuals choose between helping or not others, and are judged good or bad
for it by the rest of the population. The reputation so acquired will condition
future actions. In this model, eight strategies (referred to as 'leading
eight') enforce a high level of cooperation, generate high payoffs and are
therefore resistant to invasions by other strategies. Here we show that, by
assigning each individual one out of two labels that peers can distinguish
(e.g., political ideas, religion, skin colour...) and allowing moral and action
rules to depend on the label, intolerant behaviours can emerge within
minorities under sufficient economic stress. We analyse the sets of conditions
where this can happen and also discuss the circumstances under which tolerance
can be restored. Our results agree with empirical observations that correlate
intolerance and economic stress, and predict a correlation between the degree
of tolerance of a population and its composition and ethical stance.


  We provide a short supplement to the paper "MAPK networks and their capacity
for multistationarity due to toric steady states" by P\'erez Mill\'an and
Turjanski. We show that the capacity for toric steady states in the three
networks analyzed in that paper can be derived using the process of network
translation, which corresponds the original mass action system to a generalized
mass action system with the same steady states. In all three cases, the
translated chemical reaction network is proper, weakly reversible, and has both
a structural and kinetic deficiency of zero. This is sufficient to guarantee
toric steady states by previously established work on network translations. A
basis of the steady state ideal is then derived by consideration of the linkage
classes of the translated chemical reaction network.


  Transcription factors perform facilitated diffusion (3D diffusion in the
cytosol and 1D diffusion on the DNA) when binding to their target sites to
regulate gene expression. Here, we investigated the influence of this binding
mechanism on the noise in gene expression. Our results showed that, for
biologically relevant parameters, the binding process can be represented by a
two-state Markov model and that the accelerated target finding due to
facilitated diffusion leads to a reduction in both the mRNA and the protein
noise.


  Dynamic modelling is one of the cornerstones of systems biology. Many
research efforts are currently being invested in the development and
exploitation of large-scale kinetic models. The associated problems of
parameter estimation (model calibration) and optimal experimental design are
particularly challenging. The community has already developed many methods and
software packages which aim to facilitate these tasks. However, there is a lack
of suitable benchmark problems which allow a fair and systematic evaluation and
comparison of these contributions. Here we present BioPreDyn-bench, a set of
challenging parameter estimation problems which aspire to serve as reference
test cases in this area. This set comprises six problems including medium and
large-scale kinetic models of the bacterium E. coli, baker's yeast S.
cerevisiae, the vinegar fly D. melanogaster, Chinese Hamster Ovary cells, and a
generic signal transduction network. The level of description includes
metabolism, transcription, signal transduction, and development. For each
problem we provide (i) a basic description and formulation, (ii)
implementations ready-to-run in several formats, (iii) computational results
obtained with specific solvers, (iv) a basic analysis and interpretation. This
suite of benchmark problems can be readily used to evaluate and compare
parameter estimation methods. Further, it can also be used to build test
problems for sensitivity and identifiability analysis, model reduction and
optimal experimental design methods. The suite, including codes and
documentation, can be freely downloaded from
http://www.iim.csic.es/%7egingproc/biopredynbench/.


  Recent studies of animal social networks have significantly increased our
understanding of animal behavior, social interactions, and many important
ecological and epidemiological processes. However, most of the studies are at
low temporal and spatial resolution due to the difficulty in recording accurate
contact information. Domestic animals such as cattle have social behavior and
serve as an excellent study system because their position can be explicitly and
continuously tracked, allowing their social networks to be accurately
constructed. We used radio-frequency tags to accurately track cattle position
and analyze high-resolution cattle social networks. We tested the hypothesis of
temporal stationarity and spatial homogeneity in these high-resolution networks
and demonstrated substantial spatial-temporal heterogeneity during different
daily time periods (feeding and non-feeding) and in different areas of the pen
(grain bunk, water trough, hay bunk, and other general pen area). The social
network structure is analyzed using global network characteristics (network
density, exponential random graph model structure), subgroup clustering
(modularity), triadic property (transitivity), and dyadic interactions
(correlation coefficient from a quadratic assignment procedure). Cattle tend to
have the strongest and most consistent contacts with others around the hay bunk
during the feeding time. These results cannot be determined from data at lower
spatial (aggregated at entire pen level) or temporal (aggregated at daily
level) resolution. These results reveal new insights for real-time animal
social network structure dynamics, providing more accurate descriptions that
allow more accurate modeling of multiple (both direct and indirect) disease
transmission pathways.


  Progress in cell type reprogramming has revived the interest in Waddington's
concept of the epigenetic landscape. Recently researchers developed the
quasi-potential theory to represent the Waddington's landscape. The
Quasi-potential U(x), derived from interactions in the gene regulatory network
(GRN) of a cell, quantifies the relative stability of network states, which
determine the effort required for state transitions in a multi-stable dynamical
system. However, quasi-potential landscapes, originally developed for
continuous systems, are not suitable for discrete-valued networks which are
important tools to study complex systems. In this paper, we provide a framework
to quantify the landscape for discrete Boolean networks (BNs). We apply our
framework to study pancreas cell differentiation where an ensemble of BN models
is considered based on the structure of a minimal GRN for pancreas development.
We impose biologically motivated structural constraints (corresponding to
specific type of Boolean functions) and dynamical constraints (corresponding to
stable attractor states) to limit the space of BN models for pancreas
development. In addition, we enforce a novel functional constraint
corresponding to the relative ordering of attractor states in BN models to
restrict the space of BN models to the biological relevant class. We find that
BNs with canalyzing/sign-compatible Boolean functions best capture the dynamics
of pancreas cell differentiation. This framework can also determine the genes'
influence on cell state transitions, and thus can facilitate the rational
design of cell reprogramming protocols.


  Prion is a misfolded protein found in mammals that causes infectious diseases
of the nervous system in humans and animals. Prion diseases are invariably
fatal and highly infectious neurodegenerative diseases that affect a wide
variety of mammalian species such as sheep and goats, cattle, deer, elk and
humans etc. Recent studies have shown that rabbits have a low susceptibility to
be infected by prion diseases with respect to other animals including humans.
The present study employs molecular dynamics (MD) means to unravel the
mechanism of rabbit prion proteins (RaPrPC) based on the recently available
rabbit NMR structures (of the wild-type and its two mutants of two surface
residues). The electrostatic charge distributions on the protein surface are
the focus when analysing the MD trajectories. It is found that we can conclude
that surface electrostatic charge distributions indeed contribute to the
structural stability of wild-type RaPrPC; this may be useful for the medicinal
treatment of prion diseases.


  In the study of networks, it is often insightful to use algorithms to
determine mesoscale features such as "community structure", in which densely
connected sets of nodes constitute "communities" that have sparse connections
to other communities. The most popular way of detecting communities
algorithmically is to optimize the quality function known as modularity. When
optimizing modularity, one compares the actual connections in a (static or
time-dependent) network to the connections obtained from a random-graph
ensemble that acts as a null model. The communities are then the sets of nodes
that are connected to each other densely relative to what is expected from the
null model. Clearly, the process of community detection depends fundamentally
on the choice of null model, so it is important to develop and analyze novel
null models that take into account appropriate features of the system under
study. In this paper, we investigate the effects of using null models that take
incorporate spatial information, and we propose a novel null model based on the
radiation model of population spread. We also develop novel synthetic spatial
benchmark networks in which the connections between entities are based on
distance or flux between nodes, and we compare the performance of both static
and time-dependent radiation null models to the standard ("Newman-Girvan") null
model for modularity optimization and a recently-proposed gravity null model.
In our comparisons, we use both the above synthetic benchmarks and
time-dependent correlation networks that we construct using countrywide dengue
fever incidence data for Peru. We also evaluate a recently-proposed correlation
null model, which was developed specifically for correlation networks that are
constructed from time series, on the epidemic-correlation data.


  Translation is an important stage in gene expression. During this stage,
macro-molecules called ribosomes travel along the mRNA strand linking
amino-acids together in a specific order to create a functioning protein. An
important question is how to maximize protein production. Indeed, translation
is known to consume most of the cell's energy and it is natural to assume that
evolution shaped this process so that it maximizes the protein production rate.
If this is indeed so then one can estimate various parameters of the
translation machinery by solving an appropriate mathematical optimization
problem. The same problem also arises in the context of synthetic biology,
namely, re-engineer heterologous genes in order to maximize their translation
rate in a host organism. We consider the problem of maximizing the protein
production rate using a computational model for translation-elongation called
the ribosome flow model (RFM). This model describes the flow of the ribosomes
along an mRNA chain of length n using a set of n first-order nonlinear ODEs. It
also includes n+1 positive parameters: the ribosomal initiation rate into the
mRNA chain, and n elongation rates along the chain sites. We show that the
steady-state translation rate in the RFM is a strictly concave function of its
parameters. This means that the problem of maximizing the translation rate
under a suitable constraint always admits a unique solution, and that this
solution can be determined using highly-efficient algorithms for solving convex
optimization problems even for large values of n. Furthermore, our analysis
shows that the optimal translation rate can be computed based only on the
optimal initiation rate and the elongation rate of the codons near the
beginning of the ORF. We discuss some applications of the theoretical results
to synthetic biology, molecular evolution, and functional genomics.


  Recently we proposed an algorithm for the fast reconstruction of compact
context-specific metabolic networks (FASTCORE) that allowed dropping the
reconstruction time to the time order of seconds (Vlassis et al.,2014). This
extremely low computational demand opens new possibilities for improving the
quality of the models. Several rounds of model reconstruction, testing of the
model's predictions against real experimental data, curation steps of the input
model and the set of core reactions as well as cross-validations assays are
required to reconstruct high-quality models. These semi-automated model
curations steps are in such extend not possible with competing algorithms due
to their high computational demands. To adapt FASTCORE for the integration of
microarray data, we therefore propose a new workflow: FASTCORMICS. FASTCORMICS
requires as input microarray data and a Genome-scale reconstruction.
FASTCORMICS is devoid of heuristic parameter settings and has a low
computational demand with overall building times in the order of a few minutes.
FASTCORMICS preprocesses the microarrays data with the discretization tool
Barcode (Zillox et al, 2007). Barcode uses prior knowledge on the intensity
distribution of each probe set for a given microarray platform to segregate
between expressed genes and non-expressed genes. This preprocessing step allows
circumventing the need of setting a heuristic expression threshold, which is
critical for the output models as in response to this threshold alternative
pathways or subsystems might be included or excluded, thereby heavily changing
the functionalities of the model.
  In general, FASTCORMICS outperforms competing algorithms and allows obtaining
high-quality, robust models in a high-throughput manner. This will allow the
use of metabolic modelling as routine process for the analysis of microarray
data e.g. in the field of personalized medicine.


  Gene regulatory network inference uses genome-wide transcriptome measurements
in response to genetic, environmental or dynamic perturbations to predict
causal regulatory influences between genes. We hypothesized that evolution also
acts as a suitable network perturbation and that integration of data from
multiple closely related species can lead to improved reconstruction of gene
regulatory networks. To test this hypothesis, we predicted networks from
temporal gene expression data for 3,610 genes measured during early embryonic
development in six Drosophila species and compared predicted networks to gold
standard networks of ChIP-chip and ChIP-seq interactions for developmental
transcription factors in five species. We found that (i) the performance of
single-species networks was independent of the species where the gold standard
was measured; (ii) differences between predicted networks reflected the known
phylogeny and differences in biology between the species; (iii) an integrative
consensus network which minimized the total number of edge gains and losses
with respect to all single-species networks performed better than any
individual network. Our results show that in an evolutionarily conserved
system, integration of data from comparable experiments in multiple species
improves the inference of gene regulatory networks. They provide a basis for
future studies on the numerous multi-species gene expression datasets for other
biological processes available in the literature.


  Biological databases can be analysed as a complex network which may reveal
some its underlying biological mechanisms. Frequently, such databases are
identified as scale-free networks or as hierarchical networks depending on
connectivity distributions or clustering coefficients. Since these databases do
grow over time, one would expect that their network topology may undergo some
changes. Here, we analysed the historical versions of miRBase, a database of
microRNAs where we performed an alignment of all mature and precursor miRNAs
and calculated a pairwise similarity index. We found that the clustering
coefficient shows important changes during the growth of this database. For two
consecutive versions of the year 2009 we found a strong modification of the
network topology which we were able to associate to a technological change in
miRNA discovery. To evaluate if these changes could have happened by chance, we
performed a set of simulations of the database growth by sampling the final
version of miRBase and creating several alternative histories of miRBase. None
of the simulations were close to the actual historical evolution of this
database, which we understand as a clear indication of a very strong
construction bias.


  Emerging coherent X-ray scattering patterns of single particles have shown
dominant morphological signatures in agreement with predictions of the
scattering model used for conventional protein crystallography. The key
question is if and to what extent these scattering patterns contain volumetric
information, and what model can retrieve it. The scattering model of protein
crystallography is valid for very small crystals or those like crystalized
biomolecules with small coherent subunits. But in the general case, it fails to
model the integrated intensities of diffraction spots, and cannot even find the
size of the crystal. The more rigorous and less employed alternative is a
purely-classical crystal-specific model, which bypasses the fundamental notion
of bulk and hence the non-classical X-ray scattering from bulk. This
contribution is Part 1 out of two reports, in which we seek to clarify the
assumptions of some different regimes and models of X-ray scattering and their
implications for single particle imaging. In this part, first basic concepts
and existing models are briefly reviewed. Then the predictions of the
conventional and the rigorous models for emerging scattering patterns of
protein nanocrystals (intermediate case between conventional crystals and
single particles) are contrasted, and the terminology conflict regarding
"Diffraction Theory" is addressed. With a clearer picture of crystal
scattering, Part 2 will focus on additional concepts, limitations, correction
schemes, and alternative models relevant to single particles. Aside from such
optical details, protein crystallography is an advanced tool of analytical
chemistry and not a self-contained optical imaging technique (despite
significant instrumental role of optical data). As such, its final results can
be neither confirmed nor rejected on mere optical grounds; i.e., no
jurisdiction for optics.


  Emerging coherent X-ray scattering patterns of single-particles have shown
dominant morphological signatures in agreement with predictions of the
scattering model used for conventional protein crystallography. The key
question is if and to what extent these scattering patterns contain volumetric
information, and what model can retrieve it. This contribution is Part 2 out of
two reports, in which we seek to clarify the assumptions of some different
regimes and models of X-ray scattering and their implications for single
particle imaging. In Part 1, basic concepts and existing scattering models
along with their implications for nanocrystals, and also the misconception of
using Diffraction Theory for volumetric scattering were addressed. Here in Part
2, specific challenges ahead of single particle imaging are addressed.
Limitations of the conventional scattering model in the test case of a sphere,
schemes for improving this common scattering model, ambiguities in meaning and
properties of "density map" and non-classical effects, the crucial role of
electromagnetic boundary conditions, the uniqueness of X-ray scattering inverse
problems, and additional vulnerabilities of phase retrieval and its relation
with "resolution" are discussed. We raise concerns about the unverified use of
the common scattering model of protein crystallography for arbitrary objects,
and also leaving fundamental questions such as meaning/uniqueness/properties of
the sought 3D profiles to phase retrieval algorithms.


  The basic reproduction number R0 -- the number of individuals directly
infected by an infectious person in an otherwise susceptible population -- is
arguably the most widely used estimator of how severe an epidemic outbreak can
be. This severity can be more directly measured as the fraction people infected
once the outbreak is over, {\Omega}. In traditional mathematical epidemiology
and common formulations of static network epidemiology, there is a
deterministic relationship between R0 and {\Omega}. However, if one considers
disease spreading on a temporal contact network -- where one knows when
contacts happen, not only between whom -- then larger R0 does not necessarily
imply larger {\Omega}. In this paper, we numerically investigate the
relationship between R0 and {\Omega} for a set of empirical temporal networks
of human contacts. Among 31 explanatory descriptors of temporal network
structure, we identify those that make R0 an imperfect predictor of {\Omega}.
We find that descriptors related to both temporal and topological aspects
affect the relationship between R0 and {\Omega}, but in different ways.


  We show that the cophylogenetic distance, k-interval cospeciation, is
distinct from other metrics and accounts for global congruence between locally
incongruent trees. The growth of the neighborhood of trees which satisfy the
largest possible k-interval cospeciation with a given tree indicates that
k-interval cospeciation is useful for analyzing simulated data.


  Limitation of resources has been recently introduced as a mechanism for the
survival and coexistence of cooperators with defectors in well-mixed
populations. Here we examine the same model on a scale-free network. A
prisoner's dilemma game on a scale-free network has shown coexistence of
cooperators and defectors for the entire range of parameters. Our results show
that by introducing the network to the limited resources model, the
cooperator-dominated region in the parameter space expands comparing to the
results of well-mixed population and the coexistence region becomes narrower.
The effect of scale-free network is therefore interpreted as improving the
cooperation in the population and reducing the coexistence.


  Movement within eukaryotic cells largely originates from localized forces
exerted by myosin motors on scaffolds of actin filaments. Although individual
motors locally exert both contractile and extensile forces, large actomyosin
structures at the cellular scale are overwhelmingly contractile, suggesting
that the scaffold serves to favor contraction over extension. While this
mechanism is well understood in highly organized striated muscle, its origin in
disordered networks such as the cell cortex is unknown. Here we develop a
mathematical model of the actin scaffold's local two- or three-dimensional
mechanics and identify four competing contraction mechanisms. We predict that
one mechanism dominates, whereby local deformations of the actin break the
balance between contraction and extension. In this mechanism, contractile
forces result mostly from motors plucking the filaments transversely rather
than buckling them longitudinally. These findings sheds light on recent
$\textit{in vitro}$ experiments, and provides a new geometrical understanding
of contractility in the myriad of disordered actomyosin systems found
$\textit{in vivo}$.


  Motivation: Methods are needed to test pre-defined genomic regions such as
promoters for differential methylation in genome-wide association studies,
where the number of samples is limited and the data have large amounts of
measurement error. Results: We developed a new statistical test, the
generalized integrated functional test (GIFT), which tests for regional
differences in methylation based on differences in the functional relationship
between methylation percent and location of the CpG sites within a region. In
this method, subject-specific functional profiles are first estimated, and the
average profile within groups is compared between groups using an ANOVA-like
test. Simulations and analyses of data obtained from patients with chronic
lymphocytic leukemia indicate that GIFT has good statistical properties and is
able to identify promising genomic regions. Further, GIFT is likely to work
with multiple different types of experiments since different smoothing
functions can be used to estimate the functional relationship between
methylation percent and CpG site location. Availability and Implementation:
Matlab code for GIFT and sample data are available at
http://biostat.gru.edu/~dryu/research.html. Contact: rpodolsk@med.wayne.edu or
dryu@gru.edu


  Many ontologies have been developed in biology and these ontologies
increasingly contain large volumes of formalized knowledge commonly expressed
in the Web Ontology Language (OWL). Computational access to the knowledge
contained within these ontologies relies on the use of automated reasoning. We
have developed the Aber-OWL infrastructure that provides reasoning services for
bio-ontologies. Aber-OWL consists of an ontology repository, a set of web
services and web interfaces that enable ontology-based semantic access to
biological data and literature. Aber-OWL is freely available at
http://aber-owl.net.


  Bird flocks are a paradigmatic example of collective motion. One of the
prominent experimental traits discovered about flocks is the presence of long
range velocity correlations between individuals, which allow them to influence
each other over the large scales, keeping a high level of group coordination. A
crucial question is to understand what is the mutual interaction between birds
generating such nontrivial correlations. Here we use the Maximum Entropy (ME)
approach to infer from experimental data of natural flocks the effective
interactions between birds. Compared to previous studies, we make a significant
step forward as we retrieve the full functional dependence of the interaction
on distance and find that it decays exponentially over a range of a few
individuals. The fact that ME gives a short-range interaction even though its
experimental input is the long-range correlation function, shows that the
method is able to discriminate the relevant information encoded in such
correlations and single out a minimal number of effective parameters. Finally,
we show how the method can be used to capture the degree of anisotropy of
mutual interactions.


  The repertoire of lymphocyte receptors in the adaptive immune system protects
organisms from diverse pathogens. A well-adapted repertoire should be tuned to
the pathogenic environment to reduce the cost of infections. We develop a
general framework for predicting the optimal repertoire that minimizes the cost
of infections contracted from a given distribution of pathogens. The theory
predicts that the immune system will have more receptors for rare antigens than
expected from the frequency of encounters; individuals exposed to the same
infections will have sparse repertoires that are largely different, but
nevertheless exploit cross-reactivity to provide the same coverage of antigens;
and the optimal repertoires can be reached via the dynamics of competitive
binding of antigens by receptors, and selective amplification of stimulated
receptors. Our results follow from a tension between the statistics of pathogen
detection, which favor a broader receptor distribution, and the effects of
cross-reactivity, which tend to concentrate the optimal repertoire onto a few
highly abundant clones. Our predictions can be tested in high throughput
surveys of receptor and pathogen diversity.


  The analysis of correlations of amino acid occurrences in globular proteins
has led to the development of statistical tools that can identify native
contacts -- portions of the chains that come to close distance in folded
structural ensembles. Here we introduce a statistical coupling analysis for
repeat proteins -- natural systems for which the identification of domains
remains challenging. We show that the inherent translational symmetry of repeat
protein sequences introduces a strong bias in the pair correlations at
precisely the length scale of the repeat-unit. Equalizing for this bias reveals
true co-evolutionary signals from which local native-contacts can be
identified. Importantly, parameter values obtained for all other interactions
are not significantly affected by the equalization. We quantify the robustness
of the procedure and assign confidence levels to the interactions, identifying
the minimum number of sequences needed to extract evolutionary information in
several repeat protein families. The overall procedure can be used to
reconstruct the interactions at long distances, identifying the characteristics
of the strongest couplings in each family, and can be applied to any system
that appears translationally symmetric.


  Recent technological advances in Next Generation Sequencing tools have led to
increasing speeds of DNA sample collection, preparation, and sequencing. One
instrument can produce over 600 Gb of genetic sequence data in a single run.
This creates new opportunities to efficiently handle the increasing workload.
We propose a new method of fast genetic sequence analysis using the Dynamic
Distributed Dimensional Data Model (D4M) - an associative array environment for
MATLAB developed at MIT Lincoln Laboratory. Based on mathematical and
statistical properties, the method leverages big data techniques and the
implementation of an Apache Acculumo database to accelerate computations
one-hundred fold over other methods. Comparisons of the D4M method with the
current gold-standard for sequence analysis, BLAST, show the two are comparable
in the alignments they find. This paper will present an overview of the D4M
genetic sequence algorithm and statistical comparisons with BLAST.


  A key aim of systems biology is the reconstruction of molecular networks,
however we do not yet have networks that integrate information from all
datasets available for a particular clinical condition. This is in part due to
the limited scalability, in terms of required computational time and power, of
existing algorithms. Network reconstruction methods should also be scalable in
the sense of allowing scientists from different backgrounds to efficiently
integrate additional data. We present a network model of acute myeloid leukemia
(AML). In the current version (AML 2.1) we have used gene expression data (both
microarray and RNA-seq) from five different studies comprising a total of 771
AML samples and a protein-protein interactions dataset. Our scalable network
reconstruction method is in part based on the well-known property of gene
expression correlation among interacting molecules. The difficulty of
distinguishing between direct and indirect interactions is addressed optimizing
the coefficient of variation of gene expression, using a validated gold
standard dataset of direct interactions. Computational time is much reduced
compared to other network reconstruction methods. A key feature is the study of
the reproducibility of interactions found in independent clinical datasets. An
analysis of the most significant clusters, and of the network properties
(intraset efficiency, degree, betweenness centrality and PageRank) of common
AML mutations demonstrated the biological significance of the network. A
statistical analysis of the response of blast cells from eleven AML patients to
a library of kinase inhibitors provided an experimental validation of the
network. A combination of network and experimental data identified CDK1, CDK2,
CDK4 and CDK6 and other kinases as potential therapeutic targets in AML.


  Given the importance of non-coding RNAs to cellular regulatory functions and
rapid growth of RNA transcripts, computational prediction of RNA tertiary
structure remains highly demanded yet significantly challenging. Even for a
short RNA sequence, the space of tertiary conformations is immense; existing
methods to identify native-like conformations mostly resort to random sampling
of conformations to gain computational feasibility. However native
conformations may not be examined and prediction accuracy may be compromised
due to sampling. In particular, the state-of-the-art methods have yet to
deliver the desired prediction performance for RNAs of length beyond 50.
  This paper presents the work to tackle a key step in the RNA tertiary
structure prediction problem, the prediction of the nucleotide interactions
that constitute the desired tertiary structure. The research is established
upon a novel graph model, called backbone k-tree, to markably constrain
nucleotide interaction relationships in RNA tertiary structure. It is shown
that the new model makes it possible to efficiently predict the optimal set of
nucleotide interactions from the query sequence, including the interactions in
all recently revealed families. Evident by the preliminary results, the new
method can predict with a high accuracy the nucleotide interactions that
constitute the tertiary structure of the query sequence, thus providing a
viable solution towards ab initio prediction of RNA tertiary structure.


  Turing mechanisms can yield a large variety of patterns from noisy,
homogenous initial conditions and have been proposed as patterning mechanism
for many developmental processes. However, the molecular components that give
rise to Turing patterns have remained elusive, and the small size of the
parameter space that permits Turing patterns to emerge makes it difficult to
explain how Turing patterns could evolve. We have recently shown that Turing
patterns can be obtained with a single ligand if the ligand-receptor
interaction is taken into account. Here we show that the general properties of
ligand-receptor systems result in very large Turing spaces. Thus, the
restriction of receptors to single cells, negative feedbacks, regulatory
interactions between different ligand-receptor systems, and the clustering of
receptors on the cell surface all greatly enlarge the Turing space. We further
show that the feedbacks that occur in the FGF10/SHH network that controls lung
branching morphogenesis are sufficient to result in large Turing spaces. We
conclude that the cellular restriction of receptors provides a mechanism to
sufficiently increase the size of the Turing space to make the evolution of
Turing patterns likely. Additional feedbacks may then have further enlarged the
Turing space. Given their robustness and flexibility, we propose that
receptor-ligand based Turing mechanisms present a general mechanism for
patterning in biology.


  After a brief review of the protein folding quantum theory and a short
discussion on its experimental evidences the mechanism of glucose transport
across membrane is studied from the point of quantum conformational transition.
The structural variations among four kinds of conformations of the human
glucose transporter GLUT1 (ligand free occluded, outward open, ligand bound
occluded and inward open) are looked as the quantum transition. The comparative
studies between mechanisms of uniporter (GLUT1) and symporter (XylE and GlcP)
are given. The transitional rates are calculated from the fundamental theory.
The monosaccharide transport kinetics is proposed. The steady state of the
transporter is found and its stability is studied. The glucose (xylose)
translocation rates in two directions and in different steps are compared. The
mean transport time in a cycle is calculated and based on it the comparison of
the transport times between GLUT1,GlcP and XylE can be drawn. The non-Arrhenius
temperature dependence of the transition rate and the mean transport time is
predicted. It is suggested that the direct measurement of temperature
dependence is a useful tool for deeply understanding the transmembrane
transport mechanism.


  Here, we address the essential question of whether, in the context of
evolving populations, ecosystems attain properties that enable persistence of
the ecosystem itself. We use a simple ecosystem model describing resource,
producer, and consumer dynamics to analyze how evolution affects dynamical
stability properties of the ecosystem. In particular, we compare resilience of
the entire system after allowing the producer and consumer populations to
evolve to their evolutionarily stable strategy (ESS), to the maximum attainable
resilience. We find a substantial reduction in ecosystem resilience when
producers and consumers are allowed to evolve compared to the maximal
attainable resilience. This study illustrates the inherent difference and
possible conflict between maximizing individual-level fitness and maximizing
resilience of entire ecosystems.


  The contribution of electrostatic interactions to the free energy of binding
between model protein and a ceramic implant surface in the aqueous solvent,
considered in the framework of the non-local electrostatic model, is calculated
as a function of the implant low-frequency dielectric constant. We show that
the existence of a dynamically ordered (low-dielectric) interfacial solvent
layer at the protein-solvent and ceramic-solvent interface markedly increases
charging energy of the protein and ceramic implant, and consequently makes the
electrostatic contribution to the protein-ceramic binding energy more favorable
(attractive). Our analysis shows that the corresponding electrostatic energy
between protein and oxide ceramics depends non-monotonically on the dielectric
constant of ceramic. Obtained results indicate that protein can attract
electrostatically to the surface if ceramic material has a moderate dielectric
constant below or about 35 (in particularly ZrO2 or Ta2O5). This is in contrast
to classical (local) consideration of the solvent, which demonstrates an
unfavorable electrostatic interaction of protein with typical metal oxide
ceramic materials (dielectric constant>10). Thus, a solid implant coated by
combining oxide ceramic with a reduced dielectric constant can be beneficial to
strengthen the electrostatic binding of the protein-implant complex.


  The effects of inbreeding on human health depend critically on the number and
severity of recessive, deleterious mutations carried by individuals. In humans,
existing estimates of these quantities are based on comparisons between
consanguineous and non-consanguineous couples, an approach that confounds
socioeconomic and genetic effects of inbreeding. To circumvent this limitation,
we focused on a founder population with almost complete Mendelian disease
ascertainment and a known pedigree. By considering all recessive lethal
diseases reported in the pedigree and simulating allele transmissions, we
estimated that each haploid set of human autosomes carries on average 0.29 (95%
credible interval [0.10, 0.83]) autosomal, recessive alleles that lead to
complete sterility or severe disorders at birth or before reproductive age when
homozygous. Comparison to existing estimates of the deleterious effects of all
recessive alleles suggests that a substantial fraction of the burden of
autosomal, recessive variants is due to single mutations that lead to death
between birth and reproductive age. In turn, the comparison to estimates from
other eukaryotes points to a surprising constancy of the average number of
recessive lethal mutations across organisms with markedly different genome
sizes.


  This paper builds upon the fundamental work of Niwa et al. [34], which
provides the unique possibility to analyze the relative aggregation/folding
propensity of the elements of the entire Escherichia coli (E. coli) proteome in
a cell-free standardized microenvironment. The hardness of the problem comes
from the superposition between the driving forces of intra- and inter-molecule
interactions and it is mirrored by the evidences of shift from folding to
aggregation phenotypes by single-point mutations [10]. Here we apply several
state-of-the-art classification methods coming from the field of structural
pattern recognition, with the aim to compare different representations of the
same proteins gathered from the Niwa et al. data base; such representations
include sequences and labeled (contact) graphs enriched with chemico-physical
attributes. By this comparison, we are able to identify also some interesting
general properties of proteins. Notably, (i) we suggest a threshold around 250
residues discriminating "easily foldable" from "hardly foldable" molecules
consistent with other independent experiments, and (ii) we highlight the
relevance of contact graph spectra for folding behavior discrimination and
characterization of the E. coli solubility data. The soundness of the
experimental results presented in this paper is proved by the statistically
relevant relationships discovered among the chemico-physical description of
proteins and the developed cost matrix of substitution used in the various
discrimination systems.


  We investigated testing the likelihood of a phylogenetic tree by comparison
to its subtree pruning and regrafting (SPR) neighbors, with or without
re-optimizing branch lengths. This is inspired by aspects of Bayesian
significance tests, and the use of SPRs for heuristically finding maximum
likelihood trees. Through a number of simulations with the Jukes-Cantor model
on various topologies, it is observed that the SPR tests are informative, and
reasonably fast compared to searching for the maximum likelihood tree. This
suggests that the SPR tests would be a useful addition to the suite of existing
statistical tests, for identifying potential inaccuracies of inferred
topologies.


  Understanding how the time-complexity of evolutionary algorithms (EAs) depend
on their parameter settings and characteristics of fitness landscapes is a
fundamental problem in evolutionary computation. Most rigorous results were
derived using a handful of key analytic techniques, including drift analysis.
However, since few of these techniques apply effortlessly to population-based
EAs, most time-complexity results concern simplified EAs, such as the (1+1) EA.
  This paper describes the level-based theorem, a new technique tailored to
population-based processes. It applies to any non-elitist process where
offspring are sampled independently from a distribution depending only on the
current population. Given conditions on this distribution, our technique
provides upper bounds on the expected time until the process reaches a target
state.
  We demonstrate the technique on several pseudo-Boolean functions, the sorting
problem, and approximation of optimal solutions in combinatorial optimisation.
The conditions of the theorem are often straightforward to verify, even for
Genetic Algorithms and Estimation of Distribution Algorithms which were
considered highly non-trivial to analyse. Finally, we prove that the theorem is
nearly optimal for the processes considered. Given the information the theorem
requires about the process, a much tighter bound cannot be proved.


  We study a universal object for the genealogy of a sample in populations with
mutations: the critical birth-death process with Poissonian mutations,
conditioned on its population size at a fixed time horizon. We show how this
process arises as the law of the genealogy of a sample in a large class of
critical branching populations with mutations at birth, namely populations
converging, in a large population asymptotic, towards the continuum random
tree. We extend this model to populations with random foundation times, with
(potentially improper) prior distributions g_i: x\mapsto x^{-i}, i\in\Z_+,
including the so-called uniform (i=0) and log-uniform (i=1) priors.
  We first investigate the mutational patterns arising from these models, by
studying the site frequency spectrum of a sample with fixed size, i.e. the
number of mutations carried by k individuals in the sample. Explicit formulae
for the expected frequency spectrum of a sample are provided, in the cases of a
fixed foundation time, and of a uniform and log-uniform prior on the foundation
time. Second, we establish the convergence in distribution, for large sample
sizes, of the (suitably renormalized) tree spanned by the sample genealogy with
prior g_i on the time of origin. We finally prove that the limiting genealogies
with different priors can all be embedded in the same realization of a given
Poisson point measure.


  A Yule tree is the result of a branching process with constant birth and
death rates. Such a process serves as an instructive null model of many
empirical systems, for instance, the evolution of species leading to a
phylogenetic tree. However, often in phylogeny the only available information
is the pairwise distances between a small fraction of extant species
representing the leaves of the tree. In this article we study statistical
properties of the pairwise distances in a Yule tree. Using a method based on a
recursion, we derive an exact, analytic and compact formula for the expected
number of pairs separated by a certain time distance. This number turns out to
follow a increasing exponential function. This property of a Yule tree can
serve as a simple test for empirical data to be well described by a Yule
process. We further use this recursive method to calculate the expected number
of the $n$-most closely related pairs of leaves and the number of cherries
separated by a certain time distance. To make our results more useful for
realistic scenarios, we explicitly take into account that the leaves of a tree
may be incompletely sampled and derive a criterion for poorly sampled
phylogenies. We show that our result can account for empirical data, using two
families of birds species.


  The objective of this study was to characterize the spatiotemporal dynamics
and climatic and social-ecological risk factors associated with the largest
dengue epidemic to date in Machala, Ecuador, to inform the development of a
dengue EWS. The following data were included in analyses: neighborhood-level
georeferenced dengue cases, national census data, and entomological
surveillance data from 2010; time series of weekly dengue cases (aggregated to
the city-level) and meteorological data from 2003 to 2012. We applied LISA and
Morans I to analyze the spatial distribution of the 2010 dengue cases, and
developed multivariate logistic regression models through a multi-model
selection process to identify census variables and entomological covariates
associated with the presence of dengue at the neighborhood level. Using data
aggregated at the city-level, we conducted a time-series (wavelet) analysis of
weekly climate and dengue incidence (2003-2012) to identify significant time
periods (e.g., annual, biannual) when climate co-varied with dengue, and to
describe the climate conditions associated with the 2010 outbreak. We found
significant hotspots of dengue transmission near the center of Machala. The
best-fit model to predict the presence of dengue included older age and female
gender of the head of the household, greater access to piped water in the home,
poor housing condition, and less distance to the central hospital. Wavelet
analyses revealed that dengue transmission co-varied with rainfall and minimum
temperature at annual and biannual cycles, and we found that anomalously high
rainfall and temperatures were associated with the 2010 outbreak. Our findings
highlight the importance of geospatial information in dengue surveillance and
the potential to develop a climate-driven spatiotemporal prediction models to
inform disease prevention and control interventions.


  Low-coverage short-read resequencing experiments have the potential to expand
our understanding of Y chromosome haplogroups. However, the uncertainty
associated with these experiments mean that haplogroups must be assigned
probabilistically to avoid false inferences. We propose an efficient dynamic
programming algorithm that can assign haplogroups by maximum likelihood, and
represent the uncertainty in assignment. We apply this to both genotype and
low-coverage sequencing data, and show that it can assign haplogroups
accurately and with high resolution. The method is implemented as the program
YFitter, which can be downloaded from http://sourceforge.net/projects/yfitter/


  RNA virus populations will undergo processes of mutation and selection
resulting in a mixed population of viral particles. High throughput sequencing
of a viral population subsequently contains a mixed signal of the underlying
clones. We would like to identify the underlying evolutionary structures. We
utilize two sources of information to attempt this; within segment linkage
information, and mutation prevalence. We demonstrate that clone haplotypes,
their prevalence, and maximum parsimony reticulate evolutionary structures can
be identified, although the solutions may not be unique, even for complete sets
of information. This is applied to a chain of influenza infection, where we
infer evolutionary structures, including reassortment, and demonstrate some of
the difficulties of interpretation that arise from deep sequencing due to
artifacts such as template switching during PCR amplification.


  Ecological spatial patterns reflect the underlying processes that shape the
structure of species and communities. Mechanisms like inter and intra species
competition, dispersal and host-pathogen interactions are believed to act over
a wide range of scales, and the inference of the process from the pattern is,
despite its popularity, a challenging task. Here we call attention to a quite
unexpected phenomenon in the extensively studied tropical forest at the
Barro-Colorado Island (BCI): the spatial deployment of (almost) all tree
species is statistically equivalent, once distances are normalized by , the
typical distance between neighboring conspecific trees. Correlation function,
cluster statistics and nearest-neighbor distance distribution become
species-independent after this rescaling. Global observables (species
frequencies) and local spatial structure appear to be interrelated. This
"glocality" suggests a radical interpretation of recent experiments that show a
correlation between species' abundance and the negative feedback among
conspecifics. For the forest to be glocal, the negative feedback must govern
spatial patterns over all scales.


  This paper deals with the relations among structural, topological, and
chemical properties of the E.Coli proteome from the vantage point of the
solubility/aggregation propensity of proteins. Each E.Coli protein is initially
represented according to its known folded 3D shape. This step consists in
representing the available E.Coli proteins in terms of graphs. We first analyze
those graphs by considering pure topological characterizations, i.e., by
analyzing the mass fractal dimension and the distribution underlying both
shortest paths and vertex degrees. Results confirm the general architectural
principles of proteins. Successively, we focus on the statistical properties of
a representation of such graphs in terms of vectors composed of several
numerical features, which we extracted from their structural representation. We
found that protein size is the main discriminator for the solubility, while
however there are other factors that help explaining the solubility degree. We
finally analyze such data through a novel one-class classifier, with the aim of
discriminating among very and poorly soluble proteins. Results are encouraging
and consolidate the potential of pattern recognition techniques when employed
to describe complex biological systems.


  Markov state models (MSMs) are a widely used method for approximating the
eigenspectrum of the molecular dynamics propagator, yielding insight into the
long-timescale statistical kinetics and slow dynamical modes of biomolecular
systems. However, the lack of a unified theoretical framework for choosing
between alternative models has hampered progress, especially for non-experts
applying these methods to novel biological systems. Here, we consider
cross-validation with a new objective function for estimators of these slow
dynamical modes, a generalized matrix Rayleigh quotient (GMRQ), which measures
the ability of a rank-$m$ projection operator to capture the slow subspace of
the system. It is shown that a variational theorem bounds the GMRQ from above
by the sum of the first $m$ eigenvalues of the system's propagator, but that
this bound can be violated when the requisite matrix elements are estimated
subject to statistical uncertainty. This overfitting can be detected and
avoided through cross-validation. These result make it possible to construct
Markov state models for protein dynamics in a way that appropriately captures
the tradeoff between systematic and statistical errors.


  Anomalous diffusion, in particular subdiffusion, is frequently invoked as a
mechanism of motion in dense biological media, and may have a significant
impact on the kinetics of binding/unbinding events at the cellular level. In
this work we incorporate anomalous diffusion in a previously developed model
for FRAP experiments. Our particular implementation of subdiffusive transport
is based on a continuous time random walk (CTRW) description of the motion of
fluorescent particles, as CTRWs lend themselves particularly well to the
inclusion of binding/unbinding events. In order to model switching between
bound and unbound states of fluorescent subdiffusive particles, we derive a
fractional reaction-subdiffusion equation of rather general applicability.
Using suitable initial and boundary conditions, this equation is then
incorporated in the model describing two-dimensional kinetics of FRAP
experiments. We find that this model can be used to obtain excellent fits to
experimental data. Moreover, recovery curves corresponding to different radii
of the circular bleach spot can be fitted by a single set of parameters. While
not enough evidence has been collected to claim with certainty that CTRW is the
underlying transport mechanism in FRAP experiments, the compatibility of our
results with experimental data fuels the discussion as to whether normal
diffusion or anomalous diffusion is the appropriate model, and as to whether
anomalous diffusion effects are important to fully understand the outcomes of
FRAP experiments. On a more technical side, we derive explicit analytic
solutions of our model in certain limits.


  A simple mean field model is proposed that captures the two scales of dynamic
processes involved in the spread of the coffee rust disease. At a local level
the state variable is the fraction of coffee plants on an average farm that are
infested with the disease and at the regional level the state variable is the
fraction of farms that are infested with the disease. These two levels are
interactive with one another in obvious ways, producing qualitatively
interesting results. Specifically there are three distinct outcomes: disease
disappears, disease is potentially persistent, and disease is inevitably
persistent. The basic structure of the model generates the potential for
tipping point behavior due to the presence of a blue sky bifurcation,
suggesting that the emergence of epizootics of the disease may be
unpredictable.


  Since most analysis software for genome-wide association studies (GWAS)
currently exploit only unrelated individuals, there is a need for efficient
applications that can handle general pedigree data or mixtures of both
population and pedigree data. Even data sets thought to consist of only
unrelated individuals may include cryptic relationships that can lead to false
positives if not discovered and controlled for. In addition, family designs
possess compelling advantages. They are better equipped to detect rare
variants, control for population stratification, and facilitate the study of
parent-of-origin effects. Pedigrees selected for extreme trait values often
segregate a single gene with strong effect. Finally, many pedigrees are
available as an important legacy from the era of linkage analysis.
Unfortunately, pedigree likelihoods are notoriously hard to compute. In this
paper we re-examine the computational bottlenecks and implement ultra-fast
pedigree-based GWAS analysis. Kinship coefficients can either be based on
explicitly provided pedigrees or automatically estimated from dense markers.
Our strategy (a) works for random sample data, pedigree data, or a mix of both;
(b) entails no loss of power; (c) allows for any number of covariate
adjustments, including correction for population stratification; (d) allows for
testing SNPs under additive, dominant, and recessive models; and (e)
accommodates both univariate and multivariate quantitative traits. On a typical
personal computer (6 CPU cores at 2.67 GHz), analyzing a univariate HDL
(high-density lipoprotein) trait from the San Antonio Family Heart Study
(935,392 SNPs on 1357 individuals in 124 pedigrees) takes less than 2 minutes
and 1.5 GB of memory. Complete multivariate QTL analysis of the three
time-points of the longitudinal HDL multivariate trait takes less than 5
minutes and 1.5 GB of memory.


  The noisy dynamics of chemical systems is commonly studied using either the
chemical master equation (CME) or the chemical Fokker-Planck equation (CFPE).
The latter is a continuum approximation of the discrete CME approach. We here
show that the CFPE may fail to capture the CME's prediction of noise-induced
multistability. In particular we find a simple chemical system for which the
CME's marginal probability distribution changes from unimodal to multimodal as
the system-size decreases below a critical value, while the CFPE's marginal
probability distribution is unimodal for all physically meaningful system
sizes.


  Pedigree GWAS (Option 29) in the current version of the Mendel software is an
optimized subroutine for performing large scale genome-wide QTL analysis. This
analysis (a) works for random sample data, pedigree data, or a mix of both, (b)
is highly efficient in both run time and memory requirement, (c) accommodates
both univariate and multivariate traits, (d) works for autosomal and x-linked
loci, (e) correctly deals with missing data in traits, covariates, and
genotypes, (f) allows for covariate adjustment and constraints among
parameters, (g) uses either theoretical or SNP-based empirical kinship matrix
for additive polygenic effects, (h) allows extra variance components such as
dominant polygenic effects and household effects, (i) detects and reports
outlier individuals and pedigrees, and (j) allows for robust estimation via the
$t$-distribution. The current paper assesses these capabilities on the genetics
analysis workshop 19 (GAW19) sequencing data. We analyzed simulated and real
phenotypes for both family and random sample data sets. For instance, when
jointly testing the 8 longitudinally measured systolic blood pressure (SBP) and
diastolic blood pressure (DBP) traits, it takes Mendel 78 minutes on a standard
laptop computer to read, quality check, and analyze a data set with 849
individuals and 8.3 million SNPs. Genome-wide eQTL analysis of 20,643
expression traits on 641 individuals with 8.3 million SNPs takes 30 hours using
20 parallel runs on a cluster. Mendel is freely available at
\url{http://www.genetics.ucla.edu/software}.