  This report explores the application of machine learning techniques on short
timeseries gene expression data. Although standard machine learning algorithms
work well on longer time-series', they often fail to find meaningful insights
from fewer timepoints. In this report, we explore model-based clustering
techniques. We combine popular unsupervised learning techniques like K-Means,
Gaussian Mixture Models, Bayesian Networks, Hidden Markov Models with the
well-known Expectation Maximization algorithm. K-Means and Gaussian Mixture
Models are fairly standard, while Hidden Markov Model and Bayesian Networks
clustering are more novel ideas that suit time-series gene expression data.


  The coronavirus disease (COVID-19) pandemic has changed our lives and still
poses a challenge to science. Numerous studies have contributed to a better
understanding of the pandemic. In particular, inhalation of aerosolised
pathogens has been identified as essential for transmission. This information
is crucial to slow the spread, but the individual likelihood of becoming
infected in everyday situations remains uncertain. Mathematical models help
estimate such risks. In this study, we propose how to model airborne
transmission of SARS-CoV-2 at a local scale. In this regard, we combine
microscopic crowd simulation with a new model for disease transmission.
Inspired by compartmental models, we describe agents' health status as
susceptible, exposed, infectious or recovered. Infectious agents exhale
pathogens bound to persistent aerosols, whereas susceptible agents absorb
pathogens when moving through an aerosol cloud left by the infectious agent.
The transmission depends on the pathogen load of the aerosol cloud, which
changes over time. We propose a 'high risk' benchmark scenario to distinguish
critical from non-critical situations. Simulating indoor situations show that
the new model is suitable to evaluate the risk of exposure qualitatively and,
thus, enables scientists or even decision-makers to better assess the spread of
COVID-19 and similar diseases.


  In a stochastic reaction network setting we consider the problem of tracking
the fate of individual molecules. We show that using the classical large volume
limit results, we may approximate the dynamics of a single tracked molecule in
a simple and computationally efficient way. We give examples on how this
approach may be used to obtain various characteristics of single-molecule
dynamics (for instance, the distribution of the number of infections in a
single individual in the course of an epidemic or the activity time of a single
enzyme molecule). Moreover, we show how to approximate the overall dynamics of
species of interest in the full system with a collection of independent
single-molecule trajectories, and give explicit bounds for the approximation
error in terms of the reaction rates. This approximation, which is well defined
for all times, leads to an efficient and fully parallelizable simulation
technique for which we provide some numerical examples.


  Non-equilibrium thermodynamics has long been an area of substantial interest
to ecologists because most fundamental biological processes, such as protein
synthesis and respiration, are inherently energy-consuming. Microbial
communities are a natural system to decipher this mechanistic basis because
their interactions in the form of substrate consumption, metabolite production,
and cross-feeding can be described explicitly in thermodynamic terms. Previous
work has considered how thermodynamic constraints impact competition between
pairs of species, but restrained from analysing how this manifests in complex
dynamical systems. To address this gap, we develop a thermodynamic microbial
community model with fully reversible reaction kinetics, which allows direct
consideration of free-energy dissipation. This also allows species to interact
via products rather than just substrates, increasing the dynamical complexity,
and allowing a more nuanced classification of interaction types to emerge.
Using this model, we find that community diversity increases with substrate
lability, because greater free-energy availability allows for faster generation
of niches. Thus, more niches are generated in the time frame of community
establishment, leading to higher final species diversity. We also find that
allowing species to make use of near-to-equilibrium reactions increases
diversity in a low free-energy regime. In such a regime, two new thermodynamic
interaction types that we identify here reach comparable strengths to the
conventional (competition and facilitation) types, emphasising the key role
that thermodynamics plays in community dynamics. Our results suggest that
accounting for realistic thermodynamic constraints is vital for understanding
the dynamics of real-world microbial communities.


  Connecting the different scales of epidemic dynamics, from individuals to
communities to nations, remains one of the main challenges of disease modeling.
Here, we revisit one of the largest public health efforts deployed against a
localized epidemic: the 2014-2016 Ebola Virus Disease (EVD) epidemic in Sierra
Leone. We leverage the data collected by the surveillance and contact tracing
protocols of the Sierra Leone Ministry of Health and Sanitation, the US Centers
for Disease Control and Prevention, and other responding partners to validate a
network epidemiology framework connecting the population (incidence), community
(local forecasts), and individual (secondary infections) scales of disease
transmission. In doing so, we gain a better understanding of what brought the
EVD epidemic to an end: Reduction of introduction in new clusters (primary
cases), and not reduction in local transmission patterns (secondary
infections). We also find that the first 90 days of the epidemic contained
enough information to produce probabilistic forecasts of EVD cases; forecasts
which we show are confirmed independently by both disease surveillance and
contact tracing. Altogether, using data available two months before the start
of the international support to the local response, network epidemiology could
have inferred heterogeneity in local transmissions, the risk for superspreading
events, and probabilistic forecasts of eventual cases per community. We expect
that our framework will help connect large data collection efforts with
individual behavior, and help reduce uncertainty during health emergencies and
emerging epidemics.


  Environmental fluctuations can create population-depleted areas and even
extinct areas for the population. This effect is more severe in the presence of
the Allee effect (decreasing growth rate at low population densities).
Dispersal inside the habitat provides a rescue effect on population-depleted
areas, enhancing the population resilience to environmental fluctuations.
Habitat reduction decreases the effectiveness of the dispersal rescue
mechanism. We report here how the population resilience to environmental
fluctuations decreases when the dispersal length or the habitat size are
reduced. The resilience reduction is characterized by a decrease of the
extinction threshold for environmental fluctuations. The extinction threshold
is shown to scale with the ratio between the dispersal length and the scale of
environmental synchrony, i.e., it is the dispersal connection between
non-environmentally-correlated regions that provides resilience to
environmental fluctuations. Habitat reduction also decreases the resilience to
environmental fluctuations, when the habitat size is similar to or smaller than
the characteristic dispersal distances. The power laws of these scaling
behaviors are characterized here. Alternative scaling functions with spatial
scales of population synchrony are found to fit the simulations worse. These
results support the dispersal length as the critical scale for extinction
induced by habitat reduction.


  We consider the effect of non-constant parameters on the human-forest
interaction logistic model coupled with human technological growth introduced
in "Deforestation and world population sustainability: a quantitative
analysis"[1]. In recent years in fact, a decrease in human population growth
rate has emerged which can be measured to about 1.7% drop per year since 1960
value which coincides with latest UN projections for next decades up to year
2100 [2]. We therefore consider here the effect of decreasing human population
growth-rate on the aforementioned model and we evaluate its effect on the
probability of survival of human civilisation without going through a
catastrophic collapse in population. We find that for realistic values of the
human population carrying capacity of the earth (measured by parameter beta)
this decrease would not affect previous results leading to a low probability of
avoiding a catastrophic collapse. For larger more optimistic values of beta
instead, a decrease in growth-rate would tilt the probability in favour of a
positive outcome, i.e. from 10-20% up to even 95% likelihood of avoiding
collapse.


  This paper presents a simple mathematical model for the growth of the
Christian population in the Roman Empire during the first to fourth centuries.
The model has a subexponential growth rate of order $e^{o(t)}$, where $o$
denotes the "little-o" asymptotic bound, but still superpolynomial, and it fits
available Christian population estimates with good accuracy.


  We introduce a modified SIR model with memory for the dynamics of epidemic
spreading in a constant population of individuals. Each individual is in one of
the states susceptible (${\bf S}$), infected (${\bf I}$) or recovered (${\bf
R}$). In the state ${\bf R}$ an individual is assumed to stay immune within a
finite time interval. In the first part, we introduce a random life time or
duration of immunity which is drawn from a certain probability density
function. Once the time of immunity is elapsed an individual makes an
instantaneous transition to the susceptible state. By introducing a random
duration of immunity a memory effect is introduced into the process which
crucially determines the epidemic dynamics. In the second part, we investigate
the influence of the memory effect on the space-time dynamics of the epidemic
spreading by implementing this approach into computer simulations and employ a
multiple random walker's model. If a susceptible walker meets an infectious one
on the same site, then the susceptible one gets infected with a certain
probability. The computer experiments allow us to identify relevant parameters
for spread or extinction of an epidemic. In both parts, the finite duration of
immunity causes persistent oscillations in the number of infected individuals
with ongoing epidemic activity preventing the system from relaxation to a
steady state solution. Such oscillatory behavior is supported by real-life
observations and cannot be captured by standard SIR models.


  The reduction of computational costs for marine ecosystem models is important
for the investigation and detection of the relevant biogeochemical processes
because such models are computationally expensive. In order to lower these
computational costs by means of larger time steps we investigated the accuracy
of steady annual cycles (i.e., an annual periodic solution) calculated with
different time steps. We compared the accuracy for a hierarchy of
biogeochemical models showing an increasing complexity and computed the steady
annual cycles with offline simulations that are based on the transport matrix
approach. For each of these biogeochemical models, we obtained practically the
same solution even though larger time steps. This indicates that larger time
steps shortened the runtime with an acceptable loss of accuracy.


  Principal Component Analysis (PCA) and other multi-variate models are often
used in the analysis of "omics" data. These models contain much information
which is currently neither easily accessible nor interpretable. Here we present
an algorithmic method which has been developed to integrate this information
with existing databases of background knowledge, stored in the form of known
sets (for instance genesets or pathways). To make this accessible we have
produced a Graphical User Interface (GUI) in Matlab which allows the overlay of
known set information onto the loadings plot and thus improves the
interpretability of the multi-variate model. For each known set the optimal
convex hull, covering a subset of elements from the known set, is found through
a search algorithm and displayed. In this paper we discuss two main topics; the
details of the search algorithm for the optimal convex hull for this problem
and the GUI interface which is freely available for download for academic use.


  Several physical mechanisms have been proposed to explain allostery in
proteins. They differ by the number of internal states that they assume a
protein to occupy, leaving open the question of what controls the emergence of
these distinct physical forms of allostery. Here, we analyze a simplified model
of protein allostery under a range of physical and evolutionary constraints. We
find that two archetypal mechanisms can emerge through evolution: a
single-state mechanism where ligand binding induces a displacement along a soft
normal mode or a multi-state mechanism where ligand binding induces a switch
across an energy barrier to a different stable state. Importantly, whenever the
two mechanisms are possible, the multi-state mechanism confers a stronger
allosteric effect and thus a selective advantage. This work defines the
essential constraints on single or multi-state allostery, and sets the stage
for a physical theory of its evolutionary origins.


  As of December 2020, the COVID-19 pandemic has infected over 75 million
people, making it the deadliest pandemic in modern history. This study develops
a novel compartmental epidemiological model specific to the SARS-CoV-2 virus
and analyzes the effect of common preventative measures such as testing,
quarantine, social distancing, and vaccination. By accounting for the most
prevalent interventions that have been enacted to minimize the spread of the
virus, the model establishes a paramount foundation for future mathematical
modeling of COVID-19 and other modern pandemics. Specifically, the model
expands on the classic SIR model and introduces separate compartments for
individuals who are in the incubation period, asymptomatic, tested-positive,
quarantined, vaccinated, or deceased. It also accounts for variable infection,
testing, and death rates. I first analyze the outbreak in Santa Clara County,
California, and later generalize the findings. The results show that, although
all preventative measures reduce the spread of COVID-19, quarantine and social
distancing mandates reduce the infection rate and subsequently are the most
effective policies, followed by vaccine distribution and, finally, public
testing. Thus, governments should concentrate resources on enforcing quarantine
and social distancing policies. In addition, I find mathematical proof that the
relatively high asymptomatic rate and long incubation period are driving
factors of COVID-19's rapid spread.


  Accurately predicting the likelihood of interaction between two objects
(compound-protein sequence, user-item, author-paper, etc.) is a fundamental
problem in Computer Science. Current deep-learning models rely on learning
accurate representations of the interacting objects. Importantly, relationships
between the interacting objects, or features of the interaction, offer an
opportunity to partition the data to create multi-views of the interacting
objects. The resulting congruent and non-congruent views can then be exploited
via contrastive learning techniques to learn enhanced representations of the
objects.


  RaTG13 is the closest related coronavirus genome phylogenetically to
SARS-CoV-2, consequently understanding its provenance is of key importance to
understanding the origin of the COVID-19 pandemic. The RaTG13 NGS dataset is
attributed to a fecal swab from the intermediate horseshoe bat Rhinolophus
affinis. However, sequence analysis reveals that this is unlikely. Metagenomic
analysis using Metaxa2 shows that only 10.3 % of small subunit (SSU) rRNA
sequences in the dataset are bacterial, inconsistent with a fecal sample, which
are typically dominated by bacterial sequences. In addition, the bacterial taxa
present in the sample are inconsistent with fecal material. Assembly of
mitochondrial SSU rRNA sequences in the dataset produces a contig 98.7 %
identical to R.affinis mitochondrial SSU rRNA, indicating that the sample was
generated from this or a closely related species. 87.5 % of the NGS reads map
to the Rhinolophus ferrumequinum genome, the closest bat genome to R.affinis
available. In the annotated genome assembly, 62.2 % of mapped reads map to
protein coding genes. These results clearly demonstrate that the dataset
represents a Rhinolophus sp. transcriptome, and not a fecal swab sample.
Overall, the data show that the RaTG13 dataset was generated by the Wuhan
Institute of Virology (WIV) from a transcriptome derived from Rhinolophus sp.
tissue or cell line, indicating that RaTG13 was in live culture. This raises
the question of whether the WIV was culturing additional unreported
coronaviruses closely related to SARS-CoV-2 prior to the pandemic. The
implications for the origin of the COVID-19 pandemic are discussed.


  Machine learning shows great potential in virtual screening for drug
discovery. Current efforts on accelerating docking-based virtual screening do
not consider using existing data of other previously developed targets. To make
use of the knowledge of the other targets and take advantage of the existing
data, in this work, we apply multi-task learning to the problem of
docking-based virtual screening. With two large docking datasets, the results
of extensive experiments show that multi-task learning can achieve better
performances on docking score prediction. By learning knowledge across multiple
targets, the model trained by multi-task learning shows a better ability to
adapt to a new target. Additional empirical study shows that other problems in
drug discovery, such as the experimental drug-target affinity prediction, may
also benefit from multi-task learning. Our results demonstrate that multi-task
learning is a promising machine learning approach for docking-based virtual
screening and accelerating the process of drug discovery.


  The reconstruction of microbial genomes from large metagenomic datasets is a
critical procedure for finding uncultivated microbial populations and defining
their microbial functional roles. To achieve that, we need to perform
metagenomic binning, clustering the assembled contigs into draft genomes.
Despite the existing computational tools, most of them neglect one important
property of the metagenomic data, that is, the noise. To further improve the
metagenomic binning step and reconstruct better metagenomes, we propose a deep
Contrastive Learning framework for Metagenome Binning (CLMB), which can
efficiently eliminate the disturbance of noise and produce more stable and
robust results. Essentially, instead of denoising the data explicitly, we add
simulated noise to the training data and force the deep learning model to
produce similar and stable representations for both the noise-free data and the
distorted data. Consequently, the trained model will be robust to noise and
handle it implicitly during usage. CLMB outperforms the previous
state-of-the-art binning methods significantly, recovering the most
near-complete genomes on almost all the benchmarking datasets (up to 17\% more
reconstructed genomes compared to the second-best method). It also improves the
performance of bin refinement, reconstructing 8-22 more high-quality genomes
and 15-32 more middle-quality genomes than the second-best result.
Impressively, in addition to being compatible with the binning refiner, single
CLMB even recovers on average 15 more HQ genomes than the refiner of VAMB and
Maxbin on the benchmarking datasets. CLMB is open-source and available at
https://github.com/zpf0117b/CLMB/.


  Viral outbreaks, such as the current COVID-19 pandemic, are commonly
described by compartmental models by means of ordinary differential equation
(ODE) systems. The parameter values of these ODE models are typically unknown
and need to be estimated based on accessible data. In order to describe
realistic pandemic scenarios with strongly varying situations, these model
parameters need to be assumed as time-dependent. While parameter estimation for
the typical case of time-constant parameters does not pose larger issues, the
determination of time-dependent parameters, e.g.~the transition rates of
compartmental models, remains notoriously difficult, in particular since the
function class of these time-dependent parameters is unknown.
  In this work, we present a novel method which utilizes the Augmented Kalman
Smoother in combination with an Expectation-Maximization algorithm to
simultaneously estimate all time-dependent parameters in an SIRD compartmental
model. This approach only requires incidence data, but no prior knowledge on
model parameters or any further assumptions on the function class of the
time-dependencies. In contrast to other approaches for the estimation of the
time-dependent reproduction number, no assumptions on the parameterization of
the serial interval distribution are required. With this method, we are able to
adequately describe COVID-19 data in Germany and to give non-parametric
model-based time course estimates for the effective reproduction number.


  Biological systems have evolved to amazingly complex states, yet we do not
understand in general how evolution operates to generate increasing genetic and
functional complexity. Molecular recognition sites are short genome segments or
peptides binding a cognate recognition target of sufficient sequence
similarity. Such sites are simple, ubiquitous modules of sequence information,
cellular function, and evolution. Here we show that recognition sites, if
coupled to a time-dependent target, can rapidly evolve to complex states with
larger code length and smaller coding density than sites recognising a static
target. The underlying fitness model contains selection for recognition, which
depends on the sequence similarity between site and target, and a uniform cost
per unit of code length. Site sequences are shown to evolve in a specific
adaptive ratchet, which produces selection of different strength for code
extensions and compressions. Ratchet evolution increases the adaptive width of
evolved sites, accelerating the adaptation to moving targets and facilitating
refinement and innovation of recognition functions. We apply these results to
the recognition of fast-evolving antigens by the human immune system. Our
analysis shows how molecular complexity can evolve as a collateral to selection
for function in a dynamic environment.


  Recently, the eco-evolutionary game theory which describes the coupled
dynamics of strategies and environment have attracted great attention. At the
same time, most of the current work is focused on the classic two-player
two-strategy game. In this work, we study multi-strategy eco-evolutionary game
theory which is an extension of the framework. For simplicity, we'll focus on
the voluntary participation Prisoner's dilemma game. For the general class of
payoff-dependent feedback dynamics, we show the conditions for the existence
and stability of internal equilibrium by using the replicator dynamics,
respectively. Where internal equilibrium points, such as, two-strategy
coexistence states, three-strategy coexistence states, persistent oscillation
states and interior saddle points. These states are determined by the relative
feedback strength and payoff matrix, and are independent of the relative
feedback speed and initial state. In particular, the three-strategy coexistence
provides a new mechanism for maintaining biodiversity in biology, ecology, and
sociology. Besides, we find that this three-strategy model return to the
persistent oscillation state of the two-strategy model when there is no
defective strategy at the initial moment.


  First Few X (FFX) studies collect household-stratified data in the early
stages of a pandemic, in order to infer severity and transmissibility of an
emerging disease. We present a Bayesian method to approximately infer
population-level transmissibility for the first time from such data; previous
studies have only inferred household-level transmissibility. To do this we
perform the inference at two levels, assuming one transmission rate parameter
for within-household infection, and another transmission rate parameter for
infection between different households. We use a simplifying assumption: that
between-household infections always occur in naive households; while still
performing full joint inference on the within-household infection parameters.
In addition, a novel technique is used to remove systematic bias when the
number of new infections per day is growing or decaying, as is common in real
outbreaks, without the need for contemporaneous estimates of the serial
interval. The method is validated on simulated data and is shown to perform
well, even when the number of infected households is relatively small.


  The elastic properties of the double-stranded DNA handles used in optical
tweezers experiments on biomolecules are customarily modeled by an extensible
worm-like chain model. Fitting such model to experimental data however is no
trivial task, as the function depends on four parameters in a highly non-linear
fashion. We hereby propose a method to bypass the fitting procedure and obtain
an empirical force vs. extension curve that accurately reproduce the elasticity
of the handles.


  This chapter will provide an overview of how characterizing free-energy
landscapes can provide insights into the biophysical properties of DNA, as well
as into the behaviour of the DNA assemblies used in the field of DNA
nanotechnology. The landscapes for these complex systems are accessible through
the use of accurate coarse-grained descriptions of DNA. Particular foci will be
the landscapes associated with DNA self-assembly and mechanical deformation,
where the latter can arise from either externally imposed forces or internal
stresses.


  To a given gene tree topology $G$ and species tree topology $S$ with leaves
labeled bijectively from a fixed set $X$, one can associate a set of ancestral
configurations, each of which encodes a set of gene lineages that can be found
at a given node of the species tree. We introduce a lattice structure on
ancestral configurations, studying the directed graphs that provide graphical
representations of lattices of ancestral configurations. For a matching gene
tree topology and species tree topology $G=S$, we present a method for defining
the digraph of ancestral configurations from the tree topology by using
iterated cartesian products of graphs. We show that a specific set of paths on
the digraph of ancestral configurations is in bijection with the set of labeled
histories -- a well-known phylogenetic object that enumerates possible temporal
orderings of the coalescences of a tree. For each of a series of tree families,
we obtain closed-form expressions for the number of labeled histories by using
this bijection to count paths on associated digraphs. Finally, we prove that
our lattice construction extends to nonmatching tree pairs, and we use it to
characterize pairs $(G,S)$ having the maximal number of ancestral
configurations for a fixed $G$. We discuss how the construction provides new
methods for performing enumerations of combinatorial aspects of gene and
species trees.


  Hypertension is the leading global cause of cardiovascular disease and
premature death. Distinct hypertension subtypes may vary in their prognoses and
require different treatments. An individual's risk for hypertension is
determined by genetic and environmental factors as well as their interactions.
In this work, we studied 911 African Americans and 1,171 European Americans in
the Hypertension Genetic Epidemiology Network (HyperGEN) cohort. We built
hypertension subtype classification models using both environmental variables
and sets of genetic features selected based on different criteria. The fitted
prediction models provided insights into the genetic landscape of hypertension
subtypes, which may aid personalized diagnosis and treatment of hypertension in
the future.


  Advances in satellite imaging and GPS tracking devices have given rise to a
new era of remote sensing and geospatial analysis. In environmental science and
conservation ecology, biotelemetric data is often high-dimensional, spatially
and/or temporally, and functional in nature, meaning that there is an
underlying continuity to the biological process of interest. GPS-tracking of
animal movement is commonly characterized by irregular time-recording of animal
position, and the movement relationships between animals are prone to sudden
change. In this paper, I propose a measure of localized mutual information
(LMI) to derive a correlation function for monitoring changes in the pairwise
association between animal movement trajectories. The properties of the LMI
measure are assessed analytically and by simulation under a variety of
circumstances. Advantages and disadvantages of the LMI measure are assessed and
an alternate measure of LMI is proposed to handle potential disadvantages. The
measure of LMI is shown to be an effective tool for detecting shifts in the
correlation of animal movements, and seasonal/phasal correlatory structure.


  Predicting a structure of an antibody from its sequence is important since it
allows for a better design process of synthetic antibodies that play a vital
role in the health industry. Most of the structure of an antibody is
conservative. The most variable and hard-to-predict part is the third
complementarity-determining region of the antibody heavy chain (CDR H3).
Lately, deep learning has been employed to solve the task of CDR H3 prediction.
However, current state-of-the-art methods are not end-to-end, but rather they
output inter-residue distances and orientations to the RosettaAntibody package
that uses this additional information alongside statistical and physics-based
methods to predict the 3D structure. This does not allow a fast screening
process and, therefore, inhibits the development of targeted synthetic
antibodies. In this work, we present an end-to-end model to predict CDR H3 loop
structure, that performs on par with state-of-the-art methods in terms of
accuracy but an order of magnitude faster. We also raise an issue with a
commonly used RosettaAntibody benchmark that leads to data leaks, i.e., the
presence of identical sequences in the train and test datasets.


  There is much concern about disruption of endocrine physiology regulated by
steroid hormones in humans, other terrestrial vertebrates and fish by
industrial chemicals, such as bisphenol A, and pesticides, such as DDT. These
endocrine-disrupting chemicals influence steroid-mediated physiology in humans
and other vertebrates by competing with steroids for receptor binding sites,
disrupting diverse responses involved in reproduction, development and
differentiation. Here I discuss that due to evolution of the progesterone
receptor (PR) and mineralocorticoid receptor (MR) after ray-finned fish and
terrestrial vertebrates diverged from a common ancestor, each receptor evolved
to respond to different steroids in ray-finned fish and terrestrial
vertebrates. In elephant shark, a cartilaginous fish, ancestral to ray-finned
fish and terrestrial vertebrates, both progesterone and
17,20dihydroxyprogesterone activate the PR. During the evolution of ray-finned
fish and terrestrial vertebrates, the PR in terrestrial vertebrates continued
responding to progesterone and evolved to weakly respond to
17,20dihydroxyprogesterone. In contrast, the physiological progestin for the PR
in zebrafish and other ray-finned fish is 17,20dihydroxyprogesterone, and
ray-finned fish PR responds weakly to progesterone. The MR in fish and
terrestrial vertebrates also diverged to have different responses to
progesterone.


  Along with recent progress in structural biology and genome biology,
structural dynamics of molecular systems including nucleic acids has attracted
attention in the context of gene regulation. Structure-function relationship is
an important topic, where physicochemical properties of nucleotides are
important as well as those of amino acids in proteins. Simulation is a useful
tool for the analysis of molecular dynamics in detail, complementary to
experiments. However, molecular simulation of nucleic acids has less developed
than that of proteins yet, partly due to the physical nature of nucleic acids.
In this review, we briefly describe the current situation and future directions
of the field, hopefully as a guide to collaboration between experiment and
computation.


  Evolution occurs in populations of reproducing individuals. It is well known
that population structure can affect evolutionary dynamics. Traditionally,
natural selection is studied between mutants that differ in reproductive rate,
but are subject to the same population structure. Here we study how natural
selection acts on mutants that have the same reproductive rate, but experience
different population structures. In our framework, mutation alters population
structure, which is given by a graph that specifies the dispersal of offspring.
Reproduction can be either genetic or cultural. Competing mutants disperse
their offspring on different graphs. A more connected graph implies higher
motility. We show that enhanced motility tends to increase an invader's
fixation probability, but there are interesting exceptions. For island models,
we show that the magnitude of the effect depends crucially on the exact layout
of the additional links. Finally, we show that for low-dimensional lattices,
the effect of altered motility is comparable to that of altered fitness: in the
limit of large population size, the invader's fixation probability is either
constant or exponentially small, depending on whether it is more or less motile
than the resident.


  Numerous traits under migration-selection balance are shown to exhibit
complex patterns of genetic architecture with large variance in effect sizes.
However, the conditions under which such genetic architectures are stable have
yet to be investigated, because studying the influence of a large number of
small allelic effects on the maintenance of spatial polymorphism is
mathematically challenging, due to the high complexity of the systems that
arise. In particular, in the most simple case of a haploid population in a
two-patch environment, while it is known from population genetics that
polymorphism at a single major-effect locus is stable in the symmetric case,
there exists no analytical predictions on how this polymorphism holds when a
polygenic background also contributes to the trait. Here we propose to answer
this question by introducing a new eco-evo methodology that allows us to take
into account the combined contributions of a major-effect locus and of a
quantitative background resulting from small-effect loci, where inheritance is
encoded according to an extension to the infinitesimal model. In a regime of
small variance contributed by the quantitative loci, we justify that traits are
concentrated around the major alleles, according to a normal distribution,
using new convex analysis arguments. This allows a reduction in the complexity
of the system using a separation of time scales approach. We predict an
undocumented phenomenon of loss of polymorphism at the major-effect locus
despite strong selection for local adaptation, because the quantitative
background slowly disrupts the rapidly established polymorphism at the
major-effect locus, which is confirmed by individual-based simulations. Our
study highlights how segregation of a quantitative background can greatly
impact the dynamics of major-effect loci by provoking migrational meltdowns.


  The Gini coefficient of the life table is a concentration index that provides
information on lifespan variation. Originally proposed by economists to measure
income and wealth inequalities, it has been widely used in population studies
to investigate variation in ages at death. We focus on a complementary
indicator, Drewnowski's index, which is as a measure of equality. We study its
mathematical properties and analyze how changes over time relate to changes in
life expectancy. Further, we identify the threshold age below which mortality
improvements are translated into decreasing lifespan variation and above which
these improvements translate into increasing lifespan inequality. We illustrate
our theoretical findings simulating scenarios of mortality improvement in the
Gompertz model. Our experiments demonstrate how Drewnowski's index can serve as
an indicator of the shape of mortality patterns. These properties, along with
our analytical findings, support studying lifespan variation alongside life
expectancy trends in multiple species.


  How self-organization leads to the emergence of structure in social
populations remains a fascinating and open question in the study of complex
systems. One frequently observed structure that emerges again and again across
systems is that of self-similar community, i.e., homophily. We use a game
theoretic perspective to explore a case in which individuals choose affiliation
partnerships based on only two factors: the value they place on having social
contacts, and their risk tolerance for exposure to threat derived from social
contact (e.g., infectious disease, threatening ideas, etc.). We show how
diversity along just these two influences are sufficient to cause the emergence
of self-organizing homophily in the population. We further consider a case in
which extrinsic social factors influence the desire to maintain particular
social ties, and show the robustness of emergent homophilic patterns to these
additional influences. These results demonstrate how observable
population-level homophily may arise out of individual behaviors that balance
the value of social contacts against the potential risks associated with those
contacts. We present and discuss these results in the context of outbreaks of
infectious disease in human populations. Complementing the standard narrative
about how social division alters epidemiological risk, we here show how
epidemiological risk may deepen social divisions in human populations.


  Inferring protein-protein interactions from sequences is an important task in
computational biology. Recent methods based on Direct Coupling Analysis (DCA)
or Mutual Information (MI) allow to find interaction partners among paralogs of
two protein families. Does successful inference mainly rely on correlations
from structural contacts or from phylogeny, or both? Do these two types of
signal combine constructively or hinder each other? To address these questions,
we generate and analyze synthetic data produced using a minimal model that
allows us to control the amounts of structural constraints and phylogeny. We
show that correlations from these two sources combine constructively to
increase the performance of partner inference by DCA or MI. Furthermore, signal
from phylogeny can rescue partner inference when signal from contacts becomes
less informative, including in the realistic case where inter-protein contacts
are restricted to a small subset of sites. We also demonstrate that
DCA-inferred couplings between non-contact pairs of sites improve partner
inference in the presence of strong phylogeny, while deteriorating it
otherwise. Moreover, restricting to non-contact pairs of sites preserves
inference performance in the presence of strong phylogeny. In a natural data
set, as well as in realistic synthetic data based on it, we find that
non-contact pairs of sites contribute positively to partner inference
performance, and that restricting to them preserves performance, evidencing an
important role of phylogeny.


  Sensitivity Analysis (SA) is a useful tool to measure the impact of changes
in model parameters on the infection dynamics, particularly to quantify the
expected efficacy of disease control strategies. SA has only been applied to
epidemic models at the population level, ignoring the effect of within-host
virus-with-immune-system interactions on the disease spread. Connecting the
scales from individual to population can help inform drug and vaccine
development. Thus the value of understanding the impact of immunological
parameters on epidemiological quantities. Here we consider an
age-since-infection structured vector-host model, in which epidemiological
parameters are formulated as functions of within-host virus and antibody
densities, governed by an ODE system. We then use SA for these
immuno-epidemiological models to investigate the impact of immunological
parameters on population-level disease dynamics such as basic reproduction
number, final size of the epidemic or the infectiousness at different phases of
an outbreak. As a case study, we consider Rift Valley Fever Disease (RFVD)
utilizing parameter estimations from prior studies. SA indicates that 1%
increase in within-host pathogen growth rate can lead up to 8% increase in R0;
up to 1% increase in steady-state infected host abundance, and up to 4%
increase in infectiousness of hosts when the reproduction number R0 is larger
than one. These significant increases in population-scale disease quantities
suggest that control strategies that reduce the within-host pathogen growth can
be important in reducing disease prevalence.


  COVID-19 pandemic is an ongoing global pandemic which has caused
unprecedented disruptions in the public health sector and global economy. The
virus, SARS-CoV-2 is responsible for the rapid transmission of coronavirus
disease. Due to its contagious nature, the virus can easily infect an
unprotected and exposed individual from mild to severe symptoms. The study of
the virus effects on pregnant mothers and neonatal is now a concerning issue
globally among civilians and public health workers considering how the virus
will affect the mother and the neonates health. This paper aims to develop a
predictive model to estimate the possibility of death for a COVID-diagnosed
mother based on documented symptoms: dyspnea, cough, rhinorrhea, arthralgia,
and the diagnosis of pneumonia. The machine learning models that have been used
in our study are support vector machine, decision tree, random forest, gradient
boosting, and artificial neural network. The models have provided impressive
results and can accurately predict the mortality of pregnant mothers with a
given input.The precision rate for 3 models(ANN, Gradient Boost, Random Forest)
is 100% The highest accuracy score(Gradient Boosting,ANN) is 95%,highest
recall(Support Vector Machine) is 92.75% and highest f1 score(Gradient
Boosting,ANN) is 94.66%. Due to the accuracy of the model, pregnant mother can
expect immediate medical treatment based on their possibility of death due to
the virus. The model can be utilized by health workers globally to list down
emergency patients, which can ultimately reduce the death rate of COVID-19
diagnosed pregnant mothers.


  In [Ecological Complexity 44 (2020) Art. 100885, DOI:
10.1016/j.ecocom.2020.100885] a continuous-time compartmental mathematical
model for the spread of the Coronavirus disease 2019 (COVID-19) is presented
with Portugal as case study, from 2 March to 4 May 2020, and the local
stability of the Disease Free Equilibrium (DFE) is analysed. Here, we propose
an analogous discrete-time model and, using a suitable Lyapunov function, we
prove the global stability of the DFE point. Using COVID-19 real data, we show,
through numerical simulations, the consistence of the obtained theoretical
results.


  The need for diverse chromosomal modifications in biotechnology, synthetic
biology and basic research requires the development of new technologies. With
CRISPR SWAPnDROP, we extend the limits of genome editing to large-scale in-vivo
DNA transfer between bacterial species. Its modular platform approach
facilitates species specific adaptation to confer genome editing in various
species. In this study, we show the implementation of the CRISPR SWAPnDROP
concept for the model organism Escherichia coli and the currently fastest
growing and biotechnologically relevant organism Vibrio natriegens. We
demonstrate the excision, transfer and integration of 151kb chromosomal DNA
between E. coli strains and from E. coli to V. natriegens without size-limiting
intermediate DNA extraction. With the transfer of the E. coli MG1655 wild type
lac operon, we establish a functional lactose and galactose degradation pathway
in V. natriegens to extend its biotechnological spectrum. We also transfer the
E. coli DH5alpha lac operon and make V. natriegens capable of
alpha-complementation - a step towards an ultra-fast cloning strain.
Furthermore, CRISPR SWAPnDROP is designed to be the swiss army knife of genome
engineering. Its spectrum of application comprises scarless, marker-free,
iterative and parallel insertions and deletions, genome rearrangements, as well
as gene transfer between strains and across species. The modular character
facilitates DNA library applications and the recycling of standardized parts.
Its novel multi-color scarless co-selection system significantly improves
editing efficiency to 92% for single edits and 83% for quadruple edits and
provides visual quality controls throughout the assembly and editing process.


  Background: Recent work showed that the temporal growth of the novel
coronavirus disease (COVID-19) follows a sub-exponential power-law scaling
whenever effective control interventions are in place. Taking this into
consideration, we present a new phenomenological logistic model that is
well-suited for such power-law epidemic growth.
  Methods: We empirically develop the logistic growth model using simple
scaling arguments, known boundary conditions and a comparison with available
data from four countries, Belgium, China, Denmark and Germany, where (arguably)
effective containment measures were put in place during the first wave of the
pandemic. A non-linear least-squares minimization algorithm is used to map the
parameter space and make optimal predictions.
  Results: Unlike other logistic growth models, our presented model is shown to
consistently make accurate predictions of peak heights, peak locations and
cumulative saturation values for incomplete epidemic growth curves. We further
show that the power-law growth model also works reasonably well when
containment and lock down strategies are not as stringent as they were during
the first wave of infections in 2020. On the basis of this agreement, the model
was used to forecast COVID-19 fatalities for the third wave in South Africa,
which is currently in progress.
  Conclusions: We anticipate that our presented model will be useful for a
similar forecasting of COVID-19 induced infections/deaths in other regions as
well as other cases of infectious disease outbreaks, particularly when
power-law scaling is observed.


  Divide-and-conquer dividing by a half recurrences, of the form $x_n =a\cdot
x_{\left\lceil{n}/{2}\right\rceil}+a\cdot
x_{\left\lfloor{n}/{2}\right\rfloor}+p(n)$, $n\geq 2$, appear in many areas of
applied mathematics, from the analysis of algorithms to the optimization of
phylogenetic balance indices. The Master Theorems that solve these equations do
not provide the solution's explicit expression, only its big-$\Theta$ order of
growth. In this paper we give an explicit expression (in terms of the binary
decomposition of $n$) for the solution $x_n$ of a recurrence of this form, with
given initial condition $x_1$, when the independent term $p(n)$ is a polynomial
in $\lceil{n}/{2}\rceil$ and $\lfloor{n}/{2}\rfloor$.


  Phylogenetic trees are used to model evolution: leaves are labelled to
represent contemporary species ("taxa") and interior vertices represent extinct
ancestors. Informally, convex characters are measurements on the contemporary
species in which the subset of species (both contemporary and extinct) that
share a given state, form a connected subtree. In \cite{KelkS17} it was shown
how to efficiently count, list and sample certain restricted subfamilies of
convex characters, and algorithmic applications were given. We continue this
work in a number of directions. First, we show how combining the enumeration of
convex characters with existing parameterised algorithms can be used to speed
up exponential-time algorithms for the \emph{maximum agreement forest problem}
in phylogenetics. Second, we re-visit the quantity $g_2(T)$, defined as the
number of convex characters on $T$ in which each state appears on at least 2
taxa. We use this to give an algorithm with running time $O( \phi^{n} \cdot
\text{poly}(n) )$, where $\phi \approx 1.6181$ is the golden ratio and $n$ is
the number of taxa in the input trees, for computation of \emph{maximum
parsimony distance on two state characters}. By further restricting the
characters counted by $g_2(T)$ we open an interesting bridge to the literature
on enumeration of matchings. By crossing this bridge we improve the running
time of the aforementioned parsimony distance algorithm to $O( 1.5895^{n} \cdot
\text{poly}(n) )$, and obtain a number of new results in themselves relevant to
enumeration of matchings on at-most binary trees.


  We consider a population distributed between two habitats, in each of which
it experiences a growth rate that switches periodically between two values, $1-
\varepsilon > 0$ or $ - (1 + \varepsilon) < 0$. We study the specific case
where the growth rate is positive in one habitat and negative in the other one
for the first half of the period, and conversely for the second half of the
period, that we refer as the $(\pm 1)$ model. In the absence of migration, the
population goes to $0$ exponentially fast in each environment. In this paper,
we show that, when the period is sufficiently large, a small dispersal between
the two patches is able to produce a very high positive exponential growth rate
for the whole population, a phenomena called inflation. We prove in particular
that the threshold of the dispersal rate at which the inflation appears is
exponentially small with the period. We show that inflation is robust to random
perturbation, by considering a model where the values of the growth rate in
each patch are switched at random times: we prove, using theory of Piecewise
Deterministic Markov Processes (PDMP) that inflation occurs for low switching
rate and small dispersal. Finally, we provide some extensions to more
complicated models, especially epidemiological and density dependent models.


  In this paper we study the long term dynamics of two prey species and one
predator species. In the deterministic setting, if we assume the interactions
are of Lotka-Volterra type (competition or predation), the long term behavior
of this system is well known. However, nature is usually not deterministic. All
ecosystems experience some type of random environmental fluctuations. We
incorporate these into a natural framework as follows. Suppose the environment
has two possible states. In each of the two environmental states the dynamics
is governed by a system of Lotka-Volterra ODE. The randomness comes from
spending an exponential amount of time in each environmental state and then
switching to the other one. We show how this random switching can create very
interesting phenomena. In some cases the randomness can facilitate the
coexistence of the three species even though coexistence is impossible in each
of the two environmental states. In other cases, even though there is
coexistence in each of the two environmental states, switching can lead to the
loss of one or more species. We look into how predators and environmental
fluctuations can mediate coexistence among competing species.


  Vaccines are the most powerful pharmaceutical tool to combat the COVID-19
pandemic. While the majority (about 65%) of the German population were fully
vaccinated, incidence started growing exponentially in October 2021 with about
41% of recorded new cases aged twelve or above being symptomatic breakthrough
infections, presumably also contributing to the dynamics. At the time, it (i)
remains elusive how significant this contribution is and (ii) whether targeted
non-pharmaceutical interventions (NPIs) may stop the amplification of the
ongoing crisis. Here, we estimate that about 67%-76% of all new infections are
caused by unvaccinated individuals, implying that only 24%-33% are caused by
the vaccinated. Furthermore, we estimate 38%-51% of new infections to be caused
by unvaccinated individuals infecting other unvaccinated individuals. In total,
unvaccinated individuals are expected to be involved in 8-9 of 10 new
infections. We further show that decreasing the transmissibility of the
unvaccinated by, e.g. targeted NPIs, causes a steeper decrease in the effective
reproduction number $\mathcal{R}$ than decreasing the transmissibility of
vaccinated individuals, potentially leading to temporary epidemic control.
Furthermore, reducing contacts between vaccinated and unvaccinated individuals
serves to decrease $\mathcal R$ in a similar manner as increasing vaccine
uptake. Taken together, our results contribute to the public discourse
regarding policy changes in pandemic response and highlight the importance of
combined measures, such as vaccination campaigns and contact reduction, to
achieve epidemic control and preventing an overload of public health systems.


  The differing ability of polypeptide conformations to act as the native state
of proteins has long been rationalized in terms of differing kinetic
accessibility or thermodynamic stability. Building on the successful
applications of physical concepts and sampling algorithms recently introduced
in the study of disordered systems, in particular artificial neural networks,
we quantitatively explore how well a quantity known as the local entropy
describes the native state of model proteins. In lattice models and all-atom
representations of proteins, we are able to efficiently sample high local
entropy states and to provide a proof of concept of enhanced stability and
folding rate. Our methods are based on simple and general
statistical--mechanics arguments, and thus we expect that they are of very
general use.


  The standard approach to estimate species trees is to align a selected set of
genes, concatenate the alignments and then estimate a consensus tree. However,
individual genes contain differing levels of evolutionary information, either
supporting or conflicting with the consensus. Based on individual gene
evolutionary tree, a recent study has demonstrated that this approach may
result in incorrect solutions and developed the internode certainty (IC)
heuristic for estimating the confidence of splits made on the consensus tree.
Although an improvement, this heuristic neglects the differing rates of
molecular evolution in individual genes. Here I develop an improved version of
this method such that each gene is proportionally weighted based on its overall
signal and specifically with the imbalanced signal for each node represented
with gene tree.


  The coronavirus Covid-19 mutates quickly in the pandemic, leaves people
struggling to verify and improve the effectiveness of the vaccine based on
biochemistry. Is there any physical invariant in the variants of such kind of
pathogen that could be taken advantage to ease the tensions? To this point,
extensive numerical experiments based on continuity mechanics were carried out
to discover the vibration modes and the range of natural frequency of
coronavirus Covid-19. Such invariant could help us in developing some flexible
technique to deactivate the coronavirus, like as resonantly breaking the viral
spike by ultrasound wave. The fundamental mechanisms governing such process are
demonstrated via solving the coupled equations of acoustics and dynamics and
thereafter the technique strategies proposed to efficiently realize the
concept.


  Let $T$ be an arbitrary phylogenetic tree with $n$ leaves. It is well-known
that the average quartet distance between two assignments of taxa to the leaves
of $T$ is $\frac 23 \binom{n}{4}$. However, a longstanding conjecture of
Bandelt and Dress asserts that $(\frac 23 +o(1))\binom{n}{4}$ is also the {\em
maximum} quartet distance between two assignments. While Alon, Naves, and
Sudakov have shown this indeed holds for caterpillar trees, the general case of
the conjecture is still unresolved. A natural extension is when partial
information is given: the two assignments are known to coincide on a given
subset of taxa. The partial information setting is biologically relevant as the
location of some taxa (species) in the phylogenetic tree may be known, and for
other taxa it might not be known. What can we then say about the average and
maximum quartet distance in this more general setting? Surprisingly, even
determining the {\em average} quartet distance becomes a nontrivial task in the
partial information setting and determining the maximum quartet distance is
even more challenging, as these turn out to be dependent of the structure of
$T$. In this paper we prove nontrivial asymptotic bounds that are sometimes
tight for the average quartet distance in the partial information setting. We
also show that the Bandelt and Dress conjecture does not generally hold under
the partial information setting. Specifically, we prove that there are cases
where the average and maximum quartet distance substantially differ.


  Diffusion theory is a central tool of modern population genetics, yielding
simple expressions for fixation probabilities and other quantities that are not
easily derived from the underlying Wright-Fisher model. Unfortunately, the
textbook derivation of diffusion equations as scaling limits requires
evolutionary parameters (selection coefficients, mutation rates) to scale like
the inverse population size -- a severe restriction that does not always
reflect biological reality. Here we note that the Wright-Fisher model can be
approximated by diffusion equations under more general conditions, including in
regimes where selection and/or mutation are strong compared to genetic drift.
As an illustration, we use a diffusion approximation of the Wright-Fisher model
to improve estimates for the expected time to fixation of a strongly
deleterious allele, i.e. the rate of Muller's ratchet.


  African Swine Fever (ASF) represents the main threat to swine production,
with heavy economic consequences for both farmers and the food industry. The
spread of the virus that causes ASF through Europe raises the issues of
identifying transmission routes and assessing their relative contributions in
order to provide insights to stakeholders for adapted surveillance and control
measures. A simulation model was developed to assess ASF spread over the
commercial swine network in France. The model was designed from raw movement
data and actual farm characteristics. A metapopulation approach was used, with
transmission processes at the herd level potentially leading, through a
reaction-diffusion process, to external spread to epidemiologically connected
herds. Three transmission routes were considered: local transmission (e.g.
fomites, material exchange), movement of animals from infected to susceptible
sites, and transit of trucks without physical animal exchange. Surveillance was
based on prevalence and mortality detection thresholds, which triggered control
measures based on movement ban for detected herds and epidemiologically related
herds. The time from infection to detection varied between 8 and 21 days,
depending on the detection criteria, but was also dependent on the types of
herds in which the infection was introduced. Movement restrictions effectively
reduced the transmission between herds, but local transmission was nevertheless
observed in higher proportions highlighting the need of global awareness of all
actors of the swine industry to mitigate the risk of local spread. Raw movement
data were directly used to build a dynamic network on a realistic time-scale.
This approach allows for a rapid update of input data without any
pre-treatment, which could be important in terms of reactivity, should an
introduction occur.


  Biological signaling pathways based upon proteins binding to one another to
relay a signal for genetic expression, such as the Bone Morphogenetic Protein
(BMP) signaling pathway, can be modeled by mass action kinetics and
conservation laws that result in non-closed form polynomial equations.
Accurately determining parameters of biological pathways that represent
physically relevant features, such as binding affinity of proteins and their
associated uncertainty, presents a challenge for biological models lacking an
explicit likelihood function. Additionally, parameterizing non-closed form
biological models requires copious amounts of data from expensive
perturbation-response experiments to fit model parameters. We present an
algorithm (SBIDOEMAN) for determining optimal experiments and parameters of
systems biology models with implicit likelihoods. We evaluate our algorithm
using simulations of held-out true parameter values and demonstrate an
improvement in the rate of accurate parameter inference over random and
equidistant experimental designs when evaluated on two simple models of the BMP
signaling pathway with an implicit likelihood function.


  In the evolutionary dynamics of a rock-paper-scissor (RPS) model, the effect
of natural death plays a major role in determining the fate of the system.
Coexistence, being an unstable fixed point of the model becomes very sensitive
towards this parameter. In order to study the effect of mobility in such a
system which has explicit dependence on mortality, we perform Monte Carlo
simulation on a $2$-dimensional lattice having three cyclically-competing
species. The spatio-temporal dynamics has been studied along with the two-site
correlation function. Spatial distribution exhibits emergence of spiral
patterns in the presence of diffusion. It reveals that the joint effect of
death rate and diffusion (migration) leads to new coexistence and extinction
scenarios.


  In this work, we consider the motion of a polaron in a polynucleotide
Holstein molecular chain in a constant electric field. It is shown that the
character of the polaron motion in the chain depends not only on the chosen
parameters of the chain, but also on the initial distribution of the charge
along the chain. It is shown that for a small set value of the electric field
intensity and for fixed values of the chain parameters, changing only the
initial distribution of the charge in the chain, it is possible to observe
either a uniform movement of the charge along the chain, or an oscillatory mode
of charge movement.


  With the rapid development of high-throughput experimental technologies,
different types of omics (e.g., genomics, epigenomics, transcriptomics,
proteomics, and metabolomics) data can be produced from clinical samples. The
correlations between different omics types attracts a lot of research interest,
whereas the stduy on genome-wide omcis data translation (i.e, generation and
prediction of one type of omics data from another type of omics data) is almost
blank. Generative adversarial networks and the variants are one of the most
state-of-the-art deep learning technologies, which have shown great success in
image-to-image translation, text-to-image translation, etc. Here we proposed
OmiTrans, a deep learning framework adopted the idea of generative adversarial
networks to achieve omics-to-omics translation with promising results. OmiTrans
was able to faithfully reconstruct gene expression profiles from DNA
methylation data with high accuracy and great model generalisation, as
demonstrated in the experiments.


  We introduce a population-age-time (PAT) model which describes the temporal
evolution of the population distribution in age. The surprising result is that
the qualitative nature of the population distribution dynamics is robust with
respect to the birth rate and death rate distributions in age, and initial
conditions. When the number of children born per woman is 2, the population
distribution approaches an asymptotically steady state of a kink shape; thus
the total population approaches a constant. When the number of children born
per woman is greater than 2, the total population increases without bound; and
when the number of children born per woman is less than 2, the total population
decreases to zero.


  Binding of ligands is often crucial for function yet the effects of ligand
binding on the mechanical stability and energy landscape of proteins are
incompletely understood. Here we use a combination of single-molecule optical
tweezers and MD simulations to investigate the effect of ligand binding on the
energy landscape of acyl-coenzyme A (CoA) binding protein (ACBP). ACBP is a
topologically simple and highly conserved four-alpha-helix bundle protein that
acts as an intracellular transporter and buffer for fatty-acyl CoA and is
active in membrane assembly. We have previously described the behavior of ACBP
under tension, revealing a highly extended transition state (TS) located almost
halfway between the unfolded and native states. Here, we performed force-ramp
and force-jump experiments, in combination with advanced statistical analysis,
to show that octanoyl-CoA binding increases the activation free energy for the
unfolding reaction of ACBP without affecting the position of the transition
state along the reaction coordinate. It follows that ligand binding enhances
the mechanical resistance and thermodynamic stability of the protein, without
changing its mechanical compliance. Steered molecular dynamics simulations
allowed us to rationalize the results in terms of key interactions that
octanoyl-CoA establishes with the four alpha-helices of ACBP and showed that
the unfolding pathway is marginally affected by the ligand. The results show
that ligand-induced mechanical stabilization effects can be complex and may
prove useful for the rational design of stabilizing ligands.


  In this paper, we utilized generative models, and reformulate it for problems
in molecular dynamics (MD) simulation, by introducing an MD potential energy
component to our generative model. By incorporating potential energy as
calculated from TorchMD into a conditional generative framework, we attempt to
construct a low-potential energy route of transformation between the
helix~$\rightarrow$~coil structures of a protein. We show how to add an
additional loss function to conditional generative models, motivated by
potential energy of molecular configurations, and also present an optimization
technique for such an augmented loss function. Our results show the benefit of
this additional loss term on synthesizing realistic molecular trajectories.


  In iterated games, a player can unilaterally exert influence over the outcome
through a careful choice of strategy. A powerful class of such "payoff control"
strategies was discovered by Press and Dyson (2012). Their so-called
"zero-determinant" (ZD) strategies allow a player to unilaterally enforce a
linear relationship between both players' payoffs. It was subsequently shown by
Chen and Zinger (2014) that when the slope of this linear relationship is
positive, ZD strategies are robustly effective against a selfishly optimizing
co-player, in that all adapting paths of the selfish player lead to the maximal
payoffs for both players (at least when there are certain restrictions on the
game parameters). In this paper, we investigate the efficacy of selfish
learning against a fixed player in more general settings, for both ZD and
non-ZD strategies. We first prove that in any symmetric $2\times 2$ game, the
selfish player's final strategy must be of a certain form and cannot be fully
stochastic. We then show that there are prisoner's dilemma interactions for
which selfish optimization does not always lead to maximal payoffs against
fixed ZD strategies with positive slope. We give examples of selfish adapting
paths that lead to locally but not globally optimal payoffs, undermining the
robustness of payoff control strategies. For non-ZD strategies, these
pathologies arise regardless of the original restrictions on the game
parameters. Our results illuminate the difficulty of implementing robust payoff
control and selfish optimization, even in the simplest context of playing
against a fixed strategy.


  An effective model for protein structures is important for the study of
protein geometry, which, to a large extent, determine the functions of
proteins. There are a number of approaches for modelling; one might focus on
the conformation of the backbone or H-bonds, and the model may be based on the
geometry or the topology of the structure in focus. We focus on the topology of
H-bonds in proteins, and explore the link between the topology and the geometry
of protein structures. More specifically, we take inspiration from CASP
Evaluation of Model Accuracy and investigate the extent to which structural
similarities, via GDT_TS, can be estimated from the topology of H-bonds. We
report on two experiments; one where we attempt to mimic the computation of
GDT_TS based solely on the topology of H-bonds, and the other where we perform
linear regression where the independent variables are various scores computed
from the topology of H-bonds. We achieved an average $\Delta\text{GDT}$ of 6.45
with 54.5% of predictions inside 2 $\Delta\mathrm{GDT}$ for the first method,
and an average $\Delta\mathrm{GDT}$ of 4.41 with 72.7% of predictions inside 2
$\Delta\mathrm{GDT}$ for the second method.


  We introduce a new, simplified model of proteins, which we call protein
metastructure. The metastructure of a protein carries information about its
secondary structure and $\beta$-strand conformations. Furthermore, protein
metastructure allows us to associate an object called a fatgraph to a protein,
and a fatgraph in turn gives rise to a topological surface. It becomes thus
possible to study the topological invariants associated to a protein. We
discuss the correspondence between protein metastructures and fatgraphs, and
how one can compute topological invariants, such as genus and the number of
boundary components, from fatgraphs. We then describe an algorithm for
generating likely candidate metastructures using the information obtained from
topology of protein fatgraphs. This algorithm is further developed to predict
$\beta$-sheet topology of proteins, with a possibility to combine it with an
existing algorithm. We demonstrate the algorithm on the data from PDB, and
improve the performance of and existing algorithm by combining with it.


  Stochastic reaction networks is a powerful class of models for the
representation a wide variety of population models including biochemistry. The
control of such networks has been recently considered due to their important
implications for the control of biological systems. Their optimal control,
however, has been relatively few studied until now. The continuous-time
finite-horizon optimal control problem is formulated first and explicitly
solved in the case of unimolecular reaction networks. The problems of the
optimal sampled-data control, the continuous $H_\infty$ control, and the
sampled-data $H_\infty$ control of such networks are addressed next. The
results in the unimolecular case take the form of nonstandard Riccati
differential equations or differential Lyapunov equations coupled with
difference Riccati equations, which can all be solved numerically by
backward-in-time integration.


  A proper understanding of the links between varying gene expression levels
and complex trait adaptation is still lacking, despite recent advances in
sequencing techniques leading to new insights on their importance in some
evolutionary processes. This calls for extensions of the continuum-of-alleles
framework first introduced by Kimura (1965) that bypass the classical Gaussian
approximation. Here, we propose a novel mathematical framework to study the
evolutionary dynamics of quantitative alleles for sexually reproducing
populations under natural selection and competition through an
integro-differential equation. It involves a new reproduction operator which
allows to use techniques based on the maximum principle, unlike the
infinitesimal model operator used in other studies with sexual reproduction. In
an asymptotic regime where initially the population has a small phenotypic
variance, we analyse the long-term dynamics of the phenotypic distributions
according to the methodology of small variance (Diekmann et al. 2005). Under
some assumptions on the limit equation, we show that the population remains
monomorphic, that is the phenotypic distribution remains concentrated as a
moving Dirac mass. Moreover, in the case of a monomorphic distribution we
derive a canonical equation describing the dynamics of the dominant alleles.


  The standard models of sequence evolution on a tree determine probabilities
for every character or site pattern. A flattening is an arrangement of these
probabilities into a matrix, with rows corresponding to all possible site
patterns for one set $A$ of taxa and columns corresponding to all site patterns
for another set $B$ of taxa. Flattenings have been used to prove difficult
results relating to phylogenetic invariants and consistency and also form the
basis of several methods of phylogenetic inference. We prove that the rank of
the flattening equals $r^{\ell_T(A|B)}$, where $r$ is the number of states and
$\ell_T(A|B)$ is the parsimony length of the binary character separating $A$
and $B$. This result corrects an earlier published formula and opens up new
applications for old parsimony theorems.
  Since completing this work, we have learnt that an equivalent result has been
proved much earlier by Casanellas and Fern\'andez-S\'anchez, using a different
proof strategy.


  This study focuses on open quasispecies systems with competition and death
flow, described by modified Eigen and Crow-Kimura models. We examine the
evolutionary adaptation process as a reaction to changes in rates. One of the
fundamental assumptions, which forms the basis of our mathematical model, is
the existence of two different timescales: internal dynamics time and
evolutionary time. The latter is much slower and exhibits significant
adaptation events. These conditions allow us to represent the whole
evolutionary process through a series of steady-state equations, where all the
elements continuously depend on the evolutionary parameter.


  With the increasingly available large-scale cancer genomics datasets, machine
learning approaches have played an important role in revealing novel insights
into cancer development. Existing methods have shown encouraging performance in
identifying genes that are predictive for cancer survival, but are still
limited in modeling the distribution over genes. Here, we proposed a novel
method that can simulate the gene expression distribution at any given time
point, including those that are out of the range of the observed time points.
In order to model the irregular time series where each patient is one
observation, we integrated a neural ordinary differential equation (neural ODE)
with cox regression into our framework. We evaluated our method on eight cancer
types on TCGA and observed a substantial improvement over existing approaches.
Our visualization results and further analysis indicate how our method can be
used to simulate expression at the early cancer stage, offering the possibility
for early cancer identification.


  Protein-ligand interactions (PLIs) are fundamental to biochemical research
and their identification is crucial for estimating biophysical and biochemical
properties for rational therapeutic design. Currently, experimental
characterization of these properties is the most accurate method, however, this
is very time-consuming and labor-intensive. A number of computational methods
have been developed in this context but most of the existing PLI prediction
heavily depends on 2D protein sequence data. Here, we present a novel parallel
graph neural network (GNN) to integrate knowledge representation and reasoning
for PLI prediction to perform deep learning guided by expert knowledge and
informed by 3D structural data. We develop two distinct GNN architectures, GNNF
is the base implementation that employs distinct featurization to enhance
domain-awareness, while GNNP is a novel implementation that can predict with no
prior knowledge of the intermolecular interactions. The comprehensive
evaluation demonstrated that GNN can successfully capture the binary
interactions between ligand and proteins 3D structure with 0.979 test accuracy
for GNNF and 0.958 for GNNP for predicting activity of a protein-ligand
complex. These models are further adapted for regression tasks to predict
experimental binding affinities and pIC50 is crucial for drugs potency and
efficacy. We achieve a Pearson correlation coefficient of 0.66 and 0.65 on
experimental affinity and 0.50 and 0.51 on pIC50 with GNNF and GNNP,
respectively, outperforming similar 2D sequence-based models. Our method can
serve as an interpretable and explainable artificial intelligence (AI) tool for
predicted activity, potency, and biophysical properties of lead candidates. To
this end, we show the utility of GNNP on SARS-Cov-2 protein targets by
screening a large compound library and comparing our prediction with the
experimentally measured data.


  Structural DNA nanotechnology has advanced to the extent that extremely
complex structures can be designed. Much of this advancement has been due to
the development of automated DNA design and simulation tools. Typically, the
tools (e.g. NUPAK, cadnano, OxDNA) are created for a specific task. Ideally,
there would be an environment that can integrate all such DNA tools, also with
non-DNA tools - for example for modelling electromagnetic field along a
zero-mode waveguide made of gold nanoparticles organized on a DNA breadboard.
Such an environment would streamline design in DNA nanotechnology and enable
applying DNA nanotechnology principles to construct high performance materials
and devices from non-DNA components.
  Here we present small a programmatic tool that is a step towards building
such an environment for designing arbitrary nanostructures. In particular we
showcase how small has been used to create an integrated computational
materials engineering (ICME) framework for DNA nanotechnology, allowing the
hierarchical design, simulation and visualization of arbitrary DNA
nanostructures. Furthermore we demonstrate the design and modeling of the mode
profiles and band structure of hybrid DNA-nanoparticle materials through the
integration of small with Maxwell solvers.


  We introduce a correspondence between phylogenetic trees and Brauer diagrams,
inspired by links between binary trees and matchings described by Diaconis and
Holmes (1998). This correspondence gives rise to a range of semigroup
structures on the set of phylogenetic trees, and opens the prospect of many
applications. We furthermore extend the Diaconis-Holmes correspondence from
binary trees to non-binary trees and to forests, showing for instance that the
set of all forests is in bijection with the set of partitions of finite sets.


  We study a non-Markovian and nonstationary model of animal mobility
incorporating both exploration and memory in the form of preferential returns.
We derive exact results for the probability of visiting a given number of sites
and develop a practical WKB approximation to treat the nonstationary problem.
We further show that this model adequately describes empirical movement data of
Egyptian fruit bats (Rousettus aegyptiacus) when accounting for
inter-individual variation in the population. Finally, we study the probability
of visiting any site a given number of times and derive the corresponding
mean-field equation. Here, we find a remarkable phase transition occurring at
preferential returns which scale linearly with past visits. Following empirical
evidence, we suggest that this phase transition reflects a trade-off between
extensive and intensive foraging modes.


  Motivation: The application of PacBio HiFi and ultra-long ONT reads have
achieved huge progress in the contig-level assembly, but it is still
challenging to assemble large contigs into chromosomes with available Hi-C
scaffolding software, which all compute the contact value between contigs using
the Hi-C links from the whole contig regions. As the Hi-C links of two adjacent
contigs concentrate only at the neighbor ends of the contigs, larger contig
size will reduce the power to differentiate adjacent (signal) and non-adjacent
(noise) contig linkages, leading to a higher rate of mis-assembly.
  Results: We present a software package EndHiC, which is suitable to assemble
large contigs (> 1-Mb) into chromosomal-level scaffolds, using Hi-C links from
only the contig end regions instead of the whole contig regions. Benefiting
from the increased signal to noise ratio, EndHiC achieves much higher
scaffolding accuracy compared to existing software LACHESIS, ALLHiC, and
3D-DNA. Moreover, EndHiC has few parameters, runs 10-1000 times faster than
existing software, needs trivial memory, provides robustness evaluation, and
allows graphic viewing of the scaffold results. The high scaffolding accuracy
and user-friendly interface of EndHiC, liberate the users from labor-intensive
manual checks and revision works.
  Availability and implementation: EndHiC is written in Perl, and is freely
available at https://github.com/fanagislab/EndHiC. Contact: fanwei@caas.cn and
milrazhang@163.com Supplementary information: Supplementary data are available
at Bioinformatics online.


  Marine ecosystem models are an important tool to assess the role of the ocean
biota in climate change and to identify relevant biogeochemical processes by
validating the model outputs against observational data. For the assessment of
the marine ecosystem models, the existence and uniqueness of an annual periodic
solution (i.e., a steady annual cycle) is desirable. To analyze the uniqueness
of a steady annual cycle, we performed a larger number of simulations starting
from different initial concentrations for a hierarchy of biogeochemical models
with an increasing complexity. The numerical results suggested that the
simulations finished always with the same steady annual cycle regardless of the
initial concentration. Due to numerical instabilities, some inadmissible
approximations of the steady annual cycle, however, occurred in some cases for
the three most complex biogeochemical models. Our numerical results indicate a
unique steady annual cycle for practical applications.


  Parameter identification for marine ecosystem models is important for the
assessment and validation of marine ecosystem models against observational
data. The surrogate-based optimization (SBO) is a computationally efficient
method to optimize complex models. SBO replaces the computationally expensive
(high-fidelity) model by a surrogate constructed from a less accurate but
computationally cheaper (low-fidelity) model in combination with an appropriate
correction approach, which improves the accuracy of the low-fidelity model. To
construct a computationally cheap low-fidelity model, we tested three different
approaches to compute an approximation of the annual periodic solution (i.e., a
steady annual cycle) of a marine ecosystem model: firstly, a reduced number of
spin-up iterations (several decades instead of millennia), secondly, an
artificial neural network (ANN) approximating the steady annual cycle and,
finally, a combination of both approaches. Except for the low-fidelity model
using only the ANN, the SBO yielded a solution close to the target and reduced
the computational effort significantly. If an ANN approximating appropriately a
marine ecosystem model is available, the SBO using this ANN as low-fidelity
model presents a promising and computational efficient method for the
validation.


  Broadly neutralizing antibodies (bNAbs) are promising targets for vaccination
and therapy against HIV. Passive infusions of bNAbs have shown promise in
clinical trials as a potential alternative for anti-retroviral therapy. A key
challenge for the potential clinical application of bnAbs is the suppression of
viral escape, which is more effectively achieved with a combination of bNAbs.
However, identifying an optimal bNAb cocktail is combinatorially complex. Here,
we propose a computational approach to predict the efficacy of a bNAb therapy
trial based on the population genetics of HIV escape, which we parametrize
using high-throughput HIV sequence data from a cohort of untreated bNAb-naive
patients. By quantifying the mutational target size and the fitness cost of
HIV-1 escape from bNAbs, we reliably predict the distribution of rebound times
in three clinical trials. Importantly, we show that early rebounds are
dominated by the pre-treatment standing variation of HIV-1 populations, rather
than spontaneous mutations during treatment. Lastly, we show that a cocktail of
three bNAbs is necessary to suppress the chances of viral escape below 1%, and
we predict the optimal composition of such a bNAb cocktail. Our results offer a
rational design for bNAb therapy against HIV-1, and more generally show how
genetic data could be used to predict treatment outcomes and design new
approaches to pathogenic control.


  In this paper, we propose a novel free boundary problem to model the movement
of single species with a range boundary. The spatial movement and birth/death
processes of the species found within the range boundary are assumed to be
governed by the classic Fisher-KPP reaction-diffusion equation, while the
movement of a free boundary describing the range limit is assumed to be
influenced by the weighted total population inside the range boundary and is
described by an integro-differential equation. Our free boundary equation is a
generalization of the classical Stefan problem that allows for nonlocal
influences on the boundary movement so that range expansion and shrinkage are
both possible. In this paper we prove that the new model is well posed and
possesses steady state. We show that the spreading speed of the range boundary
is smaller than that for the equivalent problem with a Stefan condition. This
implies that the nonlocal effect of the weighted total population on the
boundary movement slows down the spreading speed of the population. While the
classical Stefan condition categorizes asymptotic behavior via a
spreading-vanishing dichotomy, the new model extends this dichotomy to a
spreading-balancing-vanishing trichotomy. We specifically analyze how habitat
boundaries expand, balance or shrink. When the model is extended to have two
free boundaries, we observe the steady state scenario, asymmetric shifts, or
even boundaries moving synchronously in the same direction. These are newly
discovered phenomena in the free boundary problems for animal movement.


  Recently, Wang and Xu [ Appl. Math. Lett. 78 (2018) 105-111] studied
thresholds and bi-stability in virus-immune dynamics. In this paper, we show
there also exist backward bifurcation and saddle node bifurcation in this
model. Our investigation demonstrates the existence of post-bifurcation
phenomenon in the system when the immune strength was selected as the
bifurcation parameter. The bifurcation of saddle knot is then analyzed. The
values of backward bifurcation point and saddle node bifurcation point are the
elite control threshold and post-treatment control threshold respectively. The
time interval between the two thresholds is the bi-stable interval.


  The capability of accurate prediction of protein functions and properties is
essential in the biotechnology industry, e.g. drug development and artificial
protein synthesis, etc. The main challenges of protein function prediction are
the large label space and the lack of labeled training data. Our method
leverages unsupervised sequence embedding and the success of deep convolutional
neural network to overcome these challenges. In contrast, most of the existing
methods delete the rare protein functions to reduce the label space.
Furthermore, some existing methods require additional bio-information (e.g.,
the 3-dimensional structure of the proteins) which is difficult to be
determined in biochemical experiments. Our proposed method significantly
outperforms the other methods on the publicly available benchmark using only
protein sequences as input. This allows the process of identifying protein
functions to be sped up.


  Many current challenges involve understanding the complex dynamical interplay
between the constituents of systems. Typically, the number of such constituents
is high, but only limited data sources on them are available. Conventional
dynamical models of complex systems are rarely mathematically tractable and
their numerical exploration suffers both from computational and data
limitations. Here we review generalized modeling, an alternative approach to
formulating dynamical models. We argue that this approach deals elegantly with
the uncertainties that exist in real world data and enables analytical insight
or highly efficient numerical investigation. We provide a survey of recent
successes of generalized modeling and a guide to the application of this
modeling approach in future studies such as complex integrative ecological
models.


  Horizontal gene transfer events partition a gene tree $T$ and thus, its leaf
set into subsets of genes whose evolutionary history is described by speciation
and duplication events alone. Indirect phylogenetic methods can be used to
infer such partitions $\mathcal{P}$ from sequence similarity or evolutionary
distances without any a priory knowledge about the underlying tree $T$. In this
contribution, we assume that such a partition $\mathcal{P}$ of a set of genes
$X$ is given and that, independently, an estimate $T$ of the original gene tree
on $X$ has been derived. We then ask to what extent $T$ and the xenology
information, i.e., $\mathcal{P}$ can be combined to determine the horizontal
transfer edges in $T$. We show that for each pair of genes $x$ and $y$ with
$x,y$ being in different parts of $\mathcal{P}$, it can be decided whether
there always exists or never exists a horizontal gene transfer in $T$ along the
path connecting $y$ and the most recent common ancestor of $x$ and $y$. This
problem is equivalent to determining the presence or absence of the directed
edge $(x,y)$ in so-called Fitch graphs; a more fine-grained version of graphs
that represent the dependencies between the sets in $\mathcal{P}$. We then
consider the generalization to insufficiently resolved gene trees and show that
analogous results can be obtained. We show that the classification of $(x,y)$
can be computed in constant time after linear-time preprocessing. Using
simulated gene family histories, we observe empirically that the vast majority
of horizontal transfer edges in the gene tree $T$ can be recovered
unambiguously.


  Cell-penetrating peptides (CPPs) have been widely used for drug-delivery
agents; however, it has not been fully understood how they translocate across
cell membranes. The Weighted Ensemble (WE) method, one of powerful and flexible
path sampling techniques, can be helpful to reveal translocation paths and free
energy barriers along those paths. Within the WE approach we show how Arg9
(nona-arginine) and Tat interact with a DOPC/DOPG (4:1) model membrane, and we
present free energy (or potential mean of forces, PMFs) profiles of
penetration, although a translocation across the membrane has not been observed
in the current simulations. Two different compositions of lipid molecules were
also tried and compared. Our approach can be applied to any CPPs interacting
with various model membranes, and it will provide useful information regarding
the transport mechanisms of CPPs.


  Summary: HTSeq 2.0 provides a more extensive API including a new
representation for sparse genomic data, enhancements in htseq-count to suit
single cell omics, a new script for data using cell and molecular barcodes,
improved documentation, testing and deployment, bug fixes, and Python 3
support. Availability and implementation: HTSeq 2.0 is released as an
open-source software under the GNU General Public Licence and available from
the Python Package Index at https://pypi.python.org/pypi/HTSeq. The source code
is available on Github at https://github.com/htseq/htseq. Contact:
fabio.zanini@unsw.edu.au


  Computational models of infectious disease can be broadly categorized into
two types: individual-based (Agent-based), or compartmental models. While
compartmental models can be structured to separate distinct sectors of a
population, they are conceptually distinct from individual-based models in
which population structure emerges from micro-scale interactions. While the
conceptual distinction is straightforward, a fair comparison of the approaches
is difficult to achieve. Here, we carry out such a comparison by building a set
of compartmental metapopulation models from an agent-based representation of a
real population. By adjusting the compartmental model to approximately match
the dynamics of the Agent-based model, we identify two key qualitative
properties of the individual-based dynamics which are lost upon aggregation
into metapopulations. These are (1) the local depletion of susceptibility to
infection, and (2) decoupling of different regional groups due to correlation
between commuting behaviors and contact rates. The first of these effects is a
general consequence of aggregating small, closely connected groups (i.e.,
families) into larger homogeneous metapopulations. The second can be
interpreted as a consequence of aggregating two distinct types of individuals:
school children, who travel short distances but have many potentially
infectious contacts, and adults, who travel further but tend to have fewer
contacts capable of transmitting infection. Our results could be generalised to
other types of correlations between the characteristics of individuals and the
behaviors that distinguish them.


  The basic social dilemma is frequently captured by a public goods game where
participants decide simultaneously whether to support a common pool or not and
after the enhanced contributions are distributed uniformly among all
competitors. What if the result of common efforts is {\it not} distributed
immediately, but it is reinvested and added to the pool for a next round? This
extension may not only result in an enhanced benefit for group members but also
opens new strategies for involved players because they may act in distinct
rounds differently. In this work we focus on the simplest case when two rounds
are considered, but the applied multiplication factors dedicated to a certain
round can be different. We show that in structured populations the winning
strategy may depend sensitively on the ratio of these factors and the last
round has a special importance to reach a fully cooperative state. We also
observe that it may pay for defectors to support the first round and after
enjoy the extra benefit of accumulated contributions. Full cooperator strategy
is only viable if the second round ensures a premium benefit of investments.


  Over a year after the start of the COVID-19 epidemics, we are still facing
the virus and it is hard to correctly predict its future spread over weeks to
come, as well as the impacts of potential political interventions. Current
epidemic models mainly fall in two approaches: compartmental models, divide the
population in epidemiological classes and rely on the mathematical resolution
of differential equations to give a macroscopic view of the epidemical
dynamics, allowing to evaluate its spread a posteriori; agent-based models are
computer models that give a microscopic view of the situation, since each human
is modelled as one autonomous agent, allowing to study the epidemical dynamics
in relation to (heterogeneous) individual behaviours. In this work, we compared
both methodologies and combined them to try and take advantage of the benefits
of each, and to overcome their limits. In particular, agent-based simulation
can be used to refine the values of the parameters of a compartmental model, or
to predict how these values evolve depending on sanitary policies applied. In
this report we discuss the conditions of such a combination of approaches, and
future improvements.


  Motivation: Modules in gene coexpression networks (GCN) can be regarded as
gene groups with individual relationships. No studies have optimized module
detection methods to extract diverse gene groups from GCN, especially for data
from clinical specimens. Results: Here, we optimized the flow from
transcriptome data to gene modules, aiming to cover diverse gene relationships.
We found the prediction accuracy of relationships in benchmark networks of
non-mammalian was not always suitable for evaluating gene relationships of
human and employed network based metrics. We also proposed a module detection
method involving a combination of graphical embedding and recursive
partitioning, and confirmed its stable and high performance in biological
plausibility of gene groupings. Analysis of differentially ex-pressed genes of
several reported cancers using the extracted modules successfully added
relational information consistent with previous reports, confirming the
usefulness of our framework.


  Rare variants are hypothesized to be largely responsible for heritability and
susceptibility to disease in humans. So rare variants association studies hold
promise for understanding disease. Conversely though, the rareness of the
variants poses practical challenges; since these variants are present in few
individuals, it can be difficult to develop data-collection and statistical
methods that effectively leverage their sparse information. In this work, we
develop a novel Bayesian nonparametric model to capture how design choices in
rare variants association studies can impact their usefulness. We then show how
to use our model to guide design choices under a fixed experimental budget in
practice. In particular, we provide a practical workflow and illustrative
experiments on simulated data.