Mutation accumulation in selfing populations under fluctuating selection

Selfing species are prone to extinction, possibly because highly selfing populations can suffer from a continuous accumulation of deleterious mutations, a process analogous to Muller's ratchet in asexual populations. However, current theory provides little insight into which types of genes are most likely to accumulate deleterious alleles and what environmental circumstances may accelerate genomic degradation. Here, we investigate temporal changes in the environment that cause fluctuations in the strength of purifying selection. We simulate selfing populations with genomes containing a mixture of loci experiencing constant selection and loci experiencing selection that fluctuates in strength (but not direction). Even when both types of loci experience the same average strength of selection, loci under fluctuating selection contribute disproportionately more to deleterious mutation accumulation. Moreover, the presence of loci experiencing fluctuating selection in the genome increases the deleterious fixation rate at loci under constant selection; under most realistic scenarios, this effect of linked selection can be attributed to a reduction in N_e. Fluctuating selection is particularly injurious when selective environments are strongly autocorrelated over time and when selection is concentrated into rare bouts of strong selection. These results imply that loci under fluctuating selection are likely important drivers of extinction in selfing species.     

RIBOSOMAL RNA GENE DIVERSITY,EFFECTIVE POPULATION SIZE,AND EVOLUTIONARY LONGEVITY IN ASEXUAL GLOMEROMYCOTA

Arbuscular mycorrhizal fungi (phylum Glomeromycota) are among the oldest and most successful symbionts of land plants. With no evidence of sexual reproduction, their evolutionary success is inconsistent with the prediction that asexual taxa are vulnerable to extinction due to accumulation of deleterious mutations. To explore why Glomeromycota defy this prediction, we studied ribosomal RNA (rRNA) gene evolution in the Claroideoglomus lineage and estimated effective population size, N_e, in C. etunicatum. We found that rRNA genes of these fungi exhibit unusual and complex patterns of molecular evolution. In C. etunicatum, these patterns can be collectively explained by an unexpectedly large N_e combined with imperfect genome‐wide and population‐level rRNA gene repeat homogenization. The mutations accumulated in rRNA gene sequences indicate that natural selection is effective at purging deleterious mutations in the Claroideoglomus lineage, which is also consistent with the large N_e of C. etunicatum. We propose that in the near absence of recombination, asexual reproduction involving massively multinucleate spores typical for Glomeromycota is responsible for the improved efficacy of selection relative to drift. We postulate that large effective population sizes contribute to the evolutionary longevity of Glomeromycota.     

Sexual selection,redundancy and survival of the most beautiful

A model is described of a highly redundant complex organism that has overlapping banks of genes such that each vital function is specified by several different genetic systems. This generates a synergistic profile linking probability of survival to the number of deleterious mutations in the genome. Computer models show that there is a dynamic interaction between the mean number of new deleterious mutations per generation (X), the mean number of deleterious mutations in the genome of the population (Y) and percentage zygote survival (Zs). IncreasedX leads to increasedY and a fall in Zs but it takes several generations before a new equilibrium is reached. If sexual attraction is influenced by the number of deleterious mutations in the genome of individuals thenY is reduced and Zs increased for any given value ofX. This fall inY and rise in Zs is more marked in polygamous than monogamous mating systems. The model is specified such that deleterious mutations can occur without any observable or measurable effect on function. Thus sexual selection, in this organism, for low levels of deleterious mutations cannot be based on assessment of performance. Instead it is based on a simple symmetrical surface pattern that is flawlessly reproduced by organisms with no deleterious mutations, but is less than perfect, and therefore less attractive, if genetic systems have been deleted. A complex vital task requires a system with a high level of redundancy that acts so that the loss of one component has no observable effect and therefore cannot be used for sexual selection. The reproduction of a beautiful surface pattern also requires a low error, high redundancy genetic system; however, in this case there is advantage if a single deleterious mutation produces a recognisable change. This leads to the conclusion that sexual selection and sexual attraction should be based on beauty rather than utility, and explains the common observation in nature that it is the most beautiful that survive.     

Fast Accumulation of Nonsynonymous Mutations on the Female-Specific W Chromosome in Birds

Following cessation of recombination during sex chromosome evolution, the nonrecombining sex chromosome is affected by a number of degenerative forces, possibly resulting in the fixation of deleterious mutations. This might take place because of weak selection against recessive or partly recessive deleterious mutations due to permanent heterozygosity of nonrecombining chromosomes. Furthermore, population genetic processes, such as selective sweeps, background selection, and Muller’s ratchet, result in a reduction in N_e, which increase the likelihood of fixation of deleterious mutations. Theory thus predicts that nonrecombining genes should show increased levels of nonsynonymous (d_N) to synonymous substitutions (d_S). We tested this in an avian system by estimating the ratio between d_N and d_S in six gametologous gene pairs located on the Z chromosome and the nonrecombining, female-specific W chromosome. In comparisons, we found a significantly higher d_N/d_S ratio for the W-linked than the Z-linked copy in three of the investigated genes. In a concatenated alignment of all six genes, the d_N/d_S ratio was six times higher for W-linked than Z-linked genes. By using human and mouse as outgroup in maximum likelihood analyses, W-linked genes were found to evolve differently compared with their Z-linked gametologues and outgroup sequences. This seems not to be a consequence of functional diversification because d_N/d_S ratios between gametologous gene copies were consistently low. We conclude that deleterious mutations are accumulating at a high rate on the avian W chromosome, probably as a result of the lack of recombination in this female-specific chromosome. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Deborah Charlesworth]     

FURTHER SIMULATION STUDIES ON EVOLUTION BY GENE DUPLICATION

In order to understand the origin of multigene families, Monte Carlo simulations were performed to see how a genetic system evolves under unequal crossing-over, mutation, random genetic drift and natural selection, starting from a single gene copy. Both haploid and diploid models were examined. Beneficial, neutral, and detrimental mutations were incorporated, and “positive” selection favors those chromosomes (haploid) or individuals (diploid) with more beneficial mutations than others. The same model for haploids was previously investigated with special reference to the evolution of gene organization, and the ratio of the numbers of beneficial genes to pseudogenes was found to be a rough indicator of the relative strengths of positive and negative (against deleterious alleles) natural selection (Ohta, 1987b). In the present paper, the evolution of gene organization and of sequence divergence among genes in the multigene family is examined. It is shown that positive selection accelerates the accumulation of arrays containing different beneficial mutations, but that total divergence including both neutral and beneficial mutations is not very sensitive to positive selection, under this model. The proportion of beneficial mutations in the total mutations accumulated is a better indicator of positive selection than is the total divergence. It is pointed out that various observed examples in which amino-acid substitutions are accelerated, as compared with synonymous substitutions in duplicated genes (Li, 1985), may reflect the effect of selection similar to the present scheme. The diploid model is shown to be more efficient for accumulating beneficial mutations in duplicated genes than the haploid one, and the relevance of this finding to the advantage of sexual reproduction is discussed.     

Fitness decline in spontaneous mutation accumulation lines of Caenorhabditis elegans with varying effective population sizes

The rate and fitness effects of new mutations have been investigated by mutation accumulation (MA) experiments in which organisms are maintained at a constant minimal population size to facilitate the accumulation of mutations with minimal efficacy of selection. We evolved 35 MA lines of Caenorhabditis elegans in parallel for 409 generations at three population sizes (N = 1, 10, and 100), representing the first spontaneous long-term MA experiment at varying population sizes with corresponding differences in the efficacy of selection. Productivity and survivorship in the N = 1 lines declined by 44% and 12%, respectively. The average effects of deleterious mutations in N = 1 lines are estimated to be 16.4% for productivity and 11.8% for survivorship. Larger populations (N = 10 and 100) did not suffer a significant decline in fitness traits despite a lengthy and sustained regime of consecutive bottlenecks exceeding 400 generations. Together, these results suggest that fitness decline in very small populations is dominated by mutations with large deleterious effects. It is possible that the MA lines at larger population sizes contain a load of cryptic deleterious mutations of small to moderate effects that would be revealed in more challenging environments.     

THE CONTRIBUTION OF NEW MUTATIONS TO GENOTYPE-ENVIRONMENT INTERACTION FOR FITNESS IN DROSOPHILA MELANOGASTER

Many studies have documented the existence of genotype-environment interaction (GEI) for traits closely related to fitness in natural populations. A type of GEI that is commonly observed is changes in the fitness ranking of genetic groups (families, clones, or inbred lines) in different environments. We refer to such changes in ranking as crossing of reaction norms for fitness. A common interpretation of crossing of reaction norms for fitness is that selection favors different alleles in the different environments (i.e., that “trade-offs” exist). If this is the case, selection could maintain genetic variation, and even lead to reproductive isolation between subpopulations using different environments. Even if the same alleles are favored in every environment, however, deleterious mutations that vary in the magnitude of their effect depending on environment could cause reaction norms for fitness to cross. If deleterious mutations with environment-dependent effects are responsible for maintaining much of the variation leading to crossing of reaction norms for fitness in natural populations, it should be possible to observe crossing of reaction norms for fitness among otherwise genetically identical lines bearing newly arisen spontaneous mutations. We examined the contribution of new mutations to GEI for fitness in Drosophila melanogaster. Eighteen lines were derived from a common, highly inbred base stock, and maintained at a population size of 10 pairs for over 200 generations, to allow them to accumulate spontaneous mutations. Because of the small population size of the lines, selection against mildly deleterious mutations should have been relatively ineffective. The lines were tested for productivity (number of surviving adult progeny from a standard number of parents) in five different environmental treatments, comprising different food media, temperatures, and levels of competition. The lines showed highly significant GEI for productivity, owing largely to considerable changes in ranking in the different environments. We conclude that mutations that are deleterious on average, but whose quantitative effects depend on environment, could be responsible for maintaining much of the variation leading to crossing of reaction norms for fitness that has been observed in samples of D. melanogaster from the wild.     

Revisiting the notion of deleterious sweeps

It has previously been shown that, conditional on its fixation, the time to fixation of a semi-dominant deleterious autosomal mutation in a randomly mating population is the same as that of an advantageous mutation. This result implies that deleterious mutations could generate selective sweep-like effects. Although their fixation probabilities greatly differ, the much larger input of deleterious relative to beneficial mutations suggests that this phenomenon could be important. We here examine how the fixation of mildly deleterious mutations affects levels and patterns of polymorphism at linked sites—both in the presence and absence of interference amongst deleterious mutations—and how this class of sites may contribute to divergence between-populations and species. We find that, while deleterious fixations are unlikely to represent a significant proportion of outliers in polymorphism-based genomic scans within populations, minor shifts in the frequencies of deleterious mutations can influence the proportions of private variants and the value of F_ST after a recent population split. As sites subject to deleterious mutations are necessarily found in functional genomic regions, interpretations in terms of recurrent positive selection may require reconsideration.     

RECURRENT AND RECENT SELECTIVE SWEEPS IN THE piRNA PATHWAY

Uncontrolled transposable element (TE) insertions and excisions can cause chromosome breaks and mutations with dramatic deleterious effects. The PIWI interacting RNA (piRNA) pathway functions as an adaptive TE silencing system during germline development. Several essential piRNA pathway proteins appear to be rapidly evolving, suggesting that TEs and the silencing machinery may be engaged in a classical “evolutionary arms race.” Using a variety of molecular evolutionary and population genetic approaches, we find that the piRNA pathway genes rhino, krimper, and aubergine show patterns suggestive of extensive recurrent positive selection across Drosophila species. We speculate that selection on these proteins reflects crucial roles in silencing unfamiliar elements during vertical and horizontal transmission of TEs into naïve populations and species, respectively.     

SEXUAL SELECTION COUNTERACTS EXTINCTION OF SMALL POPULATIONS OF THE BULB MITES

Genetic drift in small populations can increase frequency of deleterious recessives and consequently lead to inbreeding depression and population extinction. On the other hand, as homozygosity at deleterious recessives increases, they should be purged from populations more effectively by selection. Sexual selection has been postulated to strengthen selection against deleterious mutations, and should thus decrease extinction rate and intensify purging of inbreeding depression. We tested these predictions in the bulb mite Rhizoglyphus robini. We created 100 replicate lines of small populations (five males and five females) and in half of them experimentally removed sexual selection by enforcing monogamy. The lines were propagated for eight generations and then assayed for purging of inbreeding depression. We found that proportion of lines which went extinct was lower with sexual selection than without. We also found evidence for purging of inbreeding depression in the lines with sexual selection, but not in lines without sexual selection. Our results suggest that purging of inbreeding depression was more effective against mutations with relatively large deleterious effects. Thus, although our data clearly indicate a positive impact of sexual selection on short‐term survival of bottlenecked populations, long‐term consequences are less clear as they may be negatively impacted by accumulation of deleterious mutations of small effect.     

Insights into the evolutionary process from patterns of DNA sequence variability

Studies of DNA sequence diversity, particularly in Drosophila, reveal both the complexity of natural selection and the importance of the interaction of natural selection and variation in rates of recombination within genomes and between species in determining levels of sequence variability in different genes and different species. The impact of both adaptive and deleterious mutations are evident. Extension of these types of studies to other organisms has begun in earnest.     

Apparent directional selection by biased pleiotropic mutation

Pleiotropic effects of deleterious mutations are considered to be among the factors responsible for genetic constraints on evolution by long-term directional selection acting on a quantitative trait. If pleiotropic phenotypic effects are biased in a particular direction, mutations generate apparent directional selection, which refers to the covariance between fitness and the trait owing to a linear association between the number of mutations possessed by individuals and the genotypic values of the trait. The present analysis has shown how the equilibrium mean value of the trait is determined by a balance between directional selection and biased pleiotropic mutations. Assuming that genes act additively both on the trait and on fitness, the total variance-standardized directional selection gradient was decomposed into apparent and true components. Experimental data on mutation bias from the bristle traits of Drosophila and life history traits of Daphnia suggest that apparent selection explains a small but significant fraction of directional selection pressure that is observed in nature; the data suggest that changes induced in a trait by biased pleiotropic mutation (i.e., by apparent directional selection) are easily compensated for by (true) directional selection.     

A POPULATION GENETIC THEORY OF CANALIZATION

Canalization is the suppression of phenotypic variation. Depending on the causes of phenotypic variation, one speaks either of genetic or environmental canalization. Genetic canalization describes insensitivity of a character to mutations, and the insensitivity to environmental factors is called environmental canalization. Genetic canalization is of interest because it influences the availability of heritable phenotypic variation to natural selection, and is thus potentially important in determining the pattern of phenotypic evolution. In this paper a number of population genetic models are considered of a quantitative character under stabilizing selection. The main purpose of this study is to define the population genetic conditions and constraints for the evolution of canalization. Environmental canalization is modeled as genotype specific environmental variance. It is shown that stabilizing selection favors genes that decrease environmental variance of quantitative characters. However, the theoretical limit of zero environmental variance has never been observed. Of the many ways to explain this fact, two are addressed by our model. It is shown that a “canalization limit” is reached if canalizing effects of mutations are correlated with direct effects on the same character. This canalization limit is predicted to be independent of the strength of stabilizing selection, which is inconsistent with recent experimental data (Sterns et al. 1995). The second model assumes that the canalizing genes have deleterious pleiotropic effects. If these deleterious effects are of the same magnitude as all the other mutations affecting fitness very strong stabilizing selection is required to allow the evolution of environmental canalization. Genetic canalization is modeled as an influence on the average effect of mutations at a locus of other genes. It is found that the selection for genetic canalization critically depends on the amount of genetic variation present in the population. The more genetic variation, the stronger the selection for canalizing effects. All factors that increase genetic variation favor the evolution of genetic canalization (large population size, high mutation rate, large number of genes). If genetic variation is maintained by mutation-selection balance, strong stabilizing selection can inhibit the evolution of genetic canalization. Strong stabilizing selection eliminates genetic variation to a level where selection for canalization does not work anymore. It is predicted that the most important characters (in terms of fitness) are not necessarily the most canalized ones, if they are under very strong stabilizing selection (k > 0.2V_e). The rate of decrease of mutational variance V_m is found to be less than 10% of the initial V_m. From this result it is concluded that characters with typical mutational variances of about 10^–3 V_e are in a metastable state where further evolution of genetic canalization is too slow to be of importance at a microevolutionary time scale. The implications for the explanation of macroevolutionary patterns are discussed.     

Mating behavior as an indicator of quality of Drosophila subobscura males?

According to current theoretical predictions, any deleterious mutations that reduce nonsexual fitness may have a negative influence on mating success. This means that sexual selection may remove deleterious mutations from the populations. Males of good genetic quality should be more successful in mating, compared to the males of lower genetic quality. As mating success is a condition dependent trait, large fractions of the genome may be a target of sexual selection and many behavioral traits are likely to be condition dependent. We manipulated the genetic quality of Drosophila subobscura males by inducing mutations with ionizing radiation and observed the effects of the obtained heterozygous mutations on male mating behavior: courtship occurrence, courtship latency, mating occurrence, latency to mating and duration of mating. We found possible effects of mutations. Females mated more frequently with male progeny of nonirradiated males and that these males courted females faster compared to the male progeny of irradiated males. Our findings indicate a possible important role of sexual selection in purging deleterious mutations.     

Strong sexual selection in males against a mutation load that reduces offspring production in seed beetles

Theory predicts that sexual reproduction can increase population viability relative to asexual reproduction by allowing sexual selection in males to remove deleterious mutations from the population without large demographic costs. This requires that selection acts more strongly in males than females and that mutations affecting male reproductive success have pleiotropic effects on population productivity, but empirical support for these assumptions is mixed. We used the seed beetle Callosobruchus maculatus to implement a three‐generation breeding design where we induced mutations via ionizing radiation (IR) in the F₀ generation and measured mutational effects (relative to nonirradiated controls) on an estimate of population productivity in the F₁ and effects on sex‐specific competitive lifetime reproductive success (LRS) in the F₂. Regardless of whether mutations were induced via F₀ males or females, they had strong negative effects on male LRS, but a nonsignificant influence on female LRS, suggesting that selection is more efficient in removing deleterious alleles in males. Moreover, mutations had seemingly shared effects on population productivity and competitive LRS in both sexes. Thus, our results lend support to the hypothesis that strong sexual selection on males can act to remove the mutation load on population viability, thereby offering a benefit to sexual reproduction.     

Genetic Constraints on Protein Evolution

ABSTRACT

Evolution requires the generation and optimization of new traits (“adaptation”) and involves the selection of mutations that improve cellular function. These mutations were assumed to arise by selection of neutral mutations present at all times in the population. Here we review recent evidence that indicates that deleterious mutations are more frequent in the population than previously recognized and that these mutations play a significant role in protein evolution through continuous positive selection. Positively selected mutations include adaptive mutations, i.e. mutations that directly affect enzymatic function, and compensatory mutations, which suppress the pleiotropic effects of adaptive mutations. Compensatory mutations are by far the most frequent of the two and would allow potentially adaptive but deleterious mutations to persist long enough in the population to be positively selected during episodes of adaptation. Compensatory mutations are, by definition, context-dependent and thus constrain the paths available for evolution. This provides a mechanistic basis for the examples of highly constrained evolutionary landscapes and parallel evolution reported in natural and experimental populations. The present review article describes these recent advances in the field of protein evolution and discusses their implications for understanding the genetic basis of disease and for protein engineering in vitro.     

Sequence Variation in the tRNA Genes of Human Mitochondrial DNA

Recent analyses have shown that nonsynonymous variation in human mitochondrial DNA (mtDNA) contains nonneutral variants, suggesting the presence of mildly deleterious mutations. Many of the disease-causing mutations in mtDNA occur in the genes encoding the tRNAs. Nucleotide sequence variation in these genes has not been studied in human populations, nor have the structural consequences of nucleotide substitutions in tRNA molecules been examined. We therefore determined the nucleotide sequences of the 22 tRNA genes in the mtDNA of 477 Finns and, also, obtained 435 European sequences from the MitoKor database. No differences in population polymorphism indices were found between the two data sets. We assessed selective constraints against various tRNA domains by comparing allele frequencies between these domains and the synonymous and nonsynonymous sites, respectively. All tRNA domains except the variable loop were more conserved than synonymous sites, and T stem and D stem were more conserved than the respective loops. We also analyzed the energetic consequences of the 96 polymorphisms recovered in the two data sets or in the Mitomap database. The minimum free energy (ΔG) was calculated using the free energy rules as implemented in mfold version 3.1. The ΔG’s were normally distributed among the 22 wild-type tRNA genes, whereas the 96 polymorphic tRNAs departed significantly from a normal distribution. The largest differences in ΔG between the wild-type and the polymorphic tRNAs in the Finnish population tended to be in the polymorphisms that were present at low frequencies. Allele frequency distributions and minimum free energy calculations both suggested that some polymorphisms in tRNA genes are nonneutral.Reviewing Editor: Dr. Rüdiger Cerff     

Significance of Population Size on the Fixation of Nonsynonymous Mutations in Genes Under Varying Levels of Selection Pressure

Previous studies observed a higher ratio of divergences at nonsynonymous and synonymous sites (ω = d_N/d_S) in species with a small population size compared to that estimated for those with a large population size. Here we examined the theoretical relationship between ω, effective population size (N_e), and selection coefficient (s). Our analysis revealed that when purifying selection is high, ω of species with small N_e is much higher than that of species with large N_e. However the difference between the two ω reduces with the decline in selection pressure (s → 0). We examined this relationship using primate and rodent genes and found that the ω estimated for highly constrained genes of primates was up to 2.9 times higher than that obtained for their orthologous rodent genes. Conversely, for genes under weak purifying selection the ω of primates was only 17% higher than that of rodents. When tissue specificity was used as a proxy for selection pressure we found that the ω of broadly expressed genes of primates was up to 2.1-fold higher than that of their rodent counterparts and this difference was only 27% for tissue specific genes. Since most of the nonsynonymous mutations in constrained or broadly expressed genes are deleterious, fixation of these mutations is influenced by N_e. This results in a higher ω of these genes in primates compared to those from rodents. Conversely, the majority of nonsynonymous mutations in less-constrained or tissue-specific genes are neutral or nearly neutral and therefore fixation of them is largely independent of N_e, which leads to the similarity of ω in primates and rodents.     

Deleterious Mutations Accumulate Faster in Allopolyploid Than Diploid Cotton (Gossypium) and Unequally between Subgenomes

Whole-genome duplication (polyploidization) is among the most dramatic mutational processes in nature, so understanding how natural selection differs in polyploids relative to diploids is an important goal. Population genetics theory predicts that recessive deleterious mutations accumulate faster in allopolyploids than diploids due to the masking effect of redundant gene copies, but this prediction is hitherto unconfirmed. Here, we use the cotton genus (Gossypium), which contains seven allopolyploids derived from a single polyploidization event 1–2 Million years ago, to investigate deleterious mutation accumulation. We use two methods of identifying deleterious mutations at the nucleotide and amino acid level, along with whole-genome resequencing of 43 individuals spanning six allopolyploid species and their two diploid progenitors, to demonstrate that deleterious mutations accumulate faster in allopolyploids than in their diploid progenitors. We find that, unlike what would be expected under models of demographic changes alone, strongly deleterious mutations show the biggest difference between ploidy levels, and this effect diminishes for moderately and mildly deleterious mutations. We further show that the proportion of nonsynonymous mutations that are deleterious differs between the two coresident subgenomes in the allopolyploids, suggesting that homoeologous masking acts unequally between subgenomes. Our results provide a genome-wide perspective on classic notions of the significance of gene duplication that likely are broadly applicable to allopolyploids, with implications for our understanding of the evolutionary fate of deleterious mutations. Finally, we note that some measures of selection (e.g., dN/dS, π_N/π_S) may be biased when species of different ploidy levels are compared.     

Detection of selection signatures in Limousin cattle using whole-genome resequencing

Limousin, a renowned beef breed originating from central France, has been selectively bred over the last 100 years to improve economically important traits. We used whole-genome sequencing data from 10 unrelated Limousin bull calves to detect polymorphisms and identify regions under selection. A total of 13 943 766 variants were identified. Moreover, 311 852 bi-allelic SNPs and 92 229 indels located on autosomes were fixed for the alternative allele in all sequenced animals, including the previously reported missense deleterious F94L mutation in MSTN. We performed a whole-genome screen to discover genomic regions with excess homozygosity, using the pooled heterozygosity score and identified 171 different candidate selective sweeps. In total, 68 candidate genes were found in only 57 of these regions, indicating that a large fraction of the genome under selection might lie in non-coding regions and suggesting that a majority of adaptive mutations might be regulatory in nature. Many QTL were found within candidate selective sweep regions, including QTL associated with shear force or carcass weight. Among the putative selective sweeps, we located genes (MSTN, NCKAP5, RUNX2) that potentially contribute to important phenotypes in Limousin. Several candidate regions and genes under selection were also found in previous genome-wide selection scans performed in Limousin. In addition, we were able to pinpoint candidate causative regulatory polymorphisms in GRIK3 and RUNX2 that might have been under selection. Our results will contribute to improved understanding of the mechanisms and targets of artificial selection and will facilitate the interpretation of GWASs performed in Limousin.