  Understanding the emergence of biodiversity patterns in nature is a central
problem in biology. Theoretical models of speciation have addressed this
question in the macroecological scale, but little has been investigated in the
macroevolutionary context. Knowledge of the evolutionary history allows the
study of patterns underlying the processes considered in these models,
revealing their signatures and the role of speciation and extinction in shaping
macroevolutionary patterns. In this paper we introduce two algorithms to record
the evolutionary history of populations in individual-based models of
speciation, from which genealogies and phylogenies can be constructed. The
first algorithm relies on saving ancestral-descendant relationships, generating
a matrix that contains the times to the most recent common ancestor between all
pairs of individuals at every generation (the Most Recent Common Ancestor Time
matrix, MRCAT). The second algorithm directly records all speciation and
extinction events throughout the evolutionary process, generating a matrix with
the true phylogeny of species (the Sequential Speciation and Extinction Events,
SSEE). We illustrate the use of these algorithms in a spatially explicit
individual-based model of speciation. We compare the trees generated via MRCAT
and SSEE algorithms with trees inferred by methods that use only genetic
distance among extant species, commonly used in empirical studies and applied
here to simulated genetic data. Comparisons between tress are performed with
metrics describing the overall topology, branch length distribution and
imbalance of trees. We observe that both MRCAT and distance-based trees differ
from the true phylogeny, with the first being closer to the true tree than the
second.


  Phylogenetic trees illustrate the evolutionary history of genes and species.
In most cases, although genes evolve along with the species they belong to, a
species tree and gene tree are not identical, because of evolutionary events at
the gene level like duplication or transfer. These differences are handled by
phylogenetic reconciliation, which formally is a mapping between gene tree
nodes and species tree nodes and branches. We investigate models of
reconciliation with a gene transfer that replaces existing gene, which is a
biological important event but never included in reconciliation models. Also
the problem is close to a dated version of the classical subtree prune and
regraft (SPR) distance problem, where a pruned subtree has to be regrafted only
on a branch closer to the root. We prove that the reconciliation problem
including transfer and replacement is NP-hard, and that if speciations and
transfers with replacement are the only allowed evolutionary events, then it is
fixed-parameter tractable (FPT) with respect to the reconciliation's weight. We
prove that the results extend to the dated SPR problem.


  We give a non-technical introduction to convergence-divergence models, a new
modeling approach for phylogenetic data that allows for the usual divergence of
species post speciation but also allows for species to converge, i.e. become
more similar over time. By examining the $3$-taxon case in some detail we
illustrate that phylogeneticists have been "spoiled" in the sense of not having
to think about the structural parameters in their models by virtue of the
strong assumption that evolution is treelike. We show that there are not always
good statistical reasons to prefer the usual class of treelike models over more
general convergence-divergence models. Specifically we show many $3$-taxon
datasets can be equally well explained by supposing violation of the molecular
clock due to change in the rate of evolution along different edges, or by
keeping the assumption of a constant rate of evolution but instead assuming
that evolution is not a purely divergent process. Given the abundance of
evidence that evolution is not strictly treelike, our discussion is an
illustration that as phylogeneticists we often need to think clearly about the
structural form of the models we use.


  Biological network alignment (NA) aims to identify similar regions between
molecular networks of different species. NA can be local or global. Just as the
recent trend in the NA field, we also focus on global NA, which can be pairwise
(PNA) and multiple (MNA). PNA produces aligned node pairs between two networks.
MNA produces aligned node clusters between more than two networks. Recently,
the focus has shifted from PNA to MNA, because MNA captures conserved regions
between more networks than PNA (and MNA is thus considered to be more
insightful), though at higher computational complexity. The issue is that, due
to the different outputs of PNA and MNA, a PNA method is only compared to other
PNA methods, and an MNA method is only compared to other MNA methods.
Comparison of PNA against MNA must be done to evaluate whether MNA's higher
complexity is justified by its higher accuracy. We introduce a framework that
allows for this. We evaluate eight prominent PNA and MNA methods, on synthetic
and real-world biological networks, using topological and functional alignment
quality measures. We compare PNA against MNA in both a pairwise (native to PNA)
and multiple (native to MNA) manner. PNA is expected to perform better under
the pairwise evaluation framework. Indeed this is what we find. MNA is expected
to perform better under the multiple evaluation framework. Shockingly, we find
this not to always hold; PNA is often better than MNA in this framework,
depending on the choice of evaluation test.


  The theory of evolution states that the diversity of species can be explained
by descent with modification. Therefore, all living beings are related through
a common ancestor. This evolutionary process must have left traces in our
molecular composition. In this work, we present a randomization procedure in
order to determine if a group of 5 species of the primate family, namely,
macaque, guereza, orangutan, chimpanzee and human, has retained these traces in
its molecules. Firstly, we present the randomization methodology through two
toy examples, which allow to understand its logic. We then carry out a DNA data
analysis to assess if the group of primates contains phylogenetic information
which links them in a joint evolutionary history. This is carried out by
monitoring a Bayesian measure, called marginal likelihood, which we estimate by
using nested sampling. We found that it would be unusual to get the
relationship observed in the data among these primate species if they had not
shared a common ancestor. The results are in total agreement with the theory of
evolution.


  Phylogenetic comparative methods for real-valued traits usually make use of
stochastic process whose trajectories are continuous. This is despite
biological intuition that evolution is rather punctuated than gradual. On the
other hand, there has been a number of recent proposals of evolutionary models
with jump components. However, as we are only beginning to understand the
behaviour of branching Ornstein-Uhlenbeck (OU) processes the asymptotics of
branching OU processes with jumps is an even greater unknown. In this work we
build up on a previous study concerning OU with jumps evolution on a pure birth
tree. We introduce an extinction component and explore via simulations, its
effects on the weak convergence of such a process. We furthermore, also use
this work to illustrate the simulation and graphic generation possibilities of
the mvSLOUCH package.


  Models of codon evolution are commonly used to identify positive selection.
Positive selection is typically a heterogeneous process, i.e., it acts on some
branches of the evolutionary tree and not others. Previous work on DNA models
showed that when evolution occurs under a heterogeneous process it is important
to consider the property of model closure, because non-closed models can give
biased estimates of evolutionary processes. The existing codon models that
account for the genetic code are not closed; to establish this it is enough to
show that they are not linear (meaning that the sum of two codon rate matrices
in the model is not a matrix in the model). This raises the concern that a
single codon model fit to a heterogeneous process might mis-estimate both the
effect of selection and branch lengths.
  Codon models are typically constructed by choosing an underlying DNA model
(e.g., HKY) that acts identically and independently at each codon position, and
then applying the genetic code via the parameter $\omega$ to modify the rate of
transitions between codons that code for different amino acids. Here we use
simulation to investigate the accuracy of estimation of both the selection
parameter $\omega$ and branch lengths in cases where the underlying DNA process
is heterogeneous but $\omega$ is constant. We find that both $\omega$ and
branch lengths can be mis-estimated in these scenarios. Errors in $\omega$ were
usually less than 2% but could be as high as 17%. We also assessed if choosing
different underlying DNA models had any affect on accuracy, in particular we
assessed if using closed DNA models gave any advantage. However, a DNA model
being closed does not imply that the codon model constructed from it is closed,
and in general we found that using closed DNA models did not decrease errors in
the estimation of $\omega$.


  The efficiency of intracellular cargo transport from specific source to
target locations is strongly dependent upon molecular motor-assisted motion
along the cytoskeleton. Radial transport along microtubules and lateral
transport along the filaments of the actin cortex underneath the cell membrane
are characteristic for cells with a centrosome. The interplay between the
specific cytoskeleton organization and the motor performance realizes a
spatially inhomogeneous intermittent search strategy. In order to analyze the
efficiency of such intracellular search strategies we formulate a random
velocity model with intermittent arrest states. We evaluate efficiency in terms
of mean first passage times for three different, frequently encountered
intracellular transport tasks: i) the narrow escape problem, which emerges
during cargo transport to a synapse or other specific region of the cell
membrane, ii) the reaction problem, which considers the binding time of two
particles within the cell, and iii) the reaction-escape problem, which arises
when cargo must be released at a synapse only after pairing with another
particle. Our results indicate that cells are able to realize efficient search
strategies for various intracellular transport tasks economically through a
spatial cytoskeleton organization that involves only a narrow actin cortex
rather than a cell body filled with randomly oriented actin filaments.


  Helical molecules change their twist number under the effect of a mechanical
load. We study the twist-stretch relation for a set of short DNA molecules
modeled by a mesoscopic Hamiltonian. Finite temperature path integral
techniques are applied to generate a large ensemble of possible configurations
for the base pairs of the sequence. The model also accounts for the bending and
twisting fluctuations between adjacent base pairs along the molecules stack.
Simulating a broad range of twisting conformation, we compute the helix
structural parameters by averaging over the ensemble of base pairs
configurations. The method selects, for any applied force, the average twist
angle which minimizes the molecule's free energy. It is found that the chains
generally over-twist under an applied stretching and the over-twisting is
physically associated to the contraction of the average helix diameter, i.e. to
the damping of the base pair fluctuations. Instead, assuming that the maximum
amplitude of the bending fluctuations may decrease against the external load,
the DNA molecule first over-twists for weak applied forces and then untwists
above a characteristic force value. Our results are discussed in relation to
available experimental information albeit for kilo-base long molecules.


  Improvements in sequencing technologies and reduced experimental costs have
resulted in a vast number of studies generating high-throughput data. Although
the number of methods to analyze these "omics" data has also increased,
computational complexity and lack of documentation hinder researchers from
analyzing their high-throughput data to its true potential. In this chapter we
detail our data-driven, transkingdom network (TransNet) analysis protocol to
integrate and interrogate multi-omics data. This systems biology approach has
allowed us to successfully identify important causal relationships between
different taxonomic kingdoms (e.g. mammals and microbes) using diverse types of
data.


  More and more evidence shows that mating preference is a mechanism that may
lead to a reproductive isolation event. In this paper, a haploid population
living on two patches linked by migration is considered. Individuals are
ecologically and demographically neutral on the space and differ only on a
trait, $a$ or $A$, affecting both mating success and migration rate. The
special feature of this paper is to assume that the strengths of the mating
preference and the migration depend on the trait carried. Indeed, patterns of
mating preferences are generally asymmetrical between the subspecies of a
population. I prove that mating preference interacting with frequency-dependent
migration behavior can lead to a reproductive isolation. Then, I describe the
time before reproductive isolation occurs. To reach this result, I use an
original method to study the limiting dynamical system, analyzing first the
system without migration and adding migration with a perturbation method.
Finally, I study how the time before reproductive isolation is influenced by
the parameters of migration and of mating preferences, highlighting that large
migration rates tend to favor types with weak mating preferences.


  Crohn's disease (CD) is associated with an ecological imbalance of the
intestinal microbiota, consisting of hundreds of species. The underlying
complexity as well as individual differences between patients contributes to
the difficulty to define a standardized treatment. Computational modeling can
systematically investigate metabolic interactions between gut microbes to
unravel novel mechanistic insights. In this study, we integrated metagenomic
data of CD patients and healthy controls with genome-scale metabolic models
into personalized in silico microbiotas. We predicted short chain fatty acid
(SFCA) levels for patients and controls, which were overall congruent with
experimental findings. As an emergent property, low concentrations of SCFA were
predicted for CD patients and the SCFA signatures were unique to each patient.
Consequently, we suggest personalized dietary treatments that could improve
each patient's SCFA levels. The underlying modeling approach could aid clinical
practice to find novel dietary treatment and guide recovery by rationally
proposing food aliments.


  Starting from the detailed catalytic mechanism of a biocatalyst we provide a
coarse-graining procedure which, by construction, is thermodynamically
consistent. This procedure provides stoichiometries, reaction fluxes (rate
laws), and reaction forces (Gibbs energies of reaction) for the coarse-grained
level. It can treat active transporters and molecular machines, and thus
extends the applicability of ideas that originated in enzyme kinetics. Our
results lay the foundations for systematic studies of the thermodynamics of
large-scale biochemical reaction networks. Moreover, we identify the conditions
under which a relation between one-way fluxes and forces holds at the
coarse-grained level as it holds at the detailed level. In doing so, we clarify
the speculations and broad claims made in the literature about such a general
flux--force relation. As a further consequence we show that, in contrast to
common belief, the second law of thermodynamics does not require the currents
and the forces of biochemical reaction networks to be always aligned.


  The Spliced Alignment Problem (SAP) that consists in finding an optimal
semi-global alignment of a spliced RNA sequence on an unspliced genomic
sequence has been largely considered for the prediction and the annotation of
gene structures in genomes. Here, we re-visit it for the purpose of identifying
CDS ortholog groups within a set of CDS from homologous genes and for computing
multiple CDS alignments. We introduce a new constrained version of the spliced
alignment problem together with an algorithm that exploits full information on
the exon-intron structure of the input RNA and gene sequences in order to
compute high-coverage accurate alignments. We show how pairwise spliced
alignments between the CDS and the gene sequences of a gene family can be
directly used in order to clusterize the set of CDS of the gene family into a
set of ortholog groups. We also introduce an extension of the spliced alignment
problem called Multiple Spliced Alignment Problem (MSAP) that consists in
aligning simultaneously several RNA sequences on several genes from the same
gene family. We develop a heuristic algorithmic solution for the problem. We
show how to exploit multiple spliced alignments for the clustering of
homologous CDS into ortholog and close paralog groups, and for the construction
of multiple CDS alignments. An implementation of the method in Python is
available on demande to SFA@USherbrooke.ca.
  Keywords: Spliced alignment, CDS ortholog groups, Multiple CDS alignment,
Gene structure, Gene family.


  This work consists of an epidemic model with vaccination coupled with an
opinion dynamics. Our objective was to study how disease risk perception can
influence opinions about vaccination and therefore the spreading of the
disease. Differently from previous works we have considered continuous
opinions. The epidemic spreading is governed by a SIS-like model with an extra
vaccinated state. In our model individuals vaccinate with a probability
proportional to their opinions. The opinions change due to peer influence in
pairwise interactions. The epidemic feedback to the opinion dynamics acts as an
external field increasing the vaccination probability. We performed Monte Carlo
simulations in fully-connected populations. Interestingly we observed the
emergence of a first-order phase transition, besides the usual active-absorbing
phase transition presented in the SIS model. Our simulations also show that
with a certain combination of parameters, an increment in the initial fraction
of the population that is pro-vaccine has a twofold effect: it can lead to
smaller epidemic outbreaks in the short term, but it also contributes to the
survival of the chain of infections in the long term. Our results also suggest
that it is possible that more effective vaccines can decrease the long-term
vaccine coverage. This is a counterintuitive outcome, but it is in line with
empirical observations that vaccines can become a victim of their own success.


  We construct genomic predictors for heritable and extremely complex human
quantitative traits (height, heel bone density, and educational attainment)
using modern methods in high dimensional statistics (i.e., machine learning).
Replication tests show that these predictors capture, respectively, $\sim$40,
20, and 9 percent of total variance for the three traits. For example,
predicted heights correlate $\sim$0.65 with actual height; actual heights of
most individuals in validation samples are within a few cm of the prediction.
The variance captured for height is comparable to the estimated SNP
heritability from GCTA (GREML) analysis, and seems to be close to its
asymptotic value (i.e., as sample size goes to infinity), suggesting that we
have captured most of the heritability for the SNPs used. Thus, our results
resolve the common SNP portion of the "missing heritability" problem -- i.e.,
the gap between prediction R-squared and SNP heritability. The $\sim$20k
activated SNPs in our height predictor reveal the genetic architecture of human
height, at least for common SNPs. Our primary dataset is the UK Biobank cohort,
comprised of almost 500k individual genotypes with multiple phenotypes. We also
use other datasets and SNPs found in earlier GWAS for out-of-sample validation
of our results.


  Cells process external and internal signals through chemical interactions.
Cells that constitute the immune system (e.g., antigen presenting cell, T-cell,
B-cell, mast cell) can have different functions (e.g., adaptive memory,
inflammatory response) depending on the type and number of receptor molecules
on the cell surface and the specific intracellular signaling pathways activated
by those receptors. Explicitly modeling and simulating kinetic interactions
between molecules allows us to pose questions about the dynamics of a signaling
network under various conditions. However, the application of chemical kinetics
to biochemical signaling systems has been limited by the complexity of the
systems under consideration. Rule-based modeling (BioNetGen, Kappa, Simmune,
PySB) is an approach to address this complexity. In this chapter, by
application to the Fc$\varepsilon$RI receptor system, we will explore the
origins of complexity in macromolecular interactions, show how rule-based
modeling can be used to address complexity, and demonstrate how to build a
model in the BioNetGen framework. Open source BioNetGen software and
documentation are available at http://bionetgen.org.


  Introduction: Identification of blood-based metabolic changes might provide
early and easy-to-obtain biomarkers.
  Methods: We included 127 AD patients and 121 controls with
CSF-biomarker-confirmed diagnosis (cut-off tau/A$\beta_{42}$: 0.52). Mass
spectrometry platforms determined the concentrations of 53 amine, 22 organic
acid, 120 lipid, and 40 oxidative stress compounds. Multiple signatures were
assessed: differential expression (nested linear models), classification
(logistic regression), and regulatory (network extraction).
  Results: Twenty-six metabolites were differentially expressed. Metabolites
improved the classification performance of clinical variables from 74% to 79%.
Network models identified 5 hubs of metabolic dysregulation: Tyrosine,
glycylglycine, glutamine, lysophosphatic acid C18:2 and platelet activating
factor C16:0. The metabolite network for APOE $\epsilon$4 negative AD patients
was less cohesive compared to the network for APOE $\epsilon$4 positive AD
patients.
  Discussion: Multiple signatures point to various promising peripheral markers
for further validation. The network differences in AD patients according to
APOE genotype may reflect different pathways to AD.


  Background: The 2014 Ebola outbreak in West Africa was the largest on record,
resulting in over 25,000 total infections and 15,000 total deaths. Mathematical
modeling can be used to investigate the mechanisms driving transmission during
this outbreak -- in particular, burial practices appear to have been major
source of infections. Methodology/Principal Findings: We developed a
multi-stage model of Ebola virus transmission linked to a game-theoretic model
of population burial practice selection. We fit our model to cumulative
incidence and mortality data from Guinea, Liberia, and Sierra Leone from
January 2014 to March, 2016. The inclusion of behavior change substantially
improved best fit estimates and final size prediction compared to a reduced
model with fixed burials. Best fit trajectories suggest that the majority of
sanitary burial adoption occurred between July, 2014 and October, 2014.
However, these simulations also indicated that continued sanitary burial
practices waned following the resolution of the outbreak.
Conclusions/Significance: Surveillance data from the 2014 outbreak appears to
have a signal of changes in the dominant burial practices in all three
countries. Increased adoption of sanitary burials likely attenuated
transmission, but these changes occurred too late to prevent the explosive
growth of the outbreak during its early phase. For future outbreaks, explicitly
modeling behavior change and collecting data on transmission-related behaviors
may improve intervention planning and outbreak response.


  A major part of the interactions involved in the assembly and stability of
icosahedral, positive-sense single-stranded RNA (ssRNA+) viruses is
electrostatic in nature, as can be inferred from the strong $pH$- and
salt-dependence of their assembly phase diagrams. Electrostatic interactions do
not act only between the capsid coat proteins (CPs), but just as often provide
a significant contribution to the interactions of the CPs with the genomic RNA,
mediated to a large extent by positively charged, flexible N-terminal tails of
the CPs. In this work, we provide two clear and complementary definitions of an
N-terminal tail of a protein, and use them to extract the tail sequences of a
large number of CPs of ssRNA+ viruses. We examine the $pH$-dependent interplay
of charge on both tails and CPs alike, and show that -- in contrast to the
charge on the CPs -- the net positive charge on the N-tails persists even to
very basic $pH$ values. In addition, we note a limit to the length of the
wild-type genomes of those viruses which utilize positively charged tails, when
compared to viruses without charged tails and similar capsid size. At the same
time, we observe no clear connection between the charge on the N-tails and the
genome lengths of the viruses included in our study.


  The spatial dispersal of individuals is known to play an important role in
the dynamics of populations, and is central to metapopulation theory. At the
same time, local adaptation to environmental conditions creates a geographic
mosaic of evolutionary forces, where the combined drivers of selection and gene
flow interact. Although the dispersal of individuals from donor to recipient
populations provides connections within the metapopulation, promoting
demographic and evolutionary rescue, it may also introduce maladapted
individuals into habitats host to different environmental conditions,
potentially lowering the fitness of the recipient population. Here we explore a
model of the eco-evolutionary dynamics between two populations connected by
dispersal, where the productivity of each is defined by a trait complex that is
subject to local selection. Although general in nature, our model is inspired
by salmon metapopulations, where dispersal between populations is defined in
terms of the straying rate, which has been shown to be density-dependent. The
results of our model reveal that increased straying between evolving
populations leads to alternative stable states, which has large and nonlinear
effects on two measures of metapopulation robustness: the portfolio effect and
the time to recovery following an induced disturbance. We show that
intermediate levels of straying result in large gains in robustness, and that
increased habitat heterogeneity promotes robustness when straying rates are
low, and erodes robustness when straying rates are high. Finally, we show that
density-dependent straying promotes robustness, particularly when the aggregate
biomass is low and straying is correspondingly high, which has important
ramifications for the conservation of salmon metapopulations facing both
natural and anthropogenic disturbances.


  Coalescence processes have received a lot of attention in the context of
conditional branching processes with fixed population size and non-overlapping
generations. Here we focus on similar problems in the context of the standard
unconditional Bienaym\'e-Galton-Watson branching processes, either
(sub)-critical or supercritical. Using an analytical tool, we derive the
structure of some counting aspects of the ancestral genealogy of such
processes, including: the transition matrix of the ancestral count process and
an integral representation of various coalescence times distributions, such as
the time to most recent common ancestor of a random sample of arbitrary size,
including full size.
  We illustrate our results on two important examples of branching mechanisms
displaying either finite or infinite reproduction mean, their main interest
being to offer a closed form expression for their probability generating
functions at all times. Large time behaviors are investigated.


  We define and analyze a coalescent process as a recursive box-filling process
whose genealogy is given by an ancestral time-reversed, time-inhomogeneous
Bienyam\'{e}-Galton-Watson process. Special interest is on the expected size of
a typical box and its probability of being empty. Special cases leading to
exact asymptotic computations are investigated when the coalescing mechanisms
are either linear-fractional or quadratic.


  The success of a vaccination program is crucially dependent on its adoption
by a critical fraction of the population, as the resulting herd immunity
prevents future outbreaks of an epidemic. However, the effectiveness of a
campaign can engender its own undoing if individuals choose to not get
vaccinated in the belief that they are protected by herd immunity. Although
this may appear to be an optimal decision, based on a rational appraisal of
cost and benefits to the individual, it exposes the population to subsequent
outbreaks. We investigate if voluntary vaccination can emerge in a an
integrated model of an epidemic spreading on a social network of rational
agents that make informed decisions whether to be vaccinated. The information
available to each agent includes the prevalence of the disease in their local
network neighborhood and/or globally in the population, as well as the fraction
of their neighbors that are protected against the disease. Crucially, the
payoffs governing the decision of agents evolve with disease prevalence,
resulting in the co-evolution of vaccine uptake behavior with the spread of the
contagion. The collective behavior of the agents responding to local prevalence
can lead to a significant reduction in the final epidemic size, particularly
for less contagious diseases having low basic reproduction number $R_0$. Near
the epidemic threshold ($R_0\approx1$) the use of local prevalence information
can result in a dichotomous response in final vaccine coverage. The
implications of our results suggest the nature of information used by
individuals is a critical factor determining the success of public health
intervention schemes that involve mass vaccination.


  Masses of experiments have shown individual preference for fairness which
seems irrational. The reason behind it remains a focus for research. The effect
of spite (individuals are only concerned with their own relative standing) on
the evolution of fairness has attracted increasing attention from experiments,
but only has been implicitly studied in one evolutionary model. The model did
not involve high-offer rejections, which have been found in the form of
non-monotonic rejections (rejecting offers that are too high or too low) in
experiments. Here, we introduce a high offer and a non-monotonic rejection in
structured populations of finite size, and use strategy intervention to
explicitly study how spite influences the evolution of fairness: five
strategies are in sequence added into the competition of a fair strategy and a
selfish strategy. We find that spite promotes fairness, altruism inhibits
fairness, and the non-monotonic rejection can cause fairness to overcome
selfishness, which cannot happen without high-offer rejections. Particularly
for the group-structured population with seven discrete strategies, we
analytically study the effect of population size, mutation, and migration on
fairness, selfishness, altruism, and spite. A larger population size cannot
change the dominance of fairness, but it promotes altruism and inhibits
selfishness and spite. Intermediate mutation maximizes selfishness and
fairness, and minimizes spite; intermediate mutation maximizes altruism for
intermediate migration and minimizes altruism otherwise. The existence of
migration inhibits selfishness and fairness, and promotes altruism; sufficient
migration promotes spite. Our study may provide important insights into the
evolutionary origin of fairness.


  Predicting three dimensional residue-residue contacts from evolutionary
information in protein sequences was attempted already in the early 1990s.
However, contact prediction accuracies of methods evaluated in CASP experiments
before CASP11 remained quite low, typically with $<20$% true positives.
Recently, contact prediction has been significantly improved to the level that
an accurate three dimensional model of a large protein can be generated on the
basis of predicted contacts. This improvement was attained by disentangling
direct from indirect correlations in amino acid covariations or cosubstitutions
between sites in protein evolution. Here, we review statistical methods for
extracting causative correlations and various approaches to describe protein
structure, complex, and flexibility based on predicted contacts.


  A numerical method using implicit surface representations is proposed to
solve the linearized Poisson-Boltzmann equations that arise in mathematical
models for the electrostatics of molecules in solvent. The proposed method used
an implicit boundary integral formulation to derived a linear system defined on
Cartesian nodes in a narrowband surrounding the closed surface that separate
the molecule and the solvent. The needed implicit surfaces is constructed from
the given atomic description of the molecules, by a sequence of standard level
set algorithms. A fast multipole method is applied to accelerate the solution
of the linear system. A few numerical studies involving some standard test
cases are presented and compared to other existing results.


  An RNA secondary structure is designable if there is an RNA sequence which
can attain its maximum number of base pairs only by adopting that structure.
The combinatorial RNA design problem, introduced by Hale\v{s} et al. in 2016,
is to determine whether or not a given RNA secondary structure is designable.
Hale\v{s} et al. identified certain classes of designable and non-designable
secondary structures by reference to their corresponding rooted trees. We
introduce an infinite class of rooted trees containing unpaired nucleotides at
the greatest depth, and prove constructively that their corresponding secondary
structures are designable. This complements previous results for the
combinatorial RNA design problem.


  Living cells move thanks to assemblies of actin filaments and myosin motors
that range from very organized striated muscle tissue to disordered
intracellular bundles. The mechanisms powering these disordered structures are
debated, and all models studied so far predict that they are contractile. We
reexamine this prediction through a theoretical treatment of the interplay of
three well-characterized internal dynamical processes in actomyosin bundles:
actin treadmilling, the attachement-detachment dynamics of myosin and that of
crosslinking proteins. We show that these processes enable an extensive control
of the bundle's active mechanics, including reversals of the filaments'
apparent velocities and the possibility of generating extension instead of
contraction. These effects offer a new perspective on well-studied in vivo
systems, as well as a robust criterion to experimentally elucidate the
underpinnings of actomyosin activity.


  Single molecule time traces reveal the time evolution of unsynchronized
kinetic systems. Especially single molecule F\"orster resonance energy transfer
(smFRET) provides access to enzymatically important timescales, combined with
molecular distance resolution and minimal interference with the sample. Yet the
kinetic analysis of smFRET time traces is complicated by experimental
shortcomings - such as photo-bleaching and noise. Here we recapitulate the
fundamental limits of single molecule fluorescence that render the classic,
dwell-time based kinetic analysis unsuitable. In contrast, our Single Molecule
Analysis of Complex Kinetic Sequences (SMACKS) considers every data point and
combines the information of many short traces in one global kinetic rate model.
We demonstrate the potential of SMACKS by resolving the small kinetic effects
caused by different ionic strengths in the chaperone protein Hsp90. These
results show an unexpected interrelation between conformational dynamics and
ATPase activity in Hsp90.


  Local interactions, when individuals meet, can regulate collective behavior.
In a system without any central control, the rate of interaction may depend
simply on how the individuals move around. But interactions could in turn
influence movement; individuals might seek out interactions, or their movement
in response to interaction could influence further interaction rates. We
develop a general framework to address these questions, using collision theory
to establish a baseline expected rate of interaction based on proximity. We
test the models using data from harvester ant colonies. A colony uses feedback
from interactions inside the nest to regulate foraging activity. Potential
foragers leave the nest in response to interactions with returning foragers
with food. The time series of interactions and local density of ants show how
density hotspots lead to interactions that are clustered in time. A correlated
random walk null model describes the mixing of potential and returning
foragers. A model from collision theory relates walking speed and spatial
proximity with the probability of interaction. The results demonstrate that
although ants do not mix homogeneously, trends in interaction patterns can be
explained simply by the walking speed and local density of surrounding ants.


  Despite the huge interest in network resilience to stress, most of the
studies have concentrated on internal stress damaging network structure (e.g.,
node removals). Here we study how networks respond to environmental stress
deteriorating their external conditions. We show that, when regular networks
gradually disintegrate as environmental stress increases, disordered networks
can suddenly collapse at critical stress with hysteresis and vulnerability to
perturbations. We demonstrate that this difference results from a trade-off
between node resilience and network resilience to environmental stress. The
nodes in the disordered networks can suppress their collapses due to the
small-world topology of the networks but eventually collapse all together in
return. Our findings indicate that some real networks can be highly resilient
against environmental stress to a threshold yet extremely vulnerable to the
stress above the threshold because of their small-world topology.


  The Bak Sneppen (BS) model is a very simple model that exhibits all the
richness of self-organized criticality theory. At the thermodynamic limit, the
BS model converges to a situation where all particles have a fitness that is
uniformly distributed between a critical value $p_c$ and 1. The $p_c$ value is
unknown, as are the variables that influence and determine this value. Here, we
study the Bak Sneppen model in the case in which the lowest fitness particle
interacts with an arbitrary even number of $m$ nearest neighbors. We show that
$p_{c,m}$ verifies a simple local equilibrium relationship. Based on this
relationship, we can determine bounds for $p_{c,m}$.


  In the Bak-Sneppen model, the lowest fitness particle and its two nearest
neighbors are renewed at each temporal step with a uniform (0,1) fitness
distribution. The model presents a critical value that depends on the
interaction criteria (two nearest neighbors) and on the update procedure
(uniform). Here we calculate the critical value for models where one or both
properties are changed. We study models with non-uniform updates, models with
random neighbors and models with binary fitness and obtain exact results for
the average fitness and for $p_c$.


  Big data initiatives such as the Enhancing NeuroImaging Genetics through
Meta-Analysis consortium (ENIGMA), combine data collected by independent
studies worldwide to achieve more accurate estimates of effect sizes and more
reliable and reproducible outcomes. Such efforts require harmonized analyses
protocols to consistently extract phenotypes. Even so, challenges include wide
variability of fMRI protocols and scanner platforms; this leads to site-to-site
variance in quality, resolution and temporal signal-to-noise ratio (tSNR). An
effective harmonization should provide optimal measures for data of different
qualities. We developed a multi-site rsfMRI analysis pipeline to allow research
groups around the world to process rsfMRI scans in a harmonized way, to extract
consistent and quantitative measurements of connectivity and to perform
coordinated statistical tests. We used the single-modality ENIGMA rsfMRI
pipeline based on model-free Marchenko-Pastor PCA based denoising to verify and
replicate findings of significant heritability of measures from resting state
networks. We analyzed two independent cohorts, GOBS (Genetics of Brain
Structure) and HCP (the Human Connectome Project), which collected data using
conventional and connectomics oriented fMRI protocols. We used seed-based
connectivity and dual-regression approaches to show that rsfMRI signal is
consistently heritable across twenty major functional network measures.
Heritability values of 20-40% were observed across both cohorts.


  We analyze the coalescing model where a `primary' colony grows and randomly
emits secondary colonies that spread and eventually coalesce with it. This
model describes population proliferation in theoretical ecology, tumor growth
and is also of great interest for modeling the development of cities. Assuming
the primary colony to be always spherical of radius $r(t)$ and the emission
rate proportional to $r(t)^\theta$ where $\theta>0$, we derive the mean-field
equations governing the dynamics of the primary colony, calculate the scaling
exponents versus $\theta$ and compare our results with numerical simulations.
We then critically test the validity of the circular approximation and show
that it is sound for a constant emission rate ($\theta=0$). However, when the
emission rate is proportional to the perimeter, the circular approximation
breaks down and the roughness of the primary colony can not be discarded, thus
modifying the scaling exponents.


  Coarse-grained models can be of great help to address the problem of
structure prediction in nucleic acids. On one hand they can make the prediction
more efficient, while on the other hand, they can also help to identify the
essential degrees of freedom and interactions for the description of a number
of structures. With the aim to provide an all-atom representation in an
explicit solvent to the predictions of our SPlit and conQueR (SPQR)
coarse-grained model of RNA, we recently introduced a backmapping procedure
which enforces the predicted structure into an atomistic one by means of
steered Molecular Dynamics. These simulations minimize the ERMSD, a particular
metric which deals exclusively with the relative arrangement of nucleobases,
between the atomistic representation and the target structure. In this paper,
we explore the effects of this approach on the resulting interaction networks
and backbone conformations by applying it on a set of fragments using as a
target their native structure. We find that the geometry of the target
structures can be reliably recovered, with limitations in the regions with
unpaired bases such as bulges. In addition, we observe that the folding pathway
can also change depending on the parameters used in the definition of the ERMSD
and the use of other metrics such as the RMSD.


  High-throughput sequencing techniques such as metagenomic and
metatranscriptomic technologies allow cataloging of functional characteristics
of microbial community members as well as their taxonomic identity. Such
studies have found that a community's composition in terms of ecologically
relevant functional traits or guilds can be conserved more strictly across
varying settings than taxonomic composition is. I use a standard ecological
resource-consumer model to examine the dynamics of traits relevant to resource
consumption, and analyze determinants of functional composition. This model
demonstrates that interaction with essential resources can regulate the
community-wide abundance of ecologically relevant traits, keeping them at
consistent levels despite large changes in the abundances of the species
housing those traits in response to changes in the environment, and across
variation between communities in species composition. Functional composition is
shown to be able to track differences in environmental conditions faithfully
across differences in community composition. Mathematical conditions on
consumers' vital rates and functional responses sufficient to produce
conservation of functional community structure across taxonomic differences are
presented.


  Dendrograms are a way to represent evolutionary relationships between
organisms. Nowadays, these are inferred based on the comparison of genes or
protein sequences by taking into account their differences and similarities.
The genetic material of choice for the sequence alignments (all the genes or
sets of genes) results in distinct inferred dendrograms. In this work, we
evaluate differences between dendrograms reconstructed with different
methodologies and obtained for different sets of organisms chosen at random
from a much larger set. A statistical analysis is performed in order to
estimate the fluctuation between the results obtained from the different
methodologies. This analysis permit us to validate a systematic approach, based
on the comparison of the organisms' metabolic networks for inferring
dendrograms. It has the advantage that it allows the comparison of organisms
very far away in the evolutionary tree even if they have no known ortholog gene
in common.


  Phylogenomics has ushered in an age of discordance. Analyses often reveal
abundant discordances among phylogenies of different parts of genomes, as well
as incongruences between species trees obtained using different methods or data
partitions. Researchers are often left trying to make sense of such
incongruences. Interpretive ways of measuring and visualizing discordance are
needed, both among alternative species trees and gene trees, especially for
specific focal branches of a tree. Here, we introduce DiscoVista, a publicly
available tool that creates a suite of simple but interpretable visualizations.
DiscoVista helps quantify the amount of discordance and some of its potential
causes.


  The large variation in seed mass among species inspired a vast array of
theoretical and empirical research attempting to explain this variation. So
far, seed mass variation was investigated by two classes of studies: one class
focuses on species varying in seed mass within communities, while the second
focuses on variation between communities, most often with respect to resource
gradients. Here, we develop a model capable of simultaneously explaining
variation in seed mass within and between communities. The model describes
resource competition (for both soil and light resources) in annual communities
and incorporates two fundamental aspects: light asymmetry (higher light
acquisition per unit biomass for larger individuals) and growth allometry
(negative dependency of relative growth rate on plant biomass). Results show
that both factors are critical in determining patterns of seed mass variation.
In general, growth allometry increases the reproductive success of small-seeded
species while light asymmetry increases the reproductive success of
large-seeded species. Increasing availability of soil resources increases light
competition, thereby increasing the reproductive success of large-seeded
species and ultimately the community (weighted) mean seed mass. An unexpected
prediction of the model is that maximum variation in community seed mass (a
measure of functional diversity) occurs under intermediate levels of soil
resources. Extensions of the model incorporating size-dependent seed survival
and disturbance also show patterns consistent with empirical observations.
These overall results suggest that the mechanisms captured by the model are
important in determining patterns of species and functional diversity.


  Statistical estimation theory determines the optimal way of estimating
parameters of a fluctuating noisy signal. However, if the estimation is
performed on unreliable hardware, a sub-optimal estimation procedure can
outperform the previously optimal procedure. Here, we compare classes of
circadian clocks by viewing them as phase estimators for the periodic day-night
light signal. We find that continuous attractor-based free running clocks, such
as those found in the cyanobacterium Synechococcus elongatus and humans, are
nearly optimal phase estimators since their flat attractor directions
efficiently project out light intensity fluctuations due to weather patterns
(`external noise'). However, such flat directions also make these continuous
limit cycle attractors highly vulnerable to diffusive 'internal noise'. Given
such unreliable biochemical hardware, we find that point attractor-based damped
clocks, such as those found in a smaller cyanobacterium with low protein copy
number, Prochlorococcus marinus, outperform continuous attractor-based clocks.
By interpolating between the two types of clocks found in these organisms, we
demonstrate a family of biochemical phase estimation strategies that are best
suited to different relative strengths of external and internal noise.


  Diverse repertoires of hypervariable immunoglobulin receptors (TCR and BCR)
recognize antigens in the adaptive immune system. The development of
immunoglobulin receptor repertoire sequencing methods makes it possible to
perform repertoire-wide disease association studies of antigen receptor
sequences. We developed a statistical framework for associating receptors to
disease from only a small cohort of patients, with no need for a control
cohort. Our method successfully identifies previously validated Cytomegalovirus
and type 1 diabetes responsive receptors.


  This work analyses the transcriptional effects of some biochemical mechanisms
proposed in previous literature which attempts to explain the differential
spatial expression of Hedgehog target genes involved in Drosophila development.
Specifically, the expression of decapentaplegic and patched, genes whose
transcription is believed to be controlled by the activator and repressor forms
of the transcription factor Cubitus interruptus (Ci). This study is based on a
thermodynamic approach which provides binding equilibrium weighted average rate
expressions for genes controlled by transcription factors competing and
(possibly) cooperating for common binding sites, in the same way that Ci's
activator and repressor forms might do. These expressions are refined to
produce simpler equivalent formulae allowing their mathematical analysis.
Thanks to this, we can evaluate the correlation between several molecular
processes and biological features observed at tissular level. In particular, we
will focus on how high/low/differential affinity and null/total/partial
cooperation modify the activation/repression regions of the target genes or
provoke signal modulation.


  The presence of one or more species at some spatial locations but not others
is a central matter in ecology. This phenomenon is related to ecological
pattern formation. Nonlocal interactions can be considered as one of the
mechanisms causing such a phenomenon. We propose a single-species, continuous
time metapopulation model taking nonlocal interactions into account. Discrete
probability kernels are used to model these interactions in a patchy
environment. A linear stability analysis of the model shows that solutions to
this equation exhibit pattern formation if the dispersal rate of the species is
sufficiently small and the discrete interaction kernel satisfies certain
conditions. We numerically observe that traveling and stationary wave type
patterns arise near critical dispersal rate. We use weakly nonlinear analysis
to better understand the behavior of formed patterns. We show that observed
patterns arise through both supercritical and subcritical bifurcations from
spatially homogeneous steady state. Moreover, we observe that as the dispersal
rate decreases, amplitude of the patterns increases. For discontinuous
transitions to instability, we also show that there exists a threshold for the
amplitude of the initial condition, above which pattern formation is observed.


  Understanding how patterns and travelling waves form in chemical and
biological reaction-diffusion models is an area which has been widely
researched, yet is still experiencing fast development. Surprisingly enough, we
still do not have a clear understanding about all possible types of dynamical
regimes in classical reaction-diffusion models such as Lotka-Volterra
competition models with spatial dependence. In this work, we demonstrate some
new types of wave propagation and pattern formation in a classical three
species cyclic competition model with spatial diffusion, which have been so far
missed in the literature. These new patterns are characterised by a high
regularity in space, but are different from patterns previously known to exist
in reaction-diffusion models, and may have important applications in improving
our understanding of biological pattern formation and invasion theory. Finding
these new patterns is made technically possible by using an automatic adaptive
finite element method driven by a novel a posteriori error estimate which is
proven to provide a reliable bound for the error of the numerical method. We
demonstrate how this numerical framework allows us to easily explore the
dynamical patterns both in two and three spatial dimensions.


  The prisoner's dilemma describes a conflict between a pair of players, in
which defection is a dominant strategy whereas cooperation is collectively
optimal. The iterated version of the dilemma has been extensively studied to
understand the emergence of cooperation. In the evolutionary context, the
iterated prisoner's dilemma is often combined with population dynamics, in
which a more successful strategy replicates itself with a higher growth rate.
Here, we investigate the replicator dynamics of three representative
strategies, i.e., unconditional cooperation, unconditional defection, and
tit-for-tat, which prescribes reciprocal cooperation by mimicking the
opponent's previous move. Our finding is that the dynamics is self-dual in the
sense that it remains invariant when we apply time reversal and exchange the
fractions of unconditional cooperators and defectors in the population. The
duality implies that the fractions can be equalized by tit-for-tat players,
although unconditional cooperation is still dominated by defection.
Furthermore, we find that mutation among the strategies breaks the exact
duality in such a way that cooperation is more favored than defection, as long
as the cost-to-benefit ratio of cooperation is small.


  We propose a population model for HIV-TB co-infection dynamics by considering
treatments for HIV infection, active tuberculosis and co-infection. The HIV
only and TB only models are analyzed separately, as well as full model. The
basic reproduction numbers for TB ($\mathcal{R}_0^T$) and HIV
($\mathcal{R}_0^H$) and overall reproduction number for the system
$\mathcal{R}_0= \max\{\mathcal{R}_0^T, \mathcal{R}_0^H\}$ are computed. The
equilibria and their stability are studied. The main model undergoes
supercritical transcritical bifurcation at $\mathcal{R}_0^T=1$ and
$\mathcal{R}_0^H=1$ whereas the parameters $\beta^*=\beta e$ and
$\lambda^*=\lambda \sigma$ act as bifurcation parameters, respectively.
Numerical simulation claims the existence of interior equilibrium when both the
reproduction numbers are greater than unity. We explore the effect of early and
late HIV treatment on disease-induced deaths during the TB treatment course.
Mathematical analysis of our model shows that successful disease eradication
requires treatment of single disease, that is, treatment for HIV only and TB
only infected individuals with addition to co-infection treatment and in
absence of which disease eradication is extremely difficult even for
$\mathcal{R}_0<1$. When both the diseases are epidemic, the treatment for TB
only infected individuals is very effective in reducing the total infected
population and disease-induced deaths in comparison to the treatment for HIV
infected individuals while these are minimum when both the single disease
treatments are given with co-infection treatment.


  In evolutionary dynamics, well-mixed populations are almost always associated
with all-to-all interactions; mathematical models are based on complete graphs.
In most cases, these models do not predict fixation probabilities in groups of
individuals mixed by flows. We propose an analytical description in the
fast-flow limit. This approach is valid for processes with global and local
selection, and accurately predicts the suppression of selection as competition
becomes more local. It provides a modelling tool for biological or social
systems with individuals in motion.


  We investigate the effects of the relative dielectric coefficient on ionic
flows in open ion channels, using mathematical analysis of reasonably general
Poisson-Nernst-Planck type models that can include the finite sizes of ions.
The value of the relative dielectric coefficient is of course a crucial
parameter for ionic behavior in general. Using the powerful theory of
singularly perturbed problems in applied mathematics, we show that some
properties of open channels are quite insensitive to variation in the relative
dielectric coefficient, thereby explaining such effects seen unexpectedly in
simulations. The ratio between the total number of one ion species and that of
another ion species, and the ratio between the flux of one ion species and that
of another ion species do not depend significantly on the relative dielectric
coefficient.


  Determining the distribution of disease prevalence among heterogeneous
populations at the national scale is fundamental for epidemiology and public
health. Here, we use a combination of methods (spatial scan statistic,
topological data analysis, epidemic profile) to study measurable differences in
malaria intensity by regions and populations of Colombia. This study explores
three main questions: What are the regions of Colombia where malaria is
epidemic? What are the regions and populations in Colombia where malaria is
endemic? What associations exist between epidemic outbreaks between regions in
Colombia? \textit{Plasmodium falciparum} is most prevalent in the Pacific
Coast, some regions of the Amazon Basin, and some regions of the Magdalena
Basin. \textit{Plasmodium vivax} is the most prevalent parasite in Colombia,
particularly in the Northern Amazon Basin, the Caribbean, and municipalities of
Sucre, Antioquia and Cordoba. Malaria has been reported to be most common among
15-45 year old men. We find that the age-class suffering high risk of malaria
infection ranges from 20 to 30 with an acute peak at 25 years of age. Second,
this pattern was not found to be generalizable across Colombian populations,
Indigenous and Afrocolombian populations experience endemic malaria (with
household transmission). Third, clusters of epidemic malaria for
\textit{Plasmodium vivax} were detected across Southern Colombia including the
Amazon Basin and the Southern Pacific region. \textit{Plasmodium falciparum},
was is epidemic in 13 of the 1,123 municipalities (1.2\%). Some key locations
act as bridges between epidemic and endemic regions. Finally, we generate a
regional classification based on intensity and synchrony, dividing the country
into epidemic areas and bridge areas.


  The effect of phenotypic plasticity on evolution, the so-called Baldwin
effect, has been studied extensively for more than 100 years. Plasticity is
known to influence the speed of evolution towards a specific genetic
configuration, but whether it also influences what that genetic configuration
is, is still an open question. This question is investigated, in an environment
where the distribution of resources follows seasonal cycles, both analytically
and experimentally by means of an agent-based model of a foraging task.
Individuals can either specialize to foraging only one specific resource type
or generalize to foraging all resource types at a low success rate. It is found
that the introduction of learning, one instance of phenotypic plasticity,
changes what genetic configuration evolves. Specifically, the genome of
learning agents evolves a predisposition to adapt quickly to changes in the
resource distribution, under the same conditions for which non-learners would
evolve a predisposition to maximize the foraging efficiency for a specific
resource type. This paper expands the literature at the interface between
Biology and Machine Learning by identifying the Baldwin effects in
cyclically-changing environments and demonstrating that learning can change the
outcome of evolution.


  Dimensionality reduction techniques have found great success in a wide range
of fields requiring analysis of high-dimensional datasets. Time-lagged
independent components analysis (TICA), which finds independent components
(TICs) with maximal autocorrelation, is often applied to atomistic biomolecular
simulations, where the full molecular configuration can be projected onto only
a few TICs describing the slowest modes of motion. Recently, Sultan and Pande
have proposed the use of TICs as collective variables for enhanced sampling.
However, it is unclear what the best strategy for estimating the TICs of a
system is a priori. In order to evaluate the utility of TICs calculated on one
system to describe the slow dynamics of similar systems, we develop a
methodology for measuring the transferability of TICs and apply it to a wide
range of systems. We find that transferred TICs can approximate the slowest
dynamics of some systems surprisingly well, while failing to transfer between
other sets of systems, highlighting the inherent difficulties of predicting TIC
transferability. Additionally, we use two dimensional Brownian dynamics
simulations on similar potential surfaces to gain insight into the relationship
between TIC transferability and potential surface changes.


  One of the key limitations of Molecular Dynamics simulations is the
computational intractability of sampling protein conformational landscapes
associated with either large system size or long timescales. To overcome this
bottleneck, we present the REinforcement learning based Adaptive samPling
(REAP) algorithm that aims to efficiently sample conformational space by
learning the relative importance of each reaction coordinate as it samples the
landscape. To achieve this, the algorithm uses concepts from the field of
reinforcement learning, a subset of machine learning, which rewards sampling
along important degrees of freedom and disregards others that do not facilitate
exploration or exploitation. We demonstrate the effectiveness of REAP by
comparing the sampling to long continuous MD simulations and least-counts
adaptive sampling on two model landscapes (L-shaped and circular), and
realistic systems such as alanine dipeptide and Src kinase. In all four
systems, the REAP algorithm consistently demonstrates its ability to explore
conformational space faster than the other two methods when comparing the
expected values of the landscape discovered for a given amount of time. The key
advantage of REAP is on-the-fly estimation of the importance of collective
variables, which makes it particularly useful for systems with limited
structural information.


  Temporal variations in biological systems and more generally in natural
sciences are typically modelled as a set of Ordinary, Partial, or Stochastic
Differential or Difference Equations. Algorithms for learning the structure and
the parameters of a dynamical system are distinguished based on whether time is
discrete or continuous, observations are time-series or time-course, and
whether the system is deterministic or stochastic, however, there is no
approach able to handle the various types of dynamical systems simultaneously.
In this paper, we present a unified approach to infer both the structure and
the parameters of nonlinear dynamical systems of any type under the restriction
of being linear with respect to the unknown parameters. Our approach, which is
named Unified Sparse Dynamics Learning (USDL), constitutes of two steps. First,
an atemporal system of equations is derived through the application of the weak
formulation. Then, assuming a sparse representation for the dynamical system,
we show that the inference problem can be expressed as a sparse signal recovery
problem, allowing the application of an extensive body of algorithms and
theoretical results. Results on simulated data demonstrate the efficacy and
superiority of the USDL algorithm as a function of the experimental setup
(sample size, number of time measurements, number of interventions, noise
level). Additionally, USDL's accuracy significantly correlates with theoretical
metrics such as the exact recovery coefficient. On real single-cell data, the
proposed approach is able to induce high-confidence subgraphs of the signaling
pathway. USDL algorithm has been integrated in SCENERY
(\url{http://scenery.csd.uoc.gr/}); an online tool for single-cell mass
cytometry analytics.


  Assessing and managing the impact of large-scale epidemics considering only
the individual risk and severity of the disease is exceedingly difficult and
could be extremely expensive. Economic consequences, infrastructure and service
disruption, as well as the recovery speed, are just a few of the many
dimensions along which to quantify the effect of an epidemic on society's
fabric. Here, we extend the concept of resilience to characterize epidemics in
structured populations, by defining the system-wide critical functionality that
combines an individual's risk of getting the disease (disease attack rate) and
the disruption to the system's functionality (human mobility deterioration). By
studying both conceptual and data-driven models, we show that the integrated
consideration of individual risks and societal disruptions under resilience
assessment framework provides an insightful picture of how an epidemic might
impact society. In particular, containment interventions intended for a
straightforward reduction of the risk may have net negative impact on the
system by slowing down the recovery of basic societal functions. The presented
study operationalizes the resilience framework, providing a more nuanced and
comprehensive approach for optimizing containment schemes and mitigation
policies in the case of epidemic outbreaks.


  Recombinant Inbred Lines derived from divergent parental lines can display
extensive segregation distortion and long-range linkage disequilibrium (LD)
between distant loci. These genomic signatures are consistent with epistatic
selection during inbreeding. Epistatic interactions affect growth and fertility
traits or even cause complete lethality. Detecting epistasis is challenging as
multiple testing approaches are under-powered and true long-range LD is
difficult to distinguish from drift.
  Here we develop a method for reconstructing an underlying network of genomic
signatures of high-dimensional epistatic selection from multi-locus genotype
data. The network captures the conditionally dependent short- and long-range LD
structure and thus reveals "aberrant" marker-marker associations that are due
to epistatic selection rather than gametic linkage. The network estimation
relies on penalized Gaussian copula graphical models, which accounts for a
large number of markers p and a small number of individuals n.
  A multi-core implementation of our algorithm makes it feasible to estimate
the graph in high-dimensions also in the presence of significant portions of
missing data. We demonstrate the efficiency of the proposed method on simulated
datasets as well as on genotyping data in A.thaliana and maize. In addition, we
implemented the method in the R package netgwas which is freely available at
https://CRAN.R-project.org/package=netgwas.


  Linkage maps are used to identify the location of genes responsible for
traits and diseases. New sequencing techniques have created opportunities to
substantially increase the density of genetic markers. Such revolutionary
advances in technology have given rise to new challenges, such as creating
high-density linkage maps. Current multiple testing approaches based on
pairwise recombination fractions are underpowered in the high-dimensional
setting and do not extend easily to polyploid species. We propose to construct
linkage maps using graphical models either via a sparse Gaussian copula or a
nonparanormal skeptic approach. Linkage groups (LGs), typically chromosomes,
and the order of markers in each LG are determined by inferring the conditional
independence relationships among large numbers of markers in the genome.
Through simulations, we illustrate the utility of our map construction method
and compare its performance with other available methods, both when the data
are clean and contain no missing observations and when data contain genotyping
errors and are incomplete. We apply the proposed method to two genotype
datasets: barley and potato from diploid and polypoid populations,
respectively. Our comprehensive map construction method makes full use of the
dosage SNP data to reconstruct linkage map for any bi-parental diploid and
polyploid species. We have implemented the method in the R package netgwas.


  Wildlife corridors are components of landscapes, which facilitate the
movement of organisms and processes between intact habitat areas, and thus
provide connectivity between the habitats within the landscapes. Corridors are
thus regions within a given landscape that connect fragmented habitat patches
within the landscape. The major concern of designing corridors as a
conservation strategy is primarily to counter, and to the extent possible,
mitigate the effects of habitat fragmentation and loss on the biodiversity of
the landscape, as well as support continuance of land use for essential local
and global economic activities in the region of reference. In this paper, we
use game theory, graph theory, membership functions and chain code algorithm to
model and design a set of wildlife corridors with tiger (Panthera tigris
tigris) as the focal species. We identify the parameters which would affect the
tiger population in a landscape complex and using the presence of these
identified parameters construct a graph using the habitat patches supporting
tiger presence in the landscape complex as vertices and the possible paths
between them as edges. The passage of tigers through the possible paths have
been modelled as an Assurance game, with tigers as an individual player. The
game is played recursively as the tiger passes through each grid considered for
the model. The iteration causes the tiger to choose the most suitable path
signifying the emergence of adaptability. As a formal explanation of the game,
we model this interaction of tiger with the parameters as deterministic finite
automata, whose transition function is obtained by the game payoff.


  We study the dynamics of a predator-prey system in a random environment. The
dynamics evolves according to a deterministic Lotka-Volterra system for an
exponential random time after which it switches to a different deterministic
Lotka-Volterra system. This switching procedure is then repeated. The resulting
process is a Piecewise Deterministic Markov Process (PDMP). In the case when
the equilibrium points of the two deterministic Lotka--Volterra systems
coincide we show that almost surely the trajectory does not converge to the
common deterministic equilibrium. Instead, with probability one, the densities
of the prey and the predator oscillate between $0$ and $\infty$. This proves a
conjecture of Takeuchi et al (J. Math. Anal. Appl 2006).
  The proof of the conjecture is a corollary of a result we prove about linear
switched systems. Assume $(Y_t, I_t)$ is a PDMP that evolves according to
$\frac{dY_t}{dt}=A_{I_t} Y_t$ where $A_0,A_1$ are $2\times2$ matrices and $I_t$
is a Markov chain on $\{0,1\}$ with transition rates $k_0,k_1>0$. If the
matrices $A_0$ and $A_1$ are not proportional and have purely imaginary
eigenvalues, then there exists $\lambda >0$ such that $\lim_{t \to \infty}
\frac{\log \| Y_t \|}{t} = \lambda.$


  We consider the harvesting of a population in a stochastic environment whose
dynamics in the absence of harvesting is described by a one dimensional
diffusion. Using ergodic optimal control, we find the optimal harvesting
strategy which maximizes the asymptotic yield of harvested individuals. To our
knowledge, ergodic optimal control has not been used before to study harvesting
strategies. However, it is a natural framework because the optimal harvesting
strategy will never be such that the population is harvested to extinction --
instead the harvested population converges to a unique invariant probability
measure.
  When the yield function is the identity, we show that the optimal strategy
has a bang-bang property: there exists a threshold $x^*>0$ such that whenever
the population is under the threshold the harvesting rate must be zero, whereas
when the population is above the threshold the harvesting rate must be at the
upper limit. We provide upper and lower bounds on the maximal asymptotic yield,
and explore via numerical simulations how the harvesting threshold and the
maximal asymptotic yield change with the growth rate, maximal harvesting rate,
or the competition rate.
  We also show that, if the yield function is $C^2$ and strictly concave, then
the optimal harvesting strategy is continuous, whereas when the yield function
is convex the optimal strategy is of bang-bang type. This shows that one cannot
always expect bang-bang type optimal controls.


  Graphical models are a powerful tool in modelling and analysing complex
biological associations in high-dimensional data. The R-package netgwas
implements the recent methodological development on copula graphical models to
(i) construct linkage maps, (ii) infer linkage disequilibrium networks from
genotype data, and (iii) detect high-dimensional genotype-phenotype networks.
The netgwas learns the structure of networks from ordinal data and mixed
ordinal-and-continuous data. Here, we apply the functionality in netgwas to
various multivariate example datasets taken from the literature to demonstrate
the kind of insight that can be obtained from the package. We show that our
package offers a more realistic association analysis than the classical
approaches, as it discriminates between direct and induced correlations by
adjusting for the effect of all other variables while performing pairwise
associations. This feature controls for spurious interactions between variables
that can arise from conventional approaches in a biological sense. The netgwas
package uses a parallelization strategy on multi-core processors to speed-up
computations. The netgwas package is freely available at
https://cran.r-project.org/web/packages/netgwas


  We consider the task of detecting regulatory elements in the human genome
directly from raw DNA. Past work has focused on small snippets of DNA, making
it difficult to model long-distance dependencies that arise from DNA's
3-dimensional conformation. In order to study long-distance dependencies, we
develop and release a novel dataset for a larger-context modeling task. Using
this new data set we model long-distance interactions using dilated
convolutional neural networks, and compare them to standard convolutions and
recurrent neural networks. We show that dilated convolutions are effective at
modeling the locations of regulatory markers in the human genome, such as
transcription factor binding sites, histone modifications, and DNAse
hypersensitivity sites.


  Deterministic evolutionary game dynamics can lead to stable coexistences of
different types. Stochasticity, however, drives the loss of such coexistences.
This extinction is usually accompanied by population size fluctuations. We
investigate the most probable extinction trajectory under such fluctuations by
mapping a stochastic evolutionary model to a problem of classical mechanics
using the Wentzel-Kramers-Brillouin (WKB) approximation. Our results show that
more abundant types in a coexistence can be more likely to go extinct first
well agreed with previous results, and also the distance between the
coexistence and extinction point is not a good predictor of extinction.
Instead, the WKB method correctly predicts the type going extinct first.


  Blockage of some ion channels and in particular, the hERG cardiac potassium
channel delays cardiac repolarization and can induce arrhythmia. In some cases
it leads to a potentially life-threatening arrhythmia known as Torsade de
Pointes (TdP). Therefore recognizing drugs with TdP risk is essential.
Candidate drugs that are determined not to cause cardiac ion channel blockage
are more likely to pass successfully through clinical phases II and III trials
(and preclinical work) and not be withdrawn even later from the marketplace due
to cardiotoxic effects. The objective of the present study is to develop an SAR
model that can be used as an early screen for torsadogenic (causing TdP
arrhythmias) potential in drug candidates. The method is performed using
descriptors comprised of atomic NMR chemical shifts and corresponding
interatomic distances which are combined into a 3D abstract space matrix. The
method is called 3D-SDAR (3 dimensional spectral data-activity relationship)
and can be interrogated to identify molecular features responsible for the
activity, which can in turn yield simplified hERG toxicophores. A dataset of 55
hERG potassium channel inhibitors collected from Kramer et al. consisting of 32
drugs with TdP risk and 23 with no TdP risk was used for training the 3D-SDAR
model.An ANN model with multilayer perceptron was used to define collinearities
among the independent 3D-SDAR features. A composite model from 200 random
iterations with 25% of the molecules in each case yielded the following figures
of merit: training, 99.2 %; internal test sets, 66.7%; external (blind
validation) test set, 68.4%. In the external test set, 70.3% of positive TdP
drugs were correctly predicted. Moreover, toxicophores were generated from TdP
drugs. A 3D-SDAR was successfully used to build a predictive model for
drug-induced torsadogenic and non-torsadogenic drugs.


  We consider non-demographic noise in the form of uncertainty in the reaction
step size, and reveal a dramatic effect this noise may have on the stability of
self-regulating populations. Employing the reaction scheme mA->kA, but
allowing, e.g., the product number k to be a-priori unknown and sampled from a
given distribution, we show that such non-demographic noise can greatly reduce
the population's extinction risk compared to the fixed k case. Our analysis is
tested against numerical simulations, and by using empirical data of different
species, we argue that certain distributions may be more evolutionary
beneficial than others.


  A principal component analysis of a multiple sequence alignement of
hemagglutinin sequences of subtype H1 has been performed, the sequences being
encoded using the amino-acid property that maximizes the weight of the major
component. In the case of this alignment, it happens to be a well-known
hydrophobicity scale. Interestingly, sequences coming from human have large
positive amplitudes along the major component before 2009, and large negative
ones afterwards. This means that the 2009 pandemic was associated to a major
change in the hydrophobicity pattern of hemagglutinin.
  The present analysis also highlights the high variability of viral sequences
coming from swine. At a more general level, the method proposed in this paper
allows to describe a sequence coming from an alignment with a set of numbers,
the original point being that the choice of the corresponding property is
driven by the data. This approach should allow the application of numerous
methods to the study of large multiple sequence alignments.


  To estimate the time, many organisms, ranging from cyanobacteria to animals,
employ a circadian clock which is based on a limit-cycle oscillator that can
tick autonomously with a nearly 24h period. Yet, a limit-cycle oscillator is
not essential for knowing the time, as exemplified by bacteria that possess an
'hourglass': a system that when forced by an oscillatory light input exhibits
robust oscillations from which the organism can infer the time, but that in the
absence of driving relaxes to a stable fixed point. Here, using models of the
Kai system of cyanobacteria, we compare a limit- cycle oscillator with two
hourglass models, one that without driving relaxes exponentially and one that
does so in an oscillatory fashion. In the limit of low input-noise, all three
systems are equally informative on time, yet in the regime of high input-noise
the limit-cycle oscillator is far superior. The same behavior is found in the
Stuart-Landau model, indicating that our result is universal.


  Phenotypes of individuals in a population of organisms are not fixed.
Phenotypic fluctuations, which describe temporal variation of the phenotype of
an individual or individual-to-individual variation across a population, are
present in populations from microbes to higher animals. Phenotypic fluctuations
can provide a basis for adaptation and be the target of selection. Here we
present a theoretical and experimental investigation of the fate of phenotypic
fluctuations in directed evolution experiments where phenotypes are subject to
constraints. We show that selecting bacterial populations for fast migration
through a porous environment drives a reduction in cell-to-cell variation
across the population. Using sequencing and genetic engineering we study the
genetic basis for this reduction in phenotypic fluctuations. We study the
generality of this reduction by developing a simple, abstracted, numerical
simulation model of the evolution of phenotypic fluctuations subject to
constraints. Using this model we find that strong and weak selection generally
lead respectively to increasing or decreasing cell-to-cell variation as a
result of a bound on the selected phenotype under a wide range of parameters.
However, other behaviors are also possible, and we describe the outcome of
selection simulations for different model parameters and suggest future
experiments. We analyze the mechanism of the observed reduction of phenotypic
fluctuations in our experimental system, discuss the relevance of our abstract
model to the experiment and explore its broader implications for evolution.


  At time 0, start a time-continuous binary branching process, where particles
give birth to a single particle independently (at a possibly time-dependent
rate) and die independently (at a possibly time-dependent and age-dependent
rate). A particular case is the classical birth--death process. Stop this
process at time $T>0$. It is known that the tree spanned by the $N$ tips alive
at time $T$ of the tree thus obtained (called reduced tree or coalescent tree)
is a coalescent point process (CPP), which basically means that the depths of
interior nodes are iid. Now select each of the $N$ tips independently with
probability $y$ (Bernoulli sample). It is known that the tree generated by the
selected tips, which we will call Bernoulli sampled CPP, is again a CPP. Now
instead, select exactly $k$ tips uniformly at random among the $N$ tips
($k$-sample). We show that the tree generated by the selected tips is a mixture
of Bernoulli sampled CPPs with the same parent CPP, over some explicit
distribution of the sampling probability $y$. An immediate consequence is that
the genealogy of a $k$-sample can be obtained by the realization of $k$ random
variables, first the random sampling probability $Y$ and then the $k-1$ node
depths which are iid conditional on $Y=y$.


  Membrane protein transporters alternate their substrate-binding sites between
the extracellular and cytosolic side of the membrane according to the
alternating access mechanism. Inspired by this intriguing mechanism devised by
nature, we study particle transport through a channel coupled with an energy
well that oscillates its position between the two entrances of the channel. We
optimize particle transport across the channel by adjusting the oscillation
frequency. At the optimal oscillation frequency, the translocation rate through
the channel is a hundred times higher with respect to free diffusion across the
channel. Our findings reveal the effect of time dependent potentials on
particle transport across a channel and will be relevant for membrane transport
and microfluidics application.


  We study the electro-diffusion properties of a domain containing a
cusp-shaped structure in three dimensions when one ionic specie is dominant.
The mathematical problem consists in solving the steady-state
Poisson-Nernst-Planck (PNP) equation with an integral constraint for the number
of charges. A non-homogeneous Neumann boundary condition is imposed on the
boundary. We construct an asymptotic approximation for certain singular limits
that agree with numerical simulations. Finally, we analyse the consequences of
non-homogeneous surface charge density. We conclude that the geometry of
cusp-shaped domains influences the voltage profile, specifically inside the
cusp structure. The main results are summarized in the form of new
three-dimensional electrostatic laws for non-electroneutral electrolytes. We
discuss applications to dendritic spines in neuroscience.


  Natural protein molecules are exceptional polymers. Encoded in apparently
random strings of amino-acids, these objects perform clear physical tasks that
are rare to find by simple chance. Accurate folding, specific binding, powerful
catalysis, are examples of basic chemical activities that the great majority of
polypeptides do not display, and are thought to be the outcome of the natural
history of proteins. Function, a concept genuine to Biology, is at the core of
evolution and often conflicts with the physical constraints. Locating the
frustration between discrepant goals in a recurrent system leads to fundamental
insights about the chances and necessities that shape the encoding of
biological information.


  Principal component analysis is a widely-used method for the dimensionality
reduction of a given data set in a high-dimensional Euclidean space. Here we
define and analyze two analogues of principal component analysis in the setting
of tropical geometry. In one approach, we study the Stiefel tropical linear
space of fixed dimension closest to the data points in the tropical projective
torus; in the other approach, we consider the tropical polytope with a fixed
number of vertices closest to the data points. We then give approximative
algorithms for both approaches and apply them to phylogenetics, testing the
methods on simulated phylogenetic data and on an empirical dataset of
Apicomplexa genomes.


  How producers of public goods persist in microbial communities is a major
question in evolutionary biology. Cooperation is evolutionarily unstable, since
cheating strains can reproduce quicker and take over. Spatial structure has
been shown to be a robust mechanism for the evolution of cooperation. Here we
study how spatial assortment might emerge from native dynamics and show that
fluid flow shear promotes cooperative behavior. Social structures arise
naturally from our advection-diffusion-reaction model as self-reproducing
Turing patterns. We computationally study the effects of fluid advection on
these patterns as a mechanism to enable or enhance social behavior. Our central
finding is that flow shear enables and promotes social behavior in microbes by
increasing the group fragmentation rate and thereby limiting the spread of
cheating strains. Regions of the flow domain with higher shear admit high
cooperativity and large population density, whereas low shear regions are
devoid of life due to opportunistic mutations.


  Proteins are the main workhorses of biological functions in a cell, a tissue,
or an organism. Identification and quantification of proteins in a given
sample, e.g. a cell type under normal/disease conditions, are fundamental tasks
for the understanding of human health and disease. In this paper, we present
DeepNovo, a deep learning-based tool to address the problem of protein
identification from tandem mass spectrometry data. The idea was first proposed
in the context of de novo peptide sequencing [1] in which convolutional neural
networks and recurrent neural networks were applied to predict the amino acid
sequence of a peptide from its spectrum, a similar task to generating a caption
from an image. We further develop DeepNovo to perform sequence database search,
the main technique for peptide identification that greatly benefits from
numerous existing protein databases. We combine two modules de novo sequencing
and database search into a single deep learning framework for peptide
identification, and integrate de Bruijn graph assembly technique to offer a
complete solution to reconstruct protein sequences from tandem mass
spectrometry data. This paper describes a comprehensive protocol of DeepNovo
for protein identification, including training neural network models, dynamic
programming search, database querying, estimation of false discovery rate, and
de Bruijn graph assembly. Training and testing data, model implementations, and
comprehensive tutorials in form of IPython notebooks are available in our
GitHub repository (https://github.com/nh2tran/DeepNovo).


  With the maturation of metabolomics science and proliferation of biobanks,
clinical metabolic profiling is an increasingly opportunistic frontier for
advancing translational clinical research. Automated Machine Learning (AutoML)
approaches provide exciting opportunity to guide feature selection in agnostic
metabolic profiling endeavors, where potentially thousands of independent data
points must be evaluated. In previous research, AutoML using high-dimensional
data of varying types has been demonstrably robust, outperforming traditional
approaches. However, considerations for application in clinical metabolic
profiling remain to be evaluated. Particularly, regarding the robustness of
AutoML to identify and adjust for common clinical confounders. In this study,
we present a focused case study regarding AutoML considerations for using the
Tree-Based Optimization Tool (TPOT) in metabolic profiling of exposure to
metformin in a biobank cohort. First, we propose a tandem rank-accuracy measure
to guide agnostic feature selection and corresponding threshold determination
in clinical metabolic profiling endeavors. Second, while AutoML, using default
parameters, demonstrated potential to lack sensitivity to low-effect
confounding clinical covariates, we demonstrated residual training and
adjustment of metabolite features as an easily applicable approach to ensure
AutoML adjustment for potential confounding characteristics. Finally, we
present increased homocysteine with long-term exposure to metformin as a
potentially novel, non-replicated metabolite association suggested by TPOT; an
association not identified in parallel clinical metabolic profiling endeavors.
While considerations are recommended, including adjustment approaches for
clinical confounders, AutoML presents an exciting tool to enhance clinical
metabolic profiling and advance translational research endeavors.


  Recent breakthroughs in cancer research have come via the up-and-coming field
of pathway analysis. By applying statistical methods to prior known gene and
protein regulatory information, pathway analysis provides a meaningful way to
interpret genomic data. While many gene/protein regulatory relationships have
been studied, never before has such a significant amount data been made
available in organized forms of gene/protein regulatory networks and pathways.
However, pathway analysis research is still in its infancy, especially when
applying it to solve practical problems.
  In this paper we propose a new method of studying biological pathways, one
that cross analyzes mutation information, transcriptome and proteomics data.
Using this outcome, we identify routes of aberrant pathways potentially
responsible for the etiology of disease. Each pathway route is encoded as a
bayesian network which is initialized with a sequence of conditional
probabilities specifically designed to encode directionality of regulatory
relationships encoded in the pathways. Far more complex interactions, such as
phosphorylation and methylation, among others, in the pathways can be modeled
using this approach. The effectiveness of our model is demonstrated through its
ability to distinguish real pathways from decoys on TCGA mRNA-seq, mutation,
Copy Number Variation and phosphorylation data for both Breast cancer and
Ovarian cancer study. The majority of pathways distinguished can be confirmed
by biological literature. Moreover, the proportion of correctly indentified
pathways is \% higher than previous work where only mRNA-seq mutation data is
incorporated for breast cancer patients. Consequently, such an in-depth pathway
analysis incorporating more diverse data can give rise to the accuracy of
perturbed pathway detection.


  Hepatitis C virus (HCV) infection is endemic in people who inject drugs
(PWID), with prevalence estimates above 60 percent for PWID in the United
States. Previous modeling studies suggest that direct acting antiviral (DAA)
treatment can lower overall prevalence in this population, but treatment is
often delayed until the onset of advanced liver disease (fibrosis stage 3 or
later) due to cost. Lower cost interventions featuring syringe access (SA) and
medically assisted treatment (MAT) for addiction are known to be less costly,
but have shown mixed results in lowering HCV rates below current levels. Little
is known about the potential synergistic effects of combining DAA and MAT
treatment, and large-scale tests of combined interventions are rare. While
simulation experiments can reveal likely long-term effects, most prior
simulations have been performed on closed populations of model agents--a
scenario quite different from the open, mobile populations known to most health
agencies. This paper uses data from the Centers for Disease Control's National
HIV Behavioral Surveillance project, IDU round 3, collected in New York City in
2012 by the New York City Department of Health and Mental Hygiene to
parameterize simulations of open populations. Our results show that, in an open
population, SA/MAT by itself has only small effects on HCV prevalence, while
DAA treatment by itself can significantly lower both HCV and HCV-related
advanced liver disease prevalence. More importantly, the simulation experiments
suggest that cost effective synergistic combinations of the two strategies can
dramatically reduce HCV incidence. We conclude that adopting SA/MAT
implementations alongside DAA interventions can play a critical role in
reducing the long-term consequences of ongoing infection.


  How do we determine the mutational effects in exome sequencing data with
little or no statistical evidence? Can protein structural information fill in
the gap of not having enough statistical evidence? In this work, we answer the
two questions with the goal towards determining pathogenic effects of rare
variants in rare disease. We take the approach of determining the importance of
point mutation loci focusing on protein structure features. The proposed
structure-based features contain information about geometric, physicochemical,
and functional information of mutation loci and those of structural neighbors
of the loci. The performance of the structure-based features trained on 80\% of
HumDiv and tested on 20\% of HumDiv and on ClinVar datasets showed high levels
of discernibility in the mutation's pathogenic or benign effects: F score of
0.71 and 0.68 respectively using multi-layer perceptron. Combining structure-
and sequence-based feature further improve the accuracy: F score of 0.86
(HumDiv) and 0.75 (ClinVar). Also, careful examination of the rare variants in
rare diseases cases showed that structure-based features are important in
discerning importance of variant loci.


  Molecular motors play pivotal roles in organizing the interior of cells. A
motor efficient in cargo transport would move along cytoskeletal filaments with
a high speed and a minimal error in transport distance (or time) while
consuming a minimal amount of energy. The travel distance of the motor and its
variance are, however, physically constrained by the free energy being
consumed. A recently formulated \emph{thermodynamic uncertainty relation}
offers a theoretical framework for the energy-accuracy trade-off relation
ubiquitous in biological processes. According to the relation, a measure
$\mathcal{Q}$, the product between the heat dissipated from a motor and the
squared relative error in the displacement, has a minimal theoretical bound
($\mathcal{Q} \geq 2 k_B T$), which is approached when the time trajectory of
the motor is maximally regular for a given amount of free energy input. Here,
we use $\mathcal{Q}$ to quantify the transport efficiency of biological motors.
Analyses on the motility data from several types of molecular motors reveal
that $\mathcal{Q}$ is a complex function of ATP concentration and load ($f$).
For kinesin-1, $\mathcal{Q}$ approaches the theoretical bound at $f\approx 4$
pN and over a broad range of ATP concentration (1 $\mu$M - 10 mM), and is
locally minimized at [ATP] $\approx$ 200 $\mu$M. In stark contrast, this local
minimum vanishes for a mutant that has a longer neck-linker, and the value of
$\mathcal{Q}$ is significantly greater, which underscores the importance of
molecular structure. Transport efficiencies of the biological motors studied
here are semi-optimized under the cellular condition ([ATP] $\approx 1$ mM,
$f=0-1$ pN). Our study indicates that among many possible directions of
optimization, cytoskeletal motors are designed to operate at a high speed with
a minimal error while leveraging their energy resources.