The Genome Response to Artificial Selection: A Case Study in Dairy Cattle

Dairy cattle breeds have been subjected over the last fifty years to intense artificial selection towards improvement of milk production traits. In this study, we performed a whole genome scan for differentiation using 42,486 SNPs in the three major French dairy cattle breeds (Holstein, Normande and Montbéliarde) to identify the main physiological pathways and regions which were affected by this selection. After analyzing the population structure, we estimated F_ST within and across the three breeds for each SNP under a pure drift model. We further considered two different strategies to evaluate the effect of selection at the genome level. First, smoothing F_ST values over each chromosome with a local variable bandwidth kernel estimator allowed identifying 13 highly significant regions subjected to strong and/or recent positive selection. Some of them contained genes within which causal variants with strong effect on milk production traits (GHR) or coloration (MC1R) have already been reported. To go further in the interpretation of the observed signatures of selection we subsequently concentrated on the annotation of differentiated genes defined according to the F_ST value of SNPs localized close or within them. To that end we performed a comprehensive network analysis which suggested a central role of somatotropic and gonadotropic axes in the response to selection. Altogether, these observations shed light on the antagonism, at the genome level, between milk production and reproduction traits in highly producing dairy cows.     

Analysis of population differentiation in North Eurasian cattle (Bos taurus) using single nucleotide polymorphisms in three genes associated with production traits

Single nucleotide polymorphisms (SNPs) in growth hormone 1 (GH1), insulin-like growth factor 1 (IGF1) and leptin (LEP), all candidates for production traits in cattle, were characterized in North Eurasian cattle breeds. Allele frequencies of IGF1 exhibited significant (P < 0.05) deviation from neutral expectation and therefore, might be associated with divergence in North Eurasian cattle because of genetic selection. Allele frequencies and lower heterozygosity of LEP may indicate a recent introduction of an alternative allele in this geographic region. Locus F(ST) estimates were highest for IGF1 (0.151, sigma = 0.042) and lowest for GH (0.062, sigma = 0.020). Our results suggest a slightly higher population differentiation across the candidate genes (FST = 0.108) than across microsatellites (FST = 0.095), possibly because of selection and stochastic effects.     

Design of a bovine low-density SNP array optimized for imputation

The Illumina BovineLD BeadChip was designed to support imputation to higher density genotypes in dairy and beef breeds by including single-nucleotide polymorphisms (SNPs) that had a high minor allele frequency as well as uniform spacing across the genome except at the ends of the chromosome where densities were increased. The chip also includes SNPs on the Y chromosome and mitochondrial DNA loci that are useful for determining subspecies classification and certain paternal and maternal breed lineages. The total number of SNPs was 6,909. Accuracy of imputation to Illumina BovineSNP50 genotypes using the BovineLD chip was over 97% for most dairy and beef populations. The BovineLD imputations were about 3 percentage points more accurate than those from the Illumina GoldenGate Bovine3K BeadChip across multiple populations. The improvement was greatest when neither parent was genotyped. The minor allele frequencies were similar across taurine beef and dairy breeds as was the proportion of SNPs that were polymorphic. The new BovineLD chip should facilitate low-cost genomic selection in taurine beef and dairy cattle.     

Genome-wide linkage disequilibrium in two Japanese beef cattle breeds

There is little knowledge about the degree of linkage disequilibrium (LD) in beef cattle. This study aims to perform a genome-wide search for LD in Japanese Black and Japanese Brown beef cattle and to compare the level of LD between these two breeds. Parameter D' (the LD coefficient) was used as a measure of LD, and LD was tested for significance of allelic associations between syntenic and between non-syntenic marker pairs. Effects of breed, chromosome, genetic map distance and their interactions with D' were tested based on least squares analyses. Both breeds showed high levels of LD, which ranged over several tens of cM and declined as the marker distance increased for syntenic marker pairs. A rapid decline of the D' value was observed between markers that were spaced 5 and 20 cM apart. LD was significant in most cases for marker pairs <40 cM apart but was not significant between non-syntenic loci. The pattern of LD found in these two breeds was similar to that previously published for dairy cattle. The D' value between breeds was not significantly different (P > 0.05), but the interaction between breed and chromosome was highly significant (P < 0.001). Genetic selection seems to have caused the heterogeneity of the D' values among chromosomes within breed. These results indicate that LD mapping is a useful tool for fine-mapping quantitative trait loci of economically important traits in Japanese beef cattle.     

Comparison of milk protein allele frequencies in Nordic cattle breeds

Allele frequencies at four milk protein loci were studied in five modern and 17 old Nordic cattle breeds in order to reveal variants that are characteristic for these populations. The B allele of CSN3, which has been associated with improved manufacturing properties of milk, showed significantly lower frequencies in modern production breeds than in old breeds of interest for conservation purposes. Characteristic frequencies of CSN1S1 (C), CSN2 (A2) and CSN3 (B) were found in Icelandic cattle, Swedish Mountain cattle, Northern Finncattle and Western Fjord cattle, which indicate a common origin of these populations. Further comparisons of allele frequencies in old Nordic breeds suggest sorting of these breeds into two groups with a northern and southern geographic location.     

Whole-Genome Resequencing Analysis of Hanwoo and Yanbian Cattle to Identify Genome-Wide SNPs and Signatures of Selection

Over the last 30 years, Hanwoo has been selectively bred to improve economically important traits. Hanwoo is currently the representative Korean native beef cattle breed, and it is believed that it shared an ancestor with a Chinese breed, Yanbian cattle, until the last century. However, these two breeds have experienced different selection pressures during recent decades. Here, we whole-genome sequenced 10 animals each of Hanwoo and Yanbian cattle (20 total) using the Illumina HiSeq 2000 sequencer. A total of approximately 3.12 and 3.07 billion sequence reads were mapped to the bovine reference sequence assembly (UMD 3.1) at an average of approximately 10.71- and 10.53-fold coverage for Hanwoo and Yanbian cattle, respectively. A total of 17,936,399 single nucleotide polymorphisms (SNPs) were yielded, of which 22.3% were found to be novel. By annotating the SNPs, we further retrieved numerous nonsynonymous SNPs that may be associated with traits of interest in cattle. Furthermore, we performed whole-genome screening to detect signatures of selection throughout the genome. We located several promising selective sweeps that are potentially responsible for economically important traits in cattle; the PPP1R12A gene is an example of a gene that potentially affects intramuscular fat content. These discoveries provide valuable genomic information regarding potential genomic markers that could predict traits of interest for breeding programs of these cattle breeds.     

Identification of contemporary selection signatures using composite log likelihood and their associations with marbling score in Korean cattle

Positive selection not only increases beneficial allele frequency but also causes augmentation of allele frequencies of sequence variants in close proximity. Signals for positive selection were detected by the statistical differences in subsequent allele frequencies. To identify selection signatures in Korean cattle, we applied a composite log‐likelihood (CLL)‐based method, which calculates a composite likelihood of the allelic frequencies observed across sliding windows of five adjacent loci and compares the value with the critical statistic estimated by 50 000 permutations. Data for a total of 11 799 nucleotide polymorphisms were used with 71 Korean cattle and 209 foreign beef cattle. As a result, 147 signals were identified for Korean cattle based on CLL estimates (P < 0.01). The signals might be candidate genetic factors for meat quality by which the Korean cattle have been selected. Further genetic association analysis with 41 intragenic variants in the selection signatures with the greatest CLL for each chromosome revealed that marbling score was associated with five variants. Intensive association studies with all the selection signatures identified in this study are required to exclude signals associated with other phenotypes or signals falsely detected and thus to identify genetic markers for meat quality.     

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.     

奶牛重要经济性状全基因组关联分析研究进展

产奶性状是奶牛最重要的生产性状,随着平衡育种理念的提出和发展,繁殖性状、体型性状、健康性状和长寿性等功能性状也逐渐被重视并纳入育种规划中。鉴定产奶性状和功能性状主效基因或遗传标记并将之应用于奶牛标记辅助选择可望加快遗传进展。随着高密度SNP标记的高通量检测技术的发展,全基因组关联分析已成为鉴定畜禽重要经济性状基因的重要途径。文章对奶牛产奶性状和功能性状全基因组关联分析研究进展进行综述。     

奶牛重要经济性状全基因组关联分析研究进展

产奶性状是奶牛最重要的生产性状, 随着平衡育种理念的提出和发展, 繁殖性状、体型性状、健康性状和长寿性等功能性状也逐渐被重视并纳入育种规划中。鉴定产奶性状和功能性状主效基因或遗传标记并将之应用于奶牛标记辅助选择可望加快遗传进展。随着高密度SNP标记的高通量检测技术的发展, 全基因组关联分析已成为鉴定畜禽重要经济性状基因的重要途径。文章对奶牛产奶性状和功能性状全基因组关联分析研究进展进行综述。     

Genome-wide identification of runs of homozygosity islands and associated genes in local dairy cattle breeds

Runs of homozygosity (ROH) are widely used as predictors of whole-genome inbreeding levels in cattle. They identify regions that have an unfavorable effect on a phenotype when homozygous, but also identify the genes associated with traits of economic interest present in these regions. Here, the distribution of ROH islands and enriched genes within these regions in four dairy cattle breeds were investigated. Cinisara (71), Modicana (72), Reggiana (168) and Italian Holstein (96) individuals were genotyped using the 50K v2 Illumina BeadChip. The genomic regions most commonly associated with ROHs were identified by selecting the top 1% of the single nucleotide polymorphisms (SNPs) most commonly observed in the ROH of each breed. In total, 11 genomic regions were identified in Cinisara and Italian Holstein, and eight in Modicana and Reggiana, indicating an increased ROH frequency level. Generally, ROH islands differed between breeds. The most homozygous region (>45% of individuals with ROH) was found in Modicana on chromosome 6 within a quantitative trail locus affecting milk fat and protein concentrations. We identified between 126 and 347 genes within ROH islands, which are involved in multiple signaling and signal transduction pathways in a wide variety of biological processes. The gene ontology enrichment provided information on possible molecular functions, biological processes and cellular components under selection related to milk production, reproduction, immune response and resistance/susceptibility to infection and diseases. Thus, scanning the genome for ROH could be an alternative strategy to detect genomic regions and genes related to important economic traits.     

Recent artificial selection in U.S. Jersey cattle impacts autozygosity levels of specific genomic regions

Background

Genome signatures of artificial selection in U.S. Jersey cattle were identified by examining changes in haplotype homozygosity for a resource population of animals born between 1953 and 2007. Genetic merit of this population changed dramatically during this period for a number of traits, especially milk yield. The intense selection underlying these changes was achieved through extensive use of artificial insemination (AI), which also increased consanguinity of the population to a few superior Jersey bulls. As a result, allele frequencies are shifted for many contemporary animals, and in numerous cases to a homozygous state for specific genomic regions. The goal of this study was to identify those selection signatures that occurred after extensive use of AI since the 1960, using analyses of shared haplotype segments or Runs of Homozygosity. When combined with animal birth year information, signatures of selection associated with economically important traits were identified and compared to results from an extended haplotype homozygosity analysis.

Results

Overall, our results reveal that more recent selection increased autozygosity across the entire genome, but some specific regions increased more than others. A genome-wide scan identified more than 15 regions with a substantial change in autozygosity. Haplotypes found to be associated with increased milk, fat and protein yield in U.S. Jersey cattle also consistently increased in frequency.

Conclusions

The analyses used in this study was able to detect directional selection over the last few decades when individual production records for Jersey animals were available.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1500-x) contains supplementary material, which is available to authorized users.     

Use of locally weighted scatterplot smoothing (LOWESS) regression to study selection signatures in Piedmontese and Italian Brown cattle breeds

Selection is the major force affecting local levels of genetic variation in species. The availability of dense marker maps offers new opportunities for a detailed understanding of genetic diversity distribution across the animal genome. Over the last 50 years, cattle breeds have been subjected to intense artificial selection. Consequently, regions controlling traits of economic importance are expected to exhibit selection signatures. The fixation index (F_st) is an estimate of population differentiation, based on genetic polymorphism data, and it is calculated using the relationship between inbreeding and heterozygosity. In the present study, locally weighted scatterplot smoothing (LOWESS) regression and a control chart approach were used to investigate selection signatures in two cattle breeds with different production aptitudes (dairy and beef). F_st was calculated for 42 514 SNP marker loci distributed across the genome in 749 Italian Brown and 364 Piedmontese bulls. The statistical significance of F_st values was assessed using a control chart. The LOWESS technique was efficient in removing noise from the raw data and was able to highlight selection signatures in chromosomes known to harbour genes affecting dairy and beef traits. Examples include the peaks detected for BTA2 in the region where the myostatin gene is located and for BTA6 in the region harbouring the ABCG2 locus. Moreover, several loci not previously reported in cattle studies were detected.     

Discovery, characterization and validation of single nucleotide polymorphisms within 206 bovine genes that may be considered as candidate genes for beef production and quality

A large number of putative single nucleotide polymorphisms (SNPs) have been identified from the bovine genome-sequencing project. However, few of these have been validated and many will turn out to be sequencing artefacts or have low minor allele frequencies. In addition, there is little information available on SNPs within coding regions, which are likely to be responsible for phenotypic variation. Therefore, additional SNP discovery is necessary to identify and validate polymorphisms both in specific genes and genome-wide. Sequence-tagged sites within 286 genes were resequenced from a panel of animals representing a wide range of European cattle breeds. For 80 genes, no polymorphisms were identified, and 672 putative SNPs were identified within 206 genes. Fifteen European cattle breeds (436 individuals plus available parents) were genotyped with these putative SNPs, and 389 SNPs were confirmed to have minor allele frequencies above 10%. The genes containing SNPs were localized on chromosomes by radiation hybrid mapping and on the bovine genome sequence by Blast . Flanking microsatellite loci were identified, to facilitate the alignment of the genes containing the SNPs in relation to mapped quantitative trait loci. Of the 672 putative SNPs discovered in this work, only 11 were found among the validated SNPs and 100 were found among the approximately 2.3 million putative SNPs currently in dbSNP. The genes studied in this work could be considered as candidates for traits associated with beef production and the SNPs reported will help to assess the role of the genes in the genetic control of muscle development and meat quality. The allele frequency data presented allows the general utility of the SNPs to be assessed.     

Comparative selection signature analyses identify genomic footprints in Reggiana cattle,the traditional breed of the Parmigiano-Reggiano cheese production system

Reggiana is an autochthonous cattle breed reared mainly in the province of Reggio Emilia, located in the North of Italy. Reggiana cattle (originally a triple-purpose population largely diffused in the North of Italy) are characterised by a typical solid red coat colour. About 2500 cows of this breed are currently registered to its herd book. Reggiana is now considered a dual-purpose breed even if it is almost completely dedicated to the production of a mono-breed branded Protected Designation of Origin Parmigiano-Reggiano cheese, which is the main driver of the sustainable conservation of this local genetic resource. In this study, we provided the first overview of genomic footprints that characterise Reggiana and define the diversity of this local cattle breed. A total of 168 Reggiana sires (all bulls born over 35 years for which semen was available) and other 3321 sires from 3 cosmopolitan breeds (Brown, Holstein and Simmental) were genotyped with the Illumina BovineSNP50 panel. ADMIXTURE analysis suggested that Reggiana breed might have been influenced, at least in part, by the other three breeds included in this study. Selection signatures in the Reggiana genome were identified using three statistical approaches based on allele frequency differences among populations or on properties of haplotypes segregating in the populations (fixation index (F_ST); integrated haplotype score; cross-population extended haplotype homozygosity). We identified several regions under peculiar selection in the Reggiana breed, particularly on bovine chromosome (BTA) 6 in the KIT gene region, that is known to be involved in coat colour pattern distribution, and within the region of the LAP3, NCAPG and LCORL genes, that are associated with stature, conformation and carcass traits. Another already known region that includes the PLAG1 gene (BTA14), associated with conformation traits, showed a selection signature in the Reggiana cattle. On BTA18, a signal of selection included the MC1R gene that causes the red coat colour in cattle. Other selection sweeps were in regions, with high density of quantitative trait loci for milk production traits (on BTA20) and in several other large regions that might have contributed to shape and define the Reggiana genome (on BTA17 and BTA29). All these results, overall, indicate that the Reggiana genome might still contain several signs of its multipurpose and non-specialised utilisation, as already described for other local cattle populations, in addition to footprints derived by its ancestral origin and by its adaptation to the specialised Parmigiano-Reggiano cheese production system.     

Detection of selection signatures within candidate regions underlying trypanotolerance in outbred cattle populations

Breeding indigenous African taurine cattle tolerant to trypanosomosis is a straightforward approach to control costs generated by this disease. A recent study identified quantitative trait loci (QTL) underlying trypanotolerance traits in experimental crosses between tolerant N'Dama and susceptible Boran zebu cattle. As trypanotolerance is thought to result from local adaptation of indigenous cattle breeds, we propose an alternative and complementary approach to study the genetic architecture of this trait, based on the identification of selection signatures within QTL or candidate genes. A panel of 92 microsatellite markers was genotyped on 509 cattle belonging to four West African trypanotolerant taurine breeds and 10 trypanosusceptible European or African cattle breeds. Some of these markers were located within previously identified QTL regions or candidate genes, while others were chosen in regions assumed to be neutral. A detailed analysis of the genetic structure of these different breeds was carried out to confirm a priori grouping of populations based on previous data. Tests based on the comparison of the observed heterozygosities and variances in microsatellite allelic size among trypanotolerant and trypanosusceptible breeds led to the identification of two significantly less variable microsatellite markers. BM4440, one of these two outlier loci, is located within the confidence interval of a previously described QTL underlying a trypanotolerance-related trait.
Detection of selection signatures appears to be a straightforward approach for unravelling the molecular determinism of trypanosomosis pathogenesis. We expect that a whole genome approach will help confirm these results and achieve a higher resolving power.     

Assessing signatures of selection through variation in linkage disequilibrium between taurine and indicine cattle

Background

Signatures of selection are regions in the genome that have been preferentially increased in frequency and fixed in a population because of their functional importance in specific processes. These regions can be detected because of their lower genetic variability and specific regional linkage disequilibrium (LD) patterns.

Methods

By comparing the differences in regional LD variation between dairy and beef cattle types, and between indicine and taurine subspecies, we aim at finding signatures of selection for production and adaptation in cattle breeds. The VarLD method was applied to compare the LD variation in the autosomal genome between breeds, including Angus and Brown Swiss, representing taurine breeds, and Nelore and Gir, representing indicine breeds. Genomic regions containing the top 0.01 and 0.1 percentile of signals were characterized using the UMD3.1 Bos taurus genome assembly to identify genes in those regions and compared with previously reported selection signatures and regions with copy number variation.

Results

For all comparisons, the top 0.01 and 0.1 percentile included 26 and 165 signals and 17 and 125 genes, respectively, including TECRL, BT.23182 or FPPS, CAST, MYOM1, UVRAG and DNAJA1.

Conclusions

The VarLD method is a powerful tool to identify differences in linkage disequilibrium between cattle populations and putative signatures of selection with potential adaptive and productive importance.     

Genome‐wide variability and selection signatures in Italian island cattle breeds

In the present study, a sample of 88 animals belonging to four local (Modicana, Sarda, Sardo‐Bruna and Sardo‐Modicana) and one cosmopolitan (Italian Brown Swiss) cattle breeds were genotyped with a medium density SNP beadchip and compared to investigate their genetic diversity and the existence of selection signatures. A total of 43 012 SNPs distributed across all 29 autosomal chromosomes were retained after data quality control. Basic population statistics, Wright fixation index and runs of homozygosity (ROH) analyses confirmed that the Italian Brown Swiss genome was shaped mainly by selection, as underlined by the low values of heterozygosity and minor allele frequency. As expected, local cattle exhibited a large within‐breed genetic heterogeneity. The F_ST comparison revealing the largest number of significant SNPs was Sardo‐Bruna vs. Sardo‐Modicana, whereas the smallest was observed for Italian Brown Swiss vs. Sardo‐Modicana. Modicana exhibited the largest number of detected ROHs, whereas the smallest was observed for Sardo‐Modicana. Signatures of selection were detected in genomic regions that harbor genes involved in milk production traits for Italian Brown Swiss and fitness traits for local breeds. According to the results of multi‐dimensional scaling and the admixture analysis the Sardo‐Bruna is more similar to the Sarda than to the Italian Brown Swiss breed. Moreover, the Sardo‐Modicana is genetically closer to the Modicana than to the Sarda breed. Results of the present work confirm the usefulness of single nucleotide polymorphisms in deciphering the genetic architecture of livestock breeds.     

Genomic organization of the DGAT2/MOGAT gene family in cattle (Bos taurus) and other mammals

We report the cloning and initial characterization of the genes encoding DGAT2 (diacylglycerol transferase 2), MOGAT1 and MOGAT2 (monoacylglycerol transferases 1 and 2) in domestic cattle (Bos taurus). The three closely related genes belong to a gene family with at least eight members in mammals and are candidate genes for quantitative traits related to dietary fat uptake, lipid synthesis and storage. MOGAT2 and DGAT2 form a tandem and were mapped to bovine chromosome (BTA) 15q25-->q26 by fluorescence in situ hybridization. MOGAT1 was localized to BTA 2q43-->q44. The three genes were investigated for polymorphisms that might be associated with breeding values for milk fat percentage in the dairy breeds German Holstein, German Simmental and German Brown. All the detected polymorphisms were located outside exons or, with one exception, were silent. In MOGAT1, a missense mutation in exon 4 was found that causes a non-conservative substitution of cysteine170 (uncharged, hydrophobic) by lysine (positively charged, hydrophilic). However, allele frequency estimates from pooled DNA samples revealed no significant association of the observed polymorphisms with breeding values for milk fat percentage. A comparative analysis of chromosomal locations and exon-intron structure of the known members of the DGAT2/MOGAT gene family in humans, rodents and cattle indicates an ancient tandem duplication of the ancestor gene combined with an intron gain (or loss) in one copy. Further members of the family may have arisen by duplications of this gene tandem via two rounds of interchromosomal or genome duplications as well as further local (single) gene duplication and loss events.