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.     

Twenty years of artificial directional selection have shaped the genome of the Italian Large White pig breed

In this study, we investigated at the genome‐wide level if 20 years of artificial directional selection based on boar genetic evaluation obtained with a classical BLUP animal model shaped the genome of the Italian Large White pig breed. The most influential boars of this breed (n = 192), born from 1992 (the beginning of the selection program of this breed) to 2012, with an estimated breeding value reliability of >0.85, were genotyped with the Illumina Porcine SNP60 BeadChip. After grouping the boars in eight classes according to their year of birth, filtered single nucleotide polymorphisms (SNPs) were used to evaluate the effects of time on genotype frequency changes using multinomial logistic regression models. Of these markers, 493 had a P_Bonferroni < 0.10. However, there was an increasing number of SNPs with a decreasing level of allele frequency changes over time, representing a continuous profile across the genome. The largest proportion of the 493 SNPs was on porcine chromosome (SSC) 7, SSC2, SSC8 and SSC18 for a total of 204 haploblocks. Functional annotations of genomic regions, including the 493 shifted SNPs, reported a few Gene Ontology terms that might underly the biological processes that contributed to increase performances of the pigs over the 20 years of the selection program. The obtained results indicated that the genome of the Italian Large White pigs was shaped by a directional selection program derived by the application of methodologies assuming the infinitesimal model that captured a continuous trend of allele frequency changes in the boar population.     

A genome‐wide scan for signatures of recent selection in Holstein cattle

The data from the newly available 50 K SNP chip was used for tagging the genome‐wide footprints of positive selection in Holstein–Friesian cattle. For this purpose, we employed the recently described Extended Haplotype Homozygosity test, which detects selection by measuring the characteristics of haplotypes within a single population. To assess formally the significance of these results, we compared the combination of frequency and the Relative Extended Haplotype Homozygosity value of each core haplotype with equally frequent haplotypes across the genome. A subset of the putative regions showing the highest significance in the genome‐wide EHH tests was mapped. We annotated genes to identify possible influence they have in beneficial traits by using the Gene Ontology database. A panel of genes, including FABP3, CLPN3, SPERT, HTR2A5, ABCE1, BMP4 and PTGER2, was detected, which overlapped with the most extreme P‐values. This panel comprises some interesting candidate genes and QTL, representing a broad range of economically important traits such as milk yield and composition, as well as reproductive and behavioural traits. We also report high values of linkage disequilibrium and a slower decay of haplotype homozygosity for some candidate regions harbouring major genes related to dairy quality. The results of this study provide a genome‐wide map of selection footprints in the Holstein genome, and can be used to better understand the mechanisms of selection in dairy cattle breeding.     

Detection of signatures of selective sweeps in the Blonde d'Aquitaine cattle breed

Identifying recent positive selection signatures in domesticated animals could provide information on genome response to strong directional selection from domestication and artificial selection and therefore could help in identifying mutations responsible for improved traits. We used genotyping data generated using Illumina's BovineSNP50 Genotyping BeadChips to identify selection signatures in the Blonde d'Aquitaine breed, a well‐muscled French beef breed. For this purpose, we employed a hidden Markov model‐based test, which detects selection by studying local variations in the allele frequency spectrum along the genome, within a single population. Three regions containing selective sweeps were identified. Annotation of genes located within these regions revealed interesting candidate genes. For example, myostatin (also known as GDF8), a known muscle growth factor inhibitor, is located within the selection signature region found on chromosome 2. In addition, we have identified chromosomal regions that show some evidence of selection within QTL regions for economically important traits. The results of this study could help to better understand the mechanisms related to the selection of the Blonde d'Aquitaine breed.     

Long‐term selection for litter size in swine results in shifts in allelic frequency in regions involved in reproductive processes

High‐density genotype data were analyzed in three lines of swine that express substantial variation in sow fertility to uncover regions of the genome potentially influenced during selection for litter size traits. The experimental lines examined include the Nebraska Index Line (NIL), which has been subjected to long‐term selection for litter size; a control line derived from the same population that founded NIL; and a commercial Duroc × Hampshire (D × H) population, in which no selection for litter size was practiced. Regions of the genome potentially affected by selection for litter size traits in NIL were determined by multiple lines of evidence, including altered allelic frequency compared to the other lines, loss of heterozygosity and relative extended haplotype homozygosity. Additionally, a genome‐wide association study for litter size traits was conducted in a population based on NIL and commercial maternal line genetics. Several genomic regions identified as putative signatures of selection overlapped with QTL for litter size traits. One of these regions, located on SSC2 (13–14 Mb), includes the candidate gene P2X3R, which plays a role in implantation and sustained release of hormones associated with reproductive processes. Sequencing identified synonymous SNPs in P2X3R that are fixed in NIL but polymorphic with nearly equal frequencies in the D × H line, indicating a potential role of P2X3R in sow fertility. These results suggest that data derived from these lines can help to uncover and understand a portion of the genetic variance associated with fertility traits in swine.     

基因组选择技术在农业动物育种中的应用

基因组选择(genomic selection, GS)是畜禽经济性状遗传改良的重要方法。随着高密度SNP芯片和二代测序价格的下降,GS技术越来越多被应用于奶牛、猪、鸡等农业动物育种中。然而,降低全基因组SNP分型成本、提高基因组育种值(genomic estimated breeding value,GEBV)估计准确性仍然是GS研究的主要难题。本文从全基因组SNP分型策略和GEBV估计模型两个方面进行了综述,并对目前GS技术在主要畜禽品种中的应用现状进行了介绍,以期为GS在农业动物育种中的深入开展提供借鉴和参考。     

Footprints of selection in the ancestral admixture of a New World Creole cattle breed

Admixed populations represent attractive biological models to study adaptive selection. Originating from several waves of recent introduction from European (EUT), African (AFT) and zebus (ZEB) cattle, New World Creole cattle allow investigating the response to tropical environmental challenges of these three ancestries. We here provide a detailed assessment of their genetic contributions to the Creole breed from Guadeloupe (CGU). We subsequently look for footprints of selection by combining results from tests based on the extent of haplotype homozygosity and the identification of excess/deficiency of local ancestry. To tackle these issues, 140 CGU individuals and 25 Brahman zebus from Martinique were genotyped at 44 057 SNPs. These data were combined to those available on 23 populations representative of EUT, AFT or ZEB. We found average proportions of 26.1%, 36.0% and 37.9% of EUT, AFT and ZEB ancestries in the CGU genome indicating a higher level of African and zebu ancestries than suggested by historical records. We further identified 23 genomic regions displaying strong signal of selection, most of them being characterized by an excess of ZEB local ancestry. Among the candidate gene underlying these regions, several are associated with reproductive functions (RXFP2, PMEPA1, IGFBP3, KDR, PPP1R8, TBXA2R and SLC7A5) and metabolism (PDE1B and CYP46A1). Finally, two genes (CENTD3 and SAMD12) are involved in cellular signalization of immune response. This study illustrates the relevance of admixed populations to identify footprints of selection by combining several tests straightforward to implement on large data sets.     

Footprints of recent selection and variability in breed composition in the Göttingen Minipig genome

The Göttingen Minipig (GMP) developed at the University of Göttingen is a synthetic breed that is widely used in medical research and toxicology. It combines the high fertility of the Vietnamese potbellied pig, the low body weight of the Minnesota Minipig and the white coat colour of the German Landrace pig. The aim of this study was to find genomic regions that may have undergone selection since the creation of the breed in the 1960s. Therefore, the whole genome was screened for footprints of recent selection based on single nucleotide polymorphism (SNP) genotypes from the Illumina Porcine SNP60 BeadChip using two methods: the extended haplotype homozygosity (EHH) test and the estimation of the genomic proportion of the three original breeds at each SNP using a Bayesian approach. Local deviations from the average genome‐wide breed composition were tested with a permutation‐based empirical test. Results for a comprehensive whole‐genome scan for both methods are presented. Several regions showing the highest P‐values in the EHH test are related to breeding goals relevant in the GMP, such as growth (SOCS2, TXN, DDR2 and GRB10 genes) and white colour (PRLR gene). Additionally, the calculated proportion of the founder breeds diverged significantly in many regions from the pedigree‐based expectations and the genome average. The results provide a genome‐wide map of selection signatures in the GMP, which leads to a better understanding of selection that took place over the last decades in GMP breed development.     

Population genomic footprints of selection and associations with climate in natural populations of Arabidopsis halleri from the Alps

Natural genetic variation is essential for the adaptation of organisms to their local environment and to changing environmental conditions. Here, we examine genomewide patterns of nucleotide variation in natural populations of the outcrossing herb Arabidopsis halleri and associations with climatic variation among populations in the Alps. Using a pooled population sequencing (Pool‐Seq) approach, we discovered more than two million SNPs in five natural populations and identified highly differentiated genomic regions and SNPs using F_ST‐based analyses. We tested only the most strongly differentiated SNPs for associations with a nonredundant set of environmental factors using partial Mantel tests to identify topo‐climatic factors that may underlie the observed footprints of selection. Possible functions of genes showing signatures of selection were identified by Gene Ontology analysis. We found 175 genes to be highly associated with one or more of the five tested topo‐climatic factors. Of these, 23.4% had unknown functions. Genetic variation in four candidate genes was strongly associated with site water balance and solar radiation, and functional annotations were congruent with these environmental factors. Our results provide a genomewide perspective on the distribution of adaptive genetic variation in natural plant populations from a highly diverse and heterogeneous alpine environment.     

Identifying footprints of directional and balancing selection in marine and freshwater three-spined stickleback (Gasterosteus aculeatus) populations

Natural selection is expected to leave an imprint on the neutral polymorphisms at the adjacent genomic regions of a selected gene. While directional selection tends to reduce within-population genetic diversity and increase among-population differentiation, the reverse is expected under balancing selection. To identify targets of natural selection in the three-spined stickleback ( Gasterosteus aculeatus ) genome, 103 microsatellite and two indel markers including expressed sequence tags (EST) and quantitative trait loci (QTL)-associated loci, were genotyped in four freshwater and three marine populations. The results indicated that a high proportion of loci (14.7%) might be affected by balancing selection and a lower proportion (2.8%) by directional selection. The strongest signatures of directional selection were detected in a microsatellite locus and two indel markers located in the intronic regions of the Eda-gene coding for the number of lateral plates. Yet, other microsatellite loci previously found to be informative in QTL-mapping studies revealed no signatures of selection. Two novel microsatellite loci ( Stn12 and Stn90 ) located in chromosomes I and VIII, respectively, showed signals of directional selection and might be linked to genomic regions containing gene(s) important for adaptive divergence. Although the coverage of the total genomic content was relatively low, the predominance of balancing selection signals is in agreement with the contention that balancing, rather than directional selection is the predominant mode of selection in the wild.     

Genetic parameters and genome-wide associations of twinning rate in a local breed,the Maremmana cattle

This study seeks to verify the feasibility of increasing twinning in a herd of the Italian autochtonous Maremmana breed. The data set included 1260 individuals born from 1963 to 2014, 527 males and 733 females, 402 of them calving at least once from 1983 through 2015. Breeding values for twinning were estimated by a single-trait linear animal model. However, since twinning is a dichotomous trait and the frequency of twins is far smaller than the frequency of single births, breeding values were also estimated by a single-trait animal threshold model. Heritability of twinning was 0.014±0.018 and 0.062±0.093 for the linear and the threshold models, respectively. Repeatability was 0.071±0.004 and 0.286± 0.012, respectively, for the two models. Genotyping with the Illumina BovineSNP54 BeadChip was performed for cows living on farm in 2012 (119 cows) and a genome-wide association analysis was performed on the corrected phenotype of all calving during the lifespan of each cow, using the GenABEL package in R and a three step GRAMMAR-GC approach. Genomic heritability, calculated from the genomic kinship matrix estimated through genomic marker data, was 0.29±0.021. The most significant detected single nucleotide polymorphisms (Hapmap22923-BTA-129564) was located in proximity of two genes, ARHGAP8 and TMEM200C, which might be potential functional candidates for twinning rate in cattle.     

Effect of breed, gender, housing system and dietary crude protein content on performance of finishing beef cattle fed maize-silage-based diets

Maize silage-based diets with three dietary crude protein (CP) supplements were offered to 96 finishing cattle of contrasting breed (Holstein Friesian (HF) v. Simmental × HF (SHF)) and gender (bull v. steer) housed in two types of feeding system (group fed v. individually fed). The three protein supplements differed either in CP or protein degradability (degradable (LUDP) v. rumen undegradable (HUDP)) and provided CP concentrations of 142 (Con), 175 (LUDP) and 179 (HUDP) g/kg dry matter (DM) respectively, with ratios of degradable to undegradable of 3.0, 1.4 and 0.9:1 for diets Con, LUDP and HUDP, respectively. DM intakes were marginally higher (P = 0.102) for LUDP when compared with Con and HUDP. Rates of daily live-weight gain (DLWG) were higher (P = 0.005) in LUDP and HUDP when compared with Con. HF had higher DM intakes than SHF although this did not result in any improvement in HF DLWG. Bulls had significantly better DM intakes, DLWG and feed conversion efficiency than steers. Conformation scores were better in SHF than HF (P < 0.001) and fat scores lower in bulls than steers (P < 0.001). There was a number of first order interactions established between dietary treatment, breed, gender and housing system with respect to rates of gain and carcass fat scores.     

Predicted accuracy of and response to genomic selection for new traits in dairy cattle

Genomic selection relaxes the requirement of traditional selection tools to have phenotypic measurements on close relatives of all selection candidates. This opens up possibilities to select for traits that are difficult or expensive to measure. The objectives of this paper were to predict accuracy of and response to genomic selection for a new trait, considering that only a cow reference population of moderate size was available for the new trait, and that selection simultaneously targeted an index and this new trait. Accuracy for and response to selection were deterministically evaluated for three different breeding goals. Single trait selection for the new trait based only on a limited cow reference population of up to 10 000 cows, showed that maximum genetic responses of 0.20 and 0.28 genetic standard deviation (s.d.) per year can be achieved for traits with a heritability of 0.05 and 0.30, respectively. Adding information from the index based on a reference population of 5000 bulls, and assuming a genetic correlation of 0.5, increased genetic response for both heritability levels by up to 0.14 genetic s.d. per year. The scenario with simultaneous selection for the new trait and the index, yielded a substantially lower response for the new trait, especially when the genetic correlation with the index was negative. Despite the lower response for the index, whenever the new trait had considerable economic value, including the cow reference population considerably improved the genetic response for the new trait. For scenarios with a zero or negative genetic correlation with the index and equal economic value for the index and the new trait, a reference population of 2000 cows increased genetic response for the new trait with at least 0.10 and 0.20 genetic s.d. per year, for heritability levels of 0.05 and 0.30, respectively. We conclude that for new traits with a very small or positive genetic correlation with the index, and a high positive economic value, considerable genetic response can already be achieved based on a cow reference population with only 2000 records, even when the reliability of individual genomic breeding values is much lower than currently accepted in dairy cattle breeding programs. New traits may generally have a negative genetic correlation with the index and a small positive economic value. For such new traits, cow reference populations of at least 10 000 cows may be required to achieve acceptable levels of genetic response for the new trait and for the whole breeding goal.     

Use of the genomic signature in bacterial classification and identification

In this study we investigated the correlation between dinucleotide relative abundance values (the genomic signature) obtained from bacterial whole-genome sequences and two parameters widely used for bacterial classification, 16S rDNA sequence similarity and DNA-DNA hybridisation values. Twenty-eight completely sequenced bacterial genomes were included in the study. The correlation between the genomic signature and DNA-DNA hybridisation values was high and taxa that showed less than 30% DNA-DNA binding will in general not have dinucleotide relative abundance dissimilarity (delta*) values below 40. On the other hand, taxa showing more than 50% DNA-DNA binding will not have delta* values higher than 17. Our data indicate that the overall correlation between genomic signature and 16S rDNA sequence similarity is low, except for closely related organisms (16S rDNA similarity >94%). Statistical analysis of delta* values between different subgroups of the Proteobacteria indicate that the beta- and gamma-Proteobacteria are more closely related to each other than to the other subgroups of the Proteobacteria and that the alpha- and epsilon-Proteobacteria form clearly separate subgroups. Using the genomic signature we have also predicted DNA-DNA binding values for fastidious or unculturable endosymbionts belonging to the genera Rickettsia, Wigglesworthia and Buchnera.     

Microsatellite signature of ecological selection for salt tolerance in a wild sunflower hybrid species, Helianthus paradoxus

The hybrid sunflower species Helianthus paradoxus inhabits sporadic salt marshes in New Mexico and southwest Texas, USA, whereas its parental species, Helianthus annuus and Helianthus petiolaris, are salt sensitive. Previous studies identified three genomic regions - survivorship quantitative trait loci (QTLs) - that were under strong selection in experimental hybrids transplanted into the natural habitat of H. paradoxus. Here we ask whether these same genomic regions experienced significant selection during the origin and evolution of the natural hybrid, H. paradoxus. This was accomplished by comparing the variability of microsatellites linked to the three survivorship QTLs with those from genomic regions that were neutral in the experimental hybrids. As predicted if one or more selective sweeps had occurred in these regions, microsatellites linked to the survivorship QTLs exhibited a significant reduction in diversity in populations of the natural hybrid species. In contrast, no difference in diversity levels was observed between the two microsatellite classes in parental populations.     

Genetic heterogeneity and selection signature at the KIT gene in pigs showing different coat colours and patterns

Mutations in the porcine KIT gene (Dominant white locus) have been shown to affect coat colours and colour distribution in pigs. We analysed this gene in several pig breeds and populations (Sicilian black, completely black or with white patches; Cinta Senese; grey local population; Large White; Duroc; Hampshire; Pietrain; wild boar; Meishan) with different coat colours and patterns, genotyping a few polymorphisms. The 21 exons and parts of the intronic regions were sequenced in these pigs and 69 polymorphisms were identified. The grey-roan coat colour observed in a local grey population was completely associated with a 4-bp deletion of intron 18 in a single copy KIT gene, providing evidence that this mutation characterizes the I^d allele described in the early genetic literature. The white patches observed in black Sicilian pigs were not completely associated with the presence of a duplicated KIT allele (I^p), suggesting that genetic heterogeneity is a possible cause of different coat colours in this breed. Selection signature was evident at the KIT gene in two different belted pig breeds, Hampshire and Cinta Senese. The same mutation(s) may cause the belted phenotype in these breeds that originated in the 18th–19th centuries from English pigs (Hampshire) and in Tuscany (Italy) in the 14th century (Cinta Senese). Phylogenetic relationships of 28 inferred KIT haplotypes indicated two clades: one of Asian origin that included Meishan and a few Sicilian black haplotypes and another of European origin.     

Identifying footprints of selection in stocked brown trout populations: a spatio‐temporal approach

Studies of interactions between farmed and wild salmonid fishes have suggested reduced fitness of farmed strains in the wild, but evidence for selection at the genic level is lacking. We studied three brown trout populations in Denmark which have been significantly admixed with stocked hatchery trout (19–64%), along with two hatchery strains used for stocking. The wild populations were represented by contemporary samples (2000–2006) and two of them by historical samples (1943–1956). We analysed 61 microsatellite loci, nine of which showed putative functional relationships [expressed sequence tag (EST)‐linked or quantitative trait loci]. F_ST‐based outlier tests provided support for diversifying selection at chromosome regions marked by three loci, two anonymous and one EST‐linked. Patterns of differentiation suggested that the loci were candidates for being under diversifying hitch‐hiking selection in hatchery vs. wild environments. Analysis of hatchery strain admixture proportions showed that in one wild population, two of the loci showed significantly lower admixture proportions than the putatively neutral loci, implying contemporary selection against alleles introduced by hatchery strain trout. In the most strongly admixed population, however, there was no evidence for selection, possibly because of immigration by stocked trout overcoming selection against hatchery‐derived alleles or supportive breeding practices allowing hatchery strain trout to escape natural selection. To our knowledge, this is the first study demonstrating footprints of selection in wild salmonid populations subject to spawning intrusion by farmed fish.     

Genetic association with high‐resolution climate data reveals selection footprints in the genomes of barley landraces across the Iberian Peninsula

Landraces are local populations of crop plants adapted to a particular environment. Extant landraces are surviving genetic archives, keeping signatures of the selection processes experienced by them until settling in their current niches. This study intends to establish relationships between genetic diversity of barley (Hordeum vulgare L.) landraces collected in Spain and the climate of their collection sites. A high‐resolution climatic data set (5 × 5 km spatial, 1‐day temporal grid) was computed from over 2,000 temperature and 7,000 precipitation stations across peninsular Spain. This data set, spanning the period 1981–2010, was used to derive agroclimatic variables meaningful for cereal production at the collection sites of 135 barley landraces. Variables summarize temperature, precipitation, evapotranspiration, potential vernalization and frost probability at different times of the year and time scales (season and month). SNP genotyping of the landraces was carried out combining Illumina Infinium assays and genotyping‐by‐sequencing, yielding 9,920 biallelic markers (7,479 with position on the barley reference genome). The association of these SNPs with agroclimatic variables was analysed at two levels of genetic diversity, with and without taking into account population structure. The whole data sets and analysis pipelines are documented and available at https://eead-csic-compbio.github.io/barley-agroclimatic-association . We found differential adaptation of the germplasm groups identified to be dominated by reactions to cold temperature and late‐season frost occurrence, as well as to water availability. Several significant associations pointing at specific adaptations to agroclimatic features related to temperature and water availability were observed, and candidate genes underlying some of the main regions are proposed.     

Evaluation of the potential use of a meta-population for genomic selection in autochthonous beef cattle populations

This study investigated the potential application of genomic selection under a multi-breed scheme in the Spanish autochthonous beef cattle populations using a simulation study that replicates the structure of linkage disequilibrium obtained from a sample of 25 triplets of sire/dam/offspring per population and using the BovineHD Beadchip. Purebred and combined reference sets were used for the genomic evaluation and several scenarios of different genetic architecture of the trait were investigated. The single-breed evaluations yielded the highest within-breed accuracies. Across breed accuracies were found low but positive on average confirming the genetic connectedness between the populations. If the same genotyping effort is split in several populations, the accuracies were lower when compared with single-breed evaluation, but showed a small advantage over small-sized purebred reference sets over the accuracies of subsequent generations. Besides, the genetic architecture of the trait did not show any relevant effect on the accuracy with the exception of rare variants, which yielded slightly lower results and higher loss of predictive ability over the generations.