Milk fatty acid variability and association with polymorphisms in SCD1 and DGAT1 genes in White Fulani and Borgou cattle breeds

The stearoyl-CoA desaturase 1 (SCD1) A293V and acyl CoA: diacylglycerol acyltransferase 1 (DGAT1) K232A polymorphisms have been associated with significant variation in bovine milk fatty acid composition and unsaturation indices in western cattle breeds. This study aimed to estimate the milk fatty acid variability in indigenous Borgou and White Fulani cattle breeds of Benin, and the effects of the SCD1 A293V and DGAT1 K232A polymorphisms on milk and fatty acid composition and unsaturation indices. Thus, 85 Borgou and 96 White Fulani cows were genotyped for the SCD1 A293V and DGAT1 K232A polymorphisms and their milk and fatty acid composition and unsaturation indices were determined. Borgou presented milk with higher linoleic acid (P < 0.001), oleic acid (P < 0.05), C18 index (P < 0.001), total unsaturation index (P < 0.05), and lower total saturated fatty acid (SFA) compared to White Fulani. The SCD1 VV genotype was associated with higher protein and lactose contents in White Fulani (P < 0.05). In Borgou, the SCD1 AV genotype was associated with higher C14 and total unsaturation indices (P < 0.01), while the SCD1 V allele was associated with decrease in C14 index (P < 0.05). In White Fulani, the SCD1 VV genotype was associated with lower C18:1 cis-9 content (P < 0.05) while the DGAT1 K allele was associated with increased total SFA (P < 0.05), and decreased C18 index (P < 0.05), total unsaturation index (P < 0.01) and total monounsaturated fatty acid (P < 0.01). The SCD1 A293V and DGAT1 K232A may serve as genetic markers to improve milk fatty acid traits in Borgou and White Fulani breeds. Electronic supplementary material The online version of this article (10.1007/s11033-018-4331-4) contains supplementary material, which is available to authorized users.


Introduction
Milk and dairy products are important sources of energy, fat, high quality protein, vitamins, and minerals in human diets [1]. The fatty acids fraction of bovine milk fat is characterized by high amount (50-70%) of saturated fatty acids (SFA), 20-40% monounsaturated (MUFA) and low amount (1-5%) of polyunsaturated fatty acids (PUFA) [2]. Studies show that high quantities of dietary SFAs are associated with an increase in blood cholesterol and, therefore, with increased risk of atherosclerosis and coronary heart diseases in humans [3][4][5]. On the contrary, high concentrations of PUFAs in blood and tissue lipids is associated with positive benefits on cardiovascular diseases, mental health [6], exert anticancer effects [7] and play important role in preventing and/or treating various immunes disorders such as allergies [8,9].
With regards to human health aspects, increasing unsaturated fatty acids in cow's milk is an important selection goal [10]. The cow milk fatty acids are derived from two sources, the diet and the microbial activity in the rumen [11]. The fatty acids originating from the blood or from de novo fatty acid synthesis can be desaturated in the mammary gland [10]. Thus, the degree of unsaturation of milk fat is designated by a so-called unsaturation index, which is the ratio of the unsaturated product to the sum of the unsaturated product and the saturated substrate [10,12]. Several studies showed that the variability of fatty acid composition in cow milk is affected by diet, breed, genetics, parity and stage of lactation [12][13][14]. In ruminant animals, milk fatty acid synthesis is regulated by several important enzymes including stearoyl CoA desaturase 1 (SCD1), acetyl CoA carboxylase, acyltransferases and fatty acid synthase [15]. Polymorphisms in diacylglycerol acyltransferase 1 (DGAT1) and SCD1 genes coding for key enzymes in mammary lipid metabolism have been associated with significant proportion of variation in milk fatty acid composition and unsaturation indices in different cattle populations [10,[15][16][17][18][19][20][21][22].
The Stearoyl-CoA Desaturase 1 (SCD1) gene is mapped on chromosome 26 in cattle and expressed in a variety of tissues including adipose and mammary tissue [28]. Moreover, the SCD1 enzyme contributes to the desaturation of SFAs into delta-9 unsaturated fatty acids [29]. A single nucleotide polymorphism SNP in exon 5 (C878T) of the SCD1 gene leads to valine substitution by alanine at amino acid position 293 in the mature protein (Ala293Val). The SCD1 A (Ala293) allele was associated with higher cis-9 C18:1 and total monounsaturated content as well as C14:1/C14 ratio in Italian Holstein [22]. However, data on the genetic variability of milk fatty acid composition and unsaturation indices are scarce in African indigenous cattle breeds and the few studies on DGAT1 K232A polymorphism in African indigenous cattle did not show association with milk fatty acid composition [30][31][32]. In the study of Rahmatalla et al. [31], the DGAT1 K232A genotypes showed significant effects on fat content in Sudanese indigenous Kenana and Butana cattle breeds. The study of Houaga et al. [32] mentioned that the DGAT1 KK genotype was significantly associated with higher fat yield in White Fulani (P < 0.05). However, they did not investigate the milk fatty acid composition.
In Benin, indigenous White Fulani and Borgou cattle breeds are the main milk producers [32]. However, to date, there has been no study on SCD1 A293V and DGAT1 K232A influence on milk and fatty acid composition and unsaturation indices. Moreover, no detailed data on milk fatty acid composition of these breeds is available. Such data would be useful for gathering knowledge on nutritional value of indigenous cow milk for the dairy industry in Benin and for the opportunity to improve cow's milk fatty acid composition with regards to human health. Therefore, the present study aimed to estimate the milk fatty acid variation in white Fulani and Borgou cattle breeds and the effects of the DGAT1 K232A and SCD1 A293V polymorphisms on milk and fatty acid composition and unsaturation indices.

Sampling
A total of 85 Borgou and 96 White Fulani indigenous cows were sampled from state owned farms (Betecoucou, Okpara and Samiondji) and privately owned farms in Benin between May-July 2016 (raining season). Blood and milk samples were obtained from three to ten cows per herd from a total of 17 herds. Only cows in lactation were randomly selected and sampled. The map of Benin indicating the sampling sites is presented in Fig. 1. The feeding system was based solely on natural grazing without concentrate supplementation. The cows were milked once a day in the morning. Additional meta data on sampled cows obtained from the livestock keepers and herders following a written consent permitting sampling included: age, lactation stage and parity number. Blood samples were collected from the jugular vein into 10 ml EDTA vacutainer tubes and immediately transported to the laboratory in a cool box containing ice and stored at − 20 °C until further analysis. Milk samples were aseptically collected into 50 ml falcon tubes each containing one tablet of Bonopol milk preservative (Systems Plus, Canada) and sent to Valacta laboratories (Valacta Laboratories Inc., Canada, http://www.valac ta.com) for the analysis of milk components.

Milk component analysis
Test-day milk fat percent, protein percent, milk urea nitrogen and lactose percent were determined in milk samples with MilkoScan FT 6000 Series mid-range infrared Fourier transform infrared-based spectrometers (Foss, Hillerod, Denmark) by Valacta Laboratories.

Fatty acids analysis and quantification
The milk samples were prepared for analysis using alkali hydrolysis followed by methylation as described in Association of Analytical Communities (AOAC) Method 991.39 [33]. Briefly, 1 ml of cow milk was transferred into a 15 ml screwcap centrifuge tube in triplicates and 3 ml of 0.5N NaOH in methanol was added into each tube and mixed thoroughly by shaking for 30 s. Tubes were incubated in a water bath at 85 °C for 10 min and cooled at ambient temperature. In a fume hood, 1 ml of boron trifluoride (BF 3 ) was added, mixed thoroughly by shaking and incubated in a water bath at 85 °C for 10 min. Samples were allowed to cool to room temperature followed by addition of 3 ml iso-octane and 3 ml of saturated NaCl solution, mixed by shaking vigorously followed by centrifugation at 2500 rpm for 5 min. The upper layer (iso-octane) containing the fatty acid methyl esters (FAMEs) was transferred through a funnel (cotton wool was placed in the funnel and 1 g of sodium sulfate (Na 2 SO 4 ) anhydrous was added on top of the cotton wool) into a test tube with a Pasteur pipette. 2 g of anhydrous Na 2 SO 4 was added. The FAMEs were diluted to a volume of 2 ml with hexane in a volumetric flask. 1 ml of the sample was transferred into a screw cap gas chromatography (GC) sample vial and stored at − 20 °C.
The composition of individual FAMEs was analyzed by gas chromatography mass spectrometry (GC-MS). The FAMEs in hexane (1 µl) were injected in to a 7890A GC system (Agilent Technologies, USA) coupled to a 240 ion trap mass spectrometer detector (Agilent Technologies) using the Agilent 7693A automatic liquid sampler at a split ratio of 100:1. A VF5-MS (5% phenyl methylpolysiloxane), 30 m × 0.25 mm id, 0.25 µm film capillary column was used with the injector port set at 280 °C. Helium was used as carrier gas at a flow rate of 1 ml/min. The oven temperature was programmed to rise from 50 °C to 180 °C at 4 °C/min followed by an increase to 250 °C at 3 °C/min. The ion trap mass spectrometer parameters were as follow: scan range 50-540 (m/z), ionization mode EI, filament delay time 3 min and transfer line temperature, manifold temperature and trap temperature of 250 °C, 100 °C and 150 °C, respectively.
Chromatograms and spectra representing individual FAMEs were analyzed using the automated mass spectral deconvolution and identification system software (AMDIS, US). The identification of the individual FAMES was performed by comparing each of the mass spectra with the database of NIST 11 (Gaithersburg, MD, USA) and Wiley 7N (John Wiley, NY, USA) and also by comparing the calculated Kovats linear retention indices using retention times of n-alkane series against the values obtained in the NIST webbook for the same capillary column stationery phase [34]. The quantification of individual FAMEs was performed by the peak area percentage method. The fatty acid concentrations were expressed as the ratio of each individual fatty acid to the total of all fatty acids detected in the sample. The fatty acids unsaturation indices were calculated as the ratio of cis-9 unsaturated to cis-9 unsaturated + saturated for specific fatty acid pairs and multiplied by 100 [12]. We calculated the following indices: C14 index = C14:1 cis-9/(C14:1 cis-9 + C14:0) × 100 and C18 index = C18:1 cis-9/(C18:1 cis-9 + C18:0) × 100. The total index was calculated as total index = (C14:1 cis-9 + C18:1 cis-9)/(C14:1 cis-9 + C14:0 + C18:1 cis-9 + C18:0) × 100 [22].

Genotyping of DGAT1 K232A and SCD1 A293V mutations
Genomic DNA was isolated from blood samples using the phenol-chloroform method [35]. DNA quality was checked on 0.8% agarose gel and the quantity was checked using the NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies, Inc., USA). The DGAT1 K232A and SCD1 A293V genotypes were determined by the method of polymerase chain reaction-restriction fragment length polymorphisms (PCR-RFLP) [20,30].
The primers 5′-GCA CCA TCC TCT TCC TCA AG-3′ (forward) and 5′-GGA AGC GCT TTC GGATG-3′ (reverse) [30] were used to amplify a 411 bp fragment of the bovine DGAT1 gene containing the lysine/alanine substitution (exon 8) while 5′-CCC ATT CGC TCT TGT TCT GT-3 (forward) and 5′-CGT GGT CTT GCT GTG GAC T-3′ (reverse) [20] were used to amplify a 400 bp fragment containing the A293V polymorphism in exon 5 of the SCD1 gene. The PCR reactions were carried out in a 30 µl volume containing 45 ng of template DNA, 15 µl of PCR Master Mix (Bioneer, Korea) and 4.5 pmol of each primer (3 pmol/µl). The PCR conditions were as follows: an initial denaturation step at 94 °C for 3 min, 35 cycles of 94 °C for 45 s, 62 °C (DGAT1 K232A) or 64 °C (SCD1 A293V) for 60 s, 72 °C for 60 s, and a final extension step of 72 °C for 5 min. The PCR reactions were performed in the GeneAmp PCR System 9700 (Applied Biosystems, USA). Amplification was confirmed by running the PCR products on 1.8% agarose gel and visualized with GelDoc-It 2 Imager (Ultra-Violet Products Ltd., UK). The PCR products were purified with the QIA quick PCR Purification Kit (Qiagen, Germany).
Five microliter of purified PCR products containing the DGAT1 K232A were digested overnight at 37 °C with 10 U of EaeI restriction enzyme (New England Biolabs, Inc., USA). The digested PCR products were separated on 1.8% agarose gel stained with GelRed™ (Biotium, UK) resulting in two fragments of 203 and 208 bp (seen as a single band on gel) for AA genotype, two fragments of 203/208 bp and 411 bp for KA genotype and the undigested 411 bp for KK genotype (Additional file 1: Figure a).
Similarly, 5 µl of purified PCR products containing the SCD1 A293V mutation were digested overnight at 37 °C with 10 U of NcoI restriction enzyme (New England Biolabs, Inc., USA). The digested products were separated on 1.8% agarose gel stained with GelRed™ (Biotium, UK), visualized and scanned with GelDoc-It 2 Imager (Ultra-Violet Products Ltd., UK). The digestion patterns resulted in two fragments of 200 bp for the AA genotype, undigested 400 bp fragment for the VV genotype and 400 bp and 200 bp fragments for the AV genotypes (Additional file 1: Figure b).

Statistical analysis
The allele frequencies and test for Hardy-Weinberg Equilibrium (HWE) were performed with GENEPOP program version 1.2 [36]. To investigate the effect of breed, SCD1 A293V and DGAT1 K232A genotypes on milk fatty acids and milk components, the following mixed linear model with IBM SPSS version 20 software package was used: where Y ijklmn is the observed phenotype: fat%, protein%, lac-tose%, milk urea nitrogen (mg/dl), individual fatty acid%, groups of fatty acids (saturated fatty acids [SFA], monounsaturated fatty acids [MUFA], polyunsaturated fatty acids [PUFA]) and unsaturation indices (C14 index, C18 index and total index), µ is the population mean, Bi is the fixed effect of Breed, GDj is the fixed effect of DGAT1 genotypes (KK and KA), GS k is the fixed effect of SCD1 genotypes (VV, AV); (Bi × Gj) ij is the fixed interaction effect between breed and DGAT1 genotypes, (Bi × GS k ) ik is the fixed interaction effect between breed and SCD1 genotypes, GR is the geographical region (south and north), A m is the random animal effect and E ijklmn is the random residual error associated with each record. The age, lactation stage and parity number were included in the initial model and were dropped out from the final model due to the absence of significant effect. The allele substitution effect of DGAT1 (K variant) and SCD1 (V variant) was estimated following the method of Marchitelli et al. [37] by regressing the number of copies of DGAT1 K allele and SCD1 V allele against each of the milk fatty acid and milk component traits separately. The phenotypes were nested within breed to obtain breed specific estimates. The results of the different effects are presented as least squares means ± standard error. The Pearson correlation indices were calculated among the various milk and fatty acids and milk component traits using the IBM SPSS version 20 software package. The principal component analysis (PCA) of the significant milk traits (Fat, protein, C14:1 cis-9, C16:0, C18:2 cis-9, cis-12, C18:1 cis-9, C14 index, C18 index, Total index, SFA, MUFA and PUFA) between breeds was done using Minitab software version 18. Significance was declared at P < 0.05.
The principal component analysis (PCA) of the significant variables (Fat, protein, C14:1 cis-9, C16:0, C18:2 cis-9, cis-12, C18:1 cis-9, C14 index, C18 index, Total index, SFA, MUFA and PUFA) between breeds is presented in Fig. 2. The first component (Axis) clearly separated White Fulani from Borgou. The White Fulani population (red) was on the left side and the Borgou (blue) on the right side as shown in Fig. 2. The PCA analysis therefore showed that the two cattle populations are different for the studied milk traits.

Effect of geographical regions
The effect of geographical regions on milk fatty acid traits is presented in Table 2. Because the White Fulani cattle are only found in the northern part of Benin, only Borgou breed was considered in the geographical effect analysis. The Borgou cows from the South of Benin showed milk with higher contents of C12:0, C18:2 cis-9, cis-12 and PUFA (P < 0.05). However, the Borgou cows from the North of Benin produced milk with higher contents of C14:1 cis-9, C18:1 cis-9, MUFA, and higher C14 index, C18 index and total index compared to Borgou cows from the South of Benin (P < 0.05).

Effects of SCD1 A293V polymorphism on milk traits
The frequencies of 293V were 0.84 and 0.94 in Borgou and White Fulani respectively, and the genotypes were in Hardy-Weinberg equilibrium (Table 3). Only four Borgou cows and one White Fulani cow were of AA genotype. The AA genotype was therefore not included in the association analysis. The SCD1 VV genotype was associated with higher (P < 0.05) protein and lactose contents and lower (P < 0.05) C18:1 cis-9 content in White Fulani (Table 4). On the other hand, the SCD1 AV genotype was associated with higher (P < 0.01) C14 index and total index compared to the VV genotype in Borgou ( Table 4).
The estimates of the SCD1 293V allele substitution effect in Borgou breed are presented in Table 5. In comparison to the A allele, the V allele was associated with decrease in C14 index (− 5.68%, P < 0.05). However, no

Effect of breed and geographical regions
In the present study, we analyzed the milk and fatty acid composition as well as fatty acid unsaturation indices in indigenous White Fulani and Borgou cows in Benin. White Fulani produced milk with high (P < 0.001) fat content than Borgou. Significant differences between breeds were observed for individual fatty acids, fatty acid unsaturation indices and fatty acid groups. The milk of Borgou breed had higher MUFA and PUFA content and C18 and total indices than White Fulani. Breed effect on fatty acid composition was reported in the meat of Borgou and White Fulani cattle in Benin [38], in the milk of South African indigenous cattle breeds [39] and Italian Holstein-Friesian, Brown Swiss, Simmental and Alpine cattle breeds [40]. In this study, White Fulani presented the highest content of SFA (68.19 vs. 58.46%, P < 0.001) and the lowest content of PUFA (12.40 vs. 19.56%, P < 0.001) compared to Borgou, and similar MUFA content (19.41 vs. 21.98%, P > 0.05) for both breeds. These results corroborate previous studies in Benin indicating that the meat of White Fulani bulls had higher SFA content (49.68%) compared to Borgou (43.03%) and similar MUFA content (33.60% vs. 33.43%) [38]. In the current study, Borgou presented higher C18 index (P < 0.001) and total unsaturation index (P < 0.05) than White Fulani. The unsaturation or desaturation index of a specific fatty acid represents the ratio of the concentration of the monounsaturated product to the sum of the monounsaturated and the saturated substrate [22]. Considering human health aspects, increasing the amount of unsaturated fatty acids as well as unsaturation indices is an important selection objective [10]. Borgou presented lower SFA (58.46%) than the values of 64%, 63.7%; 60.9% and 71.9% reported in milk of free-ranging South African indigenous Boran, Nguni, Tuli and Afrikaner cattle breeds, respectively [39]. Moreover, Borgou and White Fulani presented lower MUFA content (21.98% and 19.41%, respectively) than South African indigenous cattle (MUFA content ranged from 25.7% in Afrikaner   [39]. On the contrary, Borgou and White Fulani produced milk with C18:2 cis-9, cis-12 of 15.84% and 9.85% respectively, much higher than the range of 1.3 to 1.7% observed in South African indigenous Boran, Nguni, Tuli and Afrikaner [39]. The difference in C18:2 cis-9, cis-12 content between Borgou/White Fulani and South African indigenous cattle may be due to the lower sample size in their study being 6 Boran, 9 Nguni, 10 Tuli and 6 Afrikaner [39]. Several factors influence milk fatty acid composition such as species, breed, individual variability, nutrition, stage of lactation, parity and season [40,41]. Studied animals were raised in the traditional system on natural grazing without concentrate supplementation and they were sampled at the same period eliminating the effect of season. The observed differences between breeds would therefore be due to their genetic background. However, the differences in fatty acid composition between the two breeds could also be due to the fatty acid composition of the forage consumed by the cows on natural grazing. The White Fulani and Borgou are raised in different agro-ecological zones with different floristic composition. The White Fulani cattle are found in the Northern part of Benin while Borgou cattle are found throughout the country. The forage species and variety, climate and stage of growth are important factors that affect fatty acid content and composition of forage [38,42] and can therefore affect the milk fatty acid composition of the cows. The Borgou cows from the North of Benin produced milk with higher C14:1 cis-9, C18:1 cis-9, C14 unsaturation index, C18 unsaturation index, total unsaturation index and MUFA contents than Borgou cows from the South of Benin. Consequently, milk from Borgou cows raised in the North of Benin seems to be healthier than from Borgou cows raised in the South of Benin due to its higher content of MUFA and higher unsaturation indices. The suggested favorable combination of bovine milk fatty acid composition for human health enhancement is ~ 30% SFA, 60% MUFA, and 10% PUFA [43]. In the present study, the milk fatty acid compositions were 58.48% SFA, 21.98% MUFA and 19.56 PUFA for Borgou and 68.19% SFA, 19.41% MUFA and 12.40% PUFA for White Fulani. It's clear that the current milk fatty acid composition of Borgou and White Fulani is far from optimal and there is need for modification towards an ideal profile. Linoleic acid (a PUFA) and oleic acid (a MUFA) have been associated with decreased serum total cholesterol and low-density lipoprotein cholesterol levels and reduced risk of coronary heart diseases in humans [44]. Oleic acid has anticancer and antiatherogenic properties [45]. Linoleic acid, an essential fatty acid in the omega-6 family is associated with reduced incidence of type 2 diabetes through its ability to improve sensibility to insulin [46]. On the other hand, palmitic acid (C16:0), considered as hypercholesterolemic is responsible for the increase in the concentration of low density lipoproteins (LDL) that are associated with coronary heart diseases in humans [47]. Borgou milk with its higher linoleic acid, oleic acid and lower total SFA contents as compared to White Fulani with higher total SFA and palmitic acid contents, may be preferred by consumers than White Fulani milk.

Effects of SCD1 A293V and DGAT1 K232A polymorphisms
The frequencies of SCD1 293V were 0.84 and 0.94 in Borgou and White Fulani, respectively. A higher frequency of the V allele (0.82) is also reported in Italian Brown cows [15] while a higher frequency of the A allele has been reported in Dutch Holstein-Friesian heifers (0.73), Italian Holsteins (0.57), and Canadian Jersey cows (0.80) [10,20,22]. The difference in SCD1 A293V allele frequencies between the studied indigenous breeds and western breeds can be explained by a breed specific effect. The SCD1 AV genotype was associated with higher C14 and total unsaturation index compared to the VV genotype in Borgou breed. This result did not agree with Conte et al. [15] who associated SCD1 VV genotype with higher C14 index in Italian Brown cows. Moreover, Kgwatalala et al. [20] showed that the AA genotype of SCD1 was associated with higher C14 index. The effect of SCD1 A293V genotypes on C14 index seems to vary from one breed to another. In the current study, SCD1 genotypes did not significantly affect protein or fat percentage in Borgou (Table 4) which is similar to the results of Schennink et al. [10] in Dutch Holstein-Friesian heifers. However, in the present study, the SCD1 V allele had significant negative association (− 5.68%, P < 0.05) with C14 index compared to the A allele in Borgou. The allele A of SCD1 is therefore significantly associated with 5.68% more C14 index in Borgou. The positive significant association between the allele A of SCD1 and C14:1 cis-9 and C14 index has been reported previously [10,20,22]. However, allele A of SCD1 did not show significant effect on C14:1 cis-9 in this study. This may be explained by the large sample sizes used in the other studies, namely 1725 Dutch Holstein-Friesian heifers [10], 297 Italian Holstein Friesian cows [22] and 525 Canadian Jersey cows [20]. Also, significant associations of the SCD1 A293V polymorphism with C10 index, C12 index, C16 index and C18 index has been reported [10,20]. However, no significant effect of SCD1 polymorphism was observed for C18 index in this study.
The frequencies of DGAT1 232K were 0.77 and 0.92 in Borgou and White Fulani breeds respectively. A higher frequency of the K allele in Borgou and White Fulani breeds in Benin [32] and in Sudanese Butana and Kenana cattle breeds have been reported previously [31]. However, a lower frequency of DGAT1 K allele (0.40) was reported in Dutch Holstein-Friesian heifers [10]. In this study, the DGAT1 K allele was associated with lower C18 index (P < 0.05), total unsaturation index (P < 0.01), and MUFA (P < 0.01), and with higher SFA (P < 0.05) in White Fulani breed. These results are similar to a reported by Schennink et al. [21], who studied 1762 Dutch Holstein Friesian cows and found that the DGAT1 232K allele was associated with more saturated fatty acid. However, no significant effect of DGAT1 232K allele on C18 and total unsaturation indices was found by Schennink et al. [21]. Similar to our data, Schennink et al. [10] showed that the DGAT1 K allele was associated with lower C18 and total unsaturation indices. The majority of milk fatty acids are present in the form of triacylglycerols and the DGAT1 enzyme plays an important role during the last step of triglyceride synthesis. The DGAT1 K232A polymorphism was reported to have significant association with milk fatty acid composition and unsaturation indices [15]. However, we did not observed significant association of DGAT1 K232A polymorphism with individual fatty acids but significant associations with SFA and MUFA was observed. This is conceivable because the effect of DGAT1 on fatty acid composition and saturation may be due to a higher activity and alteration of specificity of DGAT1 enzyme [21] which may vary between breeds. The discussion of our results on DGAT1 K232A polymorphism and fatty acid traits was limited to western dairy breeds because of the scarcity of data on African indigenous cattle breeds.

Phenotypic correlations
The Pearson correlation coefficient between total SFA and total PUFA was negative and moderate (− 0.44) in Borgou. Similar correlation (− 0.34) has been observed between total SFAs and total PUFAs in Canadian Holsteins [48]. The fat percentage showed positive correlations with C14:1 cis-9 (0.45) and C18:2 cis-9, cis-12 (0.24) and negative correlation with total PUFAs (− 0.25) in Borgou. Accordingly, an increase in the fat content of Borgou milk will lead to slightly higher C14:1 cis-9 and C18:2 cis-9, cis-12 contents and decreased total PUFA content. The fat content positively affect the price of milk in developed countries, therefore increasing PUFAs (decreasing fat content) should have negative economic impact [49]. However, in the African context, and in Benin in particular, the price of milk is not influenced by its fat content and hence decreasing total fat content for increased PUFAs in Borgou milk will be beneficial for human health and will not negatively affect farm incomes. On the contrary, increasing fat and protein percentage in White Fulani breed will lead to slightly higher increase in C14:1 cis-9 and C14 index due to the moderate positive correlation observed between the traits (Table 6). In the present study, C16:0, showed significant and negative correlations with fat percentage (− 0.23), C14:1 cis-9 (− 0.37), C18:2 cis-9, cis-12 (− 0.38) and positive correlation with total SFA (0.48) in Borgou breed. This implies that decreasing C16:0 will lead to increase in fat percentage, C14:1 cis-9, C18:2 cis-9, cis-12 and decrease in total SFAs in Borgou which would be an important selection goal.

Conclusion
This study has revealed significant differences in milk components, milk fatty acid composition and unsaturation indices between White Fulani and Borgou indigenous cattle breeds in Benin. The Borgou milk contained higher linoleic acid, higher oleic acid and lower total SFA compared to White Fulani which are beneficial traits for human health. The SCD1 AV genotype was associated with higher C14 and total indices; and the SCD1 V allele was associated with decrease in C14 index in Borgou. In White Fulani breed, the SCD1 VV genotype was associated with lower C18:1 cis-9 content while the DGAT1 K allele was associated with increased total SFA, and decreased C18 unsaturation index, total unsaturation index and total MUFA. SCD1 A293V and DGAT1 K232A polymorphisms may serve as potential genetic markers in a breeding program to improve milk fatty acids traits in indigenous cattle breeds in Benin. However, further studies with a large population of the Borgou and White Fulani breeds are needed to better understand the genetic variability of their milk fatty acids and association with genetic polymorphisms in SCD1 and DGAT1genes. and sample collection procedures were in accordance with the ethical standards and were approved by the Livestock Research Subdivision (LRS) ethics committee in Benin. The farmers gave a written consent prior to blood and milk sample collection.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creat iveco mmons .org/licen ses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.