Fatty Acid Traits Research Articles

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

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

Co-Expression Network Analysis Identifies miRNA⁻mRNA Networks Potentially Regulating Milk Traits and Blood Metabolites.

MicroRNAs (miRNA) regulate mRNA networks to coordinate cellular functions. In this study, we constructed gene co-expression networks to detect miRNA modules (clusters of miRNAs with similar expression patterns) and miRNA–mRNA pairs associated with blood (triacylglyceride and nonesterified fatty acids) and milk (milk yield, fat, protein, and lactose) components and milk fatty acid traits following dietary supplementation of cows’ diets with 5% linseed oil (LSO) (n = 6 cows) or 5% safflower oil (SFO) (n = 6 cows) for 28 days. Using miRNA transcriptome data from mammary tissues of cows for co-expression network analysis, we identified three consensus modules: blue, brown, and turquoise, composed of 70, 34, and 86 miRNA members, respectively. The hub miRNAs (miRNAs with the most connections with other miRNAs) were miR-30d, miR-484 and miR-16b for blue, brown, and turquoise modules, respectively. Cell cycle arrest, and p53 signaling and transforming growth factor–beta (TGF-β) signaling pathways were the common gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enriched for target genes of the three modules. Protein percent (p = 0.03) correlated with the turquoise module in LSO treatment while protein yield (p = 0.003) and milk yield (p = 7 × 10−04) correlated with the turquoise model, protein and milk yields and lactose percent (p < 0.05) correlated with the blue module and fat percent (p = 0.04) correlated with the brown module in SFO treatment. Several fatty acids correlated (p < 0.05) with the blue (CLA:9,11) and brown (C4:0, C12:0, C22:0, C18:1n9c and CLA:10,12) modules in LSO treatment and with the turquoise (C14:0, C18:3n3 and CLA:9,11), blue (C14:0 and C23:0) and brown (C6:0, C16:0, C22:0, C22:6n3 and CLA:10,12) modules in SFO treatment. Correlation of miRNA and mRNA data from the same animals identified the following miRNA–mRNA pairs: miR-183/RHBDD2 (p = 0.003), miR-484/EIF1AD (p = 0.011) and miR-130a/SBSPON (p = 0.004) with lowest p-values for the blue, brown, and turquoise modules, respectively. Milk yield, protein yield, and protein percentage correlated (p < 0.05) with 28, 31 and 5 miRNA–mRNA pairs, respectively. Our results suggest that, the blue, brown, and turquoise modules miRNAs, hub miRNAs, miRNA–mRNA networks, cell cycle arrest GO term, p53 signaling and TGF-β signaling pathways have considerable influence on milk and blood phenotypes following dietary supplementation of dairy cows’ diets with 5% LSO or 5% SFO.

Preliminary study of solute carrier family 23 member 3 (SLC23A3) gene as candidate marker for fatty acid traits in Kampung-Broiler crossbred chickens

SLC23A3 is one of the key genes which control the properties of the of fatty acids content in the meat. The aim of this study was to identify the association of SNP and mRNA expression of SLC23A3 as a candidate marker for fatty acid traits (FAs). Sixty two F2 Kampung × Broilers crossbred chickens population were used in this study. FA composition was measured at 12 weeks of age from thigh meat of crossbred chicken. The PCR-RFLP and qRT-PCR were used for association study and expression analysis of SLC23A3 gene with divergent fatty acid composition. A SNP in coding region c.22385690 A&gt;C of the SLC23A3 gene was associated (P&lt;0.05) with fatty acid composition including stearic acid (C18:0), elaidic acid (C18:1n9t), and linoleic acid (C18:2n6c). The SLC23A3 was detected in liver from high fatty acids (HFA) and low fatty acid (LFA) composition. However, gene expression of SLC23A3 were not differentially expressed between HFA and LFA. These results will explain better understanding of the key important role of the SLC23A3 in fatty acid traits within the liver and will propose SLC23A3 as a potential genomic selection for selection of chickens with fatty acid composition.

Genetic correlations of mid-infrared-predicted milk fatty acid groups with milk production traits

The objective of this research was to estimate the genetic correlations between milk mid-infrared-predicted fatty acid groups and production traits in first-parity Canadian Holsteins. Contents of short-chain, medium-chain, long-chain, saturated, and unsaturated fatty acid groupings in milk samples can be predicted using mid-infrared spectral data for cows enrolled in milk recording programs. Predicted fatty acid group contents were obtained for 49,127 test-day milk samples from 10,029 first-parity Holstein cows in 810 herds. Milk yield, fat and protein yield, fat and protein percentage, fat-to-protein ratio, and somatic cell score were also available for these test days. Genetic parameters were estimated for the fatty acid groups and production traits using multiple-trait random regression test day models by Bayesian methods via Gibbs sampling. Three separate 8- or 9-trait analyses were performed, including the 5 fatty acid groups with different combinations of the production traits. Posterior standard deviations ranged from <0.001 to 0.01. Average daily genetic correlations were negative and similar to each other for the fatty acid groups with milk yield (-0.62 to -0.59) and with protein yield (-0.32 to -0.25). Weak and positive average daily genetic correlations were found between somatic cell score and the fatty acid groups (from 0.25 to 0.36). Stronger genetic correlations with fat yield, fat and protein percentage, and fat-to-protein ratio were found with medium-chain and saturated fatty acid groups compared with those with long-chain and unsaturated fatty acid groups. Genetic correlations were very strong between the fatty acid groups and fat percentage, ranging between 0.88 for unsaturated and 0.99 for saturated fatty acids. Daily genetic correlations from 5 to 305 d in milk with milk, protein yield and percentage, and somatic cell score traits showed similar patterns for all fatty acid groups. The daily genetic correlations with fat yield at the beginning of lactation were decreasing for long-chain and unsaturated fatty acid groups and increasing for short-chain fatty acids. Genetic correlations between fat percentage and fatty acids were increasing at the beginning of lactation for short- and medium-chain and saturated fatty acids, but slightly decreasing for long-chain and unsaturated fatty acid groups. These results can be used in defining fatty acid traits and breeding objectives.

Genetic analysis of predicted fatty acid profiles of milk from Danish Holstein and Danish Jersey cattle populations

The objective of this study was to assess the genetic variability of the detailed fatty acid (FA) profiles of Danish Holstein (DH) and Danish Jersey (DJ) cattle populations. We estimated genetic parameters for 11 FA or groups of FA in milk samples from the Danish milk control system between May 2015 and October 2016. Concentrations of different FA and FA groups in milk samples were measured by mid-infrared spectroscopy. Data used for parameter estimation were from 132,732 first-parity DH cows and 21,966 first-parity DJ cows. We found the highest heritabilities for test day measurements in both populations for short-chain FA (DH = 0.16; DJ = 0.16) and C16:0 (DH = 0.14; DJ = 0.16). In DH, the highest heritabilities were also found for saturated FA and monounsaturated FA (both populations: 0.15). Genetic correlations between the fatty acid traits showed large differences between DH and DJ for especially short-chain FA with the other FA traits measured. Furthermore, genetic correlations of total fat with monounsaturated FA, polyunsaturated FA, short-chain FA, and C16:0 differed markedly between DH and DJ populations. In conclusion, we found genetic variation in the mid-infrared spectroscopy-predicted FA and FA groups of the DH and DJ cattle populations. This finding opens the possibility of using genetic selection to change the FA profiles of dairy cattle.

Expression and Association of SCD Gene Polymorphisms and Fatty Acid Compositions in Chicken Cross

Stearoyl-CoA desaturase (SCD) is an integral membrane protein of endoplasmic reticulum (ER) that catalyzes the rate limiting step in the monounsaturated fatty acids from saturated fatty acids. Selection for fatty acids traits based on molecular marker assisted selection is needed to increase a value of chicken meat. This study was designed to analyze expression and associations of SCD gene polymorphisms with fatty acid traits in F2 kampung-broiler chicken cross. A total of 62 F2 kampung-broiler chicken cross (29 males and 33 females) were used in this study. Fatty acid traits were measured at 26 weeks of age. Samples were divided into two groups based on fatty acid traits (the highest and the lowest). Primers in exon 2 region were designed from the genomic chicken sequence. The SNP g.37284A>G was detected and polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method was then used to genotype. The expression of SCD gene was analyzed using quantitative real time PCR (qRT-PCR). The result showed that there were three genotypes (AA, AG, and GG) found in this study. The SCD|AciI polymorphism was significantly associated with palmitoleic acid (C16:1), fatty acids total and saturated fatty acid in 26 weeks old of F2 kampung-broiler chicken cross (P<0.05). The SCD gene was expressed for polyunsaturated fatty acids in liver tissue in two groups of chickens. In conclusion, the SCD gene could be a candidate gene that affects fatty acids traits in F2 kampung-broiler chicken cross.

Phenotypic and genetic relationships of residual feed intake measures and their component traits with fatty acid composition in subcutaneous adipose of beef cattle.

Feed efficiency is of particular interest to the beef industry because feed is the largest variable cost in production and fatty acid composition is emerging as an important trait, both economically and socially, due to the potential implications of dietary fatty acids on human health. Quantifying correlations between feed efficiency and fatty acid composition will contribute to construction of optimal multiple-trait selection indexes to maximize beef production profitability. In the present study, we estimated phenotypic and genetic correlations of feed efficiency measures including residual feed intake (RFI), RFI adjusted for final ultrasound backfat thickness (RFIf); their component traits ADG, DMI, and metabolic BW; and final ultrasound backfat thickness measured at the end of feedlot test with 25 major fatty acids in the subcutaneous adipose tissues of 1,366 finishing steers and heifers using bivariate animal models. The phenotypic correlations of RFI and RFIf with the 25 individual and grouped fatty acid traits were generally low (<0.25 in magnitude). However, relatively stronger genetic correlation coefficients of RFI and RFIf with PUFA traits including the -6:-3 ratio (0.52 ± 0.29 and 0.45 ± 0.31, respectively), 18:2-6 (0.45 ± 0.18 and 0.40 ± 0.19, respectively), -6 (0.43 ± 0.18 and 0.38 ± 0.19, respectively), PUFA (0.42 ± 0.18 and 0.36 ± 0.20, respectively), and 9-16:1 (-0.43 ± 0.20 and -0.33 ± 0.22, respectively) were observed. Hence, selection for low-RFI or more efficient beef cattle will improve fatty acid profiles by lowering the content of -6 PUFA, thus reducing the ratio of -6 to -3 along with increasing the amount of 9-16:1. Moderate to moderately high genetic correlations were also observed for DMI with 9-14:1 (-0.32 ± 0.17) and the sum of CLA analyzed (SumCLA; -0.45 ± 0.21), suggesting that selection of beef cattle with lower DMI will lead to an increase amount of 9-14:1 and SumCLA in the subcutaneous adipose tissue. However, unfavorable genetic correlations were detected for ADG with 11-18:1 (-0.38 ± 0.23) and SumCLA (-0.73 ± 0.26), implying that selection of beef cattle with a better growth rate will decrease the contents of healthy fatty acids 11-18:1 and SumCLA. Therefore, it is recommended that a multiple-trait selection index be used when genetic improvements of fatty acid composition, feed efficiency, feed intake, and growth are important in the breeding objective.

Genome-wide association study confirm major QTL for backfat fatty acid composition on SSC14 in Duroc pigs

BackgroundFatty acid composition contributes importantly to meat quality and is essential to the nutritional value of the meat. Identification of genetic factors underlying levels of fatty acids can be used to breed for pigs with healthier meat. The aim of this study was to conduct genome-wide association studies (GWAS) to identify QTL regions affecting fatty acid composition in backfat from the pig breeds Duroc and Landrace.ResultsUsing data from the Axiom porcine 660 K array, we performed GWAS on 454 Duroc and 659 Landrace boars for fatty acid phenotypes measured by near-infrared spectroscopy (NIRS) technology (C16:0, C16:1n-7, C18:0, C18:1n-9, C18:2n-6, C18:3n-3, total saturated fatty acids, monounsaturated fatty acids and polyunsaturated fatty acids). Two QTL regions on SSC4 and SSC14 were identified in Duroc for the de novo synthesized fatty acids traits, whereas one QTL on SSC8 was detected in Landrace for C16:1n-7. The QTL region on SSC14 has been reported in previous studies and a putative causative mutation has been suggested in the promoter region of the SCD gene. Whole genome re-sequencing data was used for genotype imputation and to fine map the SSC14 QTL region in Norwegian Duroc. This effort confirms the location of the QTL on this chromosome as well as suggesting other putative candidate genes in the region. The most significant single nucleotide polymorphisms (SNPs) located on SSC14 explain between 55 and 76% of the genetic variance and between 27 and 54% of the phenotypic variance for the de novo synthesized fatty acid traits in Norwegian Duroc. For the QTL region on SSC8 in Landrace, the most significant SNP explained 19% of the genetic variance and 5% of the phenotypic variance for C16:1n-7.ConclusionsThis study confirms a major QTL affecting fatty acid composition on SSC14 in Duroc, which can be used in genetic selection to increase the level of fatty acid desaturation. The SSC14 QTL was not segregating in the Landrace population, but another QTL on SSC8 affecting C16:1n-7 was identified and might be used to increase the level of desaturation in meat products from this breed.

Cloning and expression of Perilla frutescens FAD2 gene and polymorphism analysis among cultivars

∆12 fatty acid desaturase (FAD2) is a key enzyme for linoleic acid and linolenic acid biosynthesis. Perilla frutescens is a special oil plant species with highest linolenic acid content. In this study, based on RACE, two alleles for one FAD2 gene were isolated from P. frutescens cultivar C2: the 3956 bp PfFAD2a and the 3959 bp PfFAD2b, both with a full-length cDNA of 1526 bp, and both encoding a 382aa basic protein. The alleles have identities of over 98%, and their encoded proteins differ only by substitution of a strongly similar residue. Saccharomyces cerevisiae heterologous expression suggested that PfFAD2a/b both encode a bio-functional FAD2 enzyme. Phylogenetic analyses indicated that PfFAD2 shows the highest homologies to FAD2 genes from dicots such as Boraginaceae and Burseraceae. PfFAD2a/b expressions are mainly restricted to developing seeds. PfFAD2a/b expression in the seedling leaf is upregulated by cold (4 °C) and repressed by heat (42 °C). Each of the eight cultivars contains two alleles for one PfFAD2 and 40 SNP sites are found. One allelic gene in cultivars C1 and P1 is pseudogene because of premature stop codon mutation in 5′ coding region. All other normal PfFAD2 genes/allelic genes encode identical or very similar proteins. PfFAD2a/b expression level in developing seeds also varies among the eight cultivars. This study provides systemic molecular and functional features of PfFAD2 and enables its application in the study of plant fatty acids traits.

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

BackgroundBovine 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).ResultsThe 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.ConclusionsAn 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).

Factors affecting variations in the detailed fatty acid profile of Mediterranean buffalo milk determined by 2-dimensional gas chromatography

Buffalo milk is the world's second most widely produced milk, and increasing attention is being paid to its composition, particularly the fatty acid profile. The objectives of the present study were (1) to characterize the fatty acid composition of Mediterranean buffalo milk, and (2) to investigate potential sources of variation in the buffalo milk fatty acid profile. We determined the profile of 69 fatty acid traits in 272 individual samples of Mediterranean buffalo milk using gas chromatography. In total, 51 individual fatty acids were identified: 24 saturated fatty acids, 13 monounsaturated fatty acids, and 14 polyunsaturated fatty acids. The major individual fatty acids in buffalo milk were in the order 16:0, 18:1 cis-9, 14:0, and 18:0. Saturated fatty acids were the predominant fraction in buffalo milk fat (70.49%); monounsaturated and polyunsaturated fatty acids were at 25.95 and 3.54%, respectively. Adopting a classification based on carbon-chain length, we found that medium-chain fatty acids (11-16 carbons) represented the greater part (53.7%) of the fatty acid fraction of buffalo milk, whereas long-chain fatty acids (17-24 carbons) and short-chain fatty acids (4-10 carbons) accounted for 32.73 and 9.72%, respectively. The n-3 and n-6 fatty acids were 0.46 and 1.77%, respectively. The main conjugated linoleic acid, rumenic acid, represented 0.45% of total milk fatty acids. Herd/test date and stage of lactation were confirmed as important sources of variation in the fatty acid profile of buffalo milk. The percentages of short-chain and medium-chain fatty acids in buffalo milk increased in early lactation (+0.6 and +3.5%, respectively), whereas long-chain fatty acids decreased (-4.2%). The only exception to this pattern was butyric acid, which linearly decreased from the beginning of lactation, confirmation that its synthesis is independent of malonyl-CoA. These results seem to suggest that in early lactation the mobilization of energy reserves may have less influence on the fatty acid profile of buffalo milk than that of cow milk, probably due to a shorter and less severe period of negative energy balance. Parity affected the profiles of a few traits and had the most significant effects on branched-chain fatty acids. This work provided a detailed overview of the fatty acid profile in buffalo milk including also those fatty acids present in small concentrations, which may have beneficial effects for human health. Our results contributed also to increase the knowledge about the effects of some of the major factors affecting buffalo production traits and fatty acid concentrations in milk, and consequently its technological and nutritional properties.

Transcriptomic Changes in Liver of Young Bulls Caused by Diets Low in Mineral and Protein Contents and Supplemented with n-3 Fatty Acids and Conjugated Linoleic Acid.

The aim of the present study was to identify transcriptional modifications and regulatory networks accounting for physiological and metabolic responses to specific nutrients in the liver of young Belgian Blue × Holstein bulls using RNA-sequencing. A larger trial has been carried out in which animals were fed with different diets: 1] a conventional diet; 2] a low-protein/low-mineral diet (low-impact diet) and 3] a diet enriched in n-3 fatty acids (FAs), conjugated linoleic acid (CLA) and vitamin E (nutraceutical diet). The initial hypothesis was that the administration of low-impact and nutraceutical diets might influence the transcriptional profiles in bovine liver and the resultant nutrient fluxes, which are essential for optimal liver function and nutrient interconversion. Results showed that the nutraceutical diet significantly reduced subcutaneous fat covering in vivo and liver pH. Dietary treatments did not affect overall liver fat content, but significantly modified the liver profile of 33 FA traits (out of the total 89 identified by gas-chromatography). In bulls fed nutraceutical diet, the percentage of n-3 and CLA FAs increased around 2.5-fold compared with the other diets, whereas the ratio of n6/n3 decreased 2.5-fold. Liver transcriptomic analyses revealed a total of 198 differentially expressed genes (DEGs) when comparing low-impact, nutraceutical and conventional diets, with the nutraceutical diet showing the greatest effects on liver transcriptome. Functional analyses using ClueGo and Ingenuity Pathway Analysis evidenced that DEGs in bovine liver were variously involved in energy reserve metabolic process, glutathione metabolism, and carbohydrate and lipid metabolism. Modifications in feeding strategies affected key transcription factors regulating the expression of several genes involved in fatty acid metabolism, e.g. insulin-induced gene 1, insulin receptor substrate 2, and RAR-related orphan receptor C. This study provides noteworthy insights into the molecular changes occurring as a result of nutrient variation in diets (aimed at reducing the environmental impact and improving human health) and broadens our understanding of the relationship between nutrients variation and phenotypic effects.

Epistatic interactions associated with fatty acid concentrations of beef from angus sired beef cattle.

BackgroundConsumers are becoming increasingly conscientious about the nutritional value of their food. Consumption of some fatty acids has been associated with human health traits such as blood pressure and cardiovascular disease. Therefore, it is important to investigate genetic variation in content of fatty acids present in meat. Previously publications reported regions of the cattle genome that are additively associated with variation in fatty acid content. This study evaluated epistatic interactions, which could account for additional genetic variation in fatty acid content.ResultsEpistatic interactions for 44 fatty acid traits in a population of Angus beef cattle were evaluated with EpiSNPmpi. False discovery rate (FDR) was controlled at 5 % and was limited to well-represented genotypic combinations. Epistatic interactions were detected for 37 triacylglyceride (TAG), 36 phospholipid (PL) fatty acid traits, and three weight traits. A total of 6,181, 7,168, and 0 significant epistatic interactions (FDR < 0.05, 50-animals per genotype combination) were associated with Triacylglyceride fatty acids, Phospholipid fatty acids, and weight traits respectively and most were additive-by-additive interactions. A large number of interactions occurred in potential regions of regulatory control along the chromosomes where genes related to fatty acid metabolism reside.ConclusionsMany fatty acids were associated with epistatic interactions. Despite a large number of significant interactions, there are a limited number of genomic locations that harbored these interactions. While larger population sizes are needed to accurately validate and quantify these epistatic interactions, the current findings point towards additional genetic variance that can be accounted for within these fatty acid traits.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-3235-8) contains supplementary material, which is available to authorized users.

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.

Investigating mutual relationship among milk fatty acids by multivariate factor analysis in dairy cows

The interpretation of mutual relationship among milk fatty acids (FA) is not easy due to the high number of FA contained in milk fat and to the high degree of correlation among them. Multivariate analysis includes different statistical approaches that could help explaining complex pattern of variables. In this study, Multivariate Factor Analysis (MFA) was used to decompose the correlation matrix of 47 FA and milk production traits (milk yield and protein and fat content) measured in 300 Italian Holstein Friesian cows reared in the North of Italy in 23 commercial dairy farms, representative of the intensive dairy system. MFA was able to extract seven latent factors with specific biologic meaning: secretion of Long Chain FA (KLCFA), mammary FA de novo synthesis (Km), rumen biohydrogenation (Kbh), synthesis of odd chain FA (Ko), synthesis of branched chain FA (Kb), mammary desaturation activity (Kd), milk yield (Kmy). According to the pattern of communalities of the factor analysis, C18:3c9c12c15 was the only FA, along with C18:2t11c15, to be uncorrelated with the other variables and it seemed to be excluded by the metabolic pattern described by the seven factors. The desaturation products of the SCD enzyme were independently associated to three latent factors, suggesting new insights in the regulation of SCD activity. Factors were considered as new quantitative phenotypes related to prominent features of milk FA profile. With the aim of evaluating the feeding regimen and animal effects, latent factors were analysed with a mixed model, which considered the fixed effect of lactation stage, parity, some feeding regimen characteristics and the random effect of bull. Lactation stage significantly affected Km and Kmy factors. In perspective, the seven factors extracted by applying MFA analysis to milk FA composition could be considered as new and more informative traits to test the effect of endogenous and exogenous variation factors.

Effects of candidate gene polymorphisms on the detailed fatty acids profile determined by gas chromatography in bovine milk

Association analyses between candidate genes and bovine milk fatty acids can improve our understanding of genetic variation in milk fatty acid profiles and reveal potential opportunities to tailor milk fat composition through selection strategies. In this work, we investigated the association of 51 single nucleotide polymorphisms (SNP) selected from 37 candidate genes using a functional and positional approach, with 47 fatty acids, 9 fatty acid groups, and 5 Δ9-desaturation indices in milk samples from Brown Swiss cows. Individual milk samples were collected from 1,158 Italian Brown Swiss cows, and gas chromatography was used to obtain detailed milk fatty acid compositions. A GoldenGate assay system (Illumina, San Diego, CA) was used to perform genotype 96 selected SNP located in 54 genes across 22 chromosomes. In total, 51 polymorphic SNP in 37 candidate genes were retained for the association analysis. A Bayesian linear animal model was used to estimate the contribution of each SNP. A total of 129 tests indicated relevant additive effects between a given SNP and a single fatty acid trait; 38 SNP belonging to 30 genes were relevant for a total of 57 fatty acid traits. Most of the studied fatty acid traits (~81%) were relevantly associated with multiple SNP. Relevantly associated SNP were mainly found in genes related to fat metabolism, linked to or contained in previously identified quantitative trait loci for fat yield or content, or associated with genes previously identified in association analyses with milk fatty acid profiles in other cow breeds. The most representative candidate genes were LEP, PRL, STAT5A, CCL3, ACACA, GHR, ADRB2, LPIN1, STAT1, FABP4, and CSN2. In particular, relevant associations with SNP located on bovine chromosome 19 (BTA19) were found. Two candidate genes on BTA19 (CCL3 and ACACA) were relevantly associated with de novo short- and medium-chain fatty acids, likely explaining the high heritability values found for these fatty acids (with the exception of C6:0). Two additional genes on BTA19 (CCL2 and GH1) showed associations with saturated and branched-chain fatty acids. Our findings provide basic information on genes and SNP affecting the milk fatty acid composition of dairy cows. These results may support the possibility of using genetic selection to modify milk fatty acid profiles to promote beneficial health-related effects.

On the role of mid-infrared predicted phenotypes in fertility and health dairy breeding programs

Fertility and health traits are of prime importance in dairy breeding programs. However, these traits are generally complex, difficult to record, and lowly heritable (<0.10), thereby hampering genetic improvement in disease resistance and fertility. Hence, indicators are useful in the prediction of genetic merit for fertility and health traits as long as they are easier to measure than direct fitness traits, heritable, and genetically correlated. Considering that changes in (fine) milk composition over a lactation reflect the physiological status of the cow, mid-infrared (MIR) analysis of milk opens the door to a wide range of potential indicator traits of fertility and health. Previous studies investigated the phenotypic and genetic relationships between fertility and MIR-predicted phenotypes, most being related to negative postpartum energy balance and body fat mobilization (e.g., fat:protein ratio, urea, fatty acids profile). Results showed that a combination of various fatty acid traits (e.g., C18:1 cis-9 and C10:0) could be used to improve fertility. Furthermore, occurrence of (sub)clinical ketosis has been related to milk-based phenotypes such as fat:protein ratio, fatty acids, and ketone bodies. Hence, MIR-predicted acetone and β-hydroxybutyrate contents in milk could be useful for breeding cows less susceptible to ketosis. Although studies investigating the genetic association among mastitis and MIR-predicted phenotypes are scarce, a wide range of traits, potentially predicted by MIR spectrometry, are worthy of consideration. These include traits related to the disease response of the cow (e.g., lactoferrin), reduced secretory activity (e.g., casein), and the alteration of the blood-milk barrier (e.g., minerals). Moreover, direct MIR prediction of fertility and health traits should be further considered. To conclude, MIR-predicted phenotypes have a role to play in the improvement of dairy cow fertility and health. However, further studies are warranted to (1) grasp underlying associations among MIR-predicted indicator and fitness traits, (2) estimate the genetic parameters, and (3) include these traits in broader breeding strategies.

Genome-wide association and genomic prediction of breeding values for fatty acid composition in subcutaneous adipose and longissimus lumborum muscle of beef cattle.

BackgroundIdentification of genetic variants that are associated with fatty acid composition in beef will enhance our understanding of host genetic influence on the trait and also allow for more effective improvement of beef fatty acid profiles through genomic selection and marker-assisted diet management. In this study, 81 and 83 fatty acid traits were measured in subcutaneous adipose (SQ) and longissimus lumborum muscle (LL), respectively, from 1366 purebred and crossbred beef steers and heifers that were genotyped on the Illumina BovineSNP50 Beadchip. The objective was to conduct genome-wide association studies (GWAS) for the fatty acid traits and to evaluate the accuracy of genomic prediction for fatty acid composition using genomic best linear unbiased prediction (GBLUP) and Bayesian methods.ResultsIn total, 302 and 360 significant SNPs spanning all autosomal chromosomes were identified to be associated with fatty acid composition in SQ and LL tissues, respectively. Proportions of total genetic variance explained by individual significant SNPs ranged from 0.03 to 11.06 % in SQ, and from 0.005 to 24.28 % in the LL muscle. Markers with relatively large effects were located near fatty acid synthase (FASN), stearoyl-CoA desaturase (SCD), and thyroid hormone responsive (THRSP) genes. For the majority of the fatty acid traits studied, the accuracy of genomic prediction was relatively low (<0.40). Relatively high accuracies (> = 0.50) were achieved for 10:0, 12:0, 14:0, 15:0, 16:0, 9c-14:1, 12c-16:1, 13c-18:1, and health index (HI) in LL, and for 12:0, 14:0, 15:0, 10 t,12c-18:2, and 11 t,13c + 11c,13 t-18:2 in SQ. The Bayesian method performed similarly as GBLUP for most of the traits but substantially better for traits that were affected by SNPs of large effects as identified by GWAS.ConclusionsFatty acid composition in beef is influenced by a few host genes with major effects and many genes of smaller effects. With the current training population size and marker density, genomic prediction has the potential to predict the breeding values of fatty acid composition in beef cattle at a moderate to relatively high accuracy for fatty acids that have moderate to high heritability.Electronic supplementary materialThe online version of this article (doi:10.1186/s12863-015-0290-0) contains supplementary material, which is available to authorized users.

Bayesian regression models outperform partial least squares methods for predicting milk components and technological properties using infrared spectral data

The aim of this study was to assess the performance of Bayesian models commonly used for genomic selection to predict “difficult-to-predict” dairy traits, such as milk fatty acid (FA) expressed as percentage of total fatty acids, and technological properties, such as fresh cheese yield and protein recovery, using Fourier-transform infrared (FTIR) spectral data. Our main hypothesis was that Bayesian models that can estimate shrinkage and perform variable selection may improve our ability to predict FA traits and technological traits above and beyond what can be achieved using the current calibration models (e.g., partial least squares, PLS). To this end, we assessed a series of Bayesian methods and compared their prediction performance with that of PLS. The comparison between models was done using the same sets of data (i.e., same samples, same variability, same spectral treatment) for each trait. Data consisted of 1,264 individual milk samples collected from Brown Swiss cows for which gas chromatographic FA composition, milk coagulation properties, and cheese-yield traits were available. For each sample, 2 spectra in the infrared region from 5,011 to 925cm−1 were available and averaged before data analysis. Three Bayesian models: Bayesian ridge regression (Bayes RR), Bayes A, and Bayes B, and 2 reference models: PLS and modified PLS (MPLS) procedures, were used to calibrate equations for each of the traits. The Bayesian models used were implemented in the R package BGLR (http://cran.r-project.org/web/packages/BGLR/index.html), whereas the PLS and MPLS were those implemented in the WinISI II software (Infrasoft International LLC, State College, PA). Prediction accuracy was estimated for each trait and model using 25 replicates of a training-testing validation procedure. Compared with PLS, which is currently the most widely used calibration method, MPLS and the 3 Bayesian methods showed significantly greater prediction accuracy. Accuracy increased in moving from calibration to external validation methods, and in moving from PLS and MPLS to Bayesian methods, particularly Bayes A and Bayes B. The maximum R2 value of validation was obtained with Bayes B and Bayes A. For the FA, C10:0 (% of each FA on total FA basis) had the highest R2 (0.75, achieved with Bayes A and Bayes B), and among the technological traits, fresh cheese yield R2 of 0.82 (achieved with Bayes B). These 2 methods have proven to be useful instruments in shrinking and selecting very informative wavelengths and inferring the structure and functions of the analyzed traits. We conclude that Bayesian models are powerful tools for deriving calibration equations, and, importantly, these equations can be easily developed using existing open-source software. As part of our study, we provide scripts based on the open source R software BGLR, which can be used to train customized prediction equations for other traits or populations.

Design of a low-density SNP chip for the main Australian sheep breeds and its effect on imputation and genomic prediction accuracy.

Genotyping sheep for genome-wide SNPs at lower density and imputing to a higher density would enable cost-effective implementation of genomic selection, provided imputation was accurate enough. Here, we describe the design of a low-density (12k) SNP chip and evaluate the accuracy of imputation from the 12k SNP genotypes to 50k SNP genotypes in the major Australian sheep breeds. In addition, the impact of imperfect imputation on genomic predictions was evaluated by comparing the accuracy of genomic predictions for 15 novel meat traits including carcass and meat quality and omega fatty acid traits in sheep, from 12k SNP genotypes, imputed 50k SNP genotypes and real 50k SNP genotypes. The 12k chip design included 12223 SNPs with a high minor allele frequency that were selected with intermarker spacing of 50-475kb. SNPs for parentage and horned or polled tests also were represented. Chromosome ends were enriched with SNPs to reduce edge effects on imputation. The imputation performance of the 12k SNP chip was evaluated using 50k SNP genotypes of 4642 animals from six breeds in three different scenarios: (1) within breed, (2) single breed from multibreed reference and (3) multibreed from a single-breed reference. The highest imputation accuracies were found with scenario 2, whereas scenario 3 was the worst, as expected. Using scenario 2, the average imputation accuracy in Border Leicester, Polled Dorset, Merino, White Suffolk and crosses was 0.95, 0.95, 0.92, 0.91 and 0.93 respectively. Imputation scenario 2 was used to impute 50k genotypes for 10396 animals with novel meat trait phenotypes to compare genomic prediction accuracy using genomic best linear unbiased prediction (GBLUP) with real and imputed 50k genotypes. The weighted mean imputation accuracy achieved was 0.92. The average accuracy of genomic estimated breeding values (GEBVs) based on only 12k data was 0.08 across traits and breeds, but accuracies varied widely. The mean GBLUP accuracies with imputed 50k data more than doubled to 0.21. Accuracies of genomic prediction were very similar for imputed and real 50k genotypes. There was no apparent impact on accuracy of GEBVs as a result of using imputed rather than real 50k genotypes, provided imputation accuracy was >90%.