Associations between the Bovine Myostatin Gene and Milk Fatty Acid Composition in New Zealand Holstein-Friesian X Jersey-Cross Cows.
1. IntroductionImproving the efficiency of cattle production systems can be achieved by selecting for fast growing animals with increased muscling that have desirable maternal reproductive traits, good milk production, and good mothering ability. To achieve this, it is important to have an understanding of the genes that underpin muscularity and adiposity.
The gene for the protein myostatin (MSTN; gene MSTN), which is also called the growth and differentiation factor 8 gene (GDF8) has pleiotropic effects. Its expression has been associated with decreased adipogenesis and increased skeletal muscle mass as a result of decreased secretion of leptin [1-3].
Sequence variation in MSTN has been associated with increases in growth and muscling traits in several species. For example, in cattle it has been connected with having increased numbers of muscle fibers (otherwise known as "double-muscling") in a number of breeds [4]. Similarly, sequence variation in the first intron has been found to influence growth and carcass traits such as the yield of leg, loin, and total lean meat in NZ Romney sheep [5].
While there is well-documented evidence describing how variation in MSTN is associated with growth and muscle traits in beef cattle breeds, there is little information about the effects of MSTN nucleotide sequence variation on milk yield and milk fatty acid (FA) composition. There are suggestions that MSTN could affect lactation by affecting the production of fatty acid in milk, and primarily through MSTN deficiency being associated with decreased adipogenesis [2,3,6].
The fat of milk from dairy cattle is approximately 70% saturated fatty acid (SFA), 25% monounsaturated fatty acid (MUFA), and 5% polyunsaturated fatty acid (PUFA) [7]. From a human health perspective, an increase in unsaturated fatty acid (UFA) content and a decrease in SFA content could be considered favorable; and equally from a physical perspective (e.g., increased spread-ability of butter), increased relative levels of UFA would be desirable. Previously, it has been identified that the concentrations of different fatty acids in milk are affected most by four parameters: the diet of the cow [8,9], genetic variation between cows within a breed [10,11], breed differences [12,13], and the number of days in milk [14]. In their report, Strucken et al. [14] established that cows in early lactation stage are characterized by negative energy balance as the dietary energy intake is unable to meet the demands of high milk production in approximately the first 60 days of lactation. To offset this balance, an alternative energy source is needed. This leads to the mobilization of body energy stores to balance the deficit between feed intake, and energy expenditure on maintenance and milk production [15]. However, as a consequence of the cow's body fat being used up for this purpose, other biological pathways are affected, resulting in a change in milk composition. The Holstein-Friesian and the Jersey breeds have one of the most notable differences in terms of the composition of milk fatty acids. Milk from Jersey cows tend to have higher concentrations of some short- and medium-chain length SFA, but lower concentrations of some UFA [16]. These differences could be capitalized upon in order to obtain the preferred fatty acid profile through cross-breeding of these breeds.
In the above context, the objective of this research was to explore the effects of MSTN variation on key milk production traits and composition of milk fatty acids in Holstein-Friesian * Jersey (HF * J)-cross (alternatively known as Kiwicross[TM]) dairy cows in Hamilton, New Zealand (NZ). This cross is now the preferred cow for dairy production in NZ and it constitutes more than half of all the dairy cows milked annually. Our working hypothesis is that variation in the gene will affect some milk traits.
2. Materials and Methods
2.1. The Dairy Cattle Investigated and the Collection of Milk Samples from Them
This research was approved by the Lincoln University Animal Ethics Committee (AEC Approval Number 521). This is a mandated and registered committee that was established, and is regularly audited, under the provisions of the New Zealand Government's Animal Welfare Act 1999 (http://www.legislation.govt.nz/act/public/1999/0142/latest/DLM49664.html). A total of 410 HF * J-cross dairy cows (alternatively known as Kiwicross[TM] cows, in Hamilton, New Zealand), of individually variable and unidentified breed ratio, and of 3 to 9 years of age were investigated. The cows were managed in two herds on the Lincoln University Dairy Farm (Canterbury, New Zealand), and all of them were grazed outdoors at all times on pasture (a blend of white clover and perennial ryegrass). All of the cows calved over the months of August-September (spring in the Southern Hemisphere) and they were then milked twice a day for up to 10 months.
Milk samples for gross milk trait analysis were collected once a month from September to February and the daily milk yield in liters was recorded using Tru-test milk meters (Tru-test Ltd., Auckland, New Zealand) and Fourier-transform infra-red spectroscopy (MilkoScan FT 120 Foss, Hillerod, Denmark) was used to analyze milk samples for fat percentage (%) and protein percentage (%). Average daily milk yield, and average protein and fat percentages were calculated over the 6 months of milk collection. The milk samples for fatty acid (FA) analysis were collected at 148 [+ or -] 19 days in milk from each cow in a single afternoon milking in mid-January (the middle of summer in the Southern Hemisphere). These were frozen at -20 [degrees]C, and then freeze-dried, prior to being individually ground to a fine powder for component analysis.
2.2. Gas Chromatography of the Fatty Acids in the Milk Samples
The FAs were methylated and extracted in n-heptane, before being analyzed by gas chromatography (GC) as FA methyl esters (FAMEs). The methylation reactions were performed in 10-mL Kimax tubes. Individual freeze-dried and powdered milk samples (0.17 g) were dissolved in 900 [micro]L of n-heptane (100%, AR grade), before 100 [micro]L of internal standard (5 mg/mL of C21:0 methyl ester in n-heptane) and 4.0 mL of 0.5 M NaOH (in 100% anhydrous methanol) were added.
The tubes were vortexed then incubated in a block heater (Ratek Instruments, Australia) at 50 [degrees]C for 15 min. After cooling to room temperature, another 2.0 mL of n-heptane and 2.0 mL of deionized water were added to each tube. After vortexing, the tubes were centrifuged (Megafuge 1.0R, Heraeus, Germany) for 5 min at 1500* g. The top layer of n-heptane was transferred into a second Kimax tube and 2.0 mL of n-heptane was added to each of the original tubes. The extraction was repeated and the n-heptane aspirates were then pooled. Anhydrous sodium sulphate (10 mg) was added to the n-heptane extracts, to remove any residual water.
The GC analysis was carried out using a Shimadzu GC-2010 Gas Chromatograph (Shimadzu Corporation, Kyoto, Japan) equipped with a flame ionization detector and an AOC-20i auto sampler. The output was analyzed with GC Solution Software (Shimadzu). For analysis, 1 [MICRO]L of the n-heptane sample extract was injected into a 100 m GC column (250 [micro]m * 0.25 [micro]m capillary column, CP-Select, Varian) with a 1:60 split ratio. The separation was undertaken with a helium carrier gas and was run for 92 min. The temperature of both the injector and detector were set at 250 [degrees]C and the thermal profile of the column consisted of 45 [degrees]C for 4 min, followed by 27 min at 175 [degrees]C (ramped at 13 [degrees]C/min), 35 min at 215 [degrees]C (ramped at 4 [degrees]C/min), and a final "bake-off" at 250 [degrees]C for 5 min (ramped at 25 [degrees]C/min). The individual FAMEs were identified by the peak retention time compared to commercially obtained external standards (ME61, ME93, BR3, BR2, ME100, GLC411, and GLC463; Laroden AB, Sweden). Quantification of the individual FAMEs was based on peak area assessment and comparison with the internal and external standards. The threshold for peak area determination on the chromatogram was a 500-unit count, with peaks that were under 500-unit count, being ignored. The calculated minimum component of an individual FAME was therefore 0.01 g per 100 g of total FA.
After their individual measurement, the FAs were arranged into various groups and indices. These groups were, short-chain length FAs (SCFA) = C4:0 + C6:0 + C8:0; medium-chain length FAs (MCFA) = C10:0 + C12:0 + C14:0; long-chain length FAs (LCFA) = C15:0 + C16:0 + C17:0 + C18:0 + C19:0 + C20:0 + C22:0 + C24:0; omega 3 FAs = C18:3 cis-9,12,15 + C20:5 cis-5,8,11,14,17 + C22:5 cis-7, 10,13,16,19; omega 6 FAs = C18:2 cis-9,12 + C18:3 cis-6, 9,12 + C20:3 cis-8,11,14 + C20:4 cis-5, 8, 11,14; monounsaturated FAs (MUFA) = C10:1 + C12:1 + C14:1 cis-9 + C15:1 + C16:1 cis-9 + C17:1 + C18:1 trans-11 + C18:1 cis-9 + C18:1 cis-(10 to 15) + C20:1 cis-5 + C20:1 cis-9 + C20:1 cis-11 + C22:1 trans-13; polyunsaturated FAs (PUFA) = C18:2 trans-9,12 + C18:2 cis-9,trans-13 + C18:2 cis-9,trans-12 + C18:2 trans-9,cis-12 + C18:2 cis-9,12 + C18:3 cis-6, 9,12 + C18:3 cis-9,12,15 + CLA + C20:3 cis-8,11, 14 + C20:4 cis-5, 8, 11, 14 + C20:5 cis-5, 8,11,14,17 + C22:5 cis-7,10,13, 16, 19; and total branched FA = C13:0 iso + C13:0 anteiso + C15:0 iso + C15:0 anteiso + C17:0 iso.
Unsaturated FA indices were also calculated as follows: C12:1 index (C12:1 divided by the sum of C12:0 andC12:1); C16:1 index (C16:1 cis-9 divided by the sum of C16:0 and C16:1 cis-9); andC18:1 index (C18:1 cis-9 divided by the sum of C18:0 and C18:1 cis-9). The method is as described by Li et al. [17], with the un-adjusted mean levels in the 430 cows being calculated and used subsequently in the statistical analyses.
2.3. Blood Sample Collection
Samples of blood were collected from each of the cows studied onto FTA[TM] cards (Whatman[TM], Middlesex, UK) by piercing the ear of the animal. This is allowed under a Code of Welfare issued by the NZ Minister of Agriculture, under Section 75 and 76 of the NZ Animal Welfare Act 1999. The samples were allowed to dry in the air and the purification of the DNA from 1.2-mm punches taken from the FTA[TM] cards was carried out using a two-step procedure described by Zhou et al. [18].
2.4. Polymerase Chain Reaction (PCR) Amplification of a Region of the Cattle Myostain Gene
The intron 1 region of MSTN was amplified using forward and reverse primers (5'-catggtactattgttgagag-3' and 5'-aaggcaaatctattccagg-3' respectively) adapted from the work of Haruna et al. [19]. The 15-[micro]L reactions contained the purified DNA on a 1.2-mm diameter disc of the FTA[TM] card, and a content of 0.25 [micro]M for each primer, 150 [micro]M for each dNTP (Eppendorf, Hamburg, Germany), 3.0 mM [Mg.sup.2+], 0.5 U of Taq DNA polymerase (Qiagen, Hilden, Germany), and 1x the buffer supplied with the DNA polymerase enzyme.
The PCR amplifications were carried out in Bio-Rad S1000 thermal cyclers (Bio-Rad, Hercules, CA, USA) and the thermal cycling parameters included an initial denaturation at 94 [degrees]C for 2 min, and then 35 repeated cycles of denaturation at 94 [degrees]C for 30 s, primer annealing at 58 [degrees]C for 30 s, and primer extension at 72 [degrees]C for 30 s. Following this process, a final extension step at 72 [degrees]C for 5 min was used.
2.5. Single Strand Conformation Polymorphism (SSCP) Analysis
An SSCP technique was employed to detect nucleotide sequence variation in the amplicons obtained from the PCR reactions. A 0.7-[micro]L aliquot of the completed reactions was added to 7 [micro]L of a solution containing 10 mM ethylenediaminetetraacetic acid (EDTA), 0.025% bromophenol blue, 0.025% xylene-cyanol, and 98% formamide. The samples were then placed on a hot plate already set at 95 [degrees]C for 5 min for denaturation, followed by immediate cooling on wet ice. They were then loaded onto 14% acrylamide: bisacrylamide (37.5:1) (Bio-Rad) gels that were 16 cm X 18 cm in size. Electrophoresis was undertaken using Protean II xi cells (Bio-Rad) for 19 h at 390 volts and 7 [degrees]C room temperature in 0.5 X Tris/Borate/EDTA running buffer.
To detect the SSCP banding patterns, the SSCP gels were stained using the silver-staining method of Byun et al. [20].
2.6. Nucleotide Sequencing and Sequence Analysis
Based on the PCR-SSCP patterns observed nucleotide sequence variation could be identified. Cows that appeared to be homozygous with unique banding patterns were subjected to direct sequencing. For cows that appeared to have heterozygous variant patterns, unique bands were excised from the wet gel and incubated in water at 69 [degrees]C for 1 h. A 1-[micro]L aliquot of the water product was pipetted into 14-[micro]L of PCR pre-mixture (as used for the original PCR reactions), re-amplified using the same thermal cycle profile, and subsequently sequenced. This approach has been described in more detail by Gong et al. [21]. The sequences were aligned, and other analyses were undertaken using version 5.2.10 of DNAMAN (Lynnon BioSoft, Vaudreuil, QC, Canada).
2.7. Statistical Analyses
All statistical analyses were performed using IBM SPSS version 22 (IBM, Armonk, NY, USA), and an alpha level of p < 0.05 was set as a threshold.
For genotypes with a frequency greater than 5% (thus having adequate sample size per group), the effect of variation in a cow's MSTN genotype on gross milk production traits, and the component levels of individual and grouped FAs was tested using general linear mixed-effects models (GLMMs) and multiple pair-wise comparisons (least significant difference tests) with Bonferroni corrections. The age of the cow expressed in an integer value of years (i.e., as a categorical variable in a range from 3 to 9 years of age), the number of days in milk for each cow (DIM; expressed as an integer value, but entered into the model as a continuous trait) and herd (to correct for herd-specific effects) were fitted to the models as fixed explanatory factors.
The model was [Y.sub.ijkl] = [micro] + [G.sub.i] + [A.sub.j] + [D.sub.k] + [H.sub.l] + [e.sub.ijkl] for genotype: where [Y.sub.ijkl] = the observed trait value in the ijklth cow; [micro] = the mean trait value for a given trait; [G.sub.i] = the fixed effect of ith MSTN genotype; [A.sub.j] = effect of age (j = 3-9 years); [D.sub.k] = effect of the number of days the cow has produced milk (DIM: k = 94-186 days); [H.sub.l] = the fixed effect of lth farm (l = 1 or 2); and [e.sub.ijkl] = random error.
The effect of sire of cow could not be included in the GLMMs, because some semen straws (sire genetics) used in NZ dairy cattle artificial insemination-based breeding approaches contain mixed-sire semen purchased from commercial semen producers. In these cases, it is impossible to ascertain individual sire identity. However, since the straws were mixed-semen straws and because different sires are used for different inseminations, in different years, it is unlikely that sire was a strongly confounding effect. Cow age and herd might also be confounded with sire, but this cannot be confirmed.
3. Results
3.1. Milk Production
Over the time the cows were studied, they produced an average of 22.53 [+ or -] 2.18 L of milk per day, with an average milk fat content of 4.94 [+ or -] 0.17% and average milk protein content of 4.03 [+ or -] 0.10%.
3.2. Identification of Nucleotide Sequence Variation in Bovine MSTN
A 406 bp fragment of the intron 1 region of bovine MSTN was amplified and analyzed using the polymerase chain reaction coupled with single strand conformation polymorphism (PCR-SSCP) analyses.
The PCR-SSCP analyses coupled with DNA sequencing revealed five banding patterns (A-E) in the region of intron 1 investigated (Figure 1). A total of seven single-nucleotide substitutions (c.373+751G/T, c.373+803T/G, c.373+877A/G, c.373+895G/C, c.374-909C/T, c.374-842G/C, c.374-812A/G) were identified, all of which have been previously reported [19].
3.3. MSTN Genotype Models
The genotypes AA (n = 149), AB (n = 87), AC (n = 48) and AD (n = 63) occurred at a frequency over 5% and were analyzed in this model. The other genotypes AE (n = 16), BB (n = 6), BC (n = 5), BD (n = 15), CC (n = 5), CD (n = 7), CE (n = 1) and DD (n = 3) had frequencies of less than 5% each, and were not included in the model.
In the genotype models (Table 1), genotype AD was associated with a lower C12:0 and C12:1 levels than genotype AA.
4. Discussion
This is the first study demonstrating association of MSTN sequence variants with the component levels of two milk FAs in NZ HF * J cross cattle.
All seven of the nucleotide substitutions identified in this study have been previously reported in a study of NZ cattle breeds, which included; Hereford, Angus, Shorthorn, Charolais, Red Poll, South Devon, Simmental, Murray Grey, HF * J cross cattle, and some composite breeds [19]. The nucleotide substitutions c.373+751G/T, c.373+803T/G, c.373+877A/G, c.373+895G/C, and c.374-909C/T were identified in all the ten aforementioned breeds, while c.374-842G/C was found in all but four breeds (Red Poll, Shorthorn, Simmental, and Composites breeds) and the c.374-812A/G was only found in Shorthorn and HF * J cross breeds.
In a previous study, Smith et al. [22] identified a mutant "mh" allele of bovine MSTN with an 11-bp deletion in the exon 3 region. This was described in British South Devon cattle, and it was the same allele first identified in the Belgian Blue breed, causing the "double-muscled" phenotype. Also, in a later investigation of MSTN in 146 British South Devon cattle, Wiener et al. [23] showed that this same mutant "mh" allele, reduced the levels of total SFA and total MUFA in the muscle (p < 0.05). Their report also revealed an increase in the ratio of PUFA: SFA in total lipid to be greater in mh/mh individuals than in the other two genotype (mh/+ and +/+) classes (p < 0.001). This suggested that the "mh" allele is associated with reduced fat levels, particularly with the levels of SFA and MUFA. The increase in the ratio of PUFA to SFA would be expected, especially when the concentration of muscle fat decreases.
In a similar report on associations between the "mh" allele carrying the 11-bp deletion with intramuscular fatty acid composition in MSTN of Belgian Blue young bulls, Raes et al. [24] revealed that animals with the +/+ (normal) genotype showed a higher relative amount of the SFA; C14:0 and C16:0 and a higher relative amount of all MUFA C16:1; C17:1 and C18:1, but the relative proportion of PUFA in total fatty acids increased with increasing mutant "mh" alleles. In the current study, the genotype model revealed that cows carrying the AA genotype showed an increase in the amount of two medium-chain length SFAs (C12:0 and C12:1) relative to AD cows. While it may be difficult to conclude that the findings in the current study are similar to the work of Weiner et al. [23] and Raes et al. [24], especially because it investigated the intronic region of MSTN in non-doubled muscled cattle breed, it is important to note that variations in the intronic regions are equally capable of influencing gene expression and/or altering the functionality of a gene [25]. In a previous investigation, He et al. [25] transformed C2C12 cell-lines with a transgene construct that contained bovine MSTN promoter (pMD-MSTNPro) and a second construct that contained the first intron of bovine MSTN (pMD-Intron1). They observed an increase in the mean fluorescence intensity of green fluorescent protein (GFP) gene and the percentage of fluorescence positive cells, and concluded from this that the presence of intron 1 of bovine MSTN increased the expression of GFP in the transformed cells.
The findings in this study are in part in agreement with the findings of Buske et al. [26]. That study revealed that one copy of the so-called "wild-type + allele" of MSTN was associated with higher milk, protein, and fat yields in dual purpose Belgian Blue (DP-BB) cows, whereas a single copy of the mutant "mh" allele (associated with double muscling) was associated with a decrease in the SFA content of milk. Even though the current study did not involve the MSTN "mh" allele, it has perhaps suggested that antagonistic effects may exist between milk and meat production traits in the context of MSTN.
Since the MSTN variations in the current study occurred in a non-expressed region of the gene, it is perhaps less obvious how they might affect the structure or function of the MSTN protein. They might however influence the rate of transcription and/or translation of the gene, as was described by Liu et al. [27] in an investigation that involved both transgenic mice and in vitro studies. In that study, intron sequences in the transgene that encoded rat growth hormone were observed to stimulate transcription by promoting assembly of an ordered nucleosome array in the vicinity of the promoter. Additionally, while the observed intron variation does not yield any amino acid change, several studies have shown that synonymous nucleotide substitutions can affect the phenotypic characteristics of the protein product by altering mRNA structure, protein stability, the electrical charge of the resulting polypeptide, and codon usage during mRNA translation [28,29].
Previous reports on the differences in the composition of milk FAs between Holstein-Friesian and the Jersey breeds suggest that milk from Jersey cows tends to have higher concentrations of some short-and medium-chain length SFA, but lower concentrations of some UFA [16]. It is therefore interesting to note that the current investigation of NZ HF * J cross cows, revealed that variant D was associated with a decrease in the amount of two medium-chain length SFAs in milk. This discovery could be of benefit in terms of its potential applicability in cross-breeding and gene-marker development, particularly in selection for decreased SFA in milk.
5. Conclusions
The findings suggest that variation in MSTN affects two milk FA traits and this may be of value in breeding dairy cattle. Notably, variant D of MSTN is associated with a decrease in two medium-chain length SFA levels. Cows with the AD genotype might therefore produce a "preferred" FA profile in milk. However, because there were insufficient cattle with the homozygous genotypes BB, CC, and DD, or the heterozygous genotypes (BC, BD, CE, and CD) in the samples investigated, further investigation involving much larger sample sizes across different farms and breeds of cattle is needed to validate this claim.
Author Contributions: Conceptualization, J.G.H.H. and H.Z.; methodology, I.L.H., U.J.E., and H.A.-N.; software, Y.L., and H.A.-N.; validation, J.G.H.H., and H.Z.; formal analysis, Y.L., H.A.-N.; investigation, I.L.H. and U.J.E.; resources, J.G.H.H. and H.Z.; data curation, H.Z.; writing--original draft preparation, I.L.H.; writing--review and editing, I.L.H., J.G.H.H., and H.Z.; visualization, H.Z.; supervision, J.G.H.H.; project administration, J.G.H.H. All authors have read and agreed to the published version of the manuscript.
Funding: This research received no external funding.
Acknowledgments: The authors would like to thank Freeman Fang for technical support.
Conflicts of Interest: The authors declare no conflict of interest.
References
[1.] McPherron, A.C.; Lawler, A.M.; Lee, S.J. Regulation of skeletal muscle mass in mice by a new TGF-b superfamily member. Nature 1997, 387, 83-90. [CrossRef] [PubMed]
[2.] Mendias, C.L.; Bakhurin, K.I.; Faulkner, J.A. Tendons of myostatin-deficient mice are small, brittle, and hypocellular. Proc. Natl. Acad. Sci. USA 2008, 105, 388-393. [CrossRef] [PubMed]
[3.] Whittemore, L.A.; Song, K.; Li, X.; Aghajanian, J.; Davies, M.; Girgenrath, S.; Hill, J.J.; Jalenak, M.; Kelley, P.; Knight, A.; et al. Inhibition of myostatin in adult mice increases skeletal muscle mass and strength. Biochem. Biophys. Res. Commun. 2003, 300, 965-971. [CrossRef]
[4.] Kambadur, R.; Sharma, M.; Smith, T.P.; Bass, J.J. Mutations in myostatin (GDF8) in double-muscled Belgian Blue and Piedmontese cattle. Genome Res. 1997, 7, 910-916. [CrossRef]
[5.] Hickford, J.G.H.; Forrest, R.H.; Zhou, H.; Fang, Q.; Han, J.; Frampton, C.M.; Horrell, A.L. Polymorphisms in the ovine myostatin gene (MSTN) and their association with growth and carcass traits in New Zealand Romney sheep. Anim. Genet. 2010, 41, 64-72. [CrossRef]
[6.] McPherron, A.C.; Lee, S.J. Suppression of body fat accumulation in myostatin-deficient mice. J. Clin. Investig. 2002, 109, 595-601. [CrossRef]
[7.] Grummer, R.R. Effects of feed on the composition of milk fat. J. Dairy Sci. 1991, 74, 3244-3257. [CrossRef]
[8.] Palmquist, D.L.; Beaulieu, A.D.; Barbano, D.M. Feed and animal factors influencing milk fat composition. J. Dairy Sci. 1993, 76, 1753-1771. [CrossRef]
[9.] Dewhurst, R.J.; Shingfield, K.J.; Lee, M.R.F.; Scollan, N.D. Increasing the concentrations of beneficial polyunsaturated fatty acids in milk produced by dairy cows in high-forage systems. Anim. Feed Sci. Technol. 2006, 131, 168-206. [CrossRef]
[10.] Soyeurt, H.; Dardenne, P.F.; Dehareng, G.; Lognay, G.; Veselko, D.; Marlier, M.; Bertozzi, C.; Mayeres, P.; Gengler, N. Estimating fatty acid content in cow milk using mid-infrared spectrometry. J. Dairy Sci. 2006, 89, 3690-3695. [CrossRef]
[11.] Stoop, W.M.; van Arendonk, J.A.M.; Heck, J.M.L.; van Valenberg, H.J.F.; Bovenhuis, H. Genetic parameters for major milk fatty acids and milk production traits of Dutch Holstein Friesians. J. Dairy Sci. 2008, 91, 385-394. [CrossRef] [PubMed]
[12.] Palladino, R.A.; Buckley, F.; Prendiville, R.; Murphy, J.J.; Callan, J.; Kenny, D.A. A comparison between Holstein-Friesian and Jersey dairy cows and their f1 hybrid on milk fatty acid composition under grazing conditions. J. Dairy Sci. 2010, 93, 2176-2184. [CrossRef] [PubMed]
[13.] Maurice-Van Eijndhoven, M.H.T.; Hiemstra, S.J.; Calus, M.P.L. Short communication: Milk fat composition of 4 cattle breeds in the Netherlands. J. Dairy Sci. 2011, 94, 1021-1025. [CrossRef] [PubMed]
[14.] Strucken, E.M.; Laurenson, Y.C.S.M.; Brockmann, G.A. Go with the flow-biology and genetics of the lactation cycle. Front. Genet. 2015, 6, 118. [CrossRef] [PubMed]
[15.] Bauman, D.E.; Bruce Currie, W. Partitioning of nutrients during pregnancy and lactation: A review of mechanisms involving homeostasis and homeorhesis. J. Dairy Sci. 1980, 63, 1514-1529. [CrossRef]
[16.] Arnould, V.M.R.; Soyeurt, H. Genetic variability of milk fatty acids. J. Appl. Genet. 2009, 209, 29-39. [CrossRef]
[17.] Li, Y.; Zhou, H.; Cheng, L.; Edwards, G.R.; Hickford, G.H.J. Effect of DGAT1 variant (K232A) on milk traits and milk fat composition in outdoor pasture-grazed dairy cattle. NZ. J. Agric. Res. 2019. [CrossRef]
[18.] Zhou, H.; Hickford, G.H.J.; Fang, Q. A two-step procedure for extracting genomic DNA from dried blood spots on filter paper for polymerase chain reaction amplification. Anal. Biochem. 2006, 354, 159-161. [CrossRef]
[19.] Haruna, I.L.; Ekegbu, U.J.; Ullah, F.; Amirpour-Najafabadi, H.; Zhou, H.; Hickford, J.G.H. Genetic variation and haplotypic diversity in the myostatin gene of New Zealand cattle breeds. Gene 2020, 740. [CrossRef]
[20.] Byun, S.O.; Fang, Q.; Zhou, H.; Hickford, J.G.H. An effective method for silver-staining DNA in large numbers of polyacrylamide gels. Anal. Biochem. 2009, 385, 174-175. [CrossRef]
[21.] Gong, H.; Zhou, H.; Hickford, J.G.H. Diversity of the glycine/tyrosine-rich keratin-associated protein 6 gene (KAP6) family in sheep. Mol. Biol. Rep. 2011, 38, 31-35. [CrossRef] [PubMed]
[22.] Smith, J.A.; Lewis, A.M.; Wiener, P.; Williams, J.L. Genetic variation in the bovine myostatin gene in UK beef cattle: Allele frequencies and haplotype analysis in the South Devon. Anim. Genet. 2000, 31, 306-309. [CrossRef] [PubMed]
[23.] Wiener, P.; Woolliams, J.A.; Frank-Lawale, A.; Ryan, M.; Richardson, R.I.; Nute, G.R.; Wood, J.D.; Homer, D.; Williams, J.L. The effects of a mutation in the myostatin gene on meat and carcass quality. Meat Sci. 2009, 83, 127-134. [CrossRef] [PubMed]
[24.] Raes, K.; de Smet, S.; Demeyer, D. Effect of double muscling in Belgian Blue young bulls on the intramuscular fatty acid composition with emphasis on conjugated linoleic acid and polyunsaturated fatty acids. Anim. Sci. 2001, 73, 253-260. [CrossRef]
[25.] He, Y.L.; Wu, Y.H.; Quan, F.S.; Liu, Y.G.; Zhang, Y. Molecular analysis of the first intron in the bovine myostatin gene. Mol. Biol. Rep. 2010, 38, 4643-4649. [CrossRef]
[26.] Buske, B.; Gengler, N.; Soyeurt, H. Short-communication: Influence of the muscle hypertrophy mutation of the myostatin gene on milk production traits and milk fatty acid composition in dual-purpose Belgian Blue dairy cattle. J. Dairy Sci. 2011, 94, 3687-3692. [CrossRef]
[27.] Liu, K.; Sandgren, E.P.; Palmiter, R.D.; Stein, A. Rat growth hormone gene introns stimulate nucleosome alignment in vitro and in transgenic mice. Proc. Natl. Acad. Sci. USA 1995, 92, 7724-7728. [CrossRef]
[28.] Sauna, Z.E.; Kimchi-Sarfaty, C. Understanding the contribution of synonymous mutations to human disease. Nat. Rev. Genet. 2011, 12, 683-691. [CrossRef]
[29.] Hunt, R.C.; Simhadri, V.L.; Iandoli, M.; Sauna, Z.E.; Kimchi-Sarfaty, C. Exposing synonymous mutations. Trends Genet. 2014, 30, 308-321. [CrossRef]
Ishaku L. Haruna (iD), Yunhai Li, Ugonna J. Ekegbu, Hamed Amirpour-Najafabadi (iD), Huitong Zhou and Jon G. H. Hickford (*) (iD)
Gene-Marker Laboratory, Department of Agricultural Sciences, Faculty of Agriculture and Life Sciences, Lincoln University, Lincoln 7647, New Zealand; Ishaku.Haruna@lincolnuni.ac.nz (I.L.H.); Yunhai.Li@Lincolnuni.ac.nz (Y.L.); Jane.Ekegbu@lincolnuni.ac.nz (U.J.E.); Hamed.Amirpour@lincolnuni.ac.nz (H.A.-N.); zhouh@lincoln.ac.nz (H.Z.)
(*) Correspondence: jon.hickford@lincoln.ac.nz; Tel.: +64-3-423-0665
Received: 26 May 2020; Accepted: 15 August 2020; Published: 19 August 2020
Animals 2020, 10, 1447; doi:10.3390/ani10091447
Table 1. Associations between milk fatty acid methyl ester levels and
myostatin genotypes in NZ HF * J cross cows (n = 347).
Individual/Grouped Mean [+ or -] Standard Error (2) (g/100 g
Milk Fatty Acids)
Fatty Acids (1) AA (n = 149)
C4:0 1.264 [+ or -] 0.012
C6:0 1.566 [+ or -] 0.010
C8:0 1.198 [+ or -] 0.009
C10:0 3.299 [+ or -] 0.035
C10:l 0.294 [+ or -] 0.004
C11:0 0.062 [+ or -] 0.002
C12:0 4.029 [+ or -] 0.046 (a)
C13:0 iso 0.084 [+ or -] 0.002
C12:l 0.096 [+ or -] 0.002 (a)
C13:0 anteiso 0.038 [+ or -] 0.000
C13:0 0.123 [+ or -] 0.002
C14:0 12.513 [+ or -] 0.079
C14:l cis-9 1.003 [+ or -] 0.020
C15:0 iso 0.296 [+ or -] 0.002
C15:0 anteiso 0.639 [+ or -] 0.009
C15:0 1.485 [+ or -] 0.016
C15:1 0.283 [+ or -] 0.003
C16:0 37.396 [+ or -] 0.288
C16:1 cis-9 1.266 [+ or -] 0.023
C17:0 iso 0.555 [+ or -] 0.006
C17:0 0.690 [+ or -] 0.006
C17:1 0.195 [+ or -] 0.002
C18:0 8.547 [+ or -] 0.119
C18:1 trans-11 2.756 [+ or -] 0.068
C18:1 cis-9 12.833 [+ or -] 0.142
C18:2 trans-9,12 0.377 [+ or -] 0.006
C18:2 cis-9, trans-13 0.288 [+ or -] 0.003
C18:2 cis-9, trans-12 0.067 [+ or -] 0.002
C18:2 trans-9, cis-12 0.475 [+ or -] 0.011
C18:2 cis-9,12 0.693 [+ or -] 0.007
C19:0 0.140 [+ or -] 0.003
C18:3 cis-6,9,12 0.075 [+ or -] 0.001
C18:3 cis-9,12,15 0.805 [+ or -] 0.010
CLA cis-9, trans-11 1.012 [+ or -] 0.028
C20:0 0.126 [+ or -] 0.002
C20:1 cis-5 0.060 [+ or -] 0.001
C20:1 cis-9 0.151 [+ or -] 0.002
C20:1 cis-11 0.077 [+ or -] 0.001
C20:3 cis-8,11,14 0.030 [+ or -] 0.001
C20:4 cis-5,8,11,14 0.035 [+ or -] 0.001
C22:0 0.064 [+ or -] 0.001
C22:1 trans-13 0.066 [+ or -] 0.001
C20:5 cis-5,8,11,14,17 0.088 [+ or -] 0.00
C24:0 0.044 [+ or -] 0.001
C22:5 cis-7,10,13,16,19 0.123 [+ or -] 0.002
SCFA 4.030 [+ or -] 0.027
MCFA 19.841 [+ or -] 0.148
LCFA 48.492 [+ or -] 0.252
Total C18:1 16.300 [+ or -] 0.168
Total C18:2 2.911 [+ or -] 0.041
Total C18:3 0.880 [+ or -] 0.011
Omega 3 1.017 [+ or -] 0.011
Omega 6 0.832 [+ or -] 0.008
MUFA 19.792 [+ or -] 0.174
PUFA 4.067 [+ or -] 0.044
Branched FA 1.611 [+ or -] 0.014
Total UFA 23.859 [+ or -] 0.206
Total SFA 72.548 [+ or -] 0.225
Unsaturated index 24.756 [+ or -] 0.218
C10:1 index 8.252 [+ or -] 0.136
C12:1 index 2.328 [+ or -] 0.034
C14:1 index 7.424 [+ or -] 0.148
C16:1 index 3.281 [+ or -] 0.052
C18:1 index 65.594 [+ or -] 0.331
CLA index 26.626 [+ or -] 0.269
Individual/Grouped
Fatty Acids (1) AB (n = 87)
C4:0 1.252 [+ or -] 0.015
C6:0 1.548 [+ or -] 0.013
C8:0 1.179 [+ or -] 0.011
C10:0 3.244 [+ or -] 0.043
C10:l 0.282 [+ or -] 0.005
C11:0 0.058 [+ or -] 0.002
C12:0 3.939 [+ or -] 0.057 (ab)
C13:0 iso 0.079 [+ or -] 0.002
C12:l 0.091 [+ or -] 0.002 (ab)
C13:0 anteiso 0.038 [+ or -] 0.001
C13:0 0.119 [+ or -] 0.003
C14:0 12.478 [+ or -] 0.098
C14:l cis-9 0.959 [+ or -] 0.025
C15:0 iso 0.290 [+ or -] 0.003
C15:0 anteiso 0.633 [+ or -] 0.011
C15:0 1.473 [+ or -] 0.019
C15:1 0.281 [+ or -] 0.004
C16:0 37.624 [+ or -] 0.354
C16:1 cis-9 1.245 [+ or -] 0.028
C17:0 iso 0.553 [+ or -] 0.008
C17:0 0.698 [+ or -] 0.008
C17:1 0.194 [+ or -] 0.003
C18:0 8.652 [+ or -] 0.146
C18:1 trans-11 2.782 [+ or -] 0.083
C18:1 cis-9 12.891 [+ or -] 0.175
C18:2 trans-9,12 0.370 [+ or -] 0.007
C18:2 cis-9, trans-13 0.281 [+ or -] 0.004
C18:2 cis-9, trans-12 0.064 [+ or -] 0.002
C18:2 trans-9, cis-12 0.483 [+ or -] 0.014
C18:2 cis-9,12 0.677 [+ or -] 0.009
C19:0 0.140 [+ or -] 0.003
C18:3 cis-6,9,12 0.073 [+ or -] 0.001
C18:3 cis-9,12,15 0.776 [+ or -] 0.013
CLA cis-9, trans-11 1.001 [+ or -] 0.034
C20:0 0.126 [+ or -] 0.002
C20:1 cis-5 0.060 [+ or -] 0.002
C20:1 cis-9 0.152 [+ or -] 0.003
C20:1 cis-11 0.075 [+ or -] 0.002
C20:3 cis-8,11,14 0.031 [+ or -] 0.001
C20:4 cis-5,8,11,14 0.035 [+ or -] 0.001
C22:0 0.067 [+ or -] 0.002
C22:1 trans-13 0.068 [+ or -] 0.002
C20:5 cis-5,8,11,14,17 0.089 [+ or -] 0.001
C24:0 0.046 [+ or -] 0.001
C22:5 cis-7,10,13,16,19 0.122 [+ or -] 0.003
SCFA 3.979 [+ or -] 0.033
MCFA 19.661 [+ or -] 0.182
LCFA 48.826 [+ or -] 0.310
Total C18:1 16.383 [+ or -] 0.206
Total C18:2 2.876 [+ or -] 0.051
Total C18:3 0.849 [+ or -] 0.013
Omega 3 0.987 [+ or -] 0.014
Omega 6 0.815 [+ or -] 0.010
MUFA 19.790 [+ or -] 0.214
PUFA 4.002 [+ or -] 0.054
Branched FA 1.593 [+ or -] 0.017
Total UFA 23.792 [+ or -] 0.254
Total SFA 72.643 [+ or -] 0.277
Unsaturated index 24.678 [+ or -] 0.269
C10:1 index 8.032 [+ or -] 0.167
C12:1 index 2.239 [+ or -] 0.042
C14:1 index 7.133 [+ or -] 0.183
C16:1 index 3.200 [+ or -] 0.063
C18:1 index 65.409 [+ or -] 0.407
CLA index 26.339 [+ or -] 0.332
Individual/Grouped p (3)
Fatty Acids (1) AC (n = 48)
C4:0 1.252 [+ or -] 0.020 0.871
C6:0 1.549 [+ or -] 0.017 0.458
C8:0 1.188 [+ or -] 0.014 0.150
C10:0 3.293 [+ or -] 0.056 0.097
C10:l 0.283 [+ or -] 0.007 0.108
C11:0 0.056 [+ or -] 0.003 0.124
C12:0 4.015 [+ or -] 0.075 (ab) 0.046
C13:0 iso 0.080 [+ or -] 0.002 0.076
C12:l 0.093 [+ or -] 0.003 (ab) 0.019
C13:0 anteiso 0.038 [+ or -] 0.001 0.450
C13:0 0.117 [+ or -] 0.004 0.364
C14:0 12.690 [+ or -] 0.128 0.266
C14:l cis-9 0.967 [+ or -] 0.033 0.390
C15:0 iso 0.299 [+ or -] 0.004 0.148
C15:0 anteiso 0.647 [+ or -] 0.014 0.814
C15:0 1.474 [+ or -] 0.026 0.922
C15:1 0.287 [+ or -] 0.005 0.714
C16:0 36.999 [+ or -] 0.465 0.711
C16:1 cis-9 1.273 [+ or -] 0.037 0.576
C17:0 iso 0.556 [+ or -] 0.010 0.987
C17:0 0.693 [+ or -] 0.010 0.816
C17:1 0.198 [+ or -] 0.004 0.639
C18:0 8.652 [+ or -] 0.192 0.918
C18:1 trans-11 2.801 [+ or -] 0.110 0.974
C18:1 cis-9 12.891 [+ or -] 0.230 0.595
C18:2 trans-9,12 0.391 [+ or -] 0.009 0.230
C18:2 cis-9, trans-13 0.291 [+ or -] 0.005 0.225
C18:2 cis-9, trans-12 0.065 [+ or -] 0.003 0.670
C18:2 trans-9, cis-12 0.492 [+ or -] 0.018 0.771
C18:2 cis-9,12 0.697 [+ or -] 0.012 0.279
C19:0 0.143 [+ or -] 0.004 0.890
C18:3 cis-6,9,12 0.074 [+ or -] 0.002 0.511
C18:3 cis-9,12,15 0.818 [+ or -] 0.017 0.144
CLA cis-9, trans-11 1.031 [+ or -] 0.045 0.920
C20:0 0.130 [+ or -] 0.003 0.530
C20:1 cis-5 0.062 [+ or -] 0.002 0.891
C20:1 cis-9 0.149 [+ or -] 0.004 0.890
C20:1 cis-11 0.076 [+ or -] 0.002 0.093
C20:3 cis-8,11,14 0.032 [+ or -] 0.001 0.137
C20:4 cis-5,8,11,14 0.034 [+ or -] 0.001 0.927
C22:0 0.066 [+ or -] 0.002 0.519
C22:1 trans-13 0.070 [+ or -] 0.002 0.190
C20:5 cis-5,8,11,14,17 0.089 [+ or -] 0.002 0.936
C24:0 0.044 [+ or -] 0.001 0.294
C22:5 cis-7,10,13,16,19 0.123 [+ or -] 0.004 0.996
SCFA 3.989 [+ or -] 0.043 0.451
MCFA 19.997 [+ or -] 0.239 0.127
LCFA 48.201 [+ or -] 0.407 0.587
Total C18:1 16.400 [+ or -] 0.271 0.634
Total C18:2 2.969 [+ or -] 0.067 0.664
Total C18:3 0.893 [+ or -] 0.017 0.123
Omega 3 1.030 [+ or -] 0.018 0.173
Omega 6 0.838 [+ or -] 0.013 0.317
MUFA 19.858 [+ or -] 0.281 0.657
PUFA 4.139 [+ or -] 0.071 0.410
Branched FA 1.620 [+ or -] 0.022 0.735
Total UFA 23.997 [+ or -] 0.334 0.663
Total SFA 72.362 [+ or -] 0.364 0.701
Unsaturated index 24.911 [+ or -] 0.353 0.677
C10:1 index 7.972 [+ or -] 0.219 0.410
C12:1 index 2.266 [+ or -] 0.055 0.243
C14:1 index 7.087 [+ or -] 0.240 0.427
C16:1 index 3.319 [+ or -] 0.083 0.435
C18:1 index 65.514 [+ or -] 0.535 0.855
CLA index 26.725 [+ or -] 0.436 0.716
(1) SCFA- short-chain length fatty acid; MCFA--medium-chain length
fatty acid; LCFA--long-chain length fatty acid; MUFA--monounsaturated
fatty acid; PUFA--polyunsaturated fatty acid; UFA--unsaturated fatty
acid; SFA--saturated fatty acid. (2) Predicted means and standard error
of those means derived from general linear mixed-effects models (GLMM).
Myostatin genotype (categorical variable), cow age (categorical), herd
(categorical), and days in milk (continuous) were fitted to the model
as fixed effects. Means within a row that do not share a superscript
letter are separated at p < 0.05. (3) 0.05 < p < 0.2 in italics, while
p < 0.05 in bold.
 Printer friendly
 Cite/link
 Email
 Feedback
| |
| Title Annotation: | Article |
|---|---|
| Author: | Haruna, Ishaku L.; Li, Yunhai; Ekegbu, Ugonna J.; Amirpour-Najafabadi, Hamed; Zhou, Huitong; Hickfor |
| Publication: | Animals (Basel) |
| Article Type: | Report |
| Geographic Code: | 8NEWZ |
| Date: | Sep 1, 2020 |
| Words: | 6863 |
| Previous Article: | Development and Evaluation of a Disease Large Animal Model for Preclinical Assessment of Renal Denervation Therapies. |
| Next Article: | A Randomized, Placebo-Controlled Clinical Trial of Famciclovir in Shelter Cats with Naturally Occurring Upper Respiratory Tract Disease. |
| Topics: | |