  Directed protein networks with only a few thousand of nodes are rather
complex and do not allow to extract easily the effective influence of one
protein to another taking into account all indirect pathways via the global
network. Furthermore, the different types of activation and inhibition actions
between proteins provide a considerable challenge in the frame work of network
analysis. At the same time these protein interactions are of crucial importance
and at the heart of cellular functioning. We develop the Google matrix analysis
of the protein-protein network from the open public database SIGNOR. The
developed approach takes into account the bi-functional activation or
inhibition nature of interactions between each pair of proteins describing it
in the frame work of Ising-spin matrix transitions. We also apply a recently
developed linear response theory for the Google matrix which highlights a
pathway of proteins whose PageRank probabilities are most sensitive with
respect to two proteins selected for the analysis. This group of proteins is
analyzed by the reduced Google matrix algorithm which allows to determine the
effective interactions between them due to direct and indirect pathways in the
global network. We show that the dominating activation or inhibition function
of each protein can be characterized by its magnetization. The results of this
Google matrix analysis are presented for three examples of selected pairs of
proteins. The developed methods work rapidly and efficiently even for networks
with several million of nodes and can be applied to various biological
networks.


  We cast the metabolism of interacting cells within a statistical mechanics
framework considering both, the actual phenotypic capacities of each cell and
its interaction with its neighbors. Reaction fluxes will be the components of
high-dimensional spin vectors, whose values will be constrained by the
stochiometry and the energy requirements of the metabolism. Within this
picture, finding the phenotypic states of the population turns out to be
equivalent to searching for the equilibrium states of a disordered spin model.
We provide a general solution of this problem for arbitrary metabolic networks
and interactions. We apply this solution to a simplified model of metabolism
and to a complex metabolic network, the central core of the \emph{E. coli}, and
demonstrate that the combination of selective pressure and interactions define
a complex phenotypic space. Cells may specialize in producing or consuming
metabolites complementing each other at the population level and this is
described by an equilibrium phase space with multiple minima, like in a
spin-glass model.


  Gene Ontology (GO) is the most important resource for gene function
annotation. It provides a way to unify biological knowledge across different
species via a dynamic and controlled vocabulary. GO is now widely represented
in the Semantic Web standard Web Ontology Language (OWL). OWL renders a rich
logic constructs to GO but also has its limitations. On the other hand, olog is
a different language for ontology and it is based on category theory. Due to
its solid mathematical background, olog can be rigorously formulated yet
considerably expressive. These advantages make ologs complementary, if not
better than, OWL. We therefore think it worthwhile to adopt ologs for GO and
took an initial step in this work.


  In the standard Susceptible-Infected-Removed (SIR) and
Susceptible-Exposed-Infected-Removed (SEIR) models, the peak of infected
individuals coincides with the in ection point of removed individuals.
Nevertheless, a survey based on the data of the 2009 H1N1 epidemic in Istanbul,
Turkey [19] displayed an unexpected time shift between the hospital referrals
and fatalities. With the motivation of investigating the underlying reason, we
use multistage SIR and SEIR models to provide an explanation for this time
shift. Numerical solutions of these models present strong evidences that the
delay is approximately half of the infection period of the epidemic disease. In
addition, graphs of the classical SIR and the multistage SIR models; and the
classical SEIR and the multistage SEIR models are compared for various epidemic
parameters. Depending on the number of stages, it is observed that the delay
varies for relatively small stage numbers whereas it does not change for large
numbers in multistage systems. One important result that follows immediately
from this observation is that this fixed delay for large numbers explains the
time shift. Additionally, depending on the stage number and the duration of the
epidemic disease, the distance between the points where each infectious stage
reaches its maximum is found approximately both graphically and qualitatively
for both systems. Variations of the time shift, the maximum point of the sum of
all infectious stages, and the in ection point of the removed stage are
observed subject to the stage number N and it is shown that these variations
stay unchanged for greater values of N.


  Cooperative behavior, the costly provision of benefits to others, is common
across all domains of life. This review article discusses cooperative behavior
in the microbial world, mediated by the exchange of extracellular products
called public goods. We focus on model species for which the production of a
public good and the related growth disadvantage for the producing cells are
well described. To unveil the biological and ecological factors promoting the
emergence and stability of cooperative traits we take an interdisciplinary
perspective and review insights gained from both mathematical models and
well-controlled experimental model systems. Ecologically, we include crucial
aspects of the microbial life cycle into our analysis and particularly consider
population structures where an ensemble of local communities (sub populations)
continuously emerge, grow, and disappear again. Biologically, we explicitly
consider the synthesis and regulation of public good production. The discussion
of the theoretical approaches includes general evolutionary concepts,
population dynamics, and evolutionary game theory. As a specific but generic
biological example we consider populations of Pseudomonas putida and its
regulation and utilization of pyoverdines, iron scavenging molecules. The
review closes with an overview on cooperation in spatially extended systems and
also provides a critical assessment of the insights gained from the
experimental and theoretical studies discussed. Current challenges and
important new research opportunities are discussed, including the biochemical
regulation of public goods, more realistic ecological scenarios resembling
native environments, cell to cell signalling, and multi-species communities.


  A major goal of molecular evolutionary biology is to identify loci or regions
of the genome under selection versus those evolving in a neutral manner.
Correct identification allows accurate inference of the evolutionary process
and thus comprehension of historical and contemporary processes driving
phenotypic change and adaptation. A fundamental difficulty lies in
distinguishing sites targeted by selection from both sites linked to these
targets and sites fully independent of selection. These three categories of
sites necessitate attention in light of the debate over the relative importance
of selection versus neutrality and the neutral theory. Modern genomic insights
have proved that complex processes such as linkage, demography, and biased gene
conversion complicate our understanding of the role of neutral versus selective
processes in evolution. In this perspective, we first highlight the importance
of the genomic and (a)biotic context of new mutations to identify the targets
of natural selection. We then present mechanisms that may constrain the
evolution of genomes and bias the inference of selection. We discuss these
mechanisms within the two critical levels that they occur: the population level
and the molecular level. We highlight that they should be taken into account to
correctly distinguish sites across the genome subject to selective or
non-selective forces and stress that a major current field-wide goal is to
quantify the absolute importance of these mechanisms.


  Here, we investigate how competition in the Eden model is affected by short
range dispersal and the requirement that site updates occur only after several
updates of the same site have been attempted previously. The latter models the
effect of tissue or media resistance to invasion. We found that the existence
of tissue intensifies 'Natural Selection' and de-accelerating 'Genetic Drift',
both to a limited extent. More interestingly, our results show that short-range
migration can eliminate Genetic demixing and conceal Natural Selection.


  Healthcare-associated infections (HAIs) remain a public health problem.
Previous work showed intensive care unit (ICU) population structure impacts
methicillin-resistant Staphylococcus aureus (MRSA) rates. Unexplored in that
work was the transient dynamics of this system. We consider the dynamics of
MRSA in an ICU in three different models: 1) a Ross-McDonald model with a
single healthcare staff type, 2) a Ross-McDonald model with nurses and doctors
considered as separate populations and 3) a meta-population model that segments
patients into smaller groups seen by a single nurse. The basic reproduction
number, R0 is derived using the Next Generation Matrix method, while the
importance of the position of patients within the meta-population model is
assessed via stochastic simulation. The single-staff model had an R0 of 0.337,
while the other two models had R0s of 0.278. The meta-population model's R0 was
not sensitive to the time nurses spent with their assigned patients vs.
unassigned patients. This suggests previous results showing that simulated
infection rates are dependent on this parameter are the result of differences
in the transient dynamics between the models, rather than differing long-term
equilibria.


  We consider the question whether a chemical reaction network preserves the
number and stability of its positive steady states upon inclusion of inflow and
outflow reactions. Often a model of a reaction network is presented without
inflows and outflows, while in fact some of the species might be degraded or
leaked to the environment, or be synthesized or transported into the system. We
provide a sufficient and easy-to-check criterion based on the stoichiometry of
the reaction network alone and discuss examples from systems biology.


  The lack of well-structured annotations in a growing amount of RNA expression
data complicates data interoperability and reusability. Commonly - used text
mining methods extract annotations from existing unstructured data descriptions
and often provide inaccurate output that requires manual curation. Automatic
data-based augmentation (generation of annotations on the base of expression
data) can considerably improve the annotation quality and has not been
well-studied. We formulate an automatic augmentation of small RNA-seq
expression data as a classification problem and investigate deep learning (DL)
and random forest (RF) approaches to solve it. We generate tissue and sex
annotations from small RNA-seq expression data for tissues and cell lines of
homo sapiens. We validate our approach on 4243 annotated small RNA-seq samples
from the Small RNA Expression Atlas (SEA) database. The average prediction
accuracy for tissue groups is 98% (DL), for tissues - 96.5% (DL), and for sex -
77% (DL). The "one dataset out" average accuracy for tissue group prediction is
83% (DL) and 59% (RF). On average, DL provides better results as compared to
RF, and considerably improves classification performance for 'unseen' datasets.


  The lack of well-structured metadata annotations complicates there-usability
and interpretation of the growing amount of publicly available RNA expression
data. The machine learning-based prediction of metadata(data augmentation) can
considerably improve the quality of expression data annotation. In this
study,we systematically benchmark deep learning (DL) and random forest
(RF)-based metadata augmentation of tissue, age, and sex using small RNA (sRNA)
expression profiles. We use 4243 annotated sRNA-Seq samples from the small RNA
expression atlas (SEA) database to train and test the augmentation performance.
In general, the DL machine learner outperforms the RF method in almost all
tested cases. The average cross-validated prediction accuracy of the DL
algorithm for tissues is 96.5%, for sex is 77%, and for age is 77.2%. The
average tissue prediction accuracy for a completely new dataset is 83.1% (DL)
and 80.8% (RF). To understand which sRNAs influence DL predictions, we employ
backpropagation-based feature importance scores using the DeepLIFT method,
which enable us to obtain information on biological relevance of sRNAs.


  According to the standard imitation protocol, a less successful player adopts
the strategy of the more successful one faithfully for future success. This is
the cornerstone of evolutionary game theory that explores the vitality of
competing strategies in different social dilemma situations. In our present
work we explore the possible consequences of two slightly modified imitation
protocols that are exaggerated and gradual learning rules. In the former case a
learner does not only adopt, but also enlarges the strategy change for the hope
of a higher income. Similarly, in the latter case a cautious learner does not
adopt the alternative behavior precisely, but takes only a smaller step towards
the other's strategy during the updating process. Evidently, both scenarios
assume that the players' propensity to cooperate may vary gradually between
zero (always defect) and one (always cooperate) where these extreme states
represent the traditional two-strategy social dilemma. We have observed that
while the usage of exaggerated learning has no particular consequence on the
final state, but gradual learning can support cooperation significantly. The
latter protocol mitigates the invasion speeds of both main strategies, but the
decline of successful defector invasion is more significant, hence the biased
impact of the modified microscopic rule on invasion processes explains our
observations.


  Commonly known as head lice, Pediculus humanus capitis are human
ectoparasites which cause infestations in children worldwide. Understanding the
life cycle of head lice is an important step in knowing how to treat lice
infestations, as the parasite behavior depends considerably on its age and
gender. In this work we propose a mathematical model for head lice population
dynamics in hosts who could be or not quarantined and treated. Considering a
lice population structured by age and gender we formulate the model as a system
of hyperbolic PDEs, which can be reduced to compartmental systems of delay or
ordinary differential equations. Besides studying fundamental properties of the
model, such as existence, uniqueness and nonnegativity of solutions, we show
the existence of (in certain cases multiple) equilibria at which the
infestation persists on the host's head. Aiming to assess the performance of
treatments against head lice infestations, by mean of computer experiments and
numerical simulations we investigate four possible treatment strategies. Our
main results can be summarized as follows: (i) early detection is crucial for
quick and efficient eradication of lice infestations; (ii) dimeticone-based
products applied every 4 days effectively remove lice in at most three
applications even in case of severe infestations and (iii) minimization of the
reinfection risk, e.g. by mean of synchronized treatments in
families/classrooms is recommended.


  Recently, the RIDL-SIT technology has been field-tested for control of Aedes
aegypti. The technique consists of releasing genetically modified mosquitoes
carrying a "lethal gene". In 2016 the World Health Organisation (WHO) and the
Pan-American Health Organization (PAHO) recommend to their constituent
countries to test the new technologies proposed to control Aedes aegypti
populations. However, issues concerning effectiveness and ecological impact
have not been thoroughly studied so far. In order to study these issues, we
develop an ecological model compatible with the information available as of
2016. It presents an interdependent dynamics of mosquito populations and food
in an homogeneous setting. Mosquito populations are described in an stochastic
compartmental setup in terms of reaction norms depending on the available food
in the environment. The development of the model allows us to indicate some
critical biological knowledge that is missing and could (should) be produced.
Hybridisation levels, release numbers during and after intervention and
population recovery time after the intervention as a function of intervention
duration and target are calculated under different hypotheses with regard to
the fitness of hybrids and compared with two field studies of actual
interventions. The minimal model should serve as a basis for detailed models
when the necessary information to construct them is produced. For the time
being, the model shows that nature will not clean the non-lethal introgressed
genes.


  We investigate calcium signaling feedback through calcium-activated potassium
channels of a dendritic spine by applying the immersed boundary method with
electrodiffusion. We simulate the stochastic gating of such ion channels and
the resulting spatial distribution of concentration, current, and membrane
voltage within the dendritic spine. In this simulation, the permeability to
ionic flow across the membrane is regulated by the amplitude of chemical
potential barriers. With spatially localized ion channels, chemical potential
barriers are locally and stochastically regulated. This regulation represents
the ion channel gating with multiple subunits, the open and closed states
governed by a continuous-time Markov process. The model simulation
recapitulates an inhibitory action on voltage-sensitive calcium channels by the
calcium-activated potassium channels in a stochastic manner as a
\emph{non-local} feedback loop. The model predicts amplified calcium influx
with more closely placed channel complexes, proposing a potential mechanism of
differential calcium handling by channel distributions. This work provides a
foundation for future computer simulation studies of dendritic spine motility
and structural plasticity.


  In recent years, functional genomics approaches combining genetic information
with bulk RNA-sequencing data have identified the downstream expression effects
of disease-associated genetic risk factors through so-called expression
quantitative trait locus (eQTL) analysis. Single-cell RNA-sequencing creates
enormous opportunities for mapping eQTLs across different cell types and in
dynamic processes, many of which are obscured when using bulk methods. The
enormous increase in throughput and reduction in cost per cell now allow this
technology to be applied to large-scale population genetics studies. Therefore,
we have founded the single-cell eQTLGen consortium (sc-eQTLGen), aimed at
pinpointing disease-causing genetic variants and identifying the cellular
contexts in which they affect gene expression. Ultimately, this information can
enable development of personalized medicine. Here, we outline the goals,
approach, potential utility and early proofs-of-concept of the sc-eQTLGen
consortium. We also provide a set of study design considerations for future
single-cell eQTL studies.


  Background: State-of-the art selection methods fail to identify weak but
cumulative effects of features found in many high-dimensional omics datasets.
Nevertheless, these features play an important role in certain diseases.
  Results: We present Netboost, a three-step dimension reduction technique.
First, a boosting-based filter is combined with the topological overlap measure
to identify the essential edges of the network. Second, sparse hierarchical
clustering is applied on the selected edges to identify modules and finally
module information is aggregated by the first principal components. The primary
analysis is than carried out on these summary measures instead of the original
data. We demonstrate the application of the newly developed Netboost in
combination with CoxBoost for survival prediction of DNA methylation and gene
expression data from 180 acute myeloid leukemia (AML) patients and show, based
on cross-validated prediction error curve estimates, its prediction superiority
over variable selection on the full dataset as well as over an alternative
clustering approach. The identified signature related to chromatin modifying
enzymes was replicated in an independent dataset of AML patients in the phase
II AMLSG 12-09 study. In a second application we combine Netboost with Random
Forest classification and improve the disease classification error in
RNA-sequencing data of Huntington's disease mice.
  Conclusion: Netboost improves definition of predictive variables for survival
analysis and classification. It is a freely available Bioconductor R package
for dimension reduction and hypothesis generation in high-dimensional omics
applications.


  Aspartate carbamoyltransferase (ATCase) is a large dodecameric enzyme with
six active sites that exhibits allostery: its catalytic rate is modulated by
the binding of various substrates at distal points from the active sites. A
recently developed method, bond-to-bond propensity analysis, has proven capable
of predicting allosteric sites in a wide range of proteins using an
energy-weighted atomistic graph obtained from the protein structure and given
knowledge only of the location of the active site. Bond-to-bond propensity
establishes if energy fluctuations at given bonds have significant effects on
any other bond in the protein, by considering their propagation through the
protein graph. In this work, we use bond-to-bond propensity analysis to study
different aspects of ATCase activity using three different protein structures
and sources of fluctuations. First, we predict key residues and bonds involved
in the transition between inactive (T) and active (R) states of ATCase by
analysing allosteric substrate binding as a source of energy perturbations in
the protein graph. Our computational results also indicate that the effect of
multiple allosteric binding is non linear: a switching effect is observed after
a particular number and arrangement of substrates is bound suggesting a form of
long range communication between the distantly arranged allosteric sites.
Second, cooperativity is explored by considering a bisubstrate analogue as the
source of energy fluctuations at the active site, also leading to the
identification of highly significant residues to the T-R transition that
enhance cooperativity across active sites. Finally, the inactive (T) structure
is shown to exhibit a strong, non linear communication between the allosteric
sites and the interface between catalytic subunits, rather than the active
site.


  Human mobility plays a key role on the transformation of local disease
outbreaks into global pandemics. Thus, the inclusion of human movements into
epidemic models has become mandatory for understanding current epidemic
episodes and to design efficient prevention policies. Following this challenge,
here we develop a Markovian framework which enables to address the impact of
recurrent mobility patterns on the epidemic onset at different temporal scales.
This formalism is validated by comparing their predictions with results from
mechanistic simulations. The fair agreement between both theory and simulations
enables to get an analytical expression for the epidemic threshold which
captures the critical conditions triggering epidemic outbreaks. Finally, by
performing an exhaustive analysis of this epidemic threshold, we reveal that
the impact of tuning human mobility on the emergence of diseases is strongly
affected by the temporal scales associated to both epidemiological and mobility
processes.


  Stochastic mechanics---the study of classical stochastic systems governed by
things like master equations and Fokker-Planck equations---exhibits striking
mathematical parallels to quantum mechanics. In this article, we make those
parallels more transparent by presenting a quantum mechanics-like formalism for
deriving a path integral description of systems described by stochastic
differential equations. Our formalism expediently recovers the usual path
integrals (the Martin-Siggia-Rose-Janssen-De Dominicis and Onsager-Machlup
forms) and is flexible enough to account for different variable domains (e.g.
real line versus compact interval), stochastic interpretations, arbitrary
numbers of variables, explicit time-dependence, dimensionful control
parameters, and more. We discuss the implications of our formalism for
stochastic biology.


  Filament formation by non-cytoskeletal enzymes has been known for decades,
yet only relatively recently has its wide-spread role in enzyme regulation and
biology come to be appreciated. This comprehensive review summarizes what is
known for each enzyme confirmed to form filamentous structures in vitro, and
for the many that are known only to form large self-assemblies within cells.
For some enzymes, studies describing both the in vitro filamentous structures
and cellular self-assembly formation are also known and described. Special
attention is paid to the detailed structures of each type of enzyme filament,
as well as the roles the structures play in enzyme regulation and in biology.
Where it is known or hypothesized, the advantages conferred by enzyme
filamentation are reviewed. Finally, the similarities, differences, and
comparison to the SgrAI system are also highlighted.


  Metagenomic studies have increasingly utilized sequencing technologies in
order to analyze DNA fragments found in environmental samples.One important
step in this analysis is the taxonomic classification of the DNA fragments.
Conventional read classification methods require large databases and vast
amounts of memory to run, with recent deep learning methods suffering from very
large model sizes. We therefore aim to develop a more memory-efficient
technique for taxonomic classification. A task of particular interest is
abundance estimation in metagenomic samples. Current attempts rely on
classifying single DNA reads independently from each other and are therefore
agnostic to co-occurence patterns between taxa. In this work, we also attempt
to take these patterns into account. We develop a novel memory-efficient read
classification technique, combining deep learning and locality-sensitive
hashing. We show that this approach outperforms conventional mapping-based and
other deep learning methods for single-read taxonomic classification when
restricting all methods to a fixed memory footprint. Moreover, we formulate the
task of abundance estimation as a Multiple Instance Learning (MIL) problem and
we extend current deep learning architectures with two different types of
permutation-invariant MIL pooling layers: a) deepsets and b) attention-based
pooling. We illustrate that our architectures can exploit the co-occurrence of
species in metagenomic read sets and outperform the single-read architectures
in predicting the distribution over taxa at higher taxonomic ranks.


  Morphogens are secreted signaling molecules that mediate tissue patterning
and growth of embryonic tissues. They are secreted in a localized region and
spread through the tissue to form a graded concentration profile. We present a
cell-based model of morphogen spreading that combines secretion in a local
source, extracellular diffusion and cellular trafficking. We introduce
hydrodynamic modes of morphogen transport and characterize the dynamics of
transport by dispersion relations of these dynamic eigenmodes. These dispersion
relations specify the characteristic relaxation time of a mode as a function of
its wavelength. In a simple model we distinguish two distinct dynamic modes
characterized by different timescales. We find that the slower mode defines the
effective diffusion and degradation as well as the shape of the concentration
profile in steady state. Using our approach we discuss mechanisms of morphogen
transport in the developing wing imaginal disc of the fruit fly
\textit{Drosophila}, distinguishing three transport regimes: transport by
extracellular diffusion, transport by transcytosis and a regime where both
transport mechanisms are combined.


  This paper describes the optimal selection of a control policy to program the
steady state of controlled nonlinear systems with hyperbolic fixed points. This
work is motivated by the field of synthetic biology, in which saddle points are
common (along with limit cycles), and the aim is to program cells to perform
both digital and analog computation, though developing genetic digital
computation has been the main focus. We frame the analog computing challenge of
generating a steady state input-output function inside living cells. To program
the steady state, a data-driven approach is taken wherein an approximation of
the Koopman operator, identified via deep dynamic mode decomposition, is used
to describe the dynamics of the system linearly. The new representation of the
dynamics are then used to solve an optimization problem for the input which
maximizes a direction in state space. Some added structure on the Koopman
operator learning process for controlled systems is given for dynamics that are
separable in the state and input. Finally, the methods are demonstrated on
simulation examples of an incoherent feedforward loop and a combinatorial
promoter system, two common network architectures seen in the field of
synthetic biology.


  Deciding whether and where a system of parametrized ordinary differential
equations displays bistability, that is, has at least two asymptotically stable
steady states for some choice of parameters, is a hard problem. For systems
modeling biochemical reaction networks, we introduce a procedure to determine,
exclusively via symbolic computations, the stability of the steady states for
unspecified parameter values. In particular, our approach fully determines the
stability type of all steady states of a broad class of networks. To this end,
we combine the Hurwitz criterion, reduction of the steady state equations to
one univariate equation, and structural reductions of the reaction network.
Using our method, we prove that bistability occurs in open regions in parameter
space for many relevant motifs in cell signaling.


  Data from discovery proteomic and phosphoproteomic experiments typically
include missing values that correspond to proteins that have not been
identified in the analyzed sample. Replacing the missing values with random
numbers, a process known as "imputation", avoids apparent infinite fold-change
values. However, the procedure comes at a cost: Imputing a large number of
missing values has the potential to significantly impact the results of the
subsequent differential expression analysis. We propose a method that
identifies differentially expressed proteins by ranking their observed changes
with respect to the changes observed for other proteins. Missing values are
taken into account by this method directly, without the need to impute them. We
illustrate the performance of the new method on two distinct datasets and show
that it is robust to missing values and, at the same time, provides results
that are otherwise similar to those obtained with edgeR which is a state-of-art
differential expression analysis method. The new method for the differential
expression analysis of proteomic data is available as an easy to use Python
package.


  Computational prediction of origin of replication (ORI) has been of great
interest in bioinformatics and several methods including GC-skew,
auto-correlation etc. have been explored in the past. In this paper, we have
extended the auto-correlation method to predict ORI location with much higher
resolution for prokaryotes and eukaryotes, which can be very helpful in
experimental validation of the computational predictions. The proposed complex
correlation method (iCorr) converts the genome sequence into a sequence of
complex numbers by mapping the nucleotides to {+1,-1,+i,-i} instead of {+1,-1}
used in the auto-correlation method (here, i is square root of -1). Thus, the
iCorr method exploits the complete spatial information about the positions of
all the four nucleotides unlike the earlier auto-correlation method which uses
the positional information of only one nucleotide. Also, the earlier
auto-correlation method required visual inspection of the obtained graphs to
identify the location of origin of replication. The proposed iCorr method does
away with this need and is able to identify the origin location simply by
picking the peak in the iCorr graph.


  Identifiability is a crucial property for a statistical model since
distributions in the model uniquely determine the parameters that produce them.
In phylogenetics, the identifiability of the tree parameter is of particular
interest since it means that phylogenetic models can be used to infer
evolutionary histories from data. In this paper we introduce a new
computational strategy for proving the identifiability of discrete parameters
in algebraic statistical models that uses algebraic matroids naturally
associated to the models. We then use this algorithm to prove that the tree
parameters are generically identifiable for 2-tree CFN and K3P mixtures. We
also show that the $k$-cycle phylogenetic network parameter is identifiable
under the K2P and K3P models.


  We introduce two algorithms for reconstructing semi-directed level-1
phylogenetic networks from their complete set of 4-leaf subnetworks, known as
quarnets. The first algorithm, the sequential method, begins with a single
quarnet and adds on one leaf at a time until all leaves have been placed. The
second algorithm, the cherry-blob method, functions similarly to cherry-picking
algorithms for phylogenetic trees by identifying exterior network structures
from the quarnets.


  Background: Genotype-phenotype maps provide a meaningful filtration of
sequence space and RNA secondary structures are particular such phenotypes.
Compatible sequences i.e.~sequences that satisfy the base pairing constraints
of a given RNA structure play an important role in the context of neutral
networks and inverse folding. Sequences satisfying the constraints of two
structures simultaneously are called bicompatible and phenotypic change,
induced by erroneously replicating populations of RNA sequences, is closely
connected to bicompatibility. Furthermore, bicompatible sequences are relevant
for riboswitch sequences, beacons of evolution, realizing two distinct
phenotypes.
  Results: We present a full loop energy model Boltzmann sampler of
bicompatible sequences for pairs of structures. The novel dynamic programming
algorithm is based on a topological framework encapsulating the relations
between loops. We utilize our sequence sampler to study the energy spectra and
density of bicompatible sequences, the rankings of the structures and key
properties for evolutionary transitions.
  Conclusion: Our analysis of riboswitch sequences shows that key properties of
bicompatible sequences depend on the particular pair of structures. While there
always exist bicompatible sequences for random structure pairs, they are less
suited to facilitate transitions. We show that native riboswitch sequences
exhibit a distinct signature with regards to the ranking of their two
phenotypes relative to the minimum free energy, suggesting a new criterion for
identifying native sequences and sequences subjected to evolutionary pressure.


  Contagion, broadly construed, refers to anything that can spread infectiously
from peer to peer. Examples include communicable diseases, rumors,
misinformation, ideas, innovations, bank failures, and electrical blackouts.
Sometimes, as in the 1918 Spanish flu epidemic, a contagion mutates at some
point as it spreads through a network. Here, using a simple
susceptible-infected (SI) model of contagion, we explore the downstream impact
of a single mutation event. Assuming that this mutation occurs at a random node
in the contact network, we calculate the distribution of the number of
"descendants," $d$, downstream from the initial "Patient Zero" mutant. We find
that the tail of the distribution decays as $d^{-2}$ for complete graphs,
random graphs, small-world networks, networks with block-like structure, and
other infinite-dimensional networks. This prediction agrees with the observed
statistics of memes propagating and mutating on Facebook, and is expected to
hold for other effectively infinite-dimensional networks, such as the global
human contact network. In a wider context, our approach suggests a possible
starting point for a mesoscopic theory of contagion. Such a theory would focus
on the paths traced by a spreading contagion, thereby furnishing an
intermediate level of description between that of individual nodes and the
total infected population. We anticipate that contagion pathways will hold
valuable lessons, given their role as the conduits through which single
mutations, innovations, or failures can sweep through a network as a whole.


  Many systems in biology and beyond employ collaborative, collective
communication strategies for improved efficiency and adaptive benefit. One such
paradigm of particular interest is the community estimation of a dynamic
signal, when, for example, an epithelial tissue of cells must decide whether to
react to a given dynamic external concentration of stress signaling molecules.
At the level of dynamic cellular communication, however, it remains unknown
what effect, if any, arises from communication beyond the mean field level.
What are the limits and benefits to communication across a network of neighbor
interactions? What is the role of Poissonian vs. super Poissonian dynamics in
such a setting? How does the particular topology of connections impact the
collective estimation and that of the individual participating cells? In this
letter we construct a robust and general framework of signal estimation over
continuous time Markov chains in order to address and answer these questions.
Our results show that in the case of Possonian estimators, the communication
solely enhances convergence speed of the Mean Squared Error (MSE) of the
estimators to their steady-state values while leaving these values unchanged.
However, in the super-Poissonian regime, MSE of estimators significantly
decreases by increasing the number of neighbors. Surprisingly, in this case,
the clustering coefficient of an estimator does not enhance its MSE while
reducing total MSE of the population.


  In the last years it has been shown that Lotka-Volterra mappings constitute a
valuable tool from both the theoretical and the applied points of view, with
developments in very diverse fields such as Physics, Population Dynamics,
Chemistry and Economy. The purpose of this work is to demonstrate that many of
the most important ideas and algebraic methods that constitute the basis of the
quasipolynomial formalism (originally conceived for the analysis of ordinary
differential equations) can be extended into the mapping domain. The extension
of the formalism into the discrete-time context is remarkable as far as the
quasipolynomial methodology had never been shown to be applicable beyond the
differential case. It will be demonstrated that Lotka-Volterra mappings play a
central role in the quasipolynomial formalism for the discrete-time case.
Moreover, the extension of the formalism into the discrete-time domain allows a
significant generalization of Lotka-Volterra mappings as well as a whole
transfer of algebraic methods into the discrete-time context. The result is a
novel and more general conceptual framework for the understanding of
Lotka-Volterra mappings, as well as a new range of possibilities that becomes
open not only for the theoretical analysis of Lotka-Volterra mappings and their
generalizations, but also for the development of new applications.


  Some of the key proteins essential for important cellular processes are
capable of recruiting other proteins from the cytosol to phospholipid
membranes. The physical basis for this cooperativity of binding is,
surprisingly, still unclear. Here, we suggest a general feedback mechanism that
explains cooperativity through mechanochemical coupling mediated by the
mechanical properties of phospholipid membranes. Our theory predicts that
protein recruitment, and therefore also protein pattern formation, involves
membrane deformation, and is strongly affected by membrane composition.


  In general, cellular phenotypes, as measured by concentrations of cellular
components, involve large degrees of freedom. However, recent measurement has
demonstrated that phenotypic changes resulting from adaptation and evolution in
response to environmental changes are effectively restricted to a
low-dimensional subspace. Thus, uncovering the origin and nature of such a
drastic dimension reduction is crucial to understanding the general
characteristics of biological adaptation and evolution. Herein, we first
formulated the dimension reduction in terms of dynamical systems theory:
considering the steady growth state of cells, the reduction is represented by
the separation of a few large singular values of the inverse Jacobian matrix
around a fixed point. We then examined this dimension reduction by numerical
evolution of cells consisting of thousands of chemicals whose concentrations
determine phenotype. As a result of the evolution, phenotypic changes due to
mutations and external perturbations were found to be mainly restricted to a
one-dimensional subspace. One singular value of the inverse Jacobian matrix at
a fixed point of concentrations was significantly larger than the others. The
major phenotypic changes due to mutations and external perturbations occur
along the corresponding left-singular vector, which leads to phenotypic
constraint, and fitness dominantly changes in the same direction. Once such
phenotypic constraint is acquired, phenotypic evolution to a novel environment
takes advantage of this restricted phenotypic direction. This results in the
convergence of phenotypic pathways across genetically different strains, as is
experimentally observed, while accelerating further evolution.


  Replicator dynamics, the continuous-time analogue of Multiplicative Weights
Updates, is the main dynamic in evolutionary game theory. In simple
evolutionary zero-sum games, such as Rock-Paper-Scissors, replicator dynamic is
periodic \cite{zeeman1980population}, however, its behavior in higher
dimensions is not well understood. We provide a complete characterization of
its behavior in zero-sum evolutionary games. We prove that, if and only if, the
system has an interior Nash equilibrium, the dynamics exhibit Poincar\'{e}
recurrence, i.e., almost all orbits come arbitrary close to their initial
conditions infinitely often. If no interior equilibria exist, then all interior
initial conditions converge to the boundary. Specifically, the strategies that
are not in the support of any equilibrium vanish in the limit of all orbits.
All recurrence results furthermore extend to a class of games that generalize
both graphical polymatrix games as well as evolutionary games, establishing a
unifying link between evolutionary and algorithmic game theory. We show that
two degrees of freedom, as in Rock-Paper-Scissors, is sufficient to prove
periodicity.


  Recent collapses of many fisheries across the globe have challenged the
mathematical approach to these systems through classic bioeconomic models.
Decimated populations did not recover as fast as predicted by these models and
depensatory effects were introduced to better fit the dynamics at low
population abundances. Alternative to depensation, modeling captures by
non-linear harvesting functions produces equivalent outcomes at small
abundances, and the dynamics undergoes a bifurcation leading to population
collapse and recovery once catching efforts are above or below certain
thresholds, respectively. The time that a population takes to undergo these
transitions has been mostly overlooked in bioeconomic contexts, though. In this
work we quantify analytically and numerically the times associated to these
collapse and recovery transitions in a model incorporating non-linear
harvesting and immigration in the presence and absence of demographic
stochasticity. Counterintuitively, although species at low abundances are prone
to extinction due to demographic stochasticity, our results show that
stochastic collapse and recovery times are upper bounded by their deterministic
estimates. This occurs over the full range of immigration rates. Our work may
have relevant quantitative implications in the context of fishery management
and rebuilding.


  We derive explicit formulas for the two first moments of he site frequency
spectrum $(SFS_{n,b})_{1\leq b\leq n-1}$ of the Bolthausen-Sznitman coalescent
along with some precise and efficient approximations, even for small sample
sizes $n$. These results provide new $L_2$-asymptotics for some values of
$b=o(n)$. We also study the length of internal branches carrying $b>n/2$
individuals. In this case we obtain the distribution function and a convergence
in law. Our results rely on the random recursive tree construction of the
Bolthausen-Sznitman coalescent.


  While generative models have shown great success in generating
high-dimensional samples conditional on low-dimensional descriptors (learning
e.g. stroke thickness in MNIST, hair color in CelebA, or speaker identity in
Wavenet), their generation out-of-sample poses fundamental problems. The
conditional variational autoencoder (CVAE) as a simple conditional generative
model does not explicitly relate conditions during training and, hence, has no
incentive of learning a compact joint distribution across conditions. We
overcome this limitation by matching their distributions using maximum mean
discrepancy (MMD) in the decoder layer that follows the bottleneck. This
introduces a strong regularization both for reconstructing samples within the
same condition and for transforming samples across conditions, resulting in
much improved generalization. We refer to the architecture as
\emph{transformer} VAE (trVAE). Benchmarking trVAE on high-dimensional image
and tabular data, we demonstrate higher robustness and higher accuracy than
existing approaches. In particular, we show qualitatively improved predictions
for cellular perturbation response to treatment and disease based on
high-dimensional single-cell gene expression data, by tackling previously
problematic minority classes and multiple conditions. For generic tasks, we
improve Pearson correlations of high-dimensional estimated means and variances
with their ground truths from 0.89 to 0.97 and 0.75 to 0.87, respectively.


  One of the most challenging problems in biomedicine and genomics is the
identification of disease biomarkers. In this study, proteomics data from seven
major cancers were used to construct two weighted protein-protein interaction
(PPI) networks i.e., one for the normal and another for the cancer conditions.
We developed rigorous, yet mathematically simple, methodology based on the
degeneracy at -1 eigenvalues to identify structural symmetry or motif
structures in network. Utilising eigenvectors corresponding to degenerate
eigenvalues in the weighted adjacency matrix, we identified structural symmetry
in underlying weighted PPI networks constructed using seven cancer data.
Functional assessment of proteins forming these structural symmetry exhibited
the property of cancer hallmarks. Survival analysis refined further this
protein list proposing BMI, MAPK11, DDIT4, CDKN2A, and FYN as putative
multi-cancer biomarkers. The combined framework of networks and spectral graph
theory developed here can be applied to identify symmetrical patterns in other
disease networks to predict proteins as potential disease biomarkers.


  One of the basic questions of phylogenomics is how gene function evolves,
whether among species or inside gene families. In this chapter, we provide a
brief overview of the problems associated with defining gene function in a
manner which allows comparisons which are both large scale and evolutionarily
relevant. The main source of functional data, despite its limitations, is
transcriptomics. Functional data provides information on evolutionary
mechanisms primarily by showing which functional classes of genes evolve under
stronger or weaker purifying or adaptive selection, and on which classes of
mutations (e.g., substitutions or duplications). However, the example of the
"ortholog conjecture" shows that we are still not at a point where we can
confidently study phylogenomically the evolution of gene function at a precise
scale.


  Synergistic and antagonistic interactions in multi-species populations - such
as resource sharing and competition - result in remarkably diverse behaviors in
populations of interacting cells, such as in soil or human microbiomes, or
clonal competition in cancer. The degree of inter- and intra-specific
interaction can often be quantified through the notion of an ecological
"niche". Typically, weakly interacting species that occupy largely distinct
niches result in stable mixed populations, while strong interactions and
competition for the same niche results in rapid extinctions of some species and
fixations of others. We investigate the transition of a deterministically
stable mixed population to a stochasticity-induced fixation as a function of
the niche overlap between the two species. We also investigate the effect of
the niche overlap on the population stability with respect to external
invasions. Our results have important implications for a number of experimental
systems.


  The commonly observed negative correlation between the number of species in
an ecological community and disease risk, typically referred to as "the
dilution effect", has received a substantial amount of attention over the past
decade. Attempts to test this relationship experimentally have revealed that,
in addition to the mean disease risk decreasing with species number, so too
does the variance of disease risk. This is referred to as the "variance
reduction effect", and has received relatively little attention in the
disease-diversity literature. Here, we set out to clarify and quantify some of
these relationships in an idealized model of a randomly assembled multi-species
community undergoing an epidemic. We specifically investigate the variance of
the community disease reproductive ratio, a multi-species extension of the
basic reproductive ratio R_0, for a family of random-parameter meta-community
SIR models, and show how the variance of community $R_0$ varies depending on
whether transmission is density or frequency-dependent. We finally outline
areas of further research on how changes in variance affect transmission
dynamics in other systems.


  The tensor-structured parametric analysis (TPA) has been recently developed
for simulating and analysing stochastic behaviours of gene regulatory networks
[Liao et. al., 2015]. The method employs the Fokker-Planck approximation of the
chemical master equation, and uses the Quantized Tensor Train (QTT) format, as
a low-parametric tensor-structured representation of classical matrices and
vectors, to approximate the high-dimensional stationary probability
distribution. This paper presents a detailed error analysis of all
approximation steps of the TPA regarding validity and accuracy, including
modelling error, artificial boundary error, discretization error, tensor
rounding error, and algebraic error. The error analysis is illustrated using
computational examples, including the death-birth process and a 50-dimensional
isomerization reaction chain.


  As proteins with similar structures often have similar functions, analysis of
protein structures can help predict protein functions and is thus important. We
consider the problem of protein structure classification, which computationally
classifies the structures of proteins into pre-defined groups. We develop a
weighted network that depicts the protein structures, and more importantly, we
propose the first graphlet-based measure that applies to weighted networks.
Further, we develop a deep neural network (DNN) composed of both convolutional
and recurrent layers to use this measure for classification. Put together, our
approach shows dramatic improvements in performance over existing
graphlet-based approaches on 36 real datasets. Even comparing with the
state-of-the-art approach, it almost halves the classification error. In
addition to protein structure networks, our weighted-graphlet measure and DNN
classifier can potentially be applied to classification of other weighted
networks in computational biology as well as in other domains.


  In a social dilemma, cooperation is collectively optimal, yet individually
each group member prefers to defect. A class of successful strategies of direct
reciprocity were recently found for the iterated prisoner's dilemma and for the
iterated three-person public-goods game: By a successful strategy, we mean that
it constitutes a cooperative Nash equilibrium under implementation error, with
assuring that the long-term payoff never becomes less than the co-players'
regardless of their strategies, when the error rate is small. Although we have
a list of actions prescribed by each successful strategy, the rationale behind
them has not been fully understood for the iterated public-goods game because
the list has hundreds of entries to deal with every relevant history of
previous interactions. In this paper, we propose a method to convert such
history-based representation into an automaton with a minimal number of states.
Our main finding is that a successful strategy for the iterated three-person
public-goods game can be represented as a $10$-state automaton by this method.
In this automaton, each state can be interpreted as the player's internal
judgement of the situation, such as trustworthiness of the co-players and the
need to redeem oneself after defection. This result thus suggests a
comprehensible way to choose an appropriate action at each step towards
cooperation based on a situational judgement, which is mapped from the history
of interactions.


  We present a simple, modular graph-based convolutional neural network that
takes structural information from protein-ligand complexes as input to generate
models for activity and binding mode prediction. Complex structures are
generated by a standard docking procedure and fed into a dual-graph
architecture that includes separate sub-networks for the ligand bonded topology
and the ligand-protein contact map. This network division allows contributions
from ligand identity to be distinguished from effects of protein-ligand
interactions on classification. We show, in agreement with recent literature,
that dataset bias drives many of the promising results on virtual screening
that have previously been reported. However, we also show that our neural
network is capable of learning from protein structural information when, as in
the case of binding mode prediction, an unbiased dataset is constructed. We
develop a deep learning model for binding mode prediction that uses docking
ranking as input in combination with docking structures. This strategy mirrors
past consensus models and outperforms the baseline docking program in a variety
of tests, including on cross-docking datasets that mimic real-world docking use
cases. Furthermore, the magnitudes of network predictions serve as reliable
measures of model confidence


  In this work, we design a type of controller that consists of adding a
specific set of reactions to an existing mass-action chemical reaction network
in order to control a target species. This set of reactions is effective for
both deterministic and stochastic networks, in the latter case controlling the
mean as well as the variance of the target species. We employ a type of network
property called absolute concentration robustness (ACR). We provide
applications to the control of a multisite phosphorylation model as well as a
receptor-ligand signaling system. For this framework, we use the so-called
deficiency zero theorem from chemical reaction network theory as well as
multiscaling model reduction methods. We show that the target species has
approximately Poisson distribution with the desired mean. We further show that
ACR controllers can bring robust perfect adaptation to a target species and are
complementary to a recently introduced antithetic feedback controller used for
stochastic chemical reactions.


  The evolution of various competing cell types in tissues, and the resulting
persistent tissue population, is studied numerically and analytically in a
particle-based model of active tissues. Mutations change the properties of
cells in various ways, including their mechanical properties. Each mutation
results in an advantage or disadvantage to grow in the competition between
different cell types. While changes in signaling processes and biochemistry
play an important role, we focus on changes in the mechanical properties by
studying the result of variation of growth force and adhesive
cross-interactions between cell types. For independent mutations of growth
force and adhesion strength, the tissue evolves towards cell types with high
growth force and low internal adhesion strength, as both increase the
homeostatic pressure. Motivated by biological evidence, we postulate a coupling
between both parameters, such that an increased growth force comes at the cost
of a higher internal adhesion strength or vice versa. This tradeoff controls
the evolution of the tissue, ranging from unidirectional evolution to very
heterogeneous and dynamic populations. The special case of two competing cell
types reveals three distinct parameter regimes: Two in which one cell type
outcompetes the other, and one in which both cell types coexist in a highly
mixed state. Interestingly, a single mutated cell alone suffices to reach the
mixed state, while a finite mutation rate affects the results only weakly.
Finally, the coupling between changes in growth force and adhesion strength
reveals a mechanical explanation for the evolution towards intra-tumor
heterogeneity, in which multiple species coexist even under a constant
evolutianary pressure.


  Motivation: Discovering functional modules in protein-protein interaction
(PPI) networks by optimization methods remains a longstanding challenge in
biology. Traditional algorithms simply consider strong protein complexes that
can be found in the original network by optimizing some metrics, which causes
obstacles for the discovery of weak and hidden complexes shielded by stronger
complexes. Also, protein complexes are not only in different density but also
in a large range of scale, making it extremely difficult to be detected. Toward
this objective, we propose a hierarchical hidden community approach to predict
protein complexes. Results: We propose a method called HirHide (Hierarchical
Hidden Community Detection), which can be combined with traditional community
detection methods to enable them to discover hierarchical hidden communities.
It is the first community detection algorithm that can find a hierarchical
structure as well as hidden structure. We compare the performance of three
traditional methods with their HirHide versions. Experimental results show that
the HirHide methods using traditional methods as the base algorithms achieve
better performance, sometimes even significantly outperform the baselines.


  Urbanization drives the epidemiology of infectious diseases to many threats
and new challenges. In this research, we study the interplay between human
mobility and dengue outbreaks in the complex urban environment of the
city-state of Singapore. We integrate both stylized and mobile phone
data-driven mobility patterns in an agent-based transmission model in which
humans and mosquitoes are represented as agents that go through the epidemic
states of dengue. We monitor with numerical simulations the system-level
response to the epidemic by comparing our results with the observed cases
reported during the 2013 and 2014 outbreaks. Our results show that human
mobility is a major factor in the spread of vector-borne diseases such as
dengue even on the short scale corresponding to intra-city distances. We
finally discuss the advantages and the limits of mobile phone data and
potential alternatives for assessing valuable mobility patterns for modeling
vector-borne diseases outbreaks in cities.


  We present BoolSi, an open-source cross-platform command line tool for
distributed simulations of deterministic Boolean networks with synchronous
update. It uses MPI standard to support execution on computational clusters, as
well as parallel processing on a single computer. BoolSi can be used to model
the behavior of complex dynamic networks, such as gene regulatory networks. In
particular, it allows for identification and statistical analysis of network
attractors. We perform a case study of the activity of a cambium cell to
demonstrate the capabilities of the tool.


  RNA interference, particularly siRNA induced gene silencing is becoming an
important avenue of modern therapeutics. The siRNA is delivered to the cells as
short double helical RNA which becomes single stranded for forming the RISC
complex. Significant experimental evidence is available for most of the steps
except the process of the separation of the two strands. We have attempted to
understand the pathway for double stranded siRNA (dsRNA) to single stranded
(ssRNA) molecules using steered molecular dynamics simulations. As the process
is completely unexplored we have applied force from all possible directions
restraining all possible residues to convert dsRNA to ssRNA. We found pulling
one strand along helical axis direction restraining far end of the other strand
demands excessive force for ssRNA formation. Pulling a central residue of one
strand, in a direction perpendicular to the helix axis, while keeping the base
paired residue fixed requires intermediate force for strand separation.
Moreover we found that in this process the force requirement is quite high for
the first bubble formation (nucleation energy) and the bubble propagation
energies are quite small. We hypothesize this could be the mode of action
adopted by the proteins in the cells.


  Indirect reciprocity is a foundational mechanism of human cooperation.
Existing models of indirect reciprocity fail to robustly support social
cooperation: image scoring models fail to provide robust incentives, while
social standing models are not informationally robust. Here we provide a new
model of indirect reciprocity based on simple, decentralized records: each
individual's record depends on their own past behavior alone, and not on their
partners' past behavior or their partners' partners' past behavior. When social
dilemmas exhibit a coordination motive (or strategic complementarity), tolerant
trigger strategies based on simple records can robustly support positive social
cooperation and exhibit strong stability properties. In the opposite case of
strategic substitutability, positive social cooperation cannot be robustly
supported. Thus, the strength of short-run coordination motives in social
dilemmas determines the prospects for robust long-run cooperation.


  During beer production, yeast generate ethanol that is exported to the
extracellular environment where it accumulates. Depending on the initial
carbohydrate concentration in the wort a high concentration of ethanol can be
achieved in beer, often higher than 20% (v/v). It is noteworthy that the
effects of elevated ethanol concentrations generated during beer fermentation
resemble those of heat shock stress, with similar responses observed in both
situations, such as the activation of proteostasis and protein quality control
mechanisms in different cell compartments, including endoplasmic reticulum
(ER), mitochondria, and cytosol. Despite the extensive published molecular and
biochemical data regarding the roles of proteostasis in different organelles of
yeast cells, little is known about how this mechanism impacts beer fermentation
and how different proteostasis mechanisms found in ER, mitochondria, and
cytosol communicate with each other during ethanol/fermentative stress.
Supporting this integrative view, transcriptome data analysis was applied using
publicly available information for a lager yeast strain grown under in beer
production conditions. The transcriptome data indicated upregulation of genes
that encode chaperones, co-chaperones, unfolded protein response elements in ER
and mitochondria, ubiquitin ligases, proteasome components, N-glycosylation
quality control pathway proteins, and components of processing bodies
(p-bodies) and stress granules (SGs) during lager beer fermentation. Thus, the
main purpose of this hypothesis and theory manuscript is to provide a concise
picture of how inter-organellar proteostasis mechanisms are connected with one
another and with biological processes that may modulate the viability and/or
vitality of yeast populations during HG/VHG beer fermentation and serial
repitching.


  Cancer is one of the leading causes of human death. Many efforts have made to
understand its mechanism and have further identified many proteins and DNA
sequence variations as suspected targets for therapy. However, drugs targeting
these targets have low success rates, suggesting the basic mechanism still
remains unclear. Here, we develop a computational software combining Cox
proportional-hazards model and stability-selection to unearth an overlooked,
yet the most important cancer drivers hidden in massive data from The Cancer
Genome Atlas (TCGA), including 11,574 RNAseq samples and clinic data.
Generally, noncoding RNAs primarily regulate cancer deaths and work as the
deadliest cancer inducers and repressors, in contrast to proteins as
conventionally thought. Especially, processed-pseudogenes serve as the primary
cancer inducers, while lincRNA and antisense RNAs dominate the repressors.
Strikingly, noncoding RNAs serves as the universal strongest regulators for all
cancer types although personal clinic variables such as alcohol and smoking
significantly alter cancer genome. Furthermore, noncoding RNAs also work as
central hubs in cancer regulatory network and as biomarkers to discriminate
cancer types. Therefore, noncoding RNAs overall serve as the deadliest cancer
regulators, which refreshes the basic concept of cancer mechanism and builds a
novel basis for cancer research and therapy. Biological functions of
pseudogenes have rarely been recognized. Here we reveal them as the most
important cancer drivers for all cancer types from big data, breaking a wall to
explore their biological potentials.


  Calculating the similarities between a pair of genomic sequences is one of
the most fundamental computational steps in genomic analysis. This step --
called sequence alignment -- is the computational bottleneck because: (1) it is
implemented using quadratic-time dynamic programming algorithms and (2) the
majority of candidate locations in the reference genome do not align with a
given read due to high dissimilarity. Calculating the alignment of such
incorrect candidate locations consumes an overwhelming majority of a modern
read mapper's execution time.
  In this thesis, we introduce four new algorithms (GateKeeper, Shouji, MAGNET,
and SneakySnake) that function as a pre-alignment step and aim to filter out
most incorrect candidate locations. The first key idea of our pre-alignment
filters is to provide high filtering accuracy by correctly detecting all
similar segments shared between two sequences. The second key idea is to
exploit the massively parallel architecture of modern FPGAs for accelerating
our filtering algorithms. We also develop an efficient CPU implementation of
the SneakySnake algorithm for commodity desktops and servers. We evaluate the
benefits and downsides of our pre-alignment filtering approach in detail using
12 real datasets. In our evaluation, we demonstrate that our hardware
pre-alignment filters show two to three orders of magnitude speedup over their
equivalent CPU implementations. We also demonstrate that integrating our
hardware pre-alignment filters with the state-of-the-art read aligners reduces
the aligner's execution time by up to 21.5x. Finally, we show that efficient
CPU implementation of pre-alignment filtering still provides significant
benefits. We show that SneakySnake on average reduces the execution time of the
best performing CPU-based read aligners Edlib and Parasail, by up to 43x and
57.9x, respectively.


  Chronic viral infections can persist in an infected person for decades. From
the perspective of the virus, a single infection can span thousands of
generations, leading to a highly diverse population of viruses with its own
complex evolutionary history. We propose a mathematical framework for
understanding how the emergence of new viral strains and phenotype within
infected persons affects the population-level control of those infections by
both non-curative treatment and chemo-prophylactic measures. We consider the
within-host emergence of new strains that lack phenotype novelty and also the
evolution of variability in contagiousness, resistance to therapy, and efficacy
of prophylaxis. Our framework balances the need for verisimilitude with our
desire to retain a model that can be approached analytically. We show how to
compute the population-level basic reproduction number accounting for the
within-host evolutionary process where new phenotypes emerge and are lost in
infected persons, which we also extend to include both treatment and
prophylactic control efforts. This allows us to make clear statements about
both the global and relative efficacy of different control efforts accounting
for within-host phenotypic evolution. Finally, we give expressions for the
endemic equilibrium of these models for certain constrained versions of the
within-host evolutionary model providing a potential method for estimating
within-host evolutionary parameters from population-level genetic sequence
data.


  We consider the problem of testing whether the points in a complex or real
variety with non-zero coordinates form a multiplicative group or, more
generally, a coset of a multiplicative group. For the coset case, we study the
notion of shifted toric varieties which generalizes the notion of toric
varieties. This requires a geometric view on the varieties rather than an
algebraic view on the ideals. We present algorithms and computations on 129
models from the BioModels repository testing for group and coset structures
over both the complex numbers and the real numbers. Our methods over the
complex numbers are based on Gr\"obner basis techniques and binomiality tests.
Over the real numbers we use first-order characterizations and employ real
quantifier elimination. In combination with suitable prime decompositions and
restrictions to subspaces it turns out that almost all models show coset
structure. Beyond our practical computations, we give upper bounds on the
asymptotic worst-case complexity of the corresponding problems by proposing
single exponential algorithms that test complex or real varieties for toricity
or shifted toricity. In the positive case, these algorithms produce generating
binomials. In addition, we propose an asymptotically fast algorithm for testing
membership in a binomial variety over the algebraic closure of the rational
numbers.


  We propose a generative model and an inference scheme for epidemic processes
on dynamic, adaptive contact networks. Network evolution is formulated as a
link-Markovian process, which is then coupled to an individual-level stochastic
SIR model, in order to describe the interplay between epidemic dynamics on a
network and network link changes. A Markov chain Monte Carlo framework is
developed for likelihood-based inference from partial epidemic observations,
with a novel data augmentation algorithm specifically designed to deal with
missing individual recovery times under the dynamic network setting. Through a
series of simulation experiments, we demonstrate the validity and flexibility
of the model as well as the efficacy and efficiency of the data augmentation
inference scheme. The model is also applied to a recent real-world dataset on
influenza-like-illness transmission with high-resolution social contact
tracking records.


  DNA sequencing, especially of microbial genomes and metagenomes, has been at
the core of recent research advances in large-scale comparative genomics. The
data deluge has resulted in exponential growth in genomic datasets over the
past years and has shown no sign of slowing down. Several recent attempts have
been made to tame the computational burden of sequence search on these terabyte
and petabyte-scale datasets, including raw reads and assembled genomes.
However, no known implementation provides both fast query and construction
time, keeps the low false-positive requirement, and offers cheap storage of the
data structure. We propose a data structure for search called RAMBO (Repeated
And Merged BloOm Filter) which is significantly faster in query time than
state-of-the-art genome indexing methods- COBS (Compact bit-sliced signature
index), Sequence Bloom Trees, HowDeSBT, and SSBT. Furthermore, it supports
insertion and query process parallelism, cheap updates for streaming inputs,
has a zero false-negative rate, a low false-positive rate, and a small index
size. RAMBO converts the search problem into set membership testing among $K$
documents. Interestingly, it is a count-min sketch type arrangement of a
membership testing utility (Bloom Filter in our case). The simplicity of the
algorithm and embarrassingly parallel architecture allows us to stream and
index a 170TB whole-genome sequence dataset in a mere 9 hours on a cluster of
100 nodes while competing methods require weeks.


  Protein-protein interactions (PPIs) are of fundamental importance for the
human body, and the knowledge of their existence can facilitate very important
tasks like drug target developing and therapy design. The high-throughput
experiments for detecting new PPIs are costly and time-consuming, stressing the
need for new computational systems able to generate high-quality PPIs
predictions. These systems have to face two main problems: the high
incompleteness of the human interactome and the lack of high-quality negative
protein-protein interactions (i.e. proteins that are known to not interact).
The latter is usually overlooked by the PPIs prediction systems, causing a
significant bias in the performances and metrics. In this work, we compare
methods for simulating negative knowledge using highly reliable training and
test sets. Moreover, we measure the performances of two state-of-the-art
systems when very reliable settings are adopted.


  One of the cornerstones in combating the HIV pandemic is being able to assess
the current state and evolution of local HIV epidemics. This remains a complex
problem, as many HIV infected individuals remain unaware of their infection
status, leading to parts of HIV epidemics being undiagnosed and under-reported.
To that end, we firstly present a method to learn epidemiological parameters
from phylogenetic trees, using approximate Bayesian computation (ABC). The
epidemiological parameters learned as a result of applying ABC are subsequently
used in epidemiological models that aim to simulate a specific epidemic.
Secondly, we continue by describing the development of a tree statistic, rooted
in coalescent theory, which we use to relate epidemiological parameters to a
phylogenetic tree, by using the simulated epidemics. We show that the presented
tree statistic enables differentiation of epidemiological parameters, while
only relying on phylogenetic trees, thus enabling the construction of new
methods to ascertain the epidemiological state of an HIV epidemic. By using
genetic data to infer epidemic sizes, we expect to enhance understanding of the
portions of the infected population in which diagnosis rates are low.


  In this paper we aim to study an optimal insurance premium level for
health-care in a deterministic and stochastic SIR models with migration fluxes
and vaccination of population. The studied model considers two standard SIR
centres connected via links and continuous migration fluxes. The premium is
calculated using the basic equivalence principle. Even in this simple setup
there are non-intuitive results that illustrate how the premium depends on
migration rates, severeness of a disease and initial distribution of healthy
and infected individuals through the centres. We investigate how the
vaccination program effects the insurance costs by comparing the savings in
benefits with the expenses for vaccination. We compare the results of
deterministic and stochastic models.


  One of the first beings affected by changes in the climate are trees, one of
our most vital resources. In this study tree species interaction and the
response to climate in different ecological environments is observed by
applying a joint species distribution model to different ecological domains in
the United States. Joint species distribution models are useful to learn
inter-species relationships and species response to the environment. The
climates' impact on the tree species is measured through species abundance in
an area. We compare the model's performance across all ecological domains and
study the sensitivity of the climate variables. With the prediction of
abundances, tree species populations can be predicted in the future and measure
the impact of climate change on tree populations.


  Vector-borne epidemics are the result of the combination of different factors
such as the crossed contagions between humans and vectors, their demographic
distribution and human mobility among others. The current availability of
information about the former ingredients demands their incorporation to current
mathematical models for vector-borne disease transmission. Here, relying on
metapopulation dynamics, we propose a framework whose results are in fair
agreement with those obtained from mechanistic simulations. This framework
allows us to derive an expression of the epidemic threshold capturing with high
accuracy the conditions leading to the onset of epidemics. Driven by these
insights, we obtain a prevalence indicator to rank the patches according to the
risk of being affected by a vector-borne disease. We illustrate the utility of
this epidemic risk indicator by reproducing the spatial distribution Dengue
cases reported in the city of Santiago de Cali (Colombia) from 2015 to 2016.


  The dynamics of a population expanding into unoccupied habitat has been
primarily studied for situations in which growth and dispersal parameters are
uniform in space or vary in one dimension. Here we study the influence of
finite-sized individual inhomogeneities and their collective effect on front
speed if randomly placed in a two-dimensional habitat. We use an
individual-based model to investigate the front dynamics for a region in which
dispersal or growth of individuals is reduced to zero (obstacles) or increased
above the background (hotspots), respectively. In a regime where front dynamics
is determined by a local front speed only, a principle of least time can be
employed to predict front speed and shape. The resulting analytical solutions
motivate an event-based algorithm illustrating the effects of several obstacles
or hotspots. We finally apply the principle of least time to large
heterogeneous environments by solving the Eikonal equation numerically.
Obstacles lead to a slow-down that is dominated by the number density and width
of obstacles, but not by their precise shape. Hotspots result in a speedup,
which we characterise as function of hotspot strength and density. Our findings
emphasise the importance of taking the dimensionality of the environment into
account.


  In 2000, Gillespie rehabilitated the chemical Langevin equation (CLE) by
describing two conditions that must be satisfied for it yield a valid
approximation of the chemical master equation (CME). In this work, we construct
an original path integral description of the CME, and show how applying
Gillespie's two conditions to it directly leads to a path integral equivalent
to the CLE. We compare this approach to the path integral equivalent of a large
system size derivation, and show that they are qualitatively different. In
particular, both approaches involve converting many sums into many integrals,
and the difference between the two methods is essentially the difference
between using the Euler-Maclaurin formula and using Riemann sums. Our results
shed light on how path integrals can be used to conceptualize coarse-graining
biochemical systems, and are readily generalizable.


  Cellular behaviors such as migration, division, and differentiation rely on
precise timing, and yet the molecular events that govern these behaviors are
highly stochastic. We investigate regulatory strategies that decrease the
timing noise of molecular events. Autoregulatory feedback increases noise. Yet,
we find that in the presence of regulation by a second species, autoregulatory
feedback decreases noise. To explain this finding, we develop a method to
calculate the optimal regulation function that minimizes the timing noise. The
method reveals that the combination of feedback and regulation minimizes noise
by maximizing the number of molecular events that must happen in sequence
before a threshold is crossed. We compute the optimal timing precision for all
two-node networks with regulation and feedback, derive a generic lower bound on
timing noise, and discuss our results in the context of neuroblast migration
during Caenorhabditis elegans development.


  Viral sequence classification is an important task in pathogen detection,
epidemiological surveys and evolutionary studies. Statistical learning methods
are widely used to classify and identify viral sequences in samples from
environments. These methods face several challenges associated with the nature
and properties of viral genomes such as recombination, mutation rate and
diversity. Also, new generations of sequencing technologies rise other
difficulties by generating massive amounts of fragmented sequences. While
linear classifiers are often used to classify viruses, there is a lack of
exploration of the accuracy space of existing models in the context of
alignment free approaches. In this study, we present an exhaustive assessment
procedure exploring the power of linear classifiers in genotyping and subtyping
partial and complete genomes. It is applied to the Hepatitis C viruses (HCV).
Several variables are considered in this investigation such as classifier types
(generative and discriminative) and their hyper-parameters (smoothing value and
regularization penalty function), the classification task (genotyping and
subtyping), the length of the tested sequences (partial and complete) and the
length of k-mer words. Overall, several classifiers perform well given a set of
precise combination of the experimental variables mentioned above. Finally, we
provide the procedure and benchmark data to allow for more robust assessment of
classification from virus genomes.


  The dynamics of infectious diseases propagating in populations depends both
on human interaction patterns, the contagion process and the pathogenesis
within hosts. The immune system follows a circadian rhythm and, consequently,
the chance of getting infected varies with the time of day an individual is
exposed to the pathogen. The movement and interaction of people also follow
24-hour cycles, which couples these two phenomena. We use a stochastic
metapopulation model informed by hourly mobility data for two medium-sized
Chinese cities. By this setup, we investigate how the epidemic risk depends on
the difference of the clocks governing the population movement and the immune
systems. In most of the scenarios we test, we observe circadian rhythms would
constrain the pace and extent of disease emergence. The three measures
(strength, outward transmission risk and introduction risk) are highly
correlated with each other. For example of the Yushu City, outward transmission
risk and introduction risk are correlated with a Pearson's correlation
coefficient of 0.83, and the risks correlate to strength with coefficients of
$-0.85$ and $-0.75$, respectively (all have $p<0.05$), in simulations with no
circadian effect and $R_0=1.5$. The relation between the circadian rhythms of
the immune system and daily routines in human mobility can affect the pace and
extent of infectious disease spreading. Shifting commuting times could mitigate
the emergence of outbreaks.


  Biochemical systems that express certain chemical species of interest at the
same level at any positive equilibrium are called "absolute concentration
robust" (ACR). These species behave in a stable, predictable way, in the sense
that their expression is robust with respect to sudden changes in the species
concentration, regardless the new positive equilibrium reached by the system.
Such a property has been proven to be fundamentally important in certain gene
regulatory networks and signaling systems. In the present paper, we
mathematically prove that a well-known class of ACR systems studied by Shinar
and Feinberg in 2010 hides an internal integral structure. This structure
confers these systems with a higher degree of robustness that what was
previously unknown. In particular, disturbances much more general than sudden
changes in the species concentrations can be rejected, and robust perfect
adaptation is achieved. Significantly, we show that these properties are
maintained when the system is interconnected with other chemical reaction
networks. This key feature enables design of insulator devices that are able to
buffer the loading effect from downstream systems - a crucial requirement for
modular circuit design in synthetic biology.


  Bacterial populations in natural conditions are expected to experience
stochastic environmental fluctuations, and in addition, environments are
affected by bacterial activities since they consume substrates and excrete
various chemicals. We here study possible outcomes of population dynamics and
evolution under the repeated cycle of substrate-rich conditions and starvation,
called the "feast-famine cycle", by a simple stochastic model with the
trade-off relationship between the growth rate and the growth yield or the
death rate. In the model, the feast (substrate-rich) period is led by a
stochastic substrate addition event, while the famine (starvation) period is
evoked because bacteria use the supplied substrate. Under the repeated
feast-famine cycle, the bacterial population tends to increase the growth rate,
even though that tends to decrease the total population size due to the
trade-off. Analysis of the model shows that the ratio between the growth rate
and the death rate becomes the effective fitness of the population. Hence, the
functional form of the trade-off between the growth and death rate determines
if the bacterial population eventually goes extinct as an evolutionary
consequence. We then show that the increase of the added substrate in the feast
period can drive the extinction faster. Overall, the model sheds light on
non-trivial possible outcomes under repeated feast-famine cycles.


  Phylogenetic networks are generalizations of phylogenetic trees that allow
the representation of reticulation events such as horizontal gene transfer or
hybridization, and can also represent uncertainty in inference. A subclass of
these, tree-based phylogenetic networks, have been introduced to capture the
extent to which reticulate evolution nevertheless broadly follows tree-like
patterns. Several important operations that change a general phylogenetic
network have been developed in recent years, and are important for allowing
algorithms to move around spaces of networks; a vital ingredient in finding an
optimal network given some biological data. A key such operation is the Nearest
Neighbor Interchange, or NNI. While it is already known that the space of
unrooted phylogenetic networks is connected under NNI, it has been unclear
whether this also holds for the subspace of tree-based networks. In this paper
we show that the space of unrooted tree-based phylogenetic networks is indeed
connected under the NNI operation. We do so by explicitly showing how to get
from one such network to another one without losing tree-basedness along the
way. Moreover, we introduce some new concepts, for instance ``shoat networks'',
and derive some interesting aspects concerning tree-basedness. Last, we use our
results to derive an upper bound on the size of the space of tree-based
networks.


  Extracting an understanding of the underlying system from high dimensional
data is a growing problem in science. Discovering informative and meaningful
features is crucial for clustering, classification, and low dimensional data
embedding. Here we propose to construct features based on their ability to
discriminate between clusters of the data points. We define a class of problems
in which linear separability of clusters is hidden in a low dimensional space.
We propose an unsupervised method to identify the subset of features that
define a low dimensional subspace in which clustering can be conducted. This is
achieved by averaging over discriminators trained on an ensemble of proposed
cluster configurations. We then apply our method to single cell RNA-seq data
from mouse gastrulation, and identify 27 key transcription factors (out of 409
total), 18 of which are known to define cell states through their expression
levels. In this inferred subspace, we find clear signatures of known cell types
that eluded classification prior to discovery of the correct low dimensional
subspace.


  In this work, a periodically-forced seasonal non-autonomous system of a
non-linear ordinary differential equation is developed that captures the
dynamics of Lassa fever transmission and seasonal variation in the birth of
mastomys rodents where time was measured in days to capture seasonality. It was
shown that the model is epidemiologically meaningful and mathematically
well-posed by using the results from the qualitative properties of the solution
of the model. It was established that in order to eliminate Lassa fever
disease, treatments with Ribavirin must be provided early to reduce mortality
and other preventive measures like an educational campaign, community hygiene,
Isolation of infected humans, and culling/destruction of rodents must be
applied to also reduce the morbidity of the disease. Finally, the obtained
results gave a primer framework for planning and designing cost-effective
strategies for good interventions in eliminating Lassa fever.


  Here we consider a game theoretic model of multilevel selection in which
individuals compete based on their payoff and groups also compete based on the
average payoff of group members. Our focus is on multilevel social dilemmas:
games in which individuals are best off cheating, while groups of individuals
do best when composed of many cooperators. We analyze the dynamics of the
two-level replicator dynamics, a nonlocal hyperbolic PDE describing
deterministic birth-death dynamics for both individuals and groups. While past
work on such multilevel dynamics has restricted attention to scenarios with
exactly solvable within-group dynamics, we use comparison principles and an
invariant property of the tail of the population distribution to extend our
analysis to all possible two-player, two-strategy social dilemmas. We
characterize the threshold level of between-group selection dividing a regime
in which the population converges to a delta function at the equilibrium of the
within-group dynamics from a regime in which between-group competition
facilitates the existence of steady-state densities supporting greater levels
of cooperation. In particular, we see that the threshold selection strength and
average payoff at steady state depend on a tug-of-war between the
individual-level incentive to be a defector in a many-cooperator group and the
group-level incentive to have many cooperators over many defectors. We also
find that lower-level selection casts a long shadow: if groups are best off
with a mix of cooperators and defectors, then there will always be fewer
cooperators than optimal at steady state, even in the limit of infinitely
strong competition between groups.


  We investigate the relation between the stationary probability distribution
of chemical reaction systems and the convective field derived from the chemical
Fokker-Planck equation (CFPE) by comparing predictions of the convective field
to the results of stochastic simulations based on Gillespie's algorithm. The
convective field takes into account the drift term of the CFPE and the reaction
bias introduced by the diffusion term. For one-dimensional systems, fixed
points and bifurcations of the convective field correspond to extrema and
phenomenological bifurcations of the stationary probability distribution
whenever the CFPE is a good approximation to the stochastic dynamics. This
provides an efficient way to calculate the effect of system size on the number
and location of probability maxima and their phenomenological bifurcations in
parameter space. For two-dimensional systems, we study models that have
saddle-node and Hopf bifurcations in the macroscopic limit. Here, the existence
of two stable fixed points of the convective field correlates either with two
peaks of the stationary probability distribution, or with a peak and a
shoulder. In contrast, a Hopf bifurcation that occurs in the convective field
for decreasing system size is not accompanied by the onset of a crater-shaped
probability distribution; decreasing system size rather destroys craters and
replaces them by local maxima.


  In this study, first we expand the cubic population model under Allee effect
with a quadratic harvest function that represents harvest effect. Then, using
four reaction equations representing the micro-interactions within the
population under the influence of demographic noise and harvest, we obtain the
mean field equations containing the effect of the harvest from the solution of
the master equation. In this way, a relationship is established between the
micro and macro parameters of the population. As a result, we calculate Mean
Time to Extinction (MTE) by using WKB approximation for single step processes
and observe the effect of harvesting.


  Recent clinical studies have shown that HIV disease pathogenesis can depend
strongly on many factors at the time of transmission, including the strength of
the initial viral load and the local availability of CD4+ T-cells. In this
article, a new within-host model of HIV infection that incorporates the
homeostatic proliferation of T-cells is formulated and analyzed. Due to the
effects of this biological process, the influence of initial conditions on the
proliferation of HIV infection is further elucidated. The identifiability of
parameters within the model is investigated and a local stability analysis,
which displays additional complexity in comparison to previous models, is
conducted. The current study extends previous theoretical and computational
work on the early stages of the disease and leads to interesting nonlinear
dynamics, including a parameter region featuring bistability of infectious and
viral clearance equilibria and the appearance of a Hopf bifurcation within
biologically relevant parameter regimes.


  In biomedical applications of machine learning, relevant information often
has a rich structure that is not easily encoded as real-valued predictors.
Examples of such data include DNA or RNA sequences, gene sets or pathways, gene
interaction or coexpression networks, ontologies, and phylogenetic trees. We
highlight recent examples of machine learning models that use structure to
constrain model architecture or incorporate structured data into model
training. For machine learning in biomedicine, where sample size is limited and
model interpretability is critical, incorporating prior knowledge in the form
of structured data can be particularly useful. The area of research would
benefit from performant open source implementations and independent
benchmarking efforts.


  We consider a stochastic susceptible-infected-recovered (SIR) model with
contact tracing on random trees and on the configuration model. On a rooted
tree, where initially all individuals are susceptible apart from the root which
is infected, we are able to find exact formulas for the distribution of the
infectious period. Thereto, we show how to extend the existing theory for
contact tracing in homogeneously mixing populations to trees. Based on these
formulas, we discuss the influence of randomness in the tree and the basic
reproduction. We find the well known results for the homogeneously mixing case
as a limit of the present model (tree-shaped contact graph). Furthermore, we
develop approximate mean field equations for the dynamics on trees, and - using
the message passing method - also for the configuration model. The
interpretation and implications of the results are discussed.


  The research reported in this paper identifies the epigenetic biomarker
(methylation beta pattern) of breast cancer. Many cancers are triggered by
abnormal gene expression levels caused by aberrant methylation of CpG sites in
the DNA. In order to develop early diagnostics of cancer-causing methylations
and to develop a treatment, it is necessary to identify a few dozen key
cancer-related CpG methylation sites out of the millions of locations in the
DNA. This research used public TCGA dataset to train a TensorFlow machine
learning model to classify breast cancer versus non-breast-cancer tissue
samples, based on over 300,000 methylation beta values in each sample. L1
regularization was applied to identify the CpG methylation sites most important
for accurate classification. It was hypothesized that CpG sites with the
highest learned model weights correspond to DNA locations most relevant to
breast cancer. A reduced model trained on methylation betas of just the 25 CpG
sites having the highest weights in the full model (trained on methylation
betas at over 300,000 CpG sites) has achieved over 94% accuracy on evaluation
data, confirming that the identified 25 CpG sites are indeed a biomarker of
breast cancer.