Reconstructing the History of Variation in Effective Population Size along Phylogenies

Genome and life-history evolution link bird diversification to the end-Cretaceous mass extinction

Complex patterns of genome evolution associated with the end-Cretaceous [Cretaceous-Paleogene (K-Pg)] mass extinction limit our understanding of the early evolutionary history of modern birds. Here, we analyzed patterns of avian molecular evolution and identified distinct macroevolutionary regimes across exons, introns, untranslated regions, and mitochondrial genomes. Bird clades originating near the K-Pg boundary exhibited numerous shifts in the mode of molecular evolution, suggesting a burst of genomic heterogeneity at this point in Earth's history. These inferred shifts in substitution patterns were closely related to evolutionary shifts in developmental mode, adult body mass, and patterns of metabolic scaling. Our results suggest that the end-Cretaceous mass extinction triggered integrated patterns of evolution across avian genomes, physiology, and life history near the dawn of the modern bird radiation.

The genetic status and rescue measure for a geographically isolated population of Amur tigers

The Amur tiger is currently confronted with challenges of anthropogenic development, leading to its population becoming fragmented into two geographically isolated groups: smaller and larger ones. Small and isolated populations frequently face a greater extinction risk, yet the small tiger population's genetic status and survival potential have not been assessed. Here, a total of 210 samples of suspected Amur tiger feces were collected from this small population, and the genetic background and population survival potentials were assessed by using 14 microsatellite loci. Our results demonstrated that the mean number of alleles in all loci was 3.7 and expected heterozygosity was 0.6, indicating a comparatively lower level of population genetic diversity compared to previously reported studies on other subspecies. The genetic estimates of effective population size (Ne) and the Ne/N ratio were merely 7.6 and 0.152, respectively, representing lower values in comparison to the Amur tiger population in Sikhote-Alin (the larger group). However, multiple methods have indicated the possibility of genetic divergence within our isolated population under study. Meanwhile, the maximum kinship recorded was 0.441, and the mean inbreeding coefficient stood at 0.0868, both of which are higher than those observed in other endangered species, such as the African lion and the grey wolf. Additionally, we have identified a significant risk of future extinction if the lethal equivalents were to reach 6.26, which is higher than that of other large carnivores. Further, our simulation results indicated that an increase in the number of breeding females would enhance the prospects of this population. In summary, our findings provide a critical theoretical basis for further bailout strategies concerning Amur tigers.

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.

Rules of teeth development align microevolution with macroevolution in extant and extinct primates

Macroevolutionary biologists have classically rejected the notion that higher-level patterns of divergence arise through microevolutionary processes acting within populations. For morphology, this consensus partly derives from the inability of quantitative genetics models to correctly predict the behaviour of evolutionary processes at the scale of millions of years. Developmental studies (evo-devo) have been proposed to reconcile micro- and macroevolution. However, there has been little progress in establishing a formal framework to apply evo-devo models of phenotypic diversification. Here we reframe this issue by asking whether using evo-devo models to quantify biological variation can improve the explanatory power of comparative models, thus helping us bridge the gap between micro- and macroevolution. We test this prediction by evaluating the evolution of primate lower molars in a comprehensive dataset densely sampled across living and extinct taxa. Our results suggest that biologically informed morphospaces alongside quantitative genetics models allow a seamless transition between the micro- and macroscales, whereas biologically uninformed spaces do not. We show that the adaptive landscape for primate teeth is corridor like, with changes in morphology within the corridor being nearly neutral. Overall, our framework provides a basis for integrating evo-devo into the modern synthesis, allowing an operational way to evaluate the ultimate causes of macroevolution.

Strong Purifying Selection in Haploid Tissue-Specific Genes of Scots Pine Supports the Masking Theory

The masking theory states that genes expressed in a haploid stage will be under more efficient selection. In contrast, selection will be less efficient in genes expressed in a diploid stage, where the fitness effects of recessive deleterious or beneficial mutations can be hidden from selection in heterozygous form. This difference can influence several evolutionary processes such as the maintenance of genetic variation, adaptation rate, and genetic load. Masking theory expectations have been confirmed in single-cell haploid and diploid organisms. However, in multicellular organisms, such as plants, the effects of haploid selection are not clear-cut. In plants, the great majority of studies indicating haploid selection have been carried out using male haploid tissues in angiosperms. Hence, evidence in these systems is confounded with the effects of sexual selection and intraspecific competition. Evidence from other plant groups is scarce, and results show no support for the masking theory. Here, we have used a gymnosperm Scots pine megagametophyte, a maternally derived seed haploid tissue, and four diploid tissues to test the strength of purifying selection on a set of genes with tissue-specific expression. By using targeted resequencing data of those genes, we obtained estimates of genetic diversity, the site frequency spectrum of 0-fold and 4-fold sites, and inferred the distribution of fitness effects of new mutations in haploid and diploid tissue-specific genes. Our results show that purifying selection is stronger for tissue-specific genes expressed in the haploid megagametophyte tissue and that this signal of strong selection is not an artifact driven by high expression levels.

What does the distribution of fitness effects of new mutations reflect? Insights from plants

The distribution of fitness effects (DFE) of new mutations plays a central role in molecular evolution. It is therefore crucial to be able to estimate it accurately from genomic data and to understand the factors that shape it. After a rapid overview of available methods to characterize the fitness effects of mutations, we review what is known on the factors affecting them in plants. Available data indicate that life history traits (e.g. mating system and longevity) have a major effect on the DFE. By contrast, the impact of demography within species appears to be more limited. These results remain to be confirmed, and methods to estimate the joint evolution of demography, life history traits, and the DFE need to be developed.

OUP accepted manuscript

The nearly neutral theory is a common framework to describe natural selection at the molecular level. This theory emphasizes the importance of slightly deleterious mutations by recognizing their ability to segregate and eventually get fixed due to genetic drift in spite of the presence of purifying selection. As genetic drift is stronger in smaller than in larger populations, a correlation between population size and molecular measures of natural selection is expected within the nearly neutral theory. However, this hypothesis was originally formulated under equilibrium conditions. As most natural populations are not in equilibrium, testing the relationship empirically may lead to confounded outcomes. Demographic nonequilibria, for instance following a change in population size, are common scenarios that are expected to push the selection–drift relationship off equilibrium. By explicitly modeling the effects of a change in population size on allele frequency trajectories in the Poisson random field framework, we obtain analytical solutions of the nonstationary allele frequency spectrum. This enables us to derive exact results of measures of natural selection and effective population size in a demographic nonequilibrium. The study of their time-dependent relationship reveals a substantial deviation from the equilibrium selection–drift balance after a change in population size. Moreover, we show that the deviation is sensitive to the combination of different measures. These results therefore constitute relevant tools for empirical studies to choose suitable measures for investigating the selection–drift relationship in natural populations. Additionally, our new modeling approach extends existing population genetics theory and can serve as foundation for methodological developments.

Quantifying the impact of changes in effective population size and expression level on the rate of coding sequence evolution

Molecular sequences are shaped by selection, where the strength of selection relative to drift is determined by effective population size (N_e). Populations with high N_e are expected to undergo stronger purifying selection, and consequently to show a lower substitution rate for selected mutations relative to the substitution rate for neutral mutations (ω). However, computational models based on biophysics of protein stability have suggested that ω can also be independent of N_e. Together, the response of ω to changes in N_e depends on the specific mapping from sequence to fitness. Importantly, an increase in protein expression level has been found empirically to result in decrease of ω, an observation predicted by theoretical models assuming selection for protein stability. Here, we derive a theoretical approximation for the response of ω to changes in N_e and expression level, under an explicit genotype-phenotype-fitness map. The method is generally valid for additive traits and log-concave fitness functions. We applied these results to protein undergoing selection for their conformational stability and corroborate out findings with simulations under more complex models. We predict a weak response of ω to changes in either N_e or expression level, which are interchangeable. Based on empirical data, we propose that fitness based on the conformational stability may not be a sufficient mechanism to explain the empirically observed variation in ω across species. Other aspects of protein biophysics might be explored, such as protein-protein interactions, which can lead to a stronger response of ω to changes in N_e.