Identifying signatures of sexual selection using genomewide selection components analysis

Sexual selection must affect the genome for it to have an evolutionary impact, yet signatures of selection remain elusive. Here we use an individual‐based model to investigate the utility of genome‐wide selection components analysis, which compares allele frequencies of individuals at different life history stages within a single population to detect selection without requiring a priori knowledge of traits under selection. We modeled a diploid, sexually reproducing population and introduced strong mate choice on a quantitative trait to simulate sexual selection. Genome‐wide allele frequencies in adults and offspring were compared using weighted F_ST values. The average number of outlier peaks (i.e., those with significantly large F_ST values) with a quantitative trait locus in close proximity (“real” peaks) represented correct diagnoses of loci under selection, whereas peaks above the F_ST significance threshold without a quantitative trait locus reflected spurious peaks. We found that, even with moderate sample sizes, signatures of strong sexual selection were detectable, but larger sample sizes improved detection rates. The model was better able to detect selection with more neutral markers, and when quantitative trait loci and neutral markers were distributed across multiple chromosomes. Although environmental variation decreased detection rates, the identification of real peaks nevertheless remained feasible. We also found that detection rates can be improved by sampling multiple populations experiencing similar selection regimes. In short, genome‐wide selection components analysis is a challenging but feasible approach for the identification of regions of the genome under selection.     

Genomic signatures reveal selection of characteristics within and between Meishan pig populations

The Chinese Meishan pig breed is well known for its high prolificacy. Moreover, this breed can be divided into three types based on their body size: big Meishan, middle Meishan (MMS) and small Meishan (SMS) pigs. Few studies have reported on the genetic signatures of Meishan pigs, particularly on a genome‐wide scale. Exploring for genetic signatures could be quite valuable for revealing the genetic architecture of phenotypic variation. Thus, we performed research in two parts based on the genome reducing and sequencing data of 143 Meishan pigs (74 MMS pigs, 69 SMS pigs). First, we detected the selection signatures among all Meishan pigs studied using the relative extended haplotype homozygosity test. Second, we detected the selection signatures between MMS and SMS pigs using the cross‐population extended haplotype homozygosity and F_ST methods. A total of 111 398 SNPs were identified from the sequenced genomes. In the population analysis, the most significant genes were associated with the mental development (RGMA), reproduction (HDAC4, FOXL2) and lipid metabolism (ACACB). From the cross‐population analysis, we detected genes related to body weight (SPDEF, PACSIN1) in both methods. We suggest that rs341373351, located within the PACSIN1 gene, might be the causal variant. This study may have achieved consistency between selection signatures and characteristics within and between Meishan pig populations. These findings can provide insight into investigating the molecular background of high prolificacy and body size in pig.     

Disentangling the effects of recombination,selection, and demography on the genetic variation at a major histocompatibility complex class II gene in the alpine chamois

The major histocompatibility complex (MHC) harbours some of the most polymorphic loci in vertebrate genomes. MHC genes are thought to be subject to some form of balancing selection, most likely pathogen‐mediated selection. Hence, MHC genes are excellent candidates for exploring adaptive processes. In this study, we investigated the genetic variation at exon 2 of the DRB class II MHC locus in 191 alpine chamois (Rupicapra rupicapra) from 10 populations in the eastern Alps of Italy. In particular, we were interested in distinguishing and estimating the relative impact of selective and demographic factors, while taking into account the confounding effect of recombination. The extremely high d_n/d_s ratio and the presence of trans‐species polymorphisms suggest that a strong long‐term balancing selection effect has been operating at this locus throughout the evolutionary history of this species. We analysed patterns of genetic variation within and between populations, and the mitochondrial D‐loop polymorphism patterns were analysed to provide a baseline indicator of the effects of demographic processes. These analyses showed that (i) the chamois experienced a demographic decline in the last 5000–30 000 years, most likely related to the postglacial elevation in temperature; (ii) this demographic process can explain the results of neutrality tests applied to MHC variation within populations, but cannot justify the much weaker divergence between populations implied by MHC as opposed to mitochondrial DNA; (iii) similar sets of divergent alleles are probably maintained with similar frequencies by balancing selection in different populations, and this mechanism is also operating in small isolated populations, which are strongly affected by drift.     

The origin of selection signatures on bovine chromosome 6

The extent and pattern of linkage disequilibrium (LD) between closely spaced markers contain information about population history, including past population size and selection history. Selection signatures can be identified by comparing the LD surrounding a putative selected allele at a locus to the putative non-selected allele. In livestock populations, locations of selection signatures identified in this way should be correlated with QTL affecting production traits, as the populations have been under strong artificial selection for these traits. We used a dense SNP map of bovine chromosome 6 to characterize the pattern of LD on this chromosome in Norwegian Red cattle, a breed which has been strongly selected for milk production. The pattern of LD was generally consistent with strong selection in regions containing QTL affecting milk production traits, including a strong selection signature in a region containing a mutation known to affect milk production. The results demonstrate that in livestock populations, the origin of selection signatures will often be QTL for livestock production traits, and illustrate the value of selection signatures in uncovering new mutations with potential effects on quantitative traits.     

Genetic structure and contrasting selection pattern at two major histocompatibility complex genes in wild house mouse populations

The mammalian major histocompatibility complex (MHC) is a tightly linked cluster of immune genes, and is often thought of as inherited as a unit. This has led to the hope that studying a single MHC gene will reveal patterns of evolution representative of the MHC as a whole. In this study we analyse a 1000-km transect of MHC variation traversing the European house mouse hybrid zone to compare signals of selection and patterns of diversification at two closely linked MHC class II genes, H-2Aa and H-2Eb. We show that although they are 0.01 cM apart (that is, recombination is expected only once in 10 000 meioses), disparate evolutionary patterns were detected. H-2Aa shows higher allelic polymorphism, faster allelic turnover due to higher mutation rates, stronger positive selection at antigen-binding sites and higher population structuring than H-2Eb. H-2Eb alleles are maintained in the gene pool for longer, including over separation of the subspecies, some H-2Eb alleles are positively and others negatively selected and some of the alleles are not expressed. We conclude that studies on MHC genes in wild-living vertebrates can give substantially different results depending on the MHC gene examined and that the level of polymorphism in a related species is a poor criterion for gene choice.     

GAL4 enhancer trap targeting of the Drosophila sex determination gene fruitless

The fru4 allele of the sex determination gene fruitless is induced by insertion of a P[lacZ,ry+] enhancer trap element. This insert also acts to disrupt expression of the fru P1 promoter derived male-specific proteins, consequently impairing male courtship behavior. fru4 maps less than 2 kb upstream of the fru P3 promoter, whose function is essential for viability. We replaced this insert with a GAL4 element, P[GAL4,w+], recovering two lines with insertions in opposite orientations at the locus, one of which demonstrated fru-specific mutant phenotypes. Reporter expression of these lines recapitulated that of P3- and P4-derived proteins which, when correlated with a developmental and tissue specific survey of fru promoters' activities, uncovered a previously unsuspected complexity of fru regulation. These novel fru alleles provide the tools for manipulation of fru-expressing cells, allowing the consequent effects to be related back to specific fru functions and the regulatory units controlling these activities.     

The effect of polymorphisms in the enhancer of split gene complex on bristle number variation in a large wild-caught cohort of Drosophila melanogaster

The Enhancer of split complex [E(spl)-C] in Drosophila encompasses a variety of functional elements controlling bristle patterning and on the basis of prior work is a strong candidate for harboring alleles having subtle effects on bristle number variation. Here we extend earlier studies identifying associations between complex phenotypes and polymorphisms segregating among inbred laboratory lines of Drosophila and test the influence of E(spl)-C on bristle number variation in a natural cohort. We describe results from an association mapping study using 203 polymorphisms spread throughout the E(spl)-C genotyped in 2000 wild-caught Drosophila melanogaster. Despite power to detect associations accounting for as little as 2% of segregating variation for bristle number, and saturating the region with single-nucleotide polymorphisms (SNPs), we identified no single SNP marker showing a significant (additive over loci) effect after correcting for multiple tests. Using a newly developed test we conservatively identify six regions of the E(spl)-C in which the insertion of transposable elements as a class contributes to variation in bristle number, apparently in a sex- or trait-limited fashion. Finally, we carry out all possible 20,503 two-way tests for epistasis and identify a slight excess of marginally significant interactions, although none survive multiple-testing correction. It may not be straightforward to extend the results of laboratory-based association studies to natural populations.     

Whole‐genome SNP data unravel population structure and signatures of selection for black plumage of indigenous chicken breeds from Jiangxi province,China

Ten indigenous chicken breeds were originally distributed in Jiangxi Province, China, and they define a critical component of Chinese chicken genetic resources. We have investigated the population genetics of seven Jiangxi chicken breeds using 600K chicken BeadChip SNP data. To provide a genome‐wide perspective for the population structure of all 10 Jiangxi chicken breeds, we herein genotyped 78 additional individuals from the seven breeds and 63 chickens from three uninvestigated breeds—Yugan Black (YG), Nancheng Black (NC) and Wanzai Yellow using 55K chicken SNP arrays. We then explored merged data of 17 101 SNPs from 235 individuals to infer the population structure of the 10 breeds. We showed that NC and YG are two regional populations of the same breed, as individuals from the two populations clustered together to form a branch separate from the other breeds in the neighbor‐joining tree, they always grouped together in multidimensional principal component analyses and they displayed an identical pattern of ancestral lineage composition. Hence, NC and YG should be considered a single breed in the state‐supported conservation scheme. Moreover, we conducted a genome scan for signatures of selection for black plumage. bayescan and hapflk analyses of two contrasting groups (three black‐feathered breeds vs. six non‐black‐feathered breeds) consistently detected 25 putative regions under selection. Nine pigmentation‐ associated genes (DCT, SLC24A5, SLC30A4, MYO5A, CYP19A1, NADK2, SLC45A2, GNAQ and DCP2) reside within these regions, and these genes are interesting candidates for black plumage and provide a starting point for further identification of causative mutations for black feathers in chicken.     

Nucleotide variability at the acetyl coenzyme A carboxylase gene and the signature of herbicide selection in the grass weed Alopecurus myosuroides (Huds.)

Acetyl coenzyme A carboxylase (ACCase) is the target of highly effective herbicides. We investigated the nucleotide variability of the ACCase gene in a sample of 18 black-grass (Alopecurus myosuroides [Huds.]) populations to search for the signature of herbicide selection. Sequencing 3,396 bp encompassing ACCase herbicide-binding domain in 86 individuals revealed 92 polymorphisms, which formed 72 haplotypes. The ratio of nonsynonymous versus synonymous substitutions was very low, in agreement with ACCase being a vital metabolic enzyme. Within black grass, most nonsynonymous substitutions were related to resistance to ACCase-inhibiting herbicides. Differentiation between populations was strong, in contrast to expectations for an allogamous, annual plant. Significant H tests revealed recent hitchhiking events within populations. These results were consistent with recent and local positive selection. We propose that, although they have only been used since at most 15 black-grass generations, ACCase-inhibiting herbicides have exerted a positive selection targeting resistant haplotypes that has been strong enough to have a marked effect upon ACCase nucleotide diversity. A minimum-spanning network of nonrecombinant haplotypes revealed multiple, independent apparitions of resistance-associated mutations. This study provides the first evidence for the signature of ongoing, recent, pesticide selection upon variation at the gene encoding the targeted enzyme in natural plant populations.     

Landscape genomic approach to detect selection signatures in locally adapted Brazilian swine genetic groups

Samples of 191 animals from 18 different Brazilian locally adapted swine genetic groups were genotyped using Illumina Porcine SNP60 BeadChip in order to identify selection signatures related to the monthly variation of Brazilian environmental variables. Using BayeScan software, 71 SNP markers were identified as F_ST outliers and 60 genotypes (58 markers) were found by Samβada software in 371 logistic models correlated with 112 environmental variables. Five markers were identified in both methods, with a Kappa value of 0.073 (95% CI: 0.011–0.134). The frequency of these markers indicated a clear north–south country division that reflects Brazilian environmental differences in temperature, solar radiation, and precipitation. Global spatial territory correlation for environmental variables corroborates this finding (average Moran's I = 0.89, range from 0.55 to 0.97). The distribution of alleles over the territory was not strongly correlated with the breed/genetic groups. These results are congruent with previous mtDNA studies and should be used to direct germplasm collection for the National gene bank.     

Genome-wide signatures of population bottlenecks and diversifying selection in European wolves

Genomic resources developed for domesticated species provide powerful tools for studying the evolutionary history of their wild relatives. Here we use 61K single-nucleotide polymorphisms (SNPs) evenly spaced throughout the canine nuclear genome to analyse evolutionary relationships among the three largest European populations of grey wolves in comparison with other populations worldwide, and investigate genome-wide effects of demographic bottlenecks and signatures of selection. European wolves have a discontinuous range, with large and connected populations in Eastern Europe and relatively smaller, isolated populations in Italy and the Iberian Peninsula. Our results suggest a continuous decline in wolf numbers in Europe since the Late Pleistocene, and long-term isolation and bottlenecks in the Italian and Iberian populations following their divergence from the Eastern European population. The Italian and Iberian populations have low genetic variability and high linkage disequilibrium, but relatively few autozygous segments across the genome. This last characteristic clearly distinguishes them from populations that underwent recent drastic demographic declines or founder events, and implies long-term bottlenecks in these two populations. Although genetic drift due to spatial isolation and bottlenecks seems to be a major evolutionary force diversifying the European populations, we detected 35 loci that are putatively under diversifying selection. Two of these loci flank the canine platelet-derived growth factor gene, which affects bone growth and may influence differences in body size between wolf populations. This study demonstrates the power of population genomics for identifying genetic signals of demographic bottlenecks and detecting signatures of directional selection in bottlenecked populations, despite their low background variability.     

Genetic diversity and genomic signatures of selection among cattle breeds from Siberia,eastern and northern Europe

Domestication in the near eastern region had a major impact on the gene pool of humpless taurine cattle (Bos taurus). As a result of subsequent natural and artificial selection, hundreds of different breeds have evolved, displaying a broad range of phenotypic traits. Here, 10 Eurasian B. taurus breeds from different biogeographic and production conditions, which exhibit different demographic histories and have been under artificial selection at various intensities, were investigated using the Illumina BovineSNP50 panel to understand their genetic diversity and population structure. In addition, we scanned genomes from eight breeds for signatures of diversifying selection. Our population structure analysis indicated six distinct breed groups, the most divergent being the Yakutian cattle from Siberia. Selection signals were shared (experimental P‐value < 0.01) with more than four breeds on chromosomes 6, 7, 13, 16 and 22. The strongest selection signals in the Yakutian cattle were found on chromosomes 7 and 21, where a miRNA gene and genes related to immune system processes are respectively located. In general, genomic regions indicating selection overlapped with known QTL associated with milk production (e.g. on chromosome 19), reproduction (e.g. on chromosome 24) and meat quality (e.g. on chromosome 7). The selection map created in this study shows that native cattle breeds and their genetic resources represent unique material for future breeding.