Lewontin did not commit Lewontin's fallacy, his critics do: Why racial taxonomy is not useful for the scientific study of human variation.

In 1972, R.C. Lewontin concluded that it follows from the fact that the large majority of human genetic variation (≈ 85%) is among individuals within local populations that racial taxonomy is unjustified. Three decades later, Edwards demonstrated that while the accuracy with which individuals may be assigned to groups is poor for a single locus, consideration of multi-locus data allows for highly accurate assignments. Edwards concluded that Lewontin's dismissal of racial taxonomy was unwarranted. Edwards misidentified the aim of Lewontin's critique, which was directed at the utility of racial classification and not at assigning individuals to groups using genetic data. Moreover, Edwards conflated distinct kinds of correlation when sketching out his argument. If we follow Edwards' argument to its natural terminus, it becomes clear that it is consideration of all of the correlation structure among local groups in human genetic data that renders racial taxonomy scientifically useless. Lewontin considers the correlation structure relevant to his analysis of racial taxonomy and does not make his eponymous misstep. Rather, critics of Lewontin who use racial taxonomies in their work are the primary offenders when it comes to committing Lewontin's fallacy.

Sociobiology on Screen. The Controversy Through the Lens of Sociobiology: Doing What Comes Naturally.

When the sociobiology debate erupted in 1975, there were almost too many contributions to the heated exchanges between sociobiologists and their critics to count. In the fall of 1976, a Canadian educational film entitled Sociobiology: Doing What Comes Naturally sparked further controversy due to its graphic visuals and outrageous narration. While critics claimed the film was a promotional tool to further the sociobiological agenda in educational settings, sociobiologists quickly distanced themselves from the film and, in turn, accused the critics of consciously misrepresenting sociobiology by organizing showings of the film. Using audio, video, archival, and published sources, this paper explores the complicated history of Sociobiology: Doing What Comes Naturally and demonstrates how the public debate about the film reflects the positions, polemics, and polarization of the sociobiology debate as a whole.

Richard Lewontin and the "complications of linkage".

During the 1960s and 1970s population geneticists pushed beyond models of single genes to grapple with the effect on evolution of multiple genes associated by linkage. The resulting models of multiple interacting loci suggested that blocks of genes, maybe even entire chromosomes or the genome itself, should be treated as a unit. In this context, Richard Lewontin wrote his famous 1974 book The Genetic Basis of Evolutionary Change, which concludes with an argument for considering the entire genome as the unit of selection as a result of linkage. Why did Lewontin and others devote so much intellectual energy to the "complications of linkage" in the 1960s and 1970s? We argue that this attention to linkage should be understood in the context of research on chromosomal inversions and co-adapted gene complexes that occupied mid-century evolutionary genetics. For Lewontin, the complications of linkage were an extension of this chromosomal focus expressed in the new language of models for linkage disequilibrium.

Persistent seed banking as eco-evolutionary determinant of plant nucleotide diversity: novel population genetics insights.

Contents Summary 725 I. Introduction 725 II. Seed banks decrease the population extinction rate 726 III. Seed banks define the effective population size 727 IV. Seed banks affect the mutation rate 728 V. Seed banks affect the effective recombination rate 728 VI. Seed banks influence the rate and signatures of natural selection 729 VII. Conclusion 729 Acknowledgements 729 References 729 SUMMARY: Long-term persistent seed banking is a common temporal bet-hedging strategy in plants to adapt to unpredictable environments. The population genomics perspective developed in this article suggests that seed banking determines plant nucleotide diversity by decreasing the rate of genetic drift and the effect of linked selection while increasing mutational input. As a result, persistent seed banks are important factors determining the magnitude of the discrepancy between the census size of the above-ground plant population and its genetic diversity, an effect known as the Lewontin paradox. The theoretical population genetics predictions presented here can be tested by combining genome-wide polymorphism data with ecological studies of dormancy.

Is genomic diversity a useful proxy for census population size? Evidence from a species-rich community of desert lizards.

Species abundance data are critical for testing ecological theory, but obtaining accurate empirical estimates for many taxa is challenging. Proxies for species abundance can help researchers circumvent time and cost constraints that are prohibitive for long-term sampling. Under simple demographic models, genetic diversity is expected to correlate with census size, such that genome-wide heterozygosity may provide a surrogate measure of species abundance. We tested whether nucleotide diversity is correlated with long-term estimates of abundance, occupancy and degree of ecological specialization in a diverse lizard community from arid Australia. Using targeted sequence capture, we obtained estimates of genomic diversity from 30 species of lizards, recovering an average of 5,066 loci covering 3.6 Mb of DNA sequence per individual. We compared measures of individual heterozygosity to a metric of habitat specialization to investigate whether ecological preference exerts a measurable effect on genetic diversity. We find that heterozygosity is significantly correlated with species abundance and occupancy, but not habitat specialization. Demonstrating the power of genomic sampling, the correlation between heterozygosity and abundance/occupancy emerged from considering just one or two individuals per species. However, genetic diversity does no better at predicting abundance than a single day of traditional sampling in this community. We conclude that genetic diversity is a useful proxy for regional-scale species abundance and occupancy, but a large amount of unexplained variation in heterozygosity suggests additional constraints or a failure of ecological sampling to adequately capture variation in true population size.

Genetic diversity is largely unpredictable but scales with museum occurrences in a species-rich clade of Australian lizards.

Genetic diversity is a fundamental characteristic of species and is affected by many factors, including mutation rate, population size, life history and demography. To better understand the processes that influence levels of genetic diversity across taxa, we collected genome-wide restriction-associated DNA data from more than 500 individuals spanning 76 nominal species of Australian scincid lizards in the genus Ctenotus To avoid potential biases associated with variation in taxonomic practice across the group, we used coalescent-based species delimitation to delineate 83 species-level lineages within the genus for downstream analyses. We then used these genetic data to infer levels of within-population genetic diversity. Using a phylogenetically informed approach, we tested whether variation in genetic diversity could be explained by population size, environmental heterogeneity or historical demography. We find that the strongest predictor of genetic diversity is a novel proxy for census population size: the number of vouchered occurrences in museum databases. However, museum occurrences only explain a limited proportion of the variance in genetic diversity, suggesting that genetic diversity might be difficult to predict at shallower phylogenetic scales.

Natural Selection, Adaptive Topographies and the Problem of Statistical Inference: The Moraba scurra Controversy Under the Microscope.

This paper gives a detailed narrative of a controversial empirical research in postwar population genetics, the analysis of the cytological polymorphisms of an Australian grasshopper, Moraba scurra. This research intertwined key technical developments in three research areas during the 1950s and 1960s: it involved Dobzhansky's empirical research program on cytological polymorphisms, the mathematical theory of natural selection in two-locus systems, and the building of reliable estimates of natural selection in the wild. In the mid-1950s the cytologist Michael White discovered an interesting case of epistasis in populations of Moraba scurra. These observations received a wide diffusion when theoretical population geneticist Richard Lewontin represented White's data on adaptive topographies. These topographies connected the information on the genetic structure of these grasshopper populations with the formal framework of theoretical population genetics. As such, they appeared at the time as the most successful application of two-locus models of natural selection to an empirical study system. However, this connection generated paradoxical results: in the landscapes, all grasshopper populations were located on a ridge (an unstable equilibrium) while they were expected to reach a peak. This puzzling result fueled years of research and triggered a controversy attracting contributors from Australia, the United States and the United Kingdom. While the original problem seemed, at first, purely empirical, the subsequent controversy affected the main mathematical tools used in the study of two-gene systems under natural selection. Adaptive topographies and their underlying mathematical structure, Wright's mean fitness equations, were submitted to close scrutiny. Suspicion eventually shifted to the statistical machinery used in data analysis, reflecting the crucial role of statistical inference in applied population genetics. In the 1950s and 1960s, population geneticists were not simply in search for new generalizations about the evolutionary process and for new evidence about genetic variation: struggling against statistical artifacts, they were searching for reliable tests and descriptors to analyze their data.

"The Theory was Beautiful Indeed": Rise, Fall and Circulation of Maximizing Methods in Population Genetics (1930-1980).

Describing the theoretical population geneticists of the 1960s, Joseph Felsenstein reminisced: "our central obsession was finding out what function evolution would try to maximize. Population geneticists used to think, following Sewall Wright, that mean relative fitness, W, would be maximized by natural selection" (Felsenstein 2000). The present paper describes the genesis, diffusion and fall of this "obsession", by giving a biography of the mean fitness function in population genetics. This modeling method devised by Sewall Wright in the 1930s found its heyday in the late 1950s and early 1960s, in the wake of Motoo Kimura's and Richard Lewontin's works. It seemed a reliable guide in the mathematical study of deterministic effects (the study of natural selection in populations of infinite size, with no drift), leading to powerful generalizations presenting law-like properties. Progress in population genetics theory, it then seemed, would come from the application of this method to the study of systems with several genes. This ambition came to a halt in the context of the influential objections made by the Australian mathematician Patrick Moran in 1963. These objections triggered a controversy between mathematically- and biologically-inclined geneticists, with affected both the formal standards and the aims of population genetics as a science. Over the course of the 1960s, the mean fitness method withered with the ambition of developing the deterministic theory. The mathematical theory became increasingly complex. Kimura re-focused his modeling work on the theory of random processes; as a result of his computer simulations, Lewontin became the staunchest critic of maximizing principles in evolutionary biology. The mean fitness method then migrated to other research areas, being refashioned and used in evolutionary quantitative genetics and behavioral ecology.

Ascertaining regions affected by GC-biased gene conversion through weak-to-strong mutational hotspots.

A major objective for evolutionary biology is to identify regions affected by positive selection. High dN/dS values for proteins and accelerated lineage-specific substitution rates for non-coding regions are considered classic signatures of positive selection. However, these could also be the result of non-adaptive phenomena, such as GC-biased gene conversion (gBGC), which favors the fixation of strong (C/G) over weak (A/T) nucleotides. Recent estimates indicate that gBGC affected up to 20% of regions with signatures of positive selection. Here we evaluate the impact of gBGC through its molecular signature of weak-to-strong mutational hotspots. We implemented specific modifications to the test proposed by Tang and Lewontin (1999) for identifying regions of differential variability and applied it to regions previously investigated for the influence of gBGC. While we found significant agreement with previous reports, our results suggest a smaller influence of gBGC than previously estimated, warranting further development of methods for its detection.

Previously unmeasured genetic diversity explains part of Lewontin's paradox in a k -mer-based meta-analysis of 112 plant species.

At the molecular level, most evolution is expected to be neutral. A key prediction of this expectation is that the level of genetic diversity in a population should scale with population size. However, as was noted by Richard Lewontin in 1974 and reaffirmed by later studies, the slope of the population size-diversity relationship in nature is much weaker than expected under neutral theory. We hypothesize that one contributor to this paradox is that current methods relying on single nucleotide polymorphisms (SNPs) called from aligning short reads to a reference genome underestimate levels of genetic diversity in many species. To test this idea, we calculated nucleotide diversity ( π ) and k -mer-based metrics of genetic diversity across 112 plant species, amounting to over 205 terabases of DNA sequencing data from 27,488 individual plants. We then compared how these different metrics correlated with proxies of population size that account for both range size and population density variation across species. We found that our population size proxies scaled anywhere from about 3 to over 20 times faster with k -mer diversity than nucleotide diversity after adjusting for evolutionary history, mating system, life cycle habit, cultivation status, and invasiveness. The relationship between k -mer diversity and population size proxies also remains significant after correcting for genome size, whereas the analogous relationship for nucleotide diversity does not. These results suggest that variation not captured by common SNP-based analyses explains part of Lewontin's paradox in plants.

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?

Taylor's power law of fluctuation scaling and the growth-rate theorem.

Taylor's law (TL), a widely verified empirical relationship in ecology, states that the variance of population density is approximately a power-law function of mean density. The growth-rate theorem (GR) states that, in a subdivided population, the rate of change of the overall growth rate is proportional to the variance of the subpopulations' growth rates. We show that continuous-time exponential change implies GR at every time and, asymptotically for large time, TL with power-law exponent 2. We also show why diverse population-dynamic models predict TL in the limit of large time by identifying simple features these models share: If the mean population density and the variance of population density are (exactly or asymptotically) non-constant exponential functions of a parameter (e.g., time), then the variance of density is (exactly or asymptotically) a power-law function of mean density.

What Determines Levels of Mitochondrial Genetic Diversity in Birds?

What determines levels of genetic diversity in mitochondrial DNA remains unresolved. We have investigated the factors that are correlated to the level of synonymous diversity of mitochondrial DNA in more than 300 bird species. We find that diversity is significantly correlated to clutch and range size, but not significantly correlated to many other variables including body mass, latitude, and longevity. The correlation between diversity and range appears to be a consequence of a correlation between range and effective population size since a measure of the effectiveness of natural selection, which is expected to be correlated to the effective population size, is also correlated to range. The slope of the relationship between diversity and range is shallow, consistent with Lewontin's paradox, and very similar to the relationship found in mammals.

L-shaped distribution of the relative substitution rate (c/µ) observed for SARS-COV-2's genome, inconsistent with the selectionist theory, the neutral theory and the nearly neutral theory but a near-neutral balanced selection theory: Implication on "neutralist-selectionist" debate.

The genomic substitution rate (GSR) of SARS-CoV-2 exhibits a molecular clock feature and does not change under fluctuating environmental factors such as the infected human population (10°-107), vaccination etc. The molecular clock feature is believed to be inconsistent with the selectionist theory (ST). The GSR shows lack of dependence on the effective population size, suggesting Ohta's nearly neutral theory (ONNT) is not applicable to this virus. Big variation of the substitution rate within its genome is also inconsistent with Kimura's neutral theory (KNT). Thus, all three existing evolution theories fail to explain the evolutionary nature of this virus. In this paper, we proposed a Segment Substitution Rate Model (SSRM) under non-neutral selections and pointed out that a balanced mechanism between negative and positive selection of some segments that could also lead to the molecular clock feature. We named this hybrid mechanism as near-neutral balanced selection theory (NNBST) and examined if it was followed by SARS-CoV-2 using the three independent sets of SARS-CoV-2 genomes selected by the Nextstrain team. Intriguingly, the relative substitution rate of this virus exhibited an L-shaped probability distribution consisting with NNBST rather than Poisson distribution predicted by KNT or an asymmetric distribution predicted by ONNT in which nearly neutral sites are believed to be slightly deleterious only, or the distribution that is lack of nearly neutral sites predicted by ST. The time-dependence of the substitution rates for some segments and their correlation with the vaccination were observed, supporting NNBST. Our relative substitution rate method provides a tool to resolve the long standing "neutralist-selectionist" controversy. Implications of NNBST in resolving Lewontin's Paradox is also discussed.

The apportionment of citations: a scientometric analysis of Lewontin 1972.

'The apportionment of human diversity' (1972) is the most highly cited research article published by geneticist Richard Lewontin in his career. This study's primary result-that most genetic diversity in humans can be accounted for by within-population differences, not between-population differences-along with Lewontin's outspoken, politically charged interpretations thereof, has become foundational to the scientific and cultural discourse pertaining to human genetic variation. The article has an unusual bibliometric trajectory in that it is much more salient in the bibliographic record today compared to the first 20 years after its publication. Here, we highlight four factors that may have played a role in shaping the paper's fame: (i) citations in influential publications across several disciplines; (ii) Lewontin's own popular books and media appearances; (iii) the renaissance of population genetics research of the early 1990s; and (iv) the serendipitous collision of scientific progress, influential books and papers, and heated controversies around the year 1994. We conclude with an analysis of Twitter data to characterize the communities and conversations that continue to keep this study at the centre of discussions about race and genetics, prompting new challenges for scientists who have inherited Lewontin's legacy. This article is part of the theme issue 'Celebrating 50 years since Lewontin's apportionment of human diversity'.

Population differentiation of polygenic score predictions under stabilizing selection.

Given the many small-effect loci uncovered by genome-wide association studies (GWAS), polygenic scores have become central to genomic medicine, and have found application in diverse settings including evolutionary studies of adaptation. Despite their promise, polygenic scores have been found to suffer from limited portability across human populations. This at first seems in conflict with the observation that most common genetic variation is shared among populations. We investigate one potential cause of this discrepancy: stabilizing selection on complex traits. Counterintuitively, while stabilizing selection constrains phenotypic evolution, it accelerates the loss and fixation of alleles underlying trait variation within populations (GWAS loci). Thus even when populations share an optimum phenotype, stabilizing selection erodes the variance contributed by their shared GWAS loci, such that predictions from GWAS in one population explain less of the phenotypic variation in another. We develop theory to quantify how stabilizing selection is expected to reduce the prediction accuracy of polygenic scores in populations not represented in GWAS samples. In addition, we find that polygenic scores can substantially overstate average genetic differences of phenotypes among populations. We emphasize stabilizing selection around a common optimum as a useful null model to connect patterns of allele frequency and polygenic score differentiation. This article is part of the theme issue 'Celebrating 50 years since Lewontin's apportionment of human diversity'.

Forensic genetics through the lens of Lewontin: population structure, ancestry and race.

In his famous 1972 paper, Richard Lewontin used 'classical' protein-based markers to show that greater than 85% of human genetic diversity was contained within, rather than between, populations. At that time, these same markers also formed the basis of forensic technology aiming to identify individuals. This review describes the evolution of forensic genetic methods into DNA profiling, and how the field has accounted for the apportionment of genetic diversity in considering the weight of forensic evidence. When investigative databases fail to provide a match to a crime-scene profile, specific markers can be used to seek intelligence about a suspect: these include inferences on population of origin (biogeographic ancestry) and externally visible characteristics, chiefly pigmentation of skin, hair and eyes. In this endeavour, ancestry and phenotypic variation are closely entangled. The markers used show patterns of inter- and intrapopulation diversity that are very atypical compared to the genome as a whole, and reinforce an apparent link between ancestry and racial divergence that is not systematically present otherwise. Despite the legacy of Lewontin's result, therefore, in a major area in which genetics coincides with issues of public interest, methods tend to exaggerate human differences and could thereby contribute to the reification of biological race. This article is part of the theme issue 'Celebrating 50 years since Lewontin's apportionment of human diversity'.

Quantifying the relationship between genetic diversity and population size suggests natural selection cannot explain Lewontin's Paradox.

Neutral theory predicts that genetic diversity increases with population size, yet observed levels of diversity across metazoans vary only two orders of magnitude while population sizes vary over several. This unexpectedly narrow range of diversity is known as Lewontin's Paradox of Variation (1974). While some have suggested selection constrains diversity, tests of this hypothesis seem to fall short. Here, I revisit Lewontin's Paradox to assess whether current models of linked selection are capable of reducing diversity to this extent. To quantify the discrepancy between pairwise diversity and census population sizes across species, I combine previously-published estimates of pairwise diversity from 172 metazoan taxa with newly derived estimates of census sizes. Using phylogenetic comparative methods, I show this relationship is significant accounting for phylogeny, but with high phylogenetic signal and evidence that some lineages experience shifts in the evolutionary rate of diversity deep in the past. Additionally, I find a negative relationship between recombination map length and census size, suggesting abundant species have less recombination and experience greater reductions in diversity due to linked selection. However, I show that even assuming strong and abundant selection, models of linked selection are unlikely to explain the observed relationship between diversity and census sizes across species.

Demography and Natural Selection Have Shaped Genetic Variation in the Widely Distributed Conifer Norway Spruce (Picea abies).

Under the neutral theory, species with larger effective population size are expected to harbor higher genetic diversity. However, across a wide variety of organisms, the range of genetic diversity is orders of magnitude more narrow than the range of effective population size. This observation has become known as Lewontin's paradox and although aspects of this phenomenon have been extensively studied, the underlying causes for the paradox remain unclear. Norway spruce (Picea abies) is a widely distributed conifer species across the northern hemisphere, and it consequently plays a major role in European forestry. Here, we use whole-genome resequencing data from 35 individuals to perform population genomic analyses in P. abies in an effort to understand what drives genome-wide patterns of variation in this species. Despite having a very wide geographic distribution and an corresponding enormous current population size, our analyses find that genetic diversity of P. abies is low across a number of populations (π = 0.0049 in Central-Europe, π = 0.0063 in Sweden-Norway, π = 0.0063 in Finland). To assess the reasons for the low levels of genetic diversity, we infer the demographic history of the species and find that it is characterized by several reoccurring bottlenecks with concomitant decreases in effective population size can, at least partly, provide an explanation for low polymorphism we observe in P. abies. Further analyses suggest that recurrent natural selection, both purifying and positive selection, can also contribute to the loss of genetic diversity in Norway spruce by reducing genetic diversity at linked sites. Finally, the overall low mutation rates seen in conifers can also help explain the low genetic diversity maintained in Norway spruce.

Evolution of Mutation Rate in Astronomically Large Phytoplankton Populations.

Genetic diversity is expected to be proportional to population size, yet, there is a well-known, but unexplained lack of genetic diversity in large populations-the "Lewontin's paradox." Larger populations are expected to evolve lower mutation rates, which may help to explain this paradox. Here, we test this conjecture by measuring the spontaneous mutation rate in a ubiquitous unicellular marine phytoplankton species Emiliania huxleyi (Haptophyta) that has modest genetic diversity despite an astronomically large population size. Genome sequencing of E. huxleyi mutation accumulation lines revealed 455 mutations, with an unusual GC-biased mutation spectrum. This yielded an estimate of the per site mutation rate µ = 5.55×10-10 (CI 95%: 5.05×10-10 - 6.09×10-10), which corresponds to an effective population size Ne ∼ 2.7×106. Such a modest Ne is surprising for a ubiquitous and abundant species that accounts for up to 10% of global primary productivity in the oceans. Our results indicate that even exceptionally large populations do not evolve mutation rates lower than ∼10-10 per nucleotide per cell division. Consequently, the extreme disparity between modest genetic diversity and astronomically large population size in the plankton species cannot be explained by an unusually low mutation rate.