The Neutral Theory in Light of Natural Selection

Identification of positions in human aldolase a that are neutral for apparent KM

According to evolutionary theory, many naturally-occurring amino acid substitutions are expected to be neutral or near-neutral, with little effect on protein structure or function. Accordingly, most changes observed in human exomes are also expected to be neutral. As such, accurate algorithms for identifying medically-relevant changes must discriminate rare, non-neutral substitutions against a background of neutral substitutions. However, due to historical biases in biochemical experiments, the data available to train and validate prediction algorithms mostly contains non-neutral substitutions, with few examples of neutral substitutions. Thus, available training sets have the opposite composition of the desired test sets. Towards improving a dataset of these critical negative controls, we have concentrated on identifying neutral positions - those positions for which most of the possible 19 amino acid substitutions have little effect on protein structure or function. Here, we used a strategy based on multiple sequence alignments to identify putative neutral positions in human aldolase A, followed by biochemical assays for 147 aldolase substitutions. Results showed that most variants had little effect on either the apparent Michaelis constant for substrate fructose-1,6-bisphosphate or its apparent cooperativity. Thus, these data are useful for training and validating prediction algorithms. In addition, we created a database of these and other biochemically characterized aldolase variants along with aldolase sequences and characteristics derived from sequence and structure analyses. This database is publicly available at https://github.com/liskinsk/Aldolase-variant-and-sequence-database.

Massively parallel experimental interrogation of natural variants in ancient signaling pathways reveals both purifying selection and local adaptation

The nature of standing genetic variation remains a central debate in population genetics, with differing perspectives on whether common variants are mostly neutral or have functional effects. We address this question by directly mapping the fitness effects of over 9,000 natural variants in the Ras/PKA and TOR/Sch9 pathways-key regulators of cell proliferation in eukaryotes-across four conditions in Saccharomyces cerevisiae. While many variants are neutral in our assay, on the order of 3,500 exhibited significant fitness effects. These non-neutral variants tend to be missense and affect conserved, more densely packed, and less solvent-exposed protein regions. They are also typically younger, occur at lower frequencies, and more often found in heterozygous states, suggesting they are subject to purifying selection. A substantial fraction of non-neutral variants showing strong fitness effects in our experiments, however, is present at high frequencies in the population. These variants show signs of local adaptation as they tend to be found specifically in domesticated strains adapted to human-made environments. Our findings support the view that while common variants are often neutral, a significant proportion have adaptive functional consequences and are driven into the population by local positive selection. This study highlights the potential to explore the functional effects of natural genetic variation on a genome scale with quantitative fitness measurements in the laboratory, bridging the gap between population genetics and functional genomics to understand evolutionary dynamics in the wild.

Phylogenetic analysis of selected species of Asteraceae on the basis of RPS 11 Gene

The Asteraceae family is a prominent group of flowering plants found across the globe, with the exception of Antarctica. The Asteraceae family is a largest flowering family pivotal group in plant evolution and diversification. Despite its importance, the genetic diversity within this family remains understudied. We focused on the rps-11 gene, a chloroplast marker previously utilized in phylogenetic studies, to investigate its potential in resolving Asteraceae relationships. The focus was on examining genetic diversity within sixteen specifically chosen species from the Asteraceae family. This assessment was based on an analysis of a chloroplast gene responsible for encoding the ribosomal protein of the smaller subunit 11 (rps 11). Nearly 417 bp of rps 11 gene was amplified, sequenced, computationally translated into amino acid sequence and the data was used for phylogenetic analysis as well as for rps 11 protein structure predictions. Based on nucleotide and amino acid sequences phylograms were drawn with the help of Molecular Evolutionary Genetic Analysis (MEGA 6), which exhibited clear genetic relationship among species under investigation. The observed genetic distance was 0.02 for Maximum likelihood tree based on nucleotide sequences whereas it was 0.05 for phylogram based on amino acid sequences. These values revealed that amino acid-based tree has demonstrated greater diversity among selected species in comparison to nucleotides-based tree. On the basis of pair wise distance calculations, genetic divergence values were found within the range of 0.015-0.309. Moreover, 3D protein modeling for rps 11 protein of sixteen selected species was also carried out by iterative threading assembly refinement (I-Tasser) software. The models exhibiting the highest C-score were picked with satisfactory plot statistics (> 90%) and structurally validated by PROCHECK. Furthermore, Ramachandran plots displayed that the rps 11 protein structures of Tagetes minuta, Xanthium strumarium, Lactuca sativa and Chrysanthemum indicum have best feature models with > 90% of residues in the allowed region and ≤ 2% in the disallowed region. The research is not enough to stand alone to validate the viability of the rps11 gene as a prospective contender for phylogenetic analysis. İn future we will focus on the maximum genetic diversity theory for phylogenetic analysis of this family.

A genome-to-proteome atlas charts natural variants controlling proteome diversity and forecasts their fitness effects

Despite abundant genomic and phenotypic data across individuals and environments, the functional impact of most mutations on phenotype remains unclear. Here, we bridge this gap by linking genome to proteome in 800 meiotic progeny from an intercross between two closely related Saccharomyces cerevisiae isolates adapted to distinct niches. Modest genetic distance between the parents generated remarkable proteomic diversity that was amplified in the progeny and captured by 6,476 genotype-protein associations, over 1,600 of which we resolved to single variants. Proteomic adaptation emerged through the combined action of numerous cis- and trans-regulatory mutations, a regulatory architecture that was conserved across the species. Notably, trans-regulatory variants often arose in proteins not traditionally associated with gene regulation, such as enzymes. Moreover, the proteomic consequences of mutations predicted fitness under various stresses. Our study demonstrates that the collective action of natural genetic variants drives dramatic proteome diversification, with molecular consequences that forecast phenotypic outcomes.

Specific amino acid changes correlate with pathogenic flavobacteria

Flavobacterium is a genus of microorganisms living in a variety of hosts and habitats across the globe. Some species are found in fish organs, and only a few, such as Flavobacterium psychrophilum and Flavobacterium columnare, cause severe disease and losses in fish farms. The evolution of flavobacteria that are pathogenic to fish is unknown, and the protein changes accountable for the selection of their colonization to fish have yet to be determined. A phylogenetic tree was constructed with the complete genomic sequences of 208 species of the Flavobacterium genus using 861 softcore genes. This phylogenetic analysis revealed clade CII comprising nine species, including five pathogenic species, and containing the most species that colonize fish. Thirteen specific amino acid changes were found to be conserved across 11 proteins within the CII clade compared with other clades, and these proteins were enriched in functions related to replication, recombination, and repair. Several of these proteins are known to be involved in pathogenicity and fitness adaptation in other bacteria. Some of the observed amino acid changes can be explained by preferential selection for certain codons and tRNA frequency. These results could help explain how species belonging to the CII clade adapt to fish environments.

Recurrent selection shapes the genomic landscape of differentiation between a pair of host-specialized haplodiploids that diverged with gene flow

Understanding the genetics of adaptation and speciation is critical for a complete picture of how biodiversity is generated and maintained. Heterogeneous genomic differentiation between diverging taxa is commonly documented, with genomic regions of high differentiation interpreted as resulting from differential gene flow, linked selection and reduced recombination rates. Disentangling the roles of each of these non-exclusive processes in shaping genome-wide patterns of divergence is challenging but will enhance our knowledge of the repeatability of genomic landscapes across taxa. Here, we combine whole-genome resequencing and genome feature data to investigate the processes shaping the genomic landscape of differentiation for a sister-species pair of haplodiploid pine sawflies, Neodiprion lecontei and Neodiprion pinetum. We find genome-wide correlations between genome features and summary statistics are consistent with pervasive linked selection, with patterns of diversity and divergence more consistently predicted by exon density and recombination rate than the neutral mutation rate (approximated by dS). We also find that both global and local patterns of FST, dXY and π provide strong support for recurrent selection as the primary selective process shaping variation across pine sawfly genomes, with some contribution from balancing selection and lineage-specific linked selection. Because inheritance patterns for haplodiploid genomes are analogous to those of sex chromosomes, we hypothesize that haplodiploids may be especially prone to recurrent selection, even if gene flow occurred throughout divergence. Overall, our study helps fill an important taxonomic gap in the genomic landscape literature and contributes to our understanding of the processes that shape genome-wide patterns of genetic variation.

Comparative genomics of dusky kob (Argyrosomus japonicus, Sciaenidae) conspecifics: Evidence for speciation and the genetic mechanisms underlying traits

Dusky kob (Argyrosomus japonicus) is a commercially important finfish, indigenous to South Africa, Australia, and China. Previous studies highlighted differences in genetic composition, life history, and morphology of the species across geographic regions. A draft genome sequence of 0.742 Gb (N50 = 5.49 Mb; BUSCO completeness = 97.8%) and 22,438 predicted protein-coding genes was generated for the South African (SA) conspecific. A comparison with the Chinese (CN) conspecific revealed a core set of 32,068 orthologous protein clusters across both genomes. The SA genome exhibited 440 unique clusters compared to 1928 unique clusters in the CN genome. Transportation and immune response processes were overrepresented among the SA accessory genome, whereas the CN accessory genome was enriched for immune response, DNA transposition, and sensory detection (FDR-adjusted p < 0.01). These unique clusters may represent an adaptive component of the species' pangenome that could explain population divergence due to differential environmental specialisation. Furthermore, 700 single-copy orthologues (SCOs) displayed evidence of positive selection between the SA and CN genomes, and globally these genomes shared only 92% similarity, suggesting they might be distinct species. These genes primarily play roles in metabolism and digestion, illustrating the evolutionary pathways that differentiate the species. Understanding these genomic mechanisms underlying adaptation and evolution within and between species provides valuable insights into growth and maturation of kob, traits that are particularly relevant to commercial aquaculture.

Selection leads to false inferences of introgression using popular methods

Detecting introgression between closely related populations or species is a fundamental objective in evolutionary biology. Existing methods for detecting migration and inferring migration rates from population genetic data often assume a neutral model of evolution. Growing evidence of the pervasive impact of selection on large portions of the genome across diverse taxa suggests that this assumption is unrealistic in most empirical systems. Further, ignoring selection has previously been shown to negatively impact demographic inferences (e.g. of population size histories). However, the impacts of biologically realistic selection on inferences of migration remain poorly explored. Here, we simulate data under models of background selection, selective sweeps, balancing selection, and adaptive introgression. We show that ignoring selection sometimes leads to false inferences of migration in popularly used methods that rely on the site frequency spectrum. Specifically, balancing selection and some models of background selection result in the rejection of isolation-only models in favor of isolation-with-migration models and lead to elevated estimates of migration rates. BPP, a method that analyzes sequence data directly, showed false positives for all conditions at recent divergence times, but balancing selection also led to false positives at medium-divergence times. Our results suggest that such methods may be unreliable in some empirical systems, such that new methods that are robust to selection need to be developed.

The Evolutionary Interplay of Somatic and Germline Mutation Rates

Novel sequencing technologies are making it increasingly possible to measure the mutation rates of somatic cell lineages. Accurate germline mutation rate measurement technologies have also been available for a decade, making it possible to assess how this fundamental evolutionary parameter varies across the tree of life. Here, we review some classical theories about germline and somatic mutation rate evolution that were formulated using principles of population genetics and the biology of aging and cancer. We find that somatic mutation rate measurements, while still limited in phylogenetic diversity, seem consistent with the theory that selection to preserve the soma is proportional to life span. However, germline and somatic theories make conflicting predictions regarding which species should have the most accurate DNA repair. Resolving this conflict will require carefully measuring how mutation rates scale with time and cell division and achieving a better understanding of mutation rate pleiotropy among cell types.

Polymorphism and the Red Queen: the selective maintenance of allelic variation in a deteriorating environment

Although allelic variation is ubiquitous in natural populations, our theoretical models are poor at predicting the existence and properties of these observed polymorphisms. In this study, inspired by Van Valen's Red Queen hypothesis, we modeled the effect of viability selection in a deteriorating environment on the properties of allelic variation in populations subject to recurrent mutation. In Monte Carlo simulations, we found that levels of polymorphism consistently built up over time. We censused the simulated populations after 10,000 generations of mutation and selection, revealing that, compared with models assuming a constant environment, the mean number of alleles was greater, as was the range of allele numbers. These results were qualitatively robust to the addition of genetic drift and to the relaxation of the assumption that the viabilities of phenogenotypes containing a new mutation are independent of each other (i.e. incorporating a model of generalized dominance). The broad range of allele numbers realized in the simulated populations-from monomorphisms to highly polymorphic populations-more closely corresponds to the observed range from numerous surveys of natural populations than previously found in theoretical studies. This match suggests that, contrary to the views of some writers, selection may actively maintain genetic variation in natural populations, particularly if the selective environment is gradually becoming harsher. Our simulations also generated many populations with heterozygote advantage, a mismatch with real data that implies that this selective property must arise extremely rarely in natural populations.

Population genomics highlights structural variations in local adaptation to saline coastal environments in woolly grape

Structural variations (SVs) are a feature of plant genomes that has been largely unexplored despite their significant impact on plant phenotypic traits and local adaptation to abiotic and biotic stress. In this study, we employed woolly grape (Vitis retordii), a species native to the tropical and subtropical regions of East Asia with both coastal and inland habitats, as a valuable model for examining the impact of SVs on local adaptation. We assembled a haplotype-resolved chromosomal reference genome for woolly grape, and conducted population genetic analyses based on whole-genome sequencing (WGS) data from coastal and inland populations. The demographic analyses revealed recent bottlenecks in all populations and asymmetric gene flow from the inland to the coastal population. In total, 1,035 genes associated with plant adaptive regulation for salt stress, radiation, and environmental adaptation were detected underlying local selection by SVs and SNPs in the coastal population, of which 37.29% and 65.26% were detected by SVs and SNPs, respectively. Candidate genes such as FSD2, RGA1, and AAP8 associated with salt tolerance were found to be highly differentiated and selected during the process of local adaptation to coastal habitats in SV regions. Our study highlights the importance of SVs in local adaptation; candidate genes related to salt stress and climatic adaptation to tropical and subtropical environments are important genomic resources for future breeding programs of grapevine and its rootstocks.

Population genetic structure of Pomacea canaliculata in China based on the COI and ITS1 genes

Comprehending the phylogeography of invasive organisms enhances our insight into their distribution dynamics, which is instrumental for the development of effective prevention and management strategies. In China, Pomacea canaliculata and Pomacea maculata are the two most widespread and damaging species of the non-native Pomacea spp.. Given this species' rapid spread throughout country, it is urgent to investigate the genetic diversity and structure of its different geographic populations, a task undertaken in the current study using the COI and ITS1 mitochondrial and ribosomal DNA genes, respectively. The result of this study, based on a nationwide systematic survey, a collection of Pomacea spp., and the identification of cryptic species, showed that there is a degree of genetic diversity and differentiation in P. canaliculata, and that all of its variations are mainly due to differences between individuals within different geographical populations. Indeed, this species contains multiple haplotypes, but none of them form a systematic geographical population structure. Furthermore, the COI gene exhibits higher genetic diversity than the ITS1 gene. Our study further clarifies the invasive pathways and dispersal patterns of P. canaliculata in China to provide a theoretical basis.

The protein domains of vertebrate species in which selection is more effective have greater intrinsic structural disorder

The nearly neutral theory of molecular evolution posits variation among species in the effectiveness of selection. In an idealized model, the census population size determines both this minimum magnitude of the selection coefficient required for deleterious variants to be reliably purged, and the amount of neutral diversity. Empirically, an "effective population size" is often estimated from the amount of putatively neutral genetic diversity and is assumed to also capture a species' effectiveness of selection. A potentially more direct measure of the effectiveness of selection is the degree to which selection maintains preferred codons. However, past metrics that compare codon bias across species are confounded by among-species variation in %GC content and/or amino acid composition. Here we propose a new Codon Adaptation Index of Species (CAIS), based on Kullback-Leibler divergence, that corrects for both confounders. We demonstrate the use of CAIS correlations, as well as the Effective Number of Codons, to show that the protein domains of more highly adapted vertebrate species evolve higher intrinsic structural disorder.

Shared evolutionary processes shape landscapes of genomic variation in the great apes

For at least the past 5 decades, population genetics, as a field, has worked to describe the precise balance of forces that shape patterns of variation in genomes. The problem is challenging because modeling the interactions between evolutionary processes is difficult, and different processes can impact genetic variation in similar ways. In this paper, we describe how diversity and divergence between closely related species change with time, using correlations between landscapes of genetic variation as a tool to understand the interplay between evolutionary processes. We find strong correlations between landscapes of diversity and divergence in a well-sampled set of great ape genomes, and explore how various processes such as incomplete lineage sorting, mutation rate variation, GC-biased gene conversion and selection contribute to these correlations. Through highly realistic, chromosome-scale, forward-in-time simulations, we show that the landscapes of diversity and divergence in the great apes are too well correlated to be explained via strictly neutral processes alone. Our best fitting simulation includes both deleterious and beneficial mutations in functional portions of the genome, in which 9% of fixations within those regions is driven by positive selection. This study provides a framework for modeling genetic variation in closely related species, an approach which can shed light on the complex balance of forces that have shaped genetic variation.

An efficient and robust ABC approach to infer the rate and strength of adaptation

Inferring the effects of positive selection on genomes remains a critical step in characterizing the ultimate and proximate causes of adaptation across species, and quantifying positive selection remains a challenge due to the confounding effects of many other evolutionary processes. Robust and efficient approaches for adaptation inference could help characterize the rate and strength of adaptation in nonmodel species for which demographic history, mutational processes, and recombination patterns are not currently well-described. Here, we introduce an efficient and user-friendly extension of the McDonald-Kreitman test (ABC-MK) for quantifying long-term protein adaptation in specific lineages of interest. We characterize the performance of our approach with forward simulations and find that it is robust to many demographic perturbations and positive selection configurations, demonstrating its suitability for applications to nonmodel genomes. We apply ABC-MK to the human proteome and a set of known virus interacting proteins (VIPs) to test the long-term adaptation in genes interacting with viruses. We find substantially stronger signatures of positive selection on RNA-VIPs than DNA-VIPs, suggesting that RNA viruses may be an important driver of human adaptation over deep evolutionary time scales.

Background Selection From Unlinked Sites Causes Nonindependent Evolution of Deleterious Mutations

Background selection describes the reduction in neutral diversity caused by selection against deleterious alleles at other loci. It is typically assumed that the purging of deleterious alleles affects linked neutral variants, and indeed simulations typically only treat a genomic window. However, background selection at unlinked loci also depresses neutral diversity. In agreement with previous analytical approximations, in our simulations of a human-like genome with a realistically high genome-wide deleterious mutation rate, the effects of unlinked background selection exceed those of linked background selection. Background selection reduces neutral genetic diversity by a factor that is independent of census population size. Outside of genic regions, the strength of background selection increases with the mean selection coefficient, contradicting the linked theory but in agreement with the unlinked theory. Neutral diversity within genic regions is fairly independent of the strength of selection. Deleterious genetic load among haploid individuals is underdispersed, indicating nonindependent evolution of deleterious mutations. Empirical evidence for underdispersion was previously interpreted as evidence for global epistasis, but we recover it from a non-epistatic model.

The contribution of gene flow, selection, and genetic drift to five thousand years of human allele frequency change

Genomic time series from experimental evolution studies and ancient DNA datasets offer us a chance to directly observe the interplay of various evolutionary forces. We show how the genome-wide variance in allele frequency change between two time points can be decomposed into the contributions of gene flow, genetic drift, and linked selection. In closed populations, the contribution of linked selection is identifiable because it creates covariances between time intervals, and genetic drift does not. However, repeated gene flow between populations can also produce directionality in allele frequency change, creating covariances. We show how to accurately separate the fraction of variance in allele frequency change due to admixture and linked selection in a population receiving gene flow. We use two human ancient DNA datasets, spanning around 5,000 y, as time transects to quantify the contributions to the genome-wide variance in allele frequency change. We find that a large fraction of genome-wide change is due to gene flow. In both cases, after correcting for known major gene flow events, we do not observe a signal of genome-wide linked selection. Thus despite the known role of selection in shaping long-term polymorphism levels, and an increasing number of examples of strong selection on single loci and polygenic scores from ancient DNA, it appears to be gene flow and drift, and not selection, that are the main determinants of recent genome-wide allele frequency change. Our approach should be applicable to the growing number of contemporary and ancient temporal population genomics datasets.

The evolution of chemodiversity in plants-From verbal to quantitative models

Plants harbour a great chemodiversity, that is diversity of specialised metabolites (SMs), at different scales. For instance, individuals can produce a large number of SMs, and populations can differ in their metabolite composition. Given the ecological and economic importance of plant chemodiversity, it is important to understand how it arises and is maintained over evolutionary time. For other dimensions of biodiversity, that is species diversity and genetic diversity, quantitative models play an important role in addressing such questions. Here, we provide a synthesis of existing hypotheses and quantitative models, that is mathematical models and computer simulations, for the evolution of plant chemodiversity. We describe each model's ingredients, that is the biological processes that shape chemodiversity, the scales it considers and whether it has been formalized as a quantitative model. Although we identify several quantitative models, not all are dynamic and many influential models have remained verbal. To fill these gaps, we outline our vision for the future of chemodiversity modelling. We identify quantitative models used for genetic variation that may be adapted for chemodiversity, and we present a flexible framework for the creation of individual-based models that address different scales of chemodiversity and combine different ingredients that bring this chemodiversity about.

Half a Century of Controversy: The Neutralist/Selectionist Debate in Molecular Evolution

The neutral and nearly neutral theories, introduced more than 50 yr ago, have raised and still raise passionate discussion regarding the forces governing molecular evolution and their relative importance. The debate, initially focused on the amount of within-species polymorphism and constancy of the substitution rate, has spread, matured, and now underlies a wide range of topics and questions. The neutralist/selectionist controversy has structured the field and influences the way molecular evolutionary scientists conceive their research.

Moderating the neutralist-selectionist debate: exactly which propositions are we debating, and which arguments are valid?

Half a century after its foundation, the neutral theory of molecular evolution continues to attract controversy. The debate has been hampered by the coexistence of different interpretations of the core proposition of the neutral theory, the 'neutral mutation-random drift' hypothesis. In this review, we trace the origins of these ambiguities and suggest potential solutions. We highlight the difference between the original, the revised and the nearly neutral hypothesis, and re-emphasise that none of them equates to the null hypothesis of strict neutrality. We distinguish the neutral hypothesis of protein evolution, the main focus of the ongoing debate, from the neutral hypotheses of genomic and functional DNA evolution, which for many species are generally accepted. We advocate a further distinction between a narrow and an extended neutral hypothesis (of which the latter posits that random non-conservative amino acid substitutions can cause non-ecological phenotypic divergence), and we discuss the implications for evolutionary biology beyond the domain of molecular evolution. We furthermore point out that the debate has widened from its initial focus on point mutations, and also concerns the fitness effects of large-scale mutations, which can alter the dosage of genes and regulatory sequences. We evaluate the validity of neutralist and selectionist arguments and find that the tested predictions, apart from being sensitive to violation of underlying assumptions, are often derived from the null hypothesis of strict neutrality, or equally consistent with the opposing selectionist hypothesis, except when assuming molecular panselectionism. Our review aims to facilitate a constructive neutralist-selectionist debate, and thereby to contribute to answering a key question of evolutionary biology: what proportions of amino acid and nucleotide substitutions and polymorphisms are adaptive?

Speciation, natural selection, and networks: three historians versus theoretical population geneticists

In 1913, the geneticist William Bateson called for a halt in studies of genetic phenomena until evolutionary fundamentals had been sufficiently addressed at the molecular level. Nevertheless, in the 1960s, the theoretical population geneticists celebrated a "modern synthesis" of the teachings of Mendel and Darwin, with an exclusive role for natural selection in speciation. This was supported, albeit with minor reservations, by historians Mark Adams and William Provine, who taught it to generations of students. In subsequent decades, doubts were raised by molecular biologists and, despite the deep influence of various mentors, Adams and Provine noted serious anomalies and began to question traditional "just-so-stories." They were joined in challenging the genetic orthodoxy by a scientist-historian, Donald Forsdyke, who suggested that a "collective variation" postulated by Darwin's young research associate, George Romanes, and a mysterious "residue" postulated by Bateson, might relate to differences in short runs of DNA bases (oligonucleotides). The dispute between a small network of historians and a large network of geneticists can be understood in the context of national politics. Contrasts are drawn between democracies, where capturing the narrative makes reversal difficult, and dictatorships, where overthrow of a supportive dictator can result in rapid reversal.

Comparative genomic analyses provide new insights into evolutionary history and conservation genomics of gorillas

Genome sequencing is a powerful tool to understand species evolutionary history, uncover genes under selection, which could be informative of local adaptation, and infer measures of genetic diversity, inbreeding and mutational load that could be used to inform conservation efforts. Gorillas, critically endangered primates, have received considerable attention and with the recently sequenced Bwindi mountain gorilla population, genomic data is now available from all gorilla subspecies and both mountain gorilla populations. Here, we reanalysed this rich dataset with a focus on evolutionary history, local adaptation and genomic parameters relevant for conservation. We estimate a recent split between western and eastern gorillas of 150,000-180,000 years ago, with gene flow around 20,000 years ago, primarily between the Cross River and Grauer's gorilla subspecies. This gene flow event likely obscures evolutionary relationships within eastern gorillas: after excluding putatively introgressed genomic regions, we uncover a sister relationship between Virunga mountain gorillas and Grauer's gorillas to the exclusion of Bwindi mountain gorillas. This makes mountain gorillas paraphyletic. Eastern gorillas are less genetically diverse and more inbred than western gorillas, yet we detected lower genetic load in the eastern species. Analyses of indels fit remarkably well with differences in genetic diversity across gorilla taxa as recovered with nucleotide diversity measures. We also identified genes under selection and unique gene variants specific for each gorilla subspecies, encoding, among others, traits involved in immunity, diet, muscular development, hair morphology and behavior. The presence of this functional variation suggests that the subspecies may be locally adapted. In conclusion, using extensive genomic resources we provide a comprehensive overview of gorilla genomic diversity, including a so-far understudied Bwindi mountain gorilla population, identify putative genes involved in local adaptation, and detect population-specific gene flow across gorilla species.

Genomic Architecture Predicts Tree Topology, Population Structuring, and Demographic History in Amazonian Birds

Geographic barriers are frequently invoked to explain genetic structuring across the landscape. However, inferences on the spatial and temporal origins of population variation have been largely limited to evolutionary neutral models, ignoring the potential role of natural selection and intrinsic genomic processes known as genomic architecture in producing heterogeneity in differentiation across the genome. To test how variation in genomic characteristics (e.g. recombination rate) impacts our ability to reconstruct general patterns of differentiation between species that cooccur across geographic barriers, we sequenced the whole genomes of multiple bird populations that are distributed across rivers in southeastern Amazonia. We found that phylogenetic relationships within species and demographic parameters varied across the genome in predictable ways. Genetic diversity was positively associated with recombination rate and negatively associated with species tree support. Gene flow was less pervasive in genomic regions of low recombination, making these windows more likely to retain patterns of population structuring that matched the species tree. We further found that approximately a third of the genome showed evidence of selective sweeps and linked selection, skewing genome-wide estimates of effective population sizes and gene flow between populations toward lower values. In sum, we showed that the effects of intrinsic genomic characteristics and selection can be disentangled from neutral processes to elucidate spatial patterns of population differentiation.

Allelic gene conversion softens selective sweeps

The prominence of positive selection, in which beneficial mutations are favored by natural selection and rapidly increase in frequency, is a subject of intense debate. Positive selection can result in selective sweeps, in which the haplotype(s) bearing the adaptive allele "sweep" through the population, thereby removing much of the genetic diversity from the region surrounding the target of selection. Two models of selective sweeps have been proposed: classical sweeps, or "hard sweeps", in which a single copy of the adaptive allele sweeps to fixation, and "soft sweeps", in which multiple distinct copies of the adaptive allele leave descendants after the sweep. Soft sweeps can be the outcome of recurrent mutation to the adaptive allele, or the presence of standing genetic variation consisting of multiple copies of the adaptive allele prior to the onset of selection. Importantly, soft sweeps will be common when populations can rapidly adapt to novel selective pressures, either because of a high mutation rate or because adaptive alleles are already present. The prevalence of soft sweeps is especially controversial, and it has been noted that selection on standing variation or recurrent mutations may not always produce soft sweeps. Here, we show that the inverse is true: selection on single-origin de novo mutations may often result in an outcome that is indistinguishable from a soft sweep. This is made possible by allelic gene conversion, which "softens" hard sweeps by copying the adaptive allele onto multiple genetic backgrounds, a process we refer to as a "pseudo-soft" sweep. We carried out a simulation study examining the impact of gene conversion on sweeps from a single de novo variant in models of human, Drosophila, and Arabidopsis populations. The fraction of simulations in which gene conversion had produced multiple haplotypes with the adaptive allele upon fixation was appreciable. Indeed, under realistic demographic histories and gene conversion rates, even if selection always acts on a single-origin mutation, sweeps involving multiple haplotypes are more likely than hard sweeps in large populations, especially when selection is not extremely strong. Thus, even when the mutation rate is low or there is no standing variation, hard sweeps are expected to be the exception rather than the rule in large populations. These results also imply that the presence of signatures of soft sweeps does not necessarily mean that adaptation has been especially rapid or is not mutation limited.

Unraveling the mitochondrial phylogenetic landscape of Thailand reveals complex admixture and demographic dynamics

The evolutionary dynamics of mitochondrial DNA within the Thai population were comprehensively explored with a specific focus on the influence of South Asian admixture. A total of 166 samples were collected through randomized sampling, ensuring a diverse representation. Our findings unveil substantial genetic and haplogroup diversity within the Thai population. We have identified 164 haplotypes categorized into 97 haplogroups, with a notable inclusion of 20 novel haplogroups. The distribution of haplogroups exhibited variations across different populations and countries. The central Thai population displayed a high diversity of haplogroups from both the M and N clades. Maternal lineage affinities were discerned between several Mainland Southeast Asia (MSEA) and South Asian populations, implying ancestral genetic connections and a substantial influence of South Asian women in establishing these relationships. f4-statistics indicates the presence of a Tibeto-Burman genetic component within the Mon population from Thailand. New findings demonstrate two phases of population expansion occurring 22,000-26,000 and 2500-3800 years ago, coinciding with the Last Glacial Maximum, and Neolithic demographic transition, respectively. This research significantly enhances our understanding of the maternal genetic history of Thailand and MSEA, emphasizing the influence of South Asian admixture. Moreover, it underscores the critical role of prior information, such as mutation rates, within the Bayesian framework for accurate estimation of coalescence times and inferring demographic history.

Comprehensive detection of structural variation and transposable element differences between wild type laboratory lineages of C. elegans

Genomic structural variations (SVs) and transposable elements (TEs) can be significant contributors to genome evolution, altered gene expression, and risk of genetic diseases. Recent advancements in long-read sequencing have greatly improved the quality of de novo genome assemblies and enhanced the detection of sequence variants at the scale of hundreds or thousands of bases. Comparisons between two diverged wild isolates of Caenorhabditis elegans, the Bristol and Hawaiian strains, have been widely utilized in the analysis of small genetic variations. Genetic drift, including SVs and rearrangements of repeated sequences such as TEs, can occur over time from long-term maintenance of wild type isolates within the laboratory. To comprehensively detect both large and small structural variations as well as TEs due to genetic drift, we generated de novo genome assemblies and annotations for each strain from our lab collection using both long- and short-read sequencing and compared our assemblies and annotations with that of other lab wild type strains. Within our lab assemblies, we annotate over 3.1Mb of sequence divergence between the Bristol and Hawaiian isolates: 337,584 SNPs, 94,503 small insertion-deletions (<50bp), and 4,334 structural variations (>50bp). Further, we define the location and movement of specific DNA TEs between N2 Bristol and CB4856 Hawaiian wild type isolates. Specifically, we find the N2 Bristol genome has 20.6% more TEs from the Tc1/mariner family than the CB4856 Hawaiian genome. Moreover, we identified Zator elements as the most abundant and mobile TE family in the genome. Using specific TE sequences with unique SNPs, we also identify 38 TEs that moved intrachromosomally and 9 TEs that moved interchromosomally between the N2 Bristol and CB4856 Hawaiian genomes. By comparing the de novo genome assembly of our lab collection Bristol isolate to the VC2010 Bristol assembly, we also reveal that lab lineages display over 2 Mb of total variation: 1,162 SNPs, 1,528 indels, and 897 SVs with 95% of the variation due to SVs. Overall, our work demonstrates the unique contribution of SVs and TEs to variation and genetic drift between wild type laboratory strains assumed to be isogenic despite growing evidence of genetic drift and phenotypic variation.

Shared evolutionary processes shape landscapes of genomic variation in the great apes

For at least the past five decades population genetics, as a field, has worked to describe the precise balance of forces that shape patterns of variation in genomes. The problem is challenging because modelling the interactions between evolutionary processes is difficult, and different processes can impact genetic variation in similar ways. In this paper, we describe how diversity and divergence between closely related species change with time, using correlations between landscapes of genetic variation as a tool to understand the interplay between evolutionary processes. We find strong correlations between landscapes of diversity and divergence in a well sampled set of great ape genomes, and explore how various processes such as incomplete lineage sorting, mutation rate variation, GC-biased gene conversion and selection contribute to these correlations. Through highly realistic, chromosome-scale, forward-in-time simulations we show that the landscapes of diversity and divergence in the great apes are too well correlated to be explained via strictly neutral processes alone. Our best fitting simulation includes both deleterious and beneficial mutations in functional portions of the genome, in which 9% of fixations within those regions is driven by positive selection. This study provides a framework for modelling genetic variation in closely related species, an approach which can shed light on the complex balance of forces that have shaped genetic variation.

An efficient and robust ABC approach to infer the rate and strength of adaptation

Inferring the effects of positive selection on genomes remains a critical step in characterizing the ultimate and proximate causes of adaptation across species, and quantifying positive selection remains a challenge due to the confounding effects of many other evolutionary processes. Robust and efficient approaches for adaptation inference could help characterize the rate and strength of adaptation in non-model species for which demographic history, mutational processes, and recombination patterns are not currently well-described. Here, we introduce an efficient and user-friendly extension of the McDonald-Kreitman test (ABC-MK) for quantifying long-term protein adaptation in specific lineages of interest. We characterize the performance of our approach with forward simulations and find that it is robust to many demographic perturbations and positive selection configurations, demonstrating its suitability for applications to non-model genomes. We apply ABC-MK to the human proteome and a set of known Virus Interacting Proteins (VIPs) to test the long-term adaptation in genes interacting with viruses. We find substantially stronger signatures of positive selection on RNA-VIPs than DNA-VIPs, suggesting that RNA viruses may be an important driver of human adaptation over deep evolutionary time scales.

Genetic diversity and ancestry of the Khmuic-speaking ethnic groups in Thailand: a genome-wide perspective

The Khmuic-speaking populations are believed to be the descendants of one of the earliest groups to settle in Mainland Southeast Asia. In Thailand, there are two agricultural Khmuic-speaking ethnic groups, the Khamu and Lua (Htin). These peoples primarily reside in scattered locations along the mountainous Thailand-Laos border in Nan province. In this study, we conducted genome-wide SNP analysis on 81 individuals from three Khamu and two Lua villages in northern Thailand. Our findings revealed that both the Khamu and Lua groups possess genetic structures that are distinct from other ethnicities in Southeast Asia, indicating a unique history of migration and settlement. Within the Khmuic group, the Khamu populations living in different locations exhibited similar genetic structures and displayed genetic affinities only with some hill-tribes and Tai-Kadai (Kra-Dai)-speaking groups in Thailand, suggesting potential intermixing or cultural exchange. Furthermore, the Lua people displayed a distinctive population structure, which could be attributed to the founder effect and endogamous marriage practices. Additionally, we discovered a relationship between the Khmuic-speaking populations in Thailand and a Neolithic ancient sample obtained from the Tham Pha Ling archaeological site in Laos. This study provides new insight into genetic substructure within the Khmuic-speaking people and their potential relationship to the indigenous inhabitants of Mainland Southeast Asia.

Population evolution of seagrasses returning to the ocean

Seagrasses are higher flowering plants that live entirely in marine environments, with the greatest habitat variation occurring from land to sea. Genetic structure or population differentiation history is a hot topic in evolutionary biology, which is of great significance for understanding speciation. Genetic information is obtained from geographically distributed subpopulations, different subspecies, or strains of the same species using next-generation sequencing techniques. Genetic variation is identified by comparison with reference genomes. Genetic diversity is explored using population structure, principal component analysis (PCA), and phylogenetic relationships. Patterns of population genetic differentiation are elucidated by combining the isolation by distance (IBD) model, linkage disequilibrium levels, and genetic statistical analysis. Demographic history is simulated using effective population size, divergence time, and site frequency spectrum (SFS). Through various population genetic analyses, the genetic structure and historical population dynamics of seagrass can be clarified, and their evolutionary processes can be further explored at the molecular level to understand how evolutionary processes contributed to the formation of early ecological species and provide data support for seagrass conservation.

Kimura's Theory of Non-Adaptive Radiation and Peto's Paradox: A Missing Link?

Karyotype diversity reflects genome integrity and stability. A strong correlation between karyotype diversity and species richness, meaning the number of species in a phylogenetic clade, was first reported in mammals over forty years ago: in mammalian phylogenetic clades, the standard deviation of karyotype diversity (KD) closely corresponded to species richness (SR) at the order level. These initial studies, however, did not control for phylogenetic signal, raising the possibility that the correlation was due to phylogenetic relatedness among species in a clade. Accordingly, karyotype diversity trivially reflects species richness simply as a passive consequence of adaptive radiation. A more recent study in mammals controlled for phylogenetic signals and established the correlation as phylogenetically independent, suggesting that species richness cannot, in itself, explain the observed corresponding karyotype diversity. The correlation is, therefore, remarkable because the molecular mechanisms contributing to karyotype diversity are evolutionarily independent of the ecological mechanisms contributing to species richness. Recently, it was shown in salamanders that the two processes generating genome size diversity and species richness were indeed independent and operate in parallel, suggesting a potential non-adaptive, non-causal but biologically meaningful relationship. KD depends on mutational input generating genetic diversity and reflects genome stability, whereas species richness depends on ecological factors and reflects natural selection acting on phenotypic diversity. As mutation and selection operate independently and involve separate and unrelated evolutionary mechanisms-there is no reason a priori to expect such a strong, let alone any, correlation between KD and SR. That such a correlation exists is more consistent with Kimura's theory of non-adaptive radiation than with ecologically based adaptive theories of macro-evolution, which are not excluded in Kimura's non-adaptive theory. The following reviews recent evidence in support of Kimura's proposal, and other findings that contribute to a wider understanding of the molecular mechanisms underlying the process of non-adaptive radiation.

The genomic basis of temporal niche evolution in a diurnal rodent

Patterns of diel activity-how animals allocate their activity throughout the 24-h daily cycle-play key roles in shaping the internal physiology of an animal and its relationship with the external environment.1,2,3,4,5 Although shifts in diel activity patterns have occurred numerous times over the course of vertebrate evolution,6 the genomic correlates of such transitions remain unknown. Here, we use the African striped mouse (Rhabdomys pumilio), a species that transitioned from the ancestrally nocturnal diel niche of its close relatives to a diurnal one,7,8,9,10,11 to define patterns of naturally occurring molecular variation in diel niche traits. First, to facilitate genomic analyses, we generate a chromosome-level genome assembly of the striped mouse. Next, using transcriptomics, we show that the switch to daytime activity in this species is associated with a realignment of daily rhythms in peripheral tissues with respect to the light:dark cycle and the central circadian clock. To uncover selection pressures associated with this temporal niche shift, we perform comparative genomic analyses with closely related rodent species and find evidence of relaxation of purifying selection on striped mouse genes in the rod phototransduction pathway. In agreement with this, electroretinogram measurements demonstrate that striped mice have functional differences in dim-light visual responses compared with nocturnal rodents. Taken together, our results show that striped mice have undergone a drastic change in circadian organization and provide evidence that the visual system has been a major target of selection as this species transitioned to a novel temporal niche.

Neutral theory and beyond: A systematic review of molecular evolution education

Molecular evolution-including the neutral theory of molecular evolution-is a major sub-discipline of evolution and is widely taught in undergraduate evolution courses. However, despite its ubiquity, there have not been any previous attempts to compile and review the molecular evolution education literature. Here, we draw upon the framework proposed in a past literature review examining the broader evolution education landscape to conduct a literature review of papers related to molecular evolution education, classifying the contributions of such papers to evolution pedagogy as well as evolution education research. We find that there remains very limited coverage of molecular evolution in the education literature, with existing papers focusing primarily on providing new instructional modules and strategies for teaching molecular evolution. Our work suggests several areas of critical need as well as opportunities to advance evolution education and evolution education research, including compiling instructional goals for the sub-discipline, developing validated assessments, and investigating student thinking related to molecular evolution. We conclude by providing general strategies, advice, and a novel curricular activity for teaching molecular evolution and the neutral theory of molecular evolution.

Timesweeper: accurately identifying selective sweeps using population genomic time series

Despite decades of research, identifying selective sweeps, the genomic footprints of positive selection, remains a core problem in population genetics. Of the myriad methods that have been developed to tackle this task, few are designed to leverage the potential of genomic time-series data. This is because in most population genetic studies of natural populations, only a single period of time can be sampled. Recent advancements in sequencing technology, including improvements in extracting and sequencing ancient DNA, have made repeated samplings of a population possible, allowing for more direct analysis of recent evolutionary dynamics. Serial sampling of organisms with shorter generation times has also become more feasible due to improvements in the cost and throughput of sequencing. With these advances in mind, here we present Timesweeper, a fast and accurate convolutional neural network-based tool for identifying selective sweeps in data consisting of multiple genomic samplings of a population over time. Timesweeper analyzes population genomic time-series data by first simulating training data under a demographic model appropriate for the data of interest, training a one-dimensional convolutional neural network on said simulations, and inferring which polymorphisms in this serialized data set were the direct target of a completed or ongoing selective sweep. We show that Timesweeper is accurate under multiple simulated demographic and sampling scenarios, identifies selected variants with high resolution, and estimates selection coefficients more accurately than existing methods. In sum, we show that more accurate inferences about natural selection are possible when genomic time-series data are available; such data will continue to proliferate in coming years due to both the sequencing of ancient samples and repeated samplings of extant populations with faster generation times, as well as experimentally evolved populations where time-series data are often generated. Methodological advances such as Timesweeper thus have the potential to help resolve the controversy over the role of positive selection in the genome. We provide Timesweeper as a Python package for use by the community.

Negative frequency dependent selection unites ecology and evolution

From genes to communities, understanding how diversity is maintained remains a fundamental question in biology. One challenging to identify, yet potentially ubiquitous, mechanism for the maintenance of diversity is negative frequency dependent selection (NFDS), which occurs when entities (e.g., genotypes, life history strategies, species) experience a per capita reduction in fitness with increases in relative abundance. Because NFDS allows rare entities to increase in frequency while preventing abundant entities from excluding others, we posit that negative frequency dependent selection plays a central role in the maintenance of diversity. In this review, we relate NFDS to coexistence, identify mechanisms of NFDS (e.g., mutualism, predation, parasitism), review strategies for identifying NFDS, and distinguish NFDS from other mechanisms of coexistence (e.g., storage effects, fluctuating selection). We also emphasize that NFDS is a key place where ecology and evolution intersect. Specifically, there are many examples of frequency dependent processes in ecology, but fewer cases that link this process to selection. Similarly, there are many examples of selection in evolution, but fewer cases that link changes in trait values to negative frequency dependence. Bridging these two well-developed fields of ecology and evolution will allow for mechanistic insights into the maintenance of diversity at multiple levels.

Evolutionarily stable strategy analysis and its links to demography and genetics through invasion fitness

Evolutionarily stable strategy (ESS) analysis pioneered by Maynard Smith and Price took off in part because it often does not require explicit assumptions about the genetics and demography of a population in contrast to population genetic models. Though this simplicity is useful, it obscures the degree to which ESS analysis applies to populations with more realistic genetics and demography: for example, how does ESS analysis handle complexities such as kin selection, group selection and variable environments when phenotypes are affected by multiple genes? In this paper, I review the history of the ESS concept and show how early uncertainty about the method lead to important mathematical theory linking ESS analysis to general population genetic models. I use this theory to emphasize the link between ESS analysis and the concept of invasion fitness. I give examples of how invasion fitness can measure kin selection, group selection and the evolution of linked modifier genes in response to variable environments. The ESSs in these examples depend crucially on demographic and genetic parameters, which highlights how ESS analysis will continue to be an important tool in understanding evolutionary patterns as new models address the increasing abundance of genetic and long-term demographic data in natural populations. This article is part of the theme issue 'Half a century of evolutionary games: a synthesis of theory, application and future directions'.

Inferring the mode and strength of ongoing selection

Genome sequence data are no longer scarce. The UK Biobank alone comprises 200,000 individual genomes, with more on the way, leading the field of human genetics toward sequencing entire populations. Within the next decades, other model organisms will follow suit, especially domesticated species such as crops and livestock. Having sequences from most individuals in a population will present new challenges for using these data to improve health and agriculture in the pursuit of a sustainable future. Existing population genetic methods are designed to model hundreds of randomly sampled sequences but are not optimized for extracting the information contained in the larger and richer data sets that are beginning to emerge, with thousands of closely related individuals. Here we develop a new method called trio-based inference of dominance and selection (TIDES) that uses data from tens of thousands of family trios to make inferences about natural selection acting in a single generation. TIDES further improves on the state of the art by making no assumptions regarding demography, linkage, or dominance. We discuss how our method paves the way for studying natural selection from new angles.

Fluctuating selection and the determinants of genetic variation

Recent studies of cosmopolitan Drosophila populations have found hundreds to thousands of genetic loci with seasonally fluctuating allele frequencies, bringing temporally fluctuating selection to the forefront of the historical debate surrounding the maintenance of genetic variation in natural populations. Numerous mechanisms have been explored in this longstanding area of research, but these exciting empirical findings have prompted several recent theoretical and experimental studies that seek to better understand the drivers, dynamics, and genome-wide influence of fluctuating selection. In this review, we evaluate the latest evidence for multilocus fluctuating selection in Drosophila and other taxa, highlighting the role of potential genetic and ecological mechanisms in maintaining these loci and their impacts on neutral genetic variation.

Selection signatures for fiber production in commercial species: A review

Natural fibers derived from diverse animal species have gained increased attention in recent years due to their favorable environmental effects, long-term sustainability benefits, and remarkable physical and mechanical properties that make them valuable raw materials used for textile and non-textile production. Domestication and selective breeding for the economically significant fiber traits play an imperative role in shaping the genomes and, thus, positively impact the overall productivity of the various fiber-producing species. These selection pressures leave unique footprints on the genome due to alteration in the allelic frequencies at specific loci, characterizing selective sweeps. Recent advances in genomics have enabled the discovery of selection signatures across the genome using a variety of methods. The increased demand for 'green products' manufactured from natural fibers necessitates a detailed investigation of the genomes of the various fiber-producing plant and animal species to identify the candidate genes associated with important fiber attributes such as fiber diameter/fineness, color, length, and strength, among others. The objective of this review is to present a comprehensive overview of the concept of selection signature and selective sweeps, discuss the main methods used for its detection, and address the selection signature studies conducted so far in the diverse fiber-producing animal species.

What Has Genomics Taught An Evolutionary Biologist?

Genomics, an interdisciplinary field of biology on the structure, function, and evolution of genomes, has revolutionized many subdisciplines of life sciences, including my field of evolutionary biology, by supplying huge data, bringing high-throughput technologies, and offering a new approach to biology. In this review, I describe what I have learned from genomics and highlight the fundamental knowledge and mechanistic insights gained. I focus on three broad topics that are central to evolutionary biology and beyond-variation, interaction, and selection-and use primarily my own research and study subjects as examples. In the next decade or two, I expect that the most important contributions of genomics to evolutionary biology will be to provide genome sequences of nearly all known species on Earth, facilitate high-throughput phenotyping of natural variants and systematically constructed mutants for mapping genotype-phenotype-fitness landscapes, and assist the determination of causality in evolutionary processes using experimental evolution.

Demography and linked selection interact to shape the genomic landscape of codistributed woodpeckers during the Ice Age

The influence of genetic drift on population dynamics during Pleistocene glacial cycles is well understood, but the role of selection in shaping patterns of genomic variation during these events is less explored. We resequenced whole genomes to investigate how demography and natural selection interact to generate the genomic landscapes of Downy and Hairy Woodpecker, species codistributed in previously glaciated North America. First, we explored the spatial and temporal patterns of genomic diversity produced by neutral evolution. Next, we tested (i) whether levels of nucleotide diversity along the genome are correlated with intrinsic genomic properties, such as recombination rate and gene density, and (ii) whether different demographic trajectories impacted the efficacy of selection. Our results revealed cycles of bottleneck and expansion, and genetic structure associated with glacial refugia. Nucleotide diversity varied widely along the genome, but this variation was highly correlated between the species, suggesting the presence of conserved genomic features. In both taxa, nucleotide diversity was positively correlated with recombination rate and negatively correlated with gene density, suggesting that linked selection played a role in reducing diversity. Despite strong fluctuations in effective population size, the maintenance of relatively large populations during glaciations may have facilitated selection. Under these conditions, we found evidence that the individual demographic trajectory of populations modulated linked selection, with purifying selection being more efficient in removing deleterious alleles in large populations. These results highlight that while genome-wide variation reflects the expected signature of demographic change during climatic perturbations, the interaction of multiple processes produces a predictable and highly heterogeneous genomic landscape.

Do Noncoding and Coding Sites in Angiosperm Chloroplast DNA Have Different Mutation Processes?

Fourfold degenerate sites within coding regions and intergenic sites have both been used as estimates of neutral evolution. In chloroplast DNA, the pattern of substitution at intergenic sites is strongly dependent on the composition of the surrounding hexanucleotide composed of the three base pairs on each side, which suggests that the mutation process is highly context-dependent in this genome. This study examines the context-dependency of substitutions at fourfold degenerate sites in protein-coding regions and compares the pattern to what has been observed at intergenic sites. Overall, there is strong similarity between the two types of sites, but there are some intriguing differences. One of these is that substitutions of G and C are significantly higher at fourfold degenerate sites across a range of contexts. In fact, A → T and T → A substitutions are the only substitution types that occur at a lower rate at fourfold degenerate sites. The data are not consistent with selective constraints being responsible for the difference in substitution patterns between intergenic and fourfold degenerate sites. Rather, it is suggested that the difference may be a result of different epigenetic modifications that result in slightly different mutation patterns in coding and intergenic DNA.

Population Genomic Time Series Data of a Natural Population Suggests Adaptive Tracking of Fluctuating Environmental Changes

Natural populations are constantly exposed to fluctuating environmental changes that negatively affect their fitness in unpredictable ways. While theoretical models show the possibility of counteracting these environmental changes through rapid evolutionary adaptations, there have been few empirical studies demonstrating such adaptive tracking in natural populations. Here, we analyzed environmental data, fitness-related phenotyping and genomic time-series data sampled over 3 years from a natural Chironomus riparius (Diptera, Insecta) population to address this question. We show that the population's environment varied significantly on the time scale of the sampling in many selectively relevant dimensions, independently of each other. Similarly, phenotypic fitness components evolved significantly on the same temporal scale (mean 0.32 Haldanes), likewise independent from each other. The allele frequencies of 367,446 SNPs across the genome showed evidence of positive selection. Using temporal correlation of spatially coherent allele frequency changes revealed 35,574 haplotypes with more than one selected SNP. The mean selection coefficient for these haplotypes was 0.30 (s.d. = 0.68). The frequency changes of these haplotypes clustered in 46 different temporal patterns, indicating concerted, independent evolution of many polygenic traits. Nine of these patterns were strongly correlated with measured environmental variables. Enrichment analysis of affected genes suggested the implication of a wide variety of biological processes. Thus, our results suggest overall that the natural population of C. riparius tracks environmental change through rapid polygenic adaptation in many independent dimensions. This is further evidence that natural selection is pervasive at the genomic level and that evolutionary and ecological time scales may not differ at all, at least in some organisms.

A common intronic single nucleotide variant modifies PKD1 expression level

Autosomal dominant polycystic kidney disease (ADPKD), caused by mutations in PKD1 and PKD2 (PKD1/2), has unexplained phenotypic variability likely affected by environmental and other genetic factors. Approximately 10% of individuals with ADPKD phenotype have no causal mutation detected, possibly due to unrecognized risk variants of PKD1/2. This study was designed to identify risk variants of PKD genes through population genetic analyses. We used Wright's F-statistics (Fst) to evaluate common single nucleotide variants (SNVs) potentially favored by positive natural selection in PKD1 from 1000 Genomes Project (1KG) and genotyped 388 subjects from the Rogosin Institute ADPKD Data Repository. The variants with >90th percentile Fst scores underwent further investigation by in silico analysis and molecular genetics analyses. We identified a deep intronic SNV, rs3874648G> A, located in a conserved binding site of the splicing regulator Tra2-β in PKD1 intron 30. Reverse-transcription PCR (RT-PCR) of peripheral blood leukocytes (PBL) from an ADPKD patient homozygous for rs3874648-A identified an atypical PKD1 splice form. Functional analyses demonstrated that rs3874648-A allele increased Tra2-β binding affinity and activated a cryptic acceptor splice-site, causing a frameshift that introduced a premature stop codon in mRNA, thereby decreasing PKD1 full-length transcript level. PKD1 transcript levels were lower in PBL from rs3874648-G/A carriers than in rs3874648-G/G homozygotes in a small cohort of normal individuals and patients with PKD2 inactivating mutations. Our findings indicate that rs3874648G > A is a PKD1 expression modifier attenuating PKD1 expression through Tra2-β, while the derived G allele advantageously maintains PKD1 expression and is predominant in all subpopulations.

Genetic slippage after sex maintains diversity for parasite resistance in a natural host population

Although parasite-mediated selection is a major driver of host evolution, its influence on genetic variation for parasite resistance is not yet well understood. We monitored resistance in a large population of the planktonic crustacean Daphnia magna over 8 years, as it underwent yearly epidemics of the bacterial pathogen Pasteuria ramosa. We observed cyclic dynamics of resistance: Resistance increased throughout the epidemics, but susceptibility was restored each spring when hosts hatched from sexual resting stages. Host resting stages collected across the year showed that largely resistant host populations can produce susceptible sexual offspring. A genetic model of resistance developed for this host-parasite system, based on multiple loci and strong epistasis, is in partial agreement with our findings. Our results reveal that, despite strong selection for resistance in a natural host population, genetic slippage after sexual reproduction can be a strong factor for the maintenance of genetic diversity of host resistance.

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.

Phylogenetic, Allometric, and Ecological Factors Affecting Morphological Variation in the Scapula and Humerus of Spiny Rats (Rodentia: Echimyidae)

Locomotion, as a fundamental function in mammals directly associated with the use of ecological resources, is expected to have anatomical structures functionally committed that evolved under intense selective pressure, possibly carrying specializations for different locomotor habits. Among caviomorph rodents, the family Echimyidae stands out for having the greatest species richness, with relatively well-resolved phylogenetic relationships, wide variation in body mass, and remarkable diversity of locomotor habits, including arboreal, scansorial, semi-aquatic, semifossorial, and terrestrial forms. Thus, Echimyidae constitutes a promising model for understanding how phylogenetic, allometric, and ecological factors affect the evolution of postcranial structures directly linked to locomotor function. We investigated the influence of these three factors on scapular and humeral morphological variation in 38 echimyid species using two-dimensional geometric morphometry and phylogenetically informed comparative methods. Scapular and humeral shape variation had a low correlation with body mass and structure size, conveying a small or negligible allometric effect. Conversely, a significant moderate to strong phylogenetic signal was detected in both structures, suggesting that an important part of their morphometric variation results from shared evolutionary history. Notably, morphological variation of the scapula was extensively structured by phylogeny, without the marked influence of locomotor habits, suggesting that its shape may be a suitable taxonomic marker. Finally, locomotor habits were important in structuring the morphological variation of the humerus. Our results suggest that the morphologies of the scapula and humerus, despite being anatomically and functionally interconnected, were differentially shaped by ecological factors associated with locomotor habits.

dN/dS-H, a New Test to Distinguish Different Selection Modes in Protein Evolution and Cancer Evolution

One of the most popular measures in the analysis of protein sequence evolution is the ratio of nonsynonymous distance (d_N) to synonymous distance (d_S). Under the assumption that synonymous substitutions in the coding region are selectively neutral, the d_N/d_S ratio can be used to statistically detect the adaptive evolution (or purifying selection) if d_N/d_S > 1 (or d_N/d_S < 1) significantly. However, due to strong structural constraints and/or variable functional constraints imposed on amino acid sites, most encoding genes in most species have demonstrated d_N/d_S < 1. Consequently, the statistical power for testing d_N/d_S = 1 may be insufficient to distinguish between different selection modes. In this paper, we propose a more powerful test, called d_N/d_S-H, in which a new parameter H, a relative measure of rate variation among sites, was introduced. Given the condition of strong purifying selections at some sites, the d_N/d_S-H model predicts d_N/d_S = 1-H for neutral evolution, d_N/d_S < 1-H for nearly neutral selection, and d_N/d_S > 1-H for adaptive evolution. The potential of this new method for resolving the neutral-adaptive debates is illustrated by the protein sequence evolution in vertebrates, Drosophila and yeasts, as well as somatic cancer evolution (specialized as the C_N/C_S-H test).

How Can We Resolve Lewontin's Paradox?

We discuss the genetic, demographic, and selective forces that are likely to be at play in restricting observed levels of DNA sequence variation in natural populations to a much smaller range of values than would be expected from the distribution of census population sizes alone-Lewontin's Paradox. While several processes that have previously been strongly emphasized must be involved, including the effects of direct selection and genetic hitchhiking, it seems unlikely that they are sufficient to explain this observation without contributions from other factors. We highlight a potentially important role for the less-appreciated contribution of population size change; specifically, the likelihood that many species and populations may be quite far from reaching the relatively high equilibrium diversity values that would be expected given their current census sizes.

Discussions of the "Not So Fit": How Ableism Limits Diverse Thought and Investigative Potential in Evolutionary Biology

AbstractEvolutionary biology and many of its foundational concepts are grounded in a history of ableism and eugenics. The field has not made a concerted effort to divest our concepts and investigative tools from this fraught history, and as a result, an ableist investigative lens has persisted in present-day evolutionary research, limiting the scope of research and harming the ability to communicate and synthesize knowledge about evolutionary processes. This failure to divest from our eugenicist and ableist history has harmed progress in evolutionary biology and allowed principles from evolutionary biology to continue to be weaponized against marginalized communities in the modern day. To rectify this problem, scholars in evolutionary research must come to terms with how the history of the field has influenced their investigations and work to establish a new framework for defining and investigating concepts such as selection and fitness.