Genome-wide association mapping for milk fat composition and fine mapping of a QTL for de novo synthesis of milk fatty acids on bovine chromosome 13

AMP-activated protein kinase mediates (-)-epigallocatechin-3-gallate (EGCG) to promote lipid synthesis in mastitis cows

Mastitis, a serious threat to the health and milk production function of dairy cows decreases milk quality. Blood from three healthy cows and three mastitis cows were collected in this study and their transcriptome was sequenced using the Illumina HiSeq platform. Differentially expressed genes (DEGs) were screened according to the |log2FoldChange| > 1 and P-value < 0.05 criteria. Pathway enrichment and functional annotation were performed through KEGG and GO analyses. Finally, the mechanism of the AMP-activated protein kinase (AMPK) mediation of (-)-epigallocatechin-3-gallate (EGCG) to promote lipid metabolism in mastitis cows was analyzed in bovine mammary epithelial cells (BMECs). Transcriptome analysis revealed a total of 825 DEGs, with 474 genes showing increased expression and 351 genes showing decreased expression. The KEGG analysis of DEGs revealed that they were mainly linked to tumour necrosis factor, nuclear factor-κB signalling pathway, and lipid metabolism-related signalling pathway, whereas GO functional annotation found that DEGs were enriched in threonine and methionine kinase activity, cellular metabolic processes, and cytoplasm. AMPK expression, which is involved in several lipid metabolism pathways, was downregulated in mastitis cows. The results of in vitro experiments showed that the inhibition of AMPK promoted the expression of lipid synthesis genes in lipopolysaccharide-induced BMECs and that EGCG could promote lipid synthesis by decreasing the expression of AMPK and downregulating the expression of inflammatory factors in inflammatory BMECs. In conclusion, our study demonstrated that AMPK mediated EGCG to inhabit of inflammatory responses and promote of lipid synthesis in inflammatory BMECs.

Genome-Wide Identification of Specific Genetic Loci Common to Sheep and Goat

Sheep and goat may become carriers of some zoonotic diseases. They are important livestock and experimental model animals for human beings. The fast and accurate identification of genetic materials originating from sheep and goat can prevent and inhibit the spread of some zoonotic diseases, monitor market product quality, and maintain the stability of animal husbandry and food industries. This study proposed a methodology for identifying sheep and goat common specific sites from a genome-wide perspective. A total of 150 specific sites were selected from three data sources, including the coding sequences of single copy genes from nine species (sheep, goat, cow, pig, dog, horse, human, mouse, and chicken), the dbSNPs for these species, and human 100-way alignment data. These 150 sites exhibited low intraspecific heterogeneity in the resequencing data of 1450 samples from five species (sheep, goat, cow, pig, and chicken) and high interspecific divergence in the human 100-way alignment data after quality control. The results were proven to be reliable at the data level. Using the process proposed in this study, specific sites of other species can be screened, and genome-level species identification can be performed using the screened sites.

Review: genetic background of milk fatty acid synthesis in bovines

Milk fat composition is an important trait for the dairy industry as it directly influences the nutritional and technological properties of milk and other dairy products. The synthesis of milk fat is a complex process regulated by a network of genes. Thus, understanding the genetic variation and molecular mechanisms regulating milk fat synthesis will help to improve the nutritional quality of dairy products. In this review, we provide an overview of milk fat synthesis in bovines along with the candidate genes involved in the pathway. We also discuss de novo synthesis of fatty acids (ACSS, ACACA, FASN), uptake of FAs (FATP, FAT, LPL), intracellular activation and channelling of FAs (ACSL, FABP), elongation (EVOLV6), desaturation (SCD, FADS), formation of triglycerides (GPAM, AGPAT, LIPIN, DGAT), and milk lipid secretion (BTN1A1, XDH, PLIN2). The genetic variability of individual fatty acids will help to develop selection strategies for obtaining a healthier milk fat profile in bovines. Thus, this review will offer a potential understanding of the molecular mechanisms that regulate milk fat synthesis in bovines.

Single-step genome-wide association for selected milk fatty acids in Dual-Purpose Belgian Blue cows

The aim of this study was to estimate genetic parameters and identify genomic regions associated with selected individual and groups of milk fatty acids (FA) predicted by milk mid-infrared spectrometry in Dual-Purpose Belgian Blue cows. The used data were 69,349 test-day records of milk yield, fat percentage, and protein percentage along with selected individual and groups FA of milk (g/dL milk) collected from 2007 to 2020 on 7,392 first-parity (40,903 test-day records), and 5,185 second-parity (28,446 test-day records) cows distributed in 104 herds in the Walloon Region of Belgium. Data of 28,466 SNPs, located on 29 Bos taurus autosomes (BTA), of 1,699 animals (639 males and 1,060 females) were used. Random regression test-day models were used to estimate genetic parameters through the Bayesian Gibbs sampling method. The SNP solutions were estimated using a single-step genomic best linear unbiased prediction approach. The proportion of genetic variance explained by each 25-SNP sliding window (with an average size of ~2 Mb) was calculated, and regions accounting for at least 1.0% of the total additive genetic variance were used to search for candidate genes. Average daily heritability estimated for the included milk FA traits ranged from 0.01 (C4:0) to 0.48 (C12:0) and 0.01 (C4:0) to 0.42 (C12:0) in the first and second parities, respectively. Genetic correlations found between milk yield and the studied individual milk FA, except for C18:0, C18:1 trans, C18:1 cis-9, were positive. The results showed that fat percentage and protein percentage were positively genetically correlated with all studied individual milk FA. Genome-wide association analyses identified 11 genomic regions distributed over 8 chromosomes [BTA1, BTA4, BTA10, BTA14 (4 regions), BTA19, BTA22, BTA24, and BTA26] associated with the studied FA traits, though those found on BTA14 partly overlapped. The genomic regions identified differed between parities and lactation stages. Although these differences in genomic regions detected may be due to the power of quantitative trait locus detection, it also suggests that candidate genes underlie the phenotypic expression of the studied traits may vary between parities and lactation stages. These findings increase our understanding about the genetic background of milk FA and can be used for the future implementation of genomic evaluation to improve milk FA profile in Dual-Purpose Belgian Blue cows.

Genetic variants associated with two major bovine milk fatty acids offer opportunities to breed for altered milk fat composition

Background

Although bovine milk is regarded as healthy and nutritious, its high content of saturated fatty acids (FA) may be harmful to cardiovascular health. Palmitic acid (C16:0) is the predominant saturated FA in milk with adverse health effects that could be countered by substituting it with higher levels of unsaturated FA, such as oleic acid (C18:1cis-9). In this work, we performed genome-wide association analyses for milk fatty acids predicted from FTIR spectroscopy data using 1811 Norwegian Red cattle genotyped and imputed to a high-density 777k single nucleotide polymorphism (SNP)-array. In a follow-up analysis, we used imputed whole-genome sequence data to detect genetic variants that are involved in FTIR-predicted levels of C16:0 and C18:1cis-9 and explore the transcript profile and protein level of candidate genes.

Results

Genome-wise significant associations were detected for C16:0 on Bos taurus (BTA) autosomes 11, 16 and 27, and for C18:1cis-9 on BTA5, 13 and 19. Closer examination of a significant locus on BTA11 identified the PAEP gene, which encodes the milk protein β-lactoglobulin, as a particularly attractive positional candidate gene. At this locus, we discovered a tightly linked cluster of genetic variants in coding and regulatory sequences that have opposing effects on the levels of C16:0 and C18:1cis-9. The favourable haplotype, linked to reduced levels of C16:0 and increased levels of C18:1cis-9 was also associated with a marked reduction in PAEP expression and β-lactoglobulin protein levels. β-lactoglobulin is the most abundant whey protein in milk and lower levels are associated with important dairy production parameters such as improved cheese yield.

Conclusions

The genetic variants detected in this study may be used in breeding to produce milk with an improved FA health-profile and enhanced cheese-making properties.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00731-9.

AGPAT3 Gene polymorphisms are associated with milk production traits in Chinese Holstein cows

The current study reports the identification of previously undiscovered single-nucleotide polymorphisms (SNPs) in the bovine AGPAT3 gene and further investigates their associations with milk production traits. Our results demonstrate that the major allele C of the SNP g.12264 C > T is positively correlated with test-day milk yield, protein percentage and 305-day milk yield. Importantly, in silico analysis showed that the C/T transition at this locus gives rise to two new transcription factor binding sites (TFBS), E2F1 and Nkx3-2. Polymorphism g.18658 G > A was the only SNP associated with milk urea nitrogen (MUN) with the G allele related to an increase in milk urea nitrogen as well as fat percentage. The GG genotype of SNP g.28731 A > G was associated with the highest fat and protein percentage and lowest 305-day milk yield and somatic cell score (SCS). The association between AGPAT3 locus and milk production traits could be utilized in marker-assisted selection for the genetic improvement of milk production traits and, probably in conjunction with other traits, for selection to improve fitness of dairy cattle.

SNPs with intermediate minor allele frequencies facilitate accurate breed assignment of Indian Tharparkar cattle

Tharparkar cattle breed is widely known for its superior milch quality and hardiness attributes. This study aimed to develop an ultra-low density breed-specific single nucleotide polymorphism (SNP) genotype panel to accurately quantify Tharparkar populations in biological samples. In this study, we selected and genotyped 72 Tharparkar animals randomly from Cattle & Buffalo Farm of IVRI, India. This Bovine SNP50 BeadChip genotypic datum was merged with the online data from six indigenous cattle breeds and five taurine breeds. Here, we used a combination of pre-selection statistics and the MAF-LD method developed in our laboratory to analyze the genotypic data obtained from 317 individuals of 12 distinct breeds to identify breed-informative SNPs for the selection of Tharparkar cattle. This methodology identified 63 unique Tharparkar-specific SNPs near intermediate gene frequencies. We report several informative SNPs in genes/QTL regions affecting phenotypes or production traits that might differentiate the Tharparkar breed.

Melatonin suppresses milk fat synthesis by inhibiting the mTOR signaling pathway via the MT1 receptor in bovine mammary epithelial cells

Milk fat content is an important criterion for assessing milk quality and is one of the main target traits of dairy cattle breeding. Recent studies have shown the importance of melatonin in regulating lipid metabolism, but the potential effects of melatonin on milk fat synthesis in bovine mammary epithelial cells (BMECs) remain unclear. Here, we showed that melatonin supplementation at 10 μmol/L significantly downregulated the mRNA expression of lipid metabolism-related genes and resulted in lower lipid droplet formation and triglyceride accumulation. Moreover, melatonin significantly upregulated melatonin receptor subtype melatonin receptor 1a (MT1) gene expression, and the negative effects of melatonin on milk fat synthesis were reversed by treatment with the nonselective MT1/melatonin receptor subtype melatonin receptor 1b (MT2) antagonist. However, a selective MT2 antagonist did not modify the negative effects of melatonin on milk fat synthesis. In addition, KEGG analysis revealed that melatonin inhibition of milk fat synthesis may occur via the mTOR signaling pathway. Further analysis revealed that melatonin significantly suppressed the activation of the mTOR pathway by restricting the phosphorylation of mTOR, 4E-BP1, and p70S6K, and the inhibition of melatonin on milk fat synthesis was reversed by mTOR activator MHY1485 in BMECs. Furthermore, in vivo experiments in Holstein dairy cows showed that exogenous melatonin significantly decreased milk fat concentration. Our data from in vitro and in vivo studies revealed that melatonin suppresses milk fat synthesis by inhibiting the mTOR signaling pathway via the MT1 receptor in BMECs. These findings lay a foundation to identify a new potential means for melatonin to modulate the fat content of raw milk in Holstein dairy cows.

A novel mutation in the promoter region of RPL8 regulates milk fat traits in dairy cattle by binding transcription factor Pax6

Ribosomal protein L8 (RPL8) was considered as a promising candidate gene for the milk fat percentage trait in dairy cattle in our previous genome-wide association studies, but the mechanism remains to be determined. Here we investigated the molecular mechanism underlying the effect of bovine RPL8 on milk fat percentage. We demonstrated that RPL8 silencing in bovine mammary epithelial cells affected the expression of genes encoding fat-related enzymes (ACACA, FASN, ACSS1, FABP3, SREBP-1, DGAT1, GPAM, PLIN2, PLIN5 and CIDEA). Furthermore, we showed here that a single nucleotide polymorphism, g.-931G > T (chr14:1508300, UMD3.1) in the putative RPL8 promoter region significantly reduced its promoter activity. Interestingly, this decrease in activity was paralleled by lower RPL8 expression in mammary gland tissues of dairy cattle with the homozygous TT genotype compared to that of cattle with the wild-type homozygous GG genotype. Importantly, we found g.-931G > T added a paired box 6 (Pax6)-binding site and this mutation located in the presumed Pax6-binding site. EMSA and co-immunoprecipitation (Co-IP) assays confirmed the interaction between RPL8 and Pax6 and the T allele exhibited a higher affinity of DNA/protein interactions than G allele, suggesting that Pax6 is an important transcription factor for RPL8 expression. In addition, lactating cows with the GG and GT genotypes presented a significant decrease in milk fat percentage compared to cows with TT genotypes. Altogether, our study indicated that g.-931G > T at RPL8 promoter altered its expression by affecting the interplay between Pax6 and RPL8, which may account for the association with milk fat traits. Findings herein first elucidated the biological function of RPL8 gene in milk fat and the identified SNP g.-931G > T may be considered as genetic makers for breeding in dairy cattle.

A post-GWAS confirming effects of PRKG1 gene on milk fatty acids in a Chinese Holstein dairy population

BACKGROUND

We previously conducted a genome-wide association study (GWAS) strategy for milk fatty acids in Chinese Holstein, and identified 83 genome-wide significant single nucleotide polymorphisms (SNPs) and 314 suggestive significant SNPs. Among them, two SNPs, BTB-01077939 and BTA-11275-no-rs associated with C10:0, C12:0, and C14 index (P = 0.000014 ~ 0.000024), were within and close to (0.85 Mb) protein kinase, cGMP-dependent, type І (PRKG1) gene on BTA26, respectively. PRKG1 gene plays a key role in lipolysis to release fatty acids and glycerol through the hydrolysis of triacyglycerol in adipocytes. We herein considered it as a promising candidate for milk fatty acids. The purpose of this study was to investigate whether PRKG1 had effects on milk fatty acids.

RESULTS

By direct sequencing the PCR products of pooled DNA, we identified a total of six SNPs, including one in 5' flanking region, four in 3' untranslated region (UTR), and one in 3' flanking region. The single-locus association analysis was carried out, and showed that the six SNPs mainly had significant associations with C6:0, C8:0 and C17:1 (P < 0.0001 ~ 0.0035). In addition, we observed a haplotype block formed by g.6903810G > A and g.6904047G > T with Haploview 4.1, and it was strongly associated with C8:0, C10:0, C16:1, C17:1, C20:0 and C16 index (P = < 0.0001 ~ 0.0123). The SNP, g.8344262A > T, was predicted to alter the binding site (BS) of transcription factor (TF) GAGA box with Genomatix software, and the subsequent luciferase assay verified that it really changed the transcriptional activity of PRKG1 gene (P = 0.0009).

CONCLUSION

In conclusion, to our best of knowledge, we are the first who identified the significant effects of PRKG1 on milk fatty acids in dairy cattle.

Genome-wide association study for milk production traits in a Brazilian Holstein population

Advances in the molecular area of selection have expanded knowledge of the genetic architecture of complex traits through genome-wide association studies (GWAS). Several GWAS have been performed so far, but confirming these results is not always possible due to several factors, including environmental conditions. Thus, our objective was to identify genomic regions associated with traditional milk production traits, including milk yield, somatic cell score, fat, protein and lactose percentages, and fatty acid composition in a Holstein cattle population producing under tropical conditions. For this, 75,228 phenotypic records from 5,981 cows and genotypic data of 56,256 SNP from 1,067 cows were used in a weighted single-step GWAS. A total of 46 windows of 10 SNP explaining more than 1% of the genetic variance across 10 Bos taurus autosomes (BTA) harbored well-known and novel genes. The MGST1 (BTA5), ABCG2 (BTA6), DGAT1 (BTA14), and PAEP (BTA11) genes were confirmed within some of the regions identified in our study. Potential novel genes involved in tissue damage and repair of the mammary gland (COL18A1), immune response (LTTC19), glucose homeostasis (SLC37A1), synthesis of unsaturated fatty acids (LTBP1), and sugar transport (SLC37A1 and MFSD4A) were found for milk yield, somatic cell score, fat percentage, and fatty acid composition. Our findings may assist genomic selection by using these regions to design a customized SNP array to improve milk production traits on farms with similar environmental conditions.

Multi-population GWAS and enrichment analyses reveal novel genomic regions and promising candidate genes underlying bovine milk fatty acid composition

Background

The power of genome-wide association studies (GWAS) is often limited by the sample size available for the analysis. Milk fatty acid (FA) traits are scarcely recorded due to expensive and time-consuming analytical techniques. Combining multi-population datasets can enhance the power of GWAS enabling detection of genomic region explaining medium to low proportions of the genetic variation. GWAS often detect broader genomic regions containing several positional candidate genes making it difficult to untangle the causative candidates. Post-GWAS analyses with data on pathways, ontology and tissue-specific gene expression status might allow prioritization among positional candidate genes.

Results

Multi-population GWAS for 16 FA traits quantified using gas chromatography (GC) in sample populations of the Chinese, Danish and Dutch Holstein with high-density (HD) genotypes detects 56 genomic regions significantly associated to at least one of the studied FAs; some of which have not been previously reported. Pathways and gene ontology (GO) analyses suggest promising candidate genes on the novel regions including OSBPL6 and AGPS on Bos taurus autosome (BTA) 2, PRLH on BTA 3, SLC51B on BTA 10, ABCG5/8 on BTA 11 and ALG5 on BTA 12. Novel genes in previously known regions, such as FABP4 on BTA 14, APOA1/5/7 on BTA 15 and MGST2 on BTA 17, are also linked to important FA metabolic processes.

Conclusion

Integration of multi-population GWAS and enrichment analyses enabled detection of several novel genomic regions, explaining relatively smaller fractions of the genetic variation, and revealed highly likely candidate genes underlying the effects. Detection of such regions and candidate genes will be crucial in understanding the complex genetic control of FA metabolism. The findings can also be used to augment genomic prediction models with regions collectively capturing most of the genetic variation in the milk FA traits.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5573-9) contains supplementary material, which is available to authorized users.

Integrated PTR-ToF-MS, GWAS and biological pathway analyses reveal the contribution of cow's genome to cheese volatilome

Volatile organic compounds (VOCs) are small molecules that contribute to the distinctive flavour of cheese which is an important attribute for consumer acceptability. To investigate whether cow's genetic background might contribute to cheese volatilome, we carried out genome-wide association studies (GWAS) and pathway-based analyses for 173 spectrometric peaks tentatively associated with several VOCs obtained from proton-transfer-reaction mass spectrometry (PTR-ToF-MS) analyses of 1,075 model cheeses produced using raw whole-milk from Brown Swiss cows. Overall, we detected 186 SNPs associated with 120 traits, several of which mapped close to genes involved in protein (e.g. CSN3, GNRHR and FAM169A), fat (e.g. AGPAT3, SCD5, and GPAM) and carbohydrate (e.g. B3GNT2, B4GALT1, and PHKB) metabolism. Gene set enrichment analysis showed that pathways connected with proteolysis/amino acid metabolism (purine and nitrogen metabolism) as well as fat metabolism (long-term potentiation) and mammary gland function (tight junction) were overrepresented. Our results provide the first evidence of a putative link between cow's genes and cheese flavour and offer new insights into the role of potential candidate loci and the biological functions contributing to the cheese volatilome.

Unravelling genetic variation underlying de novo-synthesis of bovine milk fatty acids

The relative abundance of specific fatty acids in milk can be important for consumer health and manufacturing properties of dairy products. Understanding of genes controlling milk fat synthesis may contribute to the development of dairy products with high quality and nutritional value. This study aims to identify key genes and genetic variants affecting de novo synthesis of the short- and medium-chained fatty acids C4:0 to C14:0. A genome-wide association study using 609,361 SNP markers and 1,811 animals was performed to detect genomic regions affecting fatty acid levels. These regions were further refined using sequencing data to impute millions of additional genetic variants. Results suggest associations of PAEP with the content of C4:0, AACS with the content of fatty acids C4:0-C6:0, NCOA6 or ACSS2 with the longer chain fatty acids C6:0-C14:0, and FASN mainly associated with content of C14:0. None of the top-ranking markers caused amino acid shifts but were mostly situated in putatively regulating regions and suggested a regulatory role of the QTLs. Sequencing mRNA from bovine milk confirmed the expression of all candidate genes which, combined with knowledge of their roles in fat biosynthesis, supports their potential role in de novo synthesis of bovine milk fatty acids.