Genetic Parameters of Fatty Acids in Italian Brown Swiss and Holstein Cows

Uncovering genetic parameters and environmental influences on fertility, milk production, and quality in autochthonous Reggiana cattle

Reggiana is a local cattle breed from northern Italy known for its rusticity and profitability, due to the production of branded Parmigiano Reggiano cheese. To ensure the persistence of such profitability in the long term, an adequate breeding program is required. To this aim, in the present study we estimate the genetic parameters of the main productive and reproductive traits, and we evaluate the effect of genotype by environment interaction (GxE) on these traits using 2 environmental covariates: (1) productivity and (2) temperature-humidity index (THI). Milk, fat, protein, and casein yield were considered as daily production traits, whereas protein, fat, casein percentage, casein index, and somatic cell score were considered as milk quality traits. Finally, reproductive traits such as the number of inseminations, days open, calving interval, and calving-to-first-insemination interval were evaluated. Reggiana cattle produce an average of 19 kg of milk per day with 3.7% fat and 3.4% protein content and have excellent fertility parameters. Compared with other breeds, they have slightly lower heritability for production and quality for production traits (e.g., 0.12 [0.09; 0.15] for milk yield), but similar heritability for fertility traits. Milk, protein, and fat daily yields are highly correlated but negatively correlated with the percentage of protein, fat, and casein, whereas fertility traits have an unfavorable genetic correlation with daily production traits. When considering productivity, a consistent amount of variability due to GxE was observed for all daily production traits, somatic cell count, and casein index. A modest amount of GxE was observed for fertility parameters, while the percentage of solid content showed almost no GxE effect. A similar situation occurred when considering the THI, but no GxE interaction was observed for reproduction traits. In conclusion, this study provides useful information for the implementation of accurate selection plans in this local breed, accounting for environmental plasticity measured through the consistent GxE interaction observed.

Genotype by environment interaction for fat and protein yields via reaction norms in Holstein cattle of southern Brazil

Our objective was to evaluate the genetic merit of Holstein cattle population in southern Brazil in response to variations in the regional temperature by analyzing the genotype by environment interaction using reaction norms. Fat yield (FY) and protein yield (PY) data of 67 360 primiparous cows were obtained from the database of the Paraná Holstein Breeders Association, Brazil (APCBRH). The regional average annual temperature was used as the environmental variable. A random regression model was adopted applying mixed models with Restricted Maximum Likelihood (REML) algorithm using WOMBAT software. The genetic merit of the 15 most representative bulls, depending on the temperature gradient, was evaluated. Heritability ranged from 0.21 to 0.27 for FY and from 0.14 to 0.20 for PY. The genetic correlation observed among the environmental gradients proved to be higher than 0.80 for both traits. Slight reranking of bulls for both traits was detected, demonstrating that non-relevant genotype by environment interaction for FY and PY were observed. Consequently, no inclusion of the temperature effect in the model of genetic evaluation in southern Brazilian Holstein breed is required.

Genome-Wide Association Study for Milk Fatty Acids in Holstein Cattle Accounting for the DGAT1 Gene Effect

The identification of genomic regions and candidate genes associated with milk fatty acids contributes to better understand the underlying biology of these traits and enables breeders to modify milk fat composition through genetic selection. The main objectives of this study were: (1) to perform genome-wide association analyses for five groups of milk fatty acids in Holstein cattle using a high-density (777K) SNP panel; and (2) to compare the results of GWAS accounting (or not) for the DGAT1 gene effect as a covariate in the statistical model. The five groups of milk fatty acids analyzed were: (1) saturated (SFA); (2) unsaturated (UFA); (3) short-chain (SCFA); (4) medium-chain (MCFA); and (5) long-chain (LCFA) fatty acids. When DGAT1 was not fitted as a covariate in the model, significant SNPs and candidate genes were identified on BTA5, BTA6, BTA14, BTA16, and BTA19. When fitting the DGAT1 gene in the model, only the MGST1 and PLBD1 genes were identified. Thus, this study suggests that the DGAT1 gene accounts for most of the variability in milk fatty acid composition and the PLBD1 and MGST1 genes are important additional candidate genes in Holstein cattle.

Genetic and environmental relationships of detailed milk fatty acids profile determined by gas chromatography in Brown Swiss cows

The aim of this study was to characterize the profile of 47 fatty acids, including conjugated linoleic acid (CLA), 13 fatty acid groups, and 5 Δ(9)-desaturation indices in milk samples from Brown Swiss cows. The genetic variation was assessed and the statistical relevance of the genetic background for each trait was evaluated using the Bayes factor test. The additive genetic, herd-date, and residual relationships were also estimated among all single fatty acids and groups of fatty acids. Individual milk samples were collected from 1,158 Italian Brown Swiss cows and a detailed analysis of fat percentages and milk fatty acid compositions was performed by gas chromatography. Bayesian animal models were used for (co)variance components estimation. Exploitable genetic variation was observed for most of the de novo synthesized fatty acids and saturated fatty acids, except for C4:0 and C6:0, whereas long-chain fatty acids and unsaturated fatty acids (including CLA) were mainly influenced by herd-date effects. Herd-date effect explained large portions of the total phenotypic variance for C18:2 cis-9,cis-12 (0.668), C18:3 cis-9,cis-12,cis-15 (0.631), and the biohydrogenation and elongation products of these fatty acids. The desaturation ratios showed higher heritability estimates than the individual fatty acids, except for CLA desaturation index (0.098). Among the medium-chain fatty acids, C12:0 had greater heritability than C14:0 (0.243 vs. 0.097, respectively). Both C14:0 and C16:0 showed negative additive genetic correlations with the main monounsaturated and polyunsaturated fatty acids of milk fat, suggesting that their synthesis in the mammary gland may be influenced by the presence of unsaturated fatty acids. No correlation was observed between C4:0 and the other short-chain fatty acids (except for C6:0), confirming the independence of C4:0 from de novo mammary fatty acid synthesis. Among the genetic correlations dealing with potentially beneficial fatty acids, C18:0 was positively correlated with vaccenic and rumenic acids and negatively with linoleic acid. Finally, fatty acids C6:0 through C14:0 showed relevant correlations due to unknown environmental effects, suggesting the potential existence of genetic variances in micro-environmental sensitivity. This study allowed us to acquire new knowledge about the genetic and the environmental relationships among fatty acids. Likewise, the existence of genetic variation for most of de novo synthetized fatty acids and saturated fatty acids was also observed. Overall, these results provide useful information to combine feeding with genetic selection strategies for obtaining a desirable milk fatty acids profile, depending on the origin of fatty acids in milk.

Variation of milk components in the Italian Brown cattle

The aim of this study was to evaluate the variations of protein, casein, saturated (SFA), unsaturated (UFA), monounsaturated (MUFA), polyunsaturated (PUFA) fatty acids contents and cheese yield in the milk of two groups of Italian Brown cows conventionally reared in indoor period of housing or consuming pasture during the summer months in 2008 and 2013. Milk components were obtained from samples collected during the national routine (conventionally reared) and 'extraordinary' (pasture period) milk recording scheme in herds located near Sondrio (Lombardia, Italy). Milk samples were processed with the MilkoScanTM FT6000 for the identification of milk casein, SFA, UFA, MUFA and PUFA composition. The groups were analysed separately per year and the environmental factors affecting milk protein, casein, and fatty acids contents (pasture/indoor, parity, data of sampling, days in milk, days from collection to analysis) were included in the MIXED procedure of SAS 9.3. A total of 778 milk samples were available, including 234 records from indoor and 544 observations from pasture feeding. Pasture intake affected the content of casein (%) and the proportion of fat in milk (g/100 g), enhancing milk casein levels (from 2.90 to 3) and reducing the concentration of milk SFA in milk from grazing cows (from 2.29 to 1.92). Additionally, the cheese yield was calculated as 'kg of cheese per 100 kg of milk' and resulted to be 10.4 and 12 in 2008 from milk of cows reared indoor and with pasture based diet, respectively. The dairy industry should take advantage of the milk production during grazing periods from which high quality products may be obtained.