  In genetic circuits, when the mRNA lifetime is short compared to the cell
cycle, proteins are produced in geometrically-distributed bursts, which greatly
affects the cellular switching dynamics between different metastable phenotypic
states. Motivated by this scenario, we study a general problem of switching or
escape in stochastic populations, where influx of particles occurs in groups or
bursts, sampled from an arbitrary distribution. The fact that the step size of
the influx reaction is a-priori unknown, and in general, may fluctuate in time
with a given correlation time and statistics, introduces an additional
non-demographic step-size noise into the system. Employing the probability
generating function technique in conjunction with Hamiltonian formulation, we
are able to map the problem in the leading order onto solving a stationary
Hamilton-Jacobi equation. We show that bursty influx exponentially decreases
the mean escape time compared to the "usual case" of single-step influx. In
particular, close to bifurcation we find a simple analytical expression for the
mean escape time, which solely depends on the mean and variance of the
burst-size distribution. Our results are demonstrated on several realistic
distributions and compare well with numerical Monte-Carlo simulations.


  We develop a feature allocation model for inference on genetic tumor
variation using next-generation sequencing data. Specifically, we record single
nucleotide variants (SNVs) based on short reads mapped to human reference
genome and characterize tumor heterogeneity by latent haplotypes defined as a
scaffold of SNVs on the same homologous genome. For multiple samples from a
single tumor, assuming that each sample is composed of some sample-specific
proportions of these haplotypes, we then fit the observed variant allele
fractions of SNVs for each sample and estimate the proportions of haplotypes.
Varying proportions of haplotypes across samples is evidence of tumor
heterogeneity since it implies varying composition of cell subpopulations.
Taking a Bayesian perspective, we proceed with a prior probability model for
all relevant unknown quantities, including, in particular, a prior probability
model on the binary indicators that characterize the latent haplotypes. Such
prior models are known as feature allocation models. Specifically, we define a
simplified version of the Indian buffet process, one of the most traditional
feature allocation models. The proposed model allows overlapping clustering of
SNVs in defining latent haplotypes, which reflects the evolutionary process of
subclonal expansion in tumor samples.


  Questions in computational molecular biology generate various discrete
optimization problems, such as DNA sequence alignment and RNA secondary
structure prediction. However, the optimal solutions are fundamentally
dependent on the parameters used in the objective functions. The goal of a
parametric analysis is to elucidate such dependencies, especially as they
pertain to the accuracy and robustness of the optimal solutions. Techniques
from geometric combinatorics, including polytopes and their normal fans, have
been used previously to give parametric analyses of simple models for DNA
sequence alignment and RNA branching configurations. Here, we present a new
computational framework, and proof-of-principle results, which give the first
complete parametric analysis of the branching portion of the nearest neighbor
thermodynamic model for secondary structure prediction for real RNA sequences.


  In higher organisms, all cells share the same genome, but every cell
expresses only a limited and specific set of genes that defines the cell type.
During cell division, not only the genome, but also the cell type is inherited
by the daughter cells. This intriguing phenomenon is achieved by a variety of
processes that have been collectively termed epigenetics: the stable and
inheritable changes in gene expression patterns. This article reviews the
extremely rich and exquisitely multi-scale physical mechanisms that govern the
biological processes behind the initiation, spreading and inheritance of
epigenetic states. These include not only the changes in the molecular
properties associated with the chemical modifications of DNA and histone
proteins, such as methylation and acetylation, but also less conventional ones,
such as the physics that governs the three-dimensional organization of the
genome in cell nuclei. Strikingly, to achieve stability and heritability of
epigenetic states, cells take advantage of many different physical principles,
such as the universal behavior of polymers and copolymers, the general features
of non-equilibrium dynamical systems, and the electrostatic and mechanical
properties related to chemical modifications of DNA and histones. By putting
the complex biological literature under this new light, the emerging picture is
that a limited set of general physical rules play a key role in initiating,
shaping and transmitting this crucial "epigenetic landscape". This new
perspective not only allows to rationalize the normal cellular functions, but
also helps to understand the emergence of pathological states, in which the
epigenetic landscape becomes dysfunctional.


  Context-dependent identification of orthologs customarily relies on conserved
gene order or whole-genome sequence alignment. It is shown here that
short-range context--as short as single maximal matches--also provides an
effective means to identify orthologs within whole genomes. On pristine
(un-repeatmasked) mammalian whole-genome assemblies we perform a genome
"intersection" that in general consumes less than one thirtieth of the
computation time required by commonly used methods for whole-genome alignment,
and we extract "non-embedded maximal matches," maximal matches that are not
embedded into other maximal matches, as potential orthologs. An ortholog
identified via non-embedded maximal matches is analogous to a "positional
ortholog" or a "primary ortholog" as defined in previous literature; such
orthologs constitute homologs derived from the same direct ancestor whose
ancestral positions in the genome are conserved. At the nucleotide level,
non-embedded maximal matches recapitulate most exact matches identified by a
Lastz net alignment. At the gene level, reciprocal best hits of genes
containing non-embedded maximal matches recover one-to-one orthologs annotated
by Ensembl Compara with high selectivity and high sensitivity; these reciprocal
best hits additionally include putatively novel orthologs not found in Ensembl
(e.g. over two thousand for human/chimpanzee). The method is especially
suitable for genome-wide identification of orthologs.


  Inside individual cells, expression of genes is inherently stochastic and
manifests as cell-to-cell variability or noise in protein copy numbers. Since
proteins half-lives can be comparable to the cell-cycle length, randomness in
cell-division times generates additional intercellular variability in protein
levels. Moreover, as many mRNA/protein species are expressed at low-copy
numbers, errors incurred in partitioning of molecules between the mother and
daughter cells are significant. We derive analytical formulas for the total
noise in protein levels for a general class of cell-division time and
partitioning error distributions. Using a novel hybrid approach the total noise
is decomposed into components arising from i) stochastic expression; ii)
partitioning errors at the time of cell-division and iii) random cell-division
events. These formulas reveal that random cell-division times not only generate
additional extrinsic noise but also critically affect the mean protein copy
numbers and intrinsic noise components. Counter intuitively, in some parameter
regimes noise in protein levels can decrease as cell-division times become more
stochastic. Computations are extended to consider genome duplication, where the
gene dosage is increased by two-fold at a random point in the cell-cycle. We
systematically investigate how the timing of genome duplication influences
different protein noise components. Intriguingly, results show that noise
contribution from stochastic expression is minimized at an optimal genome
duplication time. Our theoretical results motivate new experimental methods for
decomposing protein noise levels from single-cell expression data.
Characterizing the contributions of individual noise mechanisms will lead to
precise estimates of gene expression parameters and techniques for altering
stochasticity to change phenotype of individual cells.


  Recent theoretical studies have shown that demographic stochasticity can
greatly increase the tendency of asexually reproducing phenotypically diverse
organisms to spontaneously evolve into localised clusters, suggesting a simple
mechanism for sympatric speciation. Here we study the role of demographic
stochasticity in a model of competing organisms subject to assortative mating.
We find that in models with sexual reproduction, noise can also lead to the
formation of phenotypic clusters in parameter ranges where deterministic models
would lead to a homogeneous distribution. In some cases, noise can have a
sizeable effect, rendering the deterministic modelling insufficient to
understand the phenotypic distribution.


  Bovine tuberculosis (TB) is a chronic disease in cattle that causes a serious
food security challenge to the agricultural industry in terms of dairy and meat
production. In GB, Scotland has had a risk based surveillance testing policy
under which high risk herds are tested frequently, and in Sept 2009 was
officially declared as TB free. Wales have had an annual or more frequent
testing policy for all cattle herds since Jan 2010, while in England several
herds are still tested every 4 years except some high TB prevalence areas where
annual testing is applied. Time series analysis using publicly available data
for total tests on herds, total cattle slaughtered, new herd incidents, and
herds not TB free, were analysed globally for GB and locally for the
constituent regions of Wales, Scotland, West, North, and East England. After
detecting trends over time, underlying regional differences were compared with
the testing policies in the region. Total cattle slaughtered are decreasing in
Wales, Scotland and West England, but increasing in the North and East English
regions. New herd incidents, i.e., disease incidence, are decreasing in Wales,
Scotland, West English region, but increasing in North and East English
regions. Herds not TB free, are increasing in West, North, and East English
regions, while they are decreasing in Wales and Scotland. Total cattle
slaughtered were positively correlated with total tests in the West, North, and
East English regions, with high slopes of regression. There was no correlation
between total cattle slaughtered and total tests on herds in Wales indicating
that herds are tested frequent enough in order to detect all likely cases and
so control TB. The main conclusion of the analysis conducted here is that more
frequent testing is leading to lower TB infections in cattle both in terms of
TB prevalence as well as TB incidence.


  Fine-tuned regulation of the cellular nucleotide pools is indispensable for
faithful replication of DNA. The genetic information is also safeguarded by DNA
damage recognition and repair processes. Uracil is one of the most frequently
occurring erroneous base in DNA; it can arise from cytosine deamination or
thymine-replacing incorporation. Two enzyme families are primarily involved in
keeping DNA uracil-free: dUTPases that prevent thymine-replacing incorporation
and uracil-DNA glycosylases that excise uracil from DNA and initiate
uracil-excision repair. Both dUTPase and the most efficient uracil-DNA
glycosylase UNG is thought to be ubiquitous in free-living organisms. In the
present work, we have systematically investigated the genotype of deposited
fully sequenced bacterial and Archaeal genomes. Surprisingly, we have found
that in contrast to the generally held opinion, a wide number of bacterial and
Archaeal species lack the dUTPase gene(s). The dut- genotype is present in
diverse bacterial phyla indicating that loss of this (or these) gene(s) has
occurred multiple times during evolution. We have identified several survival
strategies in lack of dUTPases: i) simultaneous lack or inhibition of UNG, ii)
acquisition of a less dUTP-specific sanitizing nucleotide pyrophosphatase, and
iii) supply of dUTPase from bacteriophages. Our data indicate that several
unicellular microorganisms may efficiently cope with a dut- genotype
potentially leading to an unusual uracil-enrichment in their genomic DNA.


  In this paper, we derive and analyze a compartmental model for the control of
arboviral diseases which takes into account an imperfect vaccine combined with
individual protection and some vector control strategies already studied in the
literature. After the formulation of the model, a qualitative study based on
stability analysis and bifurcation theory reveals that the phenomenon of
backward bifurcation may occur. The stable disease-free equilibrium of the
model coexists with a stable endemic equilibrium when the reproduction number,
R 0 , is less than unity. Using Lyapunov function theory, we prove that the
trivial equilibrium is globally asymptotically stable; When the disease--
induced death is not considered, or/and, when the standard incidence is
replaced by the mass action incidence, the backward bifurcation does not occur.
Under a certain condition , we establish the global asymptotic stability of the
disease--free equilibrium of the full model. Through sensitivity analysis, we
determine the relative importance of model parameters for disease transmission.
Numerical simulations show that the combination of several control mechanisms
would significantly reduce the spread of the disease, if we maintain the level
of each control high, and this, over a long period.


  Community detection (CD) algorithms are applied to Hi-C data to discover new
communities of loci in the 3D conformation of human and mouse DNA. We find that
CD has some distinct advantages over pre-existing methods: (1) it is capable of
finding a variable number of communities, (2) it can detect communities of DNA
loci either adjacent or distant in the 1D sequence, and (3) it allows us to
obtain a principled value of k, the number of communities present. Forcing k =
2, our method recovers earlier findings of Lieberman-Aiden, et al. (2009), but
letting k be a parameter, our method obtains as optimal value k = 6,
discovering new candidate communities. In addition to discovering large
communities that partition entire chromosomes, we also show that CD can detect
small-scale topologically associating domains (TADs) such as those found in
Dixon, et al. (2012). CD thus provides a natural and flexible statistical
framework for understanding the folding structure of DNA at multiple scales in
Hi-C data.


  Minimal absent words (MAW) of a genomic sequence are subsequences that are
absent themselves but the subwords of which are all present in the sequence.
The characteristic distribution of genomic MAWs as a function of their length
has been observed to be qualitatively similar for all living organisms, the
bulk being rather short, and only relatively few being long. It has been an
open issue whether the reason behind this phenomenon is statistical or reflects
a biological mechanism, and what biological information is contained in absent
words. In this work we demonstrate that the bulk can be described by a
probabilistic model of sampling words from random sequences, while the tail of
long MAWs is of biological origin. We introduce the novel concept of a core of
a minimal absent word, which are sequences present in the genome and closest to
a given MAW. We show that in bacteria and yeast the cores of the longest MAWs,
which exist in two or more copies, are located in highly conserved regions the
most prominent example being ribosomal RNAs (rRNAs). We also show that while
the distribution of the cores of long MAWs is roughly uniform over these
genomes on a coarse-grained level, on a more detailed level it is strongly
enhanced in 3' untranslated regions (UTRs) and, to a lesser extent, also in 5'
UTRs. This indicates that MAWs and associated MAW cores correspond to
fine-tuned evolutionary relationships, and suggest that they can be more widely
used as markers for genomic complexity.


  A crop can be viewed as a complex system with outputs (e.g. yield) that are
affected by inputs of genetic, physiology, pedo-climatic and management
information. Application of numerical methods for model exploration assist in
evaluating the major most influential inputs, providing the simulation model is
a credible description of the biological system. A sensitivity analysis was
used to assess the simulated impact on yield of a suite of traits involved in
major processes of crop growth and development, and to evaluate how the
simulated value of such traits varies across environments and in relation to
other traits (which can be interpreted as a virtual change in genetic
background). The study focused on wheat in Australia, with an emphasis on
adaptation to low rainfall conditions. A large set of traits (90) was evaluated
in a wide target population of environments (4 sites x 125 years), management
practices (3 sowing dates x 2 N fertilization) and $CO_2$ (2 levels). The
Morris sensitivity analysis method was used to sample the parameter space and
reduce computational requirements, while maintaining a realistic representation
of the targeted trait x environment x management landscape ($\sim$ 82 million
individual simulations in total). The patterns of parameter x environment x
management interactions were investigated for the most influential parameters,
considering a potential genetic range of +/- 20% compared to a reference. Main
(i.e. linear) and interaction (i.e. non-linear and interaction) sensitivity
indices calculated for most of APSIM-Wheat parameters allowed the identifcation
of 42 parameters substantially impacting yield in most target environments.
Among these, a subset of parameters related to phenology, resource acquisition,
resource use efficiency and biomass allocation were identified as potential
candidates for crop (and model) improvement.


  Nonlinear effects in protein dynamics are expected to play role in function,
particularly of allosteric nature, by facilitating energy transfer between
vibrational modes. A recently proposed method focusing on the non-Gaussian
shape of the population near equilibrium projects this information onto real
space in order to identify the aminoacids relevant to function. We here apply
this method to three ancestral proteins in glucocorticoid receptor (GR) family
and show that the mutations that restrict functional activity during GR
evolution correlate significantly with locations that are highlighted by the
nonlinear contribution to the near-native configurational distribution. Our
findings demonstrate that nonlinear effects are not only indispensible for
understanding functionality in proteins, but they can also be harnessed into a
predictive tool for functional site determination.


  We study a system of $N\gg 1$ degrees of freedom coupled via a smooth
homogeneous Gaussian vector field with both gradient and divergence-free
components. In the absence of coupling, the system is exponentially relaxing to
an equilibrium with rate $\mu$. We show that, while increasing the ratio of the
coupling strength to the relaxation rate, the system experiences an abrupt
transition from a topologically trivial phase portrait with a single
equilibrium into a topologically non-trivial regime characterised by an
exponential number of equilibria, the vast majority of which are expected to be
unstable. It is suggested that this picture provides a global view on the
nature of the May-Wigner instability transition originally discovered by local
linear stability analysis.


  We study the evolution of cooperation in a finite population interacting
according to a simple model of like-with-like assortment. Evolution proceeds as
a Moran process, and payoffs from the underlying cooperator-defector game are
translated to positive fitnesses by an exponential transformation. These
evolutionary dynamics can arise, for example, in a nest-structured population
with rare migration. The use of the exponential transformation, rather than the
usual linear one, is appropriate when interactions have multiplicative fitness
effects, and allows for a tractable characterization of the effect of
assortment on the evolution of cooperation. We define two senses in which a
greater degree of assortment can favour the evolution of cooperation, the first
stronger than the second: (i) greater assortment increases, at all population
states, the probability that the number of cooperators increases, relative to
the probability that the number of defectors increases; and (ii) greater
assortment increases the fixation probability of cooperation, relative to that
of defection. We show that, by the stronger definition, greater assortment
favours the evolution of cooperation for a subset of cooperative dilemmas:
prisoners' dilemmas, snowdrift games, stag-hunt games, and some prisoners'
delight games. For other cooperative dilemmas, greater assortment favours
cooperation by the weak definition, but not by the strong definition. Our
results hold for any strength of selection.


  In this paper, we consider the bacterial point-to-point and multiple-access
molecular communications with ligand-receptors. For the point-to-point
communication, we investigate common signaling methods, namely the Level
Scenario (LS), which uses one type of a molecule with different concentration
levels, and the Type Scenario (TS), which employs multiple types of molecules
with a single concentration level. We investigate the trade-offs between the
two scenarios from the capacity point of view. We derive an upper bound on the
capacity using a Binomial Channel (BIC) model and the symmetrized
Kullback-Leibler (KL) divergence. A lower bound is also derived when the
environment noise is negligible. For the TS, we also consider the effect of
blocking of a receptor by a different molecule type. Then, we consider
multiple-access communications, for which we investigate three scenarios based
on molecule and receptor types, i.e., same types of molecules with Different
Labeling and Same types of Receptors (DLSR), Different types of Molecules and
Receptors (DMDR), and Same types of Molecules and Receptors (SMSR). We
investigate the trade-offs among the three scenarios from the total capacity
point of view. We derive some inner bounds on the capacity region of these
scenarios when the environment noise is negligible.


  Admixed populations are formed by the merging of two or more ancestral
populations, and the ancestry of each locus in an admixed genome derives from
either source. Consider a simple "pulse" admixture model, where populations A
and B merged t generations ago without subsequent gene flow. We derive the
distribution of the proportion of an admixed chromosome that has A (or B)
ancestry, as a function of the chromosome length L, t, and the initial
contribution of the A source, m. We demonstrate that these results can be used
for inference of the admixture parameters. For more complex admixture models,
we derive an expression in Laplace space for the distribution of ancestry
proportions that depends on having the distribution of the lengths of segments
of each ancestry. We obtain explicit results for the special case of a
"two-wave" admixture model, where population A contributed additional migrants
in one of the generations between the present and the initial admixture event.
Specifically, we derive formulas for the distribution of A and B segment
lengths and numerical results for the distribution of ancestry proportions. We
show that for recent admixture, data generated under a two-wave model can
hardly be distinguished from that generated under a pulse model.


  A theoretical framework for protocooperative behavior in pelotons (groups of
cyclists) is proposed. A threshold between cooperative and free-riding
behaviors in pelotons is modeled, together comprising protocooperative behavior
(different from protocooperation). Protocooperative behavior is a function of:
1. two or more cyclists coupled by drafting benefit, 2. cyclists' current power
output or speed, and 3. cyclists' maximal sustainable outputs (MSO).
Characteristics of protocooperative behavior include: 1. relatively low speed
phase in which cyclists naturally pass each other and share highest-cost front
positions, and 2. free-riding phase in which cyclists maintain speeds of those
ahead, but cannot pass. Threshold for protocooperative behavior is equivalent
to coefficient of drafting d, below which cooperative behavior occurs; above
which free-riding occurs up to a second threshold when coupled cyclists
diverge. These behaviors are also hypothesized to emerge in other biological
systems involving energy savings mechanisms. Further, the tension between
intra-group cooperation and inter-group competition is consistent with
superorganism properties.


  The majority of the human genome consists of repeated sequences. An important
type of repeated sequences common in the human genome are tandem repeats, where
identical copies appear next to each other. For example, in the sequence
$AGTC\underline{TGTG}C$, $TGTG$ is a tandem repeat, that may be generated from
$AGTCTGC$ by a tandem duplication of length $2$. In this work, we investigate
the possibility of generating a large number of sequences from a \textit{seed},
i.e.\ a small initial string, by tandem duplications of bounded length. We
study the capacity of such a system, a notion that quantifies the system's
generating power. Our results include \textit{exact capacity} values for
certain tandem duplication string systems. In addition, motivated by the role
of DNA sequences in expressing proteins via RNA and the genetic code, we define
the notion of the \textit{expressiveness} of a tandem duplication system as the
capability of expressing arbitrary substrings. We then \textit{completely}
characterize the expressiveness of tandem duplication systems for general
alphabet sizes and duplication lengths. In particular, based on a celebrated
result by Axel Thue from 1906, presenting a construction for ternary
square-free sequences, we show that for alphabets of size 4 or larger, bounded
tandem duplication systems, regardless of the seed and the bound on duplication
length, are not fully expressive, i.e. they cannot generate all strings even as
substrings of other strings. Note that the alphabet of size 4 is of particular
interest as it pertains to the genomic alphabet. Building on this result, we
also show that these systems do not have full capacity. In general, our results
illustrate that duplication lengths play a more significant role than the seed
in generating a large number of sequences for these systems.


  For dynamical systems arising from chemical reaction networks, persistence is
the property that each species concentration remains positively bounded away
from zero, as long as species concentrations were all positive in the
beginning. We describe two graphical procedures for simplifying reaction
networks without breaking known necessary or sufficient conditions for
persistence, by iteratively removing so-called intermediates and catalysts from
the network. The procedures are easy to apply and, in many cases, lead to
highly simplified network structures, such as monomolecular networks. For
specific classes of reaction networks, we show that these conditions for
persistence are equivalent to one another. Furthermore, they can also be
characterized by easily checkable strong connectivity properties of a related
graph. In particular, this is the case for (conservative) monomolecular
networks, as well as cascades of a large class of post-translational
modification systems (of which the MAPK cascade and the $n$-site futile cycle
are prominent examples). Since one of the aforementioned sufficient conditions
for persistence precludes the existence of boundary steady states, our method
also provides a graphical tool to check for that.


  Phylogenetic networks are necessary to represent the tree of life expanded by
edges to represent events such as horizontal gene transfers, hybridizations or
gene flow. Not all species follow the paradigm of vertical inheritance of their
genetic material. While a great deal of research has flourished into the
inference of phylogenetic trees, statistical methods to infer phylogenetic
networks are still limited and under development. The main disadvantage of
existing methods is a lack of scalability. Here, we present a statistical
method to infer phylogenetic networks from multi-locus genetic data in a
pseudolikelihood framework. Our model accounts for incomplete lineage sorting
through the coalescent model, and for horizontal inheritance of genes through
reticulation nodes in the network. Computation of the pseudolikelihood is fast
and simple, and it avoids the burdensome calculation of the full likelihood
which can be intractable with many species. Moreover, estimation at the
quartet-level has the added computational benefit that it is easily
parallelizable. Simulation studies comparing our method to a full likelihood
approach show that our pseudolikelihood approach is much faster without
compromising accuracy. We applied our method to reconstruct the evolutionary
relationships among swordtails and platyfishes ($Xiphophorus$: Poeciliidae),
which is characterized by widespread hybridizations.


  Split networks are a popular tool for the analysis and visualization of
complex evolutionary histories. Every collection of splits (bipartitions) of a
finite set can be represented by a split network. Here we characterize which
collection of splits can be represented using a planar split network. Our main
theorem links these collections of splits with oriented matroids and
arrangements of lines separating points in the plane. As a consequence of our
main theorem, we establish a particularly simple characterization of maximal
collections of these splits.


  Ecosystem approach to fisheries requires a thorough understanding of fishing
impacts on ecosystem status and processes as well as predictive tools such as
ecosystem models to provide useful information for management. The credibility
of such models is essential when used as decision making tools, and model
fitting to observed data is one major criterion to assess such credibility.
However, more attention has been given to the exploration of model behavior
than to a rigorous confrontation to observations, as the calibration of
ecosystem models is challenging in many ways. First, ecosystem models can only
be simulated numerically and are generally too complex for mathematical
analysis and explicit parameter estimation; secondly, the complex dynamics
represented in ecosystem models allow species-specific parameters to impact
other species parameters through ecological interactions; thirdly, critical
data about non-commercial species are often poor; lastly, technical aspects can
be impediments to the calibration with regard to the high computational cost
potentially involved and the scarce documentation published on fitting complex
ecosystem models to data. This work highlights some issues related to the
confrontation of complex ecosystem models to data and proposes a methodology
for a sequential multi-phases calibration of ecosystem models. We propose
criteria to classify the parameters of a model: model dependency and time
variability of the parameters. These criteria and the availability of
approximate initial estimates are used as decision rules to determine which
parameters need to be estimated, and their precedence order in the sequential
calibration process. The end-to-end ecosystem model ROMS-PISCES-OSMOSE applied
to the Northern Humboldt Current Ecosystem is used as an illustrative case
study.


  In this paper an algorithm is given to determine all possible structurally
different linearly conjugate realizations of a given kinetic polynomial system.
The solution is based on the iterative search for constrained dense
realizations using linear programming. Since there might exist exponentially
many different reaction graph structures, we cannot expect to have a
polynomial-time algorithm, but we can organize the computation in such a way
that polynomial time is elapsed between displaying any two consecutive
realizations. The correctness of the algorithm is proved, and possibilities of
a parallel implementation are discussed. The operation of the method is shown
on two illustrative examples.


  Florigens, the leaf-derived signals that initiate flowering, have been
described as mysterious, elusive and the Holy Grail of plant biology.


  DNA G-hairpins are potential key structures participating in folding of human
telomeric guanine quadruplexes (GQ). We examined their properties by standard
MD simulations starting from the folded state and long T-REMD starting from the
unfolded state, accumulating ~130 {\mu}s of atomistic simulations. Antiparallel
G-hairpins should spontaneously form in all stages of the folding to support
lateral and diagonal loops, with sub-{\mu}s scale rearrangements between them.
We found no clear predisposition for direct folding into specific GQ topologies
with specific syn/anti patterns. Our key prediction stemming from the T-REMD is
that an ideal unfolded ensemble of the full GQ sequence populates all 4096
syn/anti combinations of its four G-stretches. The simulations can propose
idealized folding pathways but we explain that such few-state pathways may be
misleading. In the context of the available experimental data, the simulations
strongly suggest that the GQ folding could be best understood by the kinetic
partitioning mechanism with a set of deep competing minima on the folding
landscape, with only a small fraction of molecules directly folding to the
native fold. The landscape should further include nonspecific collapse
processes where the molecules move via diffusion and consecutive random rare
transitions, which could, e.g., structure the propeller loops.


  Even closely related prokaryotes show an astounding diversity in their
ability to grow in different nutritional environments. Mechanistically, this
diversity arises predominantly through horizontal gene transfer, the exchange
of DNA between individuals from different strains. It has been hypothesized
that complex metabolic adaptations--those requiring the acquisition of multiple
distinct DNA segments--can evolve via selectively neutral intermediate steps;
an alternative explanation rests on the existence of intermediate environments
that make each individual DNA acquisition adaptive. However, it is unclear how
important changing environments are compared to neutral explorations of
phenotype space; more fundamentally, it is unknown what fraction of metabolic
adaptations are indeed complex. Here, we use metabolic network simulations to
show that all 3,363 phenotypic innovations observed in the evolutionary history
of 53 E. coli strains arose through the acquisition of a single DNA segment;
while we found no evidence for the contribution of selectively neutral
processes, 10.6% of adaptations to previously unviable environments relied on
the support of DNA acquisitions on earlier phylogenetic branches. 97.0% of all
metabolic phenotypes accessible for the E. coli pan-genome could be bestowed on
any ancestral strain by transferring a single DNA segment from one of the
extant strains. These results demonstrate an amazing ability of the E. coli
lineage to quickly adapt to previously inaccessible environments through a
single DNA acquisition, an ability likely to be mirrored in other clades of
generalist bacteria.


  Over the past 60 years, the Mycobacterium bovis bacille Calmette-Gu\'erin
(BCG) has been used worldwide to prevent tuberculosis (TB). However, BCG has
shown a very variable efficacy in different trials, showing a wide range of
protection in adults against pulmonary TB. Previous studies indicate that this
failure is related to pre-existing immune response to antigens that are common
to environmental sources of mycobacterial antigens and Mycobacterium
tuberculosis. Specifically, two different mechanisms have been hypothesized:
the masking, (previous sensitization confers some level of protection against
TB), and the blocking (previous immune response prevent vaccine taking of a new
TB vaccine), effects. In this work we introduce a series of models to
discriminate between masking and blocking mechanisms and address their relative
likelihood. The application of our models to interpret the results coming from
the BCG-REVAC clinical trials, specifically designed for the study of sources
of efficacy variability yields estimates that are consistent with high levels
of blocking (41% in Manaus -95% C.I. [14%-68%]- and 96% in Salvador -95% C.I.
[52%-100%]-), and no support for masking to play any relevant role in modifying
vaccine efficacy either alone or aside blocking. The quantification of these
effects around a plausible model constitutes a relevant step towards impact
evaluation of novel anti-tuberculosis vaccines, which are susceptible of being
affected by similar effects if applied on individuals previously exposed to
mycobacterial antigens.


  Interest in the study of in-host microbial communities has increased in
recent years due to our improved understanding of the communities' significant
role in host health. As a result, the ability to model these communities using
differential equations, for example, and analyze the results has become
increasingly relevant. The size of the models and limitations in data
collection among many other considerations require that we develop new
parameter estimation methods to address the challenges that arise when using
traditional parameter estimation methods for models of these in-host microbial
communities. In this work, we present the challenges that appear when applying
traditional parameter estimation techniques to differential equation models of
microbial communities, and we provide an original, alternative method to those
techniques. We show the derivation of our method and how our method avoids the
limitations of traditional techniques while including additional benefits. We
also provide simulation studies to demonstrate our method's viability, the
application of our method to a model of intestinal microbial communities to
demonstrate the insights that can be gained from our method, and sample code to
give readers the opportunity to apply our method to their own research.


  Wolbachia is possibly the most studied reproductive parasite of arthropod
species. It appears to be a promising candidate for biocontrol of some mosquito
borne diseases. We begin by developing a sex-structured model for a Wolbachia
infected mosquito population. Our model incorporates the key effects of
Wolbachia infection including cytoplasmic incompatibility and male killing. We
also allow the possibility of reduced reproductive output, incomplete maternal
transmission, and different mortality rates for uninfected/infected male/female
individuals. We study the existence and local stability of equilibria,
including the biologically relevant and interesting boundary equilibria. For
some biologically relevant parameter regimes there may be multiple coexistence
steady states including, very importantly, a coexistence steady state in which
Wolbachia infected individuals dominate. We also extend the model to
incorporate West Nile virus (WNv) dynamics, using an SEI modelling approach.
Recent evidence suggests that a particular strain of Wolbachia infection
significantly reduces WNv replication in Aedes aegypti. We model this via
increased time spent in the WNv-exposed compartment for Wolbachia infected
female mosquitoes. A basic reproduction number $R_0$ is computed for the WNv
infection. Our results suggest that, if the mosquito population consists mainly
of Wolbachia infected individuals, WNv eradication is likely if WNv replication
in Wolbachia infected individuals is sufficiently reduced.


  In spite of the recent interest and advances in linear controllability of
complex networks, controlling nonlinear network dynamics remains to be an
outstanding problem. We develop an experimentally feasible control framework
for nonlinear dynamical networks that exhibit multistability (multiple
coexisting final states or attractors), which are representative of, e.g., gene
regulatory networks (GRNs). The control objective is to apply parameter
perturbation to drive the system from one attractor to another, assuming that
the former is undesired and the latter is desired. To make our framework
practically useful, we consider RESTRICTED parameter perturbation by imposing
the following two constraints: (a) it must be experimentally realizable and (b)
it is applied only temporarily. We introduce the concept of ATTRACTOR NETWORK,
in which the nodes are the distinct attractors of the system, and there is a
directional link from one attractor to another if the system can be driven from
the former to the latter using restricted control perturbation. Introduction of
the attractor network allows us to formulate a controllability framework for
nonlinear dynamical networks: a network is more controllable if the underlying
attractor network is more strongly connected, which can be quantified. We
demonstrate our control framework using examples from various models of
experimental GRNs. A finding is that, due to nonlinearity, noise can
counter-intuitively facilitate control of the network dynamics.


  Biopolymers serve as one-dimensional tracks on which motor proteins move to
perform their biological roles. Motor protein phenomena have inspired
theoretical models of one-dimensional transport, crowding, and jamming.
Experiments studying the motion of Xklp1 motors on reconstituted antiparallel
microtubule overlaps demonstrated that motors recruited to the overlap walk
toward the plus end of individual microtubules and frequently switch between
filaments. We study a model of this system that couples the totally asymmetric
simple exclusion process (TASEP) for motor motion with switches between
antiparallel filaments and binding kinetics. We determine steady-state motor
density profiles for fixed-length overlaps using exact and approximate
solutions of the continuum differential equations and compare to kinetic Monte
Carlo simulations. Overlap motor density profiles and motor trajectories
resemble experimental measurements. The phase diagram of the model is similar
to the single-filament case for low switching rate, while for high switching
rate we find a new low density-high density-low density-high density phase. The
overlap center region, far from the overlap ends, has a constant motor density
as one would naively expect. However, rather than following a simple binding
equilibrium, the center motor density depends on total overlap length, motor
speed, and motor switching rate. The size of the crowded boundary layer near
the overlap ends is also dependent on the overlap length and switching rate in
addition to the motor speed and bulk concentration. The antiparallel
microtubule overlap geometry may offer a previously unrecognized mechanism for
biological regulation of protein concentration and consequent activity.


  Motivation: We introduce TRONCO (TRanslational ONCOlogy), an open-source R
package that implements the state-of-the-art algorithms for the inference of
cancer progression models from (epi)genomic mutational profiles. TRONCO can be
used to extract population-level models describing the trends of accumulation
of alterations in a cohort of cross-sectional samples, e.g., retrieved from
publicly available databases, and individual-level models that reveal the
clonal evolutionary history in single cancer patients, when multiple samples,
e.g., multiple biopsies or single-cell sequencing data, are available. The
resulting models can provide key hints in uncovering the evolutionary
trajectories of cancer, especially for precision medicine or personalized
therapy.
  Availability: TRONCO is released under the GPL license, it is hosted in the
Software section at http://bimib.disco.unimib.it/ and archived also at
bioconductor.org.
  Contact: tronco@disco.unimib.it


  Genetic oscillators, such as circadian clocks, are constantly perturbed by
molecular noise arising from the small number of molecules involved in gene
regulation. One of the strongest sources of stochasticity is the binary noise
that arises from the binding of a regulatory protein to a promoter in the
chromosomal DNA. In this study, we focus on two minimal oscillators based on
activator titration and repressor titration to understand the key parameters
that are important for oscillations and for overcoming binary noise. We show
that the rate of unbinding from the DNA, despite traditionally being considered
a fast parameter, needs to be slow to broaden the space of oscillatory
solutions. The addition of multiple, independent DNA binding sites further
expands the oscillatory parameter space for the repressor-titration oscillator
and lengthens the period of both oscillators. This effect is a combination of
increased effective delay of the unbinding kinetics due to multiple binding
sites and increased promoter ultrasensitivity that is specific for repression.
We then use stochastic simulation to show that multiple binding sites increase
the coherence of the oscillations by mitigating the binary noise. Slow values
of DNA unbinding rate are also effective in alleviating molecular noise due to
the increased distance from the bifurcation point. Our work demonstrates how
DNA binding sites and slow unbinding kinetics, which are often omitted in
biophysical models of gene circuits, can have a significant impact on the
temporal and stochastic dynamics of genetic oscillators.


  Chromatin remodeling machineries are abundant and diverse in eukaryotic
cells. They have been involved in a variety of situations such as histone
exchange and DNA repair, but their importance in gene expression remains
unclear. Although the influence of nucleosome position on the regulation of
gene expression is generally envisioned under the quasi-equilibrium
perspective, it is proposed that given the ATP-dependence of chromatin
remodeling enzymes, certain mechanisms necessitate non-equilibrium treatments.
Examination of the celebrated chromatin remodeling system of the mouse mammary
tumor virus, in which the binding of transcription factors opens the way to
other ones, reveals that breaking equilibrium offers a subtle mode of
transcription factor cooperativity, avoids molecular trapping phenomena and
allows to reconcile previously conflicting experimental data. This mechanism
provides a control lever of promoter responsiveness to transcription factor
combinations, challenging the classical view of the unilateral influence of
pioneer on secondary transcription factors.


  A recently introduced concept of "cooperative equilibrium", based on the
assumption that players have a natural attitude to cooperation, has been proven
a powerful tool in predicting human behaviour in social dilemmas. In this
paper, we extend this idea to more general game models, termed "Pareto
solvable" games, which in particular include the Nash Bargaining Problem and
the Ultimatum Game. We show that games in this class possess a unique pure
cooperative equilibrium. Furthermore, for the Ultimatum Game, this notion
appears to be strongly correlated with a suitably defined variant of the
Dictator Game. We support this observation with the results of a behavioural
experiment conducted using Amazon Mechanical Turk, which demonstrates that our
approach allows for making statistically precise predictions of average
behaviour in such settings.


  Effort has been given for the development of an analytical approach that
helps to address several sigmoidal and non-sigmoidal growth processes found in
literature. In the proposed approach, the role of specific growth rate in
different growth processes has been considered in an unified manner. It is
found that the different growth equations can be derived from the same
functional form of rate equation of specific growth rate. A common functional
form of growth and growth velocity have been derived analytically and it has
been shown that different values of the parameters involved in the description
lead to different growth function. The theta logistic growth can be explained
in this proposed framework. It is found that the competitive environment may
increase the saturation level of population size.


  We perform a geometric study of the equilibrium locus of the flow that models
the diffusion process over a circular network of cells. We prove that when
considering the set of all possible values of the parameters, the equilibrium
locus is a smooth manifold with corners, while for a given value of the
parameters, it is an embedded smooth and connected curve. For different values
of the parameters, the curves are all isomorphic.
  Moreover, we show how to build a homotopy between different curves obtained
for different values of the parameter set. This procedure allows the efficient
computation of the equilibrium point for each value of some first integral of
the system. This point would have been otherwise difficult to be computed for
higher dimensions. We illustrate this construction by some numerical
experiments.
  Eventually, we show that when considering the parameters as inputs, one can
easily bring the system asymptotically to any equilibrium point in the
reachable set, which we also easily characterize.


  Wolbachia is a genus of endosymbiotic bacteria that can infect mosquitoes and
reduce their ability to transmit dengue virus. Although the bacterium is
transmitted vertically from infected mothers to their offspring, it can be
difficult to establish an endemic infection in a wild mosquito population. We
developed and analyzed an ordinary differential equation model to investigate
the transmission dynamics of releasing Wolbachia-infected mosquitoes to
establish an endemic infection in a population of wild uninfected mosquitoes.
Our transmission model for the adult and aquatic-stage mosquitoes takes into
account Wolbachia-induced fitness change and cytoplasmic incompatibility. We
showed that, for a wide range of realistic parameter values, the basic
reproduction number is less than one. Hence, the epidemic will die out if only
a few Wolbachia-infected mosquitoes are introduced into the wild population.
Even though the basic reproduction number is less than one, an endemic
Wolbachia infection can be established if a sufficient number of infected
mosquitoes are released. This threshold effect is created by a backward
bifurcation with three coexisting equilibria: a stable zero-infection
equilibrium, an intermediate-infection unstable endemic equilibrium, and a
high-infection stable endemic equilibrium. We analyzed the impact of reducing
the wild mosquito population before introducing the infected mosquitoes and
observed that the most effective approach to establish the infection in the
wild is based on reducing mosquitoes in both the adult and aquatic stages.


  Ketamine is widely used in veterinary medicine and in medicine. Ketamine is
metabolized to its active metabolite norketamine principally by liver CYP3A.
Drug metabolism alterations during aging have severe consequences particularly
in anesthesiology and very few studies on older animals were conducted for
ketamine. The objective of the present study is to assess the influence of
aging on CYP3A metabolism of ketamine. Liver S9 fractions from 3, 6, 12 and 18
month old male Sprague Dawley rats were prepared and Michaelis-Menten
parameters were determined for primary metabolic pathways. The derived maximum
enzyme velocity (i.e. Vmax) suggests a rapid saturation of the CYP3A enzyme
active sites in liver S9 fractions of 18-month old rats. Observed Vmax for
Liver S9 fractions from 3, 6 and 12 month old male Sprague Dawley rats were
2.39 (+-0.23), 2.61 (+-0.18), and 2.07 (+-0.07) respectively compared to 0.68
(+-0.02) for Liver S9 fractions from 18 month old male Sprague Dawley rats.
Interestingly, we observed a 6 to 7 fold change in the derived Km when
comparing Liver S9 fractions from 18 month old male Sprague Dawley rats with
Liver S9 fractions from younger rats. Our results suggest that rat CYP3A enzyme
undergoes conformational changes with age particularly in our geriatric group
(e.g. 18 month rats) leading significant decrease in the rate of formation of
norketamine. Moreover, our results strongly suggest a severe impairment of
CYP3A ketamine mediated metabolism.


  Molecular noise in gene regulatory networks has two intrinsic components, one
part being due to fluctuations caused by the birth and death of protein or mRNA
molecules which are often present in small numbers and the other part arising
from gene state switching, a single molecule event. Stochastic dynamics of gene
regulatory circuits appears to be largely responsible for bifurcations into a
set of multi-attractor states that encode different cell phenotypes. The
interplay of dichotomous single molecule gene noise with the nonlinear
architecture of genetic networks generates rich and complex phenomena. In this
paper we elaborate on an approximate framework that leads to simple hybrid
multi-scale schemes well suited for the quantitative exploration of the steady
state properties of large-scale cellular genetic circuits.Through a path sum
based analysis of trajectory statistics we elucidate the connection of these
hybrid schemes to the underlying master equation and provide a rigorous
justification for using dichotomous noise based models to study genetic
networks. Numerical simulations of circuit models reveal that the contribution
of the genetic noise of single molecule origin to the total noise is
significant for a wide range of kinetic regimes.


  The genomic evolution inherent to cancer relates directly to a renewed focus
on the voluminous next generation sequencing (NGS) data, and machine learning
for the inference of explanatory models of how the (epi)genomic events are
choreographed in cancer initiation and development. However, despite the
increasing availability of multiple additional -omics data, this quest has been
frustrated by various theoretical and technical hurdles, mostly stemming from
the dramatic heterogeneity of the disease. In this paper, we build on our
recent works on "selective advantage" relation among driver mutations in cancer
progression and investigate its applicability to the modeling problem at the
population level. Here, we introduce PiCnIc (Pipeline for Cancer Inference), a
versatile, modular and customizable pipeline to extract ensemble-level
progression models from cross-sectional sequenced cancer genomes. The pipeline
has many translational implications as it combines state-of-the-art techniques
for sample stratification, driver selection, identification of
fitness-equivalent exclusive alterations and progression model inference. We
demonstrate PiCnIc's ability to reproduce much of the current knowledge on
colorectal cancer progression, as well as to suggest novel experimentally
verifiable hypotheses.


  We consider the problem of jointly estimating multiple inverse covariance
matrices from high-dimensional data consisting of distinct classes. An
$\ell_2$-penalized maximum likelihood approach is employed. The suggested
approach is flexible and generic, incorporating several other
$\ell_2$-penalized estimators as special cases. In addition, the approach
allows specification of target matrices through which prior knowledge may be
incorporated and which can stabilize the estimation procedure in
high-dimensional settings. The result is a targeted fused ridge estimator that
is of use when the precision matrices of the constituent classes are believed
to chiefly share the same structure while potentially differing in a number of
locations of interest. It has many applications in (multi)factorial study
designs. We focus on the graphical interpretation of precision matrices with
the proposed estimator then serving as a basis for integrative or meta-analytic
Gaussian graphical modeling. Situations are considered in which the classes are
defined by data sets and subtypes of diseases. The performance of the proposed
estimator in the graphical modeling setting is assessed through extensive
simulation experiments. Its practical usability is illustrated by the
differential network modeling of 12 large-scale gene expression data sets of
diffuse large B-cell lymphoma subtypes. The estimator and its related
procedures are incorporated into the R-package rags2ridges.


  Many cells contain non-centrosomal arrays of microtubules (MT), but the
assembly, organisation and function of these arrays are poorly understood. We
present the first theoretical model for the non-centrosomal MT cytoskeleton in
$Drosophila$ oocytes, in which $bicoid$ and $oskar$ mRNAs become localised to
establish the anterior-posterior body axis. Constrained by experimental
measurements, the model shows that a simple gradient of cortical MT nucleation
is sufficient to reproduce the observed MT distribution, cytoplasmic flow
patterns and localisation of $oskar$ and naive $bicoid$ mRNAs. Our simulations
exclude a major role for cytoplasmic flows in localisation and reveal an
organisation of the MT cytoskeleton that is more ordered than previously
thought. Furthermore, modulating cortical MT nucleation induces a bifurcation
in cytoskeletal organisation that accounts for the phenotypes of polarity
mutants. Thus, our three-dimensional model explains many features of the MT
network and highlights the importance of differential cortical MT nucleation
for axis formation.


  We investigate a densely packed, non-random arrangement of forty-six
chromosomes (46,XY) in human nuclei. Here, we model systems-level chromosomal
crosstalk by unifying intrinsic parameters (chromosomal length and number of
genes) across all pairs of chromosomes in the genome to derive an extrinsic
parameter called effective gene density. The hierarchical clustering and
underlying degeneracy in the effective gene density space reveal systems-level
constraints for spatial arrangement of clusters of chromosomes that were
previously unknown. Our findings corroborate experimental data on spatial
chromosomal arrangement in human nuclei, from fibroblast and lymphocyte cell
lines, thereby establishing that human genome constrains chromosomal
arrangement. We propose that this unified theory, which requires no additional
experimental input, may be extended to other eukaryotic species with annotated
genomes to infer their constrained self-organized spatial arrangement of
chromosomes.


  Species spanning the animal kingdom have evolved extravagant and costly
ornaments to attract mating partners. Zahavi's handicap principle offers an
elegant explanation for this: ornaments signal individual quality, and must be
costly to ensure honest signalling, making mate selection more efficient. Here
we incorporate the assumptions of the handicap principle into a mathematical
model and show that they are sufficient to explain the heretofore puzzling
observation of bimodally distributed ornament sizes in a variety of species.


  Rising interest in the resilience of ecological systems has spawned diverse
interpretations of the term's precise meaning. This paper classifies and
explores definitions of resilience from the ecological literature using a
dynamical systems framework. A model consisting of ordinary differential
equations is assumed to represent the ecological system. The question
"resilience of what to what?" posed by Carpenter et al. [2001] informs two
broad categories of definitions, based on resilience to state variable
perturbations and to parameter changes, respectively. Definitions of resilience
to state variable perturbations include measures of basin size (relevant to
one-time perturbations) and of basin steepness (relevant to repeated
perturbations). Resilience to parameter changes can be quantified by viewing
parameters as state variables, but could also take into account reversibility.
The dynamical systems viewpoint yields the following key insights: (i) system
properties that confer resistance to state variable perturbations can differ
from those that promote recovery from them, (ii) recovery rates, oft deemed the
purview of "engineering resilience," do matter for ecological resilience, and
(iii) resilience to repeated state variable perturbations correlates with
resilience to parameter changes due to the critical slowing down phenomenon.


  The study of network structure has uncovered signatures of the organization
of complex systems. However, there is also a need to understand how to control
them; for example, identifying strategies to revert a diseased cell to a
healthy state, or a mature cell to a pluripotent state. Two recent
methodologies suggest that the controllability of complex systems can be
predicted solely from the graph of interactions between variables, without
considering their dynamics: structural controllability and minimum dominating
sets. We demonstrate that such structure-only methods fail to characterize
controllability when dynamics are introduced. We study Boolean network
ensembles of network motifs as well as three models of biochemical regulation:
the segment polarity network in Drosophila melanogaster, the cell cycle of
budding yeast Saccharomyces cerevisiae, and the floral organ arrangement in
Arabidopsis thaliana. We demonstrate that structure-only methods both
undershoot and overshoot the number and which sets of critical variables best
control the dynamics of these models, highlighting the importance of the actual
system dynamics in determining control. Our analysis further shows that the
logic of automata transition functions, namely how canalizing they are, plays
an important role in the extent to which structure predicts dynamics.


  Network alignment (NA) aims to find regions of similarities between molecular
networks of different species. There exist two NA categories: local (LNA) or
global (GNA). LNA finds small highly conserved network regions and produces a
many-to-many node mapping. GNA finds large conserved regions and produces a
one-to-one node mapping. Given the different outputs of LNA and GNA, when a new
NA method is proposed, it is compared against existing methods from the same
category. However, both NA categories have the same goal: to allow for
transferring functional knowledge from well- to poorly-studied species between
conserved network regions. So, which one to choose, LNA or GNA? To answer this,
we introduce the first systematic evaluation of the two NA categories.
  We introduce new measures of alignment quality that allow for fair comparison
of the different LNA and GNA outputs, as such measures do not exist. We provide
user-friendly software for efficient alignment evaluation that implements the
new and existing measures. We evaluate prominent LNA and GNA methods on
synthetic and real-world biological networks. We study the effect on alignment
quality of using different interaction types and confidence levels. We find
that the superiority of one NA category over the other is context-dependent.
Further, when we contrast LNA and GNA in the application of learning novel
protein functional knowledge, the two produce very different predictions,
indicating their complementarity. Our results and software provide guidelines
for future NA method development and evaluation.


  The biological properties of proteins are uniquely determined by their
structure and dynamics. A protein in solution populates a structural ensemble
of metastable configurations around the global fold. From overall rotation to
local fluctuations, the dynamics of proteins can cover several orders of
magnitude in time scales. We propose a simulation-free coarse-grained approach
which utilizes knowledge of the important metastable folded states of the
protein to predict the protein dynamics. This approach is based upon the
Langevin Equation for Protein Dynamics (LE4PD), a Langevin formalism in the
coordinates of the protein backbone. The linear modes of this Langevin
formalism organize the fluctuations of the protein, so that more extended
dynamical cooperativity relates to increasing energy barriers to mode
diffusion. The accuracy of the LE4PD is verified by analyzing the predicted
dynamics across a set of seven different proteins for which both relaxation
data and NMR solution structures are available. Using experimental NMR
conformers as the input structural ensembles, LE4PD predicts quantitatively
accurate results, with correlation coefficient \r{ho} = .93 to NMR backbone
relaxation measurements for the seven proteins. The NMR solution structure
derived ensemble and predicted dynamical relaxation is compared with molecular
dy- namics simulation-derived structural ensembles and LE4PD predictions, and
are consistent in the timescale of the simulations. The use of the experimental
NMR conformers frees the approach from computationally demanding simulations.


  A biological system consists of a variety of polymers that are synthesized
from monomers, by catalysis that exists only for some long polymers. It is
important to elucidate the emergence and sustenance of such autocatalytic
polymerization. We analyze here the stochastic polymerization reaction
dynamics, to investigate the transition time from a state with almost no
catalysts to a state with sufficient catalysts. We found an optimal volume that
minimizes this transition time, which agrees with the inverse of the catalyst
concentration at the unstable fixed point that separates the two states, as is
theoretically explained. Relevance to the origin of life is also discussed.


  A multi-patch and multi-group modeling framework describing the dynamics of a
class of diseases driven by the interactions between vectors and hosts
structured by groups is formulated. Hosts' dispersal is modeled in terms of
patch-residence times with the nonlinear dynamics taking into account the
\textit{effective} patch-host size. The residence times basic reproduction
number $\mathcal R_0$ is computed and shown to depend on the relative
environmental risk of infection. The model is robust, that is, the disease free
equilibrium is globally asymptotically stable (GAS) if $\mathcal R_0\leq1$ and
a unique interior endemic equilibrium is shown to exist that is GAS whenever
$\mathcal R_0>1$ whenever the configuration of host-vector interactions is
irreducible. The effects of \textit{patchiness} and \textit{groupness}, a
measure of host-vector heterogeneous structure, on the basic reproduction
number $\mathcal R_0$, are explored. Numerical simulations are carried out to
highlight the effects of residence times on disease prevalence.


  The biological function of proteins is encoded in their structure and
expressed through the mediation of their dynamics. Local fluctuations are known
to initiate biologically relevant pathways as they cooperatively enhance the
dynamics in specific regions in the protein. Those biologically active regions
provide energetically-comparable conformational states that can be trapped by a
reacting partner. We analyze this mechanism as we calculate the dynamics of
monomeric and dimerized HIV protease, and free Insulin Growth Factor II
Receptor (IGF2R) domain 11 and its IGF2R:IGF2 complex. We adopt a newly
developed coarse-grained model, the Langevin Equation for Protein Dynamics
(LE4PD), which predicts dynamical relevant mechanisms with high accuracy. Both
simulation-derived and experimental NMR conformers are the input structural
ensembles for the LE4PD. The use of the experimental NMR conformers requires
minimal computational resources.


  We consider a multi-species community modelled as a complex network of
populations, where the links are given by a random asymmetric matrix J, with
fraction 1-C of zero entries, where C reflects the over-all connectivity of the
system. The non-zero elements of J are drawn from a gaussian distribution with
mean 'mu' and standard deviation . The signs of the elements J reflect the
nature of density-dependent interactions, such as predatory-prey, mutualism or
competition, and their magnitudes reflect the strength of the interaction. In
this study we try to uncover the broad features of the interspecies
interactions that determine the global robustness of this network, as indicated
by the average number of active nodes (i.e. non-extinct species) in the
network, and the total population, reflecting the biomass yield. We find that
the network transitions from a completely extinct system to one where all nodes
are active, as the mean interaction strength goes from negative to positive. We
also find that the total population, displays distinct non-monotonic scaling
behaviour with respect to the product 'mu'C, implying that survival is
dependent not merely on the number of links, but rather on the combination of
the sparseness of the connectivity matrix and the net interaction strength.
Moreover, the total population levels are optimal when the network has
intermediate net positive connection strengths. At the local level we observe
marked qualitative changes in dynamical patterns, ranging from anti-phase
clusters of period 2 cycles and chaotic bands, to fixed points, under the
variation of mean 'mu' of the interaction strengths. Lastly, we propose an
effective low dimensional map to capture the behavior of the entire network,
and this provides a broad understanding of the interplay of the local dynamical
patterns and the global robustness trends in the network.


  The stochastic nature of chemical reactions involving randomly fluctuating
population sizes has lead to a growing research interest in discrete-state
stochastic models and their analysis. A widely-used approach is the description
of the temporal evolution of the system in terms of a chemical master equation
(CME). In this paper we study two approaches for approximating the underlying
probability distributions of the CME. The first approach is based on an
integration of the statistical moments and the reconstruction of the
distribution based on the maximum entropy principle. The second approach relies
on an analytical approximation of the probability distribution of the CME using
the system size expansion, considering higher-order terms than the linear noise
approximation. We consider gene expression networks with unimodal and
multimodal protein distributions to compare the accuracy of the two approaches.
We find that both methods provide accurate approximations to the distributions
of the CME while having different benefits and limitations in applications.


  Gene products (RNAs, proteins) often occur at low molecular counts inside
individual cells, and hence are subject to considerable random fluctuations
(noise) in copy number over time. Not surprisingly, cells encode diverse
regulatory mechanisms to buffer noise. One such mechanism is the incoherent
feedforward circuit. We analyze a simplistic version of this circuit, where an
upstream regulator X affects both the production and degradation of a protein
Y. Thus, any random increase in X's copy numbers would increase both production
and degradation, keeping Y levels unchanged. To study its stochastic dynamics,
we formulate this network into a mathematical model using the Chemical Master
Equation formulation. We prove that if the functional dependence of Y's
production and degradation on X is similar, then the steady-distribution of Y's
copy numbers is independent of X. To investigate how fluctuations in Y
propagate downstream, a protein Z whose production rate only depend on Y is
introduced. Intriguingly, results show that the extent of noise in Z increases
with noise in X, in spite of the fact that the magnitude of noise in Y is
invariant of X. Such counter intuitive results arise because X enhances the
time-scale of fluctuations in Y, which amplifies fluctuations in downstream
processes. In summary, while feedforward systems can buffer a protein from
noise in its upstream regulators, noise can propagate downstream due to changes
in the time-scale of fluctuations.


  Time series data on cropping pattern at disaggregated level were analysed and
its implications on geospatial drought assessment were demonstrated. An index
of Cropping Pattern Dissimilarity (CP-DI) between a pair of years, developed in
this study, proved that the cropping pattern of a year has a higher degree of
similarity with that of recent past years only and tends to be dissimilar with
longer time difference. The temporal divergence in cropping pattern has direct
implications on geospatial approach of drought assessment, in which, time
series NDVI data are compared for drought interpretation. It was found that,
seasonal NDVI profiles of drought year and normal year did not show any anomaly
when the cropping patterns were dissimilar and two normal years having
dissimilar cropping pattern showed different NDVI profiles. Therefore, it is
suggested that such temporal comparisons of NDVI are better restricted to
recent past years to achieve more objective interpretation.


  This work gives a mathematical foundation for bifurcation from a stable
equilibrium in the genome. We construct idealized dynamics associated with the
genome. For this dynamics we investigate the two main bifurcations from a
stable equilibrium. Finally, we give a mathematical framework for the
Repressilator and the Toggle gene circuits.


  As one of the most popular statistical and machine learning models, logistic
regression with regularization has found wide adoption in biomedicine, social
sciences, information technology, and so on. These domains often involve data
of human subjects that are contingent upon strict privacy regulations.
Increasing concerns over data privacy make it more and more difficult to
coordinate and conduct large-scale collaborative studies, which typically rely
on cross-institution data sharing and joint analysis. Our work here focuses on
safeguarding regularized logistic regression, a widely-used machine learning
model in various disciplines while at the same time has not been investigated
from a data security and privacy perspective. We consider a common use scenario
of multi-institution collaborative studies, such as in the form of research
consortia or networks as widely seen in genetics, epidemiology, social
sciences, etc. To make our privacy-enhancing solution practical, we demonstrate
a non-conventional and computationally efficient method leveraging distributing
computing and strong cryptography to provide comprehensive protection over
individual-level and summary data. Extensive empirical evaluation on several
studies validated the privacy guarantees, efficiency and scalability of our
proposal. We also discuss the practical implications of our solution for
large-scale studies and applications from various disciplines, including
genetic and biomedical studies, smart grid, network analysis, etc.


  Phylogenetic networks are mathematical structures for modeling and
visualization of reticulation processes in the study of evolution. Galled
networks, reticulation visible networks, nearly-stable networks and
stable-child networks are the four classes of phylogenetic networks that are
recently introduced to study the topological and algorithmic aspects of
phylogenetic networks. We prove the following results.
  (1) A binary galled network with n leaves has at most 2(n-1) reticulation
nodes. (2) A binary nearly-stable network with n leaves has at most 3(n-1)
reticulation nodes. (3) A binary stable-child network with n leaves has at most
7(n-1) reticulation nodes.


  Environmental stress, such as oxidative or heat stress, induces the
activation of the Heat Shock Response (HSR) which leads to an increase in the
heat shock proteins (HSPs) level. These HSPs act as molecular chaperones to
maintain proteostasis. Even if the main heat shock response partners are well
known, a detailed description of the dynamical properties of the HSR network is
still missing. In this study, we derive a minimal mathematical model of
cellular response to heat shock that reproduces available experimental data
sets both on transcription factor activity and cell viability. This simplistic
model highlights the key mechanistic processes that rule the HSR network and
reveals (i) the titration of Heat Shock Factor 1 (HSF1) by chaperones as the
guiding line of the network, (ii) that protein triage governs the fate of
damaged proteins and (iii) three different temperature regimes describing
normal, acute or chronic stress.


  Two species with similar resource requirements respond in a characteristic
way to variations in their habitat -- their abundances rise and fall in
concert. We use this idea to learn how bacterial populations in the microbiota
respond to habitat conditions that vary from person-to-person across the human
population. Our mathematical framework shows that habitat fluctuations are
sufficient for explaining intra-bodysite correlations in relative species
abundances from the Human Microbiome Project. We explicitly show that the
relative abundances of phylogenetically related species are positively
correlated and can be predicted from taxonomic relationships. We identify a
small set of functional pathways related to metabolism and maintenance of the
cell wall that form the basis of a common resource sharing niche space of the
human microbiota.


  A modified ssDNA SELEX protocol was developed in order to evolve a randomized
library of imidazole modified ssDNA sequences towards sequences that mediate
the formation of gold nanoparticles from Au3+ precursor ions in aqueous
solution. Active sequences bound to nanoparticles were partitioned from
inactive sequences based on density via ultracentrifugation through a
discontinuous sucrose gradient. Colloidal gold solutions produced by the
evolved pool had a distinct absorbance spectra and produced nanoparticles with
a narrower distribution of sizes compared to colloidal gold solutions produced
by the starting, randomized pool of imidazole modified ssDNA. Sequencing data
from the evolved pool shows that conserved 5 and 6 nt motifs were shared
amongst many of the isolates, which indicates that these motifs could serve as
chelation sites for gold atoms or help stabilize colloidal gold solutions in a
base specific manner.


  Numerous experimental studies have demonstrated that the microenvironment is
a key regulator influencing the proliferative and migrative potentials of
species. Spatial and temporal disturbances lead to adverse and hazardous
microenvironments for cellular systems that is reflected in the phenotypic
heterogeneity within the system. In this paper, we study the effect of
microenvironment on the invasive capability of species, or mutants, on
structured grids under the influence of site-dependent random proliferation in
addition to a migration potential. We discuss both continuous and discrete
fitness distributions. Our results suggest that the invasion probability is
negatively correlated with the variance of fitness distribution of mutants (for
both advantageous and neutral mutants) in the absence of migration of both
types of cells. A similar behaviour is observed even in the presence of a
random fitness distribution of host cells in the system with neutral fitness
rate. In the case of a bimodal distribution, we observe zero invasion
probability until the system reaches a (specific) proportion of advantageous
phenotypes. Also, we find that the migrative potential amplifies the invasion
probability as the variance of fitness of mutants increases in the system,
which is the exact opposite in the absence of migration. Our computational
framework captures the harsh microenvironmental conditions through quenched
random fitness distributions and migration of cells, and our analysis shows
that they play an important role in the invasion dynamics of several biological
systems such as bacterial micro-habitats, epithelial dysplasia, and metastasis.
We believe that our results may lead to more experimental studies, which can in
turn provide further insights into the role and impact of heterogeneous
environments on invasion dynamics.


  We propose and study a model of hypothetical magnetosensitive ionic channels
which are long thought to be a possible candidate to explain the influence of
weak magnetic fields on living organisms ranging from magnetotactic bacteria to
fishes, birds, rats, bats and other mammals including humans. The core of the
model is provided by a short chain of magnetosomes serving as a sensor which is
coupled by elastic linkers to the gating elements of ion channels forming a
small cluster in the cell membrane. The magnetic sensor is fixed by one end on
cytoskeleton elements attached to the membrane and is exposed to viscoelastic
cytosol. Its free end can reorient stochastically and subdiffusively in
viscoelastic cytosol responding to external magnetic field changes and open the
gates of coupled ion channels. The sensor dynamics is generally bistable due to
bistability of the gates which can be in two states with probabilities which
depend on the sensor orientation. For realistic parameters, it is shown that
this model channel can operate in the magnetic field of Earth for a small
number (5 to 7) of single-domain magnetosomes constituting the sensor rod each
of which has a typical size found in magnetotactic bacteria and other
organisms, or even just one sufficiently large nanoparticle of a characteristic
size also found in nature. It is shown that due to viscoelasticity of medium
the bistable gating dynamics generally exhibits power law and stretched
exponential distributions of the residence times of the channels in their open
and closed states. This provides a generic physical mechanism for explanation
of the origin of such anomalous kinetics for other ionic channels whose sensors
move in viscoelastic environment provided by either cytosol or biological
membrane, in a quite general context, beyond the fascinating hypothesis of
magnetosensitive ionic channels we explore.


  In the absence of convincing evidence, data for Sweden and Mauritius are used
in academic publications to illustrate the Demographic Transition Theory. These
data are closely examined and found to be in clear contradiction of this
theory. Demographic Transition Theory is also contradicted by the best
available data for England. Other examples of contradicting evidence are also
discussed.


  Influenza A is a serious disease that causes significant morbidity and
mortality, and vaccines against the seasonal influenza disease are of variable
effectiveness. In this paper, we discuss use of the $p_{\rm epitope}$ method to
predict the dominant influenza strain and the expected vaccine effectiveness in
the coming flu season. We illustrate how the effectiveness of the 2014/2015
A/Texas/50/2012 [clade 3C.1] vaccine against the A/California/02/2014 [clade
3C.3a] strain that emerged in the population can be estimated via pepitope. In
addition, we show by a multidimensional scaling analysis of data collected
through 2014, the emergence of a new A/New Mexico/11/2014-like cluster [clade
3C.2a] that is immunologically distinct from the A/California/02/2014-like
strains.


  Atrial fibrillation (AF) is a leading cause of morbidity and mortality. AF
prevalence increases with age, which is attributed to pathophysiological
changes that aid AF initiation and perpetuation. Current state-of-the-art
models are only capable of simulating short periods of atrial activity at high
spatial resolution, whilst the majority of clinical recordings are based on
infrequent temporal datasets of limited spatial resolution. Being able to
estimate disease progression informed by both modelling and clinical data would
be of significant interest. In addition an analysis of the temporal
distribution of recorded fibrillation episodes AF density can provide insights
into recurrence patterns. We present an initial analysis of the AF density
measure using a simplified idealised stochastic model of a binary time series
representing AF episodes. The future aim of this work is to develop robust
clinical measures of progression which will be tested on models that generate
long-term synthetic data. These measures would then be of clinical interest in
deciding treatment strategies.


  An important occurrence in many cellular contexts is the crossing of a
prescribed threshold by a regulatory protein. The timing of such events is
stochastic as a consequence of the innate randomness in gene expression. A
question of interest is to understand how gene expression is regulated to
achieve precision in event timing. To address this, we model event timing using
the first-passage time framework - a mathe- matical tool to analyze the time
when a stochastic process first crosses a specific threshold. The protein
evolution is described via a simple stochastic model of gene expression.
Moreover, we consider the feedback regulation of protein degradation to be a
possible noise control mechanism employed to achieve the precision. Exact
analytical formulas are developed for the distribution and moments of the
first-passage time. Using these expressions, we investigate for the optimal
feedback strategy such that noise (coefficient of variation squared) in event
timing is minimized around a given fixed mean time. Our results show that the
minimum noise is achieved when the protein degradation rate is zero for all
protein levels. Lastly, the implications of this finding are discussed.


  Many fundamental cellular processes involve small numbers of molecules. When
numbers are small, fluctuations dominate, and stochastic models, which account
for these fluctuations, are required. In this chapter, we describe minimal
stochastic models of three fundamental cellular processes: gene expression,
protein modification, and polymerization. We introduce key analytic tools for
solving each model, including the generating function, eigenfunction expansion,
and operator methods, and we discuss how these tools are extended to more
complicated models. These analytic tools provide an elegant, efficient, and
often insightful alternative to stochastic simulation.


  We present a tool, keep me around (kma), a suite of python scripts and an R
package that finds retained introns in RNA-Seq experiments and incorporates
biological replicates to reduce the number of false positives when detecting
retention events. kma uses the results of existing quantification tools that
probabilistically assign multi-mapping reads, thus interfacing easily with
transcript quantification pipelines. The data is represented in a convenient,
database style format that allows for easy aggregation across introns, genes,
samples, and conditions to allow for further exploratory analysis.


  Consumer foraging behaviors are dynamic, changing in response to prey
availability, seasonality, competition, and even the consumer's physiological
state. The isotopic composition of a consumer is a product of these factors as
well as the isotopic 'landscape' of its prey, i.e. the isotopic mixing space.
Here we build a mechanistic framework that links the ecological and
physiological processes of an individual consumer to the isotopic distribution
that describes its diet, and ultimately to the isotopic composition of its own
tissues, defined as its 'isotopic niche'. By coupling these processes, we
systematically investigate under what conditions the isotopic niche of a
consumer changes as a function of both the geometric properties of its mixing
space and foraging strategies that may be static or dynamic over time. Results
of our derivations reveal general insight into the conditions impacting
isotopic niche width as a function of consumer specialization on prey, as well
as the consumer's ability to transition between diets over time. We show
analytically that moderate specialization on isotopically unique prey can serve
to maximize a consumer's isotopic niche width, while temporally dynamic diets
will tend to result in peak isotopic variance during dietary transitions. We
demonstrate the relevance of our theoretical findings by examining a marine
system composed of nine invertebrate species commonly consumed by sea otters.
In general, our analytical framework highlights the complex interplay of mixing
space geometry and consumer dietary behavior in driving expansion and
contraction of the isotopic niche. Because this approach is established on
ecological mechanism, it is well-suited for enhancing the ecological
interpretation, and uncovering the root causes, of observed isotopic data.


  Phylodynamics seeks to estimate effective population size fluctuations from
molecular sequences of individuals sampled from a population of interest. One
way to accomplish this task formulates an observed sequence data likelihood
exploiting a coalescent model for the sampled individuals' genealogy and then
integrating over all possible genealogies via Monte Carlo or, less efficiently,
by conditioning on one genealogy estimated from the sequence data. However,
when analyzing sequences sampled serially through time, current methods
implicitly assume either that sampling times are fixed deterministically by the
data collection protocol or that their distribution does not depend on the size
of the population. Through simulation, we first show that, when sampling times
do probabilistically depend on effective population size, estimation methods
may be systematically biased. To correct for this deficiency, we propose a new
model that explicitly accounts for preferential sampling by modeling the
sampling times as an inhomogeneous Poisson process dependent on effective
population size. We demonstrate that in the presence of preferential sampling
our new model not only reduces bias, but also improves estimation precision.
Finally, we compare the performance of the currently used phylodynamic methods
with our proposed model through clinically-relevant, seasonal human influenza
examples.


  Background: Synthetic lethality (SL) refers to the genetic interaction
between two or more genes where only their co-alteration (e.g. by mutations,
amplifications or deletions) results in cell death. In recent years, SL has
emerged as an attractive therapeutic strategy against cancer: by targeting the
SL partners of altered genes in cancer cells, these cells can be selectively
killed while sparing the normal cells. Consequently, a number of studies have
attempted prediction of SL interactions in human, a majority by extrapolating
SL interactions inferred through large-scale screens in model organisms.
However, these predicted SL interactions either do not hold in human cells or
do not include genes that are (frequently) altered in human cancers, and are
therefore not attractive in the context of cancer therapy.
  Results: Here, we develop a computational approach to infer SL interactions
directly from frequently altered genes in human cancers. It is based on the
observation that pairs of genes that are altered in a (significantly) mutually
exclusive manner in cancers are likely to constitute lethal combinations. Using
genomic copy-number and gene-expression data from four cancers, breast,
prostate, ovarian and uterine (total 3980 samples) from The Cancer Genome
Atlas, we identify 718 genes that are frequently amplified or upregulated, and
are likely to be synthetic lethal with six key DNA-damage response (DDR) genes
in these cancers. By comparing with published data on gene essentiality (~16000
genes) from ten DDR-deficient cancer cell lines, we show that our identified
genes are enriched among the top quartile of essential genes in these cell
lines, implying that our inferred genes are highly likely to be (synthetic)
lethal upon knockdown in these cell lines.


  Populations exhibiting partial migration consist of two groups of
individuals: Those that mi- grate between habitats, and those that remain fixed
in a single habitat. We propose several discrete-time population models to
investigate the coexistence of migrants and residents. The first class of
models is linear, and we distinguish two scenarios. In the first, there is a
single egg pool to which both populations contribute. A fraction of the eggs is
destined to become migrants, and the remainder become residents. In a second
model, there are two distinct egg pools to which the two types contribute, one
corresponding to residents and another to migrants. The asymptotic growth or
decline in these models can be phrased in terms of the value of the basic
reproduction number being larger or less than one respectively. A second class
of models incorporates density dependence effects. It is assumed that increased
densities in the various life history stages adversely affect the success of
transitioning of individuals to subsequent stages. Here too we consider models
with one or two egg pools. Although these are nonlinear models, their
asymptotic dynamics can still be classified in terms of the value of a locally
defined basic reproduction number: If it is less than one, then the entire
population goes extinct, whereas it settles at a unique fixed point consisting
of a mixture of residents and migrants, when it is larger than one. Thus, the
value of the basic reproduction number can be used to predict the stable
coexistence or collapse of populations exhibiting partial migration.


  Protein function and dynamics are closely related to its sequence and
structure. However prediction of protein function and dynamics from its
sequence and structure is still a fundamental challenge in molecular biology.
Protein classification, which is typically done through measuring the
similarity be- tween proteins based on protein sequence or physical
information, serves as a crucial step toward the understanding of protein
function and dynamics. Persistent homology is a new branch of algebraic
topology that has found its success in the topological data analysis in a
variety of disciplines, including molecular biology. The present work explores
the potential of using persistent homology as an indepen- dent tool for protein
classification. To this end, we propose a molecular topological fingerprint
based support vector machine (MTF-SVM) classifier. Specifically, we construct
machine learning feature vectors solely from protein topological fingerprints,
which are topological invariants generated during the filtration process. To
validate the present MTF-SVM approach, we consider four types of problems.
First, we study protein-drug binding by using the M2 channel protein of
influenza A virus. We achieve 96% accuracy in discriminating drug bound and
unbound M2 channels. Additionally, we examine the use of MTF-SVM for the
classification of hemoglobin molecules in their relaxed and taut forms and
obtain about 80% accuracy. The identification of all alpha, all beta, and
alpha-beta protein domains is carried out in our next study using 900 proteins.
We have found a 85% success in this identifica- tion. Finally, we apply the
present technique to 55 classification tasks of protein superfamilies over 1357
samples. An average accuracy of 82% is attained. The present study establishes
computational topology as an independent and effective alternative for protein
classification.


  Data describing the growth of the world population in the past 12,000 years
are analysed. It is shown that, if unchecked, population does not increase
exponentially but hyperbolically. This analysis reveals three
approximately-determined episodes of hyperbolic growth: 10,000-500 BC, AD
500-1200 and AD 1400-1950, representing a total of about 89% of the past 12,000
years. It also reveals three demographic transitions: 500 BC-AD 500, AD
1200-1400 and AD 1950-present, representing the remaining 11% of the past
12,000 years. The first two transitions were between sustained hyperbolic
trajectories. The current transition is to an unknown trajectory. There was
never any form of dramatic transition from stagnation to growth, described
often as a takeoff, because there was no stagnation in the growth of the world
population. Correct understanding of the historical growth of human population
is essential in the correct interpretation of the historical growth of income
per capita.


  Living organisms are non-equilibrium, fluctuating, dynamic systems containing
multi-step biological signaling cascades (BSC) with the ability to transduce
changes in the concentration of extracellular molecules such as cytokines into
changes in gene expression. Here, we derived basic equations that describe the
channel capacity and information density of BSC in terms of the average entropy
production rate deduced using the fluctuation theorem.


  Human seasonal influenza viruses evolve rapidly, enabling the virus
population to evade immunity and re-infect previously infected individuals.
Antigenic properties are largely determined by the surface glycoprotein
hemagglutinin (HA) and amino acid substitutions at exposed epitope sites in HA
mediate loss of recognition by antibodies. Here, we show that antigenic
differences measured through serological assay data are well described by a sum
of antigenic changes along the path connecting viruses in a phylogenetic tree.
This mapping onto the tree allows prediction of antigenicity from HA sequence
data alone. The mapping can further be used to make predictions about the
makeup of the future seasonal influenza virus population, and we compare
predictions between models with serological and sequence data. To make timely
model output readily available, we developed a web browser based application
that visualizes antigenic data on a continuously updated phylogeny.


  The field of information science has greatly developed, and applications in
various fields have emerged. In this paper, we evaluated the coding system in
the theory of Tsallis entropy for transmission of messages and aimed to
formulate the channel capacity by maximization of the Tsallis entropy within a
given condition of code length. As a result, we obtained a simple relational
expression between code length and code appearance probability and,
additionally, a generalized formula of the channel capacity on the basis of
Tsallis entropy statistics. This theoretical framework may contribute to data
processing techniques and other applications.