Simple Population Genetic Models Research Articles

The grapheme-valued Wright–Fisher diffusion with mutation

In Athreya et al. (2021), models from population genetics were used to define stochastic dynamics in the space of graphons arising as continuum limits of dense graphs. In the present paper we exhibit an example of a simple neutral population genetics model for which this dynamics is a Markovian diffusion that can be characterized as the solution of a martingale problem. In particular, we consider a Markov chain in the space of finite graphs that resembles a Moran model with resampling and mutation. We encode the finite graphs as graphemes, which can be represented as a triple consisting of a vertex set (or more generally, a topological space), an adjacency matrix, and a sampling (Borel) measure. We equip the space of graphons with convergence of sample subgraph densities and show that the grapheme-valued Markov chain converges to a grapheme-valued diffusion as the number of vertices goes to infinity. We show that the grapheme-valued diffusion has a stationary distribution that is linked to the Griffiths–Engen–McCloskey (GEM) distribution. In a companion paper (Greven et al. 2023), we build up a general theory for obtaining grapheme-valued diffusions via genealogies of models in population genetics.

Shared Evolutionary Path in Social Microbiomes.

Social networks can influence the ecology of gut bacteria, shaping the species composition of the gut microbiome in humans and other animals. Gut commensals evolve and can adapt at a rapid pace when colonizing healthy hosts. Here, we aimed at assessing the impact of host-to-host bacterial transmission on Escherichia coli evolution in the mammalian gut. Using an in vivo experimental evolution approach in mice, we found a transmission rate of 7% (±3% 2× standard error [2SE]) of E. coli cells per day between hosts inhabiting the same household. Consistent with the predictions of a simple population genetics model of mutation-selection-migration, the level of shared events resulting from within host evolution is greatly enhanced in cohoused mice, showing that hosts undergoing the same diet and habit are not only expected to have similar microbiome species compositions but also similar microbiome evolutionary dynamics. Furthermore, we estimated the rate of mutation accumulation of E. coli to be 3.0 × 10-3 (±0.8 × 10-3 2SE) mutations/genome/generation, irrespective of the social context of the regime. Our results reveal the impact of bacterial migration across hosts in shaping the adaptive evolution of new strains colonizing gut microbiomes.

How much does the unguarded X contribute to sex differences in life span?

Females and males often have markedly different mortality rates and life spans, but it is unclear why these forms of sexual dimorphism evolve. The unguarded X hypothesis contends that dimorphic life spans arise from sex differences in X or Z chromosome copy number (i.e., one copy in the “heterogametic” sex; two copies in the “homogametic” sex), which leads to a disproportionate expression of deleterious mutations by the heterogametic sex (e.g., mammalian males; avian females). Although data on adult sex ratios and sex‐specific longevity are consistent with predictions of the unguarded X hypothesis, direct experimental evidence remains scant, and alternative explanations are difficult to rule out. Using a simple population genetic model, we show that the unguarded X effect on sex differential mortality is a function of several reasonably well‐studied evolutionary parameters, including the proportion of the genome that is sex linked, the genomic deleterious mutation rate, the mean dominance of deleterious mutations, the relative rates of mutation and strengths of selection in each sex, and the average effect of mutations on survival and longevity relative to their effects on fitness. We review published estimates of these parameters, parameterize our model with them, and show that unguarded X effects are too small to explain observed sex differences in life span across species. For example, sex differences in mean life span are known to often exceed 20% (e.g., in mammals), whereas our parameterized models predict unguarded X effects of a few percent (e.g., 1–3% in Drosophila and mammals). Indeed, these predicted unguarded X effects fall below statistical thresholds of detectability in most experiments, potentially explaining why direct tests of the hypothesis have generated little support for it. Our results suggest that evolution of sexually dimorphic life spans is predominantly attributable to other mechanisms, potentially including “toxic Y” effects and sexual dimorphism for optimal investment in survival versus reproduction.

Dynamics of reduced genetic diversity in increasingly fragmented populations of Florida scrub jays,Aphelocoma coerulescens.

Understanding the genomic consequences of population decline is important for predicting species' vulnerability to intensifying global change. Empirical information about genomic changes in populations in the early stages of decline, especially for those still experiencing immigration, remains scarce. We used 7834 autosomal SNPs and demographic data for 288 Florida scrub jays (Aphelocoma coerulescens; FSJ) sampled in 2000 and 2008 to compare levels of genetic diversity, inbreeding, relatedness, and lengths of runs of homozygosity (ROH) between two subpopulations within dispersal distance of one another but have experienced contrasting demographic trajectories. At Archbold Biological Station (ABS), the FSJ population has been stable because of consistent habitat protection and management, while at nearby Placid Lakes Estates (PLE), the population declined precipitously due to suburban development. By the onset of our sampling in 2000, birds in PLE were already less heterozygous, more inbred, and on average more related than birds in ABS. No significant changes occurred in heterozygosity or inbreeding across the 8‐year sampling interval, but average relatedness among individuals decreased in PLE, thus by 2008 average relatedness did not differ between sites. PLE harbored a similar proportion of short ROH but a greater proportion of long ROH than ABS, suggesting one continuous population of shared demographic history in the past, which is now experiencing more recent inbreeding. These results broadly uphold the predictions of simple population genetic models based on inferred effective population sizes and rates of immigration. Our study highlights how, in just a few generations, formerly continuous populations can diverge in heterozygosity and levels of inbreeding with severe local population decline despite ongoing gene flow.

Competition dynamics in long-term propagations of Schizosaccharomyces pombe strain communities.

Experimental evolution studies with microorganisms such as bacteria and yeast have been an increasingly important and powerful tool to draw long‐term inferences of how microbes interact. However, while several strains of the same species often exist in natural environments, many ecology and evolution studies in microbes are typically performed with isogenic populations of bacteria or yeast. In the present study, we firstly perform a genotypic and phenotypic characterization of two laboratory and eight natural strains of the yeast Schizosaccharomyces pombe. We then propagated, in a rich resource environment, yeast communities of 2, 3, 4, and 5 strains for hundreds of generations and asked which fitness‐related phenotypes—maximum growth rate or relative competitive fitness—would better predict the outcome of a focal strain during the propagations. While the strain's growth rates would wrongly predict long‐term coexistence, pairwise competitive fitness with a focal strain qualitatively predicted the success or extinction of the focal strain by a simple multigenotype population genetics model, given the initial community composition. Interestingly, we have also measured the competitive fitness of the ancestral and evolved communities by the end of the experiment (≈370 generations) and observed frequent maladaptation to the abiotic environment in communities with more than three members. Overall, our results aid establishing pairwise competitive fitness as good qualitative measurement of long‐term community composition but also reveal a complex adaptive scenario when trying to predict the evolutionary outcome of those communities.

PiRNA Clusters Need a Minimum Size to Control Transposable Element Invasions.

piRNA clusters are thought to repress transposable element (TE) activity in mammals and invertebrates. Here, we show that a simple population genetics model reveals a constraint on the size of piRNA clusters: The total size of the piRNA clusters of an organism must exceed 0.2% of a genome to repress TE invasions. Moreover, larger piRNA clusters accounting for up to 3% of the genome may be necessary when populations are small, transposition rates are high, and TE insertions are recessive. If piRNA clusters are too small, the load of deleterious TE insertions that accumulate during a TE invasion may drive populations extinct before an effective piRNA-based defense against the TE can be established. Our findings are solely based on three well-supported assumptions: 1) TEs multiply within genomes, 2) TEs are mostly deleterious, and 3) piRNA clusters act as transposon traps, where a single insertion in a cluster silences all TE copies in trans. Interestingly, the piRNA clusters of some species meet our observed minimum size requirements, whereas the clusters of other species do not. Species with small piRNA clusters, such as humans and mice, may experience severe fitness reductions during invasions of novel TEs, which is possibly even threatening the persistence of some populations. This work also raises the important question of how piRNA clusters evolve. We propose that the size of piRNA clusters may be at an equilibrium between evolutionary forces that act to expand and contract piRNA clusters.

Effects of mucus trail following on the distance between individuals of opposite sex and its influence on the evolution of the trait in the Ezo abalone Haliotis discus hannai.

BackgroundAggregation affects the fertilization rate of species that utilize external fertilization. However, the process of aggregation has not been studied in detail, using either theoretical models or real world observations. We used the Ezo abalone Haliotis discus hannai as a model animal species to evaluate whether mucus trail following (MTF) facilitates aggregation by reducing the distance between members of opposite sexes. We also examined whether the fertilization rate increase with mucus-trail-following is an evolutionary driving force in mucus following behavior.MethodsWe used a y-maze to test whether H. discus hannai follows the mucus trails of other individuals. Distances between members of the opposite sex of MTF individuals were compared to non-MTF individuals using an individual-based model (IBM) consistent with the behavior of H. discus hannai. To examine whether MTF behavior evolved to reduce distances between members of the opposite sex, we constructed simple population genetic models of a diploid population with nonoverlapping, discrete generations.ResultsHaliotis discus hannai chose the y-maze arm with the mucus trail more frequently than the one without, regardless of the sex of the abalone that secreted the mucus or the reproductive season. In the IBM the distance between opposite sexes was significantly reduced by MTF behavior; however, the difference in distances between opposite sex compared to same sex individuals was only several centimeters. Simple population genetic models indicated that the aggregating effect of MTF between the opposite sex members could be an evolutionary driving force.ConclusionsThese results suggest that observed MTF behavior might have evolved as a mechanism to increase the fertilization rates although other factors could also be involved.

Imprinting sets the stage for speciation.

Sexual imprinting-a phenomenon in which offspring learn parental traits and later use them as a model for their own mate preferences-can generate reproductive barriers between species1. When the target of imprinting is a mating trait that differs among young lineages, imprinted preferences may contribute to behavioural isolation and facilitate speciation1,2. However, in most models of speciation by sexual selection, divergent natural selection is also required; the latter acts to generate and maintain variation in the sexually selected trait or traits, and in the mating preferences that act upon them3. Here we demonstrate that imprinting, in addition to mediating female mate preferences, can shape biases in male-male aggression. These biases can act similarly to natural selection to maintain variation in traits and mate preferences, which facilitates reproductive isolation driven entirely by sexual selection. Using a cross-fostering study, we show that both male and female strawberry poison frogs (Oophaga pumilio) imprint on coloration, which is a mating trait that has diverged recently and rapidly in this species4. Cross-fostered females prefer to court mates of thesame colour as their foster mother, and cross-fostered males are more aggressive towards rivals that share the colour of their foster mother. We also use a simple population-genetics model to demonstrate that when both male aggression biases and female mate preferences are formed through parental imprinting, sexual selection alone can (1) stabilize a sympatric polymorphism and (2) strengthen the trait-preference association that leads to behavioural reproductive isolation. Our study provides evidence of imprinting in an amphibian and suggests that this rarely considered combination of rival and sexual imprinting can reduce gene flow between individuals that bear divergent mating traits, which sets the stage for speciation by sexual selection.

Self-limiting population genetic control with sex-linked genome editors

In male heterogametic species the Y chromosome is transmitted solely from fathers to sons, and is selected for based only on its impacts on male fitness. This fact can be exploited to develop efficient pest control strategies that use Y-linked editors to disrupt the fitness of female descendants. With simple population genetic and dynamic models we show that Y-linked editors can be substantially more efficient than other self-limiting strategies and, while not as efficient as gene drive approaches, are expected to have less impact on non-target populations with which there is some gene flow. Efficiency can be further augmented by simultaneously releasing an autosomal X-shredder construct, in either the same or different males. Y-linked editors may be an attractive option to consider when efficient control of a species is desired in some locales but not others.

Strong hybrid male incompatibilities impede the spread of a selfish chromosome between populations of a fly.

Meiotically driving sex chromosomes manipulate gametogenesis to increase their transmission at a cost to the rest of the genome. The intragenomic conflicts they produce have major impacts on the ecology and evolution of their host species. However, their ecological dynamics remain poorly understood. Simple population genetic models predict meiotic drivers will rapidly reach fixation in populations and spread across landscapes. In contrast, natural populations commonly show spatial variation in the frequency of drivers, with drive present in clines or mosaics across species ranges. For example, Drosophila subobscura harbors a sex ratio distorting drive chromosome (SRs) at 15–25% frequency in North Africa, present at less than 2% frequency in adjacent southern Spain, and absent in other European populations. Here, we investigate the forces preventing the spread of the driver northward. We show that SRs has remained at a constant frequency in North Africa, and failed to spread in Spain. We find strong evidence that spread is impeded by genetic incompatibility between SRs and Spanish autosomal backgrounds. When we cross SRs from North Africa onto Spanish genetic backgrounds we observe strong incompatibilities specific to hybrids bearing SRs. The incompatibilities increase in severity in F2 male hybrids, leading to almost complete infertility. We find no evidence supporting an alternative hypothesis, that there is resistance to drive in Spanish populations. We conclude that the source of the stepped frequency variation is genetic incompatibility between the SRs chromosome and the genetic backgrounds of the adjacent population, preventing SRs spreading northward. The low frequency of SRs in South Spain is consistent with recurrent gene flow across the Strait of Gibraltar combined with selection against the SRs element through genetic incompatibility. This demonstrates that incompatibilities between drive chromosomes and naïve populations can prevent the spread of drive between populations, at a continental scale.

Selection on non-antigenic gene segments of seasonal influenza A virus and its impact on adaptive evolution.

Most studies on seasonal influenza A/H3N2 virus adaptation have focused on the main antigenic gene, hemagglutinin. However, there is increasing evidence that the genome-wide genetic background of novel antigenic variants can influence these variants’ emergence probabilities and impact their patterns of dominance in the population. This suggests that non-antigenic genes may be important in shaping the viral evolutionary dynamics. To better understand the role of selection on non-antigenic genes in the adaptive evolution of seasonal influenza viruses, we have developed a simple population genetic model that considers a virus with one antigenic and one non-antigenic gene segment. By simulating this model under different regimes of selection and reassortment, we find that the empirical patterns of lineage turnover for the antigenic and non-antigenic gene segments are best captured when there is both limited viral coinfection and selection operating on both gene segments. In contrast, under a scenario of only neutral evolution in the non-antigenic gene segment, we see persistence of multiple lineages for long periods of time in that segment, which is not compatible with observed molecular evolutionary patterns. Further, we find that reassortment, occurring in coinfected individuals, can increase the speed of viral adaptive evolution by primarily reducing selective interference and genetic linkage effects. Together, these findings suggest that, for influenza, with six internal or non-antigenic gene segments, the evolutionary dynamics of novel antigenic variants are likely to be influenced by the genome-wide genetic background as a result of linked selection among both beneficial and deleterious mutations.

When maladaptive gene flow does not increase selection.

Populations receiving high maladaptive gene flow are expected to experience strong directional selection-because gene flow pulls mean phenotypes away from local fitness peaks. We tested this prediction by means of a large and replicated mark-recapture study of threespine stickleback (Gasterosteus aculeatus) in two stream populations. One of the populations (outlet) experiences high gene flow from the lake population and its morphology is correspondingly poorly adapted. The other population (inlet) experiences very low gene flow from the lake population and its morphology is correspondingly well adapted. Contrary to the above prediction, selection was not stronger in the outlet than in the inlet, a result that forced us to consider potential reasons for why maladaptive gene flow might not increase selection. Of particular interest, we show by means of a simple population genetic model that maladaptive gene flow can-under reasonable conditions-decrease the strength of directional selection. This outcome occurs when immigrants decrease mean fitness in the resident population, which decreases the strength of selection against maladapted phenotypes. We argue that this previously unrecognized effect of gene flow deserves further attention in theoretical and empirical studies.

Identical twins in forensic genetics — Epidemiology and risk based estimation of weight of evidence

The increase in the number of forensic genetic loci used for identification purposes results in infinitesimal random match probabilities. These probabilities are computed under assumptions made for rather simple population genetic models. Often, the forensic expert reports likelihood ratios, where the alternative hypothesis is assumed not to encompass close relatives. However, this approach implies that important factors present in real human populations are discarded. This approach may be very unfavourable to the defendant.In this paper, we discuss some important aspects concerning the closest familial relationship, i.e., identical (monozygotic) twins, when reporting the weight of evidence. This can be done even when the suspect has no knowledge of an identical twin or when official records hold no twin information about the suspect. The derived expressions are not original as several authors previously have published results accounting for close familial relationships. However, we revisit the discussion to increase the awareness among forensic genetic practitioners and include new information on medical and societal factors to assess the risk of not considering a monozygotic twin as the true perpetrator.If accounting for a monozygotic twin in the weight of evidence, it implies that the likelihood ratio is truncated at a maximal value depending on the prevalence of monozygotic twins and the societal efficiency of recognising a monozygotic twin. If a monozygotic twin is considered as an alternative proposition, then data relevant for the Danish society suggests that the threshold of likelihood ratios should approximately be between 150,000 and 2,000,000 in order to take the risk of an unrecognised identical, monozygotic twin into consideration. In other societies, the threshold of the likelihood ratio in crime cases may reach other, often lower, values depending on the recognition of monozygotic twins and the age of the suspect. In general, more strictly kept registries will imply larger thresholds on the likelihood ratio as the monozygotic twin explanation gets less probable.

Biased introgression of mitochondrial and nuclear genes: a comparison of diploid and haplodiploid systems

Hybridization between recently diverged species, even if infrequent, can lead to the introgression of genes from one species into another. The rates of mitochondrial and nuclear introgression often differ, with some taxa showing biases for mitochondrial introgression and others for nuclear introgression. Several hypotheses exist to explain such biases, including adaptive introgression, sex differences in dispersal rates, sex-specific prezygotic isolation and sex-specific fitness of hybrids (e.g. Haldane's rule). We derive a simple population genetic model that permits an analysis of sex-specific demographic and fitness parameters and measures the relative rates of mitochondrial and nuclear introgression between hybridizing pairs. We do this separately for diploid and haplodiploid species. For diploid taxa, we recover results consistent with previous hypotheses: an excess of one sex among the hybridizing migrants or sex-specific prezygotic isolation causes a bias for one type of marker or the other; when Haldane's rule is obeyed, we find a mitochondrial bias in XY systems and a nuclear bias in ZW systems. For haplodiploid taxa, the model reveals that owing to their unique transmission genetics, they are seemingly assured of strong mitochondrial biases in introgression rates, unlike diploid taxa, where the relative fitness of male and female hybrids can tip the bias in either direction. This heretofore overlooked aspect of hybridization in haplodiploids provides what is perhaps the most likely explanation for differential introgression of mitochondrial and nuclear markers and raises concerns about the use of mitochondrial DNA barcodes for species delimitation in these taxa.

Signatures of positive selection: from selective sweeps at individual loci to subtle allele frequency changes in polygenic adaptation.

In the past 15 years, numerous methods have been developed to detect selective sweeps underlying adaptations. These methods are based on relatively simple population genetic models, including one or two loci at which positive directional selection occurs, and one or two marker loci at which the impact of selection on linked neutral variation is quantified. Information about the phenotype under selection is not included in these models (except for fitness). In contrast, in the quantitative genetic models of adaptation, selection acts on one or more phenotypic traits, such that a genotype-phenotype map is required to bridge the gap to population genetics theory. Here I describe the range of population genetic models from selective sweeps in a panmictic population of constant size to evolutionary traffic when simultaneous sweeps at multiple loci interfere, and I also consider the case of polygenic selection characterized by subtle allele frequency shifts at many loci. Furthermore, I present an overview of the statistical tests that have been proposed based on these population genetics models to detect evidence for positive selection in the genome.

Causes of natural variation in fitness: Evidence from studies of Drosophila populations

DNA sequencing has revealed high levels of variability within most species. Statistical methods based on population genetics theory have been applied to the resulting data and suggest that most mutations affecting functionally important sequences are deleterious but subject to very weak selection. Quantitative genetic studies have provided information on the extent of genetic variation within populations in traits related to fitness and the rate at which variability in these traits arises by mutation. This paper attempts to combine the available information from applications of the two approaches to populations of the fruitfly Drosophila in order to estimate some important parameters of genetic variation, using a simple population genetics model of mutational effects on fitness components. Analyses based on this model suggest the existence of a class of mutations with much larger fitness effects than those inferred from sequence variability and that contribute most of the standing variation in fitness within a population caused by the input of mildly deleterious mutations. However, deleterious mutations explain only part of this standing variation, and other processes such as balancing selection appear to make a large contribution to genetic variation in fitness components in Drosophila.

Mating advantage for rare males in wild guppy populations

To understand the processes that maintain genetic diversity is a long-standing challenge in evolutionary biology, with implications for predicting disease resistance, response to environmental change, and population persistence. Simple population genetic models are not sufficient to explain the high levels of genetic diversity sometimes observed in ecologically important traits. In guppies (Poecilia reticulata), male colour pattern is both diverse and heritable, and is arguably one of the most extreme examples of morphological polymorphism known. Negative frequency-dependent selection (NFDS), a form of selection in which genotypes are favoured when they are rare, can potentially maintain such extensive polymorphism, but few experimental studies have confirmed its operation in nature. Here we use highly replicated experimental manipulations of natural populations to show that males with rare colour patterns have higher reproductive fitness, demonstrating NFDS mediated by sexual selection. Rare males acquired more mates and sired more offspring compared to common males and, as previously reported, had higher rates of survival. Orange colour, implicated in other studies of sexual selection in guppies, did predict male reproductive success, but only in one of three populations. These data support the hypothesis that NFDS maintains diversity in the colour patterns of male guppies through two selective agents, mates and predators. Similar field-based manipulations of genotype frequencies could provide a powerful approach to reveal the underlying ecological and behavioural mechanisms that maintain genetic and phenotypic diversity.

Implications of sex-specific selection for the genetic basis of disease.

Mutation and selection are thought to shape the underlying genetic basis of many common human diseases. However, both processes depend on the context in which they occur, such as environment, genetic background, or sex. Sex has widely known effects on phenotypic expression of genotype, but an analysis of how it influences the evolutionary dynamics of disease-causing variants has not yet been explored. We develop a simple population genetic model of disease susceptibility and evaluate it using a biologically plausible empirically based distribution of fitness effects among contributing mutations. The model predicts that alleles under sex-differential selection, including sexually antagonistic alleles, will disproportionately contribute to genetic variation for disease predisposition, thereby generating substantial sexual dimorphism in the genetic architecture of complex (polygenic) diseases. This is because such alleles evolve into higher population frequencies for a given effect size, relative to alleles experiencing equally strong purifying selection in both sexes. Our results provide a theoretical justification for expecting a sexually dimorphic genetic basis for variation in complex traits such as disease. Moreover, they suggest that such dimorphism is interesting – not merely something to control for – because it reflects the action of natural selection in molding the evolution of common disease phenotypes.

The effects of inheritance in tetraploids on genetic diversity and population divergence

Polyploids are traditionally classified into allopolyploids and autopolyploids, based on their evolutionary origin and their disomic or multisomic mode of inheritance. Over the past decade it has become increasingly clear that there is a continuum between disomic and multisomic inheritance, with the rate of tetrasomy differing among species and among chromosomes within species. Here, we use a simple population genetic model to study the impact of the mode of inheritance on the genetic diversity and population divergence of tetraploids. We found that under almost strict disomic inheritance the tetraploid genome is divided into two separate subgenomes, such as found in classical allopolyploids. In those cases, assuming full tetrasomy in the analysis of polyploid genetic data will lead to an important bias in estimates of genetic diversity and population divergence. However, we found that even a low rate of allele exchange between the two subgenomes, at about one event per generation, is sufficient to homogenise the allele frequencies over the subgenomes, and the estimates become essentially unbiased. The inbreeding coefficient F(IS) can then be used to detect whether the estimates of diversity and divergence will be biased when full multisomy is assumed. Finally, we found that different summary statistics for measuring the strength of population differentiation are differentially affected by a deviation from full tetrasomy. Our model results provide several useful guidelines for the analysis of polyploid data, helping researchers to determine when their inferences are biased and which summary statistics to use.

Phylogenetic Patterns of GC-Biased Gene Conversion in Placental Mammals and the Evolutionary Dynamics of Recombination Landscapes

GC-biased gene conversion (gBGC) is a major evolutionary force shaping genomic nucleotide landscapes, distorting the estimation of the strength of selection, and having potentially deleterious effects on genome-wide fitness. Yet, a global quantitative picture, at large evolutionary scale, of the relative strength of gBGC compared with selection and random drift is still lacking. Furthermore, owing to its dependence on the local recombination rate, gBGC results in modulations of the substitution patterns along genomes and across time which, if correctly interpreted, may yield quantitative insights into the long-term evolutionary dynamics of recombination landscapes. Deriving a model of the substitution process at putatively neutral nucleotide positions from population-genetics arguments, and accounting for among-lineage and among-gene effects, we propose a reconstruction of the variation in gBGC intensity at the scale of placental mammals, and of its scaling with body-size and karyotypic traits. Our results are compatible with a simple population genetics model relating gBGC to effective population size and recombination rate. In addition, among-gene variation and phylogenetic patterns of exon-specific levels of gBGC reveal the presence of rugged recombination landscapes, and suggest that short-lived recombination hot-spots are a general feature of placentals. Across placental mammals, variation in gBGC strength spans two orders of magnitude, at its lowest in apes, strongest in lagomorphs, microbats or tenrecs, and near or above the nearly neutral threshold in most other lineages. Combined with among-gene variation, such high levels of biased gene conversion are likely to significantly impact midly selected positions, and to represent a substantial mutation load. Altogether, our analysis suggests a more important role of gBGC in placental genome evolution, compared with what could have been anticipated from studies conducted in anthropoid primates.