A post-GWAS confirming the SCD gene associated with milk medium- and long-chain unsaturated fatty acids in Chinese Holstein population.

The stearoyl-CoA desaturase (delta-9-desaturase) gene encodes a key enzyme in the cellular biosynthesis of monounsaturated fatty acids. In our initial genome-wide association study (GWAS) of Chinese Holstein cows, 19 SNPs fell in a 1.8-Mb region (20.3-22.1 Mb) on chromosome 26 underlying the SCD gene and were highly significantly associated with C14:1 or C14 index. The aims of this study were to verify whether the SCD gene has significant genetic effects on milk fatty acid composition in dairy cattle. By resequencing the entire coding region of the bovine SCD gene, a total of six variations were identified, including three coding variations (g.10153G>A, g.10213T>C and g.10329C>T) and three intronic variations (g.6926A>G, g.8646G>A and g.16158G>C). The SNP in exon 3, g.10329C>T, was predicted to result in an amino acid replacement from alanine (GCG) to valine (GTG) in the SCD protein. An association study for 16 milk fatty acids using 346 Chinese Holstein cows with accurate phenotypes and genotypes was performed using the mixed animal model with the proc mixed procedure in sas 9.2. All six detected SNPs were revealed to be associated with six medium- and long-chain unsaturated fatty acids (P = 0.0457 to P < 0.0001), specifically for C14:1 and C14 index (P = 0.0005 to P < 0.0001). Subsequently, strong linkage disequilibrium (D' = 0.88-1.00) was observed among all six SNPs in SCD and the five SNPs (rs41623887, rs109923480, rs42090224, rs42092174 and rs42091426) within the 1.8-Mb region identified in our previous GWAS, indicating that the significant association of the SCD gene with milk fatty acid content traits reduced the observed significant 1.8-Mb chromosome region in GWAS. Haplotype-based analysis revealed significant associations of the haplotypes encompassing the six SCD SNPs and one SNP (rs109923480) in a GWAS with C14:1, C14 index, C16:1 and C16 index (P = 0.0011 to P < 0.0001). In summary, our findings provide replicate evidence for our previous GWAS and demonstrate that variants in the SCD gene are significantly associated with milk fatty acid composition in dairy cattle, which provides clear evidence for an increased understanding of milk fatty acid synthesis and enhances opportunities to improve milk-fat composition in dairy cattle.

Li C ，Sun D ，Zhang S ，Liu L ，Alim MA ，Zhang Q ... -  《-》

Genetic effects of stearoyl-coenzyme A desaturase (SCD) polymorphism on milk production traits in the Chinese dairy population.

Stearoyl-CoA desaturase (SCD) is a multifunctional complex enzyme important in the cellular biosynthesis of fatty acids. The present study was to investigate the association of the SCD gene with milk production traits in dairy cattle. Two single nucleotide polymorphisms (SNPs) (g.6926A>G and g.8646A>G) in introns 3 and 4, and three SNPs (g.10153A>G, g.10213T>C and g.10329C>T) in exon 5 were identified with pooled DNA sequencing and genotyped using matrix-assisted laser desorption/ionization time of flight mass spectrometry assay in 752 Chinese Holstein cows. Polymorphism g.10329C>T was predicted to result in an amino acid replacement from alanine to valine in the SCD protein. With a mixed animal model, the significant associations of the five SNPs with 305-day milk, fat and protein yields and protein percentage were determined. We further demonstrated cows with heterozygous genotypes (A/G or C/T) had highest 305 day milk yield, fat yield, protein yield and lowest protein percentage. Heterozygous cows with genotype AG at the g.6926A>G locus showed the greatest milk yield (P < 0.0001), fat yield (P < 0.0001) and protein yield (P < 0.0001) among other heterozygous genotypes at any of the loci. Dominance effects of all identified SNPs on milk, fat and protein yields and protein percentage were significant. Moreover, significant allele substitution effects at g.6926A>G locus on milk yield and at g.10213T>C on protein yield were observed. Five-locus haplotypes and strong linkage disequilibrium (D' > 0.9) between the five SNPs were also observed. The results suggest that identified polymorphisms could be potential genetic markers to improve the production performance of Chinese Holstein.

Alim MA ，Fan YP ，Wu XP ，Xie Y ，Zhang Y ，Zhang SL ，Sun DX ，Zhang Y ，Zhang Q ，Liu L ，Guo G ... -  《-》

Genetic effects of FASN, PPARGC1A, ABCG2 and IGF1 revealing the association with milk fatty acids in a Chinese Holstein cattle population based on a post genome-wide association study.

A previous genome-wide association study deduced that one (ARS-BFGL-NGS-39328), two (Hapmap26001-BTC-038813 and Hapmap31284-BTC-039204), two (Hapmap26001-BTC-038813 and BTB-00246150), and one (Hapmap50366-BTA-46960) genome-wide significant single nucleotide polymorphisms (SNPs) associated with milk fatty acids were close to or within the fatty acid synthase (FASN), peroxisome proliferator-activated receptor gamma, coactivator 1 alpha (PPARGC1A), ATP-binding cassette, sub-family G, member 2 (ABCG2) and insulin-like growth factor 1 (IGF1) genes. To further confirm the linkage and reveal the genetic effects of these four candidate genes on milk fatty acid composition, genetic polymorphisms were identified and genotype-phenotype associations were performed in a Chinese Holstein cattle population. Nine SNPs were identified in FASN, among which SNP rs41919985 was predicted to result in an amino acid substitution from threonine (ACC) to alanine (GCC), five SNPs (rs136947640, rs134340637, rs41919992, rs41919984 and rs41919986) were synonymous mutations, and the remaining three (rs41919999, rs132865003 and rs133498277) were found in FASN introns. Only one SNP each was identified for PPARGC1A, ABCG2 and IGF1. Association studies revealed that FASN, PPARGC1A, ABCG2 and IGF1 were mainly associated with medium-chain saturated fatty acids and long-chain unsaturated fatty acids, especially FASN for C10:0, C12:0 and C14:0. Strong linkage disequilibrium was observed among ARS-BFGL-NGS-39328 and rs132865003 and rs134340637 in FASN (D´ > 0.9), and among Hapmap26001-BTC-038813 and Hapmap31284-BTC-039204 and rs109579682 in PPARGC1A (D´ > 0.9). Subsequently, haplotype-based analysis revealed significant associations of the haplotypes encompassing eight FASN SNPs (rs41919999, rs132865003, rs134340637, rs41919992, rs133498277, rs41919984, rs41919985 and rs41919986) with C10:0, C12:0, C14:0, C18:1n9c, saturated fatty acids (SFA) and unsaturated fatty acids (UFA) (P = 0.0204 to P < 0.0001). Our study confirmed the linkage between the significant SNPs in our previous genome-wide association study and variants in FASN and PPARGC1A. SNPs within FASN, PPARGC1A, ABCG2 and IGF1 showed significant genetic effects on milk fatty acid composition in dairy cattle, indicating their potential functions in milk fatty acids synthesis and metabolism. The findings presented here provide evidence for the selection of dairy cows with healthier milk fatty acid composition by marker-assisted breeding or genomic selection schemes, as well as furthering our understanding of technological processing aspects of cows' milk.

Li C ，Sun D ，Zhang S ，Yang S ，Alim MA ，Zhang Q ，Li Y ，Liu L ... -  《BMC GENETICS》

Genome-wide association and biological pathway analysis for milk-fat composition in Danish Holstein and Danish Jersey cattle.

The milk fat profile of the Danish Holstein (DH) and Danish Jersey (DJ) show clear differences. Identification of the genomic regions, genes and biological pathways underlying the milk fat biosynthesis will improve the understanding of the biology underlying bovine milk fat production and may provide new possibilities to change the milk fat composition by selective breeding. In this study a genome wide association scan (GWAS) in the DH and DJ was performed for a detailed milk fatty acid (FA) profile using the HD bovine SNP array and subsequently a biological pathway analysis based on the SNP data was performed. The GWAS identified in total 1,233 SNPs (FDR < 0.10) spread over 18 chromosomes for nine different FA traits for the DH breed and 1,122 SNPs (FDR < 0.10) spread over 26 chromosomes for 13 different FA traits were detected for the DJ breed. Of these significant SNPs, 108 SNP markers were significant in both DH and DJ (C14-index, BTA26; C16, BTA14; fat percentage (FP), BTA14). This was supported by an enrichment test. The QTL on BTA14 and BTA26 represented the known candidate genes DGAT and SCD. In addition we suggest ACSS3 to be a good candidate gene for the QTL on BTA5 for C10:0 and C15:0. In addition, genetic correlations between the FA traits within breed showed large similarity across breeds. Furthermore, the biological pathway analysis revealed that fat digestion and absorption (KEGG04975) plays a role for the traits FP, C14:1, C16 index and C16:1. There was a clear similarity between the underlying genetics of FA in the milk between DH and DJ. This was supported by the fact that there was substantial overlap between SNPs for FP, C14 index, C14:1, C16 index and C16:1. In addition genetic correlations between FA showed a similar pattern across DH and DJ. Furthermore the biological pathway analysis suggested that fat digestion and absorption KEGG04975 is important for the traits FP, C14:1, C16 index and C16:1.

Buitenhuis B ，Janss LL ，Poulsen NA ，Larsen LB ，Larsen MK ，Sørensen P ... -  《BMC GENOMICS》

Variation of Delta 9-desaturase activity in dairy cattle.

The endogenous production of unsaturated fatty acids (FA), particularly some monounsaturated FA (%MONO) and nearly all conjugated linoleic acids, is regulated by the Delta(9)-desaturase activity. The aims of this study were to assess the variation of this enzymatic activity within lactation, across dairy breeds, and to estimate its genetic parameters. The ratios of C14:1 cis-9 to C14:0, C16:1 cis-9 to C16:0, and C18:1 cis to C18:0 were calculated from FA contents predicted by mid-infrared spectrometry. Variance components and standard errors were estimated using average information REML. The multitrait mixed model included as fixed effects herd x test date x class of lactation number, class of days in milk x class of lactation number, class of age x class of lactation number, and regressions on breed composition. Four random effects were also included: animal genetic effect, 2 permanent environments (within and across lactations), and residual effect. Under the assumption that the calculated ratios are an approximate measurement of Delta(9)-desaturase activity, this study showed different sources of variation for this enzymatic activity. A slight difference was observed within lactation. The ratios of C14:1 cis-9 to C14:0 and C16:1 cis-9 to C16:0 increased as a function of days in milk. Differences across 7 dairy breeds were observed. The values of Delta(9)-desaturase indices observed for Jersey and Brown-Swiss cows were lower compared with Holstein. The opposite was observed for dual-purpose Belgian Blue cows. Values of heritability for the ratios of C14:1 cis-9 to C14:0, C16:1 cis-9 to C16:0, and C18:1 cis to C18:0 were 20, 20, and 3%, respectively. Negative genetic correlations observed between fat or protein contents and the 3 indices suggested that an increased activity of Delta(9)-desaturase could inhibit the synthesis of fat and protein in bovine milk. Negative correlations were also observed between fat or protein contents and the contents of 3 studied unsaturated FA in milk fat (C14:1 cis-9, C16:1 cis-9, and C18:1 cis). The positive genetic correlations observed between %MONO and the ratios of C14:1 cis-9 to C14:0 (0.72), C16:1 cis-9 to C16:0 (0.62), and C18:1 cis to C18:0 (0.97) showed that %MONO is linked to the Delta(9)-desaturase activity.

Soyeurt H ，Dehareng F ，Mayeres P ，Bertozzi C ，Gengler N ... -  《-》