Next Article in Journal
Bioeconomic Modelling to Assess the Impacts of Using Native Shrubs on the Marginal Portions of the Sheep and Beef Hill Country Farms in New Zealand
Previous Article in Journal
Sources of Resistance to Powdery Mildew in Barley Landraces from Turkey
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Genetic Markers Associated with Milk Production Traits in Dairy Cattle

1
State Key Laboratory of Animal Nutrition, Beijing Engineering Technology Research Center of Raw Milk Quality and Safety Control, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China
2
Faculty of Veterinary and Animal Sciences, The University of Agriculture, Dera Ismail Khan 29220, Pakistan
*
Author to whom correspondence should be addressed.
These two authors equally contributed to this work.
Submission received: 24 August 2021 / Revised: 28 September 2021 / Accepted: 14 October 2021 / Published: 18 October 2021
(This article belongs to the Special Issue Livestock Breeding and Conservation Genetics)

Abstract

:
Increasing milk production is one of the key concerns in animal production. Traditional breeding has gotten limited achievement in the improvement of milk production because of its moderate heritability. Milk production traits are controlled by many genes. Thus, identifying candidate genes associated with milk production traits may provide information that can be used to enhance the accuracy of animal selection for moderately heritable traits like milk production. The genomic selection can enhance the accuracy and intensity of selection and shortening the generation interval. The genetic progress of economically important traits can be doubled with the accuracy of selection and shortening of generation interval. Genome-wide association studies (GWAS) have made possible the screening of several single nucleotide polymorphisms (SNPs) in genes associated with milk production traits in dairy cattle. In addition, RNA-sequencing is another well-established tool used to identify genes associated with milk production in dairy cattle. Although it has been widely accepted that these three methods (GWAS, RNA-seq and DNA sequencing) are considered the first step in the screening of genes, however, the outcomes from GWAS, DNA-sequencing and RNA-seq still need further verification for the establishment of bonafide causal variants via genetic replication as well as functional validation. In the current review, we have highlighted genetic markers identified (2010-to date) for their associations with milk production traits in dairy cattle. The information regarding candidate genes associated with milk production traits provided in the current review could be helpful to select the potential genetic markers for the genetic improvement of milk production traits in dairy cattle.

1. Introduction

Milk production traits have fundamental roles in dairy development and related economy [1,2]. The bovine milk production traits such as milk yield, fat content, protein content and somatic cell score (SCS) are the essential economic traits used to measure the quality of milk [3,4]. Traditional breeding methods have achieved considerable success in many economic traits; however, milk production having moderate heritability, the gains was not fruitful with common traditional breeding [5]. Being a polygenetic trait, milk production is controlled by many genes [6,7]. Thus exploring the genetic changes underlying preferred phenotypes is the target of today’s animal producers. It has been well-established that the production of milk can be enhanced through genetic marker-assisted selection [8,9]. Various approaches such as mapping of quantitative trait loci (QTL), genome-wide association study (GWAS), RNA-sequencing, whole-genome sequencing and candidate gene analysis have been used to screen out the causal genes or their mutations associated with milk production traits [10,11,12,13,14]. So far, many candidate genes or polymorphisms within these genes have been identified that have a positive correlation with milk production traits in dairy cattle [1,15,16].
Through genomic selection, we can identify genetically superior animals at a very early age. The DNA-tested animals can get accurate genomically enhanced breeding values before they enter into sexual maturity. In addition, because of the heavier use of young, genetically superior males and females in genomic selection, the generation interval can be decreased. The intensity of selection can be enhanced because the breeders use genomic testing to identify a larger group of potentially superior animals. Altogether by enhancing the accuracy and intensity of selection and decreasing the generation interval, the rate of genetic progress for economically essential dairy traits can be almost doubled. Keeping in view the importance of genomic selection, the current review was designed to highlight the possible development on genetic markers associated with milk production traits in dairy cattle.

2. Materials and Methods

The data were collected through authentic sources, such as PubMed, ScienceDirect, Web of Science, SpringerLink, Scopus, and Google Scholar, using polymorphism, genetic markers, GWAS, RNA-seq, DNA-sequencing (Whole genome sequencing) and dairy cattle milk production traits as major keywords. All the published studies that have discussed the polymorphisms in genes and their association with milk production traits in dairy cattle were included in the current review. Moreover, we also included the published studies that reported the direct effect of genes on milk production traits in dairy cattle. Similarly, all the published articles in the English language and scientific citation index (SCI) peer-reviewed journals were incorporated for discussion in the current review. Furthermore, we considered articles (approximately 96%) published from 2010 onward in a present review article. Those genes from RNA-seq data associated with milk production traits and differentially expressed (p < 0.05, Q < 0.05) or validated through qPCR, were selected in the current review article. The present review article included all the polymorphisms in genes reported through GWAS or functional validation that were significantly associated with milk production traits in dairy cattle. Conversely, we excluded the data that was available in the form of conference papers, books, book chapters, thesis data and unpublished findings.

3. Genome-Wide Association Studies (GWAS) for Screening of Genetic Markers for Milk Production in Dairy Cattle

Genome-wide association studies have been extensively practiced to screen the polymorphism in genes associated with milk production traits in dairy cattle [17,18,19,20]. The detail of genes and their SNPs has been summarized in Table 1. Recently, Jiang et al. reported eight genes (ACSBG2, NLK, UGDH, MAP3K1, TBC1D1, RETSAT, CENPE, and FCGR2B) associated with milk production traits in 769 Chinese Holstein cows through GWAS study (Table 1) [21]. Consistently, another study documented 22 genes (SLC37A1, ALPL, MGST1, ABCG2, MEPE, PKD2, HERC3, SEPSECS, SEL1L3, DHX15, CSN1S1, CSN1S2, CSN2, CSN3, PAEP, DGAT1, RECQL4, MROH1, BOP1, ANKH, AGPAT6 and PICALM) through GWAS study which were linked to milk protein in 8080 cows (2967 Montbéliarde, 2737 Normande, and 2306 Holstein) [22]. Iung et al. performed the weighted single-step genomic BLUP (WssGBLUP) method to identify genomic regions correlated with milk production traits in the Brazilian Holstein cattle population. For this purpose 75,228 phenotypic records from 5981 cows were obtained, while genotypic data of 56,256 SNPs from 1067 cows were selected for the GWAS study [23]. Finally, ABCG2, DGAT1, MGST1, SLC37A1, LTBP1, LRRC19, PDE9A, and PAEP genes were found to be link with milk production phenotypic traits in the Brazilian Holstein cattle population [23]. By using the GWAS study, the genomic regions on BTA14 were explored for their association with milk production traits in Holstein’s population (21,068 cows) selected from four different countries (Belgium, The Netherlands, Great Britain and Denmark) [24]. Through weighted single-step genomic BLUP approach, they reported several genes (MIR2308, LOC104973955, CYHR1, ZNF34, FOXH1, COMMD5, TONSL, PPP1R16A, CPSF1, RPL8, DGAT1, ARHGAP39, GPT, LRRC14, and GML) distributed on BTA14 that were significantly correlated with milk yield (Table 2). Many of the genes such as TONSL, PPP1R16A, FOXH1, ARHGAP39, CYHR1 and ARHGAP39 have also been documented by previous studies for their association with milk production traits [25,26,27]. Furthermore, Nayeri et al. documented that SNPs in CPSF1, DGAT1, TONSL, CYHR1, FOXH1 and PPP1R16A were associated with milk yield in Canadian Holsteins [12,25]. Buitenhuis et al. [25] also showed that GML is significantly linked to milk fat and protein in dairy cattle. Interestingly, Poulsen et al. documented several genes such as ALG3, B3GALNT2, LOC520336, PIGV, MAN1C1, ST6GALNAC6, GLT6D1, GALNT14, GALNT17, COLGALT2, LFNG and SIGLEC by performing GWAS analysis in Danish Jersey and Holsteins [28].
By using the GWAS study, Ariyarathne et al. [29] reported key genetic markers that were associated with FP (MGST1, DGAT1, CEBPD, SLC52A2, GPAT4, and ACOX3), PP (CSN1S1, GOSR2, HERC6, and IGF1R) and milk urea (GMDS, E2F7, SIAH1, SLC24A4, LGMN, and ASS1) of Holstein Friesian, Jersey or crossbred cows in New Zealand [29]. Similarly, Bouwman et al. documented candidate genes (ABCG2, DGAT1 SCD1, ACLY, SREBF1, STAT5A, GH, PPARGC1A, ACSS2, AGPAT6 and FASN) that were significantly correlated with milk fatty acid traits in lactating Dutch Holstein Friesian at Netherlands [17]. Li et al. conducted a GWAS study for milk fatty acid traits in 784 Chinese Holstein cows and found that polymorphisms in some key genes showed a link with milk fatty acid traits [30]. Although we highlighted several genes documented through the GWAS study, however, the functional validation of these genes is highly warranted before adding them as genetic markers for milk improvement in dairy cattle breeding.
Table 1. GWAS study for screening genetic markers associated with milk production trait.
Table 1. GWAS study for screening genetic markers associated with milk production trait.
SNP (Gene)Production TraitsBreed and Phenotypic Traits and Method Used for AssociationCountryAuthor
rs381714237 (FCGR2B)MY, PY and PPChinese HolsteinChina[21]
ss2137349053 (CENPE)
rs385060942
ss2137349051
rs453960300
rs378415122
MY, FY and PYChinese HolsteinChina[21]
rs134985825 (RETSAT)MY, FY, PP, PYChinese HolsteinChina[21]
rs377943075 (ACSBG2)Milk FY, PPChinese HolsteinChina[21]
rs136639319 (TBC1D1)Milk FP and PPChinese HolsteinChina[21]
rs379188781 (NLK)
rs134444531
Milk PY and PPChinese HolsteinChina[21]
ss2137349058 (MAP3K1)MY, FY and PYChinese HolsteinChina[21]
ss2019489562 (UGDH)MYChinese HolsteinChina[21]
BovineHD2400007916 (CDH2)Milk FPDual-purpose Xinjiang Brown cattle
2410 individuals with 6811 reproductive records and 5441 milk records
China[31]
BTB-01731924 (GABRG2)Milk PYDual-purpose Xinjiang Brown cattleChina[31]
rs136947640 (exon10) (FASN)Milk fat traitsDual-purpose Xinjiang Brown cattleChina[31]
rs41919985 (Exon-39) (A2266T)Milk fat traitsDual-purpose Xinjiang Brown cattleChina[31]
ARS-BFGL-NGS-24998 (SAA3)
UA-IFASA-8605 (SAA3)
BFGL-NGS-119420 (TRIB3)
ARS-BFGL-NGS-69013 (SESN2)
BTA-31250-no-rs (SESN2)
Hapmap53714-rs29017586 (CHAC1)
ARS-BFGL-NGS-5790 (NR4A1)
UA-IFASA-8605 (SAA1)
UA-IFASA-8605 (M-SAA3.2)
BTA-68781-no-rs (HIST1H2AC)
BTB-00411816 (THBS4)
ARS-BFGL-NGS-85980 (FAM71A)
ARS-BFGL-NGS-72191(H4)
ARS-BFGL-NGS-29557 (PTHLH)
ARS-BFGL-NGS-81082(ARID1B)
ARS-BFGL-NGS-72191 (BoLA-DQB)
ARS-BFGL-NGS-107749 (CDH16)
ARS-BFGL-NGS-29490 (VEGFA)
ARS-BFGL-NGS-85980 (ATF3)
BTB-01766447 (RPL23A)
Milk FP, PPU.S. Holstein cows (1654 cows, Thirty one dairy traits, including 13 production, health and reproduction traits and 18 body conformation traits) were selected for this studyU.S[14]
ARS-BFGL-NGS-24998 (SAA1)
ARS-BFGL-NGS-70836 (SAA1)
ARS-BFGL-NGS-100459 (RPL23A)
ARS-BFGL-NGS-24998 (M-SAA3.2)
Hapmap49309-BTA-78604 (P4HA2)
ARS-BFGL-NGS-70836 (ATF3)
Milk PPU.S. Holstein cows (1654 cows, Thirty one dairy traits, including 13 production, health and reproduction traits and 18 body conformation traits) were selected for this studyU.S[14]
ARS-BFGL-NGS-14781 (DDIT3)Milk FPU.S. Holstein cowsU.S[14]
rs41569048 (PTHLH)Milk PPDutch HolsteinNetherlands[18]
rs41590827 (PTHLH)Milk PPDutch Holstein, 1912 first-lactation Holstein-Friesian cows from 398 commercial herds, records of milk protein,
Significance threshold used for association
Netherlands[18]
rs41640170 (HEATR7B2)Milk PPDutch HolsteinNetherlands[18]
Hapmap51303-BTA-74377 (PTHLH)Milk FPChinese Holsteins
2093 daughters as well as their 14 corresponding sires were selected to perform the current study.
The numbers of daughters of the 14 sires range from 83 to 358 daughters with an average of 150.
Transmission-disequilibrium test (TDT)-based single locus regression analyses and mixed model-based single locus regression analyses were performed for association analysis
Milk production traits such as milk yield (MY), milk fat yield (FY), milk protein yield (PY), milk fat percentage (FP) and milk protein percentage (PP) were considered for this study
China[32]
ARS-BFGL-BAC-2469 (HEATR7B2)Milk FP, PPChinese HolsteinsChina[32]
BFGL-NGS-119420 (TRIB3)Milk FP, PPCanadian Holstein
Data from 462 Canadian Holstein bulls were collected
Single locus LD regression model was used to perform association analysis
Canada[33]
rs29016156 (VEGFA)Milk PPCanadian HolsteinCanada[33]
rs41640789 (VEGFA)Milk FPCanadian HolsteinCanada[33]
rs41590827 (PTHLH)Milk PPCanadian HolsteinCanada[33]
BTB-00213370 (NUB1)
ARS-BFGL-NGS-71395 (SLC24A2)
BTB-01052867 (SLC24A2)
Milk production traitsGir cross Holstein (Girolando)
Records of 305-day milk yield of 337 dairy cows
a single-marker linear regression model was used for association
Brazil[34]
BovineHD2900015534 (SLC22A8)
BovineHD1200012381 (KLHL1)
BovineHD1200012381 (KLHL1)
BTB-00074258 (TBC1D5)
Milk fatty acids
lactation persistence
445 Chinese Holsteins
15 milk production traits were used for this study
Fixed-effect linear regression model and a mixed-effect linear model were used for association
China[35]
BovineHD2500005573 (EEF2K)
BovineHD2500005573 (EEF2K)
MYChinese HolsteinsChina[35]
rs109421300 (DGAT1)
rs109528658 (EP400)
rs42295213 (EPHA6)
rs134480235 (SLCO1A2)
PP
FY
FP
PY
Chinese Holsteins
452,920 test-day records estimated breeding values from 61,600 cows
SNP regression was performed for association annalysis
China[36]
rs211223469 (DGAT1)FY, MYKorean Holstein
2780 Korean Holsteins (926 bulls and 1854 cows) were used in the current study
Single marker regression model for association analysis while MY, FY, PY, and SCS traits were used as milk production phenotypic traits
Korea[37]
rs41596885 (PDE4B),
rs42314807 (PDE4B)
FY, MYKorean HolsteinKorea[37]
rs43454033 (ANO2)FY, MYKorean HolsteinKorea[37]
ACACA rs110562092
ADRB2 rs132839139
AGPAT6 rs110445169
CARD15 rs43710288
CSN1S1 rs43703010
CSN2 rs43703011
CSN3 rs43703015
FABP4 rs110757796
FGF2 rs110937773
GHR rs109136815
LEP rs11055965
LEP rs29004170
LEPR rs43349286
LPIN1 rs136905033
LPIN1 rs137457402
ORL1 rs135588030
PPARGC1A rs44857081
PRL rs110684599
PRL rs211032652
SCD1 rs41255693
STAT1 rs43705173
STAT1 rs43706906
STAT5A rs109578101
STAT5A rs137182814
TLR2 rs43706433
XDH rs42890834
MFAs1158 Italian Brown Swiss cows and
The bayesian linear animal mode was used for association study were considered for current study
Italy[38]
MIR2308, LOC104973955, CYHR1, ZNF34, FOXH1, COMMD5,TONSL, PPP1R16A, MFSD3, LRRC24, RPL8, C14H8orf33, KIFC2, RECQL4, ZNF7, ARHGAP39, GPT, LRRC14, C14H8orf82,LOC100141215, MIR2309, MIR1839, LOC101907640,LOC101908059, GRINA, LOC104968841, LOC104973958,LOC104973959, LOC104973960, ARC, SPATC1, LOC786966, LOC104973961, OPLAH, HGH1, LOC509114, JRK, PARP10, MAF1, SHARPIN, CYC1, GPAA1, MROH1, LOC523023, EXOSC4, PSCA, LY6K, GML, LY6D, LOC100848939, LOC101904969, LOC101905222, LYPD2, LOC104973965, LYNX1, LOC104973966, THEM6, SLURP1, LOC78762Milk yield traitsIn current study the data were records of 21,068 lactations on primiparous (9910) and multiparous (11,158) Holstein cowsBelgium, The Netherlands, Great Britain and Denmark[24]
SEMA5B, AGPAT3,DGAT1, BTN1A1 SREBF1,FASN,GHR,PRLR, LIPJ, LIPK, ECHS1, ORBS1, NFKB2, CHUK,SCD1, AGPAT6Milk fatty acids784 Chinese and 371 Danish Holstein
and 16 milk FA traits were selected for current study
SNP linear regression models were used for association analysis
China [39]
SLC37A1, MGST1, ABCG2, CSN1S1, CSN2, CSN1S2, CSN3, PAEP, DGAT1, AGPAT6, ALPL, ANKH, PICALMMilk composition traitsMontbéliarde, Normande, Holstein
848,068 test-day milk samples from 156,660 cows in the first three lactations were collected
Genotyped in 080 cows (2967 MON, 2737 NOR, and 2306 HOL)
France[22]
SLC15A2, PEPT2, SND1, LEP, CLOCK, CASR, LRRC4, DOCK1, STAT1,STAT3, ELF5Milk fat traits1256 Holstein, 624 Gir, and 477 Girolando cattle
Genomic BLUP Model was used for association analysis
The database utilized in this study was consisted of 166,628 lactations from 94,124 cows, edited for age at calving (547 to 9095 d), calving season (dry or rainy), breed composition (determined by the proportion of Holstein breed, 1/4, 3/8, 1/2, 5/8, 3/4, or 7/8), and contemporary group (determined by herd-year of calving).
Brazil[40]
rs443751026-GPATCH4 upstream
rs210886822 MGST1 intron
rs209288972 MGST1 upstream
rs208317364 DGAT1 intron
rs133931291 HSF1 intron
LGALS12 upstream
rs383292923 ANK1 intron
rs208624037 GPAT4 intron
rs467849681 ARHGEF28 intron
Milk fat Montbéliarde, Normande, Holstein
15 traits (5 linked with milk production, 2 with udder health, and 8 with udder morphology) in Montbéliarde (MON), Normande (NOR), and Holstein (HOL) cattle were selected.
The phenotypic data was collected from 2515 MON, 2203 NOR, and 6321 HOL bulls and verified in 23,926 MON, 9400 NOR, and 51,977 HOL cows
Illumina Bovine SNP50 BeadChip (50K; Illumina Inc., San Diego, CA, USA) was used for genotyping
France[20]
BTB-01603522-ACSL1
BTB-01926888-ACSL1
BTA-111275-no-rs-PRKG1
BTB-01077939-PRKG1
Hapmap26394-BTA-136497-CNTN3
Hapmap26394-BTA-136497-CNTN3
ARS-BFGL-NGS-101978-HTR1B
BTA-12468-no-rs-IGF1R
BTA-12468-no-rs-IGF1R
BTA-76414-no-rs-IGF1R
BTA-76414-no-rs-IGF1R
ARS-BFGL-NGS-40159-PLIN1
Hapmap49848-BTA-106779-CPM
BTB-01556197-HTR1B
BTB-01556197-HTR1B
ARS-BFGL-BAC-35400-FAM46A
ARS-BFGL-BAC-35400-FAM46A
ARS-BFGL-NGS-61979-UBE3D
ARS-BFGL-NGS-61979-UBE3D
ARS-BFGL-NGS-34500-ACACA
ARS-BFGL-NGS-39328-FASN
Hapmap58547-rs29023020-PRL
ARS-BFGL-NGS-111111-AGPAT3
ARS-BFGL-NGS-109493-AGPAT3
BTA-56389-no-rs-AGPAT3
ARS-BFGL-NGS-45691-FABP3
ARS-BFGL-NGS-118924-FABP3
RS-BFGL-NGS-4939-DGAT1
ARS-BFGL-NGS-118998-GHR
BTB-01373917-STAT1
ARS-BFGL-NGS-33744-STAT1
BTB-00965197-NRG1
ARS-BFGL-NGS-107403-NFKB2
ARS-BFGL-NGS-23064-SCD1
ARS-BFGL-NGS-77668-SCD1
BTB-00930925-SCD1
ARS-BFGL-NGS-39397-SCD1
BTB-00930720-SCD1
Hapmap31825-BTA-158647-SCD1
Hapmap33073-BTA-162864-SCD1
BTB-00931481-SCD1
ARS-BFGL-NGS-110077-SCD1
ARS-BFGL-NGS-108305-SCD1
BTB-00931586-SCD1
ARS-BFGL-NGS-114149-SCD1
ARS-BFGL-NGS-116481-SCD1
Hapmap24832-BTA-138805-SCD1
ARS-BFGL-NGS-6259-SCD1
BTB-00932332-SCD1
Hapmap46411-BTA-15820-CHUK
BTA-61921-no-rs-LIPJ
ARS-BFGL-NGS-21794-LIPK
ARS-BFGL-NGS-29299-SORBS1
Hapmap41595-BTA-60800-SORBS1
Hapmap58930-rs29010490-SORBS1
Hapmap28763-BTA-162328-ECHS1
Hapmap28763-BTA-162328-ECHS1
ARS-BFGL-NGS-116897-OLR1
Hapmap26001-BTC-038813-PPARGC1A
Hapmap31284-BTC-039204-PPARGC1A
Hapmap49746-BTA-76106- PPARGC1A
ARS-BFGL-NGS-12970-FADS1
ARS-BFGL-NGS-1448-AGPAT6
ARS-BFGL-NGS-85864- CYP26A1
BTB-00927439- CYP2C19
BTB-01423653-PRLR
BTB-01423676-PRLR
Hapmap30570-BTA-152778-PRLR
Milk fatty acid traitsChinese Holstein
Phenotypic data for 22 milk fatty acids in 784 Chinese Holstein cows was used
Significance threshold was considered
China[30]
MY: milk yield; FY: fat yield; FP: fat percentage; PY: protein yield; PP: protein percentage; MFAs: Milk fatty acids. Eukaryotic elongation factor 2 kinase (EEF2K); kelch like family member 1 (KLHL1); EPH receptor A6 (EPHA6); solute carrier organic anion transporter family member 1A2 (SLCO1A2); diacylglycerol O-acyltransferase 1 (DGAT1); E1A binding protein p400 (EP400); phosphodiesterase 4B (PDE4B), and anoctamin 2 (ANO2).

4. Transcriptomic Analysis for Screening of Genetic Markers Associated with Milk Production

RNA-sequencing has been a newly merged tool for screening genetic markers associated with milk production [2,10]. Besides genetic data, the gene expression profile also plays a vital role in exploring the underlying mechanism for complex traits such as milk production in dairy cattle. Constantly, Cui et al. performed the transcriptomic profiling of the bovine mammary gland of four lactating Chinese Holstein cows with high and low phenotypic milk protein and fat percentage values. They reported some promising candidate genes (TRIB3, SAA1, M-SAA3.2, SAA3, VEGFA, PTHLH, HSPD1, KRT24, and RPL23A) that were significantly correlated with milk protein and milk fat percentage [10]. Bagnato et al. also reported the significant association of HSPD1 and KRT24 genes with milk yield and protein percentage in Brown Swiss dairy cattle [41]. Furthermore, Khan et al. by using RNA-seq analysis reported the number of genes (DGAT2, ALOX5, AGPAT4, GPAT3, GGH, ALDOA, TKT, SLC11A1 and LAP3) in response to folic acid treatment that were associated with milk protein, milk yield and milk fat in dairy cattle [2]. Consistently, Ouattara et al. reported that vitamin B9 and B12 combine supplementation regulated key genes that were associated with milk production traits. Furthermore, the candidate genes (MYOM1, HP, CDK5R1, MEP1B, DLK1, PPP1R3B, GSTA5, HERC6, LOXL4, SAA3, FUT5, PYCR1 and CACNA2D1) they documented were linked to milk protein and milk fat traits [42]. The genes associated with milk production traits screened out through RNA-seq by different studies have been summarized in Table 2.
Table 2. Genes identified through the RNA-seq method.
Table 2. Genes identified through the RNA-seq method.
GenesProduction TraitsBreed CountryAuthor
LRRC73, GPX3, APOA4, HP, MFSD2, CDC42EP5, SLC13A5, SMCT1, PAQR9, SFRP2, ISG15, IFIT1, RSAD2, APOA4, MX1, MX2, USP18, LOC100298356, HERC6, ISG12(B), TLH29, RSAD2, IFIT1, FKBP5, FKBP,RXRG, ITGAD, LYZ2,HBB APOC2, ACADVL, PPP1R3B, GALE, PKLR, ANGPTL4, CDKN1A, ODC1, LPIN1, DUSP1, LMNA, APOA1, ABCG8, Kb, SAA1, PC, SDS, GADD45B, IGF-1R, CYP7A1,GK, SGLT1, FBP2MY, FY, PY, FPHolsteinChina[43]
C4BPA, SLC25A38,BMX, EIF4G3, ZC3H14, FCAMR, DNER, SAA3, HEATR7B2, TRIB3, SESN2, CHAC1, NR4A1, SAA1, ATF3, RPL23A, CDH16, VEGFA, BoLA-DQB, ARID1B, PTHLH, H4, FAM71A, THBS4, DDIT3, M-SAA3.2, HIST1H2AC, P4HA2, HSPD1, KRT24, CDKN1AFP, PPHolsteinChina[10]
GGT5, CYP2J2, ALOX12, MIF, LPL, CPT1AMFAsHolsteinChina[44]
CSN2, CSN1S1, LGB, CSN3, CSN1S2, LALBA, GLYCAM1, COX1, FASN, CLU, COX3, MT-CYB, XDH, MFGE8, EEF1A1, GPAM, ATP6, MT-ND3, ND1, MT-ND4, NADH, SPP1, SERPINA1, CNTFR, ERBB2, NEDD4L, ANG, GALE, HSPA8, LPAR6, WAP, NARS, MARS, GARS, CDO1,GATM, INSR, IGF1R, IGFBP3, CRIM1, IGFBP3Milk protein traitsHolsteinChina[45]
SLC22A1, MAPK9, PPARGC1A, FOXO1, SOCS1, SOCS2, CREB1, HNF4A, HNF4G, GADD45A, DUSP1, PDGF, SYBU, DDIT4, BAMBI, MTHFR, SLC27A2, PCK1, CPT2, SIRT3, CYP4A11, PLCB2Milk protein and fatHolsteinChina[46]
SLC27A6, ACADM, ACADs, IDH1, FABP4, CACYBP, KLHL9, UBE2B, RPS, SLC7A8Milk proteinHolsteinChina[47]
LALBA, LGB, CSN1S1, CSN1S2, CSN2, CSN3, GK,GPD1, DHCR24, COQ2, AGPAT6, GPAM, LPIN1, BTN1A1, XDH, PLIN2, SCD5,DGAT1, FADS1, FABP3, SLC22A16, ACSS1, ACSS2, ADIPOQ, HADHB, SLC1A2, SLC1A5, SLC7A4, SLC7A8, SLC38A3, SARS, PAH, ASNS, RPL,ELOVL, EIF4EBP1, INSIG1Milk protein and fatHolsteinChina[47]
LAP3, ASS1, CYP2J2, ATP6AP1, SDS, DGAT2, AGPAT4, GPAT3, ALOX5, HSD17B12, HACD4, PPT2, ELOVL6, EPHX1, LPL, GUK1, XDHMilk fat and protein yieldHolsteinChina[2]
GPC5,TECTB,IARS2,GUK1,HS3ST5, STMN4,CALB1,LALBA,GLYCAM1 GP2, LPL, SLC34A2, TUBA1C, CSN1S1, CSN2, PTHLH, BDA20, BDA20, ALOX15, STATH, BOLA-DQA1, TTC36, PAEP, SPINK4, BTN1A1, TMOD4, SCD, MYBPC1, ASB11, SLC38A3MYHolsteinChina[48]
IRF6, AGPAT6, STAT5A, XDH, B4GALT1,BCL2L11, PRLR, ALOX5, PRKAA1, NCF1, AGPAT6,STAT5A, CRYL1, GPAM, ALOX12, PGHS-2, CPT1B, SLC27A1, PRKAR2B, FADS6, PPARD, ACACA, PTGES, EHHADHMilk fat synthesisHolsteinChina[49]
MY: milk yield; FY: fat yield; FP: fat percentage; PY: protein yield; PP: protein percentage; MFAs: Milk fatty acids.

5. Whole-Genome Sequencing for Screening of Genetic Marker Associated with Milk Production in Cattle

Whole-genome sequencing is one of the next-generation sequencing methods utilized to identify a large number of SNPs more quickly and inexpensively [50]. The DNA-seq method has been widely excised in livestock genomics to identify the genetic markers associated with milk production traits [51,52]. To perform whole-genome sequencing, the data composed of about 254 k milk, fat, and protein test-day records were collected from 7522 Holstein cows calved from 2006–2016 on two dairy farms in the state of Florida [53]. They reported several genes (CDKN1B, DUSP16, HSF1, EEF1D, VPS28, TONSL, PEX16, MAPK8IP1, CREB3L1 and CRY2) on BTA5, BTA14 and BTA15 that were significantly associated with milk production in traits in dairy cattle. Interestingly, the reported genes in this study are involved in the inositol phosphate mediated signaling pathway, insulin receptor signaling pathway, JNK cascade, stress-activated MAPK cascade, and glutamine metabolic process. These pathways play a major role in maintaining milk production even under stressful condition to regulate the antioxidant system [53]. Similarly, Nanaei et al. performed DNA sequencing for screening genetic makers associated with milk production traits by using the Illumina whole genomes of 21 cattle individuals, including 3 indigenous African breeds (Ankole n = 4, Kenana n = 4 and N’Daman = 6), and two commercial breeds (Polish Holstein-Friesian n = 3 and Hereford n = 4). Importantly, they documented some key genes (IGFBP2, B4GALT1, RORA, LPIN1, ATP2B, CSN3, NME1, ACACA, PDE3A, XP-CLR, KCNIP4, GHR, NF2, ABCC9, CD44, MACF1, IL15) involved in the regulation of biological function processes such as phosphorus metabolic process, phosphate-containing compound, metabolic process, phosphorylation, protein phosphorylation” and metal ion transport that were significantly related with milk production traits [54]. Recently a study selected 45 blood samples from two-year-old animals and DNA-sequencing was carried to identify genetic markers associated with milk production traits [55]. Interestingly they found nine genes (ADCY5, CACNA1A, CREB1, INHBA, INHBB, PIK3R1, PLCB1, PRKCE, and SMAD2) distributed in the ionotropic glutamate receptor pathway, the endothelin signaling pathway, and the gonadotropin-releasing hormone receptor pathway, which are involved in the hormonal regulation of lactation [55]. Whole-genome sequencing for data obtained from 4280 progeny tested Nordic Red Cattle bulls was performed to identify the genetic markers for milk production [56]. In addition, the genes related to milk production traits including DGAT1, HSF1, TRIM26, CLEC16A, NEURL1 (Fat yield), MKL1, CPSF1, ADCK5, LAX1, GHR (milk yield), DGAT1, HSF1, UNKL, PAM16, GLIS2 (protein yield) were documented in Nordic Red Cattle [56].

6. DNA Polymorphisms and Their Association with Milk Production Traits in Dairy Cattle

The correlation of DNA polymorphisms with milk production in dairy cattle has been studied for several genes, including SCD, prolactin, DGAT1, leptin, GHR, CSN1S1, ABCG2, GH etc. In Table 3, we have summarized all the major DNA polymorphisms in genes and their association with milk production traits in dairy cattle.
Stearoyl-CoA desaturase 1 (SCD1) located on chromosome 26 has been widely studied for its association with milk production traits in dairy cattle [57,58,59,60,61,62]. Taniguchi et al. studied the polymorphisms of SCD in Holstein-Friesian, Jersey, Brown Swiss, and Japanese black cattle breeds and found their association with milk fat composition [60]. Similarly, Kgwatala et al. documented the SNP at 3-UTR of SCD and their link with milk fatty acids in Canadian Holstein and Jersey breeds [62]. Consequently, Macciotta et al. reported that the Italian Holsteins with VV genotypes produced more milk and protein than those with AA genotypes. In contrast it has been reported that cows with AA genotypes were producing more milk fat [61]. Furthermore, they highlighted that because of the involvement of SCD gene in energetic pathways, it might be the reason for their association with milk production traits such as milk, yield, and protein [62]. Mele et al. [59] studied the genotypic effect of SCD on milk fatty acids in 297 Italian Holstein Friesian cows. The genotypes in SCD were confirmed through the single-strand conformation polymorphism method. Interestingly, they found that cows having AA genotypes producing more milk fat compared to VV genotypes cows [60]. The above results were also verified by a recent study who found that heterozygous genotypes Chinese Holsteins were producing more milk than the cows of homozygous genotypes [63]. Similarly, Kesek-Wozniak et al. reported that heterozygous genotypes Polish Holstein-Friesian cows produced more milk fatty acids in milk compared to VV and AA genotypes cows [64].
Alim et al. [63] reported several SNPs in the SCD gene and their association with milk production traits in Chinese Holsteins. They documented that polymorphism g.10329C/T at exon 5 changed the amino acid alanine to valine. In addition, it was noticed that the two SNPs (g.6926A/G and g.8646A/G) at intron 3 and three polymorphisms (g.10153A/G, g.10213T/C and g.10329C/T) at exon 5 in SCD were significantly associated with milk fat, milk yield, protein yield and protein (%) in Chinese Holsteins [63]. Recently, it has been documented that A293V (c.878C/T) mutation in SCD changed the amino acid alanine to valine and is associated with milk fatty acid in Polish Holstein-Friesian cows [64]. Constantly, Bouwman et al. reported that the A allele of SNP in SCD was associated with higher milk fatty acids [17], while other studies found the less effect of V allele on milk fats in White Fulani and Borgou cattle breeds [65,66].
The polymorphism (DGAT1 K232A) in DGAT1 has been widely studied for its association with milk production traits particularly milk fatty acids in dairy cattle (Table 2) [64,66,67]. In addition, it has been documented that K allele is linked to high milk fat yield, fat content, and protein content and lower milk production protein and lactose yield [68,69]. While other studies reported that cows with AA genotypes have higher milk yield and lactose yield and low milk fat and protein contents [70,71,72]. Based on the above findings, it can be concluded that the DGAT1 K232A can be a target as a useful genetic marker for milk production improvement in dairy cattle.
Raschia et al. conducted a study in Argentina and reported several SNPs in selected candidate genes for their association with production traits in 20 Holstein and 5 Jersey cows. Furthermore they documented that the genes including ARL4A (rs43375517), SCD1 (rs41255693), GH (rs109191047; rs137651874), leptin (rs29004488), OPN (rs132812135), PRLR (rs135164815) and LTF (rs43706485) were associated with milk production traits [73]. Moreover, the SNPs (SNP1; G43737229T, SNP2; G43737229T, SNP3; G43761121A, SNP4; G43761121A, SNP5; G43761121A, SNP6) in breast cancer 1 (BRCA1) gene were significantly associated with milk yield in Karan Fries (Bos taurus × Bos indicus) cows [74]. However the stage of lactation (lactation 1, 2 and 3) and genotypes (GG, TT, TG) of cows were the key factors that affected the effect of SNPs on milk yield in Karan Fries cows [74].
Fatty acid desaturase 2 (FADS2) is another promising candidate gene that influences milk fatty acid traits and is located on bovine chromosome 29, with 16 exons encoding 359 amino acid chains [75]. The polymorphism in the FADS2 gene has been widely studied for its association with milk fatty acids (MFAs) in dairy cattle [43,76,77]. Based on published data, it can be recommended that the FADS2 can be a useful candidate marker for milk fat traits improvement in dairy cattle. The detail of FADS2 gene and their polymorphisms has been given in Table 3.
Ahmed et al. documented few milk protein genes (CSN1S1, CSN2, CSN1S2, CSN3, LALBA, and LGB) that were linked to increase milk protein traits in Sudanese Butana cattle [78]. Consistently, CSN1S1, CSN2CSN2, CSN2, CSN1S2, LALBA genes have been studied for their association with milk production traits in other Bos indicus breeds such as Sarabi, Sistani, Golpayegani and Gir in Iran and Brazil, respectively [79,80]. Moreover, Miluchová et al. proved experimentally that the CSN3 gene was significantly associated with milk production traits in the Slovakian Holstein population [81].
Haruna et al. had documented that the myostatin gene was significantly associated with increasing the amount of milk unsaturated fatty acid and decreasing the amount of saturated fatty acid in New Zealand Holstein-Friesian cross Jersey-Cross Cows. Moreover, they reported that cows with AD genotypes were linked to decreased saturated fatty acid while cows with AA genotypes correlated with increased milk unsaturated fatty acid [82]. Consequently, it has been documented that fatty acid-binding proteins (FABPs) is significantly associated with milk fatty acids synthesis in Holstein-Friesian × Jersey (HF × J) dairy cows [83].
The glycosylphosphatidylinositol-anchored high-density lipoprotein-binding protein1 (GPIHBP1) is a key gene that has been studied for its association with milk fat (%) and milk protein yield [32,84]. Consistently, another study had reported that GPIHBP1 was significantly correlated with milk fat traits in Chinese Holsteins [85]. They demonstrated that when the expression of GPIHBP1 was decreased, which decreased the LPL binding ability to GPIHBP1 and alternatively, the process of lipolysis was inhibited in mammary epithelial cells, resulting in increased fat in milk [85]. Moreover, Dong et al. [86] illustrated that the decrease of the expression of GPIHBP1, result in an increase in milk protein genes (CSN1S1, CSN1S2, CSN2, and CSN3, lactoferrin) which were associated with the regulation of milk protein biosynthesis [86].
Long-chain acyl-CoA synthetase 1 (ACSL1) is located on chromosome 27 of cattle (Bos Taurus), having 20 exons, 19 introns with 64,883 bp length [87]. The SNPs detected in the ACSL1 gene were genotyped in 992 Chinese Holstein cows and documented the significant association of these SNPs with milk production traits [87]. Consistently a study also documented the up-regulation of sic genes (ACACA, GPAM, ACSL1, FASN, LPIN1 and ACSL6) in dairy cattle during lactation [88]. Twenty candidate genes associated with milk fatty acid traits in Chinese Holstein cows were identified in a previous study, and ACSL1 was one of them [30]. Furthermore, a study had documented that mutation in the ACSL1 gene plays a key role in the milk fat enhancement of Yak [89]. Fan et al. experimentally proved that the expression (increase and decrease) of ASCL4 was significantly associated with milk fat synthesis in bovine mammary epithelial cells [90]. The interaction of ASCL4 was reported with ASCL1, FADS2, FASN, PPARD, CPT1A, FABP3 and ELOVL6 which are key genes associated with milk production traits. Based on the above findings, we can conclude that ASCL1 can be considered a key regulator of milk fat synthesis.
Acylglycerol-3-phosphate O-acyltransferase 3 (AGPAT3), located on Bos taurus autosome 1 (BTA1) having eight exons encoding 376 amino acid chains, is a crucial acyltransferase that is involved in triglyceride (TG) and phospholipid biosynthesis [91]. AGPAT3 has been identified through GWAS studies as a positional candidate gene affecting milk fatty acids in dairy cattle [30,91,92]. A study by using GWAS study documented the AGPAT3 and was found to be significantly linked with milk fatty acid traits in Chinese and Danish Holstein populations [39]. Recently, a study detected a SNP1 (g.12264 C > T) at promoter region, SNP2 (g.18852 C > T) in exon 5 and other six SNPs (g.18658 G > A, g.20046 G > A, g.23034 C > A, g.28332 C > T, g.28484 C > T, and g.28731 A > G) on intronic regions of AGPAT3 in dairy cattle [93]. All the SNPs reported by Sun et al. showed significant association with at least one phenotypic trait of milk production. Similarly, Shi et al. reported 17 SNPs in AGPAT3 that were associated with milk fatty acid traits in Chinese Holstein cows [94]. Littlejohn et al. also documented several SNPs of AGPAT3 in Holstein-Friesian × Jersey crossbreed that were associated with milk fat synthesis [95]. The detail of SNPs in AGPAT3 has been given in Table 3.
Prolactin is another key gene having an important role in lactation initiation and maintenance in mammals [96]. Several polymorphisms within prolactin have been identified which were significantly associated with milk production traits in dairy cattle [96,97]. Poglo et al. identified BACH2, E2F3 and KDM5A as key genes that are involved in the regulation of milk fat synthesis in the mammary gland of dairy cattle [98].
Table 3. The variations in genes and their association with milk production traits in dairy cattle.
Table 3. The variations in genes and their association with milk production traits in dairy cattle.
Gene (Location)Polymorphism (Location)Change in Amino Acid SequenceProduction TraitBreed CountryAuthor
SERPINA1 (BTA21)rs208607693 (5- flanking region)
rs210222822 (Exon-2)
MY, FY, PY, PPChinese HolsteinChina[99]
SERPINA1 (BTA21)rs41257068 (Exon-2)
rs207601878 (Intron-3)
MY, FY, PY, PP, FPChinese HolsteinChina[99]
SCD (BTA21)c.878C/T (Exon5)p.A293VMFAsHolstein Friesian × Jersey dairy cowsNewzealand[100]
SCD (BTA21)c.1783A/G (3-UTR) MF Holstein Friesian × Jersey dairy cowsNewzealand[100]
SCDg.10329C > T Alanine to valineMY, FY, PY, PPChinese HolsteinChina[101]
g.10153G > A
g.10213T > C
SCDg.10329C > T Alanine to valineMFACanadian Jersey cowsCanada[62]
SCDg.10329C > T Alanine to valineFY, PY, MYHolstein cowsBelgium[102]
SCDA293V SNP MFAsPolish
Holstein-Friesian
Poland[64]
SCDA293V SNP MFAsItalian HolsteinsItaly [59]
CD4 (BTA5)g.13598C > T MY, FY, PYChinese HolsteinChina[103]
STAT5B (BTA19)g.31562T > C; Exon 16 MY, PYChinese HolsteinChina[103]
DDIT3 (BTA5)g.56283814C > T;5-flanking region
g.56284880C > T;5-flanking region
FP, FY, PPChinese HolsteinChina[15]
c.*21A > G (5-UTR) MY, FY, PY China[15]
RPL23A (BTA19)g.20702088A > G 5-flanking region
g.20702122C > G 5-flanking region
g.20702782_83insG 5-flanking region
MY, FY, PYChinese HolsteinChina[15]
SESN2 (BTA19)g.125716884A > G, 5-flanking region
g.125714860_125714872del, 5-flanking region
g.125714806delinsCCCC,
5-flanking region
g.125714850A > G, 5-flanking region
g.125716686A > G,
5-flanking region
MY, FY, PYChinese HolsteinChina[15]
NR4A1 (BTA5)g.27994068A > G, 5-flanking region
g.27993737A > G, 5-flanking region
g.27992897C > T, 5-flanking region
c.*138A > G (5-UTR)
MY, FY and PYChinese HolsteinChina[15]
PTK2 (BTA14)g.4061098T > G (Exon5)p.Ile981MetMY, PY and FPChinese HolsteinChina[104]
g. 3895208T > G (Intron2) MY, PY and FP
g. 4059863 A > C (Intron13) MY, FY, PY and FP
g.3968605A > G (Intron6) MY, PY and FP
g.7012367T > C (Intron16) MY, FP and PP
UGDH (BTA6)rs61000233G/A
(Exon1)
rs60966191A/T
(Exon12)
MYChinese HolsteinChina[105]
SAA1 (BTA11)g.-1788C > T (Promoter)
g.-963C > A (Promoter)
g.-781 A > G (Promoter)
PY, MYChinese HolsteinChina[106]
c. + 2510A > G (EXON3)Gly48AspMY, FP and PY
c. + 2535C > T(EXON3)R56RPY, MY
c. + 2565G > A (EXON3)P66PPY, MY
SAA2 (BTA29)c.-22G > A (Promoter)
c.17G > C (Promoter)
c.114G > A (Promoter)
MY, FY and PYChinese HolsteinChina[107]
ACACB (BTA17)g.66218726T > C (Promoter)
g.66218117G > A(Promoter)
milk production traitsChinese HolsteinChina[108]
ERBB2 (BTA19)g.22400A > G (Intron-23)
g.22346A > T (Intron-23)
g.16431C > G (Intron-14)
g.19414A > G (Intron-14)
g.11680C > T (Intron-8)
g.10727A > G (Intron-7)
g.23650T > C (Intron-26)
g.22268T > C (Exon 23)
g.20982del (Intron-19)
Milk PPChinese HolsteinChina[109]
ERBB2g.873T > C (5-flanking region)
g.21561A > G (Exon-21)
MY, FY, FP, PY, PPChinese HolsteinChina[109]
HSPA8 (BTA15)rs132976221
g.4218T > G (Intron-3)
MY, PY, FYChinese HolsteinChina[109]
HSPA8rs136632043
g.4218T > G-exon9
MY, PY, FYChinese HolsteinChina[109]
ECHS1 (BTA6)g.25858322C > T-exon 3Leucine -phenylalanineMFAsChinese HolsteinChina[110]
ECHS1g.25857784C > T(exon 2) MFAChinese HolsteinChina[110]
FADS2 (BTA29)c. 908 C > T (Exon 7)294Ala > ValMY, PY, FY, FPChinese HolsteinChina[111]
FADS2c.1571 G > A (3-UTR) MYChinese HolsteinChina[111]
FADS2c.1571G > A (3-UTR) MFAChinese HolsteinChina[75]
FADS2rs209202414 G > A MFAHS and JerseyRomania and Poland[77]
FADS2rs211580559 G > A(exon 7)(294 Ala > Val)MFA [76]
rs42187261 G > A (exon 8) MFA [76]
rs109772589 G > A (3-UTR)
rs136261927 G > A (3-UTR) rs109772589 G > A (3-UTR)
MFACanadian HolsteinCanada[76]
THRSP (BTA29)rs42714482 (exon 7)Ala51ValMFAJersey and Polish HolsteinPoland[112]
SCD 1 (BTA26)rs41912290 G > A (Exon16)Leu/ProMFAHolsteinUS[113]
rs41255691 (Exon 5) HolsteinUS[113]
SCD5 (BTA6)ss252452201 (exon 3)
ss252452202(exon 3)
ss252452203(exon 3)
rs43687655 (exon 4)
MFAHolsteinUS[113]
INSIG1 (BTA4)ss252452218 (Exon 1)
ss252452220 (exon4)
ss252452222 (exon5)
Ser/Gly
Leu/Phe
MFAHolsteinUS[113]
INSIG2(BTA4)ss252452227 (5′UTR)
ss252452228 (5′UTR)
ss252452229 (5′UTR)
MFAHolsteinUS[113]
MBTPS1 (BTA18)ss252452238 (Exon 1) MFAHolsteinUS[113]
MBTPS2 (BTA18)ss252452240 (5′UTR) MFAHolsteinUS[113]
SCAP (BTA12)ss252452209 (5′UTR)
ss252452210 (5′UTR)
ss252452212 (exon 3)
ss252452215 (exon 7)
ss252452217 (exon 7)
rs41255691 (exon5)
MFAHolsteinUS[113]
SCAP (BTA12)ss252452215 (exon 7)Pro/SerMFAHolsteinUS[113]
MAP4K4 (BTA11)c.2061T > G (exon 18) PP, MYChinese HolsteinChina[114]
IGF2Rg.72479 G > A (exon 24) MY, PY, LCPolish HolsteinPoland[115]
Sirtuins (BTA28)g.-274C > G (Promoter) MY, FY, FP, LLChinese Red Steppe
Agerolese,
Qinchuan
Nanyang, Jiaxian, Luxi
Italy[116]
SCDSCD1-A293VA293VMFAHolsteinNetherlands[65]
FASN (BTA19)rs41919999 (Intron 12)
rs41919992 (Exon-27)
rs133498277 (Intron 28)
rs41919984 (Exon 37)
rs41919986 (Exon-42)
MCFAsChinese HolsteinChina[117]
FASNrs41919985 (Exon-39)alanine 2266threonineMCFAsChinese HolsteinChina[117]
PPARGC1A (BTA6)rs109579682-Intron-9 MCFAsChinese HolsteinChina[117]
ABCG2 (BTA6) rs137757790-Intron-7 MCFAsChinese HolsteinChina[117]
IGF1 (BTA5)rs109763947-5′-UTR MCFAsChinese HolsteinChina[117]
ABCG2ABCG2-Y581S MY, FP, PPIranian HolsteinIran[118]
LEPR (BTA3)LEPR-T945M MY, FP, PPIranian HolsteinIran[118]
SCD1SCD1-A293V MY, FP, PPIranian HolsteinIran[118]
DGAT1 (BTA14)DGAT1 K232A MFAModicana cowsItaly[119]
PRLR (BTA14)g.38948871C > T (5′
flanking region)
g.38949011G > A (5′
flanking region)
g.39115345 T > C(Exon4)
MFA [120]
PRLRg.39115344G > A (Exon4)Serine- asparagineMFA [120]
CHUK (BTA26)g.21008688G > T (5′
flanking region)
g.20966385C > G (3′ UTR)
g.20966354C > T(3′ UTR)
MFAChinese HolsteinChina[120]
MOGAT1 (BTA2)g.111599360A > G T (5′
flanking region)
g.111601747 T > A T (5′
flanking region)
MFAChinese HolsteinChina[120]
MINPP1 (BTA26)g.9206582C > T (3′ UTR)
g.9207070A > G (intron 5)
MFAChinese HolsteinChina[120]
CPM (BTA5)g.45079507A > G-5′
flanking region
rs208252716
g.45080228C > A-5′
flanking region
rs109638242
g.45080335C > G-5′
flanking region
rs136799678
g.45162113G > A 3-UTR
rs134841257
g.45163633G > T 3- flanking region
rs110822514
g.45164215A > G-3- flanking region
rs462818932
g.45164996A > G-3- flanking region
rs382501675
MFAChinese HolsteinChina[121]
CPMapmap49848-BTA-106779
(Intron 2)
MFAChinese HolsteinChina[30]
LIPK (BTA26)g.10428101G > A
rs110322221 (5-UTR)
g.10449831C > A
rs42774527 (Exon 11)
g.10214117A > C
rs41606812 (5′
flanking region)
g.10217380C > A
rs211373799 (5-UTR)
g.10247997T > C
rs42107056 (3′ UTR)
MFAsChinese HolsteinChina[122]
LIPJ (BTA26)g.10250098C > T (3- flanking
region)
rs42107122
g.10250120A > G
ss158213049726 (3- flanking
region)
g.10251075G > T
rs209219656 (3- flanking
region)
g.10251111T > C
rs42107114 (3- flanking
region)
MFAsChinese HolsteinChina[122]
DGAT1 (BTA14)DGAT1 (K232A) MFYDutch HolsteinIsrael[123]
DGAT1DGAT1 (K232A) MPTsHolsteinUSA[124]
DGAT1DGAT1 (K232A) MPTsIrish HolsteinIreland[125]
DGAT1rs109663724 PY, FY, MY USA[126]
DGAT1rs132699547
rs135423283
rs135576599
rs13675432
PY, FYHolsteinUSA[126]
DGAT1DGAT1 (K232A) MPTs + milk
coagulation properties
Italian HolsteinItaly[127]
DGAT1DGAT1 (K232A) MFY, MPY, MFAsHolsteinNetherlands[128]
DGAT1DGAT1 (K232A) MPTs, LCHolsteinNetherlands[68]
DGAT1DGAT1 (K232A) MFAsHolstein–Friesian, Jersey, Frisón
Negro, Montbeliarde, and Overo
Colorado
Chile[67]
DGAT1DGAT1 (K232A) MPTsHolsteinNetherlands[65]
DGAT1DGAT1 (K232A) MPTsBorgou and White Fulani cattleBenin[66]
DGAT1DGAT1 (K232A) MPTsHolsteinCzech Republic[129]
DGAT1DGAT1 (K232A) MFAsPolish HolsteinPoland[64]
DGAT1DGAT1 (K232A) MPTsPolish HolsteinPoland[130]
DGAT1DGAT1 (K232A) Milk metabolome and proteomeHolsteinNetherlands[131]
DGAT1DGAT1 (K232A)
MFAs
HolsteinNetherlands[132]
DGAT1DGAT1 (K232A) FY, PYHolsteinGermany[133]
DGAT1DGAT1 (K232A) MPTsPolish HolsteinPoland[134]
DGAT1DGAT1 (K232A) FP, MFAsRomanian HolsteinRomania[135]
DGAT1DGAT1 (K232A) MFAsHolsteinNetherlands[69]
DGAT1DGAT1 (K232A) MFAsHolsteinNetherlands[136]
DGAT1DGAT1 (K232A) FY, PY, MY HolsteinNetherlands[71]
DGAT1DGAT1 (K232A) MFAsHolsteinNetherlands[137]
DGAT1DGAT1 (K232A) MPTsHolstein cross Normande France [70]
DGAT1DGAT1 K232A FC, PCindigenous Ongole cattle, Indian Jersey, HolsteinIndia[138]
DGAT1rs109421300 G > A (5′flanking region) FPHolsteinChina[1]
JAK2 (BTA8)rs210148032 C > T(exon 16) PPHolsteinChina[1]
JAK2JAK2/RsaI (rs110298451) PYPolish Holstein, Montbeliarde,
Simmental, Jersey
Poland[139]
ELOVL6 (BTA6)g16379651A > G ELOVL6-Intron 3
g16458976A > G ELOVL6-Intron 3
MYHolsteinChina[3]
g16511290A > G ELOVL6-3′ UTR MY, MFC (%)HolsteinChina[3]
ACSL1 (BTA27)5′UTR-ACSL1-g.20523C > G
ACSL1-g.35651G > A- Intron 2
PC (%)HolsteinChina[140]
ACSL1ACSL1-g.35446C > T- Intron 2
ACSL1-g.35651G > A- Intron 2
g.51472C > T- ACSL1-Intron 11
TDMY (kg), FC (%), PC (%)HolsteinChina[140]
AGPAT6 (BTA27)BovineHD2700010331
g.36175805C.T (Intron 1)
ARS-BFGL-NGS-57448 g.36155097C.T (5-UTR exons)
MFYNew
Zealand Holstein-Friesian cross Jersey
Newzealand[95]
AGPAT3 (BTA1)g. 146702957G > A
rs210638665 (5′ flanking region)
g. 146704373A > G
rs209442459 (5′ flanking region)
g. 146704618A > G
rs110551271 (5′ flanking region)
g. 146704699G > A
rs110278717 (5′ flanking region)
MFAsHolsteinChina[94]
AGPAT3g.28731 A > G
g.12264 C > T
PP (%), FP (%), MYHolsteinChina[93]
ATPase 6 (BTA20)m.8308A > G (Exon 2) MFYHolsteinChina[141]
BoLA (BTA23)BoLA-DRB3.2 (exon 2) Milk microbiotaHolsteinCanada[142]
CSN3 (BTA6)g.10993T > A (exon 4)
g.10888T > C (exon 4)
g.10924C > A (exon 4)
g.10985G > A (exon 4)
Iso-Thr
Ala-Asp
Ala-Ala
PP (%), FP (%)HolsteinChina[63]
CSN3g.10944A > G(exon 4)Serine > GlycineMY, FY, FP (%),
PY
HolsteinChina[63]
CSN3g.12703T > G MY, FY, PY, PP(%)HolsteinChina[63]
GHR (BBTA20)GHR-F279Y polymorphismphenylalanine to tyrosineMY, PC (%), FC (%), LC (%)HolsteinGermany[143]
HSP90AB1 (BBTA23)SNP g.4338T > C MYFrieswal, SahiwalIndia[144]
HTR1B (BTA9)rs207969357
g.17303383G > T (Exon 1)
Alanine to SerineMFAsHolsteinChina[145]
HTR1Brs207969357
g.17307103A > T (Promoter)
rs476055046
g.17305206 T > G (Promoter)
rs476055046
g.17303761C > T (Promoter)
rs208945882
g.17303042C > G (Exon-1)
MFAsHolsteinChina[145]
IGF2 (BTA29)g.8656C > T-Exon 2
g.24507G > T-Exon 10
FY, PY, PCPolish HolsteinPoland[146]
IGF2rs42196909
IGF2.g-3815A > G
rs42196901
MY, PYIrish HolsteinIreland[147]
IGFBP2 (BTA9)rs133488718-Intron 3
rs133235938- Intron 5
MYHolsteinUK[148]
PDE9A (BTA2)c.-2012 T > C (rs42140305)
c.-2005 A > G (rs381951806)
MY, PYHolsteinChina[149]
LAP3 (BTA6)g.24564G > A (ss196003366)
g.24794T > G
g.24803T > C
g.24846T > C
g.25415T > C
MFP (%), PP (%)HolsteinChina[150]
Lipin 1 (BTA11)g.86129263C > G
rs211527179-5-flanking region
FY, FP (%), PYHolsteinChina[151]
Lipin 1c.637T > C
rs110871255-exon 5
Methionine-thrMY, FP (%),PY HolsteinChina[151]
Lipin 1c.708A > G
rs110161110-exon 5
Thr-AlaMY, PYHolsteinChina[151]
Lipin 1c.1521C > T
rs207681322-exon 8
Proline-SerineMY, FY, FP (%), PYHolsteinChina[151]
Lipin 1c.1555A > C
rs137642654-exon 8
Histidine-ProlineMY, FP (%),PYHolsteinChina[151]
Lipin 1g.86049523C > T
rs135886289-3-flanking region
g.86049389C > T
rs109039955-3-flanking region
MY, FY, PY, PP (%)HolsteinChina[151]
Leptin (BTA4)accession number MN119554 SNPp.Ala80ValMFAsHolstein Friesian × Jersey dairy cowsNewzealand[152]
MBL (BTA26)g. 2686T > C-Exon 2
g.2651G > A-Exon 2
FC (%), PC (%)Bohai Black Chinese Holstein
Luxi Yellow
China[153]
MBLg.1164 G > A (Exon 3)Proline-GlutamineFC (%), PC (%)Chinese HolsteinChina[154]
OLR 1 (BTA5)SNP10497 A > C (3′ UTR) PC (%)Israeli-HolsteinIsrael[155]
PRKG1 (BTA26)g.8344262A > T
rs109571301-5′ flanking region
g.6904047G > T
rs478962267-3′ UTR
g.6903810G > A
rs444193880-3′ UTR
g.6903365C > A
rs42630538-3′ UTR
g.6902878 T > G
rs136888798-3′ UTR
g.6901713 T > G
rs381717383-3′ flanking region
MFAsChinese HolsteinChina[156]
Prolactin (BTA23)-1043A > G (Promoter)
-402A > G (Promoter)
+ 8398G > A (Exon 4)
MY, MFC (%)Chinese HolsteinChina[96]
Prolacting.7545G > A(Intron 4) MYChinese HolsteinChina
[97]
SCAP (BTA23)ss526061914 (5-UTR/Exon 1) MY, MFAsHolsteinUS[157]
INSIG1 (BTA4)ss526061846 (Exon 4)L852PMY, MFAsHolsteinUS[157]
SREBF1 (BTA19)ss526061830 (Exon 14) MFAsHolsteinUS[157]
SLC27A1(BTA1)SNP-112T > C MYChinese HolsteinChina[158]
SLC27A6 (BTA1)g.390C > T
ss672469900-Exon1
g.15975T > C
ss672469898-Exon2
MFAsHolsteinUS[159]
SLC27A6g.242A > T
ss672469901-Exon1
Lysine81methionineMFAsHolsteinUS[159]
FABP4 (BTA14)g.3711G > C
ss672469893-Exon3
g.3691G > A
ss672469894-Exon3
Valine110MethionineMFAsHolsteinUS[159]
TLR4 (BTA8)c.-226 G > C
rs 29017188
c.2021 C > T
rs 8193069
MY, FC (%), PC (%), LCA (%)Chinese HolsteinChina[160]
TLR4c.2021 C > TThreonine-isoleucineMY, FCChinese HolsteinChina[161]
Lactoferrin (BTA22)SNPs -270 T > C
SNP-156 A > G
FC (%), PC (%), LCA (%)Chinese HolsteinChina[162]
Transferrin (BTA1)g.-1748 G > A
ss250608649-5-Flanking region
g.14037A > G
ss250608651-Exon 8
MY, PC (%)China Holstein Luxi Yellow Bohai BlackChina[163]
SCL2A12 (BTA1)g.72224078C > G--5-Flanking region MY, FC (%), PY, LY, LC (%)Polish HolsteinPoland[164]
SCL5A1 (BTA1)g.70571253A > G (Promoter) MY, FC (%), PY, LY, LC (%)Polish HolsteinPoland[164]
Leptin (BTA4)rs29004509-exon 3 MY indicine and taurine crossbred (Karan Fries) [165]
LALBA (BTA5)g.31183170T > C-promoter region LY, MYPolish HolsteinPoland[166]
PIK3R1 (BTA20)g.4453141 T > G-5′flanking region
rs207593520
c.1505G > A-3′UTR
rs208460068
g.4448024C > T-3′flanking region
rs209154772
g.4447105C > G-3′flanking region
rs210000760
MY, FC (%), PY, FY, PC (%)Chinese HolsteinChina[6]
PIK3R1 (BTA20)c.208G > A-5′UTR
rs42590258
GG 0.41
c.2776 T > C-3′UTR
rs210389799
c.2962 T > C-3′UTR
rs208819656
c.6275 T > A-3′UTR
rs41255622
g.11323546C > T-3′flanking region
rs133655926
g.11323118G > A
rs211408208-3′flanking region
MY, FC (%), PY, FY, PC (%) [6]
GH (BTA19)GHp.L127V (Exon 5) MFAs, MYModicana cowsItaly[167]
TMEM232 (BTA7)
HCN4 (BTA9)
ATP8A2 (BTA12)
LOC524642 (BTA29)
LOC524642
rs43708473
rs110025880
rs109784719
rs42169108
rs43099931
MYHolsteinGermany[168]
HAL (BTA5)ss974768522 (Promoter)
ss974768523 (Exon 1)
ss974768527 (Exon8)
ss974768525 (Exon5)
N42N
I156I
Gly228Glu
Milk production traitsChinese HolsteinChina[169]
GALE (BTA2)g.2114A > G-5′-UTR
ss1996900612
g.2037G > A-5′-UTR
ss1996900613
g.3836 G > C-Introm 9
rs211659075
MY, FY, PY, PPChinese HolsteinChina[170]
MY: milk yield; TDMY: Test-day milk yield; PC: Protein content; MFY: Milk fat yield; MFC (%): Milk fat content (%); MPTs: Milk protein traits; LY:Lactose Yield; LC: Lactose contents; FP: fat percentage, PY: protein yield, PP: protein percentage, MFA: Milk fatty acids; LC: lactose content; LL: Lactation length; milk medium-chain fatty acids (MCFAs); 3′ untranslated region (UTR); BTA21: Bos Taurus autosomal chromosome 21; Serpin peptidase inhibitor, clade A (SERPINA1); histidine ammonia-lyse gene (HAL); (Sterol regulatory element-binding protein-1 (SREBP1); UDP-galactose-4-epimerase (GALE); acetyl-CoA carboxylase beta (ACACB); SREBP cleavage-activating protein (SCAP); insulin-induced protein 1 (INSIG1);insulin-induced protein 2 (INSIG2); membrane-bound transcription factor protease, site 1 (MBTPS1); membrane-bound transcription factor protease, site 2 (MBTPS2); stearoyl coenzyme-A desaturases (SCD1 and SCD5); thyroid hormoneinducible hepatic protein gene (THRSP); UDP-glucose dehydrogenase (UDPH); Mitogen activated protein kinase kinase kinase kinase (MAP4K4); insulin-like growth factor receptor 2 (IGF2R); fatty acid synthase (FASN); peroxisome proliferator-activated receptor gamma, coactivator 1 alpha (PPARGC1A); ATP-binding cassette, sub-family G, member 2 (ABCG2); insulin-like growth factor 1 (IGF1); leptin receptor (LEPR); 1-acylglycerol-3-phosphate Oacyltransferas 3 (AGPAT3); Conserved helix-loop-helix ubiquitous kinase (CHUK); Multiple inositol-polyphosphate phosphatase 1 (MOGAT1); Multiple inositol-polyphosphate phosphatase 1 (MINPP1); carboxypeptidase M (CPM); enoyl-CoA hydratase, short chain 1 (ECHS1); Fatty acid desaturase 1 (FADS1) and 2 (FADS2); lipase family member K (LIPK); lipase family member J (LIPJ); Janus kinase 2 (JAK2); ELOVL6:Fatty Acid Elongase 6; ACSL1:acyl-CoAsynthetase 1; bovine major histocompatibility complex (BoLA); Iso:Isoleucine; Thr: Threonine; Ala: Alanine; Asp: Aspartic acid; CSN3:kappa casein; GHR: growth hormone receptor; HTR1B:hydroxytryptamine receptor 1B; IGF2: Insulin-like growth factor 2; IGFBP2:Insulin-like growth factor binding protein-2; PDE9A:Phosphodiesterase9A; MBL: Mannan-binding lectin; OLR 1: oxidized low-density lipoprotein (lectin-like) receptor 1; PRKG1: protein kinase, cGMP-dependent, type I; SREBF1: sterol regulatory element binding transcription factor 1; SLC27A1: solute carrier family 27 member 1″ protein;TLR4: Toll-like receptor 4; LALBA: Alpha-lactalbumin; phosphoinositide-3-kinase regulatory subunit 1 (PIK3R1); dual specificity phosphatase 1 (DUSP1); Growth hormone (GH);Long-chain acyl-CoA synthetase 1 (ACSL1); protein tyrosine kinase 2 (PTK2).

7. Conclusions

In the current review, we documented several genes associated with milk production traits in dairy cattle. Moreover, many SNPs within candidate genes were highlighted in the current review, which could be a useful addition to the genetic markers linked to the improvement of milk production traits in dairy cattle. There are still many candidate genes reported through GWAS studies, RNA-seq and DNA-seq need further validation in dairy cattle before selecting them as genetic markers in cattle breeding.

Author Contributions

Conceptualization, Y.M., M.Z.K. and Z.C.; writing—original draft preparation, M.Z.K., Y.M. and Z.C.; data search and collection, T.C., J.W., X.C., Z.H. and M.K.S.; editing and technical review, M.Z.K., J.X., S.L., Y.M., G.M.A. and Z.C.; visualization, Z.C.; supervision, Z.C. All authors have read and agreed to the published version of the manuscript.

Funding

The review was supported by the 2115 Talent Development Program of China Agricultural University. The funder had no role in the study design, data collection, analysis, decision to publish, and manuscript preparation.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

All the data are already provided in the main manuscript. Contact the corresponding author if further explanation is required.

Acknowledgments

We acknowledge the National Key Research and Development Program of China (2018YFD0501600), the national natural science foundation of China (U20A2062) and S & T Program of Hebei (19226625D) for their financial support. We also acknowledge the China Agricultural University, Beijing, China, for providing us with an environment of learning. Without this platform, the completion of this work would not have been an easy task.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Khan, M.Z.; Wang, D.; Liu, L.; Usman, T.; Wen, H.; Zhang, R.; Liu, S.; Shi, L.; Mi, S.; Xiao, W.; et al. Significant genetic effects of JAK2 and DGAT1 mutations on milk fat content and mastitis resistance in Holsteins. J. Dairy Res. 2019, 86, 388–393. [Google Scholar] [CrossRef]
  2. Khan, M.Z.; Liu, L.; Zhang, Z.; Khan, A.; Wang, D.; Mi, S.; Usman, T.; Liu, G.; Guo, G.; Li, X.; et al. Folic acid supplementation regulates milk production variables, metabolic associated genes and pathways in perinatal Holsteins. J. Anim. Physiol. Anim. Nutr. 2020, 104, 483–492. [Google Scholar] [CrossRef]
  3. Chen, S.; Menglin, C.; Chen, T.; Yuzhuang, L.; Tian, D.; Hui, W.; Xiaolin, L. Genetic variants of fatty acid elongase 6 in Chinese Holstein cow. Gene 2017, 670, 123–129. [Google Scholar] [CrossRef]
  4. Qian, L.; Zhao, A.; Zhang, Y.; Chen, T.; Zeisel, S.H.; Jia, W.; Cai, W. Metabolomic Approaches to Explore Chemical Di-versity of Human Breast-Milk, Formula Milk and Bovine Milk. Int. J. Mol. Sci. 2016, 17, 2128. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Ibeagha-Awemu, E.M.; Kgwatalala, P.; Zhao, X. A critical analysis of production-associated DNA polymorphisms in the genes of cattle, goat, sheep, and pig. Mamm. Genome 2008, 19, 591–617. [Google Scholar] [CrossRef]
  6. Han, B.; Yuan, Y.; Shi, L.; Li, Y.; Liu, L.; Sun, D. Identification of single nucleotide polymorphisms of PIK3R1 and DUSP1 genes and their genetic associations with milk production traits in dairy cows. J. Anim. Sci. Biotechnol. 2019, 10, 81. [Google Scholar] [CrossRef] [PubMed]
  7. Laodim, T.; Elzo, M.A.; Koonawootrittriron, S.; Suwanasopee, T.; Jattawa, D. Genomic-polygenic and polygenic predictions for milk yield, fat yield, and age at first calving in Thai multibreed dairy population using genic and functional sets of genotypes. Livest. Sci. 2019, 219, 17–24. [Google Scholar] [CrossRef]
  8. Grisart, B.; Coppieters, W.; Farnir, F.; Karim, L.; Ford, C.; Berzi, P.; Cambisano, N.; Mni, M.; Reid, S.; Simon, P.; et al. Positional Candidate Cloning of a QTL in Dairy Cattle: Identification of a Missense Mutation in the Bovine DGAT1 Gene with Major Effect on Milk Yield and Composition. Genome Res. 2002, 12, 222–231. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Grisart, B.; Farnir, F.; Karim, L.; Cambisano, N.; Kim, J.-J.; Kvasz, A.; Mni, M.; Simon, P.; Frere, J.-M.; Coppieters, W.; et al. Genetic and functional confirmation of the causality of the DGAT1 K232A quantitative trait nucleotide in affecting milk yield and composition. Proc. Natl. Acad. Sci. USA 2004, 101, 2398–2403. [Google Scholar] [CrossRef] [Green Version]
  10. Cui, X.; Yali, H.; Shaohua, Y.; Yan, X.; Shengli, Z.; Yuan, Z.; Qin, Z.; Xuemei, L.; George, E.L.; Dongxiao, S. Transcriptional profiling of mammary gland in Holstein cows with extremely different milk protein and fat percentage using RNA sequencing. BMC Genom. 2014, 15, 226. [Google Scholar] [CrossRef] [Green Version]
  11. Boutinaud, M.; Galio, L.; Lollivier, V.; Finot, L.; Wiart, S.; Esquerre, D.; Devinoy, E. Unilateral once daily milking locally induces differential gene expression in both mammary tissue and milk epithelial cells revealing mammary remodeling. Physiol. Genom. 2013, 45, 973–985. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Nayeri, S.; Sargolzaei, M.; Abo-Ismail, M.K.; May, N.; Miller, S.P.; Schenkel, F.; Moore, S.; Stothard, P. Genome-wide association for milk production and female fertility traits in Canadian dairy Holstein cattle. BMC Genet. 2016, 17, 75. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Olsen, H.G.; Hayes, B.J.; Kent, M.P.; Nome, T.; Svendsen, M.; Larsgard, A.G.; Lien, S. Genome-wide association mapping in Norwegian Red cattle identifies quantitative trait loci for fertility and milk production on BTA12. Anim. Genet. 2011, 42, 466–474. [Google Scholar] [CrossRef]
  14. Cole, J.B.; Wiggans, G.R.; Ma, L.; Sonstegard, T.S.; Lawlor, T.J., Jr.; Crooker, B.A.; Van Tassell, C.P.; Yang, J.; Wang, S.; Matukumalli, L.K.; et al. Genome-wide association analysis of thirty one production, health, reproduction and body confor-mation traits in contemporary U.S. Holstein cows. BMC Genom. 2011, 12, 408. [Google Scholar] [CrossRef] [Green Version]
  15. Li, Y.; Han, B.; Liu, L.; Zhao, F.; Liang, W.; Jiang, J.; Yang, Y.; Ma, Z.; Sun, D. Genetic association of DDIT3, RPL23A, SESN2 and NR4A1 genes with milk yield and composition in dairy cattle. Anim. Genet. 2018, 50, 123–135. [Google Scholar] [CrossRef]
  16. Chamberlain, A.J.; Hayes, B.J.; Savin, K.; Bolormaa, S.; McPartlan, H.C.; Bowman, P.J.; Van Der Jagt, C.; MacEachern, S.; Goddard, M.E. Validation of single nucleotide polymorphisms associated with milk production traits in dairy cattle. J. Dairy sci. 2012, 95, 864–875. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Bouwman, A.C.; Bovenhuis, H.; Visker, M.H.; Van Arendonk, J.A. Genome-wide association of milk fatty acids in Dutch dairy cattle. BMC Genet. 2011, 12, 43. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Schopen, G.; Visker, M.; Koks, P.; Mullaart, E.; van Arendonk, J.; Bovenhuis, H. Whole-genome association study for milk protein composition in dairy cattle. J. Dairy Sci. 2011, 94, 3148–3158. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Meredith, B.K.; Kearney, J.F.; Finlay, K.E.; Bradley, G.D.; Fahey, G.A.; Berry, P.D.; Lynn, J.D. Genome-wide associations for milk production and somatic cell score in Holstein-Friesian cattle in Ireland. BMC Genet. 2013, 13, 21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Tribout, T.; Croiseau, P.; Lefebvre, R.; Barbat, A.; Boussaha, M.; Fritz, S.; Boichard, D.; Hoze, C.; Sanchez, M.-P. Confirmed effects of candidate variants for milk production, udder health, and udder morphology in dairy cattle. Genet. Sel. Evol. 2020, 52, 55. [Google Scholar] [CrossRef]
  21. Jiang, J.; Liu, L.; Gao, Y.; Shi, L.; Li, Y.; Liang, W.; Sun, D. Determination of genetic associations between indels in 11 candidate genes and milk composition traits in Chinese Holstein population. BMC Genet. 2019, 20, 48. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Sanchez, M.; Armelle, G.; Pascal, C.; Sébastien, F.; Chris, H.; Guy, M.; Marton, P.; Barbet-Leterrier, S.; Letaief, R.; Rocha, D.; et al. Within-breed and multi-breed GWAS on imputed whole-genome sequence variants reveal candidate mutations affecting milk protein composition in dairy cattle. Genet. Sel. Evol. 2017, 49, 68. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Iung, L.H.S.; Petrini, J.; Ramírez-Díaz, J.; Salvian, M.; Rovadoscki, G.A.; Pilonetto, F.; Dauria, B.D.; Machado, P.F.; Coutinho, L.L.; Wiggans, G.R.; et al. Genome-wide association study for milk production traits in a Brazilian Holstein popula-tion. J. Dairy Sci. 2019, 102, 1–10. [Google Scholar] [CrossRef] [Green Version]
  24. Atashi, H.; Salavati, M.; De Koster, J.; Ehrlich, J.; Crowe, M.; Opsomer, G.; McLoughlin, N.; Hostens, M.; Fahey, A.; Matthews, E.; et al. Genome-wide association for milk production and lactation curve parameters in Holstein dairy cows. J. Anim. Breed. Genet. 2019, 137, 292–304. [Google Scholar] [CrossRef]
  25. Buitenhuis, B.; Janss, L.L.G.; Poulsen, N.A.; Larsen, L.B.; Larsen, M.K.; Sørensen, P. Genome-wide association and biological pathway analysis for milk-fat composition in Danish Holstein and Danish Jersey cattle. BMC Genom. 2014, 15, 1112. [Google Scholar] [CrossRef] [Green Version]
  26. Ning, C.; Wang, D.; Zheng, X.; Zhang, Q.; Zhang, S.; Mrode, R.; Liu, J.-F. Eigen decomposition expedites longitudinal ge-nome-wide association studies for milk production traits in Chinese Holstein. Genet. Sel. Evol. 2018, 50, 12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Wang, D.; Ning, C.; Liu, J.-F.; Zhang, Q.; Jiang, L. Replication of genome-wide association studies for milk production traits in Chinese Holstein by an efficient rotated linear mixed model. J. Dairy Sci. 2019, 102, 2378–2383. [Google Scholar] [CrossRef] [Green Version]
  28. Poulsen, N.A.; Robinson, R.C.; Barile, D.; Larsen, L.B.; Buitenhuis, B. A genome-wide association study reveals specific transferases as candidate loci for bovine milk oligosaccharides synthesis. BMC Genom. 2019, 20, 404. [Google Scholar] [CrossRef]
  29. Ariyarathne, H.; Correa-Luna, M.; Blair, H.; Garrick, D.; Lopez-Villalobos, N. Identification of Genomic Regions Associated with Concentrations of Milk Fat, Protein, Urea and Efficiency of Crude Protein Utilization in Grazing Dairy Cows. Genes 2021, 12, 456. [Google Scholar] [CrossRef]
  30. Li, C.; Sun, D.; Zhang, S.; Wang, S.; Wu, X.; Zhang, Q.; Liu, L.; Li, Y.; Qiao, L. Genome Wide Association Study Identifies 20 Novel Promising Genes Associated with Milk Fatty Acid Traits in Chinese Holstein. PLoS ONE 2014, 9, e96186. [Google Scholar] [CrossRef] [Green Version]
  31. Zhou, J.; Liyuan, L.; Chunpeng, C.; Menghua, Z.; Xin, L.; Zhiwu, Z.; Xixia, H.; Yuangang, S. Genome-wide association study of milk and reproductive traits in dual-purpose Xinjiang Brown cattle. BMC Genom. 2019, 20, 827. [Google Scholar] [CrossRef] [Green Version]
  32. Jiang, L.; Liu, J.-F.; Sun, D.; Ma, P.; Ding, X.; Yu, Y.; Zhang, Q. Genome Wide Association Studies for Milk Production Traits in Chinese Holstein Population. PLoS ONE 2010, 5, e13661. [Google Scholar] [CrossRef] [Green Version]
  33. Kolbehdari, D.; Wang, Z.; Grant, J.R.; Murdoch, B.; Prasad, A.; Xiu, Z.; Marques, E.; Stothard, P.; Moore, S.S. A whole genome scan to map QTL for milk production traits and somatic cell score in Canadian Holstein bulls. J. Anim. Breed. Genet. 2009, 126, 216–227. [Google Scholar] [CrossRef] [PubMed]
  34. da Cruz, A.S.; Silva, D.C.; Minasi, L.B.; Teixeira, L.K.d.F.; Rodrigues, F.M.; da Silva, C.C.; Carmo, A.S.D.; da Silva, M.V.G.B.; Utsunomiya, Y.T.; Garcia, J.F.; et al. Single-Nucleotide Polymorphism Variations Associated with Specific Genes Putatively Identified Enhanced Genetic Predisposition for 305-Day Milk Yield in the Girolando Crossbreed. Front. Genet. 2021, 11. [Google Scholar] [CrossRef] [PubMed]
  35. Yue, S.J.; Zhao, Y.Q.; Gu, X.R.; Yin, B.; Jiang, Y.L.; Wang, Z.H.; Shi, K.R. A genome-wide association study suggests new candidate genes for milk production traits in Chinese Holstein cattle. Anim. Genet. 2017, 48, 677–681. [Google Scholar] [CrossRef]
  36. Liu, L.; Jinghang, Z.; Chunpeng, C.; Juan, Z.; Wan, W.; Jia, T.; Zhiwu, Z.; Yaling, G. GWAS-Based Identification of New Loci for Milk Yield, Fat, and Protein in Holstein Cattle. Animals 2020, 10, 2048. [Google Scholar] [CrossRef] [PubMed]
  37. Kim, S.; Lim, B.; Cho, J.; Lee, S.; Dang, C.-G.; Jeon, J.-H.; Kim, J.-M.; Lee, J. Genome-Wide Identification of Candidate Genes for Milk Production Traits in Korean Holstein Cattle. Animals 2021, 11, 1392. [Google Scholar] [CrossRef]
  38. Pegolo, S.; Cecchinato, A.; Mele, M.; Conte, G.; Schiavon, S.; Bittante, G. Effects of candidate gene polymorphisms on the detailed fatty acids profile determined by gas chromatography in bovine milk. J. Dairy Sci. 2016, 99, 4558–4573. [Google Scholar] [CrossRef] [Green Version]
  39. Li, X.; Buitenhuis, A.; Lund, M.; Li, C.; Sun, D.; Zhang, Q.; Poulsen, N.; Su, G. Joint genome-wide association study for milk fatty acid traits in Chinese and Danish Holstein populations. J. Dairy Sci. 2015, 98, 8152–8163. [Google Scholar] [CrossRef] [Green Version]
  40. Otto, I.P.; Simone, E.F.; Guimarães, M.P.L.; Calus, J.V.; Marco, A.; Machado, J.C.; Panetto, C.; da Silva Marcos Vinícius, G.B. Single-step genome-wide association studies (GWAS) and post-GWAS analyses to identify genomic regions and candidate genes for milk yield in Brazilian Girolando cattle. J. Dairy Sci. 2020, 103, 10347–10360. [Google Scholar] [CrossRef]
  41. Bagnato, A.; Schiavini, F.; Rossoni, A.; Maltecca, C.; Dolezal, M.; Medugorac, I.; Sölkner, J.; Russo, V.; Fontanesi, L.; Friedmann, A.; et al. Quantitative Trait Loci Affecting Milk Yield and Protein Percentage in a Three-Country Brown Swiss Population. J. Dairy Sci. 2008, 91, 767–783. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Ouattara, B.; Bissonnette, N.; Duplessis, M.; Girard, C.L. Supplements of vitamins B9 and B12 affect hepatic and mammary gland gene expression profiles in lactating dairy cows. BMC Genom. 2016, 17, 640. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Li, Q.; Ruobing, L.; Yan, L.; Yanxia, G.; Qiufeng, L.; Dongxiao, S.; Jianguo, L. Identification of candidate genes for milk pro-duction traits by RNA sequencing on bovine liver at different lactation stages. BMC Genet. 2020, 21, 72. [Google Scholar] [CrossRef]
  44. Bai, X.; Zhuqing, Z.; Bin, L.; Xiaoyang, J.; Yongsheng, B.; Wenguang, Z. Whole blood transcriptional profiling comparison between different milk yield of Chinese Holstein cows using RNA-seq data. BMC Genom. 2016, 17, 512. [Google Scholar] [CrossRef] [Green Version]
  45. Li, C.; Wentao, C.; Chenghao, Z.; Hongwei, Y.; Ziqi, Z.; Juan, J.L.; Dongxiao, S.; Qin, Z.; Jianfeng, L.; Shengli, Z. RNA-Seq reveals 10 novel promising candidate genes affecting milk protein concentration in the Chinese Holstein population. Sci. Rep. 2016, 6, 26813. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Zhou, C.; Shen, D.; Li, C.; Cai, W.; Liu, S.; Yin, H.; Shi, S.; Cao, M.; Zhang, S. Comparative Transcriptomic and Proteomic Analyses Identify Key Genes Associated with Milk Fat Traits in Chinese Holstein Cows. Front. Genet. 2019, 10, 672. [Google Scholar] [CrossRef] [Green Version]
  47. Dai, W.; Wang, Q.; Zhao, F.; Liu, J.; Liu, H. Understanding the regulatory mechanisms of milk production using integrative transcriptomic and proteomic analyses: Improving inefficient utilization of crop by-products as forage in dairy industry. BMC Genom. 2018, 19, 403. [Google Scholar] [CrossRef]
  48. Lin, Y.; Lv, H.; Jiang, M.; Zhou, J.; Song, S.; Hou, X. Functional analysis of the dairy cow mammary transcriptome between early lactation and mid-dry period. J. Dairy Res. 2019, 86, 63–67. [Google Scholar] [CrossRef]
  49. Yang, J.; Jiang, J.; Liu, X.; Wang, H.; Guo, G.; Zhang, Q.; Jiang, L. Differential expression of genes in milk of dairy cattle during lactation. Anim. Genet. 2015, 47, 174–180. [Google Scholar] [CrossRef] [Green Version]
  50. Peterson, B.K.; Weber, J.; Kay, E.H.; Fisher, H.S.; Hoekstra, H. Double Digest RADseq: An Inexpensive Method for De Novo SNP Discovery and Genotyping in Model and Non-Model Species. PLoS ONE 2012, 7, e37135. [Google Scholar] [CrossRef] [Green Version]
  51. Weldenegodguad, M.; Popov, R.; Pokharel, K.; Ammosov, I.; Ming, Y.; Ivanova, Z.; Kantanen, J. Whole-Genome Sequencing of Three Native Cattle Breeds Originating from the Northernmost Cattle Farming Regions. Front. Genet. 2019, 9, 728. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Rosse, I.C.; Assis, J.G.; Oliveira, F.S.; Leite, L.R.; Araujo, F.; Zerlotini, A.; Volpini, A.; Dominitini, A.J.; Lopes, B.C.; Arbex, W.A.; et al. Whole genome sequencing of Guzerá cattle reveals genetic variants in candidate genes for production, disease resistance, and heat tolerance. Mamm. Genome 2017, 28, 66–80. [Google Scholar] [CrossRef] [PubMed]
  53. Sigdel, A.; Abdollahi-Arpanahi, R.; Aguilar, I.; Peñagaricano, F. Whole Genome Mapping Reveals Novel Genes and Pathways Involved in Milk Production Under Heat Stress in US Holstein Cows. Front. Genet. 2019, 10. [Google Scholar] [CrossRef]
  54. Nanaei, H.A.; Qanatqestani, M.D.; Esmailizadeh, A. Whole-genome resequencing reveals selection signatures associated with milk production traits in African Kenana dairy zebu cattle. Genomics 2019, 112, 880–885. [Google Scholar] [CrossRef] [PubMed]
  55. Ye, M.; Xu, M.; Lu, M.; Zhou, B.; El-Kader, H.A.; Alam, S.S.; Mahrous, K.F. Identification of candidate genes associated with milk yield trait in buffaloes (Bubalus bubalis) by restriction-site-associated DNA sequencing. Rev. Bras. Zootec. 2020, 49. [Google Scholar] [CrossRef]
  56. Iso-Touru, T.; Sahana, G.; Guldbrandtsen, B.; Lund, M.S.; Vilkki, J. Genome-wide association analysis of milk yield traits in Nordic Red Cattle using imputed whole genome sequence variants. BMC Genet. 2016, 17, 55. [Google Scholar] [CrossRef] [Green Version]
  57. Schennink, A.; Heck, J.M.; Bovenhuis, H.; Visker, M.H.P.W.; Van Valenberg, H.J.F.; Van Arendonk, J.A.M. Milk fatty acid unsaturation: Genetic parameters and effects of Stearoyl-CoA Desaturase (SCD1) and Acyl-CoA: Diacylglicerol Acylotrans-ferase (DGAT1). J. Dairy Sci. 2008, 91, 2135–2143. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Conte, G.; Mele, M.; Chessa, S.; Castiglioni, B.; Serra, A.; Pagnacco, G.; Secchiari, P. Diacylglycerol acyltransferase 1, stearoyl-CoA desaturase 1, and sterol regulatory element binding protein 1 gene polymorphisms and milk fatty acid composition in Italian Brown cattle. J. Dairy Sci. 2010, 93, 753–763. [Google Scholar] [CrossRef] [Green Version]
  59. Mele, M.; Conte, G.; Castiglioni, B.; Chessa, S.; Macciotta, N.P.P.; Serra, A.; Buccioni, A.; Pagnacco, G.; Secchiari, P. Stea-royl-Coenzyme A Desaturase Gene Polymorphism and Milk Fatty Acid Composition in Italian Holsteins. J. Dairy Sci. 2007, 90, 4458–4465. [Google Scholar] [CrossRef]
  60. Taniguchi, M.; Utsugi, T.; Oyama, K.; Mannen, H.; Kobayashi, M.; Tanabe, Y.; Ogino, A.; Tsuji, S. Genotype of stearoyl-CoA desaturase is associated with fatty acids composition in Japanese Black cattle. Mamm. Genom. 2004, 14, 142–148. [Google Scholar] [CrossRef]
  61. Macciotta, N.P.P.; Mele, M.; Conte, G.; Serra, A.; Cassandro, M.; Dal Zotto, R.; Borlino, A.C.; Pagnacco, G.; Secchiari, P. As-sociation between a polymorphism at the stearoyl CoA desaturase locus and milk production traits in Italian Holsteins. J. Dairy Sci. 2008, 91, 3184–3189. [Google Scholar] [CrossRef] [Green Version]
  62. Kgwatala, P.M.; Ibeagha-Awemu, M.E.; Hayes, F.J.; Zhao, X. Stearoyl-CoA desaturase 1 3′UTR SNPs and their influence on milk fatty acid composition of Canadian Holstein cows. J. Anim. Breed. Genet. 2009, 126, 394–403. [Google Scholar] [CrossRef] [PubMed]
  63. Alim, A.M.; Dong, T.; Xie, Y.; Wu, X.P.; Zhang, Y.; Shengli, Z.; Sun, X.D. Effect of polymorphisms in the CSN3 (j-casein) gene on milk production traits in Chinese Holstein Cattle. Mol. Biol. Rep. 2014, 41, 7585–7593. [Google Scholar] [CrossRef] [PubMed]
  64. Kęsek-Woźniak, M.M.; Wojtas, E.; Zielak-Steciwko, A.E. Impact of SNPs in ACACA, SCD1, and DGAT1 Genes on Fatty Acid Profile in Bovine Milk with Regard to Lactation Phases. Animals 2020, 10, 997. [Google Scholar] [CrossRef]
  65. Duchemin, S.; Bovenhuis, H.; Stoop, W.; Bouwman, A.; Van Arendonk, J.; Visker, M. Genetic correlation between composition of bovine milk fat in winter and summer, and DGAT1 and SCD1 by season interactions. J. Dairy Sci. 2013, 96, 592–604. [Google Scholar] [CrossRef] [Green Version]
  66. Houaga, I.; Muigai, A.W.T.; Ng’Ang’A, F.M.; Ibeagha-Awemu, E.M.; Kyallo, M.; Youssao, I.A.K.; Stomeo, F. Milk fatty acid variability and association with polymorphisms in SCD1 and DGAT1 genes in White Fulani and Borgou cattle breeds. Mol. Biol. Rep. 2018, 45, 1849–1862. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Carvajal, A.; Huircan, P.; Dezamour, J.; Subiabre, I.; Kerr, B.; Morales, R.; Ungerfeld, E. Milk fatty acid profile is modulated by DGAT1 and SCD1 genotypes in dairy cattle on pasture and strategic supplementation. Genet. Mol. Res. 2016, 15. [Google Scholar] [CrossRef]
  68. Bovenhuis, H.; Visker, M.; Poulsen, N.; Sehested, J.; Van Valenberg, H.; van Arendonk, J.; Larsen, L.B.; Buitenhuis, A. Effects of the diacylglycerol o-acyltransferase 1 (DGAT1) K232A polymorphism on fatty acid, protein, and mineral composition of dairy cattle milk. J. Dairy Sci. 2016, 99, 3113–3123. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Tzompa-Sosa, D.A.; van Aken, G.A.; van Hooijdonk, A.C.M.; van Valenberg, H.J.F. Influence of C16:0 and long-chain sat-urated fatty acids on normal variation of bovine milk fat triacylglycerol structure. J. Dairy Sci. 2014, 97, 4542–4551. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Vanbergue, E.; Peyraud, J.; Guinard-Flament, J.; Charton, C.; Barbey, S.; Lefebvre, R.; Gallard, Y.; Hurtaud, C. Effects of DGAT1 K232A polymorphism and milking frequency on milk composition and spontaneous lipolysis in dairy cows. J. Dairy Sci. 2016, 99, 5739–5749. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Van Gastelen, S.; Visker, M.; Edwards, J.; Antunes-Fernandes, E.; Hettinga, K.; Alferink, S.; Hendriks, W.; Bovenhuis, H.; Smidt, H.; Dijkstra, J. Linseed oil and DGAT1 K232A polymorphism: Effects on methane emission, energy and nitrogen metabolism, lactation performance, ruminal fermentation, and rumen microbial composition of Holstein-Friesian cows. J. Dairy Sci. 2017, 100, 8939–8957. [Google Scholar] [CrossRef]
  72. Mach, N.; Blum, Y.; Bannink, A.; Causeur, D.; Houee-Bigot, M.; Lagarrigue, S.; Smits, M.A. Pleiotropic effects of polymorphism of the gene diacylglycerol-O-transferase 1 (DGAT1) in the mammary gland tissue of dairy cows. J. Dairy Sci. 2012, 95, 4989–5000. [Google Scholar] [CrossRef] [Green Version]
  73. Raschia, A.M.; Juan, P.N.; Daniel, O.M.; María, J.B.; Ariel, A. Single nucleotide polymorphisms in candidate genes associated with milk yield in Argentinean Holstein and Holstein × Jersey cows. J. Anim. Sci. Technol. 2018, 60, 31. [Google Scholar] [CrossRef] [Green Version]
  74. Magotra, A.; Gupta, I.D.; Tavsief, A.; Rani, A. Polymorphism in DNA repair gene BRCA1 associated with clinical mastitis and production traits in indigenous dairy cattle. Res. Vet. Sci. 2020, 133, 194–201. [Google Scholar] [CrossRef]
  75. Li, M.; Lu, X.; Gao, Q.; Wang, M.; Arbab, A.A.I.; Sun, Y.; Chen, Z.; Zhang, H.; Karrow, N.A.; Yang, Z.; et al. A Functional 3′ UTR Polymorphism of FADS2 Affects Cow Milk Composition through Modifying Mir-744 Binding. Animals 2019, 9, 1090. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Ibeagha-Awemu, E.M.; Akwanji, K.A.; Beaudoin, F.; Zhao, X. Associations between variants of FADS genes and omega-3 and omega-6 milk fatty acids of Canadian Holstein cows. BMC Genet. 2014, 15, 25. [Google Scholar] [CrossRef] [Green Version]
  77. Proskura, W.S.; Liput, M.; Zaborski, D.; Sobek, Z.; Yu, Y.-H.; Cheng, Y.-H.; Dybus, A. The effect of polymorphism in the FADS2 gene on the fatty acid composition of bovine milk. Arch. Anim. Breed. 2019, 62, 547–555. [Google Scholar] [CrossRef] [PubMed]
  78. Ahmed, S.A.; Rahmatalla, S.R.; Bortfeldt, D.; Arends, M.; Reissmann, G.A. Milk protein polymorphisms and casein haplotypes in Butana cattle. J. Appl. Genet. 2017, 58, 261–271. [Google Scholar] [CrossRef] [PubMed]
  79. Gallinat, J.L.; Qanbari, S.; Drögemüller, C.; Pimentel, E.C.; Thaller, G.; Tetens, J. DNA-based identification of novel bovine casein gene variants. J. Dairy Sci. 2013, 96, 699–709. [Google Scholar] [CrossRef] [Green Version]
  80. Tetens, J.L.; Qanbari, S.; Drögemüller, C.; Pimentel, E.C.; Bennewitz, J.; Thaller, G.; Tetens, J. Bos indicus introgression into (peri-)alpine cattle breeds—Evidence from the analysis of bovine whey protein variants. Anim. Genet. 2014, 45, 585–588. [Google Scholar] [CrossRef]
  81. Miluchová, M.; Michal, G.; Juraj, C.; Anna, T.; Kristína, C. Association of HindIII-polymorphism in kappa-casein gene with milk, fat and protein yield in Holstein cattle. Acta Biochim. Pol. 2018, 65, 403–407. [Google Scholar] [CrossRef] [PubMed]
  82. Haruna, I.; Yunhai, L.; Ugonna, J.; Hamed, A.; Huitong, Z.; Jon, G. Associations between the Bovine Myostatin Gene and Milk Fatty Acid Composition in New Zealand Holstein-Friesian × Jersey-Cross Cows. Animal 2020, 10, 1447. [Google Scholar] [CrossRef]
  83. Zhou, H.; Cheng, P.; Azimu, W.; Hodge, S.; Edwards, G.; Hickford, J.G.H. Variation in the bovine FABP4 gene affects milk yield and milk protein content in dairy cows. Sci. Rep. 2015, 5, 10023. [Google Scholar] [CrossRef] [Green Version]
  84. Jiang, L.; Liu, X.; Yang, J.; Wang, H.; Jiang, J.; Liu, L.; He, S.; Ding, X.; Liu, J.; Zhang, Q. Targeted resequencing of GWAS loci reveals novel genetic variants for milk production traits. BMC Genom. 2014, 15, 1105. [Google Scholar] [CrossRef] [Green Version]
  85. Yang, J.; Liu, X.; Wang, D.; Ning, C.; Wang, H.; Zhang, Q.; Jiang, L. Functional validation of GPIHBP1 and identification of a functional mutation in GPIHBP1 for milk fat traits in dairy cattle. Sci. Rep. 2017, 7, 1–10. [Google Scholar] [CrossRef] [Green Version]
  86. Dong, W.; Jie, Y.; Qin, Z.; Li, J. Role of GPIHBP1 in regulating milk protein traits in dairy cattle. Anim. Biotechnol. 2018, 31, 81–85. [Google Scholar] [CrossRef] [PubMed]
  87. Bionaz, M.; Loor, J.J. ACSL1, AGPAT6, FABP3, LPINI, and SLC27A6 are the most abundant isoforms in bovine mammary tissue and their expression is affected by stage of lactation. J. Nutr. 2008, 138, 1019–1024. [Google Scholar] [CrossRef]
  88. Bionaz, M.; Loor, J.J. Gene Networks Driving Bovine Mammary Protein Synthesis during the Lactation cycle. Bioinform. Biol. Insights 2011, 5, 83–98. [Google Scholar] [CrossRef]
  89. Zhao, Z.D.; Tian, H.S.; Jiang, Y.Y.; Shi, B.G.; Liu, X.; Li, X.P.; Wang, D.Z.; Chen, J.L.; Hu, J. Association analysis of ACSL1 gene promoter polymorphism and dairy quality traits in yak. J. Agric. Biol. 2019, 27, 1596–1630. [Google Scholar]
  90. Fan, Y.; Han, Z.; Lu, X.; Zhang, H.; Idriss Arbab, A.A.; Loor, J.J.; Yang, Y.; Yang, Z. Identification of Milk Fat Metabolism-Related Pathways of the Bovine Mammary Gland during Mid and Late Lactation and Functional Verification of the ACSL4 Gene. Genes 2020, 11, 1357. [Google Scholar] [CrossRef]
  91. Olsen, H.G.; Knutsen, T.M.; Kohler, A.; Svendsen, M.; Gidskehaug, L.; Grove, H.; Nome, T.; Sodeland, M.; Sundsaasen, K.K.; Kent, M.P.; et al. 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. Genet. Sel. Evol. 2017, 49, 20. [Google Scholar] [CrossRef] [Green Version]
  92. Wang, X.; Wurmser, C.; Pausch, H.; Jung, S.; Reinhardt, F.; Tetens, J.; Thaller, G.; Fries, R. Identification and Dissection of Four Major QTL Affecting Milk Fat Content in the German Holstein-Friesian Population. PLoS ONE 2012, 7, e40711. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  93. Sun, X.; Liang, Y.; Gao, Q.; Guo, J.; Tang, C.; Shi, K.; Yang, Z.; Mao, Y. AGPAT3 Gene polymorphisms are associated with milk production traits in Chinese Holstein cows. J. Dairy Res. 2021, 88, 247–252. [Google Scholar] [CrossRef] [PubMed]
  94. Shi, L.; Wu, X.; Yang, Y.; Ma, Z.; Lv, X.; Liu, L.; Li, Y.; Zhao, F.; Han, B.; Sun, D. A post-GWAS confirming the genetic effects and functional polymorphisms of AGPAT3 gene on milk fatty acids in dairy cattle. J. Anim. Sci. Biotechnol. 2021, 12, 24. [Google Scholar] [CrossRef]
  95. Littlejohn, M.D.; Tiplady, K.; Lopdell, T.; Law, T.A.; Scott, A.; Harland, C.; Sherlock, R.; Henty, K.; Obolonkin, V. Expression Variants of the Lipogenic AGPAT6 Gene Affect Diverse Milk Composition Phenotypes in Bos taurus. PLoS ONE 2014, 9, e85757. [Google Scholar] [CrossRef]
  96. Lu, A.; Xiucai, H.; Hong, C.; Jihong, J.; Chunlei, Z.; Haixia, X.; Xueyuan, G. Single nucleotide polymorphisms in bovine PRL gene and their associations with milk production traits in Chinese Holsteins. Mol. Biol. Rep. 2010, 37, 547–551. [Google Scholar] [CrossRef]
  97. Dong, H.C.; Song, X.M.; Zhang, L.; Jiang, J.F.; Zhou, J.P.; Jiang, Y.Q. New insights into the prolactin-RsaI (PRL-RsaI) locus in Chinese Holstein cows and its effect on milk performance traits. Genet. Mol. Res. 2013, 12, 5766–5773. [Google Scholar] [CrossRef]
  98. Pegolo, S.; Dadousis, C.; Mach, N.; Ramayo-Caldas, Y.; Mele, M.; Conte, G.; Schiavon, S.; Bittante, G.; Cecchinato, A. SNP co-association and network analyses identify E2F3, KDM5A and BACH2 as key regulators of the bovine milk fatty acid profile. Sci. Rep. 2017, 7, 17317. [Google Scholar] [CrossRef] [Green Version]
  99. Li, C.; Wentao, C.; Shuli, L.; Chenghao, Z.; Hongwei, Y.; Dongxiao, S.; Shengli, Z. SERPINA1 gene identified in RNA-Seq showed strong association with milk protein concentration in Chinese Holstein cows. PeerJ 2020, 8, e8460. [Google Scholar] [CrossRef] [Green Version]
  100. Li, Y.; Huitong, Z.; Long, C.; Jenny, Z.; Jonathan, H. Variation in the stearoyl-CoA desaturase gene (SCD) and its influence on milk fatty acid composition in late-lactation dairy cattle grazed on pasture. Arch. Anim. Breed 2020, 63, 355–366. [Google Scholar]
  101. Alim, A.M.; Dong, T.; Xie, Y.; Zhang, S.L.; Zhang, D.X.; Sun, Q.; Zhang, L.; Liu, L.; Guo, G. Genetic effects of stearoyl-CoA desaturase 1 (SCD1) polymorphism on milk production traits in the Chinese dairy population. Mol. Biol. Rep. 2012, 39, 8733–8740. [Google Scholar] [CrossRef]
  102. Soyeurt, H.; Gillon, A.; Vanderick, S.; Mayeres, P.; Bertozzi, C.; Gengler, N. Estimation of Heritability and Genetic Correlations for the Major Fatty Acids in Bovine Milk. J. Dairy Sci. 2007, 90, 4435–4442. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  103. He, Y.; Chu, Q.; Ma, P.; Wang, Y.; Zhang, Q.; Sun, D.; Zhang, Y.; Yu, Y.; Zhang, Y. Association of bovine CD4 and STAT5b single nucleotide polymorphisms with somatic cell scores and milk production traits in Chinese Holsteins. J. Dairy Res. 2011, 78, 242–249. [Google Scholar] [CrossRef]
  104. Wang, H.; Jiang, L.; Liu, X.; Yang, J.; Wei, J. A Post-GWAS Replication Study Confirming the PTK2 Gene Associated with Milk Production Traits in Chinese Holstein. PLoS ONE 2013, 8, e83625. [Google Scholar] [CrossRef] [PubMed]
  105. Xu, Q.; Gui, M.; Sun, D.; Qin, Z.; Yuan, Z.; Yin, C.; Chen, H.; Ding, X.; Liu, J. Detection of genetic association and functional polymorphisms of UGDH affecting milk production trait in Chinese Holstein cattle. BMC Genom. 2012, 13, 590. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  106. Yang, S.; Gao, Y.; Zhang, S.; Zhang, Q.; Sun, D. Identification of Genetic Associations and Functional Polymorphisms of SAA1 Gene Affecting Milk Production Traits in Dairy Cattle. PLoS ONE 2016, 11, e0162195. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  107. Yang, S.; Li, C.; Xie, Y.; Cui, X.; Li, X.; Wei, J.; Zhang, Y.; Yu, Y.; Wang, Y.; Zhang, S.; et al. Detection of functional polymorphisms influencing the promoter activity of theSAA2gene and their association with milk production traits in Chinese Holstein cows. Anim. Genet. 2015, 46, 591–598. [Google Scholar] [CrossRef] [PubMed]
  108. Han, B.; Liang, W.; Liu, L.; Li, Y.; Sun, D. Genetic association of the ACACB gene with milk yield and composition traits in dairy cattle. Anim. Genet. 2018, 49, 169–177. [Google Scholar] [CrossRef] [PubMed]
  109. Li, C.; Miao, W.; Wentao, C.; Shuli, L.; Chenghao, Z.; Hongwei, Y.; Dongxiao, S.; Shengli, Z. Genetic Analyses Confirm SNPs in HSPA8 and ERBB2 are Associated with Milk Protein Concentration in Chinese Holstein Cattle. Genes 2019, 10, 104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  110. Shi, L.; Liu, L.; Ma, Z.; Lv, X.; Li, C.; Xu, L.; Han, B.; Li, Y.; Zhao, F.; Yang, Y.; et al. Identification of genetic associations of ECHS 1 gene with milk fatty acid traits in dairy cattle. Anim. Genet. 2019, 50, 430–438. [Google Scholar] [CrossRef]
  111. Li, M.; Qisong, G.; Mengqi, W.; Yan, L.; Yujia, S.; Zhi, C.; Huimin, Z.; Niel, A.K.; Zhangping, Y.; Yongjiang, M. Polymorphisms in Fatty Acid Desaturase 2 Gene Are Associated with Milk Production Traits in Chinese Holstein Cows. Animals 2020, 10, 671. [Google Scholar] [CrossRef] [Green Version]
  112. Polasik, D.; Goli’nczak, J.; Proskura, W.; Terman, A.; Dybus, A. Association between THRSP Gene Polymorphism and Fatty Acid Composition in Milk of Dairy Cows. Animals 2021, 11, 1144. [Google Scholar] [CrossRef]
  113. Rincón, G.; Islas-Trejo, A.; Castillo, A.R.; Bauman, D.E.; German, B.J.; Medrano, J.F. Polymorphisms in genes in the SREBP1 signalling pathway and SCD are associated with milk fatty acid composition in Holstein cattle. J. Dairy Res. 2012, 79, 66–75. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  114. Bhattarai, D.; Chen, X.; Rehman, Z.U.; Hao, X.; Ullah, F.; Dad, R.; Talpur, H.S.; Kadariya, I.; Cui, L.; Fan, M.; et al. Association of MAP4K4 gene single nucleotide polymorphism with mastitis and milk traits in Chinese Holstein cattle. J. Dairy Res. 2017, 84, 76–79. [Google Scholar] [CrossRef] [PubMed]
  115. Dux, M.; Magdalena, M.; Eulalia, S.; Dagmara, R.; Krzysztof, F.; Emilia, B.; Lech, Z. Association of SNP and STR polymor-phisms of insulin-like growth factor 2 receptor (IGF2R) gene with milk traits in Holstein-Friesian cows. J. Dairy Res. 2018, 85, 138–141. [Google Scholar] [CrossRef]
  116. Selvaggi, M.; Claudia, C.; Francesca, C.; Sara, A.; Giulio, A.; Tufarelli, V.; Dario, C. Determination of a possible relationship between a single nucleotide polymorphism (SNP) in the promoter region of the SIRT1 gene and production and reproduction traits in the Ag-erolese cattle breed. Arch. Anim. Breed. 2019, 62, 107–112. [Google Scholar] [CrossRef] [PubMed]
  117. Li, C.; Sun, D.; Zhang, S.; Yang, S.; Alim, M.A.; Zhang, Q.; Li, Y.; Liu, L. 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. BMC Genet. 2016, 17, 110. [Google Scholar] [CrossRef] [Green Version]
  118. Soltani-Ghombavani, M.; Ansari-Mahyari, S.; Rostami, M.; Ghanbari-Baghenoei, S.; Edriss, M. Effect of polymorphisms in the ABCG2, LEPR and SCD1 genes on milk production traits in Holstein cows. S. Afr. J. Anim. Sci. 2016, 46, 196. [Google Scholar] [CrossRef] [Green Version]
  119. Tumino, S.; Criscione, A.; Moltisanti, V.; Marletta, D.; Bordonaro, S.; Avondo, M.; Valenti, B. Feeding System Resizes the Effects of DGAT1 Polymorphism on Milk Traits and Fatty Acids Composition in Modicana Cows. Animals 2021, 11, 1616. [Google Scholar] [CrossRef]
  120. Shi, L.; Lin, L.; Xiaoqing, L.; Zhu, M.; Yuze, Y.; Yanhua, L.; Feng, Z.; Sun, D.; Han, B. Polymorphisms and genetic effects of PRLR, MOGAT1, MINPP1 and CHUK genes on milk fatty acid traits in Chinese Holstein. BMC Genet. 2019, 20, 69. [Google Scholar] [CrossRef] [PubMed]
  121. Shi, L.; Liu, L.; Lv, X.; Ma, Z.; Li, C.; Li, Y.; Zhao, F.; Sun, D.; Han, B. Identification of genetic effects and potential causal polymorphisms of CPM gene impacting milk fatty acid traits in Chinese Holstein. Anim. Genet. 2020, 51, 491–501. [Google Scholar] [CrossRef]
  122. Shi, L.; Han, B.; Lin, L.; Xiaoqing, L.; Zhu, M.; Li, C.; Lingna, X.; Li, Y.; Zhao, F.; Yang, Y.; et al. Determination of Genetic Effects of LIPK and LIPJ Genes on Milk Fatty Acids in Dairy Cattle. Genes 2019, 10, 86. [Google Scholar] [CrossRef] [Green Version]
  123. Argov-Argaman, N.; Mida, K.; Cohen, B.-C.; Visker, M.; Hettinga, K. Milk Fat Content and DGAT1 Genotype Determine Lipid Composition of the Milk Fat Globule Membrane. PLoS ONE 2013, 8, e68707. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  124. Barbosa, D.S.; Sonstegars, M.; Thallman, S.; Connor, R.; Schnabel, E.; Van Tassell, C.P. Characterization of DGAT1 allelic effects in a sample of North American Holstein cattle. Anim. Biotechnol. 2010, 21, 88–99. [Google Scholar] [CrossRef] [PubMed]
  125. Berry, D.P.; Howard, D.; O’Boyle, P.; Waters, S.; Kearney, J.; McCabe, M. Associations between the K232A polymorphism in the diacylglycerol-O-transferase 1 (DGAT1) gene and performance in Irish Holstein-Friesian dairy cattle. Ir. J. Agric. Food Res. 2010, 49, 1–9. [Google Scholar]
  126. Hill, R.; Canal, A.; Bondioli, K.; Morell, R.; Garcia, M. Molecular markers located on the DGAT1, CAST, and LEPR genes and their associations with milk production and fertility traits in Holstein cattle. Genet. Mol. Res. 2016, 15. [Google Scholar] [CrossRef]
  127. Bobbo, T.; Tiezzi, F.; Penasa, M.; De Marchi, M.; Cassandro, M. Short communication: Association analysis of diacylglycerol acyltransferase (DGAT1) mutation on chromosome 14 for milk yield and composition traits, somatic cell score, and coagulation properties in Holstein bulls. J. Dairy Sci. 2018, 101, 8087–8091. [Google Scholar] [CrossRef]
  128. Bovenhuis, H.; Visker, W.; van Valenberg, H.J.F.; Buitenhuis, A.J.; van Arendonk, J.A.M. Effects of the DGAT1 polymorphism on test-day milk production traits throughout lactation. J. Dairy Sci. 2015, 98, 6572–6582. [Google Scholar] [CrossRef] [Green Version]
  129. Kadlecova, V.; Nemeckova, D.; Jecminkova, K.; Stadnik, L. The effects of polymorphism in the dgat1 gene on energy balance and milk production traits in primiparous Holstein cows during the first six months of lactation. Bulg. J. Agric. Sci. 2014, 20, 203–209. [Google Scholar]
  130. Komisarek, J.; Kolenda, M. The effect of DGAT1 polymorphism on milk production traits in dairy cows depending on envi-ronmental temperature. Turk. J. Vet. Anim. Sci. 2016, 40, 251–254. [Google Scholar] [CrossRef]
  131. Lu, J.; Boeren, S.; van Hooijdonk, T.; Vervoort, J.; Hettinga, K. Effect of the DGAT1 K232A genotype of dairy cows on the milk metabolome and proteome. J. Dairy Sci. 2015, 98, 3460–3469. [Google Scholar] [CrossRef] [Green Version]
  132. Pacheco-Pappenheim, S.; Yener, S.; Van Valenberg, H.J.; Tzompa-Sosa, D.A.; Bovenhuis, H. The DGAT1 K232A polymorphism and feeding modify milk fat triacylglycerol composition. J. Dairy Sci. 2019, 102, 6842–6852. [Google Scholar] [CrossRef] [PubMed]
  133. Streit, M.; Neugebauer, N.; Meuwissen, T.; Bennewitz, J. Short communication: Evidence for a major gene by polygene interaction for milk production traits in German Holstein dairy cattle. J. Dairy Sci. 2011, 94, 1597–1600. [Google Scholar] [CrossRef] [Green Version]
  134. Szyda, J.; Morek-Kopeć, M.; Komisarek, J.; Żarnecki, A. Evaluating markers in selected genes for association with functional longevity of dairy cattle. BMC Genet. 2011, 12, 30. [Google Scholar] [CrossRef] [Green Version]
  135. Tabaran, A.-F.; Balteanu, V.A.; Gal, E.; Pusta, D.; Mihaiu, R.; Dan, S.D.; Mihaiu, M. Influence of DGAT1 K232A Polymorphism on Milk Fat Percentage and Fatty Acid Profiles in Romanian Holstein Cattle. Anim. Biotechnol. 2014, 26, 105–111. [Google Scholar] [CrossRef]
  136. Tzompa-Sosa, D.A.; Van Valenberg, H.; Van Aken, G.; Bovenhuis, H. Milk fat triacylglycerols and their relations with milk fatty acid composition, DGAT1 K232A polymorphism, and milk production traits. J. Dairy Sci. 2016, 99, 3624–3631. [Google Scholar] [CrossRef] [PubMed]
  137. van Gastelen, S.; Antunes-Fernandes, E.C.; Hettinga, K.A.; Dijkstra, J. The effect of linseed oil and DGAT1 K232A polymor-phism on the methane emission prediction potential of milk fatty acids. J. Dairy Sci. 2018, 101, 5599–5604. [Google Scholar] [CrossRef] [Green Version]
  138. Krovvidi, S.; Thiruvenkadan, A.K.; Murali, N.; Saravanan, R.; Vinoo, R.; Metta, M. Evaluation of non-synonym mutation in DGAT1 K232A as a marker for milk production traits in Ongole cattle and Murrah buffalo from Southern India. Trop. Anim. Health Prod. 2021, 53, 118. [Google Scholar] [CrossRef]
  139. Szewczuk, M. Association of a genetic marker at the bovine Janus kinase 2 locus (JAK2/RsaI) with milk production traits of four cattle breeds. J. Dairy Res. 2015, 82, 287–292. [Google Scholar] [CrossRef] [PubMed]
  140. Liang, Y.; Qisong, G.; Qiang, Z.; Abdelaziz, A.; Li, M.; Zhangping, Y.; Niel, K.; Mao, Y. Polymorphisms of the ACSL1 Gene Influence Milk Production Traits and Somatic Cell Score in Chinese Holstein Cows. Animals 2020, 10, 2282. [Google Scholar] [CrossRef]
  141. Qin, Y.; Chen, S.; Song-Jia, L. Polymorphisms of Mitochondrial ATPASE 8/6 Genes and Association with Milk Production Traits in Holstein Cows. Anim. Biotechnol. 2012, 23, 204–212. [Google Scholar] [CrossRef]
  142. Derakhshani, H.; Jan, P.; Jeroen, B.; Herman, B.; Ehsan, K.W. Association of bovine major histocompatibility complex (BoLA) gene polymorphism with colostrum and milk microbiota of dairy cows during the first week of lactation. Microbiome 2018, 6, 203. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  143. Rahmatalla, A.S.; Müller, U.; Strucken, E.M.; Reissmann, M.; Gudrun, A. BrockmannThe F279Y polymorphism of the GHR gene and its relation to milk production and somatic cell score in German Holstein dairy cattle. J. Appl. Genet. 2011, 52, 459–465. [Google Scholar] [CrossRef]
  144. Sajjanar, B.; Deb, R.; Singh, U.; Kumar, S.; Brahmane, M.; Nirmale, A.; Bal, S.K.; Minhas, P.S. Identification of SNP inHSP90AB1and its Association with the Relative Thermotolerance and Milk Production Traits in Indian Dairy Cattle. Anim. Biotechnol. 2014, 26, 45–50. [Google Scholar] [CrossRef] [PubMed]
  145. Cao, M.; Shi, L.; Peng, P.; Han, B.; Liu, L.; Lv, X.; Ma, Z.; Zhang, S.; Sun, D. Determination of genetic effects and functional SNPs of bovine HTR1B gene on milk fatty acid traits. BMC Genom. 2021, 22, 575. [Google Scholar] [CrossRef]
  146. Bagnicka, E.; Siadkowska, E.; Strzałkowska, N.; Żelazowska, B.; Flisikowski, K.; Krzyżewski, J.; Zwierzchowski, L. Association of polymorphisms in exons 2 and 10 of the insulin-like growth factor 2 (IGF2) gene with milk production traits in Polish Holstein-Friesian cattle. J. Dairy Res. 2010, 77, 37–42. [Google Scholar] [CrossRef] [PubMed]
  147. Berkowicz, E.W.; Magee, D.A.; Sikora, K.M.; Berry, D.P.; Howard, D.J.; Mullen, M.P.; Evans, R.D.; Spillane, C.; Machugh, D.E. Single nucleotide polymorphisms at the imprinted bovine insulin-like growth factor 2 (IGF2) locus are associated with dairy performance in Irish Holstein-Friesian cattle. J. Dairy Res. 2010, 78, 1–8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  148. Clempson, M.M.; Pollott, G.E.; Brickell, J.; Wathes, D.C. Associations Between Bovine IGFBP2 Polymorphisms with Fertility, Milk Production, and Metabolic Status in UK Dairy Cows. Anim. Biotechnol. 2012, 23, 101–113. [Google Scholar] [CrossRef]
  149. Yang, S.-H.; Bi, X.-J.; Xie, Y.; Li, C.; Zhang, S.-L.; Zhang, Q.; Sun, D.-X. Validation of PDE9A Gene Identified in GWAS Showing Strong Association with Milk Production Traits in Chinese Holstein. Int. J. Mol. Sci. 2015, 16, 26530–26542. [Google Scholar] [CrossRef] [Green Version]
  150. Zheng, X.; Ju, Z.; Wang, J.; Li, Q.; Huang, J.; Zhang, A.; Zhong, J.; Wang, C. Single nucleotide polymorphisms, haplotypes and combined gen-otypes of LAP3 gene in bovine and their association with milk production traits. Mol. Biol. Rep. 2011, 38, 4053–4061. [Google Scholar] [CrossRef]
  151. Han, B.; Yuan, Y.; Liang, R.; Li, Y.; Liu, L.; Sun, D. Genetic Effects of LPIN1 Polymorphisms on Milk Production Traits in Dairy Cattle. Genes 2019, 10, 265. [Google Scholar] [CrossRef] [Green Version]
  152. Haruna, I.L.; Zhou, H.; Hickford, J.G.H. Variation in bovine leptin gene affects milk fatty acid composition in New Zealand Holstein Friesian × Jersey dairy cows. Arch. Anim. Breed. 2021, 64, 245–256. [Google Scholar] [CrossRef]
  153. Wang, C.; Liu, M.; Li, Q.; Ju, Z.; Huang, J.; Li, J.; Wang, H.; Zhong, J. Three novel single-nucleotide polymorphisms of MBL1 gene in Chinese native cattle and their associations with milk performance traits. Vet. Immunol. Immunopathol. 2011, 139, 229–236. [Google Scholar] [CrossRef] [PubMed]
  154. Zhao, L.Z.; Wang, C.F.; Li, Q.L.; Ju, Z.H.; Huang, J.M.; Li, J.B.; Zhong, J.F.; Zhang, J.B. Novel SNPs of the mannan-binding lectin 2 gene and their association with production traits in Chinese Holsteins. Genet. Mol. Res. 2012, 11, 3744–3754. [Google Scholar] [CrossRef] [PubMed]
  155. Wang, X.; Peñagaricano, F.; Tal-Stein, R.; Lipkin, E.; Khatib, H. Short communication: Association of an OLR1 polymorphism with milk production traits in the Israeli Holstein population. J. Dairy Sci. 2012, 95, 1565–1567. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  156. Shi, L.; Xiaoqing, L.; Lin, L.; Yuze, Y.; Zhu, M.; Han, B.; Sun, D. A post-GWAS confirming effects of PRKG1 gene on milk fatty acids in a Chinese Holstein dairy population. BMC Genet. 2019, 20, 53. [Google Scholar] [CrossRef] [Green Version]
  157. Nafikov, R.A.; Schoonmaker, J.P.; Korn, K.T.; Noack, K.; Garrick, D.J.; Koehler, K.J.; Minick-Bormann, J.; Reecy, J.M.; Spurlock, D.E.; Beitz, D.C. Sterol regulatory element binding transcription factor 1 (SREBF1) polymorphism and milk fatty acid composition. J. Dairy Sci. 2013, 96, 2605–2616. [Google Scholar] [CrossRef]
  158. Lv, Y.; Wei, C.; Zhang, L.; Lu, G.; Liu, K.; Du, L. Association Between Polymorphisms in the SLC27A1 Gene and Milk Pro-duction Traits in Chinese Holstein Cattle. Anim. Biotechnol. 2011, 22, 1–6. [Google Scholar] [CrossRef]
  159. Nafikov, A.R.; Schoonmaker, P.J.; Korn, T.K.; Noack, K.; Garrick, J.D.; Koehler, K.J.; Minick-Bormann, J.; Reecy, J.M.; Spurlock, D.E.; Beitz, D.C. Association of polymorphisms in solute carrier family 27, isoform A6(SLC27A6) and fatty acid-binding protein-3 and fatty acid-binding protein-4 (FABP3 and FABP4) with fatty acid composition of bovine milk. J. Dairy Sci. 2013, 96, 6007–6021. [Google Scholar] [CrossRef] [Green Version]
  160. Wang, M.; Song, H.; Zhu, X.; Xing, S.; Zhang, M.; Zhang, H.; Wang, X.; Yang, Z.; Ding, X.; Karrow, N.A.; et al. Toll-like receptor 4 gene polymorphisms influence milk production traits in Chinese Holstein cows. J. Dairy Res. 2018, 85, 407–411. [Google Scholar] [CrossRef]
  161. Zhou, H.; Cheng, L.; Gong, H.; Byun, S.O.; Edwards, G.R.; Hickford, J.G.H. Variation in the Toll-like Receptor 4 (TLR4) gene affects milk traits in dairy cows. J. Dairy Res. 2017, 84, 426–429. [Google Scholar] [CrossRef]
  162. Mao, Y.; Zhu, X.; Xing, S.; Zhang, M.; Zhang, H.; Wang, X.; Karrow, N.; Yang, L.; Yang, Z. Polymorphisms in the promoter region of the bovine lactoferrin gene influence milk somatic cell score and milk production traits in Chinese Holstein cows. Res. Vet. Sci. 2015, 103, 107–112. [Google Scholar] [CrossRef] [PubMed]
  163. Ju, Z.; Li, Q.; Huang, J.; Hou, M.; Li, R.; Li, J.; Zhong, J.; Wang, C. Three novel SNPs of the bovine Tf gene in Chinese native cattle and their associations with milk production traits. Genet. Mol. Res. 2011, 10, 340–352. [Google Scholar] [CrossRef] [PubMed]
  164. Zwierzchowski, L.; Ostrowska, M.; Żelazowska, B.; Bagnicka, E. Single nucleotide polymorphisms in the bovine SLC2A12 and SLC5A1 glucose transporter genes—The effect on gene expression and milk traits of Holstein Friesian cows. Anim. Biotechnol. 2021, 1–11. [Google Scholar] [CrossRef] [PubMed]
  165. Dar, M.R.; Singh, M.; Thakur, S.; Verma, A. Exploring the relationship between polymorphisms of leptin and IGF-1 genes with milk yield in indicine and taurine crossbred cows. Trop. Anim. Health Prod. 2021, 53, 1–8. [Google Scholar] [CrossRef]
  166. Ostrowska, M.; Zwierzchowski, L.; Brzozowska, P.; Kawecka-Grochocka, E.; Żelazowska, B.; Bagnicka, E. The effect of single nucleotide polymorphism in the promoter region of bovine alpha-lactalbumin (LALBA) gene on LALBA expression in milk cells and milk traits of cows. J. Anim. Sci. 2021, 99, skab169. [Google Scholar] [CrossRef] [PubMed]
  167. Bordonaro, S.; Tumino, S.; Donata, M.; Anna, D.; Fortunato, D.P.; Avondo, M.; Valenti, B. Effect of GH p.L127V Polymorphism and Feeding Systems on Milk Production Traits and Fatty Acid Composition in Modicana Cows. Animals 2020, 10, 1651. [Google Scholar] [CrossRef]
  168. Friedrich, J.; Brand, B.; Ponsuksili, S.; Graunke, K.L.; Langbein, J.; Knaust, J.; Kühn, C.; Schwerin, M. Detection of genetic variants affecting cattle behaviour and their impact on milk production: A genome-wide association study. Anim. Genet. 2016, 47, 12–18. [Google Scholar] [CrossRef]
  169. Wang, H.; Jiang, L.; Wang, W.; Zhang, S.; Yin, Z.; Zhang, Q.; Liu, J.-F. Associations between variants of the HALgene and milk production traits in Chinese Holstein cows. BMC Genet. 2014, 15, 125. [Google Scholar] [CrossRef] [Green Version]
  170. Li, C.; Cai, W.; Liu, S.; Zhou, C.; Cao, M.; Yin, H.; Sun, D.; Zhang, S.; Loor, J.J. Association of UDP-galactose-4-epimerase with milk protein concentration in the Chinese Holstein population. Asian-Australas. J. Anim. Sci. 2020, 33, 1725–1731. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ma, Y.; Khan, M.Z.; Xiao, J.; Alugongo, G.M.; Chen, X.; Chen, T.; Liu, S.; He, Z.; Wang, J.; Shah, M.K.; et al. Genetic Markers Associated with Milk Production Traits in Dairy Cattle. Agriculture 2021, 11, 1018. https://0-doi-org.brum.beds.ac.uk/10.3390/agriculture11101018

AMA Style

Ma Y, Khan MZ, Xiao J, Alugongo GM, Chen X, Chen T, Liu S, He Z, Wang J, Shah MK, et al. Genetic Markers Associated with Milk Production Traits in Dairy Cattle. Agriculture. 2021; 11(10):1018. https://0-doi-org.brum.beds.ac.uk/10.3390/agriculture11101018

Chicago/Turabian Style

Ma, Yulin, Muhammad Zahoor Khan, Jianxin Xiao, Gibson Maswayi Alugongo, Xu Chen, Tianyu Chen, Shuai Liu, Zhiyuan He, Jingjun Wang, Muhammad Kamal Shah, and et al. 2021. "Genetic Markers Associated with Milk Production Traits in Dairy Cattle" Agriculture 11, no. 10: 1018. https://0-doi-org.brum.beds.ac.uk/10.3390/agriculture11101018

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop