  Many real-world complex systems are well represented as multilayer networks;
predicting interactions in those systems is one of the most pressing problems
in predictive network science. To address this challenge, we introduce two
stochastic block models for multilayer and temporal networks; one of them uses
nodes as its fundamental unit, whereas the other focuses on links. We also
develop scalable algorithms for inferring the parameters of these models.
Because our models describe all layers simultaneously, our approach takes full
advantage of the information contained in the whole network when making
predictions about any particular layer. We illustrate the potential of our
approach by analyzing two empirical datasets---a temporal network of email
communications, and a network of drug interactions for treating different
cancer types. We find that modeling all layers simultaneously does result, in
general, in more accurate link prediction. However, the most predictive model
depends on the dataset under consideration; whereas the node-based model is
more appropriate for predicting drug interactions, the link-based model is more
appropriate for predicting email communication.


  Determining the relative importance of environmental factors, biotic
interactions and stochasticity in assembling and maintaining species-rich
communities remains a major challenge in ecology. In plant communities,
interactions between individuals of different species are expected to leave a
spatial signature in the form of positive or negative spatial correlations over
distances relating to the spatial scale of interaction. Most studies using
spatial point process tools have found relatively little evidence for
interactions between pairs of species. More interactions tend to be detected in
communities with fewer species. However, there is currently no understanding of
how the power to detect spatial interactions may change with sample size, or
the scale and intensity of interactions.
  We use a simple 2-species model where the scale and intensity of interactions
are controlled to simulate point pattern data. In combination with an
approximation to the variance of the spatial summary statistics that we sample,
we investigate the power of current spatial point pattern methods to correctly
reject the null model of bivariate species independence.
  We show that the power to detect interactions is positively related to the
abundances of the species tested, and the intensity and scale of interactions.
Increasing imbalance in abundances has a negative effect on the power to detect
interactions. At population sizes typically found in currently available
datasets for species-rich plant communities we find only a very low power to
detect interactions. Differences in power may explain the increased frequency
of interactions in communities with fewer species. Furthermore, the
community-wide frequency of detected interactions is very sensitive to a
minimum abundance criterion for including species in the analyses.


  We define a new balance index for rooted phylogenetic trees based on the
symmetry of the evolutive history of every set of 4 leaves. This index makes
sense for multifurcating trees and it can be computed in time linear in the
number of leaves. We determine its maximum and minimum values for arbitrary and
bifurcating trees, and we provide exact formulas for its expected value and
variance on bifurcating trees under Ford's $\alpha$-model and Aldous'
$\beta$-model and on arbitrary trees under the $\alpha$-$\gamma$-model.


  Unequal partitioning of the molecular content at cell division has been shown
to be a source of heterogeneity in a cell population. We propose to model this
phenomenon with the help of a scalar, nonlinear impulsive differential equation
(IDE). In a first part, we consider a general autonomous IDE with fixed times
of impulse and a specific form of impulse function. We establish properties of
the solutions of that equation, most of them obtained under the hypothesis that
impulses occur periodically. In particular, we show how to investigate the
existence of periodic solutions and their stability by studying the flow of an
autonomous differential equation. A second part is dedicated to the analysis of
the convexity of this flow. Finally, we apply those results to an IDE
describing the concentration of the protein Tbet in a CD8 T-cell, where
impulses are associated to cell division, to study the effect of molecular
partitioning at cell division on the effector/memory cell-fate decision in a
CD8 T-cell lineage. We show that the degree of asymmetry in the molecular
partitioning can affect the process of cell differentiation and the
phenotypical heterogeneity of a cell population.


  Temporal evolution of a clonal bacterial population is modelled taking into
account reversible mutation and selection mechanisms. For the mutation model,
an efficient algorithm is proposed to verify whether experimental data can be
explained by this model. The selection-mutation model has unobservable fitness
parameters and, to estimate them, we use an Approximate Bayesian Computation
(ABC) algorithm. The algorithms are illustrated using in vitro data for phase
variable genes of Campylobacter jejuni.


  Forest-fire and avalanche models support the notion that frequent
catastrophes prevent the growth of very large populations and as such prevent
rare large-scale catastrophes. We show that this notion is not universal. A new
model class leads to a paradigm shift in the influence of catastrophes on the
family-size distribution of sub-populations. We study a simple population
dynamics model where individuals, as well as a whole family, may die with a
constant probability, accompanied by a logistic population growth model. We
compute the characteristics of the family-size distribution in steady-state and
the phase diagram of the steady-state distribution, and show that the family
and catastrophe size variances increase with the catastrophe frequency, which
is the opposite of common intuition. Frequent catastrophes are balanced by a
larger net-growth rate in surviving families, leading to the exponential growth
of these families. When the catastrophe rate is further increased, a second
phase transition to extinction occurs, when the rate of new families creations
is lower than their destruction rate by catastrophes.


  The activation of matrix metalloproteinase 2 (MMP2) is a crucial event during
tumor metastasis and invasion, and this pathway network consists of 3 monomers.
The pathway network of the activation obeys to a set of specified reaction
rules. According to the rules, the individual molecules localize in a
particular order and symmetrically around a homodimer following the formation
of that dimer. We generalized the homodimer pathway network obeying to similar
reaction rules, by changing the number of monomers involved in this pathway
from 3 to N. At the previous work, we found the molecules in the pathway
network are classified to some reaction groups. We derived the law of mass
conservation between the groups. Each group concentration converges to its
equilibrium solution. Using these results, we derive the concentrations of the
complexes theoretically and reveal that each complex concentration converges to
its equilibrium value. We can say the pathway network with homodimer symmetric
form complexes is asymptotic stable and identify the regulator parameter of the
target complex in the network. Our mathematical approach may help us understand
the mechanism of this type pathway network by knowing the background
mathematical laws which govern this type pathway network.


  Double-stranded DNA may contain mismatched base pairs beyond the Watson-Crick
pairs guanine-cytosine and adenine-thymine. Such mismatches bear adverse
consequences for human health. We utilize molecular dynamics and metadynamics
computer simulations to study the equilibrium structure and dynamics for both
matched and mismatched base pairs. We discover significant differences between
matched and mismatched pairs in structure, hydrogen bonding, and base flip work
profiles. Mismatched pairs shift further in the plane normal to the DNA strand
and are more likely to exhibit non-canonical structures, including the e-motif.
We discuss potential implications on mismatch repair enzymes' detection of DNA
mismatches.


  This paper builds up on T. Martinez-Coronado, A. Mir, F. Rossello and G.
Valiente's work "A balance index for phylogenetic trees based on quartets",
introducing a new balance index for trees. We show here that this balance
index, in the case of Aldous's b>=0-model, convergences weakly to a
distribution that can be characterized as the fixed point of a contraction
operator on a class of distributions.


  The Malagasy language belongs to the Greater Barito East group of the
Austronesian family, the language most closely connected to Malagasy dialects
is Maanyan (Kalimantan), but Malay as well other Indonesian and Philippine
languages are also related. The African contribution is very high in the
Malagasy genetic make-up (about 50%) but negligible in the language. Because of
the linguistic link, it is widely accepted that the island was settled by
Indonesian sailors after a maritime trek but date and place of landing are
still debated. The 50% Indonesian genetic contribution to present Malagasy
points in a different direction then Maanyan for the Asian ancestry, therefore,
the ethnic composition of the Austronesian settlers is also still debated. In
this talk I mainly review the joint research of Filippo Petroni, Dima
Volchenkov, S\"oren Wichmann and myself which tries to shed new light on these
problems. The key point is the application of a new quantitative methodology
which is able to find out the kinship relations among languages (or dialects).
New techniques are also introduced in order to extract the maximum information
from these relations concerning time and space patterns.


  The analysis of a protein-expression pattern from tissue microarray (TMA)
data will not immediately give an answer on synergistic or antagonistic effects
between the expression of the observed proteins. But contrary to apparent first
impression, it is possible to reveal those cooperative phenomena from TMA data.
We present here a largely assumption-free combinatorial analysis, related to
correlation networks but with much less arbitrary constraints. A strong focus
was put on the analysis of the basic data to analyze how the cooperative
phenomena might be imprinted in the TMA data structure. The study design was
based on two independent panels of 589 and 366 invasive breast cancer cases
from different institutions, assembled on tissue microarrays. The combinatorial
analysis generates an optimal rank ordering of protein-expression coherence.
The outcome of the analysis corresponds to all the single observations
scattered over several publications and integrates them in one context. This
means all these scattered observations can also be deduced from one TMA
experiment. A comprehensive statistical meta-analysis of the TMA data suggests
the existence of a superposition of three basic coherence situations, and
offers the opportunity to analyze these data properties with additional
real-world data and synthetic data in more detail. The presented algorithm
gives molecular pathologists a tool to extract dependency information from TMA
data. Beyond this practical benefit, some light was shed on how dependency
aspects might be imprinted into expression data. This will certainly foster the
refinement of algorithms to reconstruct dependency networks. The implementation
of the algorithm is at the moment not end-user suitable, but available on
request.


  Protein-ligand scoring is an important step in a structure-based drug design
pipeline. Selecting a correct binding pose and predicting the binding affinity
of a protein-ligand complex enables effective virtual screening. Machine
learning techniques can make use of the increasing amounts of structural data
that are becoming publicly available. Convolutional neural network (CNN)
scoring functions in particular have shown promise in pose selection and
affinity prediction for protein-ligand complexes. Neural networks are known for
being difficult to interpret. Understanding the decisions of a particular
network can help tune parameters and training data to maximize performance.
Visualization of neural networks helps decompose complex scoring functions into
pictures that are more easily parsed by humans. Here we present three methods
for visualizing how individual protein-ligand complexes are interpreted by 3D
convolutional neural networks. We also present a visualization of the
convolutional filters and their weights. We describe how the intuition provided
by these visualizations aids in network design.


  Next-generation sequencing (NGS) to profile temporal changes in living
systems is gaining more attention for deriving better insights into the
underlying biological mechanisms compared to traditional static sequencing
experiments. Nonetheless, the majority of existing statistical tools for
analyzing NGS data lack the capability of exploiting the richer information
embedded in temporal data. Several recent tools have been developed to analyze
such data but they typically impose strict model assumptions, such as
smoothness on gene expression dynamic changes. To capture a broader range of
gene expression dynamic patterns, we develop the gamma Markov negative binomial
(GMNB) model that integrates a gamma Markov chain into a negative binomial
distribution model, allowing flexible temporal variation in NGS count data.
Using Bayes factors, GMNB enables more powerful temporal gene differential
expression analysis across different phenotypes or treatment conditions. In
addition, it naturally handles the heterogeneity of sequencing depth in
different samples, removing the need for ad-hoc normalization. Efficient Gibbs
sampling inference of the GMNB model parameters is achieved by exploiting novel
data augmentation techniques. Extensive experiments on both simulated and
real-world RNA-seq data show that GMNB outperforms existing methods in both
receiver operating characteristic (ROC) and precision-recall (PR) curves of
differential expression analysis results.


  When three cyanobacterial proteins, KaiA, KaiB, and KaiC, are incubated with
ATP in vitro, the phosphorylation level of KaiC hexamers shows stable
oscillation with approximately 24 h period. In order to understand this KaiABC
clockwork, we need to analyze both the macroscopic synchronization of a large
number of KaiC hexamers and the microscopic reactions and structural changes in
individual KaiC molecules. In the present paper, we explain two coarse-grained
theoretical models, the many-molecule (MM) model and the single-molecule (SM)
model, to bridge the gap between macroscopic and microscopic understandings. In
the simulation results with these models, ATP hydrolysis drives oscillation of
individual KaiC hexamers and ATP hydrolysis is necessary for synchronizing
oscillations of a large number of KaiC hexamers. Sensitive temperature
dependence of the lifetime of the ADP bound state in the CI domain of KaiC
hexamers makes the oscillation period temperature insensitive. ATPase activity
is correlated to the frequency of phosphorylation oscillation in the single
molecule of KaiC hexamer, which should be the origin of the observed
ensemble-level correlation between the ATPase activity and the frequency of
phosphorylation oscillation. Thus, the simulation results with the MM and SM
models suggest that ATP hydrolysis randomly occurring in each CI domain of
individual KaiC hexamers is a key process for oscillatory behaviors of the
ensemble of many KaiC hexamers.


  We revisit the well-known chemostat model, considering that bacteria can be
attached together in aggregates or flocs. We distinguish explicitly free and
attached compartments in the model and give sufficient conditions for
coexistence of these two forms. We then study the case of fast attachment and
detachment and shows how it is related to density-dependent growth functions.
Finally, we give some insights concerning the cases of multi-specific flocs and
different removal rates.


  Radicalization is the process by which people come to adopt increasingly
extreme political, social or religious ideologies. When radicalization leads to
violence, radical thinking becomes a threat to national security.
De-radicalization programs are part of an effort to combat violent extremism
and terrorism. This type of initiatives attempt to alter violent extremists
radical beliefs and violent behavior with the aim to reintegrate them into
society. In this paper we introduce a simple compartmental model suitable to
describe de-radicalization programs. The population is divided into four
compartments: $ (S) $ susceptible, $ (E) $ extremists, $ (R) $ recruiters, and
$ (T) $ treatment. We calculate the basic reproduction number $ \mathcal{R}_0
$. For $ \mathcal{R}_0< 1 $ the system has one globally asymptotically stable
equilibrium where no extremist or recruiters are present. For $ \mathcal{R}_0
>1 $ the system has an additional equilibrium where extremists and recruiters
are endemic to the population. A Lyapunov function is used to show that, for $
\mathcal{R}_0 >1 $, the endemic equilibrium is globally asymptotically stable.
We use numerical simulations to support our analytical results. Based on our
model we asses strategies to counter violent extremism.


  With advances in sequencing technologies, there are now massive amounts of
genomic data from across all life, leading to the possibility that a robust
Tree of Life can be constructed. However, "gene tree heterogeneity", which is
when different genomic regions can evolve differently, is a common phenomenon
in multi-locus datasets, and reduces the accuracy of standard methods for
species tree estimation that do not take this heterogeneity into account. New
methods have been developed for species tree estimation that specifically
address gene tree heterogeneity, and that have been proven to converge to the
true species tree when the number of loci and number of sites per locus both
increase (i.e., the methods are said to be "statistically consistent"). Yet,
little is known about the biologically realistic condition where the number of
sites per locus is bounded. We show that when the sequence length of each locus
is bounded (by any arbitrarily chosen value), the most common approaches to
species tree estimation that take heterogeneity into account (i.e., traditional
fully partitioned concatenated maximum likelihood and newer approaches, called
summary methods, that estimate the species tree by combining gene trees) are
not statistically consistent, even when the heterogeneity is extremely
constrained. The main challenge is the presence of conditions such as long
branch attraction that create biased tree estimation when the number of sites
is restricted. Hence, our study uncovers a fundamental challenge to species
tree estimation using both traditional and new methods.


  Evolution is based on the assumption that competing players update their
strategies to increase their individual payoffs. However, while the applied
updating method can be different, most of previous works proposed uniform
models where players use identical way to revise their strategies. In this work
we explore how imitation-based or learning attitude and innovation-based or
myopic best response attitude compete for space in a complex model where both
attitudes are available. In the absence of additional cost the best response
trait practically dominates the whole snow-drift game parameter space which is
in agreement with the average payoff difference of basic models. When
additional cost is involved then the imitation attitude can gradually invade
the whole parameter space but this transition happens in a highly nontrivial
way. However, the role of competing attitudes is reversed in the stag-hunt
parameter space where imitation is more successful in general. Interestingly, a
four-state solution can be observed for the latter game which is a consequence
of an emerging cyclic dominance between possible states. These phenomena can be
understood by analyzing the microscopic invasion processes, which reveals the
unequal propagation velocities of strategies and attitudes.


  We study classical stochastic systems with discrete states, coupled to
switching external environments. For fast environmental processes we derive
reduced dynamics for the system itself, focusing on corrections to the
adiabatic limit of infinite time scale separation. In some cases, this leads to
master equations with negative transition `rates' or bursting events. We devise
a simulation algorithm in discrete time to unravel these master equations into
sample paths, and provide an interpretation of bursting events. Focusing on
stochastic population dynamics coupled to external environments, we discuss a
series of approximation schemes combining expansions in the inverse switching
rate of the environment, and a Kramers--Moyal expansion in the inverse size of
the population. This places the different approximations in relation to
existing work on piecewise-deterministic and piecewise-diffusive Markov
processes. We apply the model reduction methods to different examples including
systems in biology and a model of crack propagation.


  A systematic understanding of the evolution and growth dynamics of invasive
solid tumors in response to different chemotherapy strategies is crucial for
the development of individually optimized oncotherapy. Here, we develop a
hybrid three-dimensional (3D) computational model that integrates
pharmacokinetic model, continuum diffusion-reaction model and discrete cell
automaton model to investigate 3D invasive solid tumor growth in heterogeneous
microenvironment under chemotherapy. Specifically, we consider the effects of
heterogeneous environment on drug diffusion, tumor growth, invasion and the
drug-tumor interaction on individual cell level. We employ the hybrid model to
investigate the evolution and growth dynamics of avascular invasive solid
tumors under different chemotherapy strategies. Our simulations reproduce the
well-established observation that constant dosing is generally more effective
in suppressing primary tumor growth than periodic dosing, due to the resulting
continuous high drug concentration. In highly heterogeneous microenvironment,
the malignancy of the tumor is significantly enhanced, leading to inefficiency
of chemotherapies. The effects of geometrically-confined microenvironment and
non-uniform drug dosing are also investigated. Our computational model, when
supplemented with sufficient clinical data, could eventually lead to the
development of efficient in silico tools for prognosis and treatment strategy
optimization.


  Genome-wide association studies (GWAS) have identified single nucleotide
polymorphisms (SNPs) associated with trait diversity and disease
susceptibility, yet the functional properties of many genetic variants and
their molecular interactions remains unclear. It has been hypothesized that
SNPs in microRNA binding sites may disrupt gene regulation by microRNAs
(miRNAs), short non-coding RNAs that bind to mRNA and downregulate the target
gene. While a number of studies have been conducted to predict the location of
SNPs in miRNA binding sites, to date there has been no comprehensive analysis
of how SNP variants may impact miRNA regulation of genes. Here we investigate
the functional properties of genetic variants and their effects on miRNA
regulation of gene expression in cancer. Our analysis is motivated by the
hypothesis that distinct alleles may cause differential binding (from miRNAs to
mRNAs or from transcription factors to DNA) and change the expression of genes.
We previously identified pathways--systems of genes conferring specific cell
functions--that are dysregulated by miRNAs in cancer, by comparing
miRNA-pathway associations between healthy and tumor tissue. We draw on these
results as a starting point to assess whether SNPs in genes on dysregulated
pathways are responsible for miRNA dysregulation of individual genes in tumors.
Using an integrative analysis that incorporates miRNA expression, mRNA
expression, and SNP genotype data, we identify SNPs that appear to influence
the association between miRNAs and genes, which we term "regulatory QTLs
(regQTLs)": loci whose alleles impact the regulation of genes by miRNAs. We
describe the method, apply it to analyze four cancer types (breast, liver,
lung, prostate) using data from The Cancer Genome Atlas (TCGA), and provide a
tool to explore the findings.


  The evolution of antimicrobial resistance generally occurs in an environment
where antimicrobial concentration is variable, which has dramatic consequences
on the microorganisms' fitness landscape, and thus on the evolution of
resistance. We investigate the effect of these time-varying patterns of
selection within a stochastic model. We consider a homogeneous microbial
population of fixed size subjected to periodic alternations of phases of
absence and presence of an antimicrobial that stops growth. Combining
analytical approaches and stochastic simulations, we quantify how the time
necessary for fit resistant bacteria to take over the microbial population
depends on the alternation period. We demonstrate that fast alternations
strongly accelerate the evolution of resistance, reaching a plateau for
sufficiently small periods. Furthermore, this acceleration is stronger in
larger populations. For asymmetric alternations, featuring a different duration
of the phases with and without antimicrobial, we shed light on the existence of
a minimum for the time taken by the population to fully evolve resistance. The
corresponding dramatic acceleration of the evolution of antimicrobial
resistance likely occurs in realistic situations, and may have an important
impact both in clinical and experimental situations.


  Owing to its longevity and enormous information density, DNA, the molecule
encoding biological information, has emerged as a promising archival storage
medium. However, due to technological constraints, data can only be written
onto many short DNA molecules that are stored in an unordered way, and can only
be read by sampling from this DNA pool. Moreover, imperfections in writing
(synthesis), reading (sequencing), storage, and handling of the DNA, in
particular amplification via PCR, lead to a loss of DNA molecules and induce
errors within the molecules. In order to design DNA storage systems, a
qualitative and quantitative understanding of the errors and the loss of
molecules is crucial. In this paper, we characterize those error probabilities
by analyzing data from our own experiments as well as from experiments of two
different groups. We find that errors within molecules are mainly due to
synthesis and sequencing, while imperfections in handling and storage lead to a
significant loss of sequences. The aim of our study is to help guide the design
of future DNA data storage systems by providing a quantitative and qualitative
understanding of the DNA data storage channel.


  The predominant structural protein in vertebrates is collagen, which plays a
key role in extracellular matrix and connective tissue mechanics. Despite its
prevalence and physical importance in biology, the mechanical properties of
molecular collagen are far from established. The flexibility of its triple
helix is unresolved, with descriptions from different experimental techniques
ranging from flexible to semirigid. Furthermore, it is unknown how collagen
type (homo- vs. heterotrimeric) and source (tissue-derived vs. recombinant)
influence flexibility. Using SmarTrace, a chain tracing algorithm we devised,
we performed statistical analysis of collagen conformations collected with
atomic force microscopy (AFM) to determine the protein's mechanical properties.
Our results show that types I, II and III collagens - the key fibrillar
varieties - exhibit molecular flexibilities that are very similar. However,
collagen conformations are strongly modulated by salt, transitioning from
compact to extended as KCl concentration increases, in both neutral and acidic
pH. While analysis with a standard worm-like chain model suggests that the
persistence length of collagen can attain almost any value within the
literature range, closer inspection reveals that this modulation of collagen's
conformational behavior is not due to changes in flexibility, but rather arises
from the induction of curvature (either intrinsic or induced by interactions
with the mica surface). By modifying standard polymer theory to include innate
curvature, we show that collagen behaves as an equilibrated curved worm-like
chain (cWLC) in two dimensions. Analysis within the cWLC model shows that
collagen's curvature depends strongly on pH and salt, while its persistence
length does not. Thus, we find that triple-helical collagen is well described
as semiflexible, irrespective of source, type, pH and salt environment.


  The genetic regulatory network (GRN) plays a key role in controlling the
response of the cell to changes in the environment. Although the structure of
GRNs has been the subject of many studies, their large scale structure in the
light of feedbacks from the metabolic network (MN) has received relatively
little attention. Here we study the causal structure of the GRNs, namely the
chain of influence of one component on the other, taking into account feedback
from the MN. First we consider the GRNs of E. coli and B. subtilis without
feedback from MN and illustrate their causal structure. Next we augment the
GRNs with feedback from their respective MNs by including (a) links from genes
coding for enzymes to metabolites produced or consumed in reactions catalyzed
by those enzymes and (b) links from metabolites to genes coding for
transcription factors whose transcriptional activity the metabolites alter by
binding to them. We find that the inclusion of feedback from MN into GRN
significantly affects its causal structure, in particular the number of levels
and relative positions of nodes in the hierarchy, and the number and size of
the strongly connected components (SCCs). We then study the functional
significance of the SCCs. For this we identify condition specific feedbacks
from the MN into the GRN by retaining only those enzymes that are essential for
growth in specific environmental conditions simulated via the technique of flux
balance analysis (FBA). We find that the SCCs of the GRN augmented by these
feedbacks can be ascribed specific functional roles in the organism. Our
algorithmic approach thus reveals relatively autonomous subsystems with
specific functionality, or regulatory modules in the organism. This automated
approach could be useful in identifying biologically relevant modules in other
organisms for which network data is available, but whose biology is less well
studied.


  As impacts of introduced species cascade through trophic levels, they can
cause indirect and counter-intuitive effects. To investigate the impact of
invasive species at the network scale, we use a generalized food web model,
capable of propagating changes through networks with a series of ecologically
realistic criteria. Using data from a small British offshore island, we
quantify the impacts of four virtual invasive species (an insectivore, a
herbivore, a carnivore and an omnivore whose diet is based on a rat) and
explore which clusters of species react in similar ways. We find that the
predictions for the impacts of invasive species are ecologically plausible,
even for large networks robust predictions for the impacts of invasive species
can be obtained. Species in the same taxonomic group are similarly impacted by
a virtual invasive species. However, interesting differences within a given
taxonomic group can occur. The results suggest that some native species may be
at risk from a wider range of invasives than previously believed. The
implications of these results for ecologists and land managers are discussed.


  This paper studies the dynamics of a network-based SIRS epidemic model with
vaccination and a nonmonotone incidence rate. This type of nonlinear incidence
can be used to describe the psychological or inhibitory effect from the
behavioral change of the susceptible individuals when the number of infective
individuals on heterogeneous networks is getting larger. Using the analytical
method, epidemic threshold $R_0$ is obtained. When $R_0$ is less than one, we
prove the disease-free equilibrium is globally asymptotically stable and the
disease dies out, while $R_0$ is greater than one, there exists a unique
endemic equilibrium. By constructing a suitable Lyapunov function, we also
prove the endemic equilibrium is globally asymptotically stable if the
inhibitory factor $\alpha$ is sufficiently large. Numerical experiments are
also given to support the theoretical results. It is shown both theoretically
and numerically a larger $\alpha$ can accelerate the extinction of the disease
and reduce the level of disease.


  Investigations of a possible connection between population density and the
propagation and magnitude of epidemics have so far led to mixed and
unconvincing results. There are three reasons for that. (i) Previous studies
did not focus on the appropriate density interval. (ii) For the density to be a
meaningful variable the population must be distributed as uniformly as
possible. If an area has towns and cities where a majority of the population is
concentrated its average density is meaningless. (iii) In the propagation of an
epidemic the initial proportion of susceptibles (that is to say persons who
have not developed an immunity) is an essential, yet usually unknown, factor.
The assumption that most of the population is susceptible holds only for new
strain of diseases.
  It will be shown that when these requirements are taken care of, the size of
epidemics is indeed closely connected with the population density. This
empirical observation comes as a welcome confirmation of the classical KMK
(Kermack-McKendrick 1927) model. Indeed, one of its key predictions is that the
size of the epidemic increases strongly (and in a non linear way) with the
initial density of susceptibles.
  An interesting consequence is that, contrary to common beliefs, in sparsely
populated territories, like Alaska, Australia or the west coast of the United
states the size of epidemics among native populations must have been limited by
the low density even for diseases for which the natives had no immunity (i.e.,
were susceptibles).


  The totally asymmetric simple exclusion process (TASEP) is a paradigmatic
stochastic model for non-equilibrium physics, and has been successfully applied
to describe active transport of molecular motors along cytoskeletal filaments.
Building on this simple model, we consider a two-lane lattice-gas model that
couples directed transport (TASEP) to diffusive motion in a semi-closed
geometry, and simultaneously accounts for spontaneous growth and
particle-induced shrinkage of the system's size. This particular extension of
the TASEP is motivated by the question of how active transport and diffusion
might influence length regulation in confined systems. Surprisingly, we find
that the size of our intrinsically stochastic system exhibits robust temporal
patterns over a broad range of growth rates. More specifically, when particle
diffusion is slow relative to the shrinkage dynamics, we observe quasi-periodic
changes in length. We provide an intuitive explanation for the occurrence of
these self-organized temporal patterns, which is based on the imbalance between
the diffusion and shrinkage speed in the confined geometry. Finally, we
formulate an effective theory for the oscillatory regime, which explains the
origin of the oscillations and correctly predicts the dependence of key
quantities, as for instance the oscillation frequency, on the growth rate.


  In this paper, we affirm our earlier findings of evidence for a limit to
human lifespan. In particular, we assess the analyses in extreme value theory
(EVT) performed by Rootz\'en and Zholud. We find that their criticisms of our
work are unfounded and that their analyses are contradicted by several other
papers using EVT. Furthermore, we find that even if we completely accept the
conclusions about late-life human mortality reached by Rootz\'en and Zholud,
their results do not actually contradict the findings presented in our original
paper: whether unbounded or not, human lifespan is unlikely to greatly exceed
120 years, and the improbability of longer survival---whether it is exactly
zero or merely astronomically small---acts as a de facto limit. In order to
eliminate the confusion surrounding the issue, we propose the adoption of the
term "limit" to denote the age at which the chance of survival is exactly zero
and the term "effective limit" to denote the age at which the change of
survival falls below a given threshold. Once this distinction is made, it can
be demonstrated that the final result of Rootz\'en and Zholud is essentially a
recapitulation of the main conclusion of our paper. Ultimately, much of the
controversy surrounding the issue of a limit to human lifespan can be avoided
by carefully reading the literature and applying statistics to practical human
scales.


  Spectropolarimetry of intact plant leaves allows to probe the molecular
architecture of vegetation photosynthesis in a non-invasive and non-destructive
way and, as such, can offer a wealth of physiological information. In addition
to the molecular signals due to the photosynthetic machinery, the cell
structure and its arrangement within a leaf can create and modify polarization
signals. Using Mueller matrix polarimetry with rotating retarder modulation, we
have visualized spatial variations in polarization in transmission around the
chlorophyll a absorbance band from 650 nm to 710 nm. We show linear and
circular polarization measurements of maple leaves and cultivated maize leaves
and discuss the corresponding Mueller matrices and the Mueller matrix
decompositions, which show distinct features in diattenuation, polarizance,
retardance and depolarization. Importantly, while normal leaf tissue shows a
typical split signal with both a negative and a positive peak in the induced
fractional circular polarization and circular dichroism, the signals close to
the veins only display a negative band. The results are similar to the negative
band as reported earlier for single macrodomains. We discuss the possible role
of the chloroplast orientation around the veins as a cause of this phenomenon.
Systematic artefacts are ruled out as three independent measurements by
different instruments gave similar results. These results provide better
insight into circular polarization measurements on whole leaves and options for
vegetation remote sensing using circular polarization.


  Bacterial communities have rich social lives. A well-established interaction
involves the exchange of a public good in Pseudomonas populations, where the
iron-scavenging compound pyoverdine, synthesized by some cells, is shared with
the rest. Pyoverdine thus mediates interactions between producers and
non-producers and can constitute a public good. This interaction is often used
to test game theoretical predictions on the "social dilemma" of producers. Such
an approach, however, underestimates the impact of specific properties of the
public good, for example consequences of its accumulation in the environment.
Here, we experimentally quantify costs and benefits of pyoverdine production in
a specific environment, and build a model of population dynamics that
explicitly accounts for the changing significance of accumulating pyoverdine as
chemical mediator of social interactions. The model predicts that, in an
ensemble of growing populations (metapopulation) with different initial
producer fractions (and consequently pyoverdine contents), the global producer
fraction initially increases. Because the benefit of pyoverdine declines at
saturating concentrations, the increase need only be transient. Confirmed by
experiments on metapopulations, our results show how a changing benefit of a
public good can shape social interactions in a bacterial population.


  In order to solve a system of nonlinear rate equations one can try to use
some soliton methods. The procedure involves three steps: (1) Find a `Lax
representation' where all the kinetic variables are combined into a single
matrix $\rho$, all the kinetic constants are encoded in a matrix $H$; (2) find
a Darboux-Backund dressing transformation for the Lax representation $i\dot
\rho=[H,f(\rho)]$, where $f$ models a time-dependent environment; (3) find a
class of seed solutions $\rho=\rho[0]$ that lead, via a nontrivial chain of
dressings $\rho[0]\to \rho[1]\to \rho[2]\to\dots$ to new solutions, difficult
to find by other methods. The latter step is not a trivial one since a
non-soliton method has to be employed to find an appropriate initial $\rho[0]$.
Procedures that lead to a correct $\rho[0]$ have been discussed in the
literature only for a limited class of $H$ and $f$. Here, we develop a
formalism that works for practically any $H$, and any explicitly time-dependent
$f$. As a result, we are able to find exact solutions to a system of equations
describing an arbitrary number of species interacting through (auto)catalytic
feedbacks, with general time dependent parameters characterizing the
nonlinearity. Explicit examples involve up to 42 interacting species.


  Recent progress in experimental techniques has enabled us to quantitatively
study stochastic and flexible behavior of biological systems. For example, gene
regulatory networks perform stochastic information processing and their
functionalities have been extensively studied. In gene regulatory networks,
there are specific subgraphs called network motifs that occur at frequencies
much higher than those found in randomized networks. Further understanding of
the designing principle of such networks is highly desirable. In a different
context, information thermodynamics has been developed as a theoretical
framework that generalizes non-equilibrium thermodynamics to stochastically
fluctuating systems with information. Here we systematically characterize gene
regulatory networks on the basis of information thermodynamics. We model
three-node gene regulatory patterns by a stochastic Boolean model, which
receive one or two input signals that carry external information. For the case
of a single input, we found that all the three-node patterns are classified
into four types by using information-thermodynamic quantities such as
dissipation and mutual information, and reveal to which type each network motif
belongs. Next, we consider the case where there are two inputs, and evaluate
the capacity of logical operation of the three-node patterns by using
tripartite mutual information, and argue the reason why patterns with fewer
edges are preferred in natural selection. This result might also explain the
difference of the occurrence frequencies among different types of
feedforward-loop network motifs.


  When plated onto substrates, cell morphology and even stem cell
differentiation are influenced by the stiffness of their environment. Stiffer
substrates give strongly spread (eventually polarized) cells with strong focal
adhesions, and stress fibers; very soft substrates give a less developed
cytoskeleton, and much lower cell spreading. The kinetics of this process of
cell spreading is studied extensively, and important universal relationships
are established on how the cell area grows with time. Here we study the
population dynamics of spreading cells, investigating the characteristic
processes involved in cell response to the substrate. We show that unlike the
individual cell morphology, this population dynamics does not depend on the
substrate stiffness. Instead, a strong activation temperature dependence is
observed. Different cell lines on different substrates all have long-time
statistics controlled by the thermal activation over a single energy barrier
dG=19 kcal/mol, while the early-time kinetics follows a power law $t^5$. This
implies that the rate of spreading depends on an internal process of
adhesion-mechanosensing complex assembly and activation: the operational
complex must have 5 component proteins, and the last process in the sequence
(which we believe is the activation of focal adhesion kinase) is controlled by
the binding energy dG.


  Computational chemists typically assay drug candidates by virtually screening
compounds against crystal structures of a protein despite the fact that some
targets, like the $\mu$ Opioid Receptor and other members of the GPCR family,
traverse many non-crystallographic states. We discover new conformational
states of $\mu OR$ with molecular dynamics simulation and then machine learn
ligand-structure relationships to predict opioid ligand function. These
artificial intelligence models identified a novel $\mu$ opioid chemotype.


  Clusters of differentiation (CD) are cell surface biomarkers that denote key
biological differences between cell types and disease state. CD-targeting
therapeutic monoclonal antibodies (mAB) afford rich trans-disease repositioning
opportunities. Within a compendium of systemic lupus erythematous (SLE)
patients, we applied the Integrated machine learning pipeline for aberrant
biomarker enrichment (i-mAB) to profile de novo gene expression features
affecting CD20, CD22 and CD30 gene aberrance. First, a novel relief-based
algorithm identified interdependent features(p=681) predicting
treatment-na\"ive SLE patients (balanced accuracy=0.822). We then compiled
CD-associated expression profiles using regularized logistic regression and
pathway enrichment analyses. On an independent general cell line model system
data, we replicated associations (in silico) of BCL7A(padj=1.69e-9) and
STRBP(padj=4.63e-8) with CD22; NCOA2(padj=7.00e-4), ATN1(padj=1.71e-2), and
HOXC4(padj=3.34e-2) with CD30; and PHOSPHO1, a phosphatase linked to bone
mineralization, with both CD22(padj=4.37e-2) and CD30(padj=7.40e-3). Utilizing
carefully aggregated secondary data and leveraging a priori hypotheses, i-mAB
fostered robust biomarker profiling among interdependent biological features.


  The human neurotensin 1 receptor (hNTS1) is a G protein-coupled receptor
involved in many physiological functions, including analgesia, hypothermia, and
hypotension. To gain a better understanding of which signaling pathways or
combination of pathways are linked to NTS1 activation and function, we
investigated the ability of activated hNTS1, which was stably expressed by
CHO-K1 cells, to directly engage G proteins, activate second messenger cascades
and recruit \b{eta}-arrestins. Using BRET-based biosensors, we found that
neurotensin (NT), NT(8-13) and neuromedin N (NN) activated the G{\alpha}q-,
G{\alpha}i1-, G{\alpha}oA-, and G{\alpha}13-protein signaling pathways as well
as the recruitment of \b{eta}-arrestins 1 and 2. Using pharmacological
inhibitors, we further demonstrated that all three ligands stimulated the
production of inositol phosphate and modulation of cAMP accumulation along with
ERK1/2 activation. Interestingly, despite the functional coupling to
G{\alpha}i1 and G{\alpha}oA, NT was found to produce higher levels of cAMP in
the presence of pertussis toxin, supporting that hNTS1 activation leads to cAMP
accumulation in a G{\alpha}s-dependent manner. Additionally, we demonstrated
that the full activation of ERK1/2 required signaling through both a
PTX-sensitive Gi/o-c-Src signaling pathway and PLCb-DAG-PKC-Raf-1- dependent
pathway downstream of Gq. Finally, the whole-cell integrated signatures
monitored by the cell-based surface plasmon resonance and changes in the
electrical impedance of a confluent cell monolayer led to identical phenotypic
responses between the three ligands. The characterization of the hNTS1-mediated
cellular signaling network will be helpful to accelerate the validation of
potential NTS1 biased ligands with an improved therapeutic/adverse effect
profile.


  The apelinergic system is an important player in the regulation of both
vascular tone and cardiovascular function, making this physiological system an
attractive target for drug development for hypertension, heart failure and
ischemic heart disease. Indeed, apelin exerts a positive inotropic effect in
humans whilst reducing peripheral vascular resistance. In this study, we
investigated the signaling pathways through which apelin exerts its hypotensive
action. We synthesized a series of apelin-13 analogs whereby the C-terminal
Phe13 residue was replaced by natural or unnatural amino acids. In HEK293 cells
expressing APJ, we evaluated the relative efficacy of these compounds to
activate G{\alpha}i1 and G{\alpha}oA G-proteins, recruit \b{eta}-arrestins 1
and 2 (\b{eta}arrs), and inhibit cAMP production. Calculating the transduction
ratio for each pathway allowed us to identify several analogs with distinct
signaling profiles. Furthermore, we found that these analogs delivered i.v. to
Sprague-Dawley rats exerted a wide range of hypotensive responses. Indeed, two
compounds lost their ability to lower blood pressure, while other analogs
significantly reduced blood pressure as apelin-13. Interestingly, analogs that
did not lower blood pressure were less effective at recruiting \b{eta}arrs.
Finally, using Spearman correlations, we established that the hypotensive
response was significantly correlated with \b{eta}arr recruitment but not with
G protein- dependent signaling. In conclusion, our results demonstrated that
the \b{eta}arr recruitment potency is involved in the hypotensive efficacy of
activated APJ.


  Pathogenic Gram-negative bacteria have developed resistance to antibiotics
due to their ability in creating an envelope on the outer layer of
lipooligosaccharides (LOS). The cationic phosphoethanolamine (PEA) decoration
of LOS lipid A is regulated by lipid A-PEA transferase (EptA) which may serve
as a prominent target for developing new antibiotics. The structural
characterization of Neisserial EptA has provided a structural basis to its
catalytic mechanisms and ligand recognition that are crucial for inhibitor
development. In this study, a combination of pharmacophore- and ligand-based
approach has been employed to explore novel potent EptA inhibitors among
millions of commercially-available compounds and approved drugs. A total of
8166 hit molecules obtained from ZincPharmer pharmacophore-based screening and
PubMed ligand similarity search were further examined through individual
two-step semi-flexible docking simulation performed in MOE. Best hits were
therefore selected based on their docking score and consensus of the two
docking validations. Free energy of binding calculation suggests that the best
20 consensus compounds have a stronger binding affinity than EptA natural
substrate PEA. Further interaction analyses of selected eight ligands
demonstrate that these ligands have overall more effective interactions with
catalytically-essential residues and metal cofactors of EptA. Selected hits can
be further analyzed in vitro and examined through a pre-clinical trial. This
study provides an insight into drug repurposing which may serve as an initial
step to develop novel potent EptA inhibitors to combat the virulence of
multi-drug resistant Gram-negative bacteria.


  Logical models offer a simple but powerful means to understand the complex
dynamics of biochemical regulation, without the need to estimate kinetic
parameters. However, even simple automata components can lead to collective
dynamics that are computationally intractable when aggregated into networks. In
previous work we demonstrated that automata network models of biochemical
regulation are highly canalizing, whereby many variable states and their
groupings are redundant (Marques-Pita and Rocha, 2013). The precise charting
and measurement of such canalization simplifies these models, making even very
large networks amenable to analysis. Moreover, canalization plays an important
role in the control, robustness, modularity and criticality of Boolean network
dynamics, especially those used to model biochemical regulation (Gates and
Rocha, 2016; Gates et al., 2016; Manicka, 2017). Here we describe a new
publicly-available Python package that provides the necessary tools to extract,
measure, and visualize canalizing redundancy present in Boolean network models.
It extracts the pathways most effective in controlling dynamics in these
models, including their effective graph and dynamics canalizing map, as well as
other tools to uncover minimum sets of control variables.


  Probing the reaction mechanisms of supramolecular processes in soft- and
biological matter, such as protein aggregation, is inherently challenging.
These processes emerge from the simultaneous action of multiple molecular
mechanisms, each of which is associated with the rearrangement of a large
number of weak bonds, resulting in a complex free energy landscape with many
kinetic barriers. Reaction rate measurements of supramolecular processes at
different temperatures can offer unprecedented insights into the underlying
molecular mechanisms and their thermodynamic properties. However, to be able to
interpret such measurements in terms of the underlying microscopic mechanisms,
a key challenge is to establish which properties of the complex free energy
landscapes are probed by the reaction rate. Here, we present a reaction rate
theory for supramolecular kinetics based on Kramers rate theory for diffusive
reactions over multiple kinetic barriers, and apply the results to protein
aggregation. Using this framework and Monte Carlo simulations, we show that
reaction rates for protein aggregation are of the Arrhenius-Eyring type and
that the associated activation energies probe only one relevant barrier along
the respective free energy landscapes. We apply this advancement to interpret,
both in experiments and in coarse-grained computer simulations, reaction rate
measurements of amyloid aggregation kinetics in terms of the underlying
molecular mechanisms and associated thermodynamic signatures. Our results
establish a general platform for probing the mechanisms and energetics of
supramolecular phenomena in soft- and biological matter using the framework of
chemical kinetics.


  Constraint-based stoichiometric models are ubiquitous in metabolic research,
with Flux Balance Analysis (FBA) being the most widely used method to describe
metabolic phenotypes of cells growing in steady-state. Of the many variants of
constrain-based modelling methods published throughout the years, only few have
focused on thermodynamic issues, in particular the elimination of non-physical
and non-physiological cyclic fluxes. In this work, we revisit two of these
methods, namely thermodynamic FBA and loopless FBA, and analyze the strengths
and weaknesses of each one. Finally, we suggest a compromise denoted
semi-thermodynamic FBA (st-FBA) which imposes stronger thermodynamic constrains
on the flux polytope compared to loopless FBA, without requiring a large set of
thermodynamic parameters as in the case of thermodynamic FBA. We show that
st-FBA is a useful and simple way to eliminate thermodynamically infeasible
cycles that generate ATP.


  The single nucleotide polymorphism (SNP) is the most widely studied type of
genetic variation. A haplotype is defined as the sequence of alleles at SNP
sites on each haploid chromosome. Haplotype information is essential in
unravelling the genome-phenotype association. Haplotype assembly is a
well-known approach for reconstructing haplotypes, exploiting reads generated
by DNA sequencing devices. The Minimum Error Correction (MEC) metric is often
used for reconstruction of haplotypes from reads. However, problems with the
MEC metric have been reported. Here, we investigate the MEC approach to
demonstrate that it may result in incorrectly reconstructed haplotypes for
devices that produce error-prone long reads. Specifically, we evaluate this
approach for devices developed by Illumina, Pacific BioSciences and Oxford
Nanopore Technologies. We show that imprecise haplotypes may be reconstructed
with a lower MEC than that of the exact haplotype. The performance of MEC is
explored for different coverage levels and error rates of data. Our simulation
results reveal that in order to avoid incorrect MEC-based haplotypes, a
coverage of 25 is needed for reads generated by Pacific BioSciences RS systems.


  We explore the application of probability generating functions (PGFs) to
invasive processes, focusing on infectious disease introduced into large
populations. Our goal is to acquaint the reader with applications of PGFs,
moreso than to derive new results. PGFs help predict a number of properties
about early outbreak behavior while the population is still effectively
infinite, including the probability of an epidemic, the size distribution after
some number of generations, and the cumulative size distribution of
non-epidemic outbreaks. We show how PGFs can be used in both discrete-time and
continuous-time settings, and discuss how to use these results to infer disease
parameters from observed outbreaks. In the large population limit for
susceptible-infected-recovered (SIR) epidemics PGFs lead to survival-function
based models that are equivalent the the usual mass-action SIR models but with
fewer ODEs. We use these to explore properties such as the final size of
epidemics or even the dynamics once stochastic effects are negligible. We
target this tutorial to biologists and public health researchers who want to
learn how to apply PGFs to invasive diseases, but it could also be used in an
introductory mathematics course on PGFs. We include many exercises to help
demonstrate concepts and to give practice applying the results. We summarize
our main results in a few tables. Additionally we provide a small python
package which performs many of the relevant calculations.


  The delimitation of bioregions helps to understand historical and ecological
drivers of species distribution. In this work, we performed a network analysis
of the spatial distribution patterns of plants in south of France
(Languedoc-Roussillon and Provence-Alpes-C\^ote d'Azur) to analyze the
biogeographical structure of the French Mediterranean flora at different
scales. We used a network approach to identify and characterize biogeographical
regions, based on a large database containing 2.5 million of geolocalized plant
records corresponding to more than 3500 plant species. This methodology is
performed following five steps, from the biogeographical bipartite network
construction, to the identification of biogeographical regions under the form
of spatial network communities, the analysis of their interactions and the
identification of clusters of plant species based on the species contribution
to the biogeographical regions. First, we identified two sub-networks that
distinguish Mediterranean and temperate biota. Then, we separated eight
statistically significant bioregions that present a complex spatial structure.
Some of them are spatially well delimited, and match with particular geological
entities. On the other hand fuzzy transitions arise between adjacent bioregions
that share a common geological setting, but are spread along a climatic
gradient. The proposed network approach illustrates the biogeographical
structure of the flora in southern France, and provides precise insights into
the relationships between bioregions. This approach sheds light on ecological
drivers shaping the distribution of Mediterranean biota: the interplay between
a climatic gradient and geological substrate shapes biodiversity patterns.
Finally this work exemplifies why fragmented distributions are common in the
Mediterranean region, isolating groups of species that share a similar
eco-evolutionary history.


  We propose a new model for the time evolution of livestock commodities which
exhibits endogenous deterministic stochastic behaviour. The model is based on
the Yoccoz-Birkeland integral equation, a model first developed for studying
the time-evolution of single species with high average fertility, a relatively
short mating season and density dependent reproduction rates. This equation is
then coupled with a differential equation describing the price of a livestock
commodity driven by the unbalance between its demand and supply. At its birth
the cattle population is split into two parts: reproducing females and cattle
for butchery. The relative amount of the two is determined by the spot price of
the meat. We prove the existence of an attractor and we investigate numerically
its properties: the strange attractor existing for the original
Yoccoz-Birkeland model is persistent but its chaotic behaviour depends also
from the price evolution in an essential way.


  Genome-scale metabolic models have become a fundamental tool for examining
metabolic principles. However, metabolism is not solely characterized by the
underlying biochemical reactions and catalyzing enzymes, but also affected by
regulatory events. Since the pioneering work of Covert and co-workers as well
as Shlomi and co-workers it is debated, how regulation and metabolism
synergistically characterize a coherent cellular state. The first approaches
started from metabolic models which were extended by the regulation of the
encoding genes of the catalyzing enzymes. By now, bioinformatics databases in
principle allow addressing the challenge of integrating regulation and
metabolism on a system-wide level. Collecting information from several
databases we provide a network representation of the integrated gene regulatory
and metabolic system for Escherichia coli, including major cellular processes,
from metabolic processes via protein modification to a variety of regulatory
events. Besides transcriptional regulation, we also take into account
regulation of translation, enzyme activities and reactions. Our network model
provides novel topological characterizations of system components based on
their positions in the network. We show that network characteristics suggest a
representation of the integrated system as three network domains (regulatory,
metabolic and interface networks) instead of two. This new three-domain
representation reveals the structural centrality of components with known high
functional relevance. This integrated network can serve as a platform for
understanding coherent cellular states as active subnetworks and to elucidate
crossover effects between metabolism and gene regulation.


  Translation of mRNA into protein is a fundamental yet complex biological
process with multiple factors that can potentially affect its efficiency. Here,
we study a stochastic model describing the traffic flow of ribosomes along the
mRNA (namely, the inhomogeneous $\ell$-TASEP), and identify the key parameters
that govern the overall rate of protein synthesis, sensitivity to initiation
rate changes, and efficiency of ribosome usage. By analyzing a continuum limit
of the model, we obtain closed-form expressions for stationary currents and
ribosomal densities, which agree well with Monte Carlo simulations.
Furthermore, we completely characterize the phase transitions in the system,
and by applying our theoretical results, we formulate design principles that
detail how to tune the key parameters we identified to optimize translation
efficiency. Using ribosome profiling data from S. cerevisiae, we shows that its
translation system is generally consistent with these principles. Our
theoretical results have implications for evolutionary biology, as well as
synthetic biology.


  We adapt a fitness function from evolutionary game theory as a mechanism for
aggregation and dispersal in a partial differential equation (PDE) model of two
interacting populations, described by density functions $u$ and $v$. We
consider a spatial model where individuals migrate up local fitness gradients,
seeking out locations where their given traits are more advantageous. The
resulting system of fitness gradient equations is a degenerate system having
spatially structured, smooth, steady state solutions characterized by constant
fitness throughout the domain. When populations are viewed as predator and
prey, our model captures prey aggregation behavior consistent with Hamilton's
selfish herd hypothesis. We also present weak steady state solutions in 1d that
are continuous but in general not smooth everywhere, with an associated fitness
that is discontinuous, piecewise constant. We give numerical examples of
solutions that evolve toward such weak steady states. We also give an example
of a spatial Lotka--Volterra model, where a fitness gradient flux creates
instabilities that lead to spatially structured steady states. Our results also
suggest that when fitness has some dependence on local interactions, a
fitness-based dispersal mechanism may act to create spatial variation across a
habitat.


  A multilabeled tree (or MUL-tree) is a rooted tree in which every leaf is
labelled by an element from some set, but in which more than one leaf may be
labelled by the same element of that set. In phylogenetics, such trees are used
in biogeographical studies, to study the evolution of gene families, and also
within approaches to construct phylogenetic networks. A multilabelled tree in
which no leaf-labels are repeated is called a phylogenetic tree, and one in
which every label is the same is also known as a tree-shape. In this paper, we
consider the complexity of computing metrics on MUL-trees that are obtained by
extending metrics on phylogenetic trees. In particular, by restricting our
attention to tree shapes, we show that computing the metric extension on
MUL-trees is NP complete for two well-known metrics on phylogenetic trees,
namely, the path-difference and Robinson Foulds distances. We also show that
the extension of the Robinson Foulds distance is fixed parameter tractable with
respect to the distance parameter. The path distance complexity result allows
us to also answer an open problem concerning the complexity of solving the
quadratic assignment problem for two matrices that are a Robinson similarity
and a Robinson dissimilarity, which we show to be NP-complete. We conclude by
considering the maximum agreement subtree (MAST) distance on phylogenetic trees
to MUL-trees. Although its extension to MUL-trees can be computed in polynomial
time, we show that computing its natural generalization to more than two
MUL-trees is NP-complete, although fixed-parameter tractable in the maximum
degree when the number of given trees is bounded.


  Most biological molecules contain acido-basic groups that modulate their
structure and interactions. A consequence is that pH gradients, local
heterogeneities and dynamic variations are used by cells and organisms to drive
or regulate specific biological functions including energetic metabolism,
vesicular traffic, migration and spatial patterning of tissues in development.
While the direct or regulatory role of pH in protein function is well
documented, the role of hydrogen and hydroxyl ions in modulating the properties
of lipid assemblies such as bilayer membranes is only beginning to be
understood. Here, we review approaches using artificial lipid vesicles that
have been instrumental in providing an understanding of the influence of pH
gradients and local variations on membrane vectorial motional processes:
migration, membrane curvature effects promoting global or local deformations,
crowding generation by segregative polarization processes. In the case of pH
induced local deformations, an extensive theoretical framework is given and an
application to a specific biological issue, namely the structure and stability
of mitochondrial cristae, is described.


  The intracellular environment is a dynamic space filled with various
organelles moving in all directions. Included in this diverse group of
organelles are vesicles, which are involved in transport of molecular cargo
throughout the cell. Vesicles move in either a directed or non-directed
fashion, often depending on interactions with cytoskeletal proteins such as
microtubules, actin filaments, and molecular motors. How these proteins affect
the local fluctuations of vesicles in the cytoplasm is not clear since they
have the potential to both facilitate and impede movement. Here we show that
vesicle mobility is significantly affected by myosin-II, even though it is not
a cargo transport motor. We find that myosin-II activity increases the
effective diffusivity of vesicles and its inhibition facilitates longer states
of non-directed motion. Our study suggests that altering myosin-II activity in
the cytoplasm of cells can modulate the mobility of vesicles, providing a
possible mechanism for cells to dynamically tune the cytoplasmic environment in
space and time.


  We study that the breakdown of epidemic depends on some parameters, that is
expressed in epidemic reproduction ratio number. It is noted that when $R_0 $
exceeds 1, the stochastic model have two different results. But, eventually the
extinction will be reached even though the major epidemic occurs. The question
is how long this process will reach extinction. In this paper, we will focus on
the Markovian process of SIS model when major epidemic occurs. Using the
approximation of quasi--stationary distribution, the expected mean time of
extinction only occurs when the process is one step away from being extinct.
Combining the theorm from Ethier and Kurtz, we use CLT to find the
approximation of this quasi distribution and successfully determine the
asymptotic mean time to extinction of SIS model without demography.


  The evolution of multicellularity was a major transition in the history of
life on earth. Conditions under which multicellularity is favored have been
studied theoretically and experimentally. But since the construction of a
multicellular organism requires multiple rounds of cell division, a natural
question is whether these cell divisions should be synchronous or not. We study
a simple population model in which there compete simple multicellular organisms
that grow either by synchronous or asynchronous cell divisions. We demonstrate
that natural selection can act differently on synchronous and asynchronous cell
division, and we offer intuition for why these phenotypes are generally not
neutral variants of each other.


  Background: Networks in different domains are characterized by similar global
characteristics while differing in local structures. To further extend this
concept, we investigated network regularities on a fine scale in order to
examine the functional impact of recurring motifs in signed oriented biological
networks. In this work we generalize to signaling net works some considerations
made on feedback and feed forward loops and extend them by adding a close
scrutiny of Linear Triplets, which have not yet been investigate in detail.
Results: We studied the role of triplets, either open or closed (Loops or
linear events) by enumerating them in different biological signaling networks
and by comparing their significance profiles. We compared different data
sources and investigated the fine topology of protein networks representing
causal relationships based on transcriptional control, phosphorylation,
ubiquitination and binding. Not only were we able to generalize findings that
have already been reported but we also highlighted a connection between
relative motif abundance and node function. Furthermore, by analyzing for the
first time Linear Triplets, we highlighted the relative importance of nodes
sitting in specific positions in closed signaling triplets. Finally, we tried
to apply machine learning to show that a combination of motifs features can be
used to derive node function. Availability: The triplets counter used for this
work is available as a Cytoscape App and as a standalone command line Java
application. http://apps.cytoscape.org/apps/counttriplets Keywords: Graph
theory, graph analysis, graph topology, machine learning, cytoscape


  The recent rise of cryo-EM and X-ray high-throughput techniques is providing
a wealth of new structures trapped in different conformations. Understanding
how proteins transition between different conformers, and how they relate to
each other in terms of function is not straightforward, and highly depends on
the choice of the right set of degrees of freedom. Here we present eBDIMS
server, an online tool and software for automatic classification of structural
ensembles and reconstruction of transition pathways using coarse-grained (CG)
simulations. The server generates CG-pathways between two protein conformations
along with a representation in a simplified 2D-motion landscape based on the
Principal Components (PCs) from experimental structures. For a conformationally
rich ensemble, the PCs provide powerful reaction coordinates for automatic
structure classification, detection of on-pathway intermediates and validation
of in silico pathways. When the number of available structures is low or
sampling is limited, Normal Modes (NMs) provide alternative motion axes for
trajectory analysis. The path-generation eBDIMS method is available at a
user-friendly website: https://login.biophysics.kth.se/eBDIMS/ or as standalone
software. The server incorporates a powerful interactive graphical interface
for simultaneous visualization of transition pathways in 2D-motion space and
3D-molecular graphics, which greatly facilitates the exploration of the
relationships between different conformations.


  Distributions of anthropogenic signatures (impacts and activities) are
mathematically analysed. The aim is to understand the Anthropocene and to see
whether anthropogenic signatures could be used to determine its beginning. A
total of 23 signatures were analysed and results are presented in 31 diagrams.
Some of these signatures contain undistinguishable natural components but most
of them are of purely anthropogenic origin. Great care was taken to identify
abrupt accelerations, which could be used to determine the beginning of the
Anthropocene. Results of the analysis can be summarised in three conclusions.
1. Anthropogenic signatures cannot be used to determine the beginning of the
Anthropocene. 2. There was no abrupt Great Acceleration around 1950 or around
any other time. 3. Anthropogenic signatures are characterised by the Great
Deceleration in the second half of the 20th century. The second half of the
20th century does not mark the beginning of the Anthropocene but most likely
the beginning of the end of the strong anthropogenic impacts, maybe even the
beginning of a transition to a sustainable future. The Anthropocene is a unique
stage in human experience but it has no clearly marked beginning and it is
probably not a new geological epoch.


  A discrete electrostatic/diffusion model has been developed to describe the
selective permeation of ion channels, based on ionic Coulomb blockade (ICB) and
quantised dehydration (QD). It has been applied to describe selectivity
phenomena measured in the bacterial NaChBac sodium channel and some of its
mutants. Site-directed mutagenesis and the whole-cell patch-clamp technique
were used to investigate how the value $Q_f$ of the fixed charge at the
selectivity filter (SF) affected both valence and alike-charge selectivity. The
new ICB/QD model predicts that increasing ${Q_f}$ should lead to a shift of
selectivity sequences towards larger ion sizes and charges, a result that
agrees with the present experiments and with earlier work. Comparison of the
model with experimental data provides evidence for an {\it effective charge}
$Q_f^*$ at the SF that is smaller in magnitude than the nominal $Q_f$
corresponding to the charge on the isolated protein residues. Furthermore,
$Q_f^*$ was different for aspartate and glutamate charged rings and also
depended on their position within the SF. It is suggested that protonation of
the residues within the restricted space is an important factor in
significantly reducing the effective charge of the EEEE ring. Values of $Q_f^*$
derived from experiments on the anomalous mole fraction effect (AMFE) agree
well with expectations based on the ICB/QD model and have led to the first
clear demonstration of the expected ICB oscillations in Ca$^{2+}$ conduction as
a function of the fixed charge. Pilot studies of the dependence of Ca$^{2+}$
conduction on pH are consistent with the predictions of the model.


  In the present study, we continue investigation of the high-contrast vision
in the inverted retina of the vertebrates eyes. We report a method of
separation and purification of porcine (Sus scrofa domestica) intermediate
filaments (IFs), extracted from the retinal M\"uller cells (MCs). We also
report experimental and theoretical methods of measurements and calculations of
the reduced resistivity and light transmission by the IFs and single-wall
carbon nanotubes (SWCNTs). The measured reduced resistivity values were
(3.1+-0.3)*10^-4 and (2.8+-0.2)*10^-4 Ohm m^-1 cm^2, respectively, being quite
close to those of typical metals. We report a method for measuring the light
energy transmission by the intermediate filaments and single-wall carbon
nanotubes. We found that these structures efficiently transfer light energy
along its axis, with the light reemitted at the other end of the structure. We
also report spectral selectivity of the IFs. The reported results demonstrate
that the assumptions we made in deducing the theory of high-contrast vision in
an inverted retina were correct and fully supported by the presently reported
experimental results.


  While the mechanisms and quantitative details are complex, few analysts would
doubt that antibiotic use increases the prevalence of drug-resistant bacterial
pathogens among all bacteria causing disease in a population. The causal
connection between antibiotic use and the total incidence of severe bacterial
infections, possibly mediated by antibiotic resistance, is less clearly
established. The increasing burden of severe bacterial infections and their
sequelae, particularly sepsis, in the United States and other countries, calls
out for an explanation. In this Perspective we consider the evidence bearing on
the hypothesis that prevalence of antibiotic resistance and levels of
antibiotic use are important contributors to the rates of sepsis
hospitalizations and other outcomes with bacterial etiology. In the process, we
discuss the consequences of resistance to/use of commonly prescribed
antibiotics, including fluoroquinolones, and provide a comparison of the
epidemiology of antibiotic use/resistance and severe outcomes associated with
infections with bacterial pathogens, particularly Escherichia coli,
methicillin-resistant Staphylococcus aureus (MRSA), and Clostridium Difficile
in the UK vs. US.


  Coevolutionary game dynamics is the result of players that may change their
strategies and their network of interaction. For such games, and based on
interpreting strategies as configurations, strategy-to-payoff maps can be
defined for every interaction network, which opens up to derive game
landscapes. This paper presents an analysis of these game landscapes by their
information content. By this analysis, we particularly study the effect of a
rescaled payoff matrix generalizing social dilemmas and differences between
well-mixed and structured populations.


  Microtubules are filamentous tubular protein polymers which are essential for
a range of cellular behaviour, and are generally straight over micron length
scales. However, in some gliding assays, where microtubules move over a carpet
of molecular motors, individual microtubules can also form tight arcs or rings,
even in the absence of crosslinking proteins. Understanding this phenomenon may
provide important explanations for similar highly curved microtubules which can
be found in nerve cells undergoing neurodegeneration. We propose a model for
gliding assays where the kinesins moving the microtubules over the surface
induce ring formation through differential binding, substantiated by recent
findings that a mutant version of the motor protein kinesin applied in solution
is able to lock-in microtubule curvature. For certain parameter regimes, our
model predicts that both straight and curved microtubules can exist
simultaneously as stable steady-states, as has been seen experimentally.
Additionally, unsteady solutions are found, where a wave of differential
binding propagates down the microtubule as it glides across the surface, which
can lead to chaotic motion. Whilst this model explains two-dimensional
microtubule behaviour in an experimental gliding assay, it has the potential to
be adapted to explain pathological curling in nerve cells.


  Observed bimodal tree cover distributions at particular environmental
conditions and theoretical models indicate that some areas in the tropics can
be in either of the alternative stable vegetation states forest or savanna.
However, when including spatial interaction in nonspatial differential equation
models of a bistable quantity, only the state with the lowest potential energy
remains stable. Our recent reaction-diffusion model of Amazonian tree cover
confirmed this and was able to reproduce the observed spatial distribution of
forest versus savanna satisfactorily when forced by heterogeneous environmental
and anthropogenic variables, even though bistability was underestimated. These
conclusions were solely based on simulation results. Here, we perform an
analytical and numerical analysis of the model. We derive the Maxwell point
(MP) of the homogeneous reaction-diffusion equation without savanna trees as a
function of rainfall and human impact and show that the front between forest
and nonforest settles at this point as long as savanna tree cover near the
front remains sufficiently low. For parameters resulting in higher savanna tree
cover near the front, we also find irregular forest-savanna cycles and
woodland-savanna bistability, which can both explain the remaining observed
bimodality.


  Shifting ecosystem disturbance patterns due to climate change (e.g. storms,
droughts, wildfires) or direct human interference (e.g. harvests, nutrient
loading) highlight the importance of quantifying and strengthening the
resilience of desired ecological regimes. Although existing metrics capture
resilience to isolated shocks, gradual parameter changes, and continuous noise,
quantifying resilience to repeated, discrete disturbances requires novel
analytical tools. Here we introduce a flow-kick framework that quantifies
resilience to disturbances explicitly in terms of their magnitude and
frequency. We present a resilience boundary between disturbances that cause
either escape from a basin of attraction or stabilization within it, and use
the resilience boundary to build resilience metrics tailored to repeated,
discrete disturbances. The flow-kick model suggests that the
distance-to-threshold resilience metric overestimates resilience in the context
of repeated disturbances. It also reveals counterintuitive triggers for regime
shifts, such as increasing recovery times between disturbances, or increasing
disturbance magnitude and recovery times proportionately.


  We study how the complexity of evolutionary dynamics in the classic MacArthur
consumer-resource model depends on resource uptake and utilization rates. The
traditional assumption in such models is that the utilization rate of the
consumer is proportional to the uptake rate. More generally, we show that if
these two rates are related through a power law (which includes the traditional
assumption as a special case), then the resulting evolutionary dynamics in the
consumer is necessarily a simple hill-climbing process leading to an
evolutionary equilibrium, regardless of the dimension of phenotype space. When
utilization and uptake rates are not related by a power law, more complex
evolutionary trajectories can occur, including the chaotic dynamics observed in
previous studies for high-dimensional phenotype spaces. These results draw
attention to the importance of distinguishing between utilization and uptake
rates in consumer-resource models.


  The Turing patterning mechanism is believed to underly the formation of
repetitive structures in development, such as zebrafish stripes and mammalian
digits, but it has proved difficult to isolate the specific biochemical species
responsible for pattern formation. Meanwhile, synthetic biologists have
designed Turing systems for implementation in cell colonies, but none have yet
led to visible patterns in the laboratory. In both cases, the relationship
between underlying chemistry and emergent biology remains mysterious. To help
resolve the mystery, this article asks the question: what kinds of biochemical
systems can generate Turing patterns? We find general conditions for Turing
pattern inception -- the ability to generate unstable patterns from random
noise -- which may lead to the ultimate formation of stable patterns, depending
on biochemical non-linearities. We find that a wide variety of systems can
generate stable Turing patterns, including several which are currently unknown,
such as two-species systems composed of two self-activators, and systems
composed of a short-range inhibitor and a long-range activator. We furthermore
find that systems which are widely believed to generate stable patterns may in
fact only generate unstable patterns, which ultimately converge to
spatially-homogeneous concentrations. Our results suggest that a much wider
variety of systems than is commonly believed could be responsible for observed
patterns in development, or could be good candidates for synthetic patterning
networks.


  Herbaria worldwide are housing a treasure of 100s of millions of herbarium
specimens, which are increasingly being digitized in recent years and thereby
made more easily accessible to the scientific community. At the same time, deep
learning algorithms are rapidly improving pattern recognition from images and
these techniques are more and more being applied to biological objects. We are
using digital images of herbarium specimens in order to identify taxa and
traits of these collection objects by applying convolutional neural networks
(CNN). Images of the 1000 species most frequently documented by herbarium
specimens on GBIF have been downloaded and combined with morphological trait
data, preprocessed and divided into training and test datasets for species and
trait recognition. Good performance in both domains is promising to use this
approach in future tools supporting taxonomy and natural history collection
management.


  Interaction patterns among individuals play vital roles in spreading
infectious diseases. Understanding these patterns and integrating their impact
in modeling diffusion dynamics of infectious diseases are important for
epidemiological studies. Current network-based diffusion models assume that
diseases transmit through interactions where both infected and susceptible
individuals are co-located at the same time. However, there are several
infectious diseases that can transmit when a susceptible individual visits a
location after an infected individual has left. Recently, we introduced a
diffusion model called same place different time (SPDT) transmission to capture
the indirect transmissions that happen when an infected individual leaves
before a susceptible individual's arrival along with direct transmissions. In
this paper, we demonstrate how these indirect transmission links significantly
enhance the emergence of infectious diseases simulating airborne disease
spreading on a synthetic social contact network. We denote individuals having
indirect links but no direct links during their infectious periods as hidden
spreaders. Our simulation shows that indirect links play similar roles of
direct links and a single hidden spreader can cause large outbreak in the SPDT
model which causes no infection in the current model based on direct link. Our
work opens new direction in modeling infectious diseases.


  As a living information and communications system, the genome encodes
patterns in single nucleotide polymorphisms (SNPs) reflecting human adaption
that optimizes population survival in differing environments. This paper
mathematically models environmentally induced adaptive forces that quantify
changes in the distribution of SNP frequencies between populations. We make
direct connections between biophysical methods (e.g. minimizing genomic free
energy) and concepts in population genetics. Our unbiased computer program
scanned a large set of SNPs in the major histocompatibility complex region, and
flagged an altitude dependency on a SNP associated with response to oxygen
deprivation. The statistical power of our double-blind approach is demonstrated
in the flagging of mathematical functional correlations of SNP
information-based potentials in multiple populations with specific
environmental parameters. Furthermore, our approach provides insights for new
discoveries on the biology of common variants. This paper demonstrates the
power of biophysical modeling of population diversity for better understanding
genome-environment interactions in biological phenomenon.


  We introduce a model for the global optimization problem of nectar harvesting
by flower visitors, e.g., nectar-feeding bats, as a generalization of the
(multiple) traveling-salesperson problem (TSP). The model includes multiple
independent animals and many flowers with time-dependent content. This provides
an ensemble of realistic combinatorial optimization problems, in contrast to
previously studied models like random Satisfiability or standard TSP.
  We numerically studied the optimum harvesting of these foragers, with
parameters obtained from experiments, by using genetic algorithms. For the
distribution of travel distances, we find a power-law (or L\'evy) distribution,
as often found for natural foragers. Note, in contrast to many models, we make
no assumption about the nature of the flight-distance distribution, the power
law just emerges. This is in contrast to the TSP, where we find in the present
study an exponential tail.
  Furthermore, the optimization problem exhibits a {phase transition}, similar
to the TSP, at a critical value for the amount of nectar which can be
harvested. This phase transition coincides with a dramatic increase in the
typical running time of the optimization algorithm. For the value of the
critical exponent nu, describing the divergence of the correlation length, we
find nu=1.7(4), which is on the other hand compatible with the value found for
the TSP.
  Finally, we also present data from field experiments in Costa Rica for the
resource use for freely visiting flower bats. We found that the temporal
patterns in experiments and model agree remarkably, confirming our model. Also
the data show that the bats are able to memorize the positions of food sources
and optimize, at least partially, their routes.


  Many theoretical works have attempted to coarse grain gene expression at the
level of transcription and translation via frameworks based on exclusion
processes. Usually in these models the three-dimensional conformation of the
substrates (DNA and mRNA) is neglected, and particles move on a static
unidimensional lattice in contact to an infinite reservoir. In this work we
generalise the paradigmatic exclusion process and study the transport of
particles along a unidimensional polymer-like flexible lattice immersed in a
three-dimensional particle reservoir. We study the recycling of particles in
the reservoir, how the transport is influenced by the global conformation of
the lattice and, in turn, how particle density dictates the structure of the
polymer.


  Multipartite viruses replicate through a puzzling evolutionary strategy.
Their genome is segmented into two or more parts, and encapsidated in separate
particles that appear to propagate independently. Completing the replication
cycle, however, requires the full genome, so that a systemic infection of a
host requires the concurrent presence of several particles. This represents an
apparent evolutionary drawback of multipartitism, while its advantages remain
unclear. A transition from monopartite to multipartite viral forms has been
described in vitro under conditions of high multiplicity of infection,
suggesting that cooperation between defective mutants is a plausible
evolutionary pathway towards multipartitism. However, it is unknown how the
putative advantages that multipartitism might enjoy at the microscopic level
affect its epidemiology, or if an explicit advantange is needed to explain its
ecological persistence. To disentangle which mechanisms might contribute to the
rise and fixation of multipartitism, we investigate the interaction between
viral spreading dynamics and host population structure. We set up a
compartmental model of the spread of a virus in its different forms and explore
its epidemiology using both analytical and numerical techniques. We uncover
that the impact of host contact structure on spreading dynamics entails a rich
phenomenology of ecological relationships that includes cooperation,
competition, and commensality. We find that multipartitism might rise to
fixation even in the absence of explicit microscopic advantages. Multipartitism
allows the virus to colonize environments that could not be invaded by the
monopartite form, facilitated by homogeneous contacts among hosts. We
conjecture that these features might have led to an increase in the diversity
and prevalence of multipartite viral forms concomitantly with the expansion of
agricultural practices.


  Behavioural differences may arise in the absence of genetic or environmental
variation. Chaotic dynamics may influence behavioural development, and so this
among-individual variation. We discuss methods and experimental designs to test
this idea. Ultimately, nonlinear and chaotic behavioural development may
explain much of natural variation.


  The evolution of microbial and viral organisms often generates clonal
interference, a mode of competition between genetic clades within a population.
In this paper, we show that interference strongly constrains the genetic and
phenotypic complexity of evolving systems. Our analysis uses biophysically
grounded evolutionary models for an organism's quantitative molecular
phenotypes, such as fold stability and enzymatic activity of genes. We find a
generic mode of asexual evolution called phenotypic interference with strong
implications for systems biology: it couples the stability and function of
individual genes to the population's global speed of evolution. This mode
occurs over a wide range of evolutionary parameters appropriate for microbial
populations. It generates selection against genome complexity, because the
fitness cost of mutations increases faster than linearly with the number of
genes. Recombination can generate a distinct mode of sexual evolution that
eliminates the superlinear cost. We show that positive selection can drive a
transition from asexual to facultative sexual evolution, providing a specific,
biophysically grounded scenario for the evolution of sex. In a broader context,
our analysis suggests that the systems biology of microbial organisms is
strongly intertwined with their mode of evolution.


  These lecture notes introduce key concepts of mathematical population
genetics within the most elementary setting and describe a few recent
applications to microbial evolution experiments. Pointers to the literature for
further reading are provided, and some of the derivations are left as exercises
for the reader.


  Background: The chromatin remodelers of the SWI/SNF family are critical
transcriptional regulators. Recognition of lysine acetylation through a
bromodomain (BRD) component is key to SWI/SNF function; in most eukaryotes,
this function is attributed to SNF2/Brg1.
  Results: Using affinity purification coupled to mass spectrometry (AP-MS) we
identified members of a SWI/SNF complex (SWI/SNFTt) in Tetrahymena thermophila.
SWI/SNFTt is composed of 11 proteins, Snf5Tt, Swi1Tt, Swi3Tt, Snf12Tt, Brg1Tt,
two proteins with potential chromatin interacting domains and four proteins
without orthologs to SWI/SNF proteins in yeast or mammals. SWI/SNFTt subunits
localize exclusively to the transcriptionally active macronucleus (MAC) during
growth and development, consistent with a role in transcription. While
Tetrahymena Brg1 does not contain a BRD, our AP-MS results identified a
BRD-containing SWI/SNFTt component, Ibd1 that associates with SWI/SNFTt during
growth but not development. AP-MS analysis of epitope-tagged Ibd1 revealed it
to be a subunit of several additional protein complexes, including putative
SWRTt, and SAGATt complexes as well as a putative H3K4-specific histone methyl
transferase complex. Recombinant Ibd1 recognizes acetyl-lysine marks on
histones correlated with active transcription. Consistent with our AP-MS and
histone array data suggesting a role in regulation of gene expression, ChIP-Seq
analysis of Ibd1 indicated that it primarily binds near promoters and within
gene bodies of highly expressed genes during growth.
  Conclusions: Our results suggest that through recognizing specific histones
marks, Ibd1 targets active chromatin regions of highly expressed genes in
Tetrahymena where it subsequently might coordinate the recruitment of several
chromatin remodeling complexes to regulate the transcriptional landscape of
vegetatively growing Tetrahymena cells.


  We use a stochastic birth-death model for a population of cells to estimate
the normal tissue complication probability (NTCP) under a particular
radiotherapy protocol. We specifically allow for interaction between cells, via
a nonlinear logistic growth model. To capture some of the effects of intrinsic
noise in the population we develop several approximations of NTCP, using
Kramers-Moyal expansion techniques. These approaches provide an approximation
to the first and second moments of a general first-passage time problem in the
limit of large, but finite populations. We use this method to study NTCP in a
simple model of normal cells and in a model of normal and damaged cells. We
also study a combined model of normal tissue cells and tumour cells. Based on
existing methods to calculate tumour control probabilities, and our procedure
to approximate NTCP, we estimate the probability of complication free tumour
control.