Variation in Fatty Acid Contents of Milk and Milk Fat Within and Across Breeds

Overexpression of ATGL impairs lipid droplet accumulation by accelerating lipolysis in goat mammary epithelial cells

Adipose triglyceride lipase (ATGL) is the key enzyme for the degradation of triacylglycerols (TAGs). It functions in concert with other enzymes to mobilize TAG and supply fatty acids (FAs) for energy production. Dysregulated lipolysis leads to excess concentrations of circulating FAs, which may lead to destructive and lipotoxic effects to the organism. To understand the role of ATGL in mammary lipid metabolism, ATGL was overexpressed in goat mammary epithelial cells (GMECs) by using a recombinant adenovirus system. ATGL overexpression decreased lipid droplet (LD) accumulation and cellular TG content (p < 0.05) along with a decrease in the expression of the key enzyme that catalyzes the final step of TG synthesis (DGAT). Significant increases were observed in the expression of genes related to lipolysis (hormone-sensitive lipase [HSL]) and FA desaturation (SCD) by ATGL overexpression. Genes responsible for FA oxidation (PPARα), LD formation and secretion (ADRP and BTN1A1), and long-chain FA uptake (CD36) were all decreased by ATGL overexpression (p < 0.05). The primary products of TAG lipolysis, free FAs (FFAs), were notably increased in the ATGL-overexpressing cells. Taken together, our results demonstrated that ATGL activation impairs lipid formation partially through accelerating lipolysis in GMECs.

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.

Correlation between breeding values for milk fatty acids and Nordic Total Merit index traits for Danish Holstein and Danish Jersey

Milk fatty acid composition is gaining interest in the Danish dairy industry both to develop new dairy products and as a management tool. To be able to implement milk fatty acid (FA) composition in the breeding program, it is important to know the correlations with the traits in the breeding goal. To estimate these correlations, we measured milk fat composition in Danish Holstein (DH) and Danish Jersey (DJ) cattle breeds using mid-infrared spectroscopy. Breeding values were estimated for specific FA and for groups of FA. Correlations with the estimated breeding values (EBV) underlying the Nordic Total Merit index (NTM) were calculated within breed. For both DH and DJ, we showed that FA EBV had moderate correlations with the NTM and production traits. For both DH and DJ, the correlation of FA EBV and NTM were in the same direction, except for C16:0 (0 in DH, 0.23 in DJ). A few correlations differed between DH and DJ. The correlation between claw health index and C18:0 was negative in DH (-0.09) but positive in DJ (0.12). In addition, some correlations were not significant in DH but were significant in DJ. The correlations between udder health index and long-chain FA, trans FA, C16:0, and C18:0 were not significant in DH (-0.05 to 0.02), but were significant in DJ (-0.17, -0.15, 0.14, and -0.16, respectively). For both DH and DJ, the correlations between FA EBV and nonproduction traits were low. This implies that it is possible to breed for a different fat composition in the milk without affecting the nonproduction traits in the breeding goal.

Association of Single Nucleotide Polymorphism in the DGAT1 Gene with the Fatty Acid Composition of Cows Milked Once and Twice a Day

A single nucleotide polymorphism (SNP) rs109421300 of the diacylglycerol acyltransferase 1 (DGAT1) on bovine chromosome 14 is associated with fat yield, fat percentage, and protein percentage. This study aimed to investigate the effect of SNP rs109421300 on production traits and the fatty acid composition of milk from cows milked once a day (OAD) and twice a day (TAD) under New Zealand grazing conditions. Between September 2020 and March 2021, 232 cows from a OAD herd and 182 cows from a TAD herd were genotyped. The CC genotype of SNP rs109421300 was associated with significantly (p < 0.05) higher fat yield, fat percentage, and protein percentage, and lower milk and protein yields in both milking frequencies. The CC genotype was also associated with significantly (p < 0.05) higher proportions of C16:0 and C18:0, higher predicted solid fat content at 10 °C (SFC₁₀), and lower proportions of C4:0 and C18:1 cis-9 in both milking frequencies. The association of SNP with fatty acids was similar in both milking frequencies, with differences in magnitudes. The SFC₁₀ of cows milked OAD was lower than cows milked TAD for all three SNP genotypes suggesting the suitability of OAD milk for producing easily spreadable butter. These results demonstrate that selecting cows with the CC genotype is beneficial for New Zealand dairy farmers with the current payment system, however, this would likely result in less spreadable butter.

Health-Promoting Ingredients in Goat's Milk and Fermented Goat's Milk Drinks

The present study aimed to determine the content of health-promoting compounds, and fatty acids, with particular emphasis on the content of cis9trans11 C18:2 (CLA) acid, selected minerals, folates in organic and commercial goat's milk and fermented goat's milk drinks. The analyzed milk and yoghurts had various contents of particular groups of fatty acids, CLA, minerals, and folates. Raw organic goat's milk had a significantly (p < 0.05) higher content of CLA (3.26 mg/g fat) compared to commercial milk (2.88 mg/g fat and 2.54 mg/g fat). Among the analyzed fermented goat's milk drinks, the highest CLA content (4.39 mg/g fat) was determined in commercial natural yoghurts, while the lowest one was in organic natural yoghurts (3.28 mg/g fat). The highest levels of calcium (1322.9-2324.4 µg/g), phosphorus (8148.1-11,309.9 µg/g), and copper (0.072-0.104 µg/g) were found in all commercial products and those of manganese (0.067-0.209 µg/g) in organic products. The contents of the other assayed elements (magnesium, sodium, potassium, iron, and zinc) did not depend on the production method, but only on the product type, i.e., the degree of goat's milk processing. The highest folate content in the analyzed milks was found in the organic sample (3.16 µg/100 g). Organic Greek yoghurts had a several times higher content of folates, reaching 9.18 µg/100 g, compared to the other analyzed fermented products.

Variety of rumen microbial populations involved in biohydrogenation related to individual milk fat percentage of dairy cows

Our objective was to investigate the contribution of the rumen microbiome on the individual milk fat percentage (MFP) of Holstein dairy cows under the same nutritional and management conditions. From 92 early lactation dairy cows, the top 10 with the highest MFP (HF; n = 10) and the last 10 with the lowest MFP (LF; n = 10) were selected for the study. As a result, the milk trans-10, cis-12 C18:2 content was significant lower in the HF group than that in the LF group (P < 0.001). The rumen acetate to propionate ratio was significant higher in the HF group than that in the LF group (P = 0.035). According to the results of 16S rRNA gene sequencing, a minor but significant difference existed between the groups (P = 0.040). Three genera of the family Lachnospiraceae and four genera of the order Bacteroidales were identified to be the biomarkers for the LF group and HF group in the LEfSe analysis, respectively. Three microbial modules enriched by the family Lachnospiraceae were positively related to the milk trans-10, cis-12 C18:2 content (r _s > 0.60, P < 0.05). According to the results of shotgun metagenome sequencing, three kinds of linoleic acid (LA) isomerase genes were present in the gene pools of the rumen microbiome. Among them, the relative abundance of Bifidobacterium LA isomerase (BBI) was higher in the HF group than that in the LF group (P = 0.007). Three metagenome-assembled genomes (MAGs) with LA isomerase genes were positively correlated to the milk trans-10, cis-12 C18:2 content (r _s > 0.40, P < 0.05). Furthermore, all of these three MAGs were found to be able to produce lactate. Taken together, these results indicate that the increased relative abundance of microbial population with the trans-10 biohydrogenation pathway within the rumen microbiome contributes to the decrease of MFP via the increase of rumen trans-10, cis-12 C18:2 production. This study provides a new perspective for the development of measures for improving the milking performance of dairy cows.

Polymorphism of a Highly Asymmetrical Triacylglycerol in Milk Fat: 1-Butyryl 2-Stearoyl 3-Palmitoyl-glycerol

Milk fat has more than 200 triacylglycerols (TAGs), which play a pivotal role in its crystallization behavior. Asymmetrical TAGs containing short butyryl chains contribute to a significant portion of milk fat TAGs. This work aims to elucidate the crystallization behavior of asymmetrical milk fat TAGs by employing the pure compound of 1-butyryl 2-stearoyl 3-palmitoyl-glycerol (BuSP). The structural evolution of BuSP after being cooled down to 20 °C from the melt is evaluated by small- and wide-angle X-ray scattering (SAXS and WAXS) and differential scanning calorimetry (DSC). The temporal structural observation shows that BuSP crystallizes into the α-form with short and long spacings of 4.10 and 56.9 Å, respectively, during the first hour of isothermal hold at 20 °C. The polymorphic transformation of the α to β' phase occurred after 4 h of isothermal hold, and the β'- to α-form fraction ratio was about 70:30 at the end of the isothermal experiment (18 h). Pure β'-form X-ray patterns are obtained from the BuSP powder with short spacings of 4.33, 4.14, and 3.80 Å, while the long spacing of 51.2 Å depicts a three-chain-length lamellar structure with a tilt angle of 32°. Corresponding DSC measurements display that BuSP crystallizes from the melt at 29.1 °C, whereas the melting of α- and β'-forms was recorded at 30.3 and 47.8 °C, respectively. In the absence of the β-form, the β'-polymorph is the most stable observed form in BuSP. This work exemplarily explains the crystallization behavior of asymmetrical milk fat TAGs and thus provides new insights into their role in overall milk fat crystallization.

Effect of Breed on the Fatty Acid Composition of Milk from Dairy Cows Milked Once and Twice a Day in Different Stages of Lactation

The objective of this study was to evaluate the effect of breed on the overall composition and fatty acid composition of milk from cows milked once a day (OAD) and twice a day (TAD) in different stages of lactation. Milk samples were taken from 39 Holstein-Friesian (F), 27 Jersey (J), and 34 Holstein-Friesian × Jersey (F × J) crossbred cows from a OAD milking herd and 104 F and 83 F × J cows from a TAD milking herd in early (49 ± 15 days in milk), mid (129 ± 12 days in milk), and late (229 ± 13 days in milk) lactation. Calibration equations to predict the concentrations of individual fatty acids were developed using mid-infrared (MIR) spectroscopy. There was a significant interaction between breed within the milking frequency and stage of lactation for the production traits and composition traits. Holstein-Friesian cows milked OAD produced milk with lower concentrations of C18:0 in early and mid lactations compared to F × J and J cows. Holstein-Friesian cows milked TAD produced lower concentrations of C18:0 in early lactation and lower concentrations of C16:0 and C18:0 in late lactation compared to F × J. Lower concentrations of these fatty acids would reduce the hardness of the butter when the milk is processed. In the OAD milking herd, F cows were superior for daily milk yield compared to J cows, but Jersey cows produced significantly (p < 0.05) higher percentages of fat and a higher concentration of C18:0 fatty acid. The relative concentrations of C18:0 and C18 cis-9 in F and J cows milked OAD imply there is no breed effect on the activity of delta-9-desaturase, whereas stages of lactation likely have an effect. These results can be used to assist with selecting breeds and cows that are suitable for either OAD or TAD milking, allowing closer alignment with milk processing needs.

Oilseed Supplementation Improves Milk Composition and Fatty Acid Profile of Cow Milk: A Meta-Analysis and Meta-Regression

Oilseed supplementation is a strategy to improve milk production and milk composition in dairy cows; however, the response to this approach is inconsistent. Thus, the aim of this study was to evaluate the effect of oilseed supplementation on milk production and milk composition in dairy cows via a meta-analysis and meta-regression. A comprehensive and structured search was performed using the following electronic databases: Google Scholar, Primo-UAEH and PubMed. The response variables were: milk yield (MY), atherogenic index (AI), Σ omega-3 PUFA, Σ omega-6 PUFA, fat, protein, lactose, linoleic acid (LA), linolenic acid (LNA), oleic acid (OA), vaccenic acid (VA), conjugated linoleic acid (CLA), unsaturated fatty acid (UFA) and saturated fatty acid (SFA) contents. The explanatory variables were breed, lactation stage (first, second, and third), oilseed type (linseed, soybean, rapeseed, cottonseed, and sunflower), way (whole, extruded, ground, and roasted), dietary inclusion level, difference of the LA, LNA, OA, forage and NDF of supplemented and control rations, washout period and experimental design. A meta-analysis was performed with the “meta” package of the statistical program R. A meta-regression analysis was applied to explore the sources of heretogeneity. The inclusion of oilseeds in dairy cow rations had a positive effect on CLA (+0.27 g 100 g−1 fatty acids (FA); p < 0.0001), VA (+1.03 g 100 g−1 FA; p < 0.0001), OA (+3.44 g 100 g−1 FA; p < 0.0001), LNA (+0.28 g 100 g−1 FA; p < 0.0001) and UFA (+8.32 g 100 g−1 FA; p < 0.0001), and negative effects on AI (−1.01; p < 0.0001), SFA (−6.51; p < 0.0001), fat milk (−0.11%; p < 0.001) and protein milk (−0.04%; p < 0.007). Fat content was affected by animal breed, lactation stage, type and processing of oilseed and dietary NDF and LA contents. CLA, LA, OA and UFA, desirable FA milk components, were affected by type, processing, and the intake of oilseed; additionally, the concentrations of CLA and VA are affected by washout and design. Oilseed supplementation in dairy cow rations has a positive effect on desirable milk components for human consumption. However, animal response to oilseed supplementation depends on explanatory variables related to experimental design, animal characteristics and the type of oilseed.

Genetic variants associated with two major bovine milk fatty acids offer opportunities to breed for altered milk fat composition

Background

Although bovine milk is regarded as healthy and nutritious, its high content of saturated fatty acids (FA) may be harmful to cardiovascular health. Palmitic acid (C16:0) is the predominant saturated FA in milk with adverse health effects that could be countered by substituting it with higher levels of unsaturated FA, such as oleic acid (C18:1cis-9). In this work, we performed genome-wide association analyses for milk fatty acids predicted from FTIR spectroscopy data using 1811 Norwegian Red cattle genotyped and imputed to a high-density 777k single nucleotide polymorphism (SNP)-array. In a follow-up analysis, we used imputed whole-genome sequence data to detect genetic variants that are involved in FTIR-predicted levels of C16:0 and C18:1cis-9 and explore the transcript profile and protein level of candidate genes.

Results

Genome-wise significant associations were detected for C16:0 on Bos taurus (BTA) autosomes 11, 16 and 27, and for C18:1cis-9 on BTA5, 13 and 19. Closer examination of a significant locus on BTA11 identified the PAEP gene, which encodes the milk protein β-lactoglobulin, as a particularly attractive positional candidate gene. At this locus, we discovered a tightly linked cluster of genetic variants in coding and regulatory sequences that have opposing effects on the levels of C16:0 and C18:1cis-9. The favourable haplotype, linked to reduced levels of C16:0 and increased levels of C18:1cis-9 was also associated with a marked reduction in PAEP expression and β-lactoglobulin protein levels. β-lactoglobulin is the most abundant whey protein in milk and lower levels are associated with important dairy production parameters such as improved cheese yield.

Conclusions

The genetic variants detected in this study may be used in breeding to produce milk with an improved FA health-profile and enhanced cheese-making properties.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00731-9.

Performance and milk quality parameters of Jersey crossbreds in low-input dairy systems

Previous work has demonstrated some benefit from alternative breeds in low-input dairying, although there has been no systematic analysis of the simultaneous effect of Jersey crossbreeding on productivity, health, fertility parameters or milk nutritional quality. This work aimed to understand the effects of, and interactions/interrelations between, dairy cow genotypes (Holstein-Friesian (HF), Holstein-Friesian × Jersey crossbreds (HF × J)) and season (spring, summer, autumn) on milk yield; basic composition; feed efficiency, health, and fertility parameters; and milk fatty acid (FA) profiles. Milk samples (n = 219) and breed/diet data were collected from 74 cows in four UK low-input dairy farms between March and October 2012. HF × J cows produced milk with more fat (+ 3.2 g/kg milk), protein (+ 2.9 g/kg milk) and casein (+ 2.7 g/kg milk); and showed higher feed, fat, and protein efficiency (expressed as milk, fat and protein outputs per kg DMI) than HF cows. Milk from HF × J cows contained more C4:0 (+ 2.6 g/kg FA), C6:0 (+ 1.9 g/kg FA), C8:0 (+ 1.3 g/kg FA), C10:0 (+ 3.0 g/kg FA), C12:0 (+ 3.7 g/kg FA), C14:0 (+ 4.6 g/kg FA) and saturated FA (SFA; + 27.3 g/kg milk) and less monounsaturated FA (MUFA; -23.7 g/kg milk) and polyunsaturated FA (− 22.3 g/kg milk). There was no significant difference for most health and fertility parameters, but HF × J cows had shorter calving interval (by 39 days). The superior feed, fat and protein efficiency of HF × J cows, as well as shorter calving interval can be considered beneficial for the financial sustainability of low-input dairy farms; and using such alternative breeds in crossbreeding schemes may be recommended. Although statistically significant, it is difficult to determine if differences observed between HF and HF × J cows in fat composition are likely to impact human health, considering average population dairy fat intakes and the relatively small difference. Thus, the HF × J cow could be used in low-input dairying to improve efficiency and productivity without impacting milk nutritional properties.

Chemical Composition, Fatty Acid Profile, and Lipid Quality Indices in Commercial Ripening of Cow Cheeses from Different Seasons

The aim of the study was to compare and demonstrate whether commercial rennet ripening cheeses available on the market in summer and winter differ in their chemical composition, fatty acid profile, content of cis9trans11 C18:2 (CLA) acid and other trans isomers of C18:1 and C18:2 acid and whether they are characterized by different values of lipid quality assessment indices. The experimental material consisted of rennet ripening of cheeses produced from cow's milk available in the Polish market. The first batch contained cheeses produced in winter and purchased from the market between May and June. The second batch contained cheeses produced in summer and purchased between November and December. Chemical composition was analyzed by FoodScan apparatus. The gas chromatography (GC) method was used to determine the content of fatty acids. Results obtained in the presented study indicate that the chemical composition, content of fatty acids trans isomers, and lipid quality indices varied between summer and winter cheeses. The summer cheeses were richer sources of MUFA and PUFA compared to winter cheeses. Summer cheeses were also characterized by lower content of SFA, higher content n - 3, lower n - 6/n - 3 ratio, and higher content of DFA. Higher contents of CLA and trans C18:1 and C18:2 were found in summer cheeses.

Whole Genome DNA Methylation Variations in Mammary Gland Tissues from Holstein Cattle Producing Milk with Various Fat and Protein Contents

Milk fat and protein contents are among key elements of milk quality, and they are attracting more attention in response to consumers' demand for high-quality dairy products. To investigate the potential regulatory roles of DNA methylation underlying milk component yield, whole genome bisulfite sequencing was employed to profile the global DNA methylation patterns of mammary gland tissues from 17 Canada Holstein cows with various milk fat and protein contents. A total of 706, 2420 and 1645 differentially methylated CpG sites (DMCs) were found between high vs. low milk fat (HMF vs. LMF), high vs. low milk protein (HMP vs. LMP), and high vs. low milk fat and protein (HMFP vs. LMFP) groups, respectively (q value < 0.1). Twenty-seven, 56 and 67 genes harboring DMCs in gene regions (denoted DMC genes) were identified for HMF vs. LMF, HMP vs. LMP and HMFP vs. LMFP, respectively. DMC genes from HMP vs. LMP and HMFP vs. LMFP comparisons were significantly overrepresented in GO terms related to aerobic electron transport chain and/or mitochondrial ATP (adenosine triphosphate) synthesis coupled electron transport. A total of 83 (HMF vs. LMF), 708 (HMP vs. LMP) and 408 (HMFP vs. LMFP) DMCs were co-located with 87, 147 and 158 quantitative trait loci (QTL) for milk component and yield traits, respectively. In conclusion, the identified methylation changes are potentially involved in the regulation of milk fat and protein yields, as well as the variation in reported co-located QTLs.

Association of DGAT1 With Cattle, Buffalo, Goat, and Sheep Milk and Meat Production Traits

Milk fatty acids are essential for many dairy product productions, while intramuscular fat (IMF) is associated with the quality of meat. The triacylglycerols (TAGs) are the major components of IMF and milk fat. Therefore, understanding the polymorphisms and genes linked to fat synthesis is important for animal production. Identifying quantitative trait loci (QTLs) and genes associated with milk and meat production traits has been the objective of various mapping studies in the last decade. Consistently, the QTLs on chromosomes 14, 15, and 9 have been found to be associated with milk and meat production traits in cattle, goat, and buffalo and sheep, respectively. Diacylglycerol O-acyltransferase 1 (DGAT1) gene has been reported on chromosomes 14, 15, and 9 in cattle, goat, and buffalo and sheep, respectively. Being a key role in fat metabolism and TAG synthesis, the DGAT1 has obtained considerable attention especially in animal milk production. In addition to milk production, DGAT1 has also been a subject of interest in animal meat production. Several polymorphisms have been documented in DGAT1 in various animal species including cattle, buffalo, goat, and sheep for their association with milk production traits. In addition, the DGAT1 has also been studied for their role in meat production traits in cattle, sheep, and goat. However, very limited studies have been conducted in cattle for association of DGAT1 with meat production traits in cattle. Moreover, not a single study reported the association of DGAT1 with meat production traits in buffalo; thus, further studies are warranted to fulfill this huge gap. Keeping in view the important role of DGAT1 in animal production, the current review article was designed to highlight the major development and new insights on DGAT1 effect on milk and meat production traits in cattle, buffalo, sheep, and goat. Moreover, we have also highlighted the possible future contributions of DGAT1 for the studied species.

The Impact of Seasonality in Pasture-Based Production Systems on Milk Composition and Functionality

Seasonal calving, pasture-based dairy systems are widely practiced in countries with a temperate climate and plentiful rainfall such as Ireland and New Zealand. This approach maximizes milk production from pasture and, consequently, is a low-cost, low-input dairy production system. On the other hand, the majority of global milk supply is derived from high input indoor total mixed ration systems where seasonal calving is not practiced due to the dependence on ensiled silages, grains and concentrated feeds, which are available year-round. Synchronous changes in the macro and micronutrients in milk are much more noticeable as lactation progresses through early, mid and late stages in seasonal systems compared to non-seasonal systems-which can have implications on the processability and functionality of milk.

Meeting Breeding Potential in Organic and Low-Input Dairy Farming

Low-input (LI) dairy farming, relying heavily on grazing, is increasing in popularity for perceived sustainability, welfare, and milk nutritional quality benefits. However, there is little research into the breed suitability for these systems. The popular Holstein–Friesians are not well-suited to LI production as, to achieve their potential high yields, they require high levels of concentrate intakes and veterinary inputs. Holstein–Friesians were traditionally bred for high milk yields, which often correlate negatively with functional traits, such as fertility and health. This drives the need for alternative breed choices, and UK dairy farmers use several crossbreeding practices. Additionally, classic measures of production efficiency (kilogram feed per liter of milk) are not the sole priority in LI systems, which also aim for improved health, fertility, forage conversion, and milk quality. This study aimed to explore the effect of breeding strategy on LI and organic production in dairy systems, collecting data from 17 farms throughout England and Wales: 7 organic and 10 low-input conventional systems with both purebred and crossbred cows from different breeds. Records from 1,070 cows were collected, including background data, health, fertility, breeding, and parity. Additionally, milk was analyzed on four occasions (autumn 2011 and winter, spring, and summer 2012). Principal components analysis was used to visualize the effect of management, Farm ID, and stage of lactation on LI production. The analysis clustered cows by Farm ID, showing that individual management practice on each farm had the greatest impact on various production traits. Cows were allocated a composite score based on their yield, health records, and milk fatty acid profile, and a linear mixed-effects model indicated (p < 0.01) that crossbred New Zealand Friesian cows scored highest, whereas Dairy Shorthorn cows scored the lowest. This paper highlights weaknesses in current breeding programs for LI and organic farms in the UK, in terms of the alignment of breeds with husbandry practices. Additional research is needed to explore any gene by environment interactions to meet the true potential of individual cows and certain breeds under LI and organic management.

Integration of Wet-Lab Measures, Milk Infrared Spectra, and Genomics to Improve Difficult-to-Measure Traits in Dairy Cattle Populations

The objective of this study was to evaluate the contribution of Fourier-transformed infrared spectroscopy (FTIR) data for dairy cattle breeding through two different approaches: (i) estimating the genetic parameters for 30 measured milk traits and their FTIR predictions and investigating the additive genetic correlation between them and (ii) evaluating the effectiveness of FTIR-derived phenotyping to replicate a candidate bull’s progeny testing or breeding value prediction at birth. Records were available from 1,123 cows phenotyped using gold standard laboratory methodologies (LAB data). This included phenotypes related to fine milk composition and milk technological characteristics, milk acidity, and milk protein fractions. The dataset used to generate FTIR predictions comprised 729,202 test-day records from 51,059 Brown Swiss cows (FIELD data). A first approach consisted of estimating genetic parameters for phenotypes available from LAB and FIELD datasets. To do so, a set of bivariate animal models were run, and genetic correlations between LAB and FIELD phenotypes were estimated using FIELD information obtained at the population level. Heritability estimates were generally higher for FIELD predictions than for the corresponding LAB measures. The additive genetic correlations (r_a) between LAB and FIELD phenotypes had different magnitudes across traits but were generally strong. Overall, these results demonstrated the potential of using FIELD information as indicator traits for the indirect genetic improvement of LAB measures. In the second approach, we included genotype information for 1,011 cows from the LAB dataset, 1,493 cows from the FIELD dataset, 181 sires with daughters in both LAB and FIELD datasets, and 540 sires with daughters in the FIELD dataset only. Predictions were obtained using the single-step GBLUP method. A four fold cross-validation was used to assess the predictive ability of the different models, assessed as the ability to predict masked LAB records from daughters of progeny testing bulls. The correlation between observed and predicted LAB measures in validation was averaged over the four training-validation sets. Different sets of phenotypic information were used sequentially in cross-validation schemes: (i) LAB cows from the training set; (ii) FIELD cows from the training set; and (iii) FIELD cows from the validation set. Models that included FIELD records showed an improvement for the majority of traits. This study suggests that breeding programs for difficult-to-measure traits could be implemented using FTIR information. While these programs should use progeny testing, acceptable values of accuracy can be achieved also for bulls without phenotyped progeny. Robust calibration equations are, deemed as essential.

The Effect of Mixing Milk of Different Species on Chemical, Physicochemical, and Sensory Features of Cheeses: A Review

The yield and quality of cheese are associated with the composition, physicochemical, sensory, rheological, and microbiological properties of milk and with the technology applied to the milk before and/or during cheese processing. This review describes the most important research on cheeses obtained from processing mixtures of different milk species and discusses the effect of milk mixtures (i.e., species and mixture ratios) on composition, physicochemical, sensory, rheological, and microbiological properties of cheeses. More specifically, the present review paper will gather and focus only on studies that have provided a clear comparison between cheeses produced from a mixture of two milk species to cheeses produced from only one species.

A comparative study on milk composition of Jersey and Holstein dairy cows during the early lactation

Recently, Jersey cattle was introduced and produced by embryo transfer to Korea. This study was conducted to investigate the differences of milk compositions between Jersey and Holstein cows and the relationship between days in milk (DIM) and milk compositions during early lactation. Data were collected from twelve lactating cows from Department of Animal Resources Development at National Institute of Animal Science. Cows in parity 1 were used, and calved at spring from April to March of 2017. All cows were housed in two sections within a free-stall barn, which divided into six from each breed, and received a basal total mixed ration. Milk samples of each cow were collected at 3 DIM and 30 DIM for analyzing the milk compositions, including fatty acids (FA), amino acids and minerals. Total solids, citrate, and milk urea nitrogen level were differed between the breeds (p < 0.05). As DIM went from 3 to 30, milk protein, total solids, and somatic cell count decreased (p < 0.05), but lactose increased in all breed milk (p < 0.05). Citrate and free fatty acid (FFA) elevated in Jersey milk (p < 0.05), whereas reduced in Holstein milk (p < 0.05). Proportions of some individual FA varied from the breeds. Myristic (C14:0), palmitic (C16:0), and arachidonic acid (C20:4) in milk from all cows were higher at 3 DIM than at 30 DIM (p < 0.05). Also, stearic (C18:0) and oleic acid (C18:1) were lower at 3 DIM than at 30 DIM (p < 0.05), and the C18:1 to C18:0 ratio was significantly differed in DIM × breed interactions (p < 0.05). The contents of the individual amino acids did not differ from the breeds. Calcium, phosphorous, magnesium, and zinc (Zn) contents was significantly increased in Holstein milk than Jersey milk at 3 DIM. Also, K and Zn concentrations were higher in Holstein milk than in Jersey milk at 30 DIM (p < 0.05). It was concluded that Jersey cows would produce more effective milk in processing dairy products and more proper energy status compared with Holstein cows in early lactation under the same environmental and nutritional conditions.

Genetic correlations of milk fatty acid contents predicted from milk mid-infrared spectra in New Zealand dairy cattle

The objective of this study was to estimate genetic correlations among milk fatty acid (FA) concentrations in New Zealand dairy cattle. Concentrations of each of the most common FA, expressed as a percentage of the total FA, were determined by gas chromatography on a specific cohort of animals. Using this data set, prediction equations were derived using mid-infrared (MIR) spectroscopy data collected from the same samples. These prediction equations were applied to a large data set of MIR measurements in 34,141 milk samples from 3,445 Holstein-Friesian, 2,935 Jersey, and 3,609 crossbred Holstein-Friesian × Jersey cows, sampled an average of 3.42 times during the 2007-2008 season. Data were analyzed using univariate and bivariate repeatability animal models. Heritability of predicted FA concentration in milk fat ranged from 0.21 to 0.42, indicating that genetic selection could be used to change the FA composition of milk. The de novo synthesized FA (C6:0, C8:0, C10:0, C12:0, and C14:0) showed strong positive genetic correlations with each other, ranging from 0.24 to 0.99. Saturated FA were negatively correlated with unsaturated (-0.93) and polyunsaturated (-0.84) FA. The saturated FA were positively correlated with milk fat yield and fat percentage, whereas the unsaturated FA were negatively associated with fat yield and fat percentage. Our results indicate that bovine milk FA composition can be changed through genetic selection using MIR as a phenotypic proxy.

Polymorphisms in Fatty Acid Desaturase 2 Gene are Associated with Milk Production Traits in Chinese Holstein Cows

Simple Summary

Searching for causative polymorphisms underlying the variability of milk production traits and then incorporating them into breeding programs are very effective ways to improve the efficiency and reliability of conventional dairy cattle breeding. Fatty acid desaturase 2 (FADS2) plays a pivotal role in the biosynthesis of polyunsaturated fatty acids. Previous studies provided evidence that FADS2 was one of the most downregulated genes during negative energy balance in the liver of postpartum dairy cattle. Genes involved in the energetic pathways may influence other production traits, such as protein, fat and milk yields. Therefore, in the present study, we investigated the common genetic variants of the FADS2 gene in Chinese Holstein cows. Our results provided direct evidence that FADS2 was an interesting candidate for selection to increase milk production and improve resistance against mastitis.

Abstract

This study investigated the single nucleotide polymorphisms (SNPs) of Fatty acid desaturase 2 (FADS2) gene and further explored their genetic effects on conventionally collected milk production traits in Chinese Holstein cows using 18,264 test-day records of 841 cows. One missense mutation c. 908 C > T (SNP site in the complementary DNA sequence), which caused an amino acid change from alanine to valine (294Ala > Val), and two 3’ untranslated region (UTR) SNPs, c.1571 G > A and c.2776 A > G were finally identified. The SNP c.908 C > T was significantly associated with test-day milk yield, fat percentage and 305-day milk, fat and protein yield. In particular, the T allele of the SNP c.908 C > T showed a significant association with decreased somatic cell score (SCS) in the investigated population. Significant relationship between the SNP c.1571 G > A and 305-day milk yield showed that genotype GG was linked to the highest milk yield. Substituting the allele G for A at the c.2776 A > G locus resulted in a decrease of protein percentage. Our results demonstrated that FADS2 was an interesting candidate for selection to increase milk production and improve resistance against mastitis.

Effects of FABP4 variation on milk fatty-acid composition for dairy cattle grazed on pasture in late lactation

This research communication describes associations between variation in the fatty acid binding protein 4 gene (FABP4) and milk fat composition in New Zealand Holstein-Friesian × Jersey cross dairy cows. After correcting for the effect of the amino acid substitution p.K232A in diacylglycerol O-acyltransferase 1 (DGAT1), which is associated with variation in many milk fatty acid (FA) component levels, the effect of FABP4 c.328A/G on milk FA levels was typically small. For the five genotypes analysed, the AB cows produced more medium-chain fatty acids than CC cows (P < 0.05), and more C14:0 FA than AA and AC cows (P < 0.05). The AA and AC cows produced less C22:0 FA (P < 0.01) than the BC cows, and the AC cows produced more C24:0 FA (P < 0.05) than was produced by the BC cows. Cows of genotype CC produce more long-chain fatty acids than cows of genotype BC (P < 0.05).

Genome-wide association study on Fourier transform infrared milk spectra for two Danish dairy cattle breeds

Background

Infrared spectral analysis of milk is cheap, fast, and accurate. Infrared light interacts with chemical bonds present inside the milk, which means that Fourier transform infrared milk spectra are a reflection of the chemical composition of milk. Heritability of Fourier transform infrared milk spectra has been analysed previously. Further genetic analysis of Fourier transform infrared milk spectra could give us a better insight in the genes underlying milk composition. Breed influences milk composition, yet not much is known about the effect of breed on Fourier transform infrared milk spectra. Improved understanding of the effect of breed on Fourier transform infrared milk spectra could enhance efficient application of Fourier transform infrared milk spectra. The aim of this study is to perform a genome wide association study on a selection of wavenumbers for Danish Holstein and Danish Jersey. This will improve our understanding of the genetics underlying milk composition in these two dairy cattle breeds.

Results

For each breed separately, fifteen wavenumbers were analysed. Overall, more quantitative trait loci were observed for Danish Jersey compared to Danish Holstein. For both breeds, the majority of the wavenumbers was most strongly associated to a genomic region on BTA 14 harbouring DGAT1. Furthermore, for both breeds most quantitative trait loci were observed for wavenumbers that interact with the chemical bond C-O. For Danish Jersey, wavenumbers that interact with C-H were associated to genes that are involved in fatty acid synthesis, such as AGPAT3, AGPAT6, PPARGC1A, SREBF1, and FADS1. For wavenumbers which interact with –OH, associations were observed to genomic regions that have been linked to alpha-lactalbumin.

Conclusions

The current study identified many quantitative trait loci that underlie Fourier transform infrared milk spectra, and thus milk composition. Differences were observed between groups of wavenumbers that interact with different chemical bonds. Both overlapping and different QTL were observed for Danish Holstein and Danish Jersey.

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.

Impact of crossbreeding Holstein and Brown Swiss cows on milk yield, composition, and fatty acid profiles in subtropics

Milk yield (MY), composition, and fatty acid profiles of purebred Holstein (HO) cows, Brown Swiss (BS), and their F₁ crossbreds (HS) were compared under subtropical climate in Egypt. Pure HO had significantly greater 305-MY, total-MY, and daily-MY (p < 0.05) than pure BS and crossbred HS. Furthermore, HO and HS had significantly higher peak-MY (44.2 and 43.3 kg, respectively) than BS (36.1 kg). In comparison with HO, BS had significantly higher milk fat, protein, total solids, and solid-not-fat % (p < 0.05); however, no significant differences were observed between BS and HS for milk fat %. The milk fat of BS had higher concentrations of saturated (C17:0 and C18:0) fatty acids (FA) than that of the HO and HS (p = 0.001 and 0.008, respectively). Content of C4:0, C6:0, C8:0, C12:0, and C15.0 FA did not differ between genotypes. Milk from HO and BS had significantly higher concentrations of unsaturated (C20:1 and C20:5) FA than that from HS (p < 0.05). C14:1 (myristoleic), C16:1 (palmitoleic), and C18:2 (linoleic) FA contents were similar for all genotypes. In conclusion, BS surpassed HO for all milk composition traits and some FA components; however, HS had comparable milk fat percentage with BS. There appears to be an opportunity to modify the concentration of certain FA by breeding.

Relationships between milk mid-IR predicted gastro-enteric methane production and the technical and financial performance of commercial dairy herds

Considering economic and environmental issues is important in ensuring the sustainability of dairy farms. The objective of this study was to investigate univariate relationships between lactating dairy cow gastro-enteric methane (CH4) production predicted from milk mid-IR (MIR) spectra and technico-economic variables by the use of large scale and on-farm data. A total of 525 697 individual CH4 predictions from milk MIR spectra (MIR-CH4 (g/day)) of milk samples collected on 206 farms during the Walloon milk recording scheme were used to create a MIR-CH4 prediction for each herd and year (HYMIR-CH4). These predictions were merged with dairy herd accounting data. This allowed a simultaneous study of HYMIR-CH4 and 42 technical and economic variables for 1024 herd and year records from 2007 to 2014. Pearson correlation coefficients (r) were used to assess significant relationships (P<0.05). Low HYMIR-CH4 was significantly associated with, amongst others, lower fat and protein corrected milk (FPCM) yield (r=0.18), lower milk fat and protein content (r=0.38 and 0.33, respectively), lower quantity of milk produced from forages (r=0.12) and suboptimal reproduction and health performance (e.g. longer calving interval (r=-0.21) and higher culling rate (r=-0.15)). Concerning economic results, low HYMIR-CH4 was significantly associated with lower gross margin per cow (r=0.19) and per litre FPCM (r=0.09). To conclude, this study suggested that low lactating dairy cow gastro-enteric CH4 production tended to be associated with more extensive or suboptimal management practices, which could lead to lower profitability. The observed low correlations suggest complex interactions between variables due to the use of on-farm data with large variability in technical and management practices.

Milk fatty acid composition as an indicator of energy status in Holstein dairy cows

Abstract. Transition dairy cows often enter a stage of negative energy balance during which the utilization of energy reserves is reflected in the milk fatty acid (FA) composition. In this study, metabolic status was evaluated by measuring milk FA, ruminal short-chain FA (RSCFA), and serum biochemical parameters in Holstein cows. Samples (milk, rumen contents, and blood) were collected around days 30 (early) and 150 (middle) of lactation, and rumen contents and blood samples were collected 30 days before calving (dry). Fatty acids were extracted and FA composition was determined. Glucose, triacylglycerols, total cholesterol, low-density lipoprotein cholesterol (LDL-C) and high-density lipoprotein cholesterol (HDL-C), beta-hydroxybutyric acid (BHB), and non-esterified fatty acid (NEFA) concentrations were determined in serum samples. Lower percentages of saturated FA in milk and higher percentages of monounsaturated FA, polyunsaturated FA, and C18:1n-9 were observed in early lactation compared to mid-lactation. In rumen higher concentrations of propionic acid were determined at mid-lactation compared to early lactation. Acetic and butyric acid concentrations showed no significant differences between sampling intervals. In serum higher glucose concentrations were observed during the dry period and mid-lactation than during early lactation. Lower BHB and higher NEFA concentrations were noted during early lactation compared to mid-lactation and the dry period. Total cholesterol, LDL-C, HDL-C and triacylglycerols showed no significant differences between sampling intervals. The results of the present study suggest that determination of milk FA is a potential indicator of energy status in dairy cows.

Genetic and nongenetic factors associated with milk color in dairy cows

Milk color is one of the sensory properties that can influence consumer choice of one product over another and it influences the quality of processed dairy products. This study aims to quantify the cow-level genetic and nongenetic factors associated with bovine milk color traits. A total of 136,807 spectra from Irish commercial and research herds (with multiple breeds and crosses) were used. Milk lightness (Lˆ*), red-green index (aˆ*) and yellow-blue index (bˆ*) were predicted for individual milk samples using only the mid-infrared spectrum of the milk sample. Factors associated with milk color were breed, stage of lactation, parity, milking-time, udder health status, pasture grazing, and seasonal calving. (Co)variance components for Lˆ*,aˆ*, and bˆ* were estimated using random regressions on the additive genetic and within-lactation permanent environmental effects. Greater bˆ* value (i.e., more yellow color) was evident in milk from Jersey cows. Milk Lˆ* increased consistently with stage of lactation, whereas aˆ* increased until mid lactation to subsequently plateau. Milk bˆ* deteriorated until 31 to 60 DIM, but then improved thereafter until the end of lactation. Relative to multiparous cows, milk yielded by primiparae was, on average, lighter (i.e., greater Lˆ*), more red (i.e., greater aˆ*), and less yellow (i.e., lower bˆ*). Milk from the morning milk session had lower Lˆ*,aˆ*, and bˆ* Heritability estimates (±SE) for milk color varied between 0.15 ± 0.02 (30 DIM) and 0.46 ± 0.02 (210 DIM) for Lˆ*, between 0.09 ± 0.01 (30 DIM) and 0.15 ± 0.02 (305 DIM) for aˆ*, and between 0.18 ± 0.02 (21 DIM) and 0.56 ± 0.03 (305 DIM) for bˆ* For all the 3 milk color features, the within-trait genetic correlations approached unity as the time intervals compared shortened and were generally <0.40 between the peripheries of the lactation. Strong positive genetic correlations existed between bˆ* value and milk fat concentration, ranging from 0.82 ± 0.19 at 5 DIM to 0.96 ± 0.01 at 305 DIM and confirming the observed phenotypic correlation (0.64, SE = 0.01). Results of the present study suggest that breeding strategies for the enhancement of milk color traits could be implemented for dairy cattle populations. Such strategies, coupled with the knowledge of milk color traits variation due to nongenetic factors, may represent a tool for the dairy processors to reduce, if not eliminate, the use of artificial pigments during milk manufacturing.

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.

Detecting concentrations of milk components by an iterative optical technique

This paper introduces a theoretical and practical model for reconstructing the scattering properties of a participating media. Our theory is based on a robust generalization of the Gerchberg-Saxton (G-S) algorithm. At the end of this algorithm the reduced scattering coefficient μ's of a given substance, can be estimated from the standard deviation (STD) of the retrieved phase of the remitted light. We use the theory to compute the phase's STD that directly correlated to the optical properties for different types of milk components, and we derive a novel appearance model for milk parameterized by the lactose and protein contents. Our results show that we are able to detect the possibility of lactose and milk proteins' quantitative signature by the G-S optical tool, en route to the design of a novel milk-content-monitoring tool. Sketch of the experimental setup for light intensity measurements and reduced scattering coefficient reconstruction. The samples were prepared from various milk components: whey protein, sodium casienate and lactose, at different concentrations.

Overlap in genomic variation associated with milk fat composition in Holstein Friesian and Dutch native dual-purpose breeds

The aim of this study was to identify if genomic variations associated with fatty acid (FA) composition are similar between the Holstein-Friesian (HF) and native dual-purpose breeds used in the Dutch dairy industry. Phenotypic and genotypic information were available for the breeds Meuse-Rhine-Yssel (MRY), Dutch Friesian (DF), Groningen White Headed (GWH), and HF. First, the reliability of genomic breeding values of the native Dutch dual-purpose cattle breeds MRY, DF, and GWH was evaluated using single nucleotide polymorphism (SNP) effects estimated in HF, including all SNP or subsets with stronger associations in HF. Second, the genomic variation of the regions associated with FA composition in HF (regions on Bos taurus autosome 5, 14, and 26), were studied in the different breeds. Finally, similarities in genotype and allele frequencies between MRY, DF, GWH, and HF breeds were assessed for specific regions associated with FA composition. On average across the traits, the highest reliabilities of genomic prediction were estimated for GWH (0.158) and DF (0.116) when the 8 to 22 SNP with the strongest association in HF were included. With the same set of SNP, GEBV for MRY were the least reliable (0.022). This indicates that on average only 2 (MRY) to 16% (GWH) of the genomic variation in HF is shared with the native Dutch dual-purpose breeds. The comparison of predicted variances of different regions associated with milk and milk fat composition showed that breeds clearly differed in genomic variation within these regions. Finally, the correlations of allele frequencies between breeds across the 8 to 22 SNP with the strongest association in HF were around 0.8 between the Dutch native dual-purpose breeds, whereas the correlations between the native breeds and HF were clearly lower and around 0.5. There was no consistent relationship between the reliabilities of genomic prediction for a specific breed and the correlation between the allele frequencies of this breed and HF. In conclusion, most of the genomic variation associated with FA composition in the Dutch dual-purpose breeds appears to be breed-specific. Furthermore, the minor allele frequencies of genes having an effect on the milk FA composition in HF were shown to be much smaller in the breeds MRY, DF, and GWH, especially for the MRY breed.

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.