  Exact results for product-form stationary distributions of Markov chains are
of interest in different fields. In stochastic reaction networks (CRNs),
stationary distributions are mostly known in special cases where they are of
product-form. However, there is no full characterization of the classes of
networks whose stationary distributions have product-form. We develop an
algebraic approach to product-form stationary distributions in the framework of
CRNs. Under certain hypotheses on linearity and decomposition of the state
space for conservative ergodic CRNs, this gives sufficient and necessary
algebraic conditions for product-form stationary distributions. Correspondingly
we obtain a semialgebraic subset of the parameter space that captures rates
where, under the corresponding hypotheses, CRNs have product-form. We employ
the developed theory to CRNs and some models of statistical mechanics, besides
sketching the pertinence in other models from applied probability.


  Background: The inception of next generations sequencing technologies have
exponentially increased the volume of biological sequence data. Protein
sequences, being quoted as the `language of life', has been analyzed for a
multitude of applications and inferences.
  Motivation: Owing to the rapid development of deep learning, in recent years
there have been a number of breakthroughs in the domain of Natural Language
Processing. Since these methods are capable of performing different tasks when
trained with a sufficient amount of data, off-the-shelf models are used to
perform various biological applications. In this study, we investigated the
applicability of the popular Skip-gram model for protein sequence analysis and
made an attempt to incorporate some biological insights into it.
  Results: We propose a novel $k$-mer embedding scheme, Align-gram, which is
capable of mapping the similar $k$-mers close to each other in a vector space.
Furthermore, we experiment with other sequence-based protein representations
and observe that the embeddings derived from Align-gram aids modeling and
training deep learning models better. Our experiments with a simple baseline
LSTM model and a much complex CNN model of DeepGoPlus shows the potential of
Align-gram in performing different types of deep learning applications for
protein sequence analysis.


  Recent advancements in transfer learning have made it a promising approach
for domain adaptation via transfer of learned representations. This is
especially when relevant when alternate tasks have limited samples of
well-defined and labeled data, which is common in the molecule data domain.
This makes transfer learning an ideal approach to solve molecular learning
tasks. While Adversarial reprogramming has proven to be a successful method to
repurpose neural networks for alternate tasks, most works consider source and
alternate tasks within the same domain. In this work, we propose a new
algorithm, Representation Reprogramming via Dictionary Learning (R2DL), for
adversarially reprogramming pretrained language models for molecular learning
tasks, motivated by leveraging learned representations in massive state of the
art language models. The adversarial program learns a linear transformation
between a dense source model input space (language data) and a sparse target
model input space (e.g., chemical and biological molecule data) using a k-SVD
solver to approximate a sparse representation of the encoded data, via
dictionary learning. R2DL achieves the baseline established by state of the art
toxicity prediction models trained on domain-specific data and outperforms the
baseline in a limited training-data setting, thereby establishing avenues for
domain-agnostic transfer learning for tasks with molecule data.


  We describe a mathematical model for the aggregation of starved first-stage C
elegans larvae (L1s). We propose that starved L1s produce and respond
chemotactically to two labile diffusible chemical signals, a short-range
attractant and a longer range repellent. This model takes the mathematical form
of three coupled partial differential equations, one that describes the
movement of the worms and one for each of the chemical signals. Numerical
solution of these equations produced a pattern of aggregates that resembled
that of worm aggregates observed in experiments. We also describe the
identification of a sensory receptor gene, srh-2, whose expression is induced
under conditions that promote L1 aggregation. Worms whose srh-2 gene has been
knocked out form irregularly shaped aggregates. Our model suggests this
phenotype may be explained by the mutant worms slowing their movement more
quickly than the wild type.


  Countries around the world implement nonpharmaceutical interventions (NPIs)
to mitigate the spread of COVID-19. Design of efficient NPIs requires
identification of the structure of the disease transmission network. We here
identify the key parameters of the COVID-19 transmission network for time
periods before, during, and after the application of strict NPIs for the first
wave of COVID-19 infections in Germany combining Bayesian parameter inference
with an agent-based epidemiological model. We assume a Watts-Strogatz
small-world network which allows to distinguish contacts within clustered
cliques and unclustered, random contacts in the population, which have been
shown to be crucial in sustaining the epidemic. In contrast to other works,
which use coarse-grained network structures from anonymized data, like cell
phone data, we consider the contacts of individual agents explicitly. We show
that NPIs drastically reduced random contacts in the transmission network,
increased network clustering, and resulted in a change from an exponential to a
constant regime of newcases. In this regime, the disease spreads like a wave
with a finite wave speed that depends on the number of contacts in a nonlinear
fashion, which we can predict by mean field theory. Our analysis indicates that
besides the well-known transitionbetween exponential increase and exponential
decrease in the number of new cases, NPIs can induce a transition to another,
previously unappreciated regime of constant new cases.


  Computational methods that operate on three-dimensional molecular structure
have the potential to solve important questions in biology and chemistry. In
particular, deep neural networks have gained significant attention, but their
widespread adoption in the biomolecular domain has been limited by a lack of
either systematic performance benchmarks or a unified toolkit for interacting
with molecular data. To address this, we present ATOM3D, a collection of both
novel and existing benchmark datasets spanning several key classes of
biomolecules. We implement several classes of three-dimensional molecular
learning methods for each of these tasks and show that they consistently
improve performance relative to methods based on one- and two-dimensional
representations. The specific choice of architecture proves to be critical for
performance, with three-dimensional convolutional networks excelling at tasks
involving complex geometries, graph networks performing well on systems
requiring detailed positional information, and the more recently developed
equivariant networks showing significant promise. Our results indicate that
many molecular problems stand to gain from three-dimensional molecular
learning, and that there is potential for improvement on many tasks which
remain underexplored. To lower the barrier to entry and facilitate further
developments in the field, we also provide a comprehensive suite of tools for
dataset processing, model training, and evaluation in our open-source atom3d
Python package. All datasets are available for download from
https://www.atom3d.ai .


  In the past decade, advances in genome sequencing have allowed researchers to
uncover the history of hybridization in diverse groups of species, including
our own. Although the field has made impressive progress in documenting the
extent of natural hybridization, both historical and recent, there are still
many unanswered questions about its genetic and evolutionary consequences.
Recent work has suggested that the outcomes of hybridization in the genome may
be in part predictable, but many open questions about the nature of selection
on hybrids and the biological variables that shape such selection have hampered
progress in this area. We discuss what is known about the mechanisms that drive
changes in ancestry in the genome after hybridization, highlight major
unresolved questions, and discuss their implications for the predictability of
genome evolution after hybridization.


  The movement of single kinesin molecules was observed while applying noisy
external forces that mimic intracellular active fluctuations. We found kinesin
accelerates under noise, especially when a large hindering load is added. The
behavior quantitatively conformed to a theoretical model that describes the
kinesin movement with simple two-state reactions. The universality of the
kinetic theory suggests that intracellular enzymes share a similar
noise-induced acceleration mechanism, i.e. active fluctuations in cells are not
just noise but are utilized to promote various physiological processes.


  Living organisms must respond to environmental changes. Generally, accurate
and rapid responses are provided by simple, unidirectional networks that
connect inputs with outputs. Besides accuracy and speed, biological responses
should also be robust to environmental or intracellular noise and mutations.
Furthermore, cells must also respond to unforeseen environmental changes that
have not previously been experienced, to avoid extinction prior to the
evolutionary rewiring of their networks, which takes numerous generations. We
have investigated gene regulatory networks that mutually activate or inhibit,
and have demonstrated that complex entangled networks can make appropriate
input-output relationships that satisfy the robust and adaptive responses
required for unforeseen challenges. Such entangled networks function for sloppy
and unreliable responses with low Hill coefficient reactions for the expression
of each gene. To compensate for such sloppiness, several detours in the
regulatory network exist. By taking advantage of the averaging over such
detours, the network shows a higher robustness to environmental and
intracellular noise as well as to mutations in the network, when compared to
simple unidirectional circuits. Furthermore, the appropriate response to
unforeseen challenges, allowing for functional outputs, is achieved as many
genes exhibit similar dynamic expression responses, irrespective of inputs, as
confirmed by applying dynamic time warping and dynamic mode decomposition. As
complex entangled networks are common in gene regulatory networks and global
gene expression responses are observed in microbial experiments, the present
results provide a novel design principle for cellular networks.


  Accurate prediction of the transmission of epidemic diseases such as COVID-19
is crucial for implementing effective mitigation measures. In this work, we
develop a tensor method to predict the evolution of epidemic trends for many
regions simultaneously. We construct a 3-way spatio-temporal tensor (location,
attribute, time) of case counts and propose a nonnegative tensor factorization
with latent epidemiological model regularization named STELAR. Unlike standard
tensor factorization methods which cannot predict slabs ahead, STELAR enables
long-term prediction by incorporating latent temporal regularization through a
system of discrete-time difference equations of a widely adopted
epidemiological model. We use latent instead of location/attribute-level
epidemiological dynamics to capture common epidemic profile sub-types and
improve collaborative learning and prediction. We conduct experiments using
both county- and state-level COVID-19 data and show that our model can identify
interesting latent patterns of the epidemic. Finally, we evaluate the
predictive ability of our method and show superior performance compared to the
baselines, achieving up to 21% lower root mean square error and 25% lower mean
absolute error for county-level prediction.


  We investigate the celebrated mathematical SICA model but using fractional
differential equations in order to better describe the dynamics of HIV-AIDS
infection. The infection process is modelled by a general functional response
and the memory effect is described by the Caputo fractional derivative.
Stability and instability of equilibrium points are determined in terms of the
basic reproduction number. Furthermore, a fractional optimal control system is
formulated and the best strategy for minimizing the spread of the disease into
the population is determined through numerical simulations based on the derived
necessary optimality conditions.


  In physics of living systems, a search for relationships of a few macroscopic
variables that emerge from many microscopic elements is a central issue. We
evolved gene regulatory networks so that the expression of target genes
(partial system) is insensitive to environmental changes. Then, we found the
expression levels of the remaining genes autonomously increase as a plastic
response. Negative proportionality was observed between the average changes in
target and remnant genes, reflecting reciprocity between the macroscopic
robustness of homeostatic genes and plasticity of regulator genes. This
reciprocity follows the lever principle, which was satisfied throughout the
evolutionary course, imposing an evolutionary constraint.


  Despite increasing data from population-wide sequencing studies, the risk for
recessive disorders in consanguineous partnerships is still heavily debated. An
important aspect that has not sufficiently been investigated theoretically, is
the influence of inbreeding on mutation load and incidence rates when the
population sizes change. We therefore developed a model to study these dynamics
for a wide range of growth and mating conditions. In the phase of population
expansion and shortly afterwards, our simulations show that there is a drop of
diseased individuals at the expense of an increasing mutation load for random
mating, while both parameters remain almost constant in highly consanguineous
partnerships. This explains the empirical observation in present times that a
high degree of consanguinity is associated with an increased risk of autosomal
recessive disorders. However, it also states that the higher frequency of
severe recessive disorders with developmental delay in inbred populations is a
transient phenomenon before a mutation-selection balance is reached again.


  Interaction with its immediate environment determines the ecology of an
organism. Species present in any habitat, wild or urban, may face extreme
pressure due to sudden perturbations. When such disturbances are unpredictable,
it becomes more challenging to tackle. Implementation of specific strategies is
therefore essential for different species to overcome adverse situations.
Numerous biotic and abiotic factors can alter the dynamics of a species.
Anthropogenic disturbance is one such factor that has considerable implications
and also the potential to impact species living in the proximity of human
habitats. We investigated the response of an urban adapted species to a sudden
surge in human footfall or overcrowding. Dogs (Canis lupus familiaris) living
freely in the streets of developing countries experience tremendous
anthropogenic pressure. It is known that human movement in an area can predict
the behaviour of these dogs by largely influencing their personalities. In the
current study, we observed a strong effect of high and sudden human footfall on
the abundance and behavioural activity of dogs. A decline in both the abundance
of dogs and behavioural activities was seen with the increase in human
movement. Further investigation over a restricted temporal scale revealed
reinstated behavioural activity but non-restoration of population abundance.
This provides important evidence on the extent to which humans influence the
behaviour of free-ranging dogs in urban environments.


  Bacteria, in contrast to eukaryotic cells contain two types of genes:
chromosomal genes that are fixed to the cell, and plasmids that are mobile
genes, easily shared to other cells. The sharing of plasmid genes between
individual bacteria and between bacterial lineages has contributed vastly to
bacterial evolution, allowing specialized traits to `jump ship' between one
lineage or species and the next. The benefits of this generosity from the point
of view of both recipient and plasmid are generally understood, but come at the
expense of chromosomal genes in the donor cell, which share potentially
advantageous genes with their competition while receiving no benefit. Using
both continuous models and agent based simulations, we demonstrate that
`secretive' genes which restrict horizontal gene transfer are favored over wide
range of models and parameter values. Our findings lead to a peculiar paradox:
given the obvious benefits of keeping secrets, why do bacteria share
information so freely?


  A simplified method to compute $R_t$, the Effective Reproduction Number, is
presented. The method relates the value of $R_t$ to the estimation of the
doubling time performed with a local exponential fit. The condition $R_t = 1$
corresponds to a growth rate equal to zero or equivalently an infinite doubling
time. Different assumptions on the probability distribution of the generation
time are considered. A simple analytical solution is presented in case the
generation time follows a gamma distribution.


  Understanding epistasis (genetic interaction) may shed some light on the
genomic basis of common diseases, including disorders of maximum interest due
to their high socioeconomic burden, like schizophrenia. Distance correlation is
an association measure that characterises general statistical independence
between random variables, not only the linear one. Here, we propose distance
correlation as a novel tool for the detection of epistasis from case-control
data of single-nucleotide polymorphisms (SNPs). On the methodological side, we
highlight the derivation of the explicit asymptotic distribution of the test
statistic. We show that this is the only way to obtain enough computational
speed for the method to be used in practice, in a scenario where the resampling
techniques found in the literature are impractical. Our simulations show
satisfactory calibration of significance, as well as comparable or better power
than existing methodology. We conclude with the application of our technique to
a schizophrenia genetics dataset, obtaining biologically sound insights.


  Fitting Susceptible-Infected-Recovered (SIR) models to incidence data is
problematic when not all infected individuals are reported. Assuming an
underlying SIR model with general but known distribution for the time to
recovery, this paper derives the implied differential-integral equations for
observed incidence data when a fixed fraction of newly infected individuals are
not observed. The parameters of the resulting system of differential equations
are identifiable. Using these differential equations, we develop a stochastic
model for the conditional distribution of current disease incidence given the
entire past history of reported cases. We estimate the model parameters using
Bayesian Markov Chain Monte-Carlo sampling of the posterior distribution. We
use our model to estimate the transmission rate and fraction of asymptomatic
individuals for the current Coronavirus 2019 outbreak in eight American
Countries: the United States of America, Brazil, Mexico, Argentina, Chile,
Colombia, Peru, and Panama, from January 2020 to May 2021. Our analysis reveals
that consistently, about 40-60% of the infections were not observed in the
American outbreaks. The two exception are Mexico and Peru, with acute
under-reporting in Mexico.


  Background and Objective: The massive parallel sequencing technology
facilitates new discoveries in terms of transcript differential analysis;
however, all the new findings must be validated, since the diversity of
transcript expression may impair the identification of the most relevant ones.
  Methods: The proposed RLeave algorithm (implemented in the R environment)
utilizes a combination of conventional analysis (classic edgeR) together with
other mathematical methods (Leave-one-out sample technique and Decision Trees
validation) to identify more relevant candidates to be in vitro or in silico
validated.
  Results: The RLeave protocol was tested using miRNome expression analysis of
two sample groups (diabetes mellitus and acute lymphoblastic leukemia), and
both had their most important differentially expressed miRNA confirmed by
RT-qPCR.
  Conclusion: This protocol is applicable in RNA-SEQ research, highlighting the
most relevant transcripts for in silico and/or in vitro validation.


  Systems Biology is a fundamental field and paradigm that introduces a new era
in Biology. The crux of its functionality and usefulness relies on metabolic
networks that model the reactions occurring inside an organism and provide the
means to understand the underlying mechanisms that govern biological systems.
Even more, metabolic networks have a broader impact that ranges from resolution
of ecosystems to personalized medicine.The analysis of metabolic networks is a
computational geometry oriented field as one of the main operations they depend
on is sampling uniformly points from polytopes; the latter provides a
representation of the steady states of the metabolic networks. However, the
polytopes that result from biological data are of very high dimension (to the
order of thousands) and in most, if not all, the cases are considerably skinny.
Therefore, to perform uniform random sampling efficiently in this setting, we
need a novel algorithmic and computational framework specially tailored for the
properties of metabolic networks.We present a complete software framework to
handle sampling in metabolic networks. Its backbone is a Multiphase Monte Carlo
Sampling (MMCS) algorithm that unifies rounding and sampling in one pass,
obtaining both upon termination. It exploits an improved variant of the
Billiard Walk that enjoys faster arithmetic complexity per step. We demonstrate
the efficiency of our approach by performing extensive experiments on various
metabolic networks. Notably, sampling on the most complicated human metabolic
network accessible today, Recon3D, corresponding to a polytope of dimension 5
335 took less than 30 hours. To our knowledge, that is out of reach for
existing software.


  The accurate prediction of biological features from genomic data is paramount
for precision medicine and sustainable agriculture. For decades, neural network
models have been widely popular in fields like computer vision, astrophysics
and targeted marketing given their prediction accuracy and their robust
performance under big data settings. Yet neural network models have not made a
successful transition into the medical and biological world due to the
ubiquitous characteristics of biological data such as modest sample sizes,
sparsity, and extreme heterogeneity. Here, we investigate the robustness,
generalization potential and prediction accuracy of widely used convolutional
neural network and natural language processing models with a variety of
heterogeneous genomic datasets. Mainly, recurrent neural network models
outperform convolutional neural network models in terms of prediction accuracy,
overfitting and transferability across the datasets under study. While the
perspective of a robust out-of-the-box neural network model is out of reach, we
identify certain model characteristics that translate well across datasets and
could serve as a baseline model for translational researchers.


  The main focus of this chapter is on public health control strategies which
are currently the main way to mitigate COVID-19 pandemic. We introduce and
compare compartmental models of increasing complexity for COVID-19 transmission
to describe dynamics of the disease spread. We begin by considering an SEAIR
model including basic characteristics related to COVID-19. Next, we shall pay
attention to age-structure modeling to emphasis the role of age-group
individuals on the disease spread. A Model with constant delay is also
formulated to show the impact of the latency period on the severity of
COVID-19. Since there is evidence that for COVID-19 disease, important
relationships exist between what is happening in the host and what is occurring
at the population level, we shall link the basic model to in-host dynamics
through the so-called threshold-type delay models. Finally, we will include
demographic effects to the most complex models and we will conduct rigorous
bifurcation analysis to quantify possible factors responsible for disease
progression.


  Autocatalytic systems are very often incorporated in the "origin of life"
models, a connection that has been analyzed in the context of the classical
hypercycles introduced by Manfred Eigen. We investigate the dynamics of certain
networks called bimolecular autocatalytic systems. In particular, we consider
the dynamics corresponding to the relative populations in these networks, and
show that they can be analyzed by studying well-chosen autonomous polynomial
dynamical systems. Moreover, we find that one can use results from reaction
network theory to prove persistence and permanence of several types of
bimolecular autocatalytic systems called autocatalytic recombination networks.


  Modern genomics experiments measure functional behaviors for many thousands
of DNA sequences. We suggest that, especially when these sequences are chosen
at random, it is natural to compute correlation functions between sequences and
measured behaviors. In simple models for the dependence of DNA flexibility on
sequence, for example, correlation functions can be interpreted directly as
interaction parameters. Analysis of recent experiments shows that this is
surprisingly effective, leading directly to extraction of distinct features for
DNA flexibility and predictions that are as accurate as more complex models.
This approach follows the conventional use of correlation functions in
statistical physics and connects the search for relevant DNA sequence features
to the search for relevant stimulus features in the analysis of sensory
neurons.


  SARS-CoV-2 (COVID-19), a positive single stranded RNA virus, member of corona
virus family, is spreading its tentacles across the world due to lack of drugs
at present. Being associated with cough, fever, and respiratory distress, this
disease caused more than 15 % mortality worldwide. Due to its vital role in
virus replication, Mpro/3CLpro has recently been regarded as a suitable target
for drug design. The current study focused on the inhibitory activity of
Calotropin, a component from milk of Calotropis gigantean, against Mpro protein
from SARS-CoV-2. Till date there is no work is undertaken on in-silico analysis
of this compound against Mpro of COVID-19 protein. In the present study,
molecular docking studies were conducted by using Patchdock tool. Protein
Interactions tool was used for protein interactions. The calculated parameters
such as docking score indicated effective binding of Calotropin to Mpro
protein. Interactions results indicated that, Mpro/ Calotropin complexes forms
hydrophobic interactions. Therefore, Calotropin may represent potential herbal
treatment to act as COVID-19 Mpro inhibitor. However, further research is
necessary to investigate their potential medicinal use.


  Cells have evolved a metabolic control of DNA replication to respond to a
wide range of nutritional conditions. Accumulating data suggest that this
poorly understood control depends, at least in part, on Central Carbon
Metabolism (CCM). In Bacillus subtilis , the glycolytic pyruvate kinase (PykA)
is intricately linked to replication. This 585 amino-acid-long enzyme comprises
a catalytic (Cat) domain that binds to phosphoenolpyruvate (PEP) and ADP to
produce pyruvate and ATP, and a C-terminal domain of unknown function.
Interestingly, the C-terminal domain termed PEPut interacts with Cat and is
homologous a domain that, in other metabolic enzymes, are phosphorylated at a
conserved TSH motif at the expense of PEP and ATP to drive sugar import and
catalytic or regulatory activities. To gain insights into the role of PykA in
replication, DNA synthesis was analyzed in various Cat and PEPut mutants grown
in a medium where the metabolic activity of PykA is dispensable for growth.
Measurements of replication parameters ( ori/ter ratio, C period and fork
speed) and of the pyruvate kinase activity showed that PykA mutants exhibit
replication defects resulting from side chain modifications in the PykA protein
rather than from a reduction of its metabolic activity. Interestingly, Cat and
PEPut have distinct commitments in replication: while Cat impacts positively
and negatively replication fork speed, PEPut stimulates initiation through a
process depending on Cat-PEPut interaction and growth conditions. Residues
binding to PEP and ADP in Cat, stabilizing the Cat-PEPut interaction and
belonging to the TSH motif of PEPut were found important for the commitment of
PykA in replication. In vitro , PykA affects the activities of replication
enzymes (the polymerase DnaE, helicase DnaC and primase DnaG) essential for
initiation and elongation and genetically linked to pykA . Our results thus
connect replication initiation and elongation to CCM metabolites (PEP, ATP and
ADP), critical Cat and PEPut residues and to multiple links between PykA and
the replication enzymes DnaE, DnaC and DnaG. We propose that PykA is endowed
with a moonlighting activity that senses the concentration of signaling
metabolites and interacts with replication enzymes to convey information on the
cellular metabolic state to the replication machinery and adjust replication
initiation and elongation to metabolism. This defines a new type of replication
regulator proposed to be part of the metabolic control that gates replication
in the cell cycle.


  We introduce and analyze a spatial Lotka-Volterra competition model with
local and nonlocal interactions. We study two alternative classes of nonlocal
competition that differ in how each species' characteristics determine the
range of the nonlocal interactions. In both cases, nonlocal interactions can
create spatial patterns of population densities in which highly populated
clumps alternate with unpopulated regions. This non-populated regions provide
spatial niches for a weaker competitor to establish in the community and
persist in conditions in which local models predict competitive exclusion.
Moreover, depending on the balance between local and nonlocal competition
intensity, the clumps of the weaker competitor vary from M-like structures with
higher densities of individuals accumulating at the edges of each clump to
triangular structures with most individuals occupying their centers. These
results suggest that long-range competition, through the creation of spatial
patterns in population densities, might be an important driving force behind
the rich diversity of species observed in real ecological communities.


  The phenomenological mathematical model of COVID-19 spreading is proposed for
large countries with geographical differentiation of population density.
According to the model COVID-19 spreading takes the form of several
spatio-temporal waves developing almost independently and simultaneously in
areas with different population density. The intensity of each wave is
described by a power-law dependence. The parameters of dependence are
determined by real statistical data at the initial stage of the disease spread.
The results of the model simulation were verified using statistical data for
the Republic of Belarus. Based on the developed model, a forecast calculation
was made at the end of May, 2020. The accuracy of forecasting the total number
of cases for a period of 3 months in the proposed approach was about 3 percent.


  Climate change, which is largely linked to human activities, is already
having a considerable impact on our societies. Based on current trends, climate
change is expected to accelerate in the coming decades. Beyond its impact on
the pace of natural disasters (floods, hurricanes, etc.), climate change may
have catastrophic consequences for human life and health. One of the concerns
is the increase in the transmission of viruses spread by mosquitoes. Indeed,
rising temperatures have a direct positive impact on the viability of
mosquitoes in ecosystems, leading to their abundance and thus the risk of
exposure of human populations to these pathogens. This study quantifies the
consequences of global warming on the risk of epidemics of viruses transmitted
by the Aedes Albopictus mosquito in metropolitan France. This mosquito, which
is a vector for the Dengue, Chikungunya and Zika viruses, among others, arrived
in mainland France in 2004 and has since spread throughout the country. Thanks
to the association previously established between the probability of the
presence of the mosquito and the average temperature combined with a
mathematical model, the probability of an epidemic and the number of people who
could be infected and die during a season in each department are estimated. If
there is a high degree of heterogeneity in metropolitan France, nearly 2,000
deaths per year could be expected by 2040.


  In apparent contradiction to competition theory, the number of known,
co-existing plankton species far exceeds their explicable biodiversity - a
discrepancy termed the Paradox of the Plankton. We introduce a new
game-theoretic model for competing micro-organisms in which one player consists
of all organisms of one species. The stable points for the population dynamics
in our model, known as strategic behavior distributions (SBDs), are probability
distributions of behaviors across all organisms which imply a stable population
of the species as a whole. We find that intra-specific variability is the key
characteristic that ultimately allows co-existence because the outcomes of
competitions between individuals with variable competitive abilities is
unpredictable. Our simulations based on the theoretical model show that up to
100 species can coexist for at least 10000 generations, and that even small
population sizes or species with inferior competitive ability can survive when
there is intra-specific variability. In nature, this variability can be
observed as niche differentiation, variability in environmental and ecological
factors, and variability of individual behaviors or physiology. Therefore
previous specific explanations of the paradox are consistent with and provide
specific examples of our suggestion that individual variability is the
mechanism which solves the paradox.


  Using game theory we provide mathematical proof that if a species of
asexually reproducing microbes does not possess maximum variability in
competitive abilities amongst its individual organisms, then that species is
vulnerable to replacement by competitors. Furthermore, we prove that such
maximally variable species are neutral towards each other in competition for
limited resources; they coexist. Our proof is constructive: given one species
which does not possess maximum variability, we construct a species with the
same (or lower) mean competitive ability which can invade, in the sense that
its expected value in competition is positive whereas the expected value of the
non-maximally variable species is negative. Our results point towards the
mechanistic underpinnings for the frequent observations that (1) microbes are
characterized by large intra-specific variability and that (2) the number of
extant microbe species is very large.


  Why, contrary to theoretical predictions, do marine microbe communities
harbor tremendous phenotypic heterogeneity? How can so many marine microbe
species competing in the same niche coexist? We discovered a unifying
explanation for both phenomena by investigating a non-cooperative game that
interpolates between individual-level competitions and species-level outcomes.
We identified all equilibrium strategies of the game. These strategies are
characterized by maximal phenotypic heterogeneity. They are also neutral
towards each other in the sense that an unlimited number of species can
co-exist while competing according to the equilibrium strategies. Whereas prior
theory predicts that natural selection would minimize trait variation around an
optimum value, here we obtained a rigorous mathematical proof that species with
maximally variable traits are those that endure. This discrepancy may reflect a
disparity between predictions from models developed for larger organisms in
contrast to our microbe-centric model. Rigorous mathematics proves that
phenotypic heterogeneity is itself a mechanistic underpinning of microbial
diversity. This discovery has fundamental ramifications for microbial ecology
and may represent an adaptive reservoir sheltering biodiversity in changing
environmental conditions.


  This article introduces a novel protein structure alignment method (named
TALI) based on the protein backbone torsion angle instead of the more
traditional distance matrix. Because the structural alignment of the two
proteins is based on the comparison of two sequences of numbers (backbone
torsion angles), we can take advantage of a large number of well-developed
methods such as Smith-Waterman or Needleman-Wunsch. Here we report the result
of TALI in comparison to other structure alignment methods such as DALI, CE,
and SSM ass well as sequence alignment based on PSI-BLAST. TALI demonstrated
great success over all other methods in application to challenging proteins.
TALI was more successful in recognizing remote structural homology. TALI also
demonstrated an ability to identify structural homology between two proteins
where the structural difference was due to a rotation of internal domains by
nearly 180$^\circ$.


  We consider the SEIRS compartment epidemiology model suitable for predicting
the evolution of the COVID-19 pandemy in the extreme limiting case of no
acquired immunity. The disease-free and endemic fixed points are found and
their stability is analysed. The expression for the basic reproduction ratio is
obtained and discussed, emphasizing on its dependence on the model parameters.
The threshold contact ratio is found which determines the possibility for a
stable disease-free fixed point existence. Numeric solution for the pandemy
evolution is also undertaken together with the approximate analytic solutions
for the early stage of the disease spread as well as as for its decay after the
rapid measures are undertaken. We analysed several possible scenarios for
introducing and relaxing the quarantine measures. The cyclic "quarantine on"
and "quarantine off" strategy at fixed identification and isolation ratios fail
to reduce the lowering of the second and the consecutive waves, whereas this
goal is possible to achieve if the flexible increase of the identification and
isolation ratios is also involved.


  A key question in evolution is how likely a mutant is to take over. This
depends on natural selection and on stochastic fluctuations. Population spatial
structure can impact mutant fixation probabilities. We introduce a model for
structured populations on graphs that generalizes previous ones by making
migrations independent of birth and death. We demonstrate that by tuning
migration asymmetry, the star graph transitions from amplifying to suppressing
natural selection. The results from our model are universal in the sense that
they do not hinge on a modeling choice of microscopic dynamics or update rules.
Instead, they depend on migration asymmetry, which can be experimentally tuned
and measured.


  A reasonable prediction of infectious diseases transmission process under
different disease control strategies is an important reference point for policy
makers. Here we established a dynamic transmission model via Python and
realized comprehensive regulation of disease control measures. We classified
government interventions into three categories and introduced three parameters
as descriptions for the key points in disease control, these being
intraregional growth rate, interregional communication rate, and detection rate
of infectors. Our simulation predicts the infection by COVID-19 in the UK would
be out of control in 73 days without any interventions; at the same time, herd
immunity acquisition will begin from the epicentre. After we introduced
government interventions, single intervention is effective in disease control
but at huge expense while combined interventions would be more efficient, among
which, enhancing detection number is crucial in control strategy of COVID-19.
In addition, we calculated requirements for the most effective vaccination
strategy based on infection number in real situation. Our model was programmed
with iterative algorithms, and visualized via cellular automata, it can be
applied to similar epidemics in other regions if the basic parameters are
inputted, and is able to synthetically mimick the effect of multiple factors in
infectious disease control.


  By the end of 2020, a year since the first cases of infection by the Covid-19
virus have been reported, there is a light at the end of the tunnel. Several
pharmaceutical companies made significant progress in developing effective
vaccines against the Covid-19 virus that has claimed the lives of more than a
million people over the world. On the other hand, there is growing evidence of
re-infection by the virus, which can cause further outbreaks. In this paper, we
apply statistical physics tools to examine the vaccination rate required to
control the pandemic for three different vaccine efficiency scenarios. Also, we
study the effect of temporal restrictions/reliefs on the pandemic's outbreak,
assuming that re-infection is possible. When examining the efficiency of the
vaccination rate of the general population in preventing an additional outbreak
of the disease, we find that a high vaccination rate (where at least 0.3% of
the population is vaccinated daily, which is equivalent to ~ 1 million vaccine
doses in the United States daily) is required to gain control over the spread
of the virus without further restrictions. Due to feasible limitations on the
vaccination rate, the vaccination process should be accompanied by the
prevailing restrictions until most of the population is vaccinated.


  Phylogenetic trees are leaf-labelled trees used to model the evolution of
species. Here we explore the practical impact of kernelization (i.e. data
reduction) on the NP-hard problem of computing the TBR distance between two
unrooted binary phylogenetic trees. This problem is better-known in the
literature as the maximum agreement forest problem, where the goal is to
partition the two trees into a minimum number of common, non-overlapping
subtrees. We have implemented two well-known reduction rules, the subtree and
chain reduction, and five more recent, theoretically stronger reduction rules,
and compare the reduction achieved with and without the stronger rules. We find
that the new rules yield smaller reduced instances and thus have clear
practical added value. In many cases they also cause the TBR distance to
decrease in a controlled fashion. Next, we compare the achieved reduction to
the known worst-case theoretical bounds of 15k-9 and 11k-9 respectively, on the
number of leaves of the two reduced trees, where k is the TBR distance,
observing in both cases a far larger reduction in practice. As a by-product of
our experimental framework we obtain a number of new insights into the actual
computation of TBR distance. We find, for example, that very strong lower
bounds on TBR distance can be obtained efficiently by randomly sampling certain
carefully constructed partitions of the leaf labels, and identify instances
which seem particularly challenging to solve exactly. The reduction rules have
been implemented within our new solver Tubro which combines kernelization with
an Integer Linear Programming (ILP) approach. Tubro also incorporates a number
of additional features, such as a cluster reduction and a practical
upper-bounding heuristic, and it can leverage combinatorial insights emerging
from the proofs of correctness of the reduction rules to simplify the ILP.


  Drug Discovery is a fundamental and ever-evolving field of research. The
design of new candidate molecules requires large amounts of time and money, and
computational methods are being increasingly employed to cut these costs.
Machine learning methods are ideal for the design of large amounts of potential
new candidate molecules, which are naturally represented as graphs. Graph
generation is being revolutionized by deep learning methods, and molecular
generation is one of its most promising applications. In this paper, we
introduce a sequential molecular graph generator based on a set of graph neural
network modules, which we call MG^2N^2. At each step, a node or a group of
nodes is added to the graph, along with its connections. The modular
architecture simplifies the training procedure, also allowing an independent
retraining of a single module. Sequentiality and modularity make the generation
process interpretable. The use of graph neural networks maximizes the
information in input at each generative step, which consists of the subgraph
produced during the previous steps. Experiments of unconditional generation on
the QM9 and Zinc datasets show that our model is capable of generalizing
molecular patterns seen during the training phase, without overfitting. The
results indicate that our method is competitive, and outperforms challenging
baselines for unconditional generation.


  We present an evaluation of the effectiveness of manual contact tracing
compared to bulletin board contact tracing. We show that bulletin board contact
tracing gives comparable results in terms of the reproductive number, duration,
prevalence and incidence but is less resource intensive, easier to implement
and offers a wider range of privacy options. Classical contact tracing focuses
on contacting individuals whom an infectious person has been in proximity to. A
bulletin board approach focuses on identifying locations visited by an
infectious person, and then contacting those who were at those locations. We
present results comparing their effects on the overall reproductive number as
well as the incidence and prevalence of disease. We evaluate them by building a
new discrete time stochastic model based on the Susceptible Exposed Infectious
and Recovered (SEIR) framework for disease spread. We conduct simulation
experiments to quantify the effectiveness of these two models of contact
tracing by calibrating the model to be compatible with SARS-CoV-2. Our
experiments show that location-based bulletin board contact tracing can improve
manual contact tracing.


  Coupled dynamical systems in ecology are known to respond to the seasonal
forcing of their parameters with multiple dynamical behaviours, ranging from
seasonal cycles to chaos. Seasonal forcing is predominantly modelled as a sine
wave but the transition between seasons is often more sudden. Some studies
mentioned the robustness of their results to the shape of the forcing signal,
without detailed analyses. Therefore, whether and how the shape of seasonal
forcing affects the dynamics of coupled dynamical systems remains unclear,
while abrupt seasonal transitions are widespread in ecological systems. To
provide some answers, we conduct a numerical analysis of the dynamical response
of predator--prey communities to the shape of the forcing signal by exploring
the joint effect of two features of seasonal forcing: the magnitude of the
signal, which is classically the only one studied, and the shape of the signal,
abrupt or sinusoidal. We consider both linear and saturating functional
responses, and focus on seasonal forcing of the predator's discovery rate,
which fluctuates with changing environmental conditions and prey's ability to
escape predation. Our numerical results highlight that a more abrupt seasonal
forcing mostly alters the magnitude of population fluctuations and triggers
period--doubling bifurcations, as well as the emergence of chaos, at lower
forcing strength than for sine waves. Controlling the variance of the forcing
signal mitigates this trend but does not fully suppress it, which suggests that
the variance is not the only feature of the shape of seasonal forcing that acts
on community dynamics. Although theoretical studies may predict correctly the
sequence of bifurcations using sine waves as a representation of seasonality,
there is a rationale for applied studies to implement as realistic seasonal
forcing as possible to make precise predictions of community dynamics.


  Explaining why humans cooperate in anonymous contexts is a major goal of
human behavioral ecology, cultural evolution, and related fields. What predicts
cooperation in anonymous contexts is inconsistent across populations, levels of
analysis, and games. For instance, market integration is a key predictor across
ethnolinguistic groups but has inconsistent predictive power at the individual
level. We adapt an idea from 19th-century sociology: people in societies with
greater overlap in ties across domains among community members (Durkheim's
"mechanical" solidarity) will cooperate more with their network partners and
less in anonymous contexts than people in societies with less overlap
("organic" solidarity). This hypothesis, which can be tested at the individual
and community level, assumes that these two types of societies differ in the
importance of keeping existing relationships as opposed to recruiting new
partners. Using multiplex networks, we test this idea by comparing cooperative
tendencies in both anonymous experimental games and real-life communal labor
tasks across 9 Makushi villages in Guyana that vary in the degree of
within-village overlap. Average overlap in a village predicts both real-world
cooperative and anonymous interactions in the predicted direction; individual
overlap also has effects in the expected direction. These results reveal a
consistent patterning of cooperative tendencies at both individual and local
levels and contribute to the debate over the emergence of norms for cooperation
among humans. Multiplex overlap can help us understand inconsistencies in
previous studies of cooperation in anonymous contexts and is an unexplored
dimension with explanatory power at multiple levels of analysis.


  The outer hair cell (OHC) membrane harbors a voltage-dependent protein,
prestin (SLC26a5), in high density, whose charge movement is evidenced as a
nonlinear capacitance (NLC). NLC is bell-shaped, with its peak occurring at a
voltage, Vh, where sensor charge is equally distributed across the plasma
membrane. Thus, Vh provides information on the conformational state of prestin.
Vh is sensitive to membrane tension, shifting to positive voltage as tension
increases and is the basis for considering prestin piezoelectric (PZE). NLC can
be deconstructed into real and imaginary components that report on charge
movements in phase or 90 degrees out of phase with AC voltage. Here we show in
membrane macro-patches of the OHC that there is a partial trade-off in the
magnitude of real and imaginary components as interrogation frequency
increases, as predicted by a recent PZE model (Rabbitt, 2020). However, we find
similar behavior in a simple kinetic model of prestin that lacks piezoelectric
coupling, the meno presto model. At a particular frequency, Fis, the complex
component magnitudes intersect. Using this metric, Fis, which depends on the
frequency response of each complex component, we find that initial Vh
influences Fis; thus, by categorizing patches into groups of different Vh,
(above and below -30 mV) we find that Fis is lower for the negative Vh group.
We also find that the effect of membrane tension on complex NLC is dependent,
but differentially so, on initial Vh. Whereas the negative group exhibits
shifts to higher frequencies for increasing tension, the opposite occurs for
the positive group. Despite complex component trade-offs, the low-pass roll-off
in absolute magnitude of NLC, which varies little with our perturbations and is
indicative of diminishing total charge movement, poses a challenge for a role
of voltage-driven prestin in cochlear amplification at very high frequencies.


  Deterministic compartmental models have been used extensively in modeling
epidemic propagation. These models are required to fit available data and
numerical procedures are often implemented to this end. But not every model
architecture is able to fit the data because the structure of the model imposes
hard constraints on the solutions. We investigate in this work two such
situations: first the distribution of transition times from a compartment to
another may impose a variable number of intermediary states; secondly, a
non-linear relationship between time-dependent measures of compartments sizes
may indicate the need for structurations (i.e., considering several groups of
individuals of heterogeneous characteristics).


  Due to recent breakthroughs in state-of-the-art DNA sequencing technology,
genomics data sets have become ubiquitous. The emergence of large-scale data
sets provides great opportunities for better understanding of genomics,
especially gene regulation. Although each cell in the human body contains the
same set of DNA information, gene expression controls the functions of these
cells by either turning genes on or off, known as gene expression levels. There
are two important factors that control the expression level of each gene: (1)
Gene regulation such as histone modifications can directly regulate gene
expression. (2) Neighboring genes that are functionally related to or interact
with each other that can also affect gene expression level. Previous efforts
have tried to address the former using Attention-based model. However,
addressing the second problem requires the incorporation of all potentially
related gene information into the model. Though modern machine learning and
deep learning models have been able to capture gene expression signals when
applied to moderately sized data, they have struggled to recover the underlying
signals of the data due to the nature of the data's higher dimensionality. To
remedy this issue, we present SimpleChrome, a deep learning model that learns
the latent histone modification representations of genes. The features learned
from the model allow us to better understand the combinatorial effects of
cross-gene interactions and direct gene regulation on the target gene
expression. The results of this paper show outstanding improvements on the
predictive capabilities of downstream models and greatly relaxes the need for a
large data set to learn a robust, generalized neural network. These results
have immediate downstream effects in epigenomics research and drug development.


  The COVID-19 pandemic began over two years ago, yet schools, businesses, and
other organizations are still struggling to keep the risk of disease outbreak
low while returning to (near) normal functionality. Observations from these
past years suggest that this goal can be achieved through the right balance of
mitigation strategies, which may include some combination of mask use,
vaccinations, viral testing, and contact tracing. The choice of mitigation
measures will be uniquely based on the needs and available resources of each
organization. This article presents practical guidance for creating these
policies based on an analytical model of disease spread that captures the
combined effects of each of these interventions. The resulting guidance is
tested through simulation across a wide range of parameters and used to discuss
the spread of disease on college campuses.


  A primary quantity of interest in the study of infectious diseases is the
average number of new infections that an infected person produces. This
so-called reproduction number has significant implications for the disease
progression. There has been increasing literature suggesting that
superspreading, the significant variability in number of new infections caused
by individuals, plays an important role in the spread of SARS-CoV-2. In this
paper, we consider the effect that such superspreading has on the estimation of
the reproduction number and subsequent estimates of future cases. Accordingly,
we employ a simple extension to models currently used in the literature to
estimate the reproduction number and present a case-study of the progression of
COVID-19 in Austria. Our models demonstrate that the estimation uncertainty of
the reproduction number increases with superspreading and that this improves
the performance of prediction intervals. Of independent interest is the
derivation of a transparent formula that connects the extent of superspreading
to the width of credible intervals for the reproduction number. This serves as
a valuable heuristic for understanding the uncertainty surrounding diseases
with superspreading.


  Several implicit methods to infer Horizontal Gene Transfer (HGT) focus on
pairs of genes that have diverged only after the divergence of the two species
in which the genes reside. This situation defines the edge set of a graph, the
later-divergence-time (LDT) graph, whose vertices correspond to genes colored
by their species. We investigate these graphs in the setting of relaxed
scenarios, i.e., evolutionary scenarios that encompass all commonly used
variants of duplication-transfer-loss scenarios in the literature. We
characterize LDT graphs as a subclass of properly vertex-colored cographs, and
provide a polynomial-time recognition algorithm as well as an algorithm to
construct a relaxed scenario that explains a given LDT. An edge in an LDT graph
implies that the two corresponding genes are separated by at least one HGT
event. The converse is not true, however. We show that the complete xenology
relation is described by an rs-Fitch graph, i.e., a complete multipartite graph
satisfying constraints on the vertex coloring. This class of vertex-colored
graphs is also recognizable in polynomial time. We finally address the question
"how much information about all HGT events is contained in LDT graphs" with the
help of simulations of evolutionary scenarios with a wide range of duplication,
loss, and HGT events. In particular, we show that a simple greedy graph editing
scheme can be used to efficiently detect HGT events that are implicitly
contained in LDT graphs.


  Analysis of 1H-NMR spectra is often hindered by large variations that occur
during the collection of these spectra. Large solvent and standard peaks, base
line drift and negative peaks (due to improper phasing) are among some of these
variations. Furthermore, some instrument dependent alterations, such as
incorrect shimming, are also embedded in the recorded spectrum. The
unpredictable nature of these alterations of the signal has rendered the
automated and instrument independent computer analysis of these spectra
unreliable. In this paper, a novel method of extracting the information content
of a signal (in this paper, frequency domain 1H-NMR spectrum), called the
frequency-information transformation (FIT), is presented and compared to a
previously used method (SPUTNIK). FIT can successfully extract the relevant
information to a pattern matching task present in a signal, while discarding
the remainder of a signal by transforming a Fourier transformed signal into an
information spectrum (IS). This technique exhibits the ability of decreasing
the inter-class correlation coefficients while increasing the intra-class
correlation coefficients. Different spectra of the same molecule, in other
words, will resemble more to each other while the spectra of different
molecules will look more different from each other. This feature allows easier
automated identification and analysis of molecules based on their spectral
signatures using computer algorithms.


  Expression level is known to be a strong determinant of a protein's rate of
evolution. But the converse can also be true: evolutionary dynamics can affect
expression levels of proteins. Having implications in both directions fosters
the possibility of a feedback loop, where higher expressed systems are more
likely to improve and be expressed even higher, while those that are expressed
less are eventually lost to drift. Using a minimal model to study this in the
context of a changing environment, we demonstrate that one unexpected
consequence of such a feedback loop is that a slow switch to a new environment
can allow genotypes to reach higher fitness sooner than a direct exposure to
it.


  Coinfection is the process by which a host that is infected with a pathogen
becomes infected by a second pathogen at a later point in time. An
immunosuppressant host response to a primary disease can facilitate spreading
of a subsequent emergent pathogen among the population. Social contact patterns
within the substrate populace can be modelled using complex networks and it has
been shown that contact patterns vastly influence the emergent disease
dynamics. In this paper, we consider the effect of contact clustering on the
coinfection dynamics of two pathogens spreading over a network. We use the
generating function formulation to describe the expected outbreak sizes of each
pathogen and numerically study the threshold criteria that permit the
coexistence of each strain among the network. We find that the effects of
clustering on the levels of coinfection are governed by the details of the
contact topology.


  The reconstruction of possible histories given a sample of genetic data in
the presence of recombination and recurrent mutation is a challenging problem,
but can provide key insights into the evolution of a population. We present
KwARG, which implements a parsimony-based greedy heuristic algorithm for
finding plausible genealogical histories (ancestral recombination graphs) that
are minimal or near-minimal in the number of posited recombination and mutation
events. Given an input dataset of aligned sequences, KwARG outputs a list of
possible candidate solutions, each comprising a list of mutation and
recombination events that could have generated the dataset; the relative
proportion of recombinations and recurrent mutations in a solution can be
controlled via specifying a set of 'cost' parameters. We demonstrate that the
algorithm performs well when compared against existing methods. The software is
made available on GitHub.


  Multi-task learning (MTL) is frequently used in settings where a target task
has to be learnt based on limited training data, but knowledge can be leveraged
from related auxiliary tasks. While MTL can improve task performance overall
relative to single-task learning (STL), these improvements can hide negative
transfer (NT), where STL may deliver better performance for many individual
tasks. Asymmetric multitask feature learning (AMTFL) is an approach that tries
to address this by allowing tasks with higher loss values to have smaller
influence on feature representations for learning other tasks. Task loss values
do not necessarily indicate reliability of models for a specific task. We
present examples of NT in two orthogonal datasets (image recognition and
pharmacogenomics) and tackle this challenge by using aleatoric homoscedastic
uncertainty to capture the relative confidence between tasks, and set weights
for task loss. Our results show that this approach reduces NT providing a new
approach to enable robust MTL.


  The relationship between sequence variation and phenotype is poorly
understood. Here we use metabolomic analysis to elucidate the molecular
mechanism underlying the filamentous phenotype of E. coli strains that carry
destabilizing mutations in the Dihydrofolate Reductase (DHFR). We find that
partial loss of DHFR activity causes SOS response indicative of DNA damage and
cell filamentation. This phenotype is triggered by an imbalance in deoxy
nucleotide levels, most prominently a disproportionate drop in the
intracellular dTTP. We show that a highly cooperative (Hill coefficient 2.5) in
vivo activity of Thymidylate Kinase (Tmk), a downstream enzyme that
phosphorylates dTMP to dTDP, is the cause of suboptimal dTTP levels. dTMP
supplementation in the media rescues filamentation and restores in vivo Tmk
kinetics to almost perfect Michaelis-Menten, like its kinetics in vitro.
Overall, this study highlights the important role of cellular environment in
sculpting enzymatic kinetics with system level implications for bacterial
phenotype.


  Cataract is one of the most prevalent protein aggregation disorders and still
the biggest cause of vision loss worldwide. The human lens, in its core region,
lacks turnover of any cells or cellular components; it has therefore evolved
remarkable mechanisms for resisting protein aggregation for a lifetime. We now
report that one such mechanism relies on an unusually abundant metabolite,
myo-inositol, to suppress light-scattering aggregation of lens proteins. We
quantified aggregation suppression by in vitro turbidimetry and characterized
both macroscopic and microscopic mechanisms of myo-inositol action using
negative-stain electron microscopy, differential scanning fluorometry, and a
thermal scanning Raman spectroscopy apparatus. Given recent metabolomic
evidence that it is dramatically depleted in human cataractous lenses compared
to age-matched controls, we suggest that maintaining or restoring healthy
levels of myo-inositol in the lens may be a simple, safe, and widely available
strategy for reducing the global burden of cataract.


  Protein complexes involved in DNA mismatch repair appear to diffuse along
dsDNA in order to locate a hemimethylated incision site via a dissociative
mechanism. Here, we study the probability that these complexes locate a given
target site via a semi-analytic, Monte Carlo calculation that tracks the
association and dissociation of the complexes. We compare such probabilities to
those obtained using a non-dissociative diffusive scan, and determine that for
experimentally observed diffusion constants, search distances, and search
durations $\textit{in vitro}$, there is neither a significant advantage nor
disadvantage associated with the dissociative mechanism in terms of probability
of successful search, and that both search mechanisms are highly efficient for
a majority of hemimethylated site distances. Furthermore, we examine the space
of physically realistic diffusion constants, hemimethylated site distances, and
association lifetimes and determine the regions in which dissociative searching
is more or less efficient than non-dissociative searching. We conclude that the
dissociative search mechanism is advantageous in the majority of the physically
realistic parameter space.


  Finding ways to overcome the temptation to exploit one another is still a
challenge in behavioural sciences. In the framework of evolutionary game
theory, punishing strategies are frequently used to promote cooperation in
competitive environments. Here, we introduce altruistic punishers in the
spatial public goods game. This strategy acts as a cooperator in the absence of
defectors, otherwise it will punish all defectors in their vicinity while
bearing a cost to do so. We observe three distinct behaviours in our model: i)
in the absence of punishers, cooperators (who don't punish defectors) are
driven to extinction by defectors for most parameter values; ii) clusters of
punishers thrive by sharing the punishment costs when these are low iii) for
higher punishment costs, punishers, when alone, are subject to exploitation but
in the presence of cooperators can form a symbiotic spatial structure that
benefits both. This last observation is our main finding since neither
cooperation nor punishment alone can survive the defector strategy in this
parameter region and the specificity of the symbiotic spatial configuration
shows that lattice topology plays a central role in sustaining cooperation.
Results were obtained by means of Monte Carlo simulations on a square lattice
and subsequently confirmed by a pairwise comparison of different strategies'
payoffs in diverse group compositions, leading to a phase diagram of the
possible states.


  When studying the dynamics of trait distribution of populations in a
heterogeneous environment, classical models from quantitative genetics choose
to look at its system of moments, specifically the first two ones.
Additionally, in order to close the resulting system of equations, they often
assume that the trait distribution is Gaussian (see for instance Ronce and
Kirkpatrick 2001). The aim of this paper is to introduce a mathematical
framework that follows the whole trait distribution (without prior assumption)
to study evolutionary dynamics of sexually reproducing populations.
Specifically, it focuses on complex traits, whose inheritance can be encoded by
the infinitesimal model of segregation (Fisher 1919). We show that it allows us
to derive a regime in which our model gives the same dynamics as when assuming
a Gaussian trait distribution. To support that, we compare the stationary
problems of the system of moments derived from our model with the one given in
Ronce and Kirkpatrick 2001 and show that they are equivalent under this regime
and do not need to be otherwise. Moreover, under this regime of equivalence, we
show that a separation bewteen ecological and evolutionary time scales arises.
A fast relaxation toward monomorphism allows us to reduce the complexity of the
system of moments, using a slow-fast analysis. This reduction leads us to
complete, still in this regime, the analytical description of the bistable
asymmetrical equilibria numerically found in Ronce and Kirkpatrick 2001. More
globally, we provide explicit modelling hypotheses that allow for such local
adaptation patterns to occur.


  The recent availability of huge, many-dimensional data sets, like those
arising from genome-wide association studies (GWAS), provides many
opportunities for strengthening causal inference. One popular approach is to
utilize these many-dimensional measurements as instrumental variables
(instruments) for improving the causal effect estimate between other pairs of
variables. Unfortunately, searching for proper instruments in a
many-dimensional set of candidates is a daunting task due to the intractable
model space and the fact that we cannot directly test which of these candidates
are valid, so most existing search methods either rely on overly stringent
modeling assumptions or fail to capture the inherent model uncertainty in the
selection process. We show that, as long as at least some of the candidates are
(close to) valid, without knowing a priori which ones, they collectively still
pose enough restrictions on the target interaction to obtain a reliable causal
effect estimate. We propose a general and efficient causal inference algorithm
that accounts for model uncertainty by performing Bayesian model averaging over
the most promising many-dimensional instrumental variable models, while at the
same time employing weaker assumptions regarding the data generating process.
We showcase the efficiency, robustness and predictive performance of our
algorithm through experimental results on both simulated and real-world data.


  The use of genetic variants as instrumental variables - an approach known as
Mendelian randomization - is a popular epidemiological method for estimating
the causal effect of an exposure (phenotype, biomarker, risk factor) on a
disease or health-related outcome from observational data. Instrumental
variables must satisfy strong, often untestable assumptions, which means that
finding good genetic instruments among a large list of potential candidates is
challenging. This difficulty is compounded by the fact that many genetic
variants influence more than one phenotype through different causal pathways, a
phenomenon called horizontal pleiotropy. This leads to errors not only in
estimating the magnitude of the causal effect but also in inferring the
direction of the putative causal link. In this paper, we propose a Bayesian
approach called BayesMR that is a generalization of the Mendelian randomization
technique in which we allow for pleiotropic effects and, crucially, for the
possibility of reverse causation. The output of the method is a posterior
distribution over the target causal effect, which provides an immediate and
easily interpretable measure of the uncertainty in the estimation. More
importantly, we use Bayesian model averaging to determine how much more likely
the inferred direction is relative to the reverse direction.


  One of the most common mutations in the serine protease inhibitor Kazal type
1 (SPINK1) gene is the N34S variant which is strongly associated with chronic
pancreatitis. Although it is assumed that N34S mutation constitutes a high-risk
factor, the underlying pathologic mechanism is still unknown. In the present
study, we investigated the impact of physiological stress factors on SPINK1
protein structure and trypsin inhibitor function using biophysical methods. Our
circular dichroism spectroscopy data revealed differences in the secondary
structure of SPINK1 and N34S mutant suggesting protein structural changes
induced by the mutation as an impairment that could be disease-relevant. We
further confirmed that both SPINK1 (KD of 0.15 +/- 0.06 nM) and its N34S
variant (KD of 0.08 +/- 0.02 nM) have similar binding affinity and inhibitory
effect towards trypsin as shown by surface plasmon resonance and trypsin
inhibition assay studies, respectively. We found that stress conditions such as
altered ion concentrations (i.e. potassium, calcium), temperature shifts, as
well as environmental pH lead to insignificant differences in trypsin
inhibition between SPINK1 and N34S mutant. However, we have shown that the
environmental pH induces structural changes in both SPINK1 constructs in a
different manner. Our findings suggest protein structural changes in the N34S
variant as an impairment of SPINK1 and environmental pH shift as a trigger that
could play a role in disease progression of pancreatitis.


  So far most of the analysis of coronavirus 2020 epidemic data has been
focusing on a short-time window and consequently a quantitative test of
statistical physical laws of Coronavirus Epidemics with Containment Measures
(CEwCM) is currently lacking. Here we report a quantitative analysis of CEwCM
over 230 days, covering the full-time lapse of the first epidemic wave. We use
a 3D phase diagram tracking the simultaneous evolution of the doubling time
Td(t) and reproductive number Rt(t) showing that this expanded parameter space
is needed for biological physics of CEwCP. We have verified that in the
supercritical [Rt(t)>1, Td(t)<40 days] regime i) the curve Z(t) of total
infected cases follows the growth rate called Ostwald law; ii) the doubling
time follows the exponential law Td(t)=A exp((t-t0)/s) as a function of time
and iii) the power law Td(t)=C(Rt(t)-1)^-n is verified with the exponent n
depending on the definition of Rt(t). The log-log plots Td(t) versus (Rt-1) of
the second 2020 epidemic wave unveil in the subcritical regime [Td(t)>100 days]
arrested metastable phases with Rt>1 where Td(t) was kept constant followed by
its explosion and its containment following the same power law as in the first
wave


  Background
  SAMtools and BCFtools are widely used programs for processing and analysing
high-throughput sequencing data.
  Findings
  The first version appeared online twelve years ago and has been maintained
and further developed ever since, with many new features and improvements added
over the years. The SAMtools and BCFtools packages represent a unique
collection of tools that have been used in numerous other software projects and
countless genomic pipelines.
  Conclusion
  Both SAMtools and BCFtools are freely available on GitHub under the
permissive MIT licence, free for both non-commercial and commercial use. Both
packages have been installed over a million times via Bioconda. The source code
and documentation are available from http://www.htslib.org.


  Vivid structural colors in birds are a conspicuous and vital part of their
phenotype. They are produced by a rich diversity of integumentary photonic
nanostructures in skin and feathers. Unlike pigmentary coloration, whose
molecular genetic basis is being elucidated, little is known regarding the
pathways underpinning organismal structural coloration. Here, we review
available data on the development of avian structural colors. In particular,
feather photonic nanostructures are understood to be intracellularly
self-assembled by physicochemical forces typically seen in soft colloidal
systems. We identify promising avenues for future research that can address
current knowledge gaps, which is also highly relevant for the sustainable
engineering of advanced bioinspired and biomimetic materials.


  Sensitivity studies the network response to perturbations. We consider local
perturbations of the concentrations of metabolites at an equilibrium. We
investigate the responses in the network, both of the metabolite concentrations
and of the reaction fluxes. Our approach is purely qualitative, rather than
quantitative. In fact, our analysis is based, solely, on the stoichiometry of
the reaction network and we do not require any quantitative information on the
reaction rates. Instead, the description is done only in algebraic terms, and
the only data required is the network structure. Biological applications
include metabolic control and network reconstruction.


  We study the correlation between the codon usage bias of genetic sequences
and the network features of protein-protein interaction (PPI) in bacterial
species. We use PCA techniques in the space of codon bias indices to show that
genes with similar patterns of codon usage have a significantly higher
probability that their encoded proteins are functionally connected and
interacting. Importantly, this signal emerges when multiple aspects of codon
bias are taken into account at the same time. The present study extends our
previous observations on E.Coli over a wide set of 34 bacteria. These findings
could allow for future investigations on the possible effects of codon bias on
the topology of the PPI network, with the aim of improving existing
bioinformatics methods for predicting protein interactions.


  The estrogen receptor is a nuclear hormone receptor activated by the natural
steroid hormone 17$\beta$-estradiol (E2). Fragment molecular orbital (FMO)
calculations were performed which allowed us to obtain the interaction energy
($E_{int}$) between E2, 17$\alpha$-estradiol (17$\alpha$-E2) and the human
estrogen receptor $\alpha$ ligand-binding domain. In aqueous media the
MP2/6-31G(d) $E_{int}$ was of -88.52 kcal/mol for E2 and -78.73 kcal/mol for
17$\alpha$-E2. Attractive dispersion interactions were observed between ligands
and all surrounding hydrophobic residues. Water molecules were found at the
binding site and strong attractive electrostatic interactions were observed
between the ligands and the Glu 353 and His 524 residues. The essential
dynamics revealed that E2 adapts to the binding site and its motion, in a
sense, synchronizes with the whole receptor; while 17$\alpha$-E2, with its
motion of greater amplitude compared to E2, disturbs the binding site. Perhaps
this feature of the normal substrate is a necessary condition for biological
function. Another important requirement relates to the number of water
molecules at the binding site. Therefore, negative values in $E_{int}$ is a
necessary but not sufficient condition since, it is also necessary to consider
the conformers population that fulfill all the requirements that ensure a
biological response.


  Widespread growth signatures in COVID-19 confirmed case counts are reported,
with sharp transitions between three distinct dynamical regimes (exponential,
superlinear and sublinear). Through analytical and numerical analysis, a novel
framework is developed that exploits information in these signatures. An
approach well known to physics is applied, where one looks for common dynamical
features, independently from differences in other factors. These features and
associated scaling laws are used as a powerful tool to pinpoint regions where
analytical derivations are effective, get an insight into qualitative changes
of the disease progression, and infer the key infection parameters. The
developed framework for joint analytical and numerical analysis of empirically
observed COVID-19 growth patterns can lead to a fundamental understanding of
infection progression under strong control measures, applicable to outbursts of
both COVID-19 and other infectious diseases.


  The spatial scale of population synchrony gives the characteristic distance
at which the population fluctuations are correlated. Therefore, it gives also
the characteristic size of the regions of simultaneous population depletion, or
even extinction. Single-species previous results imply that the spatial scale
of population synchrony is equal or greater (due to dispersion) than the
spatial scale of synchrony of environmental fluctuations. Theoretical results
on multispecies ecosystems points that interspecies interactions modify the
spatial scale of population synchrony. In particular, recent results on two
species ecosystems, for two competitors and for predator-prey, point that the
spatial scale of population synchrony generally increases as the fluctuations
propagates through the food web, i.e., the species more directly affected by
environmental fluctuations presents the smaller spatial scale of population
synchrony. Here, we found that this behaviour is generally true for a two
species ecosystem. The exception to this behaviour are the particular cases
where the population fluctuations of one of the species does not damp by its
own, but requires a strong transfer of the fluctuation to the other species to
be damped. These analytical results illustrate the importance of applying an
ecosystem rather than a single-species perspective when developing sustainable
harvestings or assessing the extinction risk of endangered species.


  How to produce expressive molecular representations is a fundamental
challenge in AI-driven drug discovery. Graph neural network (GNN) has emerged
as a powerful technique for modeling molecular data. However, previous
supervised approaches usually suffer from the scarcity of labeled data and have
poor generalization capability. Here, we proposed a novel Molecular
Pre-training Graph-based deep learning framework, named MPG, that leans
molecular representations from large-scale unlabeled molecules. In MPG, we
proposed a powerful MolGNet model and an effective self-supervised strategy for
pre-training the model at both the node and graph-level. After pre-training on
11 million unlabeled molecules, we revealed that MolGNet can capture valuable
chemistry insights to produce interpretable representation. The pre-trained
MolGNet can be fine-tuned with just one additional output layer to create
state-of-the-art models for a wide range of drug discovery tasks, including
molecular properties prediction, drug-drug interaction, and drug-target
interaction, involving 13 benchmark datasets. Our work demonstrates that MPG is
promising to become a novel approach in the drug discovery pipeline.


  Supergenes are genomic regions containing sets of tightly linked loci that
control multi-trait phenotypic polymorphisms under balancing selection. Recent
advances in genomics have uncovered significant variation in both the genomic
architecture as well as the mode of origin of supergenes across diverse
organismal systems. Although the role of genomic architecture for the origin of
supergenes has been much discussed, differences in the genomic architecture
also subsequently affect the evolutionary trajectory of supergenes and the rate
of degeneration of supergene haplotypes. In this review, we synthesize recent
genomic work and historical models of supergene evolution, highlighting how the
genomic architecture of supergenes affects their evolutionary fate. We discuss
how recent findings on classic supergenes involved in governing ant colony
social form, mimicry in butterflies, and heterostyly in flowering plants relate
to theoretical expectations. Furthermore, we use forward simulations to
demonstrate that differences in genomic architecture affect the degeneration of
supergenes. Finally, we discuss implications of the evolution of supergene
haplotypes for the long-term fate of balanced polymorphisms governed by
supergenes.


  When designing new molecules with particular properties, it is not only
important what to make but crucially how to make it. These instructions form a
synthesis directed acyclic graph (DAG), describing how a large vocabulary of
simple building blocks can be recursively combined through chemical reactions
to create more complicated molecules of interest. In contrast, many current
deep generative models for molecules ignore synthesizability. We therefore
propose a deep generative model that better represents the real world process,
by directly outputting molecule synthesis DAGs. We argue that this provides
sensible inductive biases, ensuring that our model searches over the same
chemical space that chemists would also have access to, as well as
interpretability. We show that our approach is able to model chemical space
well, producing a wide range of diverse molecules, and allows for unconstrained
optimization of an inherently constrained problem: maximize certain chemical
properties such that discovered molecules are synthesizable.


  We present a unified framework for modelling genomes and their rearrangements
in a genome algebra, as elements that simultaneously incorporate all physical
symmetries. Building on previous work utilising the group algebra of the
symmetric group, we explicitly construct the genome algebra for the case of
unsigned circular genomes with dihedral symmetry and show that the maximum
likelihood estimate (MLE) of genome rearrangement distance can be validly and
more efficiently performed in this setting. We then construct the genome
algebra for a more general case, that is, for genomes that may be represented
by elements of an arbitrary group and symmetry group, and show that the MLE
computations can be performed entirely within this framework. There is no
prescribed model in this framework; that is, it allows any choice of
rearrangements that preserve the set of regions, along with arbitrary weights.
Further, since the likelihood function is built from path probabilities -- a
generalisation of path counts -- the framework may be utilised for any distance
measure that is based on path probabilities.


  Evolution gave rise to creatures that are arguably more sophisticated than
the greatest human-designed systems. This feat has inspired computer scientists
since the advent of computing and led to optimization tools that can evolve
complex neural networks for machines -- an approach known as "neuroevolution".
After a few successes in designing evolvable representations for
high-dimensional artifacts, the field has been recently revitalized by going
beyond optimization: to many, the wonder of evolution is less in the perfect
optimization of each species than in the creativity of such a simple iterative
process, that is, in the diversity of species. This modern view of artificial
evolution is moving the field away from microevolution, following a fitness
gradient in a niche, to macroevolution, filling many niches with highly
different species. It already opened promising applications, like evolving gait
repertoires, video game levels for different tastes, and diverse designs for
aerodynamic bikes.


  One of the last chapters in the long course of human evolution was the shift
from hunting and gathering to the production of food or strategies of
subsistence based on farming and the herding of animals. In Southwest Asia, the
first steps towards the origins of agriculture began some 12,000 years ago and
then spread over most regions of Europe during the span of time from about
10,400 years ago (the start of the so-called PPNB on the island of Cyprus)
through around 6,000 years ago. The aim of this chapter is to provide an
overview on the research that we have done on the question of the Neolithic
transition in Europe, which began when Luca CavalliSforza, a leading figure in
the field of human population genetics, and I began to work in collaboration at
the University of Pavia in November of 1970. This draft forms the basis of my
paper as part of the Festschrift prepared for the 45th anniversary of Ryszard
Grygiel and Peter Bogucki's scientific cooperation.


  Graphs as a type of data structure have recently attracted significant
attention. Representation learning of geometric graphs has achieved great
success in many fields including molecular, social, and financial networks. It
is natural to present proteins as graphs in which nodes represent the residues
and edges represent the pairwise interactions between residues. However, 3D
protein structures have rarely been studied as graphs directly. The challenges
include: 1) Proteins are complex macromolecules composed of thousands of atoms
making them much harder to model than micro-molecules. 2) Capturing the
long-range pairwise relations for protein structure modeling remains
under-explored. 3) Few studies have focused on learning the different
attributes of proteins together. To address the above challenges, we propose a
new graph neural network architecture to represent the proteins as 3D graphs
and predict both distance geometric graph representation and dihedral geometric
graph representation together. This gives a significant advantage because this
network opens a new path from the sequence to structure. We conducted extensive
experiments on four different datasets and demonstrated the effectiveness of
the proposed method.


  Background: Coronavirus disease 2019 (COVID-19) and Influenza A are common
disease caused by viral infection. The clinical symptoms and transmission
routes of the two diseases are similar. However, there are no relevant studies
on laboratory diagnostic models to discriminate COVID-19 and influenza A. This
study aims at establishing a signature of laboratory findings to tell patients
with COVID-19 apart from those with influenza A perfectly. Materials: In this
study, 56 COVID-19 patients and 54 influenza A patients were included.
Laboratory findings, epidemiological characteristics and demographic data were
obtained from electronic medical record databases. Elastic network models,
followed by a stepwise logistic regression model were implemented to identify
indicators capable of discriminating COVID-19 and influenza A. A nomogram is
diagramed to show the resulting discriminative model. Results: The majority of
hematological and biochemical parameters in COVID-19 patients were
significantly different from those in influenza A patients. In the final model,
albumin/globulin (A/G), total bilirubin (TBIL) and erythrocyte specific volume
(HCT) were selected as predictors. Using an external dataset, the model was
validated to perform well. Conclusion: A diagnostic model of laboratory
findings was established, in which A/G, TBIL and HCT were included as highly
relevant indicators for the segmentation of COVID-19 and influenza A, providing
a complimentary means for the precise diagnosis of these two diseases.


  We employ the epidemic Renormalization Group (eRG) framework to understand,
reproduce and predict the COVID-19 pandemic diffusion across the US. The human
mobility across different geographical US divisions is modelled via open source
flight data alongside the impact of social distancing for each such division.
We analyse the impact of the vaccination strategy on the current pandemic wave
dynamics in the US. We observe that the ongoing vaccination campaign will not
impact the current pandemic wave and therefore strict social distancing
measures must still be enacted. To curb the current and the next waves our
results indisputably show that vaccinations alone are not enough and strict
social distancing measures are required until sufficient immunity is achieved.
Our results are essential for a successful vaccination strategy in the US.


  We report three pre-registered studies (total N=1,799) exploring the effect
of nudging personal and injunctive norms in decisions that involve a trade-off
between objective equality and efficiency. The first two studies provide
evidence that: (i) nudging the personal norm has a similar effect to nudging
the injunctive norm; (ii) when both norms are nudged towards the same
direction, there is no additive effect; (iii) when the personal norm and the
injunctive norm are nudged towards opposite directions, some people tend to
follow the personal norm, while others tend to follow the injunctive norm.
Study 3 tests whether these two classes of people, those who tend to follow the
injunctive norm and those who tend to follow the personal norm, map onto the
two sub-dimensions of Aquino and Reed's moral identity scale. We find partial
evidence of this hypothesis: people higher in the symbolisation dimension are
weakly more likely to follow the injunctive norm; however, we do not find any
evidence that people higher in the internalisation dimension are more likely to
follow the personal norm.


  Branching morphogenesis generates epithelial trees which facilitate gas
exchange, filtering, as well as secretion processes with their large surface to
volume ratio. In this review, we focus on the developmental mechanisms that
control the early stages of lung branching morphogenesis. Lung branching
morphogenesis involves the stereotypic, recurrent definition of new branch
points, subsequent epithelial budding, and lung tube elongation. We discuss
current models and experimental evidence for each of these steps. Finally, we
discuss the role of the mesenchyme in determining the organ-specific shape.


  Faced with the 2020 SARS-CoV2 epidemic, public health officials have been
seeking models that could be used to predict not only the number of new cases
but also the levels of hospitalisation, critical care and deaths. In this paper
we present a stochastic compartmental model capable of real-time monitoring and
forecasting of the pandemic incorporating multiple streams of real-world data,
reported cases, testing intensity, deaths, hospitalisations and critical care
occupancy. Model parameters are estimated via a Bayesian particle filtering
technique. The model successfully tracks the key variables (reported cases,
critical care and deaths) throughout the two waves (March-June and
September-November 2020) of the COVID-19 outbreak in Scotland. The model
hospitalisation predictions in Summer 2020 are consistently lower than the
recorded data, but consistent with the change to the reporting criteria by the
Health Protection Scotland on 15th September. Most parameter estimates were
constant over the two waves, but the infection rate and consequently the
reproductive number decrease in the later stages of the first wave and increase
again from July 2020. The death rates are initially high but decrease over
Summer 2020 before rising again in November. The model can also be used to
provide short-term predictions. We show that the 2-week predictability is very
good for the period from March to June 2020, even at early stages of the
pandemic. The model has been slower to pick up the increase in the case numbers
in September 2020 but forecasting improves again in the later stages of the
epidemic.


  The infection fatality rate (IFR) of COVID-19 is one of the measures of
disease impact that can be of importance for policy making. Here we show that
many of the studies on which these estimates are based are scientifically
flawed for reasons which include: nonsensical equations, unjustified
assumptions, small sample sizes, non-representative sampling (systematic
biases), incorrect definitions of symptomatic and asymptomatic cases
(identified and unidentified cases), typically assuming that cases which are
asymptomatic at the time of testing are the same as completely asymptomatic
(never symptomatic) cases. Moreover, a widely cited meta-analysis misrepresents
some of the IFR values in the original studies, and makes inappropriate
duplicate use of studies, or the information from studies, so that the results
that are averaged are not independent from each other. The lack of validity of
these research papers is of particular importance in view of their influence on
policies that affect lives and well-being in confronting a worldwide pandemic.


  There is an urgent and well-recognized need to extend genetic studies to
diverse populations, but several obstacles continue to be prohibitive,
including (but not limited to) the difficulty of recruiting individuals from
diverse populations in large numbers and the lack of representation in
available genomic references. These obstacles notwithstanding, studying
multiple diverse populations would provide informative, population-specific
insights. Using Native Hawaiians as an example of an understudied population
with a unique evolutionary history, I will argue that by developing key genomic
resources and integrating evolutionary thinking into genetic epidemiology, we
will have the opportunity to efficiently advance our knowledge of the genetic
risk factors, ameliorate health disparity, and improve healthcare in this
underserved population.


  Many of the building blocks of life such as amino acids and nucleotides are
chiral, i.e., different from their mirror image. Contemporary life selects and
synthesizes only one of two possible handednesses. In an abiotic environment,
however, there are usually equally many left- and right-handed molecules. If
homochirality was a prerequisite of life, there must have been physical or
chemical circumstances that led to the selection of a certain preference.
Conversely, if it was a consequence of life, we must identify possible pathways
for accomplishing a transition from a racemic to a homochiral chemistry. After
a discussion of the observational evidence, I will review ideas where
homochirality of any handedness could emerge as a consequence of the first
polymerization events of nucleotides in an emerging RNA world. These mechanisms
are not limited to nucleotides, but can also occur for peptides, as a precursor
to the RNA world. The question of homochirality is, in this sense, intimately
tied to the origin of life. Future Mars missions may be able to detect
biomolecules of extant or extinct life. We will therefore also discuss possible
experimental setups for determining the chirality of primitive life forms in
situ on Mars.