A Portable Infrared Attenuated Total Reflection Spectrometer for Food Analysis

The analytical performance of a compact infrared attenuated total reflection spectrometer using a pyroelectric detector array has been evaluated and compared to a conventional laboratory Fourier transform infrared system for applications in food analysis. Analytical characteristics including sensitivity, repeatability, linearity of the calibration functions, signal-to-noise ratio, and spectral resolution have been derived for both approaches. Representative analytes of relevance in food industries (i.e., organic solvents, fatty acids, and mycotoxins) have been used for the assessment of the performance of the device and to discuss the potential of this technology in food and feed analysis.

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

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

Metagenomic, microbiological, chemical and sensory profiling of Caciocavallo Podolico Lucano cheese

In this study, Caciocavallo Podolico Lucano (CPL) cheese was deeply characterized for its bacterial community, chemical composition and sensory aspects. The entire cheese making process (from milk collection to ripened cheese) was performed by strictly applying the traditional protocol for CPL production in four dairy factories (A-D) representative of the production area. The vat made of wood represents the main transformation tool for CPL cheese production and the biofilms hosted onto the internal surfaces of all vats analyzed in this study were dominated by lactic acid bacteria. Total mesophilic microorganisms present in bulk milk (4.7-5.0 log CFU/ml) increased consistently after contact with the wooden vat surfaces (5.4-6.4 log CFU/ml). The application of Illumina sequencing technology identified barely 18 taxonomic groups among processed samples; streptococci and lactobacilli constituted the major groups of the wooden vat biofilms [94.74-99.70 % of relative abundance (RA)], while lactobacilli dominated almost entirely (94.19-100 % of total RA) the bacterial community of ripened cheeses. Except coagulase positive staphylococci, undesirable bacteria were undetectable. Among chemical parameters, significant variations were registered for unsaturated, monounsaturated, polyunsaturated fatty acids and antioxidant properties (significantly lower for CPL cheeses produced in factory B). The cheeses from factories A, C and D were characterized by a higher lactic acid and persistence smell attributes than factory B. This work indicated that the strict application of CPL cheese making protocol harmonized the main microbiological, physicochemical and sensory parameters of the final cheeses produced in the four factories investigated.

Heat stress on breeding value prediction for milk yield and composition of a Brazilian Holstein cattle population

Due to the high milk production of Holstein cows, many countries have chosen to import semen to improve local dairy herds. This strategy would be more effective if this semen was used in the same environment conditions in which the bulls were selected. If the effect of genotype by environment (G × E) interaction is not considered, the estimated breeding values (EBVs) may vary, potentially reducing the selection response. We evaluate the impact of heat stress on selection for milk yield and composition of Holstein cows using random regression models. To verify the interference of heat stress in milk yield (MY) and composition traits (fat, protein, total saturated, and total unsaturated fatty acids content in milk), temperature-humidity index (THI) on test-day milk records was used. The threshold value to divide the environments using test-day information from Brazilian Holstein cows was 72 units of THI, i.e., < 72 represented no heat stress and > 72 represented heat stress. Legendre polynomials of second-order (Leg 2) model and two lactation points (33 and 122 DIM) were used to estimate heritabilities and EBVs for five important dairy traits. The heritabilities of milk components and fatty acids were low (0.09-0.29), regardless of lactation period and degree of heat stress, with the exception of protein content (0.30-0.35). Fat content was the only milk component that was reduced according to the degree of heat stress and lactation period. The EBVs tended to decrease in heat stress conditions, thus animals with high genetic potential demonstrated evidence of G × E interaction. However, acclimatization of dairy cows to heat stress in the farm production systems may have been responsible for the low differences among genetic parameters and EBVs with and without heat stress found in this study.

Seasonal and lactational variations in fatty acid profile of milk in indigenous cattle

The present study was conducted on quality and bioactive components in milk of Kankrej and Sahiwal indigenous cattle breeds at Livestock Research Station, Kodemdesar and College of Veterinary and Animal Science, Bikaner situated in hot-arid region of Rajasthan. Representative milk samples were collected and brought to the laboratory. The variance analysis observed significant effect of season on SFA of Kankrej and Sahiwal, and revealed that the total MUFA in Kankrej differ significantly. Significant effect of season on PUFA in Sahiwal and Kankrej cattle was also observed. Significant effect of season on Omega fatty acids in Sahiwal and Kankrej cattle except Omega-9 fatty acid in Sahiwal cattle was seen. The analysis of variance found highly significant effect of season on SCFA whereas non-significant effect on LCFA in Kankrej and Sahiwal cattle. Effect of season on MCFA was found significant in Sahiwal whereas non-significant effect was found in Kankrej cattle. The analysis of variance observed significant effect of various lactation stage on MUFA, Omega-9 and MCFA fatty acids for Kankrej and Sahiwal cattle and on SCFA in Sahiwal, and non-significant effect in Kankrej whereas significant effect on LCFA in Kankrej and nonsignificant effect in Sahiwal cattle.

Random regression test-day models to describe milk production and fatty acid traits in first lactation Walloon Holstein cows

We investigated the use of different Legendre polynomial orders to estimate genetic parameters for milk production and fatty acid (FA) traits in the first lactation Walloon Holstein cows. The data set comprised 302,684 test-day records of milk yield, fat and protein contents, and FAs generated by mid-infrared (MIR) spectroscopy, C16:0 (palmitic acid), C18:1 cis-9 (oleic acid), LCFAs (long-chain FAs), SFAs (saturated FAs) and UFAs (unsaturated FAs) were studied. The models included random regression coefficients for herd-year of calving (h), additive genetic (a) and permanent environment (p) effects. The selection of the best random regression model (RRM) was based on the deviance information criterion (DIC), and genetic parameters were estimated via a Bayesian approach. For all analysed random effects, DIC values decreased as the order of the Legendre polynomials increased. Best-fit models had fifth-order (degree 4) for the p effect and ranged from second- to fifth-order (degree 1-4) for the a and h effects (LEGhap: LEG555 for milk yield and protein content; LEG335 for fat content and SFA; LEG545 for C16:0 and UFA; and LEG535 for C18:1 cis-9 and LCFA). Based on the best-fit models, an effect of overcorrection was observed in early lactation (5-35 days in milk [DIM]). On the contrary, third-order (LEG333; degree 2) models showed flat residual trajectories throughout lactation. In general, the estimates of genetic variance tended to increase over DIM, for all traits. Heritabilities for milk production traits ranged from 0.11 to 0.58. Milk FA heritabilities ranged from low-to-high magnitude (0.03-0.56). High Spearman correlations (>0.90 for all bulls and >0.97 for top 100) were found among breeding values for 155 and 305 DIM between the best RRM and LEG333 model. Therefore, third-order Legendre polynomials seem to be most parsimonious and sufficient to describe milk production and FA traits in Walloon Holstein cows.

Genome-wide associations for heat stress response suggest potential candidate genes underlying milk fatty acid composition in dairy cattle

Contents of milk fatty acids (FA) display remarkable alterations along climatic gradients. Detecting candidate genes underlying such alterations might be beneficial for the exploration of climate sensitivity in dairy cattle. Consequently, we aimed on the definition of FA heat stress indicators, considering FA breeding values in response to temperature-humidity index (THI) alterations. Indicators were used in GWAS, in ongoing gene annotations and for the estimation of chromosome-wide variance components. The phenotypic data set consisted of 39,600 test-day milk FA records from 5,757 first-lactation Holstein dairy cows kept in 16 large-scale German cooperator herds. The FA traits were C18:0, polyunsaturated fatty acids (PUFA), saturated fatty acids (SFA), and unsaturated fatty acids (UFA). After genotype quality control, 40,523 SNP markers from 3,266 cows and 930 sires were considered. Meteorological data from the weather station in closest herd distance were used for the calculation of maximum hourly daily THI, which were allocated to 10 different THI classes. The same FA from 3 stages of lactation were considered as different, but genetically correlated traits. Consequently, a 3-trait reaction norm model was used to estimate genetic parameters and breeding values for FA along THI classes, considering either pedigree (A) or genomic (G) relationship matrices. De-regressed proofs and genomic estimated breeding values at the intermediate THI class 5 and at the extreme THI class 10 were used as pseudophenotypes in ongoing genomic analyses for thermoneutral (TNC) and heat stress conditions (HSC), respectively. The differences in de-regressed proofs and in genomic estimated breeding values from both THI classes were pseudophenotypes for heat stress response (HSR). Genetic correlations between the same FA under TNC and HSC were smallest in the first lactation stage and ranged from 0.20 for PUFA to 0.87 for SFA when modeling with the A matrix, and from 0.35 for UFA to 0.86 for SFA when modeling with the G matrix. In the first lactation stage, larger additive genetic variances under HSC compared with TNC indicate climate sensitivity for C18:0, PUFA, and UFA. Climate sensitivity was also reflected by pronounced chromosome-wide genetic variances for HSR of PUFA and UFA in the first stage of lactation. For all FA under TNC, HSC, and HSR, quite large genetic variance proportions were explained by BTA14. In GWAS, 30 SNP (within or close to 38 potential candidate genes) overlapped for HSR of the different FA. One unique potential candidate gene (AMFR) was detected for HSR of PUFA, 15 for HSR of SFA (ADGRB1, DENND3, DUSP16, EFR3A, EMP1, ENSBTAG00000003838, EPS8, MGP, PIK3C2G, STYK1, TMEM71, GSG1, SMARCE1, CCDC57, and FASN) and 3 for HSR of UFA (ENSBTAG00000048091, PAEP, and EPPK1). The identified unique genes play key roles in milk FA synthesis and are associated with disease resistance in dairy cattle. The results suggest consideration of FA in combination with climatic responses when inferring genetic mechanisms of heat stress in dairy cows.

Climate sensitivity of milk production traits and milk fatty acids in genotyped Holstein dairy cows

The aim of this study was the evaluation of climate sensitivity via genomic reaction norm models [i.e., to infer cow milk production and milk fatty acid (FA) responses on temperature-humidity index (THI) alterations]. Test-day milk traits were recorded between 2010 and 2016 from 5,257 first-lactation genotyped Holstein dairy cows. The cows were kept in 16 large-scale cooperator herds, being daughters of 344 genotyped sires. The longitudinal data consisted of 47,789 test-day records for the production traits milk yield (MY), fat yield (FY), and protein yield (PY), and of 20,742 test-day records for 6 FA including C16:0, C18:0, saturated fatty acids (SFA), unsaturated fatty acids (UFA), monounsaturated fatty acids (MUFA), and polyunsaturated fatty acids (PUFA). After quality control of the genotypic data, 41,057 SNP markers remained for genomic analyses. Meteorological data from the weather station in closest herd distance were used for the calculation of maximum hourly daily THI. Genomic reaction norm models were applied to estimate genetic parameters in a single-step approach for production traits and FA in dependency of THI at different lactation stages, and to evaluate the model stability. In a first evaluation strategy (New_sire), all phenotypic records from daughters of genotyped sires born after 2010 were masked, to mimic a validation population. In the second strategy (New_env), only daughter records of the new sires recorded in the most extreme THI classes were masked, aiming at predicting sire genomic estimated breeding values (GEBV) under heat stress conditions. Model stability was the correlation between GEBV of the new sires in the reduced data set with respective GEBV estimated from all phenotypic data. Among all test-day production traits, PY responded as the most sensitive to heat stress. As observed for the remaining production traits, genetic variances were quite stable across THI, but genetic correlations between PY from temperate climates with PY from extreme THI classes dropped to 0.68. Genetic variances in dependency of THI were very similar for C16:0 and SFA, indicating marginal climatic sensitivity. In the early lactation stage, genetic variances for C18:0, MUFA, PUFA, and UFA were significantly larger in the extreme THI classes compared with the estimates under thermoneutral conditions. For C18:0 and MUFA, PUFA, and UFA in the middle THI classes, genetic correlations in same traits from the early and the later lactation stages were lower than 0.50, indicating strong days in milk influence. Interestingly, within lactation stages, genetic correlations for C18:0 and UFA recorded at low and high THI were quite large, indicating similar genetic mechanisms under stress conditions. The model stability was improved when applying the New_env instead of New_sire strategy, especially for FA in the first stage of lactation. Results indicate moderately accurate genomic predictions for milk traits in extreme THI classes when considering phenotypic data from a broad range of remaining THI. Phenotypically, thermal stress conditions contributed to an increase of UFA, suggesting value as a heat stress biomarker. Furthermore, the quite large genetic variances for UFA at high THI suggest the consideration of UFA in selection strategies for improved heat stress resistance.

A Functional 3' UTR Polymorphism of FADS2 Affects Cow Milk Composition through Modifying Mir-744 Binding

Simple Summary

Fatty acid desaturase 2 (FADS2) is the rate-limiting enzyme involved in the synthesis of long-chain polyunsaturated fatty acids (LC-PUFAs). Many studies have suggested that polymorphisms in the FADS2 gene can modify its delta-6 desaturase activity. However, much remains unknown in regards to the regulatory mechanisms interpreting how DNA variants influence the function of FADS2. A previous study has suggested that c.1571G>A is located within the miR-744 binding site, indicating that this single nucleotide polymorphism (SNP) may be functional. Therefore, the aim of the present study was to determine the association of SNP c.1571G>A with milk PUFAs content and to further elucidate how this SNP contributes to regulating FADS2 expression. Our results support the use of SNP c.1571G>A as a potential genetic marker in the selective breeding of cattle to increase beneficial FAs content in milk.

Abstract

This study determined the associations of FADS2 c.1571G>A with milk FAs content and revealed that cows with the GG genotype had improved levels of delta-6 desaturase substrates (linoleic acid, C18:2n-6; p < 0.001) and decreased levels of desaturase products (gamma-linolenic acid, C18:3n-6; p < 0.001), indicating a reduction in FADS2 expression or delta-6 desaturase activity caused by this polymorphism. Computer alignment demonstrated that c.1571G>A occurred within a potential miR-744 binding site. When the c.1571G allele was present, the luciferase activity of reporter constructs was significantly suppressed by miR-744, while no such effect was observed with the A allele. Overexpression of miR-744 in bovine mammary epithelial cells (with the 1571GG genotype) downregulated FADS2 expression at both mRNA and protein levels. In contrast, inhibition of endogenous miR-744 with a specific inhibitor dramatically upregulated FADS2 expression. Taken together, these lines of evidence indicated that the c.1571A minor allele abolished the ability of miR-744 to bind FADS2, with a consequent increase in FADS2 expression levels and synthesis of omega-6 LC-PUFAs.

Co-Expression Networks Reveal Potential Regulatory Roles of miRNAs in Fatty Acid Composition of Nelore Cattle

Fatty acid (FA) content affects the sensorial and nutritional value of meat and plays a significant role in biological processes such as adipogenesis and immune response. It is well known that, in beef, the main FAs associated with these biological processes are oleic acid (C18:1 cis9, OA) and conjugated linoleic acid (CLA-c9t11), which may have beneficial effects on metabolic diseases such as type 2 diabetes and obesity. Here, we performed differential expression and co-expression analyses, weighted gene co-expression network analysis (WGCNA) and partial correlation with information theory (PCIT), to uncover the complex interactions between miRNAs and mRNAs expressed in skeletal muscle associated with FA content. miRNA and mRNA expression data were obtained from skeletal muscle of Nelore cattle that had extreme genomic breeding values for OA and CLA. Insulin and MAPK signaling pathways were identified by WGCNA as central pathways associated with both of these fatty acids. Co-expression network analysis identified bta-miR-33a/b, bta-miR-100, bta-miR-204, bta-miR-365-5p, bta-miR-660, bta-miR-411a, bta-miR-136, bta-miR-30-5p, bta-miR-146b, bta-let-7a-5p, bta-let-7f, bta-let-7, bta-miR 339, bta-miR-10b, bta-miR 486, and the genes ACTA1 and ALDOA as potential regulators of fatty acid synthesis. This study provides evidence and insights into the molecular mechanisms and potential target genes involved in fatty acid content differences in Nelore beef cattle, revealing new candidate pathways of phenotype modulation that could positively benefit beef production and human consumption.

Comparative milk metabolite profiling for exploring superiority of indigenous Indian cow milk over exotic and crossbred counterparts

This study was planned to identify differences in the milk metabolite composition of Indian (Sahiwal), exotic (Holstein-Friesian) and their crossbred cows in intensive system of management. To mimic the management system of ancient India, indigenous cattle under extensive system (zero input) were also included. Holstein-Friesian (HF) had significantly higher amount of saturated fatty acids (SFA, 76.3%) as compared to the crossbred (73.3%) and Sahiwal (68.0%). HF had the highest concentration (42.7%) of hypercholesterolemic fatty acids and the maximum value (68.5) of athrogenecity index (AI). Sahiwal had the highest proportion (32.1%) of total unsaturated fatty acids (UFA). Mineral, vitamin, n-3 fatty acids and total amount of essential amino acids did not vary across the three groups. Milk of indigenous cattle maintained only on grazing had more favorable nutrient profile. It had low SFA (61.4%), high UFA (38.6%) and higher concentrations of both monounsaturated fatty acids (31.4%) and polyunsaturated fatty acids (7.2%). The n-6/n-3 ratio (2.7) and the AI (33.9) were significantly lower. Significantly higher concentrations of minerals (Zn, Fe, P and Cu) and vitamins except vitamin B5 were recorded in their milk. The study revealed that milk metabolite characteristics can be used to promote indigenous cattle.

A 100-Year Review: Identification and genetic selection of economically important traits in dairy cattle

Over the past 100 yr, the range of traits considered for genetic selection in dairy cattle populations has progressed to meet the demands of both industry and society. At the turn of the 20th century, dairy farmers were interested in increasing milk production; however, a systematic strategy for selection was not available. Organized milk performance recording took shape, followed quickly by conformation scoring. Methodological advances in both genetic theory and statistics around the middle of the century, together with technological innovations in computing, paved the way for powerful multitrait analyses. As more sophisticated analytical techniques for traits were developed and incorporated into selection programs, production began to increase rapidly, and the wheels of genetic progress began to turn. By the end of the century, the focus of selection had moved away from being purely production oriented toward a more balanced breeding goal. This shift occurred partly due to increasing health and fertility issues and partly due to societal pressure and welfare concerns. Traits encompassing longevity, fertility, calving, health, and workability have now been integrated into selection indices. Current research focuses on fitness, health, welfare, milk quality, and environmental sustainability, underlying the concentrated emphasis on a more comprehensive breeding goal. In the future, on-farm sensors, data loggers, precision measurement techniques, and other technological aids will provide even more data for use in selection, and the difficulty will lie not in measuring phenotypes but rather in choosing which traits to select for.

Role of Fatty Acids in Milk Fat and the Influence of Selected Factors on Their Variability-A Review

Fatty acids (FAs) of milk fat are considered to be important nutritional components of the diets of a significant portion of the human population and substantially affect human health. With regard to dairy farming, the FA profile is also seen as an important factor in the technological quality of raw milk. In this sense, making targeted modifications to the FA profile has the potential to significantly contribute to the production of dairy products with higher added value. Thus, FAs also have economic importance. Current developments in analytical methods and their increasing efficiency enable the study of FA profiles not only for scientific purposes but also in terms of practical technological applications. It is important to study the sources of variability of FAs in milk, which include population genetics, type of farming, and targeted animal nutrition. It is equally important to study the health and technological impacts of FAs. This review summarizes current knowledge in the field regarding sources of FA variability, including the impact of factors such as: animal nutrition, seasonal feed changes, type of animal farming (conventional and organic), genetic parameters (influence of breed), animal individuality, lactation, and milk yield. Potential practical applications (to improve food technology and consumer health) of FA profile information are also reviewed.

Fine-mapping of BTA17 using imputed sequences for associations with de novo synthesized fatty acids in bovine milk

The aim of this study was to fine-map a genomic region associated with milk fatty acids (FA) on Bos taurus autosome (BTA) 17. This genomic region has been discovered with 50,000 (50k) single nucleotide polymorphisms (SNP) imputed to 777,000 (777k) SNP. In this study, high-density genotypes were imputed to whole-genome sequences level to identify candidate gene(s) associated with milk FA composition on BTA17. Phenotypes and genotypes were available for 1,640 cows sampled in winter, and for 1,581 cows sampled in summer. Phenotypes consisted of gas chromatography measurements in winter and in summer milk samples of 6 individual FA and the indicator of de novo synthesis, C6:0-C14:0. Genotypes consisted of imputed 777k SNP, and 89 sequenced ancestors of the population of genotyped cows. In addition, 450 whole-genome sequences from the 1,000 Bull Genome Consortium were available. Using 495 Holstein-Friesian sequences as a reference population, the 777k SNP genotypes of the cows were imputed to sequence level. We then applied single-variant analyses with an animal model, and identified thousands of significant associations with C6:0, C8:0, C10:0, C12:0, C14:0, and C6:0-C14:0. For C8:0 in summer milk samples, the genomic region located between 29 and 34 Mbp on BTA17 revealed a total of 646 significant associations. The most significant associations [-log₁₀(P-value) = 7.82] were 8 SNP in perfect linkage disequilibrium. After fitting one of these 8 SNP as a fixed effect in the model, and re-running the single-variant analyses, no further significant associations were found for any of the 6 FA or C6:0-C14:0. These findings suggest that one polymorphism underlying this QTL on BTA17 influences multiple de novo synthesized milk FA. Thirteen genes in the QTL region were identified and analyzed carefully. Six out of the 8 SNP that showed the strongest associations were located in the La ribonucleoprotein domain family, member 1B (LARP1B) gene, and we suggest LARP1B as a primary candidate gene. Another gene of interest for this QTL region might be PKL4. None of these suggested candidate genes have previously been associated with milk fat synthesis or milk FA composition.

Genetic analysis of groups of mid-infrared predicted fatty acids in milk

The objective of this study was to investigate genetic variability of mid-infrared predicted fatty acid groups in Canadian Holstein cattle. Genetic parameters were estimated for 5 groups of fatty acids: short-chain (4 to 10 carbons), medium-chain (11 to 16 carbons), long-chain (17 to 22 carbons), saturated, and unsaturated fatty acids. The data set included 49,127 test-day records from 10,029 first-lactation Holstein cows in 810 herds. The random regression animal test-day model included days in milk, herd-test date, and age-season of calving (polynomial regression) as fixed effects, herd-year of calving, animal additive genetic effect, and permanent environment effects as random polynomial regressions, and random residual effect. Legendre polynomials of the third degree were selected for the fixed regression for age-season of calving effect and Legendre polynomials of the fourth degree were selected for the random regression for animal additive genetic, permanent environment, and herd-year effect. The average daily heritability over the lactation for the medium-chain fatty acid group (0.32) was higher than for the short-chain (0.24) and long-chain (0.23) fatty acid groups. The average daily heritability for the saturated fatty acid group (0.33) was greater than for the unsaturated fatty acid group (0.21). Estimated average daily genetic correlations were positive among all fatty acid groups and ranged from moderate to high (0.63-0.96). The genetic correlations illustrated similarities and differences in their origin and the makeup of the groupings based on chain length and saturation. These results provide evidence for the existence of genetic variation in mid-infrared predicted fatty acid groups, and the possibility of improving milk fatty acid profile through genetic selection in Canadian dairy cattle.

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

Background

Bovine milk is widely regarded as a nutritious food source for humans, although the effects of individual fatty acids on human health is a subject of debate. Based on the assumption that genomic selection offers potential to improve milk fat composition, there is strong interest to understand more about the genetic factors that influence the biosynthesis of bovine milk and the molecular mechanisms that regulate milk fat synthesis and secretion. For this reason, the work reported here aimed at identifying genetic variants that affect milk fatty acid composition in Norwegian Red cattle. Milk fatty acid composition was predicted from the nation-wide recording scheme using Fourier transform infrared spectroscopy data and applied to estimate heritabilities for 36 individual and combined fatty acid traits. The recordings were used to generate daughter yield deviations that were first applied in a genome-wide association (GWAS) study with 17,343 markers to identify quantitative trait loci (QTL) affecting fatty acid composition, and next on high-density and sequence-level datasets to fine-map the most significant QTL on BTA13 (BTA for Bos taurus chromosome).

Results

The initial GWAS revealed 200 significant associations, with the strongest signals on BTA1, 13 and 15. The BTA13 QTL highlighted a strong functional candidate gene for de novo synthesis of short- and medium-chained saturated fatty acids; acyl-CoA synthetase short-chain family member 2. However, subsequent fine-mapping using single nucleotide polymorphisms (SNPs) from a high-density chip and variants detected by resequencing showed that the effect was more likely caused by a second nearby gene; nuclear receptor coactivator 6 (NCOA6). These findings were confirmed with results from haplotype studies. NCOA6 is a nuclear receptor that interacts with transcription factors such as PPARγ, which is a major regulator of bovine milk fat synthesis.

Conclusions

An initial GWAS revealed a highly significant QTL for de novo-synthesized fatty acids on BTA13 and was followed by fine-mapping of the QTL within NCOA6. The most significant SNPs were either synonymous or situated in introns; more research is needed to uncover the underlying causal DNA variation(s).

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-017-0294-5) contains supplementary material, which is available to authorized users.

Amiata Donkey Milk Chain: Animal Health Evaluation and Milk Quality

This study presents an investigation of Amiata donkey health and quality of milk for human consumption. Thirty-one lactating dairy jennies were examined. The following samples were collected: faecal samples from the rectum of animals for parasitological examination; cervical swabs for the detection of bacteria causing reproductive disorders; and blood samples for serological diagnosis of main zoonotic (Brucella spp., Leptospira spp.) and donkey abortion agents (Brucella spp., Leptospira spp., Salmonella abortus equi, Equine viral arterithis virus, Equine herpesvirus type 1). In addition, individual milk samples were collected and analysed for mastitis-causing pathogens and milk quality. Regarding animal health, we detected a high prevalence of strongyle parasites in donkeys. It is very important to tackle parasitic diseases correctly. Selective control programmes are preferable in order to reduce anthelmintic drug use. For dairy donkeys, withdrawal periods from anthelmintic drugs need to be carefully managed, in accordance with EU and national regulations. The isolation of Staphylococcus aureus in milk highlights the importance of preventing contamination during milking, by adopting appropriate hygiene and safety practices at a farm level. Lysozyme activity was high compared to cow’s milk, contributing to the inhibitory activity against certain bacteria. Donkey milk was characterised by high lactose content, low caseins, low fat, higher levels of unsaturated fatty acids compared to ruminant milks. Unsaturated fatty acids and omega 3 fatty acids in particular have become known for their beneficial health effect, which is favourable for human diet. These characteristics make it suitable for infants and children affected by food intolerance/allergies to bovine milk proteins and multiple food allergies as well as for adults with dyslipidemias. It is also recommended to prevent cardiovascular diseases.

Changes throughout lactation in phenotypic and genetic correlations between methane emissions and milk fatty acid contents predicted from milk mid-infrared spectra

The aim of this study was to estimate phenotypic and genetic correlations between methane production (Mp) and milk fatty acid contents of first-parity Walloon Holstein cows throughout lactation. Calibration equations predicting daily Mp (g/d) and milk fatty acid contents (g/100 dL of milk) were applied on milk mid-infrared spectra related to Walloon milk recording. A total of 241,236 predictions of Mp and milk fatty acids were used. These data were collected between 5 and 305 d in milk in 33,555 first-parity Holstein cows from 626 herds. Pedigree data included 109,975 animals. Bivariate (i.e., Mp and a fatty acid trait) random regression test-day models were developed to estimate phenotypic and genetic parameters of Mp and milk fatty acids. Individual short-chain fatty acids (SCFA) and groups of saturated fatty acids, SCFA, and medium-chain fatty acids showed positive phenotypic and genetic correlations with Mp (from 0.10 to 0.16 and from 0.23 to 0.30 for phenotypic and genetic correlations, respectively), whereas individual long-chain fatty acids (LCFA), and groups of LCFA, monounsaturated fatty acids, and unsaturated fatty acids showed null to positive phenotypic and genetic correlations with Mp (from -0.03 to 0.13 and from -0.02 to 0.32 for phenotypic and genetic correlations, respectively). However, these correlations changed throughout lactation. First, de novo individual and group fatty acids (i.e., C4:0, C6:0, C8:0, C10:0, C12:0, C14:0, SCFA group) showed low phenotypic or genetic correlations (or both) in early lactation and higher at the end of lactation. In contrast, phenotypic and genetic correlations between Mp and C16:0, which could be de novo synthetized or derived from blood lipids, were more stable during lactation. This fatty acid is the most abundant fatty acid of the saturated fatty acid and medium-chain fatty acid groups of which correlations with Mp showed the same pattern across lactation. Phenotypic and genetic correlations between Mp and C17:0 and C18:0 were low in early lactation and increased afterward. Phenotypic and genetic correlations between Mp and C18:1 cis-9 originating from the blood lipids were negative in early lactation and increased afterward to become null from 18 wk until the end of lactation. Correlations between Mp and groups of LCFA, monounsaturated fatty acids, and unsaturated fatty acids showed a similar or intermediate pattern across lactation compared with fatty acids that compose them. Finally, these results indicate that correlations between Mp and milk fatty acids vary following lactation stage of the cow, a fact still often ignored when trying to predict Mp from milk fatty acid profile.

A post-GWAS confirming the SCD gene associated with milk medium- and long-chain unsaturated fatty acids in Chinese Holstein population

The stearoyl-CoA desaturase (delta-9-desaturase) gene encodes a key enzyme in the cellular biosynthesis of monounsaturated fatty acids. In our initial genome-wide association study (GWAS) of Chinese Holstein cows, 19 SNPs fell in a 1.8-Mb region (20.3-22.1 Mb) on chromosome 26 underlying the SCD gene and were highly significantly associated with C14:1 or C14 index. The aims of this study were to verify whether the SCD gene has significant genetic effects on milk fatty acid composition in dairy cattle. By resequencing the entire coding region of the bovine SCD gene, a total of six variations were identified, including three coding variations (g.10153G>A, g.10213T>C and g.10329C>T) and three intronic variations (g.6926A>G, g.8646G>A and g.16158G>C). The SNP in exon 3, g.10329C>T, was predicted to result in an amino acid replacement from alanine (GCG) to valine (GTG) in the SCD protein. An association study for 16 milk fatty acids using 346 Chinese Holstein cows with accurate phenotypes and genotypes was performed using the mixed animal model with the proc mixed procedure in sas 9.2. All six detected SNPs were revealed to be associated with six medium- and long-chain unsaturated fatty acids (P = 0.0457 to P < 0.0001), specifically for C14:1 and C14 index (P = 0.0005 to P < 0.0001). Subsequently, strong linkage disequilibrium (D' = 0.88-1.00) was observed among all six SNPs in SCD and the five SNPs (rs41623887, rs109923480, rs42090224, rs42092174 and rs42091426) within the 1.8-Mb region identified in our previous GWAS, indicating that the significant association of the SCD gene with milk fatty acid content traits reduced the observed significant 1.8-Mb chromosome region in GWAS. Haplotype-based analysis revealed significant associations of the haplotypes encompassing the six SCD SNPs and one SNP (rs109923480) in a GWAS with C14:1, C14 index, C16:1 and C16 index (P = 0.0011 to P < 0.0001). In summary, our findings provide replicate evidence for our previous GWAS and demonstrate that variants in the SCD gene are significantly associated with milk fatty acid composition in dairy cattle, which provides clear evidence for an increased understanding of milk fatty acid synthesis and enhances opportunities to improve milk-fat composition in dairy cattle.

Genetic parameters for milk fatty acids, milk yield and quality traits of a Holstein cattle population reared under tropical conditions

Information about genetic parameters is essential for selection decisions and genetic evaluation. These estimates are population specific; however, there are few studies with dairy cattle populations reared under tropical and sub-tropical conditions. Thus, the aim was to obtain estimates of heritability and genetic correlations for milk yield and quality traits using pedigree and genomic information from a Holstein population maintained in a tropical environment. Phenotypic records (n = 36 457) of 4203 cows as well as the genotypes for 57 368 single nucleotide polymorphisms from 755 of these cows were used. Covariance components were estimated using the restricted maximum likelihood method under a mixed animal model, considering a pedigree-based relationship matrix or a combined pedigree-genomic matrix. High heritabilities (around 0.30) were estimated for lactose and protein content in milk whereas moderate values (between 0.19 and 0.26) were obtained for percentages of fat, saturated fatty acids and palmitic acid in milk. Genetic correlations ranging from -0.38 to -0.13 were determined between milk yield and composition traits. The smaller estimates compared to other similar studies can be due to poor environmental conditions, which may reduce genetic variability. These results highlight the importance in using genetic parameters estimated in the population under evaluation for selection decisions.

Fatty acid profile differs between organic and conventionally produced cow milk independent of season or milking time

Differing amounts of fresh forage and concentrates fed, and level of input contributes to the differences reported in fatty acid (FA) composition of organic and conventionally produced cow milk. In many previous studies designed to investigate this phenomenon, comparisons were made between grazed organic cows and housed conventional cows. In the present study, we have investigated differences between organic and conventional milk produced using year-round pasture grazing, as practiced in New Zealand. The FA composition was determined in milk sampled at morning and evening milking in both spring and autumn. Samples were taken from 45 cows from the Massey University organic herd and compared with 50 cows from the corresponding conventional herd grazed and managed similarly at the same location. Forty-three out of 51 analyzed FA were influenced by season, whereas 28 were different between production systems. In addition, one-half were also different due to time of milking. Levels of linoleic acid and α-linolenic acid were higher in organic milk, whereas conjugated linoleic acid (CLA) and vaccenic acid were higher in conventional milk. The first 3 FA (linoleic acid, α-linolenic acid, and CLA) were more abundant in milk harvested during autumn, and the CLA concentration was also significantly influenced by time of milking. Our results confirm reports that the FA profile is affected by season and time of milking, and we also showed an effect due to the production system, when both sets of cows were kept continuously on pasture, even after taking milking time and seasonal effect into account.

Invited review: organic and conventionally produced milk-an evaluation of factors influencing milk composition

Consumer perception of organic cow milk is associated with the assumption that organic milk differs from conventionally produced milk. The value associated with this difference justifies the premium retail price for organic milk. It includes the perceptions that organic dairy farming is kinder to the environment, animals, and people; that organic milk products are produced without the use of antibiotics, added hormones, synthetic chemicals, and genetic modification; and that they may have potential benefits for human health. Controlled studies investigating whether differences exist between organic and conventionally produced milk have so far been largely equivocal due principally to the complexity of the research question and the number of factors that can influence milk composition. A main complication is that farming practices and their effects differ depending on country, region, year, and season between and within organic and conventional systems. Factors influencing milk composition (e.g., diet, breed, and stage of lactation) have been studied individually, whereas interactions between multiple factors have been largely ignored. Studies that fail to consider that factors other than the farming system (organic vs. conventional) could have caused or contributed to the reported differences in milk composition make it impossible to determine whether a system-related difference exists between organic and conventional milk. Milk fatty acid composition has been a central research area when comparing organic and conventional milk largely because the milk fatty acid profile responds rapidly and is very sensitive to changes in diet. Consequently, the effect of farming practices (high input vs. low input) rather than farming system (organic vs. conventional) determines milk fatty acid profile, and similar results are seen between low-input organic and low-input conventional milks. This confounds our ability to develop an analytical method to distinguish organic from conventionally produced milk and provide product verification. Lack of research on interactions between several influential factors and differences in trial complexity and consistency between studies (e.g., sampling period, sample size, reporting of experimental conditions) complicate data interpretation and prevent us from making unequivocal conclusions. The first part of this review provides a detailed summary of individual factors known to influence milk composition. The second part presents an overview of studies that have compared organic and conventional milk and discusses their findings within the framework of the various factors presented in part one.

Genomic additive and dominance variance of milk performance traits

Milk performance traits are likely influenced by both additive and non-additive (e.g. dominance) genetic effects. Genetic variation can be partitioned using genomic information. The objective of this study was to estimate genetic variance components of production and milk component traits (e.g. acetone, fatty acids), which are particularly important for milk processing or which can provide information on the health status of cows. A genomic relationship approach was applied to phenotypic and genetic information of 1295 Holstein cows for estimating additive genetic and dominance variance components. Most of the 17 investigated traits were mainly affected by additive genetic effects, but protein content and casein content also showed a significant contribution of dominance. The ratio of dominance to additive variance was estimated as 0.64 for protein content and 0.56 for casein content. This ratio was highest for SCS (1.36) although dominance was not significant. Dominance effects were negligible in other moderately heritable milk traits.

Associations between variants of FADS genes and omega-3 and omega-6 milk fatty acids of Canadian Holstein cows

BACKGROUND

Fatty acid desaturase 1 (FADS1) and 2 (FADS2) genes code respectively for the enzymes delta-5 and delta-6 desaturases which are rate limiting enzymes in the synthesis of polyunsaturated omega-3 and omega-6 fatty acids (FAs). Omega-3 and-6 FAs as well as conjugated linoleic acid (CLA) are present in bovine milk and have demonstrated positive health effects in humans. Studies in humans have shown significant relationships between genetic variants in FADS1 and 2 genes with plasma and tissue concentrations of omega-3 and-6 FAs. The aim of this study was to evaluate the extent of sequence variations within these two genes in Canadian Holstein cows as well as the association between sequence variants and health promoting FAs in milk.

RESULTS

Thirty three SNPs were detected within the studied regions of genes including a synonymous mutation (FADS1-07, rs42187261, 306Tyr > Tyr) in exon 8 of FADS1, a non-synonymous mutation (FADS2-14, rs211580559, 294Ala > Val) within FADS2 exon 7, a splice site SNP (FADS2-05, rs211263660), a 3'UTR SNP (FADS2-23, rs109772589), and another 3'UTR SNP with an effect on a microRNA binding site within FADS2 gene (FADS2-19, rs210169303). Association analyses showed significant relations between three out of seven tested SNPs and several FAs. Significant associations (FDR P < 0.05) were recorded between FADS2-23 (rs109772589) and two omega-6 FAs (dihomogamma linolenic acid [C20:3n6] and arachidonic acid [C20:4n6]), FADS1-07 (rs42187261) and one omega-3 FA (eicosapentaenoic acid, C20:5n3) and tricosanoic acid (C23:0), and one intronic SNP, FADS1-01 (rs136261927) and C20:3n6.

CONCLUSION

Our study has demonstrated positive associations between three SNPs within FADS1 and FADS2 genes (a SNP within the 3'UTR, a synonymous SNP and an intronic SNP), with three milk PUFAs of Canadian Holstein cows thus suggesting possible involvement of synonymous and non-coding region variants in FA synthesis. These SNPs may serve as potential genetic markers in breeding programs to increase milk FAs that are of benefit to human health.

Concordance analysis between estimation methods of milk fatty acid content

Considering the milk fatty acid influence on human health, the aim of this study was to compare gas chromatography (GC) and Fourier transform infrared (FTIR) spectroscopy for the determination of these compounds. Fatty acid content (g/100g of fat) were obtained by both methods and compared through Pearson's correlation, linear Bayesian regression, and the Bland-Altman method. Despite the high correlations between the measurements (r=0.60-0.92), the regression coefficient values indicated higher measures for palmitic acid, oleic acid, unsaturated and monounsaturated fatty acids and lower values for stearic acid, saturated and polyunsaturated fatty acids estimated by GC in comparison to FTIR results. This inequality was confirmed in the Bland-Altman test, with an average bias varying from -8.65 to 6.91g/100g of fat. However, the inclusion of 94% of the samples into the concordance limits suggested that the variability of the differences between the methods was constant throughout the range of measurement. Therefore, despite the inequality between the estimates, the methods displayed the same pattern of milk fat composition, allowing similar conclusions about the milk samples under evaluation.

Genetic ancestry modifies fatty acid concentrations in different adipose tissue depots and milk fat

The objective of this study was to determine the effect of cow genetic strain on fatty acid (FA) profiles in adipose tissue and milk. Adipose samples from two subcutaneous (shoulder and tail-head) and three visceral (kidney channel, mesenteric and omental) depots were obtained post mortem from New Zealand (NZ; n = 8) and North American (NA; n = 8) Holstein-Friesian cows. At the time of slaughter cows were in similar body condition (NZ: 4.0 ± 0.03, NA: 4.0 ± 0.02; ± SD) and stage of lactation (NZ: 90 ± 11.2 d; NA: 83 ± 4.3 d; ± SD). Milk was collected during the a.m. milking prior to slaughter and milk fat was extracted. Adipose and milk fat FA were quantified using gas chromatography. NZ cows had a lower proportion of saturated FA in shoulder, tail-head and omental adipose tissue and a greater proportion of mono-unsaturated FA and an elevated Δ9-desaturase index in shoulder and tail-head adipose tissue. The proportions of individual FA differed between adipose depots, with proportions of de-novo FA greater in subcutaneous compared with visceral adipose depots. Milk from NZ cows contained greater concentrations of short chain FA (C8 : 0-12 : 0) and CLA, and less cis-9 18 : 1 than milk from NA cows. Regression analysis identified moderate associations between milk FA and shoulder adipose tissue FA for 18 : 2 (R(2) = 0.24), 18 : 3 n - 3 (R(2) = 0.39), and polyunsaturated fatty acids (R(2) = 0.38). Results from this study support the hypothesis that genetic strain dictates FA profiles in adipose tissue and milk and may alter the metabolic status of the various adipose depots differently. The data further support the premise that genetic strain affects the metabolic status of the various adipose depots differently. Elucidating the mechanisms that regulate the different adipose depots in the NZ and NA strains will increase our understanding of tissue mobilization and replenishment.

Variation in nutritionally relevant components in retail Jersey and Guernsey whole milk

This study investigates the quality of retail milk labelled as Jersey & Guernsey (JG) when compared with milk without breed specifications (NS) and repeatability of differences over seasons and years. 16 different brands of milk (4 Jersey & Guernsey, 12 non specified breed) were sampled over 2 years on 4 occasions. JG milk was associated with both favourable traits for human health, such as the higher total protein, total casein, α-casein, β-casein, κ-casein and α-tocopherol contents, and unfavourable traits, such as the higher concentrations of saturated fat, C12:0, C14:0 and lower concentrations of monounsaturated fatty acids. In summer, JG milk had a higher omega-3:omega-6 ratio than had NS milk. Also, the relative increase in omega-3 fatty acids and α-tocopherol, from winter to summer, was greater in JG milk. The latter characteristic could be of use in breeding schemes and farming systems producing niche dairy products. Seasonality had a more marked impact on the fatty acid composition of JG milk than had NS milk, while the opposite was found for protein composition. Potential implication for the findings in human health, producers, industry and consumers are considered.

New phenotypes for new breeding goals in dairy cattle

Cattle production faces new challenges regarding sustainability with its three pillars - economic, societal and environmental. The following three main factors will drive dairy cattle selection in the future: (1) During a long period, intensive selection for enhanced productivity has deteriorated most functional traits, some reaching a critical point and needing to be restored. This is especially the case for the Holstein breed and for female fertility, mastitis resistance, longevity and metabolic diseases. (2) Genomic selection offers two new opportunities: as the potential genetic gain can be almost doubled, more traits can be efficiently selected; phenotype recording can be decoupled from selection and limited to several thousand animals. (3) Additional information from other traits can be used, either from existing traditional recording systems at the farm level or from the recent and rapid development of new technologies and precision farming. Milk composition (i.e. mainly fatty acids) should be adapted to better meet human nutritional requirements. Fatty acids can be measured through a new interpretation of the usual medium infrared spectra. Milk composition can also provide additional information about reproduction and health. Modern milk recorders also provide new information, that is, on milking speed or on the shape of milking curves. Electronic devices measuring physiological or activity parameters can predict physiological status like estrus or diseases, and can record behavioral traits. Slaughterhouse data may permit effective selection on carcass traits. Efficient observatories should be set up for early detection of new emerging genetic defects. In the near future, social acceptance of cattle production could depend on its capacity to decrease its ecological footprint. The first solution consists in increasing survival and longevity to reduce replacement needs and the number of nonproductive animals. At the individual level, selection on rumen activity may lead to decreased methane production and concomitantly to improved feed efficiency. A major effort should be dedicated to this new field of research and particularly to rumen flora metagenomics. Low input in cattle production is very important and tomorrow's cow will need to adapt to a less intensive production environment, particularly lower feed quality and limited care. Finally, global climate change will increase pathogen pressure, thus more accurate predictors for disease resistance will be required.

Short communication: influence of the muscle hypertrophy mutation of the myostatin gene on milk production traits and milk fatty acid composition in dual-purpose Belgian Blue dairy cattle

The objective of this study was to estimate genetic effects for the muscle hypertrophy mutation (mh) of the myostatin gene for conventional milk production traits and for milk fatty acid composition in dual-purpose Belgian Blue dairy cows. For the present study, only cows from a single herd, in which genotype frequencies were as balanced as possible (0.266 for +/+, 0.523 for mh/+, and 0.211 for mh/mh), were chosen to avoid confounding between herd and genotype effects. A total of 109 cows with 3,190 test-day records for fat, protein, and milk yields and 1,064 test-day records for saturated and monounsaturated fatty acids were used for the calculations. Variance component and gene effect estimations were performed via expectation-maximization REML and BLUP methods, respectively, using a multi-trait mixed test-day model with an additional fixed regression on the muscle hypertrophy genotype. Results showed that one copy of the wild-type "+" allele led to a significant additive effect of 26.35 g/d for fat yield. Significant dominance effects of 23.22 g/d for protein yield and 30.28 g/d for fat yield were also observed. In contrast, a nonsignificant trend was observed in favor of lower saturated fatty acid contents in milk for one copy of the mutant "mh" allele. Concerning milk, fat, and protein yields, our results confirmed literature results indicating a superior effect of the "+" allele compared with the mutant allele. Therefore, the selection of the "+" allele has the potential to increase conventional milk production traits in the dual-purpose Belgian Blue breed. However, when focus is given to milk fatty acid profile, a possible antagonistic effect between the benefit of the "+" allele for higher milk production and the "mh" allele for reduced saturated fatty acid content in milk should be confirmed in further studies.

Mid-infrared prediction of bovine milk fatty acids across multiple breeds, production systems, and countries

Increasing consumer concern exists over the relationship between food composition and human health. Because of the known effects of fatty acids on human health, the development of a quick, inexpensive, and accurate method to directly quantify the fatty acid (FA) composition in milk would be valuable for milk processors to develop a payment system for milk pertinent to their customer requirements and for farmers to adapt their feeding systems and breeding strategies accordingly. The aim of this study was (1) to confirm the ability of mid-infrared spectrometry (MIR) to quantify individual FA content in milk by using an innovative procedure of sampling (i.e., samples were collected from cows belonging to different breeds, different countries, and in different production systems); (2) to compare 6 mathematical methods to develop robust calibration equations for predicting the contents of individual FA in milk; and (3) to test interest in using the FA equations developed in milk as basis to predict FA content in fat without corrections for the slope and the bias of the developed equations. In total, 517 samples selected based on their spectral variability in 3 countries (Belgium, Ireland, and United Kingdom) from various breeds, cows, and production systems were analyzed by gas chromatography (GC). The samples presenting the largest spectral variability were used to calibrate the prediction of FA by MIR. The remaining samples were used to externally validate the 28 FA equations developed. The 6 methods were (1) partial least squares regression (PLS); (2) PLS+repeatability file (REP); (3) first derivative of spectral data+PLS; (4) first derivative+REP+PLS; (5) second derivative of spectral data+PLS; and (6) second derivative+REP+PLS. Methods were compared on the basis of the cross-validation coefficient of determination (R2cv), the ratio of standard deviation of GC values to the standard error of cross-validation (RPD), and the validation coefficient of determination (R2v). The third and fourth methods had, on average, the highest R2cv, RPD, and R2v. The final equations were built using all GC and the best accuracy was observed for the infrared predictions of C4:0, C6:0, C8:0, C10:0, C12:0, C14:0, C16:0, C18:0, C18:1 trans, C18:1 cis-9, C18:1 cis, and for some groups of FA studied in milk (saturated, monounsaturated, unsaturated, short-chain, medium-chain, and long-chain FA). These equations showed R2cv greater than 0.95. With R2cv equal to 0.85, the MIR prediction of polyunsaturated FA could be used to screen the cow population. As previously published, infrared predictions of FA in fat are less accurate than those developed from FA content in milk (g/dL of milk) and no better results were obtained by using milk FA predictions if no corrections for bias and slope based on reference milk samples with known contents of FA were used. These results indicate the usefulness of equations with R2cv greater than 95% in milk payment systems and the usefulness of equations with R2cv greater than 75% for animal breeding purposes.

Predicting the fatty acid composition of milk: a comparison of two Fourier transform infrared sampling techniques

In the present study a novel approach for Fourier transform infrared (FT-IR) characterization of the fatty acid composition of milk based on dried film measurements has been presented and compared to a standard FT-IR approach based on liquid milk measurements. Two hundred and sixty-two (262) milk samples were obtained from a feeding experiment, and the samples were measured with FT-IR as dried films as well as liquid samples. Calibrations against the most abundant fatty acids, CLA (i.e., 18:2cis-9, trans-11), 18:3cis-9, cis-12, cis-15, and summed fatty acid parameters were obtained for both approaches. The estimation errors obtained in the dried film calibrations were overall lower than the corresponding liquid sample calibrations. Similar and good calibrations (i.e., R(2) ranges from 0.82 to 0.94 (liquid samples) and from 0.88 to 0.97 (dried films)) for short-chain fatty acids (6:0-14:0), 18:1cis-9, SAT, MUFA, and iodine value were obtained by both approaches. However, the dried film approach was the only approach for which feasible calibrations (i.e., R(2) ranges from 0.78 to 0.93) were obtained for the major saturated fatty acids 16:0 and 18:0, the minor fatty acid features 4:0, CLA (i.e., 18:2cis-9, trans-11), PUFA, and the summed 18:1 trans isomers. For the dried film approach, logical spectral features were found to dominate the respective fatty acid calibration models. The preconcentration step of the dried film approach could be expected to account for a major part of the prediction improvements going from predictions in liquid milk to predictions in dried films. The dried film approach has a significant potential for use in high-throughput applications in industrial environments and might also serve as a valuable supplement for determination of genetic and breeding factors within research communities.