Email updates

Keep up to date with the latest news and content from BMC Genetics and BioMed Central.

Open Access Research article

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

Eveline M Ibeagha-Awemu1*, Kingsley A Akwanji2, Frédéric Beaudoin1 and Xin Zhao2*

Author Affiliations

1 Dairy and Swine Research and Development Centre, Agriculture and Agri-Food Canada, 2000 Rue College, Sherbrooke, Quebec J1M 0C8, Canada

2 Department of Animal Science, McGill University, 21,111 Lakeshore Road, Ste-Anne-De-Bellevue, Quebec H9X 3V9, Canada

For all author emails, please log on.

BMC Genetics 2014, 15:25  doi:10.1186/1471-2156-15-25


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2156/15/25


Received:5 June 2013
Accepted:12 February 2014
Published:17 February 2014

© 2014 Ibeagha-Awemu et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.

Abstract

Background

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

Results

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

Conclusion

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

Keywords:
FADS1; FADS2; SNP; Omega-3 fatty acids; Omega-6 fatty acids; Milk fatty acids; Canadian Holstein cows

Background

The polyunsaturated fatty acid (PUFA) concentrations in blood and tissue lipids are known to be closely related to several positive health outcomes on cardiovascular disease morbidity and mortality, early visual, cognitive and motor development, mental health and psychiatric disorders, and early growth and development during pregnancy and early childhood [1]. In addition, PUFAs exert anticancer effects [2] and play beneficial roles in preventing and/or treating various inflammatory immune disorders such as allergies [3,4] and in male fertility and spermatogenesis [5]. Most of these effects are primarily facilitated by the long chain PUFAs (LC-PUFAs) (have 20 carbon atoms or more) arachidonic acid (AA, 20:4n-6), eicosapentaenoic acid (EPA, 20:5n-3) and docosahexaenoic acid (DHA, 22:6n-3). PUFAs are also known to play crucial roles in various cellular functions, such as acting as components of cellular membranes and as precursors of important mediators of inflammation (e.g. EPA, [6]). PUFAs of 18 carbon atoms are mainly composed of isomers of conjugated linoleic acid (CLA). They too like LC-PUFAs have demonstrated health related benefits such as antiadipogenic, anticarcinogenic, antiatherogenic, antidiabetogenic and antiinflammatory properties (review by [7,8]). LC-PUFAs are classified in three principal families, omega-6 (n-6), omega-3 (n-3) or omega-9 (n-9) fatty acids (FAs), according to the position of the terminal double bond from the methyl end of the FA [9]. The precursor FAs of omega-3 and -6 families are respectively linoleic acid (LA) and alpha linolenic acid (ALA). They cannot be synthesised by mammals so they must be provided by the diet and are therefore defined as essential FAs. Sources of unsaturated fatty acids (USFAs) including monounsaturated fatty acids (MUFAs) and PUFAs in human diets are plant oils, oil seeds, and animal products including fish, meat, eggs and milk. The major source of CLA is from ruminant products, mainly milk and dairy products. Rumenic acid is the major form of CLA in milk and it is endogenously synthesized through the activity of stearoyl-CoA desaturase 1 (SCD1) enzyme or delta 9 desaturase. SCD1 desaturates vaccenic acid, an intermediate in rumen biohydrogenation to form CLA [10]. Mammals or humans can also synthesize LC-PUFAs endogenously starting from the precursor essential FAs, LA and ALA. In western diets, bovine milk is therefore a major contributor of CLA, omega-3 and omega-6 FAs to humans. However, the proportion of these beneficial FAs in bovine milk is generally low.

Endogenous synthesis of omega-3 and -6 FAs is through desaturation of LA and ALA by the activities of fatty acid synthase 1 (FADS1) and 2 (FADS2), and elongases. FADS1 and FADS2 code, respectively, for the enzymes Δ-5 and Δ-6 desaturases considered rate limiting enzymes in the synthesis of LC-PUFAs [11,12]. FADS1 and FADS2 add double bonds at the Δ-5 and Δ-6 position, respectively, of LC-PUFAs. These genes are clustered in a head-on-head direction on bovine chromosome 29 (BTA 29). A third desaturase gene, FADS3 is also clustered at the same location but its role in FA metabolism in cows is not yet clear. The concentration of PUFAs in bovine milk is known to be influenced by dietary factors [13,14] and genetics/breed [15-17]. Furthermore, low to medium heritability values for USFAs recorded for many breeds [18,19] indicate the possibility of improvement through genetic selection. In humans, several studies suggest that plasma and tissue concentrations of omega-3 and -6 FAs are associated with several single nucleotide polymorphisms (SNPs) in the FADS1 and FADS2 genes [20-23]. Also, recent gene and genome-wide association studies (GWAS) in humans have highlighted the influence of variations in the FADS1 and FADS2 gene cluster on lipid metabolism, glucose metabolism, and other quantitative traits such as total cholesterol and low-density lipoprotein [24-26], and disease conditions such as development of metabolic syndrome [27], coronary artery disease [28], myocardial infarction [29], and dyslipidemia [30]. In these studies, changes in the levels of individual FAs related to FADS enzyme activity were noted, thus suggesting a relationship between these conditions and a deregulation in desaturase activity. In a recent study, Merino et al. [31] demonstrated that genetic variants in FADS1 and FADS2 could alter desaturase activities in subjects of Caucasians and Asian descent. Even though several sequence variations including SNPs are listed for these genes for bovine in the SNP database (dbSNP, http://www.ncbi.nlm.nih.gov/projects/SNP/ webcite), no association with milk FAs contents has been attempted. Increasing the proportion of USFAs, especially PUFAs in milk may greatly influence human health and wellness. Dietary manipulations through feeding diets rich in USFAs have proven successful but come with significant extra cost. Exploiting animal’s inherent abilities through searching for FADS gene variants with increased desaturase activities and including them in breeding programs may lead to increased milk PUFA content.

The aim of this study was to characterize the FADS1 and FADS2 genes of Canadian Holstein cows for sequence variations and to determine possible associations with beneficial FAs content in milk.

Methods

Animal ethics

Animal care use and protocols were approved by the McGill University Animal Care Committee.

Screening for polymorphisms in FADS genes

DNA from 40 randomly selected Holstein cows from the Macdonald Campus Dairy Unit, McGill University were used to screen for polymorphisms within the FADS genes. DNA was isolated from leucocytes using a modified method of Montgomery and Sise [32]. To screen for polymorphisms, only the regulatory regions (promoters, 5′ and 3′UTRs) and coding regions (exons including about 200 bp of surrounding intronic sequences) were targeted for selective amplification by PCR and Sanger sequencing. Primers targeting the regions of interest were designed using primer3 program and the reference sequences NC_007330.4 (region 41911653-41925160, FADS1) and NC_007330.4 (region 41998252-42036271, FADS2). After optimization, each PCR reaction was carried out in a final volume of 40 μL. PCR reaction mix contained 30 ng of genomic DNA, 0.25 mM dNTPs, 1.5 to 3.0 mM MgCl2 (Additional file 1: Table S1), 10 μM each primer, 1× Crimpson Taq buffer and 2 units Crimpson Tag DNA polymerase (NEB Canada, Pickering, ON, Canada). PCR cycling conditions with an Eppendorf Mastercycler® pro (Eppendorf North America, Hauppauge, NY, USA) included an initial denaturation step at 94°C for 2 min followed by 30 cycles of 94°C for 30 sec, 56 to 62°C (Additional file 1: Table S1) for 30 sec and 72°C for 1 min, plus a final elongation step of 72°C for 5 min. Successful PCR amplification was confirmed by running 5 μL of PCR products on 1% agarose gels and visualization under UV rays. Thirty μL of PCR products were sent to McGill University/Genome Quebec Innovation Centre for Sanger sequencing using the big dye termination technique and an ABI 3700 sequencer. Sequences were processed for polymorphisms with ChromasPro version 1.7.3 (http://technelysium.com.au/ webcite) and ABI Prism SeqScape® version 2.1 (Life Technologies Inc., Burlington, ON, Canada) software. Effects of regulatory region polymorphisms on microRNA binding sites were searched with TargetScan Release 6.0 (http://www.targetscan.org/vert_60/ webcite).

Additional file 1: Table S1. Primer sequences used in charactering FADS genes by Sanger sequencing. Reference sequences used in primer design with primer3 program were NC_007330.4 (region 41911653-41925160, FADS1) and NC_007330.4 (region 41998252-42036271, FADS2).

Format: DOCX Size: 21KB Download fileOpen Data

Milk sampling and fatty acid analysis

Milk samples were collected from 450 cows from five herds in the province of Quebec. Cows were in mid-lactation and their parities ranged from 1 to 4. Milk sample collection was coordinated by Valacta, the dairy production center of expertise for Quebec and the Atlantic provinces (http://www.valacta.com webcite). Milk samples were centrifuged (12000 × g at 4°C for 30 min) immediately upon arrival at the laboratory to separate the fat from the whey and somatic cell fractions. The fat portion was used for FA analysis while the somatic cells were used for DNA isolation. The different FA isomers were separated by capillary gas chromatography on a Varian CP-3900 gas chromatograph equipped with a Varian CP-8400 auto-sampler and auto-injector, column oven and a flame ionization detector (Varian Inc., Walnut Creek, CA, USA) according to O’Fallon et al. [33].

Genomic DNA purification from milk somatic cells

Genomic DNA was isolated from milk somatic cells using NucleoSpin® Blood QuickPure kit (MJS Biolynx, Ontario, Canada) with some modifications. Briefly, somatic cells were washed twice with phosphate buffered saline (PBS) and resuspended in 200 µL of TE buffer (10 mMTris-HCl and 1 mM EDTA pH 7.6). Three hundred µL of 0.5 M EDTA was added and mixed by shaking vigorously for 45 min followed by centrifugation for 10 min at 3000 × g. The pellet was resuspended in 200 µL of PBS. Twenty five μL proteinase K solution (10 mg/mL) and 200 µL of lysis buffer BQ1 (NucleoSpin® kit) were added to each tube before incubation for 15 min at 70°C. Seven hundred µL of chloroform (Sigma, Ontario, Canada) was added and tubes were vigorously shaken (15 secs). After 10 min of centrifugation at 14000 × g, the aqueous phase was transferred to a 1.5 mL tube containing 210 µL of anhydrous ethanol, mixed briefly and loaded on the NucleoSpin® blood column. The standard protocol of the NucleoSpin® kit was followed from this point except that three washes were done instead of two. Genomic DNA concentration was measured with a NanoDrop® spectrophotometer (NanoDrop Technologies, Inc., Wilmington, DE, USA).

SNP genotyping

Five potentially functional SNPs (includes one synonymous, one non-synonymous, one splice site and two 3′UTR mutations) and two randomly selected intronic mutations within FADS1 and FADS2 genes (Table 1), out of 33 detected in this study were genotyped in 450 samples and used in association analyses. Genotyping was accomplished with TaqMan probes and a StepOnePlus™ Real-Time PCR system (Applied Biosystems, Foster City, CA, USA). All reagents required for the TaqMan assay including universal master mix, amplifying primers and probes were obtained from Applied Biosystems. One allelic probe was labeled with FAM dye and the other with the fluorescent VIC dye (Additional file 2: Table S2). PCR was run in the TaqMan universal master mix at a probe concentration of 200 nM. The reaction was performed in a 96-well format in a total reaction volume of 10 µL using 10 ng of gDNA. The real-time PCR cycling condition was as follows: reaction plates were heated for 30 secs at 60°C and then for 10 minutes at 95°C, followed by 35 cycles of 92°C for 15 secs and 60°C for 1 minute. The fluorescence intensity of each well in the TaqMan assay plate was subsequently read and the SNP call from the system (StepOne Software v2.2.2) was used to determine the genotype of each animal.

Table 1. Genotyped markers, frequencies of alleles and genotypes, and test for Hardy-Weinberg equilibrium

Additional file 2: Table S2. Primers and probes used in the TaqMan assay for SNP genotyping.

Format: DOCX Size: 16KB Download fileOpen Data

Statistical analyses

Statistical analyses were performed with SAS version 9.3 software (SAS Institute Inc., Cary, NC, USA). ALLELE procedure was used to compute allele frequencies, genotype frequencies and tests for HWE employing Fisher’s exact test. Associations between genotypes and the different FA profiles were evaluated using PROC MIXED with genotype as fixed effect (model 1).

Model 1: Yijk = μ + Hi + Gj + eijk

Where Yijk = trait measured; μ = overall mean; Hi = fixed effect of herd; Gj = fixed effect of genotype (j = 3, each SNP was analyzed separately with the same model) and eijk = residual error. A false discovery rate (FDR) was applied on the raw p-values of the global genotype effect using the MULTTEST procedure of SAS. For the FAs with significant results, least square means multiple comparisons were done with Tukey adjustments and considered to be significantly different at P ≤0.05. Spearman correlations between FAs were calculated over all data with the CORR procedure of SAS.

Allele substitution effects were estimated using PROC MIXED and the following model (Model 2).

Model 2: Yij = μ + Hi + ASE + eij

Where Yij = trait measured; μ = overall mean; Hi = fixed effect of herd; ASE = allele substitution effects (alleles at the same locus were coded 0, 1 or 2 [where 0 = CC or GG, 1 = CT or AG and 2 = TT or AA]) and treated as a linear covariant and eij = residual error. FDR was also applied on the raw p-values of the estimates.

Results

SNPs within the coding and regulatory regions of FADS1 and FADS2 genes in Canadian Holstein cows

A total of 33 SNPs were detected within the studied regions, 9 within FADS1 and 24 within FADS2, in Canadian Holstein cows (Additional file 3: Table S3).The SNPs have been submitted to the Single Nucleotide Polymorphism Database (dbSNP). Two SNPs are located within the coding regions, a synonymous mutation (FADS1-07, rs42187261) in exon 8 of FADS1 with no effect on amino acid 306 (306Tyr > Tyr) and a non-synonymous mutation (FADS2-14, rs211580559) within FADS2 exon 7 that causes a change in amino acid 294 from alanine to valine (294Ala > Val). One SNP (FADS2-05, rs211263660) occurred within a splice site. Seven SNPs were identified within the 3′UTR of FADS2 gene and Target Scan analysis indicated that FADS2-19 (rs210169303) is situated within the binding site for bta-miR-744. The presence of the mutated allele (T) may disrupt binding of this microRNA. Frequencies of genotyped alleles are shown in Table 1 and are all in agreement with Hardy-Weinberg equilibrium. The distribution of genotypes and their frequencies are shown in Table 1.

Additional file 3: Table S3. Identified SNPs within studied regions of FADS1 and FADS2 genes and functional annotation.

Format: DOCX Size: 22KB Download fileOpen Data

FADS SNP association with milk polyunsaturated fatty acids (PUFAs)

Results of significant SNP marker associations with milk PUFAs are shown in Table 2, while allele substitution results are presented in Table 3. Complete results, including significant and non-significant associations as well as substitution effects are shown in Additional file 4: Table S4 and Additional file 5: Table S5. Our results indicate that three markers (FADS1-01, FADS1-07 and FADS2-23) were significantly associated with three individual milk PUFAs (dihomogamma linolenic acid [C20:3n6], arachidonic acid [C20:4n6] and eicosapentaenoic acid [C20:5n3]) (FDR P < 0.05, Table 2). In addition, there was a significant association between FADS1-07 and tricosanoic acid (C23:0), a saturated FA. Out of these associations, significant relationship between FADS1-07, a synonymous mutation within FADS1 gene, and one omega-3 FA (C20:5n3, eicosapentanoic acid) showed that genotype TT was linked to the highest increase in eicosapentanoic acid. Allele substitution at this locus (allele T substituted for C) indicated that allele T increased milk C20:5n3 by 0.00420 ± 0.0008 (Table 3). A significant association was also recorded between one of the evaluated intronic mutations within FADS1 gene (FADS1-01[rs136261927]) and the omega-6 FA, C20:3n6. In this case, genotype CC was associated with the highest increase in milk C20:3n6 content. This result is supported by allele substitution analysis whereby substituting allele T for C at FADS1-01 associated a decrease of -0.009 ± 0.002 of milk C20:3n6 content to allele T (Table 3). FADS2-23, a 3′UTR variant, was significantly associated with the omega-6 FAs, C20:3n6 and C20:4n6, with genotype GG showing higher increases in the affected FAs (Table 2). Similarly, allele substitution effects shows that allele A at this locus decreased milk C20:3n6 and C20:4n6 significantly (Table 3).

Table 2. Results of significant1 marker associations with milk fatty acid isomers (g/100 g of total fat) of Canadian Holstein cows

Table 3. Results of significant1 allele substitution effects P-value

Additional file 4: Table S4. Complete marker-trait association results. Significant P-values (P ≤ 0.05) are highlighted yellow.

Format: RTF Size: 1.5MB Download fileOpen Data

Additional file 5: Table S5. Complete results on allele substitution effect. The reference genotype is either CC or GG. Significant p-values (P ≤ 0.05) are highlighted yellow.

Format: RTF Size: 763KB Download fileOpen Data

FADS SNP association with milk monounsaturated fatty acids (MUFAs)

Significant associations (raw p-values) recorded between two studied SNPs and two MUFA isomers as well as total MUFA (Tables 2 and 3) disappeared after FDR correction. Complete association results and allele substitution effects are shown in Additional file 4: Table S4 and Additional file 5: Table S5, respectively.

FADS SNP association with milk saturated fatty acids (SFAs)

Only one significant association was recorded between FADS1-07 and one saturated FA (C23:0, tracosanoic acid). Genotype TT was linked to the highest increases in tracosanoic acid. Allele substitution at this locus (allele T substituted for C) indicated that allele T increased milk C23:0 by 0.0021 ± 0.0008. Complete association results and allele substitution effects are shown in Tables 2 and 3, Additional file 4: Table S4 and Additional file 5: Table S5, respectively.

Correlations between fatty acids

Correlations between total SFAs and total MUFAs was positive but very low (0.099), and between total SFAs and total PUFAs was moderate but negative (-0.343) (Table 4). SFAs of 10 to 16 chain lengths showed very high but negative correlations with total MUFAs (-0.903) and low but also negative with total PUFAs (-0.268). Between total MUFAs and total PUFAs, results showed a low but positive relationship (0.234). Details of Pearson correlation coefficients between individual FAs are shown in Additional file 6: Table S6.

Table 4. Pearson correlation coefficients between total PUFAs, MUFAs, SFAs and different SFA classes

Additional file 6: Table S6. Pearson correlation coefficients between studied fatty acids.

Format: DOCX Size: 49KB Download fileOpen Data

Discussion

Our study is the first to report SNP characterization within bovine FADS1 and FADS2 genes in Canadian Holstein cows, and association with milk SFAs, MUFAs and PUFAs. FADS1 and FADS2 genes code for Δ-5 and Δ-6 desaturases, respectively and these enzymes are crucial in the endogenous synthesis of LC-PUFAs from the precursor essential FAs, LA and ALA, obtained from the diet. Numerous human studies have indicated that, changes in the levels of individual FAs are related to FADS enzyme activity [27-30] and recently, Merino et al [31] reported that genetic variants in FADS1 and FADS2 genes in humans can alter desaturase activity. Furthermore, it has been shown in humans that genetic variants of FADS1 and FADS2 influence blood lipid and breast milk essential FAs in pregnancy and lactation [34]. In this study, mutations were detected in the coding as well as non-coding (introns and 3′UTR) regions of studied genes. Mutations in these regions have been variously shown to associate significantly with economically important traits in farm animals, affect gene expression or function of resultant protein products [35,36].

The reported non-synonymous mutation within FADS2 gene has either valine or alanine at position 294 of the protein. These two amino acids are both nonpolar amino acids with aliphatic side chain groups. This change in amino acid does not seem to affect the structure and function of the resultant protein which may explain why no significant association was recorded between this SNP and any of studied FAs.

Significant associations were recorded between one PUFA (C20:5n3), and the synonymous FADS1 (FADS1-07, rs42187261) mutation. According to allele substitution results, allele T, the major allele with frequency of 57% in studied individuals, was superior in increasing the content of C20:5n3. It is worthwhile to point out that positive influence of this SNP on oleic acid and other health promoting FAs was significant before FDR correction, but not after FDR correction. Stearic acid (C18:0), a SFA and one of the main products of extensive biohydrogenation of USFAs in the ruminal environment is desaturated by the SCD1 enzyme [37] to oleic acid. The influence of this SNP (FADS1-07, rs42187261), together with favorable associations with SNPs detected within the SCD1 gene [38-40] may contribute to desaturate products of ruminal biohydrogenation or SFAs to much needed health promoting FAs. Furthermore, this mutation (rs42187261) also favored increase in the content of one individual SFA (C23:0). This study shows only a moderate correlation between C23:0 and C20:5n3 (r2 = 0.459). Furthermore, C23:0 is not among SFAs (C16:0, C14:0 and C12:0) which were associated with raised plasma low density lipoprotein cholesterol levels [41]. Since most SFAs in milk are synthesized de novo in the mammary gland through the activities of enzymes other than FADS1 and FADS2 while USFAs arise mostly from dietary sources, the effect of this mutation (rs42187261) may be more relevant for USFAs. Furthermore, very low and in most cases negative correlations were recorded between PUFAs positively influenced in this study and SFAs; implying that the effect of the rs42187261, rs109772589 and rs136261927 mutations can be useful in the positive management of these PUFAs.

Synonymous mutations have generally been regarded as silent mutations, that is, have no effect on protein structure and function. Recent evidence, however, indicates that silent mutations are able to cause changes in protein expression, conformation and function [36]. Numerous recent GWAS studies have revealed substantial contributions of synonymous SNPs to human disease risk and other complex traits [42-45]. Evidence suggest that synonymous SNPs can result in aberrant mRNA splicing [46], affect mRNA stability and thus protein expression and enzymatic activity [44], affect protein conformation and have functional and clinical consequences [47]. Consequently, Chamary and Hurst [48] estimated that 5-10% of human genes contain at least one region in which silent mutations could be harmful. Synonymous SNPs are generally not included in most association studies between SNPs and milk fat traits in ruminants. Therefore, important information that can assist in breeding programs could be missed. Our result on the synonymous FADS1 (FADS1-07, rs42187261) mutation support the concept of using both non-synonymous and synonymous mutations in animal breeding.

Our study also recorded a significant association between one intronic SNP (FADS1-01, rs136261927) and one MUFA (C20:3n6). Reports of significant associations between intronic SNPs in several lipogenic genes and milk fat traits as well as other production traits have emerged [35]. For example, intronic SNPs in the SCD1 (g.6926A > G) and FASN genes (g.8948C > T, ss491228481) were recently shown to be significantly associated with milk yield, fat yield and protein yield in Chinese Holstein cows [49]. Mutations in intronic regions are known to play important regulatory roles in mammalian gene expression, including roles in transcription regulation, polyadenylation, mRNA export, translational efficiency, and the rate of mRNA decay [50]. Our study further points to the importance of non-coding region SNPs in the modulation of gene expression and resultant effects on milk yield traits.

Conclusion

Our study has uncovered positive associations between three out of seven tested SNPs and health promoting FAs in the milk of Canadian Holstein cows. These SNPs may serve as potential genetic markers that can be used in breeding programs for increased milk beneficial FAs.

Competing interests

The authors have no competing interests to declare.

Authors’ contributions

EMIA and XZ designed the study and wrote the paper. KAA isolated DNA from milk somatic cells, prepared samples for Sanger sequencing and performed gas chromatographic analysis of milk fatty acids, FB performed SNP genotyping and contributed to manuscript writing. All authors read and approved the final manuscript.

Acknowledgements

We thank Celine Lesuffleur for help with DNA isolation from milk somatic cells, Valacta for facilitating collection of milk samples from participating herds, Dr. Arif Mustapha for help with FA profile analysis and Steve Méthot for statistical management of data. Financial support was provided by DairyGen (Dairy Cattle Genetics Research and Development Council of Canadian Dairy Network), NSERC (The Natural Sciences and Engineering Research Council of Canada) and Agriculture and Agri-Food Canada.

References

  1. Glaser C, Heinrich J, Koletzko B: Role of FADS1 and FADS2 polymorphisms in polyunsaturated fatty acid metabolism.

    Metabolism 2010, 59(7):993-999. PubMed Abstract | Publisher Full Text OpenURL

  2. Yang T, Fang S, Zhang H-X, Xu L-X, Zhang Z-Q, Yuan K-T, Xue C-L, Yu H-L, Zhang S, Li Y-F, et al.: N-3 PUFAs have antiproliferative and apoptotic effects on human colorectal cancer stem-like cells in vitro.

    J Nutr Biochem 2013, 24(5):744-753. PubMed Abstract | Publisher Full Text OpenURL

  3. Park B-K, Park S, Park J-B, Park MC, Min TS, Jin M: Omega-3 fatty acids suppress Th2-associated cytokine gene expressions and GATA transcription factors in mast cells.

    J Nutr Biochem 2013, 24(5):868-876. PubMed Abstract | Publisher Full Text OpenURL

  4. Prescott SL, Calder PC: N-3 polyunsaturated fatty acids and allergic disease.

    Curr Opin Clin Nutr Metab Care 2004, 7(2):123-129. PubMed Abstract | Publisher Full Text OpenURL

  5. Roqueta-Rivera M, Stroud CK, Haschek WM, Akare SJ, Segre M, Brush RS, Agbaga M-P, Anderson RE, Hess RA, Nakamura MT: Docosahexaenoic acid supplementation fully restores fertility and spermatogenesis in male delta-6 desaturase-null mice.

    J Lipid Res 2010, 51(2):360-367. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Jump DB: Dietary polyunsaturated fatty acids and regulation of gene transcription.

    Curr Opin Lipidol 2002, 13(2):155-164. PubMed Abstract | Publisher Full Text OpenURL

  7. Bhattacharya A, Banu J, Rahman M, Causey J, Fernandes G: Biological effects of conjugated linoleic acids in health and disease.

    J Nutr Biochem 2006, 17(12):789-810. PubMed Abstract | Publisher Full Text OpenURL

  8. Parodi PW: Conjugated linoleic acid and other anticarcinogenic agents of bovine milk fat.

    J Dairy Sci 1999, 82(6):1339-1349. PubMed Abstract | Publisher Full Text OpenURL

  9. Glaser C, Lattka E, Rzehak P, Steer C, Koletzko B: Genetic variation in polyunsaturated fatty acid metabolism and its potential relevance for human development and health.

    Matern Child Nutr 2011, 7:27-40. OpenURL

  10. Bauman DE, Lock AL: Conjugated linoleic acid: biosynthesis and nutritional significance. In Advanced Dairy Chemistry Volume 2 Lipids. Edited by Fox PF, McSweeney PLH. US: Springer; 2006:93-136. OpenURL

  11. Nakamura MT, Nara TY: Structure, function, and dietary regulation of δ6, δ5, and δ9 desaturases.

    Annu Rev Nutr 2004, 24(1):345-376. PubMed Abstract | Publisher Full Text OpenURL

  12. Sprecher H, Luthria DL, Mohammed BS, Baykousheva SP: Reevaluation of the pathways for the biosynthesis of polyunsaturated fatty acids.

    J Lipid Res 1995, 36(12):2471-2477. PubMed Abstract | Publisher Full Text OpenURL

  13. Chilliard Y, Glasser F, Ferlay A, Bernard L, Rouel J, Doreau M: Diet, rumen biohydrogenation and nutritional quality of cow and goat milk fat.

    Eur J Lipid Sci Technol 2007, 109(8):828-855. Publisher Full Text OpenURL

  14. Dewhurst RJ, Shingfield KJ, Lee MRF, Scollan ND: Increasing the concentrations of beneficial polyunsaturated fatty acids in milk produced by dairy cows in high-forage systems.

    Anim Feed Sci Technol 2006, 131(3-4):168-206. Publisher Full Text OpenURL

  15. Marchitelli C, Contarini G, De Matteis G, Crisà A, Pariset L, Scatà MC, Catillo G, Napolitano F, Moioli B: Milk fatty acid variability: effect of some candidate genes involved in lipid synthesis.

    J Dairy Res 2013, 80(02):165-173. PubMed Abstract | Publisher Full Text OpenURL

  16. Stoop WM, Schennink A, Visker MHPW, Mullaart E, van Arendonk JAM, Bovenhuis H: Genome-wide scan for bovine milk-fat composition: I: quantitative trait loci for short- and medium-chain fatty acids.

    J Dairy Sci 2009, 92(9):4664-4675. PubMed Abstract | Publisher Full Text OpenURL

  17. White SL, Benson GA, Washburn SP, Green JT Jr: Milk production and economic measures in confinement or pasture systems using seasonally calved Holstein and Jersey cows.

    J Dairy Sci 2002, 85(1):95-104. PubMed Abstract | Publisher Full Text OpenURL

  18. Bastin C, Soyeurt H, Gengler N: Genetic parameters of milk production traits and fatty acid contents in milk for Holstein cows in parity 1-3.

    J Anim Breed Genet 2013, 130:118-127. PubMed Abstract | Publisher Full Text OpenURL

  19. Soyeurt H, Dardenne P, Dehareng F, Bastin C, Gengler N: Genetic parameters of saturated and monounsaturated fatty acid content and the ratio of saturated to unsaturated fatty acids in bovine milk.

    J Dairy Sci 2008, 91(9):3611-3626. PubMed Abstract | Publisher Full Text OpenURL

  20. Gillingham LG, Harding SV, Rideout TC, Yurkova N, Cunnane SC, Eck PK, Jones PJH: Dietary oils and FADS1-FADS2 genetic variants modulate [13C]α-linolenic acid metabolism and plasma fatty acid composition.

    Am J Clin Nutr 2013, 97(1):195-207. PubMed Abstract | Publisher Full Text OpenURL

  21. Malerba G, Schaeffer L, Xumerle L, Klopp N, Trabetti E, Biscuola M, Cavallari U, Galavotti R, Martinelli N, Guarini P, et al.: SNPs of the FADS gene cluster are associated with polyunsaturated fatty acids in a cohort of patients with cardiovascular disease.

    Lipids 2008, 43(4):289-299. PubMed Abstract | Publisher Full Text OpenURL

  22. Rzehak P, Heinrich J, Klopp N, Schaeffer L, Hoff S, Wolfram G, Illig T, Linseisen J: Evidence for an association between genetic variants of the fatty acid desaturase 1 fatty acid desaturase 2 (FADS1 FADS2) gene cluster and the fatty acid composition of erythrocyte membranes.

    Br J Nutr 2009, 101(1):20-26. PubMed Abstract | Publisher Full Text OpenURL

  23. Schaeffer L, Gohlke H, Müller M, Heid IM, Palmer LJ, Kompauer I, Demmelmair H, Illig T, Koletzko B, Heinrich J: Common genetic variants of the FADS1 FADS2 gene cluster and their reconstructed haplotypes are associated with the fatty acid composition in phospholipids.

    Hum Mol Genet 2006, 15(11):1745-1756. PubMed Abstract | Publisher Full Text OpenURL

  24. Illig T, Gieger C, Zhai G, Romisch-Margl W, Wang-Sattler R, Prehn C, Altmaier E, Kastenmuller G, Kato BS, Mewes H-W, et al.: A genome-wide perspective of genetic variation in human metabolism.

    Nat Genet 2010, 42(2):137-141. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Kathiresan S, Melander O, Guiducci C, Surti A, Burtt NP, Rieder MJ, Cooper GM, Roos C, Voight BF, Havulinna AS, et al.: Six new loci associated with blood low-density lipoprotein cholesterol, high-density lipoprotein cholesterol or triglycerides in humans.

    Nat Genet 2008, 40(2):189-197. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  26. Tanaka T, Shen J, Abecasis GR, Kisialiou A, Ordovas JM, Guralnik JM, Singleton A, Bandinelli S, Cherubini A, Arnett D, et al.: Genome-wide association study of plasma polyunsaturated fatty acids in the InCHIANTI study.

    PLoS Genet 2009, 5(1):e1000338. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Truong H, DiBello JR, Ruiz-Narvaez E, Kraft P, Campos H, Baylin A: Does genetic variation in the Δ6-desaturase promoter modify the association between α-linolenic acid and the prevalence of metabolic syndrome?

    Am J Clin Nutr 2009, 89(3):920-925. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Martinelli N, Girelli D, Malerba G, Guarini P, Illig T, Trabetti E, Sandri M, Friso S, Pizzolo F, Schaeffer L, et al.: FADS genotypes and desaturase activity estimated by the ratio of arachidonic acid to linoleic acid are associated with inflammation and coronary artery disease.

    Am J Clin Nutr 2008, 88(4):941-949. PubMed Abstract | Publisher Full Text OpenURL

  29. Baylin A, Ruiz-Narvaez E, Kraft P, Campos H: α-Linolenic acid, Δ6-desaturase gene polymorphism, and the risk of nonfatal myocardial infarction.

    Am J Clin Nutr 2007, 85(2):554-560. PubMed Abstract | Publisher Full Text OpenURL

  30. Lu Y, Feskens EJ, Dolle ME, Imholz S, Verschuren WM, Muller M, Boer JM: Dietary n - 3 and n - 6 polyunsaturated fatty acid intake interacts with FADS1 genetic variation to affect total and HDL-cholesterol concentrations in the Doetinchem Cohort Study.

    Am J Clin Nutr 2010, 92(1):258-265. PubMed Abstract | Publisher Full Text OpenURL

  31. Merino DM, Johnston H, Clarke S, Roke K, Nielsen D, Badawi A, El-Sohemy A, Ma DWL, Mutch DM: Polymorphisms in FADS1 and FADS2 alter desaturase activity in young Caucasian and Asian adults.

    Mol Genet Metab 2011, 103(2):171-178. PubMed Abstract | Publisher Full Text OpenURL

  32. Montgomery GW, Sise JA: Extraction of DNA from sheep white blood cells.

    N Z J Agric Res 1990, 33(3):437-441. Publisher Full Text OpenURL

  33. O’Fallon JV, Busboom JR, Nelson ML, Gaskins CT: A direct method for fatty acid methyl ester synthesis: application to wet meat tissues, oils, and feedstuffs.

    J Anim Sci 2007, 85(6):1511-1521. PubMed Abstract | Publisher Full Text OpenURL

  34. Xie L, Innis SM: Genetic variants of the FADS1 FADS2 gene cluster are associated with altered (n-6) and (n-3) essential fatty acids in plasma and erythrocyte phospholipids in women during pregnancy and in breast milk during lactation.

    J Nutr 2008, 138(11):2222-2228. PubMed Abstract | Publisher Full Text OpenURL

  35. Ibeagha-Awemu E, Kgwatalala P, Zhao X: A critical analysis of production-associated DNA polymorphisms in the genes of cattle, goat, sheep, and pig.

    Mamm Genome 2008, 19(9):591-617. PubMed Abstract | Publisher Full Text OpenURL

  36. Sauna ZE, Kimchi-Sarfaty C: Understanding the contribution of synonymous mutations to human disease.

    Nat Rev Genet 2011, 12(10):683-691. PubMed Abstract | Publisher Full Text OpenURL

  37. Bauman DE, Lock AL: Milk fatty acid composition: challenges and opportunities related to human health. World BUIATRICS Congress; 2010:278-289. OpenURL

  38. Kgwatalala PM, Ibeagha-Awemu EM, Mustafa AF, Zhao X: Influence of stearoyl-coenzyme Adesaturase 1 genotype and stage of lactation on fatty acid composition of Canadian Jersey cows.

    J Dairy Sci 2009, 92(3):1220-1228. PubMed Abstract | Publisher Full Text OpenURL

  39. Schennink A, Heck JML, Bovenhuis H, Visker MHPW, van Valenberg HJF, van Arendonk JAM: Milk fatty acid unsaturation: genetic parameters and effects of stearoyl-CoA Desaturase (SCD1) and Acyl CoA: DiacylglycerolAcyltransferase 1 (DGAT1).

    J Dairy Sci 2008, 91(5):2135-2143. PubMed Abstract | Publisher Full Text OpenURL

  40. Taniguchi M, Utsugi T, Oyama K, Mannen H, Kobayashi M, Tanabe Y, Ogino A, Tsuji S: Genotype of stearoyl-CoA desaturase is associated with fatty acid composition in Japanese Black cattle.

    Mamm Genome 2004, 15(2):142-148. PubMed Abstract | Publisher Full Text OpenURL

  41. Ohlsson L: Dairy products and plasma cholesterol levels.

    Food Nutr Res 2010, 54:5124. OpenURL

  42. Bartoszewski RA, Jablonsky M, Bartoszewska S, Stevenson L, Dai Q, Kappes J, Collawn JF, Bebok Z: A synonymous single nucleotide polymorphism in ΔF508 CFTR alters the secondary structure of the mRNA and the expression of the mutant protein.

    J Biol Chem 2010, 285(37):28741-28748. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  43. Macaya D, Katsanis SH, Hefferon TW, Audlin S, Mendelsohn NJ, Roggenbuck J, Cutting GR: A synonymous mutation in TCOF1 causes Treacher Collins syndrome due to mis-splicing of a constitutive exon.

    Am J Med Genet A 2009, 149A(8):1624-1627. PubMed Abstract | Publisher Full Text OpenURL

  44. Nackley AG, Shabalina SA, Tchivileva IE, Satterfield K, Korchynskyi O, Makarov SS, Maixner W, Diatchenko L: Human catechol-O-methyltransferase haplotypes modulate protein expression by altering mRNA secondary structure.

    Science 2006, 314(5807):1930-1933. PubMed Abstract | Publisher Full Text OpenURL

  45. Ramser J, Ahearn ME, Lenski C, Yariz KO, Hellebrand H, von Rhein M, Clark RD, Schmutzler RK, Lichtner P, Hoffman EP, et al.: Rare missense and synonymous variants in UBE1 are associated with X-linked infantile spinal muscular atrophy.

    Am J Hum Genet 2008, 82(1):188-193. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  46. Cartegni L, Chew SL, Krainer AR: Listening to silence and understanding nonsense: exonic mutations that affect splicing.

    Nat Rev Genet 2002, 3(4):285-298. PubMed Abstract | Publisher Full Text OpenURL

  47. Kimchi-Sarfaty C, Oh JM, Kim I-W, Sauna ZE, Calcagno AM, Ambudkar SV, Gottesman MM: A “silent” polymorphism in the MDR1 gene changes substrate specificity.

    Science 2007, 315(5811):525-528. PubMed Abstract | Publisher Full Text OpenURL

  48. Chamary JV, Hurst LD: The price of silent mutations.

    Sci Am 2009, 300(6):46-53. Publisher Full Text OpenURL

  49. Alim MA, Wang P, Wu XP, Li C, Cui XG, Zhang SL, Zhang Q, Zhang Y, Sun DX: Effect of FASN gene on milk yield and milk composition in the Chinese Holstein dairy population.

    Anim Genet 2014, 45(1):111-113. Publisher Full Text OpenURL

  50. Nott A, Meislin SH, Moore MJ: A quantitative analysis of intron effects on mammalian gene expression.

    RNA 2003, 9(5):607-617. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL