Genome-wide association study of milk fatty acid composition in Italian Simmental and Italian Holstein cows using single nucleotide polymorphism arrays.

Bovine milk is important for human nutrition, but its fat content is often criticized as a risk factor in cardiovascular disease. Selective breeding programs could be used to alter the fatty acid (FA) composition of bovine milk to improve the healthiness of dairy products for human consumption. Here, we performed a genome-wide association study (GWAS) on bovine milk to identify genomic regions or specific genes associated with FA profile and to investigate genetic differences between the Italian Simmental (IS) and Italian Holstein (IH) breeds. To achieve this, we first characterized milk samples from 416 IS cows and 436 IH cows for their fat profile by gas chromatography. Subjects were genotyped with single nucleotide polymorphism array and a single-marker regression model for GWAS was performed. Our findings confirm previously reported quantitative trait loci strongly associated with bovine milk fat composition. More specifically, our GWAS results revealed significant signals on chromosomes Bos taurus autosome 19 and 26 for milk FA. Further analysis using a gene-centric approach and pathway meta-analysis identified not only some well-known genes underlying quantitative trait loci for milk FA components, such as FASN, SCD, and DGAT1, but also other significant candidate genes, including some with functional roles in pathways related to "Lipid metabolism." Highlighted genes related to FA profile include ECI2, PCYT2, DCXR, G6PC3, PYCR1, and ALG12 in IS, and CYP17A1, ACO2, PI4K2A, GOT1, GPT, NT5C2, PDE6G, POLR3H, and COX15 in IH. Overall, the breed-specific association outcomes reflect differences in the genetic backgrounds of the IS and IH breeds and their selective breeding histories.

Palombo V ，Milanesi M ，Sgorlon S ，Capomaccio S ，Mele M ，Nicolazzi E ，Ajmone-Marsan P ，Pilla F ，Stefanon B ，D'Andrea M ... -  《-》

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》

Fine-mapping of BTA17 using imputed sequences for associations with de novo synthesized fatty acids in bovine milk.

The aim of this study was to fine-map a genomic region associated with milk fatty acids (FA) on Bos taurus autosome (BTA) 17. This genomic region has been discovered with 50,000 (50k) single nucleotide polymorphisms (SNP) imputed to 777,000 (777k) SNP. In this study, high-density genotypes were imputed to whole-genome sequences level to identify candidate gene(s) associated with milk FA composition on BTA17. Phenotypes and genotypes were available for 1,640 cows sampled in winter, and for 1,581 cows sampled in summer. Phenotypes consisted of gas chromatography measurements in winter and in summer milk samples of 6 individual FA and the indicator of de novo synthesis, C6:0-C14:0. Genotypes consisted of imputed 777k SNP, and 89 sequenced ancestors of the population of genotyped cows. In addition, 450 whole-genome sequences from the 1,000 Bull Genome Consortium were available. Using 495 Holstein-Friesian sequences as a reference population, the 777k SNP genotypes of the cows were imputed to sequence level. We then applied single-variant analyses with an animal model, and identified thousands of significant associations with C6:0, C8:0, C10:0, C12:0, C14:0, and C6:0-C14:0. For C8:0 in summer milk samples, the genomic region located between 29 and 34 Mbp on BTA17 revealed a total of 646 significant associations. The most significant associations [-log(P-value) = 7.82] were 8 SNP in perfect linkage disequilibrium. After fitting one of these 8 SNP as a fixed effect in the model, and re-running the single-variant analyses, no further significant associations were found for any of the 6 FA or C6:0-C14:0. These findings suggest that one polymorphism underlying this QTL on BTA17 influences multiple de novo synthesized milk FA. Thirteen genes in the QTL region were identified and analyzed carefully. Six out of the 8 SNP that showed the strongest associations were located in the La ribonucleoprotein domain family, member 1B (LARP1B) gene, and we suggest LARP1B as a primary candidate gene. Another gene of interest for this QTL region might be PKL4. None of these suggested candidate genes have previously been associated with milk fat synthesis or milk FA composition.

Duchemin SI ，Bovenhuis H ，Megens HJ ，Van Arendonk JAM ，Visker MHPW ... -  《-》

Effects of candidate gene polymorphisms on the detailed fatty acids profile determined by gas chromatography in bovine milk.

Association analyses between candidate genes and bovine milk fatty acids can improve our understanding of genetic variation in milk fatty acid profiles and reveal potential opportunities to tailor milk fat composition through selection strategies. In this work, we investigated the association of 51 single nucleotide polymorphisms (SNP) selected from 37 candidate genes using a functional and positional approach, with 47 fatty acids, 9 fatty acid groups, and 5 Δ(9)-desaturation indices in milk samples from Brown Swiss cows. Individual milk samples were collected from 1,158 Italian Brown Swiss cows, and gas chromatography was used to obtain detailed milk fatty acid compositions. A GoldenGate assay system (Illumina, San Diego, CA) was used to perform genotype 96 selected SNP located in 54 genes across 22 chromosomes. In total, 51 polymorphic SNP in 37 candidate genes were retained for the association analysis. A Bayesian linear animal model was used to estimate the contribution of each SNP. A total of 129 tests indicated relevant additive effects between a given SNP and a single fatty acid trait; 38 SNP belonging to 30 genes were relevant for a total of 57 fatty acid traits. Most of the studied fatty acid traits (~81%) were relevantly associated with multiple SNP. Relevantly associated SNP were mainly found in genes related to fat metabolism, linked to or contained in previously identified quantitative trait loci for fat yield or content, or associated with genes previously identified in association analyses with milk fatty acid profiles in other cow breeds. The most representative candidate genes were LEP, PRL, STAT5A, CCL3, ACACA, GHR, ADRB2, LPIN1, STAT1, FABP4, and CSN2. In particular, relevant associations with SNP located on bovine chromosome 19 (BTA19) were found. Two candidate genes on BTA19 (CCL3 and ACACA) were relevantly associated with de novo short- and medium-chain fatty acids, likely explaining the high heritability values found for these fatty acids (with the exception of C6:0). Two additional genes on BTA19 (CCL2 and GH1) showed associations with saturated and branched-chain fatty acids. Our findings provide basic information on genes and SNP affecting the milk fatty acid composition of dairy cows. These results may support the possibility of using genetic selection to modify milk fatty acid profiles to promote beneficial health-related effects.

Pegolo S ，Cecchinato A ，Mele M ，Conte G ，Schiavon S ，Bittante G ... -  《-》

SNP co-association and network analyses identify E2F3, KDM5A and BACH2 as key regulators of the bovine milk fatty acid profile.

The fatty acid (FA) profile has a considerable impact on the nutritional and technological quality of milk and dairy products. The molecular mechanism underlying the regulation of fat metabolism in bovine mammary gland have been not completely elucidated. We conducted genome-wide association studies (GWAS) across 65 milk FAs and fat percentage in 1,152 Brown Swiss cows. In total, we identified 175 significant single nucleotide polymorphism (SNPs) spanning all chromosomes. Pathway analyses revealed that 12:0 was associated with the greatest number of overrepresented categories/pathways (e.g. mitogen-activated protein kinase (MAPK) activity and protein phosphorylation), suggesting that it might play an important biological role in controlling milk fat composition. An Associated Weight Matrix approach based on SNP co-associations predicted a network of 791 genes related to the milk FA profile, which were involved in several connected molecular pathways (e.g., MAPK, lipid metabolism and hormone signalling) and undetectable through standard GWAS. Analysis of transcription factors and their putative target genes within the network identified BACH2, E2F3 and KDM5A as key regulators of milk FA metabolism. These findings contribute to increasing knowledge of FA metabolism and mammary gland functionality in dairy cows and may be useful in developing targeted breeding practices to improve milk quality.

Pegolo S ，Dadousis C ，Mach N ，Ramayo-Caldas Y ，Mele M ，Conte G ，Schiavon S ，Bittante G ，Cecchinato A ... -  《Scientific Reports》