  Several biological tissues undergo changes in their geometry and in their
bulk material properties by modelling and remodelling processes. Modelling
synthesises tissue in some regions and removes tissue in others. Remodelling
overwrites old tissue material properties with newly formed, immature tissue
properties. As a result, tissues are made up of different "patches", i.e.,
adjacent tissue regions of different ages and different material properties,
within evolving boundaries. In this paper, generalised equations governing the
spatio-temporal evolution of such tissues are developed within the continuum
model. These equations take into account nonconservative, discontinuous surface
mass balance due to creation and destruction of material at moving interfaces,
and bulk balance due to tissue maturation. These equations make it possible to
model patchy tissue states and their evolution without explicitly maintaining a
record of when/where resorption and formation processes occurred. The time
evolution of spatially averaged tissue properties is derived systematically by
integration. These spatially-averaged equations cannot be written in closed
form as they retain traces that tissue destruction is localised at tissue
boundaries.
  The formalism developed in this paper is applied to bone tissues, which
exhibit strong material heterogeneities due to their slow mineralisation and
remodelling processes. Evolution equations are proposed in particular for
osteocyte density and bone mineral density. Effective average equations for
bone mineral density (BMD) and tissue mineral density (TMD) are derived using a
mean-field approximation. The error made by this approximation when remodelling
patchy tissue is investigated. The specific time signatures of BMD or TMD
during remodelling events may provide a way to detect these events occurring at
lower, unseen spatial resolutions from microCT scans.


  Crystallization is a key step in macromolecular structure determination by
crystallography. While a robust theoretical treatment of the process is
available, due to the complexity of the system, the experimental process is
still largely one of trial and error. In this article, efforts in the field are
discussed together with a theoretical underpinning using a solubility phase
diagram. Prior knowledge has been used to develop tools that computationally
predict the crystallization outcome and define mutational approaches that
enhance the likelihood of crystallization. For the most part these tools are
based on binary outcomes (crystal or no crystal), and the full information
contained in an assembly of crystallization screening experiments is lost. The
potential of this additional information is illustrated by examples where new
biological knowledge can be obtained and where a target can be sub-categorized
to predict which class of reagents provides the crystallization driving force.
Computational analysis of crystallization requires complete and correctly
formatted data. While massive crystallization screening efforts are under way,
the data available from many of these studies are sparse. The potential for
this data and the steps needed to realize this potential are discussed.


  How exponentially growing cells maintain size homeostasis is an important
fundamental problem. Recent single-cell studies in prokaryotes have uncovered
the adder principle, where cells on average, add a fixed size (volume) from
birth to division. Interestingly, this added volume differs considerably among
genetically-identical newborn cells with similar sizes suggesting a stochastic
component in the timing of cell-division. To mechanistically explain the adder
principle, we consider a time-keeper protein that begins to get stochastically
expressed after cell birth at a rate proportional to the volume. Cell-division
time is formulated as the first-passage time for protein copy numbers to hit a
fixed threshold. Consistent with data, the model predicts that while the mean
cell-division time decreases with increasing size of newborns, the noise in
timing increases with size at birth. Intriguingly, our results show that the
distribution of the volume added between successive cell-division events is
independent of the newborn cell size. This was dramatically seen in
experimental studies, where histograms of the added volume corresponding to
different newborn sizes collapsed on top of each other. The model provides
further insights consistent with experimental observations: the distributions
of the added volume and the cell-division time when scaled by their respective
means become invariant of the growth rate. Finally, we discuss various
modifications to the proposed model that lead to deviations from the adder
principle. In summary, our simple yet elegant model explains key experimental
findings and suggests a mechanism for regulating both the mean and fluctuations
in cell-division timing for size control.


  Understanding the processes that give rise to quantitative measurements
associated with molecular sequence data remains an important issue in
statistical phylogenetics. Examples of such measurements include geographic
coordinates in the context of phylogeography and phenotypic traits in the
context of comparative studies. A popular approach is to model the evolution of
continuously varying traits as a Brownian diffusion process. However, standard
Brownian diffusion is quite restrictive and may not accurately characterize
certain trait evolutionary processes. Here, we relax one of the major
restrictions of standard Brownian diffusion by incorporating a nontrivial
estimable drift into the process. We introduce a relaxed drift diffusion model
for the evolution of multivariate continuously varying traits along a
phylogenetic tree via Brownian diffusion with drift. Notably, the relaxed drift
model accommodates branch-specific variation of drift rates while preserving
model identifiability. We implement the relaxed drift model in a Bayesian
inference framework to simultaneously reconstruct the evolutionary histories of
molecular sequence data and associated multivariate continuous trait data, and
provide tools to visualize evolutionary reconstructions. We illustrate our
approach in three viral examples. In the first two, we examine the
spatiotemporal spread of HIV-1 in central Africa and West Nile virus in North
America and show that a relaxed drift approach uncovers a clearer, more
detailed picture of the dynamics of viral dispersal than standard Brownian
diffusion. Finally, we study antigenic evolution in the context of HIV-1
resistance to three broadly neutralizing antibodies. Our analysis reveals
evidence of a continuous drift at the HIV-1 population level towards enhanced
resistance to neutralization by the VRC01 monoclonal antibody over the course
of the epidemic.


  Wide-angle x-ray scattering (WAXS) is emerging as a powerful tool for
increasing the resolution of solution structure measurements of biomolecules.
Compared to its better known complement, small angle x-ray scattering (SAXS),
WAXS targets higher scattering angles and can enhance structural studies of
molecules by accessing finer details of solution structures. Although the
extension from SAXS to WAXS is easy to implement experimentally, the
computational tools required to fully harness the power of WAXS are still under
development. Currently, WAXS is employed to study structural changes and ligand
binding in proteins; however the methods are not as fully developed for nucleic
acids. Here, we show how WAXS can qualitatively characterize nucleic acid
structures as well as the small but significant structural changes driven by
the addition of multivalent ions. We show the potential of WAXS to test
all-atom molecular dynamics (MD) simulations and to provide insight in
understanding how the trivalent ion cobalt(III) hexammine (CoHex) affects the
structure of RNA and DNA helices. We find that MD simulations capture the RNA
structural change that occurs due to addition of CoHex.


  As cells grow and divide under a given environment, they become crowded and
resources are limited, as seen in bacterial biofilms and multicellular
aggregates. These cells often show strong interactions through exchanging
chemicals, as in quorum sensing, to achieve mutualism. Here, to achieve stable
division of labor, three properties are required. First, isogenous cells
differentiate into several types. Second, this aggregate of distinct cell types
shows better growth than that of isolated cells, by achieving division of
labor. Third, this cell aggregate is robust in the number distribution of
differentiated cell types. We here address how cells acquire the ability of
cell differentiation and division of labor simultaneously, which is also
connected with the robustness of a cell society. For this purpose, we developed
a dynamical-systems model of cells consisting of chemical components with
intracellular catalytic reaction dynamics. The reactions convert external
nutrients into internal components for cellular growth, and the divided cells
interact via chemical diffusion. We found that cells sharing an identical
catalytic network spontaneously differentiate via induction from cell-cell
interactions, and then achieve division of labor, enabling a higher growth rate
than that in the unicellular case. This symbiotic differentiation emerged for a
class of reaction networks with limited resources and strong cell-cell
interactions. Then, robustness in the cell type distribution was achieved,
while instability of collective growth could emerge even among the cooperative
cells when the internal reserves of products were dominant. The present
mechanism is simple and general as a natural result of interacting cells with
resource limitation, and is consistent with the observed behaviors and forms of
several aggregates of unicellular organisms.


  Proteins are essential components of living systems, capable of performing a
huge variety of tasks at the molecular level, such as recognition, signalling,
copy, transport, ... The protein sequences realizing a given function may
largely vary across organisms, giving rise to a protein family. Here, we
estimate the entropy of those families based on different approaches, including
Hidden Markov Models used for protein databases and inferred statistical models
reproducing the low-order (1-and 2-point) statistics of multi-sequence
alignments. We also compute the entropic cost, that is, the loss in entropy
resulting from a constraint acting on the protein, such as the fixation of one
particular amino-acid on a specific site, and relate this notion to the escape
probability of the HIV virus. The case of lattice proteins, for which the
entropy can be computed exactly, allows us to provide another illustration of
the concept of cost, due to the competition of different folds. The relevance
of the entropy in relation to directed evolution experiments is stressed.


  Multivariate regression is a widespread computational technique that may give
meaningless results if the explanatory variables are too numerous or highly
collinear. Tikhonov regularization, or ridge regression, is a popular approach
to address this issue. We reveal here a formal analogy between ridge regression
and statistical mechanics, where the objective function is comparable to a free
energy and the ridge parameter plays the role of temperature. This analogy
suggests two new criteria to select a suitable ridge parameter: the
specific-heat (Cv) and the maximum penalty (MP) fits. We apply these methods to
the calibration of the force constant in elastic network models (ENM). This key
parameter determines the amplitude of the predicted atomic fluctuations, and is
commonly obtained by fitting crystallographic B-factors. However, rigid-body
motions are often partially neglected in such fits, even though their
importance has been repeatedly stressed. Considering the full set of rigid-body
and internal degrees of freedom bears significant risks of overfitting, due to
strong correlations between explanatory variables, and requires thus careful
regularization. Using simulated data, we show that ridge regression with the Cv
or MP criterion markedly reduces the error of the estimated force constant, its
across-protein variation, and the number of proteins with unphysical values of
the fit parameters, in comparison with popular regularization schemes such as
generalized cross-validation. When applied to protein crystals, the new methods
provide a more robust calibration of ENM force constants, even though
rigid-body motions account on average for more than 80% of the amplitude of
B-factors. While MP emerges as the optimal choice for fitting crystallographic
B-factors, the Cv fit is more robust to the nature of the data, and is thus an
interesting candidate for other applications.


  Gradient sensing underlies important biological processes including
morphogenesis, polarization, and cell migration. The precision of gradient
sensing increases with the length of a detector (a cell or group of cells) in
the gradient direction, since a longer detector spans a larger range of
concentration values. Intuition from analyses of concentration sensing suggests
that precision should also increase with detector length in the direction
transverse to the gradient, since then spatial averaging should reduce the
noise. However, here we show that, unlike for concentration sensing, the
precision of gradient sensing decreases with transverse length for the simplest
gradient sensing model, local excitation--global inhibition (LEGI). The reason
is that gradient sensing ultimately relies on a subtraction of measured
concentration values. While spatial averaging indeed reduces the noise in these
measurements, which increases precision, it also reduces the covariance between
the measurements, which results in the net decrease in precision. We
demonstrate how a recently introduced gradient sensing mechanism, regional
excitation--global inhibition (REGI), overcomes this effect and recovers the
benefit of transverse averaging. Using a REGI-based model, we compute the
optimal two- and three-dimensional detector shapes, and argue that they are
consistent with the shapes of naturally occurring gradient-sensing cell
populations.


  The structure and function of biological molecules are strongly influenced by
the water and dissolved ions that surround them. This aqueous solution
(solvent) exerts significant electrostatic forces in response to the
biomolecule's ubiquitous atomic charges and polar chemical groups. In this
work, we investigate a simple approach to numerical calculation of this model
using boundary-integral equation (BIE) methods and boundary-element methods
(BEM). Traditional BEM discretizes the protein--solvent boundary into a set of
boundary elements, or panels, and the approximate solution is defined as a
weighted combination of basis functions with compact support. The resulting BEM
matrix then requires integrating singular or near singular functions, which can
be slow and challenging to compute. Here we investigate the accuracy and
convergence of a simpler representation, namely modeling the unknown surface
charge distribution as a set of discrete point charges on the surface. We find
that at low resolution, point-based BEM is more accurate than panel-based
methods, due to the fact that the protein surface is sampled directly, and can
be of significant value for numerous important calculations that require only
moderate accuracy, such as the preliminary stages of rational drug design and
protein engineering.


  Electrostatic forces play many important roles in molecular biology, but are
hard to model due to the complicated interactions between biomolecules and the
surrounding solvent, a fluid composed of water and dissolved ions. Continuum
model have been surprisingly successful for simple biological questions, but
fail for important problems such as understanding the effects of protein
mutations. In this paper we highlight the advantages of boundary-integral
methods for these problems, and our use of boundary integrals to design and
test more accurate theories. Examples include a multiscale model based on
nonlocal continuum theory, and a nonlinear boundary condition that captures
atomic-scale effects at biomolecular surfaces.


  Minimizing time delays in manufacturing vaccines appropriate to rapidly
mutating viruses is the key step for improving vaccine effectiveness. The
vaccine for the H3N2 flu type has failed for the last two years (~ 15%
effective). Here we summarize the state of the predictive art and report the
most current results for H3N2 flu vaccine design. Using a 2006 model of
dimensional reduction of viral mutational complexity, we show that this model
can reduce vaccine time delays by a year or more in some cases.


  Although known that single domain proteins fold and unfold by parallel
pathways, demonstration of this expectation has been difficult to establish in
experiments. Unfolding rate, $k_\mathrm{u}(f)$, as a function of force $f$,
obtained in single molecule pulling experiments on src SH3 domain, exhibits
upward curvature on a $\log k_\mathrm{u}(f)$ plot. Similar observations were
reported for other proteins for the unfolding rate $k_\mathrm{u}([C])$. These
findings imply unfolding in these single domain proteins involves a switch in
the pathway as $f$ or $[C]$ is increased from a low to a high value. We provide
a unified theory demonstrating that if $\log k_\mathrm{u}$ as a function of a
perturbation ($f$ or $[C]$) exhibits upward curvature then the underlying
energy landscape must be strongly multidimensional. Using molecular simulations
we provide a structural basis for the switch in the pathways and dramatic
shifts in the transition state ensemble (TSE) in src SH3 domain as $f$ is
increased. We show that a single point mutation shifts the upward curvature in
$\log k_\mathrm{u}(f)$ to a lower force, thus establishing the malleability of
the underlying folding landscape. Our theory, applicable to any perturbation
that affects the free energy of the protein linearly, readily explains movement
in the TSE in a $\beta$-sandwich (I27) protein and single chain monellin as the
denaturant concentration is varied. We predict that in the force range
accessible in laser optical tweezer experiments there should be a switch in the
unfolding pathways in I27 or its mutants.


  A spatially explicit model previously developed to apply to a population of
arboreal nesting ants in southern Mexico is explored in detail, including
notions of criticality and robust criticality and expected power law
adjustment.


  The diverse, specialized genes in today's lifeforms evolved from a common
core of ancient, elementary genes. However, these genes did not evolve
individually: gene expression is controlled by a complex network of
interactions, and alterations in one gene may drive reciprocal changes in its
proteins' binding partners. We show that the topology of a leukemia gene
regulatory network is strongly coupled with evolutionary properties.
Slowly-evolving ("cold"), old genes tend to interact with each other, as do
rapidly-evolving ("hot"), young genes, causing genes to evolve in clusters. We
argue that gene duplication placed old, cold genes at the center of the
network, and young, hot genes on the periphery, and demonstrate this with
single-node centrality measures and two new measures of efficiency. Integrating
centrality measures with evolutionary information, we define a
medically-relevant "cancer network core," strongly enriched for common cancer
mutations ($p=2\times 10^{-14}$). This could aid in identifying driver
mutations and therapeutic targets.


  The common understanding of protein evolution has been that neutral or
slightly deleterious mutations are fixed by random drift, and evolutionary rate
is determined primarily by the proportion of neutral mutations. However, recent
studies have revealed that highly expressed genes evolve slowly because of
fitness costs due to misfolded proteins. Here we study selection maintaining
protein stability.
  Protein fitness is taken to be $s = \kappa \exp(\beta\Delta G) (1 -
\exp(\beta\Delta\Delta G))$, where $s$ and $\Delta\Delta G$ are selective
advantage and stability change of a mutant protein, $\Delta G$ is the folding
free energy of the wild-type protein, and $\kappa$ represents protein abundance
and indispensability. The distribution of $\Delta\Delta G$ is approximated to
be a bi-Gaussian function, which represents structurally slightly- or
highly-constrained sites. Also, the mean of the distribution is negatively
proportional to $\Delta G$.
  The evolution of this gene has an equilibrium ($\Delta G_e$) of protein
stability, the range of which is consistent with experimental values. The
probability distribution of $K_a/K_s$, the ratio of nonsynonymous to synonymous
substitution rate per site, over fixed mutants in the vicinity of the
equilibrium shows that nearly neutral selection is predominant only in
low-abundant, non-essential proteins of $\Delta G_e > -2.5$ kcal/mol. In the
other proteins, positive selection on stabilizing mutations is significant to
maintain protein stability at equilibrium as well as random drift on slightly
negative mutations, although the average $\langle K_a/K_s \rangle$ is less than
1. Slow evolutionary rates can be caused by high protein
abundance/indispensability, which produces positive shifts of $\Delta\Delta G$
through decreasing $\Delta G_e$, and by strong structural constraints, which
directly make $\Delta\Delta G$ more positive.


  In spite of their many facets, the phenomena of autoimmunity and
immunodeficiency seem to be related to each other through the subtle links
connecting retroviral mutation and action to immune response and adaptation. In
a previous work, we introduced a network model of how a set of interrelated
genotypes (called a quasispecies, in the stationary state) and a set of
interrelated idiotypes (an idiotypic network) interact. That model, which does
not cover the case of a retroviral quasispecies, was instrumental for the study
of quasispecies survival when confronting the immune system and led to the
conclusion that, unlike what happens when a quasispecies is left to evolve by
itself, letting genotypes mutate too infrequently leads to the destruction of
the quasispecies. Here we extend that genotype-idiotype interaction model by
the addition of a further parameter ($\nu$) to account for the action of
retroviruses (i.e., the destruction of idiotypes by genotypes). We give
simulation results within a suitable parameter niche, highlighting the issues
of quasispecies survival and of the onset of autoimmunity through the
appearance of the so-called pathogenic idiotypes. Our main findings refer to
how $\nu$ and $\lambda$, a parameter describing the rate at which idiotypes get
stimulated, relate to each other. While for $\nu>\lambda$ the quasispecies
survives at the expense of weakening the immune system significantly or even
destroying it, for $\nu<\lambda$ the fittest genotypes of the quasispecies
become mimicked inside the immune system as pathogenic idiotypes. The latter is
in agreement with the current understanding of the HIV quasispecies.


  Analyses of a set of 47 fossil and 4 modern skulls using phylogenetic
geometric morphometric methods corroborate and refine earlier results. These
include evidence that the African Iwo Eleru skull, only about 12,000 years old,
indeed represents a new species of near human. In contrast, the earliest known
anatomically modern human skull, Qafzeh 9, the skull of Eve from
Israel/Palestine, is validated as fully modern in form. Analyses clearly show
evidence of archaic introgression into Gravettian, pre_Gravettian, Qafzeh, and
Upper Cave (China) populations of near modern humans, and in about that order
of increasing archaic content. The enigmatic Saldahna (Elandsfontein) skull
emerges as a probable first representative of that lineage, which exclusive of
Neanderthals that, eventually lead to modern humans. There is also evidence
that the poorly dated Kabwe (Broken Hill) skull represents a much earlier
distinct lineage. The clarity of the results bode well for quantitative
statistical phylogenetic methods making significant inroads in the stalemates
of paleoanthropology.


  Malfunction of the system which regulates the bloodflow in the brain is a
major cause of stroke and dementia, costing many lives and many billions of
pounds each year in the UK alone. This regulatory system, known as cerebral
autoregulation, has been the subject of much experimental and mathematical
investigation yet our understanding of it is still quite limited. One area in
which our understanding is particularly lacking is that of the role of nitric
oxide, understood to be a potent vasodilator. The interactions of nitric oxide
with the better understood myogenic response remain un-modelled and poorly
understood. In this thesis we present a novel model of the arteriolar control
mechanism, comprising a mixture of well-established and new models of
individual processes, brought together for the first time. We show that this
model is capable of reproducing experimentally observed behaviour very closely
and go on to investigate its stability in the context of the vasculature of the
whole brain. In conclusion we find that nitric oxide, although it plays a
central role in determining equilibrium vessel radius, is unimportant to the
dynamics of the system and its responses to variation in arterial blood
pressure. We also find that the stability of the system is very sensitive to
the dynamics of Ca$^{2+}$ within the muscle cell, and that self-sustaining
Ca$^{2+}$ waves are not necessary to cause whole-vessel radius oscillations
consistent with vasomotion.


  Epidemic spread in single-host systems strongly depends on the population's
contact network. However, little is known regarding the spread of epidemics
across networks representing populations of multiple hosts. We explored
cross-species transmission in a multilayer network where layers represent
populations of two distinct hosts, and disease can spread across intralayer
(within-host) and interlayer (between-host) edges. We developed an analytic
framework for the SIR epidemic model to examine the effect of (i) source of
infection and (ii) between-host asymmetry in infection probabilities, on
disease risk. We measured risk as outbreak probability and outbreak size in a
focal host, represented by one network layer. Numeric simulations were used to
validate the analytic formulations. We found that outbreak probability is
determined by a complex interaction between source of infection and
between-host infection probabilities, whereas outbreak size is mainly affected
by the non-focal host to focal host infection probability alone. Hence,
inter-specific asymmetry in infection probabilities shapes disease dynamics in
multihost networks. These results expand current theory of monolayer networks,
where outbreak size and probability are considered equal, highlighting the
importance of considering multiple measures of disease risk. Our study advances
understanding of multihost systems and non-biological systems with asymmetric
flow rates.


  A three population system with a top predator population, i.e. the
herbivores, and two prey populations, grass and trees, is considered to model
the interaction of herbivores with natural resources. We apply the model for
four natural mountain parks in Northern Italy, three located in the Eastern
Alps, two of which in the Dolomites and one in the Julian Alps, and one in the
Maritime Alps, Northwest Italy. The simulations, based on actual data gathered
from contacts with rangers and parks administrators, field samplings and
published material, provide useful information on the behavior of the
vegetation-wild herbivores interactions and the possible medium-long term
evolution of these ecosystems. At the same time they show that these ecosystems
are in a very delicate situation, for which the animal populations could become
extinguished in case of adverse environmental conditions. The determination of
the so called sensitivity surfaces support our findings and indicate some
possible preventive measures to the park admistrators.


  We study two specific measures of quality of chemical reaction networks,
Precision and Sensitivity. The two measures arise in the study of sensory
adaptation, in which the reaction network is viewed as an input-output system.
Given a step change in input, Sensitivity is a measure of the magnitude of the
response, while Precision is a measure of the degree to which the system
returns to its original output for large time. High values of both are
necessary for high-quality adaptation.
  We focus on reaction networks without dissipation, which we interpret as
detailed-balance, mass-action networks. We give various upper and lower bounds
on the optimal values of Sensitivity and Precision, characterized in terms of
the stoichiometry, by using a combination of ideas from matroid theory and
differential-equation theory.
  Among other results, we show that this class of non-dissipative systems
contains networks with arbitrarily high values of both Sensitivity and
Precision. This good performance does come at a cost, however, since certain
ratios of concentrations need to be large, the network has to be extensive, or
the network should show strongly different time scales.


  The self-assembly of proteins into $\beta$-sheet-rich amyloid fibrils has
been observed to occur with sigmoidal kinetics, indicating that the system
initially is trapped in a metastable state. Here, we use a minimal
lattice-based model to explore the thermodynamic forces driving amyloid
formation in a finite canonical ($NVT$) system. By means of
generalized-ensemble Monte Carlo techniques and a semi-analytical method, the
thermodynamic properties of this model are investigated for different sets of
intersheet interaction parameters. When the interactions support lateral growth
into multi-layered fibrillar structures, an evaporation/condensation transition
is observed, between a supersaturated solution state and a thermodynamically
distinct state where small and large fibril-like species exist in equilibrium.
Intermediate-size aggregates are statistically suppressed. These properties do
not hold if aggregate growth is one-dimensional.


  This paper studies how patterns derived from a system of reaction-diffusion
equations may vary significantly depending upon boundary and initial
conditions, as well as in the spatial dependence of the coefficients involved.
From an extensive numerical study of the BVAM model, we demonstrate that the
geometric pattern of a reaction-diffusion system is not uniquely determined by
the value of the parameters in the equation. From this result, we suggest that
the variability of patterns among individuals of the same species may have its
roots in this sensitivity. Furthermore, this study analyzes briefly how the
inclusion of the advection and the space dependency in the parameters of the
model influences the forms of a specific pattern. The results of this study are
compared to the skin patterns that appear in Pseudoplatystom} fishes.


  One of the most important tasks in understanding the complex spatial
organization of the genome consists in extracting information about this
spatial organization, the function and structure of chromatin topological
domains from existing experimental data, in particular, from genome
colocalization (Hi-C) matrices. Here we present an algorithm allowing to reveal
the underlying hierarchical domain structure of a polymer conformation from
analyzing the modularity of colocalization matrices. We also test this
algorithm on several model polymer structures: equilibrium globules, random
fractal globules and regular fractal (Peano) conformations. We define what we
call a spectrum of cluster borders, and show that these spectra behave
strikingly differently for equilibrium and fractal conformations, allowing us
to suggest an additional criterion to identify fractal polymer conformations.


  Popular online enrichment analysis tools from the field of molecular systems
biology provide users with the ability to submit their experimental results as
gene sets for individual analysis. Such queries are kept private, and have
never before been considered as a resource for integrative analysis. By
harnessing gene set query submissions from thousands of users, we aim to
discover biological knowledge beyond the scope of an individual study. In this
work, we investigated a large collection of gene sets submitted to the tool
Enrichr by thousands of users. Based on co-occurrence, we constructed a global
gene-gene association network. We interpret this inferred network as providing
a summary of the structure present in this crowdsourced gene set library, and
show that this network recapitulates known protein-protein interactions and
functional associations between genes. This finding implies that this network
also offers predictive value. Furthermore, we visualize this gene-gene
association network using a new edge-pruning algorithm that retains both the
local and global structures of large-scale networks. Our ability to make
predictions for currently unknown gene associations, that may not be captured
by individual researchers and data sources, is a demonstration of the potential
of harnessing collective knowledge from users of popular tools in the field of
molecular systems biology.


  Phylogenetic networks are a generalization of phylogenetic trees that are
used to represent reticulate evolution. Unrooted phylogenetic networks form a
special class of such networks, which naturally generalize unrooted
phylogenetic trees. In this paper we define two operations on unrooted
phylogenetic networks, one of which is a generalization of the well-known
nearest-neighbor interchange (NNI) operation on phylogenetic trees. We show
that any unrooted phylogenetic network can be transformed into any other such
network using only these operations. This generalizes the well-known fact that
any phylogenetic tree can be transformed into any other such tree using only
NNI operations. It also allows us to define a generalization of tree space and
to define some new metrics on unrooted phylogenetic networks. To prove our main
results, we employ some fascinating new connections between phylogenetic
networks and cubic graphs that we have recently discovered. Our results should
be useful in developing new strategies to search for optimal phylogenetic
networks, a topic that has recently generated some interest in the literature,
as well as for providing new ways to compare networks.


  Although social network analysis is a promising tool to study the structure
and dynamics of wildlife communities, the current methods require costly and
detailed network data, which often are not available over long time periods
(e.g. decades). The present work aims to resolve this issue by developing a new
methodology that requires much less detailed data and it relies on
well-established Lotka-Volterra models. Using the long-term abundance data
(1973-2003) of northeastern Kansas rodents, the changes in the magnitude and
direction of interactions (e.g., changes from cooperative behavior to
competitive behavior or changes in the magnitude of competition) are
quantified.


  Lifetimes of bound states of protein complexes or biomolecule folded states
typically decrease when subject to mechanical force. However, a plethora of
biological systems exhibit the counter-intuitive phenomenon of catch bonding,
where non-covalent bonds become stronger under externally applied forces. The
quest to understand the origin of catch-bond behavior has lead to the
development of phenomenological and microscopic theories that can
quantitatively recapitulate experimental data. Here, we assess the successes
and limitations of such theories in explaining experimental data. The most
widely applied approach is a phenomenological two-state model, which fits all
of the available data on a variety of complexes: actomyosin,
kinetochore-microtubule, selectin-ligand, and cadherin-catenin binding to
filamentous actin. With a primary focus on the selectin family of cell-adhesion
complexes, we discuss the positives and negatives of phenomenological models
and the importance of evaluating the physical relevance of fitting parameters.
We describe a microscopic theory for selectins, which provides a structural
basis for catch bonds and predicts a crucial allosteric role for residues
Asn82--Glu88. We emphasize the need for new theories and simulations that can
mimic experimental conditions, given the complex response of cell adhesion
complexes to force and their potential role in a variety of biological
contexts.


  Despite the central role that antibodies play in the adaptive immune system
and in biotechnology, much remains unknown about the quantitative relationship
between an antibody's amino acid sequence and its antigen binding affinity.
Here we describe a new experimental approach, called Tite-Seq, that is capable
of measuring binding titration curves and corresponding affinities for
thousands of variant antibodies in parallel. The measurement of titration
curves eliminates the confounding effects of antibody expression and stability
that arise in standard deep mutational scanning assays. We demonstrate Tite-Seq
on the CDR1H and CDR3H regions of a well-studied scFv antibody. Our data shed
light on the structural basis for antigen binding affinity and suggests a role
for secondary CDR loops in establishing antibody stability. Tite-Seq fills a
large gap in the ability to measure critical aspects of the adaptive immune
system, and can be readily used for studying sequence-affinity landscapes in
other protein systems.


  The US and Hungarian statistical records of the years 1900 and 1896,
respectively, before the dramatic medical advances, show 32% and 27% deaths
attributable to infections, whereas only 5% and 2% due to cancer. These data
can be interpreted to mean that (i) the immune system evolved for purging
nascent selfish cells, which establish natural chimerism littering the soma and
the germline by conspecific alien cells and (ii) defense against pathogens that
represent xenogeneic aliens appeared later in evolution.
  `Liberating' T cells from the semantic trap of immunity and the shackles of
the `two-signal' model of T cell activation, we point out theoretical grounds
that the immune response to cancer is conceptually different from the immune
response to infection. We argue for a one-signal model (with stochastic
influences) as the explanation for T cell activation in preference to the
widely accepted two-signal model of co-stimulation. Convincing evidence for our
one-signal model emerged from the widespread autoimmune adverse events in 64.2%
of advanced melanoma patients treated with the anti-CTLA-4 antibody
(ipilimumab) that blocks an immune checkpoint. Harnessing the unleashed
autoimmune power of T cells could be rewarding to defeat cancer. Assuming that
immunization against isogeneic tumors also would be effective is a fallacy.


  We explore the range of probabilistic behaviours that can be engineered with
Chemical Reaction Networks (CRNs). We show that at steady state CRNs are able
to "program" any distribution with finite support in $\mathbb{N}^m$, with $m
\geq 1$. Moreover, any distribution with countable infinite support can be
approximated with arbitrarily small error under the $L^1$ norm. We also give
optimized schemes for special distributions, including the uniform
distribution. Finally, we formulate a calculus to compute on distributions that
is complete for finite support distributions, and can be compiled to a
restricted class of CRNs that at steady state realize those distributions.


  The three-dimensional organisation of chromosomes can be probed using methods
such as Capture-C. However it is unclear how such population level data relates
to the organisation within a single cell, and the mechanisms leading to the
observed interactions are still largely obscure. We present a polymer modelling
scheme based on the assumption that chromosome architecture is maintained by
protein bridges which form chromatin loops. To test the model we perform FISH
experiments and also compare with Capture-C data. Starting merely from the
locations of protein binding sites, our model accurately predicts the
experimentally observed chromatin interactions, revealing a population of 3D
conformations.


  Spatial self-organization emerges in distributed systems exhibiting local
interactions when nonlinearities and the appropriate propagation of signals are
at work. These kinds of phenomena can be modeled with different frameworks,
typically cellular automata or reaction-diffusion systems. A different class of
dynamical processes involves the correlated movement of agents over space,
which can be mediated through chemotactic movement or minimization of cell-cell
interaction energy. A classic example of the latter is given by the formation
of spatially segregated assemblies when cells display differential adhesion.
Here we consider a new class of dynamical models, involving cell adhesion among
two stochastically exchangeable cell states as a minimal model capable of
exhibiting well-defined, ordered spatial patterns. Our results suggest that a
whole space of pattern-forming rules is hosted by the combination of physical
differential adhesion and the value of probabilities modulating cell phenotypic
switching, showing that Turing-like patterns can be obtained without resorting
to reaction-diffusion processes. If the model is expanded allowing cells to
proliferate and die in an environment where diffusible nutrient and toxic waste
are at play, different phases are observed, characterized by regularly spaced
patterns. The analysis of the parameter space reveals that certain phases reach
higher population levels than other modes of organization. A detailed
exploration of the mean-field theory is also presented. Finally we let
populations of cells with different adhesion matrices compete for reproduction,
showing that, in our model, structural organization can improve the fitness of
a given cell population. The implications of these results for ecological and
evolutionary models of pattern formation and the emergence of multicellularity
are outlined.


  Contig assembly is the first stage that most assemblers solve when
reconstructing a genome from a set of reads. Its output consists of contigs --
a set of strings that are promised to appear in any genome that could have
generated the reads. From the introduction of contigs 20 years ago, assemblers
have tried to obtain longer and longer contigs, but the following question was
never solved: given a genome graph $G$ (e.g. a de Bruijn, or a string graph),
what are all the strings that can be safely reported from $G$ as contigs? In
this paper we finally answer this question, and also give a polynomial time
algorithm to find them. Our experiments show that these strings, which we call
omnitigs, are 66% to 82% longer on average than the popular unitigs, and 29% of
dbSNP locations have more neighbors in omnitigs than in unitigs.


  The chemical master equation (CME) is the exact mathematical formulation of
chemical reactions occurring in a dilute and well-mixed volume. The
reaction-diffusion master equation (RDME) is a stochastic description of
reaction-diffusion processes on a spatial lattice, assuming well-mixing only on
the length scale of the lattice. It is clear that, for the sake of consistency,
the solution of the RDME of a chemical system should converge to the solution
of the CME of the same system in the limit of fast diffusion: indeed, this has
been tacitly assumed in most literature concerning the RDME. We show that, in
the limit of fast diffusion, the RDME indeed converges to a master equation,
but not necessarily the CME. We introduce a class of propensity functions, such
that if the RDME has propensities exclusively of this class then the RDME
converges to the CME of the same system; while if the RDME has propensities not
in this class then convergence is not guaranteed. These are revealed to be
elementary and non-elementary propensities respectively. We also show that
independent of the type of propensity, the RDME converges to the CME in the
simultaneous limit of fast diffusion and large volumes. We illustrate our
results with some simple example systems, and argue that the RDME cannot be an
accurate description of systems with non-elementary rates.


  Proper balance between exploitation and exploration is what makes good
decisions, which achieve high rewards like payoff or evolutionary fitness. The
Infomax principle postulates that maximization of information directs the
function of diverse systems, from living systems to artificial neural networks.
While specific applications are successful, the validity of information as a
proxy for reward remains unclear. Here, we consider the multi-armed bandit
decision problem, which features arms (slot-machines) of unknown probabilities
of success and a player trying to maximize cumulative payoff by choosing the
sequence of arms to play. We show that an Infomax strategy (Info-p) which
optimally gathers information on the highest mean reward among the arms
saturates known optimal bounds and compares favorably to existing policies. The
highest mean reward considered by Info-p is not the quantity actually needed
for the choice of the arm to play, yet it allows for optimal tradeoffs between
exploration and exploitation.


  The yeast galactose network has provided many insights into how eukaryotic
gene circuits regulate metabolic function. However, there is currently no
consensus model of the network that incorporates protein dilution due to
cellular growth. We address this by adapting a well-known model and having it
account for growth benefit and burden due to expression of the network
proteins. Modifying the model to incorporate galactose transport and basal
Gal1p production allows us to better reproduce experimental observations.
Incorporating the growth rate effect demonstrates how the native network can
optimize growth in different galactose environments. These findings advance our
quantitative understanding of this gene network, and implement a general
approach for analysing the balance between growth costs and benefits in a range
of metabolic control networks.


  The solution space of genome-scale models of cellular metabolism provides a
map between physically viable flux configurations and cellular metabolic
phenotypes described, at the most basic level, by the corresponding growth
rates. By sampling the solution space of E. coli's metabolic network, we show
that empirical growth rate distributions recently obtained in experiments at
single-cell resolution can be explained in terms of a trade-off between the
higher fitness of fast-growing phenotypes and the higher entropy of
slow-growing ones. Based on this, we propose a minimal model for the evolution
of a large bacterial population that captures this trade-off. The scaling
relationships observed in experiments encode, in such frameworks, for the same
distance from the maximum achievable growth rate, the same degree of growth
rate maximization, and/or the same rate of phenotypic change. Being grounded on
genome-scale metabolic network reconstructions, these results allow for
multiple implications and extensions in spite of the underlying conceptual
simplicity.


  We study a large data set of protein structure ensembles of very diverse
sizes determined by nuclear magnetic resonance. By examining the
distance-dependent correlations in the displacement of residues pairs and
conducting finite size scaling analysis it was found that the correlations and
susceptibility behave as in systems near a critical point implying that, at the
native state, the motion of each amino acid residue is felt by every other
residue up to the size of the protein molecule. Furthermore certain protein's
shapes corresponding to maximum susceptibility were found to be more probable
than others. Overall the results suggest that the protein's native state is
critical, implying that despite being posed near the minimum of the energy
landscape, they still preserve their dynamic flexibility.


  Over the last decade new technologies for making large numbers of
fine-grained measurements have led to the surprising discovery that many
biological systems sit near a critical point. These systems are potentially
more adaptive in that small changes to component behavior can induce
large-scale changes in aggregate structure and function. Accounting for
criticality remains a challenge as sensitivity to perturbation suggests a lack
of robustness. Furthermore, change induced by perturbation may not be adaptive.
Complicating matters further critical phenomena can result from
history-dependent stochastic processes. A question central to distinguishing
among these conflicting views of criticality is to what degree criticality can
be controlled by the components of the system. We address the control of
criticality using data on conflict dynamics and fight sizes from an animal
society model system (Macaca nemestrina, n=48). The system is fundamentally
finite so we operationalize criticality in information theoretic terms using
Fisher information and a measure of instability. We analyze criticality using
empirically-grounded equilibrium (maximum entropy) and dynamic (branching
process) models of the monkeys' fight-joining behavior. We find that (1) this
heterogeneous, socially organized system, like homogeneous, spatial systems
(flocks and schools), sits near a critical point, (2) the contributions
individuals make to how critical the system is can be quantified and vary, and
(3) the distance from the critical point (DFC) can be controlled through
biologically plausible mechanisms operating on this heterogeneity. These
mechanisms include third-party policing, which dampens fight participation of
the individuals with the largest effect on DFC. Control of DFC allows
biological systems to balance the tradeoff between robustness and need for
rapid change.


  Circadian clocks are the central timekeepers of life, allowing cells to
anticipate changes between day and night. Experiments in recent years have
revealed that circadian clocks can be highly stable, raising the question how
reliably they can be read out. Here, we combine mathematical modeling with
information theory to address the question how accurately a cell can infer the
time from an ensemble of protein oscillations, which are driven by a circadian
clock. We show that the precision increases with the number of oscillations and
their amplitude relative to their noise. Our analysis also reveals that their
exists an optimal phase relation that minimizes the error in the estimate of
time, which depends on the relative noise levels of the protein oscillations.
Lastly, our work shows that cross-correlations in the noise of the protein
oscillations can enhance the mutual information, which suggests that
cross-regulatory interactions between the proteins that read out the clock can
be beneficial for temporal information transmission.


  Hand-hygiene compliance and contacts of health-care workers largely determine
the potential paths of pathogen transmission in hospital wards. We explored how
the combination of data collected by two automated infrastructures based on
wearable sensors and recording (i) use of hydro-alcoholic solution and (ii)
contacts of health-care workers provide an enhanced view of the risk of
transmission events in the ward. We perform a proof-of-concept observational
study. Detailed data on contact patterns and hand-hygiene compliance of
health-care workers were collected by wearable sensors over 12 days in an
infectious disease unit of a hospital in Marseilles, France. 10837 contact
events among 10 doctors, 4 nurses, 4 nurses' aids and 4 housekeeping staff were
recorded during the study. Most contacts took place among medical doctors.
Aggregate contact durations were highly heterogeneous and the resulting contact
network was highly structured. 510 visits of health-care workers to patients'
rooms were recorded, with a low rate of hand-hygiene compliance. Both data sets
were used to construct histories and statistics of contacts informed by the use
of hydro-alcoholic solution, or lack thereof, of the involved health-care
workers. Hand-hygiene compliance data strongly enrich the information
concerning contacts among health-care workers, by assigning a 'safe' or
'at-risk' value to each contact. The global contact network can thus be divided
into 'at-risk' and 'safe' contact networks. The combined data could be of high
relevance for outbreak investigation and to inform data-driven models of
nosocomial disease spread.


  The isolation-with-migration (IM) model is a common tool to make inferences
about the presence of gene flow during speciation, using polymorphism data.
However, Becquet and Przeworski (2009) report that the parameter estimates
obtained by fitting the IM model are very sensitive to the model's assumptions,
including the assumption of constant gene flow until the present. This paper is
concerned with the isolation-with-initial-migration (IIM) model of
Wilkinson-Herbots (2012), which drops precisely this assumption. In the IIM
model, one ancestral population divides into two descendant subpopulations,
between which there is an initial period of gene flow and a subsequent period
of isolation. We derive a fast method of fitting an extended version of the IIM
model, which allows for asymmetric gene flow and unequal subpopulation sizes.
This is a maximum-likelihood method, applicable to observations on the number
of different sites between pairs of DNA sequences from a large number of
independent loci. In addition to obtaining parameter estimates, our method can
also be used to distinguish between alternative models representing different
evolutionary cenarios, by means of likelihood ratio tests. We illustrate the
procedure on pairs of Drosophila sequences from approximately 30,000 loci. The
computing time needed to fit the most complex version of the model to this data
set is only a couple of minutes. The code to fit the IIM model can be found in
the supplementary files of this paper.


  A plethora of natural, artificial and social systems exist which do not
belong to the Boltzmann-Gibbs (BG) statistical-mechanical world, based on the
standard additive entropy $S_{BG}$ and its associated exponential BG factor.
Frequent behaviors in such complex systems have been shown to be closely
related to $q$-statistics instead, based on the nonadditive entropy $S_q$ (with
$S_1=S_{BG}$), and its associated $q$-exponential factor which generalizes the
usual BG one. In fact, a wide range of phenomena of quite different nature
exist which can be described and, in the simplest cases, understood through
analytic (and explicit) functions and probability distributions which exhibit
some universal features. Universality classes are concomitantly observed which
can be characterized through indices such as $q$. We will exhibit here some
such cases, namely concerning the distribution of inter-occurrence (or
inter-event) times in the areas of finance, earthquakes and genomes.


  We introduce an heterogeneous nonlinear $q$-voter model with zealots and two
types of susceptible voters, and study its non-equilibrium properties when the
population is finite and well mixed. In this two-opinion model, each individual
supports one of two parties and is either a zealot or a susceptible voter of
type $q_1$ or $q_2$. While here zealots never change their opinion, a
$q_i$-susceptible voter ($i=1,2$) consults a group of $q_i$ neighbors at each
time step, and adopts their opinion if all group members agree. We show that
this model violates the detailed balance whenever $q_1 \neq q_2$ and has
surprisingly rich properties. Here, we focus on the characterization of the
model's non-equilibrium stationary state (NESS) in terms of its probability
distribution and currents in the distinct regimes of low and high density of
zealotry. We unveil the NESS properties in each of these phases by computing
the opinion distribution and the circulation of probability currents, as well
as the two-point correlation functions at unequal times (formally related to a
"probability angular momentum"). Our analytical calculations obtained in the
realm of a linear Gaussian approximation are compared with numerical results.


  We consider the extension of the `Moran effect', where correlated noise
generates synchrony between isolated single species populations, to the study
of synchrony between populations embedded in multi-species communities. In
laboratory experiments on complex microbial communities, comprising both
predators (protozoa) and prey (bacteria), we observe synchrony in abundances
between isolated replicates. A breakdown in synchrony occurs for both predator
and prey as the reactor dilution rate increases, which corresponds to both an
increased rate of input of external resources and an increased effective
mortality though washout. The breakdown is more rapid, however, for the lower
trophic level. We can explain this phenomenon using a mathematical framework
for determining synchrony between populations in multi-species communities at
equilibrium. We assume that there are multiple sources of environmental noise
with different degrees of correlation that affect the individual species
population dynamics differently. The deterministic dynamics can then influence
the degree of synchrony between species in different communities. In the case
of a stable equilibrium community synchrony is controlled by the eigenvalue
with smallest negative real part. Intuitively fluctuations are minimally damped
in this direction. We show that the experimental observations are consistent
with this framework but only for multiplicative noise.


  Applying the constraints dictated by the principle of detailed balance, we
analyze a recent proposal for spontaneous mirror symmetry breaking (SMSB) based
on enantioselective autocatalysis coupled to a linear decay of the enantiomers
and in the presence of reaction noise. We find the racemic state is the final
stable outcome for both deterministic as well as for stochastic dynamics, and
for both well-mixed and small spatially-coupled systems. The racemic outcome
results even when the autocatalytic cycles are driven irreversibly by external
reagents, in manifestly non-equilibrium conditions. Our findings suggest that
first-order autocatalysis coupled to reactions involving \textit{non-linear}
heterochiral dynamics is a necessary pre-condition for any mechanism purporting
to lead to molecular homochirality.


  A major challenge in molecular simulations is to describe
denaturant-dependent folding of proteins order to make direct comparisons with
{\it in vitro} experiments. We use the molecular transfer model, which is
currently the only method that accomplishes this goal albeit
phenomenologically, to quantitatively describe urea-dependent folding of PDZ
domain, which plays a significant role in molecular recognition and signaling.
Experiments show that urea-dependent unfolding rates of the PDZ2 domain exhibit
a downward curvature at high urea concentrations, which has been interpreted by
invoking the presence of a sparsely populated high energy intermediate.
Simulations using the MTM and a coarse-grained model of PDZ2 are used to show
that the intermediate, which has some native-like character, is present in
equilibrium both in the presence and absence of urea. The free energy profiles
show that there are two barriers separating the folded and unfolded states.
Structures of the transition state ensembles, ($TSE1$ separating the unfolded
and $I_{EQ}$ and $TSE2$ separating $I_{EQ}$ and the native state), determined
using the $P_{fold}$ method, show that $TSE1$ is expanded; $TSE2$ and
native-like. Folding trajectories reveal that PDZ2 folds by parallel routes. In
one pathway folding occurs exclusively through $I_1$, which resembles $I_{EQ}$.
In a fraction of trajectories, constituting the second pathway, folding occurs
through a combination of $I_{1}$ and a kinetic intermediate. The radius of
gyration ($R_g^{U}$) of the unfolded state is more compact (by $\sim$ 9\%)
under native conditions. Decrease in $R_g^{U}$ occurs on the time scale on the
order of utmost $\sim$ 20 $\mu s$. The modest decrease in $R_g^{U}$ and the
rapid collapse suggest that high spatial and temporal resolution are needed to
detect compaction in finite-sized proteins.


  The search to understand how the HIV virus spreads inside the human body and
how the immune response works to control it has motivated studies related to
Mathematical Immunology. Actually, researches include the idea of mathematical
models representing the dynamics of healthy and infected cell populations and
focusing on mechanisms used by HIV to invade target-host cells, viral
dissemination (which leads to depletion of the T-cell pool and collapse of the
immune system), and impairment of immune response. In this work, we show the
importance of specific cells of immune response, as infection vectors involved
in the dynamics of viral proliferation within an untreated patient, by using an
ordinary differential equation model in which we considered that the virus
infected target-cells such as macrophages, dendritic cells, and lymphocytes
TCD4 and TCD8 populations. In conclusion, we demonstrate the importance of each
cell-host and the threshold of viral establishment and posterior spread based
on the presence of infected macrophages and dendritic cells in
antigen-presenting processing, which leads to new infections (by wild or mutant
virions) after immune response activation episodes. We presented an $R_0$
expression that provides the major parameters of HIV infection. Additionally,
we suggest some possibilities of new targets for functional vaccines.


  Dual graphs have been applied to model RNA secondary structures. The purpose
of the paper is two-fold: we present new graph-theoretic properties of dual
graphs to validate the further analysis and classification of RNAs using these
topological representations; we also present a linear-time algorithm to
partition dual graphs into topological components called {\it blocks} and
determine if each block contains a {\it pseudoknot} or not. We show that a
block contains a pseudoknot if and only if the block has a vertex of degree $3$
or more; this characterization allows us to efficiently isolate smaller RNA
fragments and classify them as pseudoknotted or pseudoknot-free regions, while
keeping these sub-structures intact. Even though non-topological techniques to
detect and classify pseudoknots have been efficiently applied, structural
properties of dual graphs provide a unique perspective for the further analysis
of RNAs. Applications to RNA design can be envisioned since modular building
blocks with intact pseudoknots can be combined to form new
constructs.\end{abstract}


  Many models of genome rearrangement involve operations (e.g. inversions and
translocations) that are self-inverse, and hence generate a group acting on the
space of genomes. This gives a correspondence between genome arrangements and
the elements of a group, and consequently, between evolutionary paths and walks
on the Cayley graph. Many common methods for phylogeny reconstruction rely on
calculating the minimal distance between two genomes; this omits much of the
other information available from the Cayley graph. In this paper we begin an
exploration of some of this additional information, in particular describing
the phylogeny as a Steiner tree within the Cayley graph, and exploring the
"interval" between two genomes. While motivated by problems in systematic
biology, many of these ideas are of independent group-theoretic interest.


  Comprehension of metabolic pathways is considerably enhanced by metabolic
flux analysis (MFA-ILE) in isotope labeling experiments. The balance equations
are given by hundreds of algebraic (stationary MFA) or ordinary differential
equations (nonstationary MFA), and reducing the number of operations is
therefore a crucial part of reducing the computation cost. The main bottleneck
for deterministic algorithms is the computation of derivatives, particularly
for nonstationary MFA. In this article we explain how the overall
identification process may be speeded up by using the adjoint approach to
compute the gradient of the residual sum of squares. The proposed approach
shows significant improvements in terms of complexity and computation time when
it is compared with the usual (direct) approach. Numerical results are obtained
for the central metabolic pathways of Escherichia coli and are validated
against reference software in the stationary case. The methods and algorithms
described in this paper are included in the sysmetab software package
distributed under an Open Source license at
http://forge.scilab.org/index.php/p/sysmetab/.


  Stochastic phenotype switching has been suggested to play a beneficial role
in microbial populations by leading to the division of labour among cells, or
ensuring that at least some of the population survives an unexpected change in
environmental conditions. Here we use a computational model to investigate an
alternative possible function of stochastic phenotype switching - as a way to
adapt more quickly even in a static environment. We show that when a genetic
mutation causes a population to become less fit, switching to an alternative
phenotype with higher fitness (growth rate) may give the population enough time
to develop compensatory mutations that increase the fitness again. The
possibility of switching phenotypes can reduce the time to adaptation by orders
of magnitude if the "fitness valley" caused by the deleterious mutation is deep
enough. Our work has important implications for the emergence of
antibiotic-resistant bacteria. In line with recent experimental findings we
hypothesise that switching to a slower growing but less sensitive phenotype
helps bacteria to develop resistance by exploring a larger set of beneficial
mutations while avoiding deleterious ones.


  Many DNA profiles recovered from crime scene samples are of a quality that
does not allow them to be searched against, nor entered into, databases. We
propose a method for the comparison of profiles arising from two DNA samples,
one or both of which can have multiple donors and be affected by low DNA
template or degraded DNA. We compute likelihood ratios to evaluate the
hypothesis that the two samples have a common DNA donor, and hypotheses
specifying the relatedness of two donors. Our method uses a probability
distribution for the genotype of the donor of interest in each sample. This
distribution can be obtained from a statistical model, or we can exploit the
ability of trained human experts to assess genotype probabilities, thus
extracting much information that would be discarded by standard interpretation
rules. Our method is compatible with established methods in simple settings,
but is more widely applicable and can make better use of information than many
current methods for the analysis of mixed-source, low-template DNA profiles. It
can accommodate uncertainty arising from relatedness instead of or in addition
to uncertainty arising from noisy genotyping. We describe a computer program
GPMDNA, available under an open source license, to calculate LRs using the
method presented in this paper.


  Rooted phylogenetic networks are used to describe evolutionary histories that
contain non-treelike evolutionary events such as hybridization and horizontal
gene transfer. In some cases, such histories can be described by a phylogenetic
base-tree with additional linking arcs, which can for example represent gene
transfer events. Such phylogenetic networks are called tree-based. Here, we
consider two possible generalizations of this concept to nonbinary networks,
which we call tree-based and strictly-tree-based nonbinary phylogenetic
networks. We give simple graph-theoretic characterizations of tree-based and
strictly-tree-based nonbinary phylogenetic networks. Moreover, we show for each
of these two classes that it can be decided in polynomial time whether a given
network is contained in the class. Our approach also provides a new view on
tree-based binary phylogenetic networks. Finally, we discuss two examples of
nonbinary phylogenetic networks in biology and show how our results can be
applied to them.


  Effective population size characterizes the genetic variability in a
population and is a parameter of paramount importance in population genetics.
Kingman's coalescent process enables inference of past population dynamics
directly from molecular sequence data, and researchers have developed a number
of flexible coalescent-based models for Bayesian nonparametric estimation of
the effective population size as a function of time. A major goal of
demographic reconstruction is understanding the association between the
effective population size and potential explanatory factors. Building upon
Bayesian nonparametric coalescent-based approaches, we introduce a flexible
framework that incorporates time-varying covariates through Gaussian Markov
random fields. To approximate the posterior distribution, we adapt efficient
Markov chain Monte Carlo algorithms designed for highly structured Gaussian
models. Incorporating covariates into the demographic inference framework
enables the modeling of associations between the effective population size and
covariates while accounting for uncertainty in population histories.
Furthermore, it can lead to more precise estimates of population dynamics. We
apply our model to four examples. We reconstruct the demographic history of
raccoon rabies in North America and find a significant association with the
spatiotemporal spread of the outbreak. Next, we examine the effective
population size trajectory of the DENV-4 virus in Puerto Rico along with viral
isolate count data and find similar cyclic patterns. We compare the population
history of the HIV-1 CRF02_AG clade in Cameroon with HIV incidence and
prevalence data and find that the effective population size is more reflective
of incidence rate. Finally, we explore the hypothesis that the population
dynamics of musk ox during the Late Quaternary period were related to climate
change.


  Copy number alterations (CNAs) are thought to account for 85% of the
variation in gene expression observed among breast tumours. The expression of
cis-associated genes is impacted by CNAs occurring at proximal loci of these
genes, whereas the expression of trans-associated genes is impacted by CNAs
occurring at distal loci. While a majority of these CNA-driven genes
responsible for breast tumourigenesis are cis-associated, trans-associated
genes are thought to further abet the development of cancer and influence
disease outcomes in patients. Here we present a network-based approach that
integrates copy-number and expression profiles to identify putative cis- and
trans-associated genes in breast cancer pathogenesis. We validate these cis-
and trans-associated genes by employing them to subtype a large cohort of
breast tumours obtained from the METABRIC consortium, and demonstrate that
these genes accurately reconstruct the ten subtypes of breast cancer. We
observe that individual breast cancer subtypes are driven by distinct sets of
cis- and trans-associated genes. Among the cis-associated genes, we recover
several known drivers of breast cancer (e.g. CCND1, ERRB2, MDM2 and ZNF703) and
some novel putative drivers (e.g. BRF2 and SF3B3). siRNA-mediated knockdown of
BRF2 across a panel of breast cancer cell lines showed significant reduction
specifically in cell proliferation in HER2+ lines, thereby indicating that BRF2
could be a context-dependent oncogene and potentially targetable in these
lines. Among the trans-associated genes, we identify modules of immune-response
(CD2, CD19, CD38 and CD79B), mitotic/cell-cycle kinases (e.g. AURKB, MELK, PLK1
and TTK), and DNA-damage response genes (e.g. RFC4 and FEN1).


  For an eco-epidemic model with disease in the prey and periodic coefficients
it is conjectured in [Xingge Niu, Tailei Zhang, Zhidong Teng, The asymptotic
behavior of a nonautonomous eco-epidemic model with disease in the prey,
Applied Mathematical Modelling 35, 457-470 (2011)] that, when the infected prey
is permanent, there is a positive periodic orbit that is globally
asymptotically stable in the interior of $(\R^+)^3$. In this paper we prove the
existence part of the conjecture.


  Comparative transcriptomics has gained increasing popularity in genomic
research thanks to the development of high-throughput technologies including
microarray and next-generation RNA sequencing that have generated numerous
transcriptomic data. An important question is to understand the conservation
and differentiation of biological processes in different species. We propose a
testing-based method TROM (Transcriptome Overlap Measure) for comparing
transcriptomes within or between different species, and provide a different
perspective to interpret transcriptomic similarity in contrast to traditional
correlation analyses. Specifically, the TROM method focuses on identifying
associated genes that capture molecular characteristics of biological samples,
and subsequently comparing the biological samples by testing the overlap of
their associated genes. We use simulation and real data studies to demonstrate
that TROM is more powerful in identifying similar transcriptomes and more
robust to stochastic gene expression noise than Pearson and Spearman
correlations. We apply TROM to compare the developmental stages of six
Drosophila species, C. elegans, S. purpuratus, D. rerio and mouse liver, and
find interesting correspondence patterns that imply conserved gene expression
programs in the development of these species. The TROM method is available as
an R package on CRAN (http://cran.r-project.org/) with manuals and source codes
available at http://www.stat.ucla.edu/ jingyi.li/software-and-data/trom.html.


  We formulate models for the spread of infection on networks that are amenable
to analysis in the large population limit. We distinguish three different
levels: (1) binding sites, (2) individuals, and (3) the population. In the
tradition of Physiologically Structured Population Models, the formulation
starts on the individual level. Influences from the `outside world' on an
individual are captured by environmental variables. These environmental
variables are population level quantities. A key characteristic of the network
models is that individuals can be decomposed into a number of conditionally
independent components: each individual has a fixed number of `binding sites'
for partners. The Markov chain dynamics of binding sites are described by only
a few equations. In particular, individual-level probabilities are obtained
from binding-site-level probabilities by combinatorics while population-level
quantities are obtained by averaging over individuals in the population. Thus
we are able to characterize population-level epidemiological quantities, such
as $R_0$, $r$, the final size, and the endemic equilibrium, in terms of the
corresponding variables.


  The foundational concepts behind the persistence of ecological communities
have been based on two ecological properties: dynamical stability and
feasibility. The former is typically regarded as the capacity of a community to
return to an original equilibrium state after a perturbation in species
abundances, and is usually linked to the strength of interspecific
interactions. The latter is the capacity to sustain positive abundances on all
its constituent species, and is linked to both interspecific interactions and
species demographic characteristics. Over the last 40 years, theoretical
research in ecology has emphasized the search for conditions leading to the
dynamical stability of ecological communities, while the conditions leading to
feasibility have been overlooked. However, thus far, we have no evidence of
whether species interactions are more conditioned by the community's need to be
stable or feasible. Here, we introduce novel quantitative methods and use
empirical data to investigate the consequences of species interactions on the
dynamical stability and feasibility of mutualistic communities. First, we
demonstrate that the more nested the species interactions in a community are,
the lower the mutualistic strength that the community can tolerate without
losing dynamical stability. Second, we show that high feasibility in a
community can be reached either with high mutualistic strength or with highly
nested species interactions. Third, we find that during the assembly process of
a seasonal pollinator community located at The Zackenberg Research Station (NE
Greenland), a high feasibility is reached through the nested species
interactions established between newcomer and resident species. Our findings
imply that nested mutualistic communities promote feasibility over stability,
which may suggest that the former can be key for community persistence.


  The advent of accessible ancient DNA technology now allows the direct
ascertainment of allele frequencies in ancestral populations, thereby enabling
the use of allele frequency time series to detect and estimate natural
selection. Such direct observations of allele frequency dynamics are expected
to be more powerful than inferences made using patterns of linked neutral
variation obtained from modern individuals. We develop a Bayesian method to
make use of allele frequency time series data and infer the parameters of
general diploid selection, along with allele age, in non-equilibrium
populations. We introduce a novel path augmentation approach, in which we use
Markov chain Monte Carlo to integrate over the space of allele frequency
trajectories consistent with the observed data. Using simulations, we show that
this approach has good power to estimate selection coefficients and allele age.
Moreover, when applying our approach to data on horse coat color, we find that
ignoring a relevant demographic history can significantly bias the results of
inference. Our approach is made available in a C++ software package.


  Diffusion-based network models are widely used for protein function
prediction using protein network data and have been shown to outperform
neighborhood- and module-based methods. Recent studies have shown that
integrating the hierarchical structure of the Gene Ontology (GO) data
dramatically improves prediction accuracy. However, previous methods usually
either used the GO hierarchy to refine the prediction results of multiple
classifiers, or flattened the hierarchy into a function-function similarity
kernel. No study has taken the GO hierarchy into account together with the
protein network as a two-layer network model.
  We first construct a Bi-relational graph (Birg) model comprised of both
protein-protein association and function-function hierarchical networks. We
then propose two diffusion-based methods, BirgRank and AptRank, both of which
use PageRank to diffuse information on this two-layer graph model. BirgRank is
an application of traditional PageRank with fixed decay parameters. In
contrast, AptRank uses an adaptive mechanism to improve the performance of
BirgRank. We evaluate both methods in predicting protein function on yeast,
fly, and human datasets, and compare with four previous methods: GeneMANIA,
TMC, ProteinRank and clusDCA. We design three validation strategies: missing
function prediction, de novo function prediction, and guided function
prediction to comprehensively evaluate all six methods. We find that both
BirgRank and AptRank outperform the others, especially in missing function
prediction when using only 10% of the data for training.
  AptRank combines protein-protein associations and the GO function-function
hierarchy into a two-layer network model without flattening the hierarchy into
a similarity kernel. Introducing an adaptive mechanism to the traditional,
fixed-parameter model of PageRank greatly improves the accuracy of protein
function prediction.


  We present a stochastic framework to study signal transmission in a generic
two-step cascade $\text{S} \rightarrow \text{X} \rightarrow \text{Y}$. Starting
from a set of Langevin equations obeying Gaussian noise processes we calculate
the variance and covariance while considering both linear and nonlinear
production terms for different biochemical species of the cascade. These
quantities are then used to calculate the net synergy within the purview of
partial information decomposition. We show that redundancy in information
transmission is essentially an important consequence of Markovian property of
the two-step cascade motif. We also show that redundancy increases fidelity of
the signalling pathway.


  In this study we present an extension of the replicator equation with
diffusion to multiplex graphs. We derive an exact formula for the diffusion
term, which shows that, while diffusion is linear for numbers of agents, it is
necessary to account for non-linear terms when working with fractions of
individuals. We also derive the transition probabilities that give rise to such
macroscopic behavior, completing the bottom-up description. Finally, it shown
that the usual assumption of constant population sizes induces a hidden
selective pressure in the system.


  Identifying the critical domain size necessary for a population to persist is
an important question in ecology. Both demographic and environmental
stochasticity impact a population's ability to persist. Here we explore ways of
including this variability. We study populations which have traditionally been
modelled using a deterministic integrodifference equation (IDE) framework, with
distinct dispersal and sedentary stages. Individual based models (IBMs) are the
most intuitive stochastic analogues to IDEs but yield few analytic insights. We
explore two alternate approaches; one is a scaling up to the population level
using the Central Limit Theorem, and the other a variation on both
Galton-Watson branching processes and branching processes in random
environments. These branching process models closely approximate the IBM and
yield insight into the factors determining the critical domain size for a given
population subject to stochasticity.


  Because genomes are products of natural processes rather than intelligent
design, all genomes contain functional and nonfunctional parts. The fraction of
the genome that has no biological function is called rubbish DNA. Rubbish DNA
consists of junk DNA, i.e., the fraction of the genome on which selection does
not operate, and garbage DNA, i.e., sequences that lower the fitness of the
organism, but exist in the genome because purifying selection is neither
omnipotent nor instantaneous. In this chapter, I (1) review the concepts of
genomic function and functionlessness from an evolutionary perspective, (2)
present a precise nomenclature of genomic function, (3) discuss the evidence
for the existence of vast quantities of junk DNA within the human genome, (4)
discuss the mutational mechanisms responsible for generating junk DNA, (5)
spell out the necessary evolutionary conditions for maintaining junk DNA, (6)
outline various methodologies for estimating the functional fraction within the
genome, and (7) present a recent estimate for the functional fraction of our
genome.


  The two-parameter Poisson--Dirichlet diffusion, introduced in 2009 by Petrov,
extends the infinitely-many-neutral-alleles diffusion model, related to
Kingman's one-parameter Poisson--Dirichlet distribution and to certain
Fleming--Viot processes. The additional parameter has been shown to regulate
the clustering structure of the population, but is yet to be fully understood
in the way it governs the reproductive process. Here we shed some light on
these dynamics by formulating a $K$-allele Wright--Fisher model for a
population of size $N$, involving a uniform mutation pattern and a specific
state-dependent migration mechanism. Suitably scaled, this process converges in
distribution to a $K$-dimensional diffusion process as $N\to\infty$. Moreover,
the descending order statistics of the $K$-dimensional diffusion converge in
distribution to the two-parameter Poisson--Dirichlet diffusion as $K\to\infty$.
The choice of the migration mechanism depends on a delicate balance between
reinforcement and redistributive effects. The proof of convergence to the
infinite-dimensional diffusion is nontrivial because the generators do not
converge on a core. Our strategy for overcoming this complication is to prove
\textit{a priori} that in the limit there is no "loss of mass", i.e., that, for
each limit point of the sequence of finite-dimensional diffusions (after a
reordering of components by size), allele frequencies sum to one.


  The regeneration of forest resources is one of the main objectives of modern
forest management. The present book is based on my personal PhD thesis and aims
to identify detailed features of the regeneration process and the
characteristics of sapling interaction. The research involved evaluation of
regeneration structure, dimensional differentiation of saplings, spatial
pattern analysis - sapling aggregation and species association, and sapling
competition using distance-dependent methods. A new distance-dependent
dimensional differentiation index (IDIV) and a new geometrical criterion for
selecting neighbouring competitor saplings were created and applied.
Appropriate methods and techniques were identified for investigating
regeneration and computer technology was used intensively. Nine stand-alone
software programs were produced in order to achieve the objectives of the
thesis.


  A two-patch mathematical model of Dengue virus type 2 (DENV-2) that accounts
for vectors' vertical transmission and between patches human dispersal is
introduced. Dispersal is modeled via a Lagrangian approach. A host-patch
residence-times basic reproduction number is derived and conditions under which
the disease dies out or persists are established. Analytical and numerical
results highlight the role of hosts' dispersal in mitigating or exacerbating
disease dynamics. The framework is used to explore dengue dynamics using, as a
starting point, the 2002 outbreak in the state of Colima, Mexico.


  We develop theory and computational methods to investigate particle
inclusions embedded within curved lipid bilayer membranes. We consider the case
of spherical lipid vesicles where inclusion particles are coupled through (i)
intramembrane hydrodynamics, (ii) traction stresses with the external and
trapped solvent fluid, and (iii) intermonolayer slip between the two leaflets
of the bilayer. We investigate relative to flat membranes how the membrane
curvature and topology augment hydrodynamic responses. We show how both the
translational and rotational mobility of protein inclusions are effected by the
membrane curvature, ratio of intramembrane viscosity to solvent viscosity, and
inter-monolayer slip. For general investigations of many-particle dynamics, we
also discuss how our approaches can be used to treat the collective diffusion
and hydrodynamic coupling within spherical bilayers.


  Calculating solvent accessible surface areas (SASA) is a run-of-the-mill
calculation in structural biology. Although there are many programs available
for this calculation, there are no free-standing, open-source tools designed
for easy tool-chain integration. FreeSASA is an open source C library for SASA
calculations that provides both command-line and Python interfaces in addition
to its C API. The library implements both Lee and Richards' and Shrake and
Rupley's approximations, and is highly configurable to allow the user to
control molecular parameters, accuracy and output granularity. It only depends
on standard C libraries and should therefore be easy to compile and install on
any platform. The source code is freely available from
http://freesasa.github.io/. The library is well-documented, stable and
efficient. The command-line interface can easily replace closed source legacy
programs, with comparable or better accuracy and speed, and with some added
functionality.


  Dynamic models - often deterministic in nature - were used to estimate the
basic reproductive number, R_0, of the 2014-5 Ebola virus disease (EVD)
epidemic outbreak in West Africa. Estimates of R_0 were then used to project
the likelihood for large outbreak sizes, e.g., exceeding hundreds of thousands
of cases. Yet fitting deterministic models can lead to over-confidence in the
confidence intervals of the fitted R_0, and, in turn, the type and scope of
necessary interventions. In this manuscript we propose a hybrid
stochastic-deterministic method to estimate R_0 and associated confidence
intervals (CIs). The core idea is that stochastic realizations of an underlying
deterministic model can be used to evaluate the compatibility of candidate
values of R_0 with observed epidemic curves. The compatibility is based on
comparing the distribution of expected epidemic growth rates with the observed
epidemic growth rate given "process noise", i.e., arising due to stochastic
transmission, recovery and death events. By applying our method to reported EVD
case counts from Guinea, Liberia and Sierra Leone, we show that prior estimates
of R_0 based on deterministic fits appear to be more confident than analysis of
stochastic trajectories suggests should be possible. Moving forward, we
recommend including a hybrid stochastic-deterministic fitting procedure when
quantifying the full R_0 CI at the onset of an epidemic due to multiple sources
of noise.


  We study the spatio-temporal patterns of the proportion of influenza B out of
laboratory confirmations of both influenza A and B, with data from 139
countries and regions downloaded from the FluNet compiled by the World Health
Organization, from January 2006 to October 2015, excluding 2009. We restricted
our analysis to 34 countries that reported more than 2000 confirmations for
each of types A and B over the study period. We find that Pearson's correlation
is 0.669 between effective distance from Mexico and influenza B proportion
among the countries from January 2006 to October 2015. In the United States,
influenza B proportion in the pre-pandemic period (2003-2008) negatively
correlated with that in the post-pandemic era (2010-2015) at the regional
level. Our study limitations are the country-level variations in both
surveillance methods and testing policies. Influenza B proportion displayed
wide variations over the study period. Our findings suggest that even after
excluding 2009's data, the influenza pandemic still has an evident impact on
the relative burden of the two influenza types. Future studies could examine
whether there are other additional factors. This study has potential
implications in prioritizing public health control measures.


  Observation times of freely diffusing single molecules in solution are
limited by the photophysics of the attached fluorescence markers and by a small
observation volume in the femtolitre range that is required for a sufficient
signal-to-background ratio. To extend diffusion-limited observation times
through a confocal detection volume, A. E. Cohen and W. E. Moerner have
invented and built the ABELtrap - a microfluidic device to actively counteract
Brownian motion of single nanoparticles with an electrokinetic trap. Here we
present a version of an ABELtrap with a laser focus pattern generated by
electro-optical beam deflectors and controlled by a programmable FPGA chip.
This ABELtrap holds single fluorescent nanoparticles for more than 100 seconds,
increasing the observation time of fluorescent nanoparticles compared to free
diffusion by a factor of 10000. To monitor conformational changes of individual
membrane proteins in real time, we record sequential distance changes between
two specifically attached dyes using F\"orster resonance energy transfer
(smFRET). Fusing the a-subunit of the FoF1-ATP synthase with mNeonGreen results
in an improved signal-to-background ratio at lower laser excitation powers.
This increases our measured trap duration of proteoliposomes beyond 2 s.
Additionally, we observe different smFRET levels attributed to varying
distances between the FRET donor (mNeonGreen) and acceptor (Alexa568)
fluorophore attached at the a- and c-subunit of the FoF1-ATP synthase
respectively.


  We discuss how information encoded in a template polymer can be
stochastically copied into a copy polymer. We consider four different
stochastic copy protocols of increasing complexity, inspired by building blocks
of the mRNA translation pathway. In the first protocol, monomer incorporation
occurs in a single stochastic transition. We then move to a more elaborate
protocol in which an intermediate step can be used for error correction.
Finally, we discuss the operating regimes of two kinetic proofreading
protocols: one in which proofreading acts from the final copying step, and one
in which it acts from an intermediate step. We review known results for these
models and, in some cases, extend them to analyze all possible combinations of
energetic and kinetic discrimination. We show that, in each of these protocols,
only a limited number of these combinations leads to an improvement of the
overall copying accuracy.


  We study an individual-based model in which two spatially-distributed
species, characterized by different diffusivities, compete for resources. We
consider three different ecological settings. In the first, diffusing faster
has a cost in terms of reproduction rate. In the second case, resources are not
uniformly distributed in space. In the third case, the two species are
transported by a fluid flow. In all these cases, at varying the parameters, we
observe a transition from a regime in which diffusing faster confers an
effective selective advantage to one in which it constitutes a disadvantage. We
analytically estimate the magnitude of this advantage (or disadvantage) and
test it by measuring fixation probabilities in simulations of the
individual-based model. Our results provide a framework to quantify
evolutionary pressure for increased or decreased dispersal in a given
environment.


  Species tree reconstruction is complicated by effects of Incomplete Lineage
Sorting (ILS), commonly modeled by the multi-species coalescent model. While
there has been substantial progress in developing methods that estimate a
species tree given a collection of gene trees, less attention has been paid to
fast and accurate methods of quantifying support. In this paper, we propose a
fast algorithm to compute quartet-based support for each branch of a given
species tree with regard to a given set of gene trees. We then show how the
quartet support can be used in the context of the multi-species coalescent
model to compute i) the local posterior probability that the branch is in the
species tree and ii) the length of the branch in coalescent units. We evaluate
the precision and recall of the local posterior probability on a wide set of
simulated and biological data, and show that it has very high precision and
improved recall compared to multi-locus bootstrapping. The estimated branch
lengths are highly accurate when gene trees have little error, but are
underestimated when gene tree estimation error increases. Computation of both
branch length and local posterior probability is implemented as a new feature
in ASTRAL.


  Runs of homozygosity (ROH) are contiguous stretches of homozygous genome
which likely reflect transmission from common ances- tors and can be used to
track the inheritance of haplotypes of interest. In the present paper, ROH were
extracted from 50K SNPs and used to detect regions of the genome associated
with susceptibility to diseases in a population of 468 Holstein-Frisian cows.
Diagnosed diseases were categorised as infectious diseases, metabolic
syndromes, mastitis, reproductive diseases and locomotive disorders. ROH
associated with infectious diseases, mastitis and locomotive disorders were
found on BTA 12. A long region of homozygosity linked with metabolic syndromes,
infectious and reproductive diseases was detected on BTA 15, disclosing complex
relationships between immunity, metabolism and functional disorders. ROH
associated with infectious and reproductive diseases, mastitis and metabolic
syndromes were observed on chromosomes 3, 5, 7, 13 and 18. Previous studies
reported QTLs for milk production traits on all of these regions, thus
substantiating the known negative relationship between selection for milk
production and health in dairy cattle.


  The prevalent technique for DNA sequencing consists of two main steps:
shotgun sequencing, where many randomly located fragments, called reads, are
extracted from the overall sequence, followed by an assembly algorithm that
aims to reconstruct the original sequence. There are many different
technologies that generate the reads: widely-used second-generation methods
create short reads with low error rates, while emerging third-generation
methods create long reads with high error rates. Both error rates and error
profiles differ among methods, so reconstruction algorithms are often tailored
to specific shotgun sequencing technologies. As these methods change over time,
a fundamental question is whether there exist reconstruction algorithms which
are robust, i.e., which perform well under a wide range of error distributions.
  Here we study this question of sequence assembly from corrupted reads. We
make no assumption on the types of errors in the reads, but only assume a bound
on their magnitude. More precisely, for each read we assume that instead of
receiving the true read with no errors, we receive a corrupted read which has
edit distance at most $\epsilon$ times the length of the read from the true
read. We show that if the reads are long enough and there are sufficiently many
of them, then approximate reconstruction is possible: we construct a simple
algorithm such that for almost all original sequences the output of the
algorithm is a sequence whose edit distance from the original one is at most
$O(\epsilon)$ times the length of the original sequence.