Effects of DGAT1 on milk performance in Sudanese Butana × Holstein crossbred cattle

The improvement of milk production of indigenous Sudanese cattle such as Bos indicus Butana and its cross with Holstein is a major goal of the Sudanese government to ensure sufficient healthy nutrition in the country. In this study, we investigated the K232A polymorphism of diacylglycerol acyltransferase (DGAT1), a well-known modulator of milk production in other breeds. We determined allele frequencies and the allele effects on milk production. Therefore, 93 purebred Butana and 203 Butana × Holstein crossbred cattle were genotyped using competitive allele-specific PCR assays. Association analysis was performed using a linear mixed model in R. In purebred Butana cattle, the lysine DGAT1 protein variant K232, which is found to be associated with higher fat and protein contents, as well as higher fat yield was highly frequent at 0.929, while its frequency in Butana × Holstein crossbred cattle was 0.394. Significant effects were found on milk yield (P = 7.6 × 10−20), fat yield (P = 2.2 × 10−17), protein yield (P = 2.0 × 10−19) and lactose yield (P = 4.0 × 10−18) in crossbred cattle. As expected, the protein variant K232 was disadvantageous since it was decreasing milk, protein, and lactose yields by 1.741 kg, 0.063 kg and 0.084 kg, respectively. No significant effects were found for milk fat, protein, and lactose contents. The high frequency of the lysine DGAT1 protein variant K232 in Butana cattle could contribute to their high milk fat content in combination with low milk yield. In Butana × Holstein crossbred cattle, the DGAT1 marker can be used for effective selection and thus genetic improvement of milk production. Supplementary Information The online version contains supplementary material available at 10.1007/s11250-022-03141-7.


Introduction
The improvement of milk production is a major goal of the Sudanese government to ensure sufficient and healthy nutrition and to meet the increasing demand for milk and dairy products in Sudan. One of the indigenous breeds in Sudan is the Bos indicus Zebu breed Butana which is known for its good milk yield and adaption to the extreme environmental conditions. To increase the productivity of Butana cattle, they are crossed with the high-yielding dairy breed Holstein-Friesian. Moreover, the productivity of Butana as well as Butana × Holstein cattle could be improved by developing sustainable management systems, breeding programs, and also genetic selection for higher milk yield.
The DGAT1 (diacylglycerol acyltransferase) gene is located on bovine chromosome 14. The gene affects milk production traits in a wide-range of cattle breeds. Especially, the non-synonymous mutation K232A in exon 8 of DGAT1 which causes an amino acid change from lysine to alanine has a strong effect on milk yield and composition (Grisart et al., 2002;Spelman et al., 2002;Winter et al., 2002). The lysine DGAT1 protein variant K232 was found to be associated with higher fat and protein contents, as well as higher fat yield, but lower milk and protein yields (Bobbo et al., 2018;Bovenhuis et al., 2016;Li et al., 2021;Grisart et al., 2002;Spelman et al., 2002). Because of the strong negative relationship between lactose and fat content in milk across species (Fox, 2009), it is assumed that the lysine variant K232 also decreases lactose content. The frequency of the lysine variant K232 is high in Bos indicus cattle ranging from 0.80 to 1 (Abu Sara et al., 2015;Rahamattalla et al., 2015;Tantia et al., 2006;Kaupe et al., 2004;Spelman et al., 2002), but low in Holstein-Friesian cattle ranging from 0.21 to 0.63 (Grisart et al., 2002;Winter et al., 2002;Molee et al., 2012;Patel et al., 2009).
In this study, we investigated the K232A protein polymorphism of DGAT1 using the DNA marker rs109234250 (14:611,019 G/A) in Sudanese purebred Butana cattle as well as Butana × Holstein crossbred cattle. We determined its allele frequencies and performed association tests with milk production traits to test the suitability of the DGAT1 polymorphism as a marker for increasing milk yield and food security in Sudan. . DNA was extracted from blood using a salting out procedure (Miller et al., 1988). The K232A polymorphism of DGAT1 was investigated using the marker rs109234250 (14:611,019 G/A), whereby the allele A of this marker corresponds to the protein variant K232 with lysine. Genotypes were determined using allele-specific genotyping with competitive allele-specific PCR assays (KASP) (Kreuzer et al., 2013) (Supplementary Table 1). The genomic position is given with regard to the Bos taurus ARS_UCD1.2 genome assembly.

Phenotypic data
One thousand one hundred twenty-seven test-day records for 203 Butana × Holstein crossbred cattle from five farms and 134 test-day records for 37 Butana cattle from Atbara Research Station were used in this study. Daily milk yield (MY) and milk composition were recorded monthly in 2017 and 2018. Samples of fresh milk were analysed using the milk analyser Lactoscan for fat (FC), protein (PC) and lactose (LC) contents in percentages. Fat yield (FY), protein yield (PY) and lactose yield (LY) in kilograms were calculated by multiplying the respective contents with MY. Based on the calving and sampling dates, the age at first calving, lactation number and days in milk were calculated. The total MY per lactation and the lactation length were recorded from the latest finished lactation of each cow. On average, 5.25 measurements per cow were recorded. Outliers were defined as values outside the mean ± 3 standard deviations and removed from the data set.

Statistical analysis
Association analysis was performed only for Butana × Holstein crossbred cattle because the investigated DGAT1 marker was not segregating among the 37 Butana cattle, which had phenotypic data available. Association analysis between phenotypic and genotypic data was performed with a linear mixed model using the lmer function implemented in R language for statistical computing (version 4.0.3). The model for testing the additive effect of the DGAT1 marker was: where Y is the test-day record of the respective trait, X is the matrix of significantly estimated fixed effects, β is the vector of the fixed-effects regression coefficients, DGAT1 is the genotype at the DGAT1 marker rs109234250 and (1|animal) is the animal as a random effect to compensate for repeated measurements of the animals followed by the error term. The fixed effects farm, lactation number, year of sampling, season of sampling, birth year, birth season, calving season, age at first calving in days and days in milk and were only included in the model, if they significantly contributed to the model tested using ANOVA (Supplementary Table 2).

Overview of the investigated traits
The average MY per day was 2.4 times lower in Butana cattle (4.15 ± 1.40 kg) compared to Butana × Holstein cattle (10.12 ± 1.10 kg) (
The association study provides evidence for DGAT1 gene effects on milk production traits in the investigated Butana × Holstein crossbred cattle. The direction of effects of the lysine variant K232 on MY, PY, and LY was the same as in previous studies, decreasing them (Grisart et al., 2002(Grisart et al., , 2004Spelman et al., 2002;Gautier et al., 2007;Schennink et al., 2007;Streit et al., 2011;Cardoso et al., 2011). We observed that the lysine variant K232 decreased FY, in contrast to what has been reported in literature. However, this might be due to the overall low MY in the crossbreed animals causing also a slight decrease in FY. Furthermore, the very high correlations between the associated traits in our study (r > 0.89) may justify the same direction of effects for MY, FY, PY and LY.
Different levels of Holstein genomic background in the different crossbred populations and differences in the management of the herds influence the MY of individual animals. We observed that cows at Abbas Alreef Farm produced around twice as much milk (on average 17.16 kg/ day) as the cows on the other Butana × Holstein farms (Table 1). One reason for this are the excellent conditions with regard to farm management and feeding at Abbas Alreef Farm, which allows to make use of the genetic potential of Butana × Holstein crossbred cattle despite the harsh environmental conditions in Sudan. That is not to be taken for granted in Sudan considering the costs of feeding and management challenges at keeping high-yielding cattle. Thus, high variation of milk production traits and lactation properties between different farms is normal in Sudan. However, the average values reported in this study were in similar ranges compared to previously reported values for Butana × Holstein crossbred cattle (Supplementary Table 4) and Butana cattle (Supplementary Table 5). Another reason for the high MY is the higher amount of Holstein genomic background of crossbred cattle at the Abbas Alreef Farm compared to the other Butana × Holstein farms. This could be seen also in the allele frequency of the investigated DGAT1 marker for the K232A polymorphism, where alanine variant 232A is the favourite variant in Holstein cattle: At Abbas Alreef Farm, where the Holstein genomic background is the highest, the frequency of the alanine variant 232A, which is beneficial for MY, is the highest. All in all, a negative trend between the frequency of lysine variant K232, and the daily MY in the five investigated farms is visible. It would be interesting to know to which extent the DGAT1 marker could serve as a marker for the Holstein genomic background in crossbreeds with Holstein, where it is clear that this SNP was introduced from the Holstein breed. Furthermore, it would be desirable to consider the confounding factor of the Holstein genomic background. Unfortunately, neither pedigree data nor whole-genome data (e.g. from 50K SNP chip) are routinely collected in Sudan, so that we could not consider or investigate this factor in our current analyses. The only possibility to compensate for the farm-effect was to include "farm" as covariate in the association model. This farm-effect captured all differences between farms including farm management as well as Holstein genomic background.
Despite the small sample size in this study, we could confirm significant effects of DGAT1 on Butana × Holstein crossbred cattle. With a higher sample size, it might be also possible to obtain significant associations for fat and protein contents in Butana × Holstein crossbred cattle.
In conclusion, the high observed frequency of the lysine DGAT1 protein variant K232 in Butana cattle could explain their higher milk FC and lower MY. In Butana × Holstein crossbred cattle, the DGAT1 marker can be used for selection and based on our results should lead to an increase in MY. Nonetheless, other factors such as the farm management and feed supply have to be improved in parallel to increase production efficiency.
Author contribution GAB and SE designed the study. SE collected the samples and phenotypic data and did the laboratory work. MR helped with the laboratory work. SE and PK performed all the statistical analyses. SE and PK drafted the manuscript. GAB and DA helped draft the manuscript. All the authors read and approved the final manuscript.
Funding Open Access funding enabled and organized by Projekt DEAL. SE was supported by the Sudanese Ministry of Higher Education and Scientific Research and the Fiat Panis Foundation. The collection of data and samples was supported by the Ministry of Higher Education and Scientific Research Sudan.
Data availability All phenotypic and genotypic input data used in this study can be found in Supplementary Table 6.

Declarations
Ethics approval Blood samples of animals were collected based on routine procedures on farms, according to the animal protection law in Sudan.

Consent to participate
Animal owners agreed to the participation of their animals in this study.

Conflict of interest The authors declare no competing interests.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.