Limits to selection on standing variation in an asexual population.

We consider how a population of N haploid individuals responds to directional selection on standing variation, with no new variation from recombination or mutation. Individuals have trait values z1,…,zN, which are drawn from a distribution ψ; the fitness of individual i is proportional to [Formula: see text] . For illustration, we consider the Laplace and Gaussian distributions, which are parametrised only by the variance V0, and show that for large N, there is a scaling limit which depends on a single parameter NV0. When selection is weak relative to drift (NV0≪1), the variance decreases exponentially at rate 1/N, and the expected ultimate gain in log fitness (scaled by V0), is just NV0, which is the same as Robertson's (1960) prediction for a sexual population. In contrast, when selection is strong relative to drift (NV0≫1), the ultimate gain can be found by approximating the establishment of alleles by a branching process in which each allele competes independently with the population mean and the fittest allele to establish is certain to fix. Then, if the probability of survival to time t∼1/V0 of an allele with value z is P(z), with mean P¯, the winning allele is the fittest of NP¯ survivors drawn from a distribution ψP/P¯. The expected ultimate change is ∼2log(1.15NV0) for a Gaussian distribution, and ∼-12log0.36NV0-log-log0.36NV0 for a Laplace distribution. This approach also predicts the variability of the process, and its dynamics; we show that in the strong selection regime, the expected genetic variance decreases as ∼t-3 at large times. We discuss how these results may be related to selection on standing variation that is spread along a linear chromosome.

What is reproductive isolation?

Reproductive isolation (RI) is a core concept in evolutionary biology. It has been the central focus of speciation research since the modern synthesis and is the basis by which biological species are defined. Despite this, the term is used in seemingly different ways, and attempts to quantify RI have used very different approaches. After showing that the field lacks a clear definition of the term, we attempt to clarify key issues, including what RI is, how it can be quantified in principle, and how it can be measured in practice. Following other definitions with a genetic focus, we propose that RI is a quantitative measure of the effect that genetic differences between populations have on gene flow. Specifically, RI compares the flow of neutral alleles in the presence of these genetic differences to the flow without any such differences. RI is thus greater than zero when genetic differences between populations reduce the flow of neutral alleles between populations. We show how RI can be quantified in a range of scenarios. A key conclusion is that RI depends strongly on circumstances-including the spatial, temporal and genomic context-making it difficult to compare across systems. After reviewing methods for estimating RI from data, we conclude that it is difficult to measure in practice. We discuss our findings in light of the goals of speciation research and encourage the use of methods for estimating RI that integrate organismal and genetic approaches.

Accumulation and maintenance of information in evolution.

Selection accumulates information in the genome-it guides stochastically evolving populations toward states (genotype frequencies) that would be unlikely under neutrality. This can be quantified as the Kullback-Leibler (KL) divergence between the actual distribution of genotype frequencies and the corresponding neutral distribution. First, we show that this population-level information sets an upper bound on the information at the level of genotype and phenotype, limiting how precisely they can be specified by selection. Next, we study how the accumulation and maintenance of information is limited by the cost of selection, measured as the genetic load or the relative fitness variance, both of which we connect to the control-theoretic KL cost of control. The information accumulation rate is upper bounded by the population size times the cost of selection. This bound is very general, and applies across models (Wright-Fisher, Moran, diffusion) and to arbitrary forms of selection, mutation, and recombination. Finally, the cost of maintaining information depends on how it is encoded: Specifying a single allele out of two is expensive, but one bit encoded among many weakly specified loci (as in a polygenic trait) is cheap.

Inversions and parallel evolution.

Local adaptation leads to differences between populations within a species. In many systems, similar environmental contrasts occur repeatedly, sometimes driving parallel phenotypic evolution. Understanding the genomic basis of local adaptation and parallel evolution is a major goal of evolutionary genomics. It is now known that by preventing the break-up of favourable combinations of alleles across multiple loci, genetic architectures that reduce recombination, like chromosomal inversions, can make an important contribution to local adaptation. However, little is known about whether inversions also contribute disproportionately to parallel evolution. Our aim here is to highlight this knowledge gap, to showcase existing studies, and to illustrate the differences between genomic architectures with and without inversions using simple models. We predict that by generating stronger effective selection, inversions can sometimes speed up the parallel adaptive process or enable parallel adaptation where it would be impossible otherwise, but this is highly dependent on the spatial setting. We highlight that further empirical work is needed, in particular to cover a broader taxonomic range and to understand the relative importance of inversions compared to genomic regions without inversions. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

The response of a metapopulation to a changing environment.

A species distributed across diverse environments may adapt to local conditions. We ask how quickly such a species changes its range in response to changed conditions. Szép et al. (Szép E, Sachdeva H, Barton NH. 2021 Polygenic local adaptation in metapopulations: a stochastic eco-evolutionary model. Evolution75, 1030-1045 (doi:10.1111/evo.14210)) used the infinite island model to find the stationary distribution of allele frequencies and deme sizes. We extend this to find how a metapopulation responds to changes in carrying capacity, selection strength, or migration rate when deme sizes are fixed. We further develop a 'fixed-state' approximation. Under this approximation, polymorphism is only possible for a narrow range of habitat proportions when selection is weak compared to drift, but for a much wider range otherwise. When rates of selection or migration relative to drift change in a single deme of the metapopulation, the population takes a time of order m-1 to reach the new equilibrium. However, even with many loci, there can be substantial fluctuations in net adaptation, because at each locus, alleles randomly get lost or fixed. Thus, in a finite metapopulation, variation may gradually be lost by chance, even if it would persist in an infinite metapopulation. When conditions change across the whole metapopulation, there can be rapid change, which is predicted well by the fixed-state approximation. This work helps towards an understanding of how metapopulations extend their range across diverse environments. This article is part of the theme issue 'Species' ranges in the face of changing environments (Part II)'.

Genetic load and extinction in peripheral populations: the roles of migration, drift and demographic stochasticity.

We analyse how migration from a large mainland influences genetic load and population numbers on an island, in a scenario where fitness-affecting variants are unconditionally deleterious, and where numbers decline with increasing load. Our analysis shows that migration can have qualitatively different effects, depending on the total mutation target and fitness effects of deleterious variants. In particular, we find that populations exhibit a genetic Allee effect across a wide range of parameter combinations, when variants are partially recessive, cycling between low-load (large-population) and high-load (sink) states. Increased migration reduces load in the sink state (by increasing heterozygosity) but further inflates load in the large-population state (by hindering purging). We identify various critical parameter thresholds at which one or other stable state collapses, and discuss how these thresholds are influenced by the genetic versus demographic effects of migration. Our analysis is based on a 'semi-deterministic' analysis, which accounts for genetic drift but neglects demographic stochasticity. We also compare against simulations which account for both demographic stochasticity and drift. Our results clarify the importance of gene flow as a key determinant of extinction risk in peripheral populations, even in the absence of ecological gradients. This article is part of the theme issue 'Species' ranges in the face of changing environments (part I)'.

Discovery of GSK251: A Highly Potent, Highly Selective, Orally Bioavailable Inhibitor of PI3Kδ with a Novel Binding Mode.

Optimization of a previously reported lead series of PI3Kδ inhibitors with a novel binding mode led to the identification of a clinical candidate compound 31 (GSK251). Removal of an embedded Ames-positive heteroaromatic amine by reversing a sulfonamide followed by locating an interaction with Trp760 led to a highly selective compound 9. Further optimization to avoid glutathione trapping, to enhance potency and selectivity, and to optimize an oral pharmacokinetic profile led to the discovery of compound 31 (GSK215) that had a low predicted daily dose (45 mg, b.i.d) and a rat toxicity profile suitable for further development.

Design and Development of a Macrocyclic Series Targeting Phosphoinositide 3-Kinase δ.

A macrocyclization approach has been explored on a series of benzoxazine phosphoinositide 3-kinase δ inhibitors, resulting in compounds with improved potency, permeability, and in vivo clearance while maintaining good solubility. The thermodynamics of binding was explored via surface plasmon resonance, and the binding of lead macrocycle 19 was found to be almost exclusively entropically driven compared with progenitor 18, which demonstrated both enthalpic and entropic contributions. The pharmacokinetics of macrocycle 19 was also explored in vivo, where it showed reduced clearance when compared with the progenitor 18. This work adds to the growing body of evidence that macrocyclization could provide an alternative and complementary approach to the design of small-molecule inhibitors, with the potential to deliver differentiated properties.

Rapid identification of highly potent human anti-GPCR antagonist monoclonal antibodies.

Complex cellular targets such as G protein-coupled receptors (GPCRs), ion channels, and other multi-transmembrane proteins represent a significant challenge for therapeutic antibody discovery, primarily because of poor stability of the target protein upon extraction from cell membranes. To assess whether a limited set of membrane-bound antigen formats could be exploited to identify functional antibodies directed against such targets, we selected a GPCR of therapeutic relevance (CCR1) and identified target binders using an in vitro yeast-based antibody discovery platform (AdimabTM) to expedite hit identification. Initially, we compared two different biotinylated antigen formats overexpressing human CCR1 in a 'scouting' approach using a subset of the antibody library. Binders were isolated using streptavidin-coated beads, expressed as yeast supernatants, and screened using a high-throughput binding assay and flow cytometry on appropriate cell lines. The most suitable antigen was then selected to isolate target binders using the full library diversity. This approach identified a combined total of 183 mAbs with diverse heavy chain sequences. A subset of clones exhibited high potencies in primary cell chemotaxis assays, with IC50 values in the low nM/high pM range. To assess the feasibility of any further affinity enhancement, full-length hCCR1 protein was purified, complementary-determining region diversified libraries were constructed from a high and lower affinity mAb, and improved binders were isolated by fluorescence-activated cell sorting selections. A significant affinity enhancement was observed for the lower affinity parental mAb, but not the high affinity mAb. These data exemplify a methodology to generate potent human mAbs for challenging targets rapidly using whole cells as antigen and define a route to the identification of affinity-matured variants if required.

Bone retouchers and technological continuity in the Middle Stone Age of North Africa.

Evidence for specialised bone tools has recently been reported for the Middle Stone Age of North Africa [one], which complements similar finds of slightly younger age in South Africa [two, three]. However, until now scant reference has been made to lesser known tools also made of bone ('bone retouchers') that were employed specifically as intermediaries for working or refining stone artefacts, that are sometimes present in these assemblages. In this paper we describe 20 bone retouchers from the cave of Grotte des Pigeons at Taforalt in north-east Morocco. This is the largest stratified assemblage of bone retouchers from a North African MSA site, and the biggest single collection so far from the African Continent. A total of 18 bone retouchers was recovered in securely dated archaeological levels spanning a period from ~ 84.5 ka to 24 ka cal BP. A further two bone retouchers were found in a layer at the base of the deposits in association with Aterian artefacts dating to around 85,000 BP and so far represent the earliest evidence of this type of tool at Taforalt. In this paper we present a first, detailed description of the finds and trace the stages of their production, use and discard (chaîne opératoire). At the same time, we assess if there were diachronic changes in their form and function and, finally, explore their presence in relation to stone tools from the same occupation layers of the cave.

Optimization of Orally Bioavailable PI3Kδ Inhibitors and Identification of Vps34 as a Key Selectivity Target.

Optimization of a lead series of PI3Kδ inhibitors based on a dihydroisobenzofuran core led to the identification of potent, orally bioavailable compound 19. Selectivity profiling of compound 19 showed similar potency for class III PI3K, Vps34, and PI3Kδ, and compound 19 was not well-tolerated in a 7-day rat toxicity study. Structure-based design led to an improvement in selectivity for PI3Kδ over Vps34 and, a focus on oral phramacokinetics properties resulted in the discovery of compound 41, which showed improved toxicological outcomes at similar exposure levels to compound 19.

Free-Wilson Analysis of Comprehensive Data on Phosphoinositide-3-kinase (PI3K) Inhibitors Reveals Importance of N-Methylation for PI3Kδ Activity.

Phosphoinositide-3-kinase δ (PI3Kδ) is a critical regulator of cell growth and transformation and has been explored as a therapeutic target for a range of diseases. Through the exploration of the thienopyrimidine scaffold, we have identified a ligand-efficient methylation that leads to remarkable selectivity for PI3Kδ over the closely related isoforms. Interrogation through the Free-Wilson analysis highlights the innate selectivity the thienopyrimidine scaffold has for PI3Kδ and provides a predictive model for the activity against the PI3K isoforms.

Infant funerary behavior and kinship in Pleistocene hunter-gatherers from Morocco.

Infant remains are relatively uncommon in the late Pleistocene (Upper Palaeolithic) archaeological record. Funerary treatment is considered indicative of social status and mirrors cultural attitudes toward the deceased or the group they represent. Here we report on the burials of six infants, including three who died at birth or shortly thereafter, from Later Stone Age (Iberomaurusian) levels at Grotte des Pigeons, Taforalt, in Morocco dating to ∼14,500 cal BP. Funerary treatment of the infants was equivalent to that of older individuals within the community, indicating an inclusive social status. The burials of two of the six infants, shown by previous aDNA analysis to be brother and sister, were overlain by ochre stained grinding stones that may have served as grave markers. In this case, a uniquely shared funerary treatment mirrored a close biological relationship, suggesting that kinship contributed to the patterning of funerary behavior within this Pleistocene burial assemblage.

Why structure matters.

Great care is needed when interpreting claims about the genetic basis of human variation based on data from genome-wide association studies.

Discovery of Potent, Efficient, and Selective Inhibitors of Phosphoinositide 3-Kinase δ through a Deconstruction and Regrowth Approach.

A deconstruction of previously reported phosphoinositide 3-kinase δ (PI3Kδ) inhibitors and subsequent regrowth led to the identification of a privileged fragment for PI3Kδ, which was exploited to deliver a potent, efficient, and selective lead series with a novel binding mode observed in the PI3Kδ crystal structure.

90,000 year-old specialised bone technology in the Aterian Middle Stone Age of North Africa.

The question of cognitive complexity in early Homo sapiens in North Africa is intimately tied to the emergence of the Aterian culture (~145 ka). One of the diagnostic indicators of cognitive complexity is the presence of specialised bone tools, however significant uncertainty remains over the manufacture and use of these artefacts within the Aterian techno-complex. In this paper we report on a bone artefact from Aterian Middle Stone Age (MSA) deposits in Dar es-Soltan 1 cave on the Atlantic coast of Morocco. It comes from a layer that can be securely dated to ~90 ka. The typological characteristics of this tool, which suggest its manufacture and use as a bone knife, are comparatively similar to other bone artefacts from dated Aterian levels at the nearby site of El Mnasra and significantly different from any other African MSA bone technology. The new find from Dar es-Soltan 1 cave combined with those from El Mnasra suggest the development of a bone technology unique to the Aterian.

Pleistocene North African genomes link Near Eastern and sub-Saharan African human populations.

North Africa is a key region for understanding human history, but the genetic history of its people is largely unknown. We present genomic data from seven 15,000-year-old modern humans, attributed to the Iberomaurusian culture, from Morocco. We find a genetic affinity with early Holocene Near Easterners, best represented by Levantine Natufians, suggesting a pre-agricultural connection between Africa and the Near East. We do not find evidence for gene flow from Paleolithic Europeans to Late Pleistocene North Africans. The Taforalt individuals derive one-third of their ancestry from sub-Saharan Africans, best approximated by a mixture of genetic components preserved in present-day West and East Africans. Thus, we provide direct evidence for genetic interactions between modern humans across Africa and Eurasia in the Pleistocene.

Identifying Bird Remains Using Ancient DNA Barcoding.

Bird remains that are difficult to identify taxonomically using morphological methods, are common in the palaeontological record. Other types of challenging avian material include artefacts and food items from endangered taxa, as well as remains from aircraft strikes. We here present a DNA-based method that enables taxonomic identification of bird remains, even from material where the DNA is heavily degraded. The method is based on the amplification and sequencing of two short variable parts of the 16S region in the mitochondrial genome. To demonstrate the applicability of this approach, we evaluated the method on a set of Holocene and Late Pleistocene postcranial bird bones from several palaeontological and archaeological sites in Europe with good success.

Range shifts or extinction? Ancient DNA and distribution modelling reveal past and future responses to climate warming in cold-adapted birds.

Global warming is predicted to cause substantial habitat rearrangements, with the most severe effects expected to occur in high-latitude biomes. However, one major uncertainty is whether species will be able to shift their ranges to keep pace with climate-driven environmental changes. Many recent studies on mammals have shown that past range contractions have been associated with local extinctions rather than survival by habitat tracking. Here, we have used an interdisciplinary approach that combines ancient DNA techniques, coalescent simulations and species distribution modelling, to investigate how two common cold-adapted bird species, willow and rock ptarmigan (Lagopus lagopus and Lagopus muta), respond to long-term climate warming. Contrary to previous findings in mammals, we demonstrate a genetic continuity in Europe over the last 20 millennia. Results from back-casted species distribution models suggest that this continuity may have been facilitated by uninterrupted habitat availability and potentially also the greater dispersal ability of birds. However, our predictions show that in the near future, some isolated regions will have little suitable habitat left, implying a future decrease in local populations at a scale unprecedented since the last glacial maximum.