Next Article in Journal
Modest Ag-Extension and Access to Seeds of Aromatic Rice Can Boost Returns of Smallholder Farmers in Uganda, A Case Study
Next Article in Special Issue
Comparative Transcriptome Analysis Reveals the Key Genes Involved in Lipid Deposition in Pekin Ducks (Anas platyrhynchos domesticus)
Previous Article in Journal
An Algorithm for Calculating Apple Picking Direction Based on 3D Vision
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Single Nucleotide Polymorphisms in the Fatty Acid Binding Protein 4, Fatty Acid Synthase and Stearoyl-CoA Desaturase Genes Influence Carcass Characteristics of Tropical Crossbred Beef Steers

by
Felista W. Mwangi
1,
Shedrach B. Pewan
1,2,
John R. Otto
1,
Oyelola A. Adegboye
3,
Edward Charmley
4,
Christopher P. Gardiner
1,
Bunmi S. Malau-Aduli
5,
Robert T. Kinobe
1 and
Aduli E. O. Malau-Aduli
1,*
1
Animal Genetics and Nutrition, Veterinary Sciences Discipline, College of Public Health, Medical and Veterinary Sciences, Division of Tropical Health and Medicine, James Cook University, Townsville, QLD 4811, Australia
2
National Veterinary Research Institute, Vom Private Mail Bag 01, Plateau State, Nigeria
3
Public Health and Tropical Medicine Discipline, College of Public Health, Medical and Veterinary Sciences, Division of Tropical Health and Medicine, James Cook University, Townsville, QLD 4811, Australia
4
Commonwealth Scientific and Industrial Research Organisation, Agriculture and Food, Australian Tropical Sciences and Innovation Precinct, James Cook University, Townsville, QLD 4811, Australia
5
College of Medicine and Dentistry, Division of Tropical Health and Medicine, James Cook University, Townsville, QLD 4811, Australia
*
Author to whom correspondence should be addressed.
Agriculture 2022, 12(8), 1171; https://doi.org/10.3390/agriculture12081171
Submission received: 13 June 2022 / Revised: 2 August 2022 / Accepted: 4 August 2022 / Published: 6 August 2022
(This article belongs to the Special Issue Advances in Molecular Genetics in Domestic Animals)

Abstract

:
This study explored the identification of single nucleotide polymorphisms (SNP) in fatty acid binding protein 4 (FABP4), stearoyl-CoA desaturase (SCD), and fatty acid synthase (FASN) genes that may influence the carcass traits of tropical crossbred beef cattle. The hypothesis tested was that SNP in the FABP4, SCD, and FASN genes are associated with chiller-assessed carcass traits of tropically adapted northern Australian crossbred beef cattle. Fifty Bos indicus and Bos taurus crossbred steers were backgrounded on either buffel grass only, or buffel grass and desmanthus mixed pastures for 147 days and finished in a commercial feedlot for 110 days. Steers were slaughtered within 48 h of leaving the feedlot within a lairage period not exceeding 12 h and carcasses graded 12 h after slaughter. Next-generation sequencing of the FASN, FABP4, and SCD genes identified multiple SNP loci that were correlated and significantly associated with carcass traits. The FABP4 g.44677205A>G locus was significantly associated with hump height and correlated with loin eye muscle area (EMA; p < 0.05). Polymorphism in the SCD gene g.21275851C>A locus was associated with subcutaneous fat depth and marbling score (p < 0.05). The CC genotype had a higher subcutaneous fat depth and marbling score (p < 0.05) than the AA genotype. Significant correlations were observed between carcass marbling score and subcutaneous fat depth within the FASN SNP locus (p < 0.05). Therefore, the hypothesis that SNP in the FABP4, SCD, and FASN genes are associated with chiller-assessed carcass traits of tropically adapted northern Australian crossbred beef cattle was accepted. These findings suggest that SNP in the FABP4, SCD, and FASN genes may be used in carcass grading and meat quality improvement through marker-assisted selection of northern Australian crossbred beef cattle.

1. Introduction

Observable and measurable carcass traits are important in assigning meat quality values. The postmortem measurements of subcutaneous fat depth and marbling are useful in carcass grading as indirect meat quality indicators [1,2]. Characteristics such as loin eye muscle area (EMA) and subcutaneous fat depth are indicative of the amount of saleable carcass, while marbling score is a vital indicator of meat eating quality [3]. High beef quality standards are vital for consumer satisfaction, influencing the decision to re-purchase and maintain or increase the global market share of export-dependent beef industries [4,5]. Tropical beef cattle breeds (Bos indicus) are predominant in northern Australia (where over half of the Australian beef cattle herd is reared), due to their ability to adapt and survive in harsh environments characterized by poor feed quality, high ambient temperatures, and high parasite load [6]. However, tropical beef cattle have a low growth rate and produce comparatively tougher meat with lower fat content than the taurine breeds. Thus, they are crossed with the taurine breeds to improve growth rate and meat quality without compromising their ability to survive in harsh environments [7]. Therefore, the measure of tropical breed content, as determined by the hump height, is critical during carcass assessment [8].
Carcass traits are influenced by age, genetics, and management; hence, selective breeding to achieve long-term enhancement of economically important carcass traits is a relevant tool [3]. However, it is difficult to attain efficient genetic gain using traditional breeding methods since most meat quality measurements are obtained after animal slaughter, thus making it difficult to determine meat quality in a living animal [9,10]. This creates a substantial genetic lag when progeny phenotypic carcass traits are used to select breeding sires [3]. Advances in molecular genetics have led to the identification of candidate genes influencing meat quality traits [11], and single nucleotide polymorphisms (SNP) loci in genes have been evaluated for their potential use as genetic markers for marker-assisted selection aimed at improving meat quality [12,13]. The fatty acid binding protein 4 (FABP4), fatty acid synthase (FASN), and stearoyl-CoA desaturase (SCD) genes are among the candidate genes associated with carcass traits. The bovine FABP4 gene encodes a cytoplasmic protein that binds long-chain fatty acids (FA) and other hydrophobic ligands [14]. The protein is involved in the regulation of lipid hydrolysis and intracellular fatty acids trafficking [14]; as a result, FABP4 is considered a functional and positional candidate gene for fat content and distribution rate in the muscle [15,16]. Associations of SNP in the FABP4 gene with marbling score and hot carcass weight in Holstein bulls [10], carcass weight, marbling score, and meat quality grade in Hanwoo cattle [17,18], and marbling score in Yanbian yellow cattle [19] have been reported. FASN encodes a multifunctional polypeptide of enzymes associated with fatty acid biosynthesis in the cytosol of animal cells [20] and is considered a candidate gene for marker-assisted selection and breeding for intramuscular fat (IMF) for improved meat quality and carcass grades [21]. Several FASN SNP have been associated with backfat thickness, loin eye muscle area (EMA), and IMF content in Qinchuan cattle [3] and with hot carcass weight in Beefmaster, Brangus, Bonsmara, Romosinuano, Hereford, and Angus beef cattle [22]. Polymorphisms of the SCD gene have been reported to be associated with marbling score, but not backfat thickness in Wagyu x Limousin cattle [23]. In Chinese Simmental cattle, SNP were associated with IMF, but not marbling score [24]. In contrast, neither backfat thickness nor IMF in Spanish commercial beef cattle [25] were associated with any SNP. Barendse et al. [16] reported that SCD g.2502C>G SNP was associated with muscle IMF in Angus, but not Hereford, Shorthorn, Murray Grey, Belmont Red, Brahman, and Santa Gertrudis cattle. These findings indicate that SNP-based genetic markers need to be tested and confirmed for the specific breeds of interest [26], hence the need for a targeted search for SNP in FABP4, SCD, and FASN genes in northern Australian crossbred tropical beef cattle and to determine their association with carcass traits. Therefore, this study aimed to detect SNP in the exons and introns of FABP4, SCD, and FASN genes of northern Australian tropical crossbred beef cattle and investigate their associations with carcass traits. It was hypothesised that SNP in the FABP4, SCD, and FASN genes are associated with chiller-assessed carcass traits of tropically adapted northern Australian crossbred beef cattle.

2. Materials and Methods

This study followed the Australian code of practice for the care and use of animals for scientific purposes [27] and was approved by the James Cook University Animal Ethics Committee (Approval Number 2639).

2.1. Animal Management

Fifty tropical steers of crossbred Bos indicus and Bos taurus breeds were randomly selected and used for this study based on an a priori G-Power analysis to achieve 80% power with a 4.0 critical F-value for a large effect size and α of 0.05. Steers were backgrounded on buffel grass only or buffel grass and desmanthus pastures for 147 days as described previously [28] and finished in a commercial feedlot for 110 days. The steers were managed and transported according to approved Meat Standards Australia (MSA) protocols [29] and slaughtered in a commercial abattoir within 48 h of leaving the feedlot with a lairage period not exceeding 12 h.

2.2. Carcass Evaluation

Carcasses were graded according to the Aus-Meat and MSA (Meat Standards Australia) grading standards [30] by an accredited meat grader at a commercial export-licensed abattoir in Queensland, Australia. The traits of interest were hot standard carcass weight (HCW), hump height, marbling score, subcutaneous rib fat (back fat) thickness, rump fat thickness at the P8 site, eye muscle area (EMA), meat colour, and MSA index. Carcass assessment was carried out at the 12th–13th rib interface 12 h after slaughter. The hump height, an indicator of tropical breed content of the carcass, was measured at the position of the greatest hump width by placing a ruler parallel to the surface of the sawn chine perpendicular to the first thoracic vertebrae and recorded in millimetres. Carcass marbling indicates the level of intramuscular fat content and was scored using the 100 (devoid) to 1100+ (abundant) reference standards, while subcutaneous fat thickness at the 12th rib and P8 site were measured in millimetres using a ruler grid. The EMA was the M. longissimus dorsi surface area at the 12th rib calculated in square centimetres and measured using a plastic grid. The meat colour was scored at the bloomed loin eye muscle using the 1 (light) to 7 (dark) scale colour reference standards. The MSA Index was calculated as the sum of the predicted eating quality scores for 39 MSA cuts weighted by their relative proportion of total carcass weight.

2.3. Blood Collection and Genomic DNA Extraction

Blood was collected at the start of the backgrounding phase via caudal venipuncture into 4 mL EDTA-containing vacutainer tubes (BD, Sydney, Australia) and stored at −80 °C prior to analysis. Blood samples were later thawed at room temperature and an aliquot of 2 mL was used for genomic DNA extraction using NucleoSpin Blood Kit (Macherey-Nagel GmbH and Co. KG, Duren, Germany) according to the manufacturer’s protocol. DNA yield and purity were quantified using a NanoDrop ND-1000 (Thermo Fisher Scientific Australia Pty Ltd., Victoria, Australia) and diluted to 50 ng/μL with nuclease-free water.

2.4. Primer Design

All primers were designed using the Geneious Prime software version 2020.2.5 (Biomatters Ltd., Auckland, New Zealand) and synthesised at Integrated DNA Technologies Pte Ltd. (Queenstown, Singapore). Bos taurus (Hereford) sequences of FASN, FABP4, and SCD gene were obtained from the National Center for Biotechnology Information database (GenBank) and used as references. The 15 kb SCD sequence (NC_ 037353.1) was split into three fragments and the 18 kb long FASN (NC_ 037346.1) sequence was split into two overlapping fragments, while the 4 kb FABP4 (NC_ 037341.1) gene was not fragmented before amplification (Table 1, Supplementary Figure S1).

2.5. Target Gene Amplification

Gene amplification, library preparation, normalization, sequencing, and data analysis procedures were carried out as described previously [31]. In brief, three diverse long-range polymerase chain reaction (PCR) approaches were tested to amplify the SCD, FASN, and FABP4 genes, and the optimum combinations were selected. Platinum SuperFi PCR Master Mix and Phusion Hot Start II High-Fidelity PCR Master Mix (Thermo Fisher Scientific, Scoresby, Australia) were used to amplify the FABP4 and SCD genes, respectively, under similar PCR conditions. The amplification reactions were executed in a SimpliAmp Thermal Cycler (Thermofisher Scientific, Scoresby, Australia) in a total volume of 50 µL consisting of 25 µL of PCR master mix, 100 ng of DNA template, and 0.5 µM of each primer in a 3-step procedure. The conditions were as follows: single initial denaturation at 98 °C for 1 min, 35 cycles of denaturation, annealing, and extension at 98 °C for 15 s, 60 °C for 15 s, and 72 °C for 9 min, respectively, followed by a final extension at 72 °C for 9 min and a 4 °C hold. The FASN gene fragments were amplified using PrimeSTAR GXL Master Mix (TaKaRa Bio Inc., Kusatsu, Shiga, Japan) in a 2-step protocol. The amplification reaction mix consisted of 1.25 units of polymerase, 10 µL of 5× buffer, 0.2 µM of each primer, 200 µM of dNTP mixture, and 100 ng of DNA template in a total volume of 50 µL. The amplification reaction conditions included initial denaturation for 1 min at 98 °C and 30 cycles of denaturation and annealing combined with extension at 98 °C for 10 s and 68 °C for 10 min, respectively. The amplification products were visualised using gel electrophoresis in 0.8% agarose gel to evaluate success of the assay.

2.6. PCR Products Clean-Up

All PCR products were cleaned using the Sera-Mag SpeedBeads in a Zephyr NGS Workstation (PerkinElmer, Waltham, MA, USA) and quantified using the Promega dsDNA Quantifluor System Kit (Ref: E2670, 00002484139) on the PerkinElmer Enspire Workstation. The six amplification products from the three genes were pooled at 0.4 nM, normalised to 2 ng/µL and then diluted to 0.2 ng/µL with 10 mM Tris-HCl at pH 8.0 for library preparation.

2.7. Library Preparation and Sequencing

Libraries were prepared with the Nextera XT DNA Library Prep kit (Illumina, CA, USA) according to the manufacturer’s instructions. Libraries were purified using 0.6× Sera-Mag SpeedBeads (Merck KGaA, Darmstadt, Germany) and washed twice with 80% ethanol to remove unincorporated adapters and fragments shorter than 250 bp length. The DNA libraries were analysed for fragment size and concentration using Agilent High Sensitivity D5000 reagents and ScreenTape on the Tape Station 4200 (Agilent Technologies, Santa Clara, CA, USA) and further concentration determined using the QuantiFluor® dsDNA System (Promega, Madison, WI, USA). The fragment size and concentration data were used to normalise each library to 2 nM with 10 mM and pH 8.5 Tris-HCl, and samples of the different steers were then pooled together. A 10 pM input of the pooled samples and 10% PhiX spike-in were used for sequencing on the Illumina MiSeq benchtop sequencer with a 500-cycle MiSeq Reagent Nano Kit v2 (Illumina, Inc, San Diego, CA, USA).

2.8. Sequence Data Analysis and Calculations

Sequence data were downloaded from the Illumina Dashboard-BaseSpace Sequence Hub (https://basespace.illumina.com/dashboard, accessed on 27 July 2021) and analysed using the Geneious Prime software with NC_ 037353.1, NC_ 037346.1 and NC_ 037341.1 as the SCD, FASN, and FABP4 reference sequences, respectively. The retrieved reads were trimmed, and adapters were removed using the BBDuk trimmer. The Phred quality score was set at 20 to increase the probability of calling true SNP to 99%. All short reads of 20 bp or less were discarded and low-coverage regions were excluded during SNP calling. The allele and genotype frequencies were determined by direct counting, while the polymorphism information content (PIC) was determined using the GDIcall online calculator (http://www.msrcall.com/Gdicall.aspx) (accessed on 11 January 2022). Hardy–Weinberg equilibrium (HWE) and expected heterozygosity (He) were calculated as described previously [32].

2.9. Statistical Analysis

Statistical data analyses were carried out using the R software version 4.0.2 (R Foundation for Statistical Computing, Vienna, Austria). Initial summary statistics were computed for means, standard deviations, distribution, and range to identify any outliers. The chi-square test was used to determine if the SNP were in HWE. Distance-based hierarchical clustering of the SNP loci was used to examine the degree of linkage disequilibrium between each pair of loci [33] and the results were presented as dendrograms and heatmaps. Spearman’s correlation analysis was used to estimate the strength of relationships between genomic variants and carcass traits. Linear models were fitted to investigate whether a gene SNP was associated with carcass traits using generalized least squares linear procedures to compute group means, while significant differences between least-square means were compared using the Tukey-adjusted multiple comparisons test. Differences were declared significant at p < 0.05.

3. Results

3.1. Genetic Variants and Population Diversity

In a targeted evaluation of the associations between SNP loci in the FABP4, SCD, and FASN genes with chiller-assessed carcass traits of tropical crossbred beef cattle, a total of 65 SNP loci were identified. These comprise 11, 27, and 27 SNP for the FABP4, SCD, and FASN genes, respectively (Supplementary Table S1). Twenty-eight SNP were considered novel since they were not found in the Bovine Genome Variation Database (BGVD) (http://animal.nwsuaf.edu.cn/code/index.php/BosVar, accessed on 28 January 2022). Only seven SNP were predicted to be non-synonymous amino acid substitutions (Table 2). The SNP were in HWE except for g.50786977A>G and g.50790973C>A in the FASN gene (Supplementary Table S1; p ≤ 0.04). The MAF, He, and PIC for the SNP ranged between 0.11–0.47, 0.20–0.49, and 0.18–0.37, respectively, for all three genes. Correlation coefficients were used to measure the degree of linkage disequilibrium between each pair of loci depicted in Supplementary Figures S2–S4 and showed the presence of local patterns forming four clusters for each gene.

3.2. Correlations between SNP and Carcass Traits

The clustering patterns of the SNP loci are shown in the top panels of Figure 1A, Figure 2A and Figure 3A. The level of SNP heterozygosity and homozygosity varied widely between closely related individuals in the SCD gene compared to both FABP4 and FASN genes. The correlation coefficients between the SNP and carcass traits are presented in the bottom panels of Figure 1B, Figure 2B and Figure 3B. Most of the SNP loci in the FABP4 gene were in linkage disequilibrium (p < 0.05), but the correlations between SNP and carcass traits were observed in only three SNP. The g.44677205A>G (rs109388335), g.44677611G>C (rs41729172), and g.44679833A>G (rs133333024) SNP were all positively correlated with hump height (p < 0.05), while g.44677205A>G (rs109388335) was negatively correlated with EMA (p < 0.05). All SNP identified in the SCD gene were in linkage disequilibrium (p < 0.05), except g.21267896C>T (rs136334180) that was not correlated with g.21274479G>A (rs382184952) (Figure 2B). Polymorphisms on the g.21271392G>A (rs211294052) and g.21274479G>A (rs382184952) loci were positively correlated with HCW and EMA, respectively (p < 0.05). Moreover, g.21275851C>A (novel) was negatively correlated with P8 fat thickness, backfat thickness, and fat class (p < 0.01), while g.21273692T>C (rs208058585), g.21276141C>T (rs41255697), and g.21276672A>G (rs41255698) were negatively correlated with carcass marbling score (p < 0.05). Most SNP detected in the FASN gene were in linkage disequilibrium (Figure 3B; p < 0.05). Positive correlations were observed in g.50784533C>G (rs481622676) and g.50792445C>T (novel) SNP with carcass P8 fat and backfat thickness, while g.50786221A>G (rs518879624) and g.50787362C>T (novel) were negatively correlated with carcass marbling (p < 0.05). Several carcass traits were significantly correlated with each other. For instance, HCW was positively correlated (≥30) with P8 fat, backfat thickness, EMA, and MSA index, but negatively correlated with meat colour scores (p < 0.05). A highly positive correlation (≥80) between P8 fat with back fat and fat class was observed (p < 0.001). Hump height was negatively correlated with MSA index (p < 0.001), but there was no correlation between hump height and other carcass traits.

3.3. Associations between SNP and Carcass Traits

Associations between FABP4 g.44677205A>G (rs109388335), SCD g.21275851C>A (novel) and FASN g.50784533C>G (rs481622676) SNP with carcass traits are presented in Table 3. The FABP4 g.44677205A>G SNP had a significant association with hump height (p < 0.01) and tended to be associated with EMA (p = 0.05). Multiple genotype comparisons at the g.44677205A>G locus showed that the GG genotype had a higher hump height (125.8 ± 10.68 mm) compared to the AA (109.3 ± 13.39 mm) genotype (p ≤ 0.01), but similar to the AG genotype (p = 0.08). The hump height difference between the AA and AG genotypes was insignificant (p = 0.64). EMA was significantly higher for the AA (88.4 ± 8.71 cm2) than the GG (81.3 ± 3.88 cm2) and AG (86.7 ± 4.93 cm2) genotypes (Figure 4; p ≤ 0.03), but genotypic differences between AA and AG were not significant (p = 0.81). The AA genotype of SCD g.21275851C>A SNP had the lowest P8 fat, marbling, backfat, and fat class, while the highest scores were observed for the CC genotype (p < 0.02). Multiple comparisons between SCD g.21275851C>A SNP variants and carcass traits indicated that P8 fat, marbling, backfat, and fat class scores were significantly different between the CC and AA genotypes (Figure 5; p < 0.03). No significant associations were observed between the FASN g.50784533C>G variants and the carcass traits measured (p ≥ 0.10).

4. Discussion

Chiller-assessed carcass quality measurements are obtained after slaughter, hence the data cannot be used to make management decisions such as culling the inferior and breeding the superior performing animals [9,10]. Regions on chromosomes 1, 9, 14, 16, 19, 23, 26, and 29, had been previously identified to be associated with fatty acid composition and carcass traits in a cohort of Bos indicus, Bos taurus, and tropical composite beef cattle using high-density data [34]. Other studies reported that single allelic substitutions influence carcass traits and can be used to predict carcass quality in living animals [17,23]. In this study, a targeted next-generation sequencing technique was used to identify SNP in the FABP4, SCD, and FASN genes that may be used as molecular markers for carcass quality selection in northern Australian tropically adapted beef cattle. A marker PIC is one of the indicators of marker quality and it reflects the ability of a marker to detect polymorphisms between individuals in a population [35]. Markers with PIC values above 0.5 are deemed very informative, 0.25–0.50 are a bit informative, and below 0.25 are minimally informative [36]. In this study, all the markers were informative with the exception of only one marker at the g.50783803G>A locus that had a PIC below 0.25.

4.1. FABP4 Gene Polymorphisms

The FABP4 gene is a primary metabolic indicator of IMF deposition as it is located within the quantitative trait loci region that contributes to serum leptin, a protein involved in body fat regulation, and FABP4 also encodes a protein involved in intracellular fatty acids trafficking [37,38]. Damon et al. [39] reported that FABP4 protein level was positively correlated with fat cell count and lipid content in porcine. Genetic variants of the FABP4 gene in cattle and their associations with meat and carcass traits have been reported [14,16,18,40]. The FABP4 g.3691G>A SNP was associated with marbling score, while g.2834C>G was reported to be associated with HCW in Holstein bulls [10]. The g.3473T>A SNP was reported to have a significant effect on carcass weight, while g.3631A>G significantly influenced marbling score in Hanwoo cattle [17]. In another study, FABP4 g.3691G>A SNP was reported to influence marbling score and meat quality grade in Hanwoo cattle [18]. Five SNP (g.3496A>C, g.3745T>C, g.3533A>T, g.3767T>C, and g.3711G>C) were associated with marbling scores in Yanbian yellow cattle [19]. However, these polymorphisms were not observed in this study. Furthermore, there were no correlations or associations between the identified FABP4 SNP in the present study and carcass subcutaneous fat depth (P8 fat and backfat) or marbling score. These findings agree with a previous study [38] that reported no significant effect of g.44677959T>C (c.220) with subcutaneous fat nor marbling score in Japanese Black cattle. In contrast, Cho et al. [41] reported an SNP association with backfat thickness in Korean cattle, while Goszczynski et al. [2] reported an additive effect in cattle of varying breeds. Three SNP (g.44677205A>G, g.44677611G>C, and g44679833A<G) were positively correlated with hump height in this study, with the homozygous GG genotype having the highest and AA the lowest measurements. Hump height is used to estimate the amount of tropical breed content of the carcass, with high values indicating higher tropical breed content and, subsequently, lower meat eating quality [42]. Therefore, the correlations indicate that genetic variants of the three loci are influenced by breed, in agreement with previous studies [18,43]. Both EMA and subcutaneous fat depth indicate the amount of saleable meat from the carcass [3]. The negative correlation between g.44677205A>G and EMA, and the lack of correlation with subcutaneous fat depth, indicate that selecting for the A allele may increase the amount of carcass saleable meat without increasing the subcutaneous fat that ends up being trimmed off [44].

4.2. SCD Gene Polymorphisms

The SCD gene is highly expressed in lipogenic tissues and encodes for an enzyme that desaturates SFA by introducing a cis double bond at the 9th and 10th carbon interface of the fatty acid [45]. The unsaturation of the fatty acid chain is a key determining factor of the melting temperature of triacylglycerols, thus influencing meat fat hardness [44]. Alleles associated with high SCD enzyme activity have been reported to be associated with increased marbling scores in Wagyu × Limousin cattle [23]. Similarly, three SNP were correlated with carcass marbling score in this study, where the C, T, and G alleles of the g.21273692T>C, g.21276141C>T, and g.21276672A>G loci, respectively, had favourably higher marbling scores. The missense mutation that causes substitution of valine (type V) to alanine (type A) at PCS position 878 (g.21272422C>T) had no effect on carcass traits measured in this study. These findings align with previous studies that reported no effect on backfat thickness in Spanish commercial beef and Chinese Simmental cattle [24,25]. In addition, the SNP was not associated with HCW, EMA, backfat thickness, and marbling score of Japanese black cattle reported by Ohsaki et al. [46]. More studies are required to determine the mechanism behind the effect of SCD SNP reported in the Wagyu x Limousin cattle [23] and lack of it in other cattle breeds [24,25,46]. Accordingly, the study of Li et al. [11] on the effect of SNP on carcass traits of beef cattle populations indicated no significant effect on marbling score and meat colour soon after slaughter. Genes responsible for fat deposition and metabolism are reported to influence meat colour since high fat levels accelerate lipid myoglobin oxidation and, subsequently, meat discoloration, since marbling influences visually assessed meat colour [47]. The lack of effect of SCD SNP on meat colour observed in this study may be due to the short period between slaughter and carcass assessments.

4.3. FASN Gene Polymorphisms

The FASN gene encodes an essential homodimeric cytosolic enzyme critical for de novo lipogenesis by catalysing palmitic acid synthesis from acetyl-CoA and malonyl-CoA [48,49]. Palmitic acid is a predominant fatty acid in beef and is used as a substrate for the synthesis of other fatty acids through elongation and desaturation [5,50]. Fatty acid composition is reported to influence carcass traits [51,52]. The majority of association studies on FASN gene polymorphisms have focused on fatty acid composition [53,54,55,56] with little emphasis on carcass traits [3,40]. Polymorphism in the FASN thioesterase domain was found to have a significant effect on beef grade by influencing fat deposition in Korean cattle [48]. Oh et al. [57] examined the effect of five missense SNP in the FASN gene and reported an association of all the SNP with carcass backfat thickness and marbling score with no influence on carcass weight in Korean cattle. In another study, Raza et al. [3] reported associations between g.13192T>C with backfat thickness and EMA, and g.13232C>T with IMF in Qinchuan cattle. Polymorphism in FASN influenced HCW, but no associations with subcutaneous fat thickness and marbling score were observed in crossbred beef cattle [22], or with marbling score in the loin eye muscle of purebred American Angus bulls [58]. Matsuhashi et al. [56] reported no association between FASN g.16024A>G SNP with carcass weight and subcutaneous fat thickness of Japanese black cattle. In this study, polymorphisms in the g.50784533C>G and g.50792445C>T loci were correlated with subcutaneous fat depth, while g.50786221A>G and g.50787362C>T were correlated with marbling scores. The loci g.50787362C>T and g.50792445C>T are located in the intron, hence they may influence carcass fat content through regulation of alternative splicing or gene expression [59]. Furthermore, the allelic mutations in g.50784533C>G locus in exon 13 resulted in amino acid substitution from alanine to glycine at PCS position 2066 that encodes for malonyl-CoA-/acetyl-CoA-acyl carrier protein-transacylase (AT/MT) domain [48]. The AT/MT domain catalyses the transfer of acetyl-CoA to the acyl carrier protein and transacylates the malonyl-CoA to the acyl carrier protein, which is the initial step in de novo fatty acids synthesis [60]; thus, the SNP may influence fatty acids synthesis. The differences between studies may be due to breed effect or in the SNP studied [61,62]. Effects of genetic markers on phenotypic outcomes are often breed-specific, and may not be extrapolated to all cattle breeds [18,25,43]. Since marbling and MSA index scores were highly correlated, these findings indicate that SNP in the FASN gene may be used in marker-assisted selection for improved carcass grades and beef eating quality in northern Australian tropical crossbred cattle.

5. Conclusions

This study aimed to examine SNP present in the exons and introns of FABP4, SCD, and FASN genes of northern Australian tropical crossbred beef cattle and examine their associations with the chiller-assessed carcass traits. The results showed significant correlations between SNP in FABP4 with hump height and loin eye muscle area, SCD SNP with carcass weight, loin eye muscle area, and marbling, while FASN SNP were correlated with subcutaneous fat thickness and marbling score. Therefore, the hypothesis that SNP in the FABP4, SCD, and FASN genes are associated with chiller-assessed carcass traits of tropically adapted northern Australian crossbred beef cattle was accepted. These findings indicate the involvement of these genes in carcass traits of the tropical crossbred beef cattle population of northern Australia and the potential to use SNP in carcass grade and meat quality improvement through marker-assisted selection.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/agriculture12081171/s1, Figure S1: Gel image of the amplification products of the three target genes visualised in 0.8% agarose gel; Table S1: SNP genetic variants identified in the FABP4, SCD, and FASN genes; Figure S2: Correlation coefficients between all SNP in FABP4 gene. The rectangle represents distance-based clustering of SNP loci; Figure S3: Correlation coefficients between all SNP in SCD gene. The rectangle represents distance-based clustering of SNP loci; Figure S4: Correlation coefficients between all SNP in FASN gene. The rectangle represents distance-based clustering of SNP loci.

Author Contributions

Conceptualisation, A.E.O.M.-A., C.P.G., E.C., B.S.M.-A., R.T.K. and F.W.M.; methodology, A.E.O.M.-A., C.P.G., B.S.M.-A., R.T.K., E.C., O.A.A., S.B.P., J.R.O. and F.W.M.; software, A.E.O.M.-A. and O.A.A.; validation, A.E.O.M.-A., C.P.G., R.T.K., E.C. and B.S.M.-A.; formal analysis, F.W.M. and O.A.A.; investigation, F.W.M.; resources, A.E.O.M.-A., C.P.G., E.C., R.T.K. and B.S.M.-A.; data curation and writing—original draft preparation, F.W.M.; writing—reviewing and editing, A.E.O.M.-A., C.P.G., E.C., R.T.K., O.A.A., S.B.P., J.R.O. and B.S.M.-A.; supervision, A.E.O.M.-A., C.P.G., E.C., R.T.K. and B.S.M.-A.; project administration, A.E.O.M.-A., E.C. and C.P.G.; funding acquisition, A.E.O.M.-A., C.P.G. and E.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Cooperative Research Centre Projects (CRC-P) [grant number CRC P-58599] from the Australian Government’s Department of Industry, Innovation and Science, and a Ph.D. scholarship funded by CRC-P and the College of Public Health, Medical and Veterinary Sciences, James Cook University, Queensland, Australia, awarded to the first author.

Institutional Review Board Statement

The study was conducted in accordance with the Australian code of practice for the care and use of animals for scientific purposes (8th Edition, 2013), and approved by the James Cook University Animal Ethics Committee (Approval Number 2639 issued on 5 July 2019).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available from the corresponding author on request.

Acknowledgments

The authors gratefully acknowledge the College of Public Health, Medical and Veterinary Sciences, James Cook University, Townsville, QLD, the Cooperative Research Centre Projects, Canberra, ACT, the Commonwealth Scientific and Industrial Research Organization, Agrimix Pty Ltd and North Australian Pastoral Company, QLD. We are also grateful to the feedlot and abattoir staff for cattle management, handling, slaughter, and grading. We appreciate the financial support for this research through the Cooperative Research Centre Projects (CRC-P) [grant number CRC P-58599] research grant from the Australian Government’s Department of Industry, Innovation and Science, and a Ph.D. scholarship funded by CRC-P and the College of Public Health, Medical and Veterinary Sciences, James Cook University, Queensland, Australia, awarded to the first-named author.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Testa, M.L.; Grigioni, G.; Panea, B.; Pavan, E. Color and Marbling as Predictors of Meat Quality Perception of Argentinian Consumers. Foods 2021, 10, 1465. [Google Scholar] [CrossRef] [PubMed]
  2. Goszczynski, D.E.; Papaleo-Mazzucco, J.; Ripoli, M.V.; Villarreal, E.L.; Rogberg-Muñoz, A.; Mezzadra, C.A.; Melucci, L.M.; Giovambattista, G. Genetic Variation in FABP4 and Evaluation of Its Effects on Beef Cattle Fat Content. Anim. Biotechnol. 2017, 28, 211–219. [Google Scholar] [CrossRef] [PubMed]
  3. Raza, S.H.A.; Gui, L.; Khan, R.; Schreurs, N.M.; Xiaoyu, W.; Wu, S.; Mei, C.; Wang, L.; Ma, X.; Wei, D.; et al. Association between FASN gene polymorphisms ultrasound carcass traits and intramuscular fat in Qinchuan cattle. Gene 2018, 645, 55–59. [Google Scholar] [CrossRef] [PubMed]
  4. Coleman, L.W.; Hickson, R.E.; Schreurs, N.M.; Martin, N.P.; Kenyon, P.R.; Lopez-Villalobos, N.; Morris, S.T. Carcass characteristics and meat quality of Hereford sired steers born to beef-cross-dairy and Angus breeding cows. Meat Sci. 2016, 121, 403–408. [Google Scholar] [CrossRef]
  5. Scollan, N.; Hocquette, J.-F.; Nuernberg, K.; Dannenberger, D.; Richardson, I.; Moloney, A. Innovations in beef production systems that enhance the nutritional and health value of beef lipids and their relationship with meat quality. Meat Sci. 2006, 74, 17–33. [Google Scholar] [CrossRef]
  6. Davis, G. Genetic parameters for tropical beef cattle in northern Australia: A review. Aust. J. Agric. Res. 1993, 44, 179–198. [Google Scholar] [CrossRef]
  7. Greenwood, P.L.; Gardner, G.E.; Ferguson, D.M. Current situation and future prospects for the Australian beef industry—A review. Asian-Australas. J. Anim. Sci. 2018, 31, 992–1006. [Google Scholar] [CrossRef]
  8. Polkinghorne, R.; Thompson, J.M.; Watson, R.; Gee, A.; Porter, M. Evolution of the Meat Standards Australia (MSA) beef grading system. Aust. J. Exp. Agric. 2008, 48, 1351–1359. [Google Scholar] [CrossRef]
  9. Tait, R.G.; Shackelford, S.D.; Wheeler, T.L.; King, D.A.; Keele, J.W.; Casas, E.; Smith, T.P.L.; Bennett, G.L. CAPN1, CAST, and DGAT1 genetic effects on preweaning performance, carcass quality traits, and residual variance of tenderness in a beef cattle population selected for haplotype and allele equalization. J. Anim. Sci. 2014, 92, 5382–5393. [Google Scholar] [CrossRef]
  10. Ardicli, S.; Samli, H.; Alpay, F.; Dincel, D.; Soyudal, B.; Balci, F. Association of Single Nucleotide Polymorphisms in the FABP4 Gene with Carcass Characteristics and Meat Quality in Holstein Bulls. Ann. Anim. Sci. 2017, 17, 117–130. [Google Scholar] [CrossRef]
  11. Li, X.; Ekerljung, M.; Lundström, K.; Lundén, A. Association of polymorphisms at DGAT1, leptin, SCD1, CAPN1 and CAST genes with color, marbling and water holding capacity in meat from beef cattle populations in Sweden. Meat Sci. 2013, 94, 153–158. [Google Scholar] [CrossRef] [PubMed]
  12. Gill, J.L.; Bishop, S.C.; McCorquodale, C.; Williams, J.L.; Wiener, P. Association of selected SNP with carcass and taste panel assessed meat quality traits in a commercial population of Aberdeen Angus-sired beef cattle. Genet. Sel. Evol. 2009, 41, 36. [Google Scholar] [CrossRef] [PubMed]
  13. Curi, R.A.; Palmieri, D.A.; Suguisawa, L.; Ferraz, A.L.J.; de Oliveira, H.N.; Furlan, L.R.; Silveira, A.C.; Lopes, C.R. Effects of GHR gene polymorphisms on growth and carcass traits in Zebu and crossbred beef cattle. Livest. Sci. 2006, 101, 94–100. [Google Scholar] [CrossRef]
  14. Michal, J.J.; Zhang, Z.W.; Gaskins, C.T.; Jiang, Z. The bovine fatty acid binding protein 4 gene is significantly associated with marbling and subcutaneous fat depth in Wagyu x Limousin F2 crosses. Anim. Genet. 2006, 37, 400–402. [Google Scholar] [CrossRef] [PubMed]
  15. Mannen, H. Identification and utilization of genes associated with beef qualities. Anim. Sci. J. 2011, 82, 1–7. [Google Scholar] [CrossRef]
  16. Barendse, W.; Bunch, R.J.; Thomas, M.B.; Harrison, B.E. A splice site single nucleotide polymorphism of the fatty acid binding protein 4 gene appears to be associated with intramuscular fat deposition in longissimus muscle in Australian cattle. Anim. Genet. 2009, 40, 770–773. [Google Scholar] [CrossRef] [PubMed]
  17. Lee, S.H.; Van Der Werf, J.H.J.; Lee, S.H.; Park, E.W.; Oh, S.J.; Gibson, J.P.; Thompson, J.M. Genetic polymorphisms of the bovineFatty acid binding protein 4gene are significantly associated with marbling and carcass weight in Hanwoo (Korean Cattle). Anim. Genet. 2010, 41, 442–444. [Google Scholar] [CrossRef]
  18. Shin, S.-C.; Heo, J.-P.; Chung, E.-R. Genetic variants of the FABP4 gene are associated with marbling scores and meat quality grades in Hanwoo (Korean cattle). Mol. Biol. Rep. 2012, 39, 5323–5330. [Google Scholar] [CrossRef]
  19. Yin, B.-Z.; Fang, J.-C.; Zhang, J.-S.; Zhang, L.-M.; Xu, C.; Xu, H.-Y.; Shao, J.; Xia, G.-J. Correlations between single nucleotide polymorphisms in FABP4 and meat quality and lipid metabolism gene expression in Yanbian yellow cattle. PLoS ONE 2020, 15, e0234328. [Google Scholar] [CrossRef]
  20. Kneeland, J.; Li, C.; Basarab, J.; Snelling, W.M.; Benkel, B.; Murdoch, B.; Hansen, C.; Moore, S.S. Identification and fine mapping of quantitative trait loci for growth traits on bovine chromosomes 2, 6, 14, 19, 21, and 23 within one commercial line of Bos taurus. J. Anim. Sci. 2004, 82, 3405–3414. [Google Scholar] [CrossRef]
  21. Hillgartner, F.B.; Salati, L.M.; Goodridge, A.G. Physiological and molecular mechanisms involved in nutritional regulation of fatty acid synthesis. Physiol. Rev. 1995, 75, 47–76. [Google Scholar] [CrossRef] [PubMed]
  22. Rempel, L.A.; Casas, E.; Shackelford, S.D.; Wheeler, T.L. Relationship of polymorphisms within metabolic genes and carcass traits in crossbred beef cattle. J. Anim. Sci. 2012, 90, 1311–1316. [Google Scholar] [CrossRef] [PubMed]
  23. Jiang, Z.; Michal, J.J.; Tobey, D.J.; Daniels, T.F.; Rule, D.C.; MacNeil, M.D. Significant associations of stearoyl-CoA desaturase (SCD1) gene with fat deposition and composition in skeletal muscle. Int. J. Biol. Sci. 2008, 4, 345–351. [Google Scholar] [CrossRef] [PubMed]
  24. Wu, X.X.; Yang, Z.P.; Shi, X.K.; Li, J.Y.; Ji, D.J.; Mao, Y.J.; Chang, L.L.; Gao, H.J. Association of SCD1 and DGAT1 SNP with the intramuscular fat traits in Chinese Simmental cattle and their distribution in eight Chinese cattle breeds. Mol. Biol. Rep. 2012, 39, 1065–1071. [Google Scholar] [CrossRef]
  25. Avilés, C.; Polvillo, O.; Peña, F.; Juárez, M.; Martínez, A.L.; Molina, A. Associations between DGAT1, FABP4, LEP, RORC, and SCD1 gene polymorphisms and fat deposition in Spanish commercial beef. J. Anim. Sci. 2013, 91, 4571–4577. [Google Scholar] [CrossRef]
  26. Hocquette, J.-F.; Botreau, R.; Picard, B.; Jacquet, A.; Pethick, D.W.; Scollan, N.D. Opportunities for predicting and manipulating beef quality. Meat Sci. 2012, 92, 197–209. [Google Scholar] [CrossRef]
  27. National Health and Medical Research Council. Australian Code of Practice for the Care and Use of Animals for Scientific Purposes, 8th ed.; National Health and Medical Research Council: Canberra, ACT, Australia, 2013; ISBN 186-496-5-975. [Google Scholar]
  28. Mwangi, F.W.; Gardiner, C.P.; Walker, G.; Hall, T.J.; Malau-Aduli, B.S.; Kinobe, R.T.; Malau-Aduli, A.E.O. Growth Performance and Plasma Metabolites of Grazing Beef Cattle Backgrounded on Buffel or Buffel-Desmanthus Mixed Pastures. Animals 2021, 11, 2355. [Google Scholar] [CrossRef]
  29. Meat & Livestock Australia Meat Standards Australia Standards Manual. Section 5: Livestock Supply; Meat & Livestock Australia: Sydney, Australia, 2020. [Google Scholar]
  30. Meat Standards Australia Standards Manual. Section 7: Processor; Meat & Livestock Australia Limited: Sydney, Australia, 2020. [Google Scholar]
  31. Pewan, S.B.; Otto, J.R.; Huerlimann, R.; Budd, A.M.; Mwangi, F.W.; Edmunds, R.C.; Holman, B.W.B.; Henry, M.L.E.; Kinobe, R.T.; Adegboye, O.A.; et al. Next Generation Sequencing of Single Nucleotide Polymorphic DNA-Markers in Selecting for Intramuscular Fat, Fat Melting Point, Omega-3 Long-Chain Polyunsaturated Fatty Acids and Meat Eating Quality in Tattykeel Australian White MARGRA Lamb. Foods 2021, 10, 2288. [Google Scholar] [CrossRef]
  32. Nei, M.; Roychoudhury, A.K. Sampling variances of heterozygosity and genetic distance. Genetics 1974, 76, 379–390. [Google Scholar] [CrossRef]
  33. Lin, C.-Y.; Xing, G.; Xing, C. Measuring linkage disequilibrium by the partial correlation coefficient. Heredity 2012, 109, 401–402. [Google Scholar] [CrossRef]
  34. Kelly, M.J.; Tume, R.K.; Fortes, M.; Thompson, J.M. Whole-genome association study of fatty acid composition in a diverse range of beef cattle breeds. J. Anim. Sci. 2014, 92, 1895–1901. [Google Scholar] [CrossRef] [PubMed]
  35. Serrote, C.M.L.; Reiniger, L.R.S.; Silva, K.B.; Rabaiolli, S.M.D.S.; Stefanel, C.M. Determining the Polymorphism Information Content of a molecular marker. Gene 2020, 726, 144175. [Google Scholar] [CrossRef] [PubMed]
  36. Botstein, D.; White, R.L.; Skolnick, M.; Davis, R.W. Construction of a genetic linkage map in man using restriction fragment length polymorphisms. Am. J. Hum. Genet. 1980, 32, 314–331. [Google Scholar] [PubMed]
  37. Jurie, C.; Cassar-Malek, I.; Bonnet, M.; Leroux, C.; Bauchart, D.; Boulesteix, P.; Pethick, D.W.; Hocquette, J.-F. Adipocyte fatty acid-binding protein and mitochondrial enzyme activities in muscles as relevant indicators of marbling in cattle. J. Anim. Sci. 2007, 85, 2660–2669. [Google Scholar] [CrossRef] [PubMed]
  38. Hoashi, S.; Hinenoya, T.; Tanaka, A.; Ohsaki, H.; Sasazaki, S.; Taniguchi, M.; Oyama, K.; Mukai, F.; Mannen, H. Association between fatty acid compositions and genotypes of FABP4 and LXR-alpha in Japanese Black cattle. BMC Genet. 2008, 9, 84. [Google Scholar] [CrossRef]
  39. Damon, M.; Louveau, I.; Lefaucheur, L.; Lebret, B.; Vincent, A.; Leroy, P.; Sanchez, M.P.; Herpin, P.; Gondret, F. Number of intramuscular adipocytes and fatty acid binding protein-4 content are significant indicators of intramuscular fat level in crossbred Large White × Duroc pigs. J. Anim. Sci. 2006, 84, 1083–1092. [Google Scholar] [CrossRef]
  40. Oh, D.-Y.; Lee, Y.-S.; La, B.-M.; Yeo, J.-S. Identification of the SNP (Single Nucleotide Polymorphism) for Fatty Acid Composition Associated with Beef Flavor-related FABP4 (Fatty Acid Binding Protein 4) in Korean Cattle. Asian-Australas. J. Anim. Sci. 2012, 25, 913–920. [Google Scholar] [CrossRef]
  41. Cho, S.-A.; Park, T.S.; Yoon, D.-H.; Cheong, H.S.; Namgoong, S.; Park, B.L.; Lee, H.W.; Han, C.S.; Kim, E.M.; Cheong, I.-C.; et al. Identification of genetic polymorphisms in FABP3 and FABP4 and putative association with back fat thickness in Korean native cattle. BMB Rep. 2008, 41, 29–34. [Google Scholar] [CrossRef]
  42. Pringle, T.D.; Williams, S.E.; Lamb, B.S.; Johnson, D.D.; West, R.L. Carcass characteristics, the calpain proteinase system, and aged tenderness of Angus and Brahman crossbred steers. J. Anim. Sci. 1997, 75, 2955–2961. [Google Scholar] [CrossRef]
  43. Bartoň, L.; Bureš, D.; Kott, T.; Řehák, D. Associations of polymorphisms in bovine DGAT1, FABP4, FASN, and PPARGC1A genes with intramuscular fat content and the fatty acid composition of muscle and subcutaneous fat in Fleckvieh bulls. Meat Sci. 2016, 114, 18–23. [Google Scholar] [CrossRef]
  44. Kelly, M.J.; Tume, R.K.; Newman, S.; Thompson, J.M. Genetic variation in fatty acid composition of subcutaneous fat in cattle. Anim. Prod. Sci. 2013, 53, 129–133. [Google Scholar] [CrossRef]
  45. Bernard, A.; Leroux, C.; Chilliard, Y. Expression and nutritional regulation of Stearoyl-CoA desaturase genes in the ruminant mammary gland: Relationship with milk fatty acid composition. In Stearoyl-CoA Desaturase Genes in Lipid Metabolism; Ntambi, J.M., Ed.; Springer Science+Business Media: Clermont-Ferrand, France, 2013; ISBN 9781461479697. [Google Scholar]
  46. Ohsaki, H.; Tanaka, A.; Hoashi, S.; Sasazaki, S.; Oyama, K.; Taniguchi, M.; Mukai, F.; Mannen, H. Effect of SCD and SREBP genotypes on fatty acid composition in adipose tissue of Japanese Black cattle herds. Anim. Sci. J. 2009, 80, 225–232. [Google Scholar] [CrossRef]
  47. Fiems, L.O.; Campeneere, S.D.; De Smet, S.; Van de Voorde, G.; Vanacker, J.M.; Boucqué, C.V. Relationship between fat depots in carcasses of beef bulls and effect on meat colour and tenderness. Meat Sci. 2000, 56, 41–47. [Google Scholar] [CrossRef]
  48. Oh, D.-Y.; Nam, I.; Hwang, S.; Kong, H.; Lee, H.; Ha, J.; Baik, M.; Oh, M.H.; Kim, S.; Han, K.; et al. In vivo evidence on the functional variation within fatty acid synthase gene associated with lipid metabolism in bovine longissimus dorsi muscle tissue. Genes Genom. 2018, 40, 289–294. [Google Scholar] [CrossRef] [PubMed]
  49. Roy, R.; Taourit, S.; Zaragoza, P.; Eggen, A.; Rodellar, C. Genomic structure and alternative transcript of bovine fatty acid synthase gene (FASN): Comparative analysis of the FASN gene between monogastric and ruminant species. Cytogenet. Genome Res. 2005, 111, 65–73. [Google Scholar] [CrossRef] [PubMed]
  50. Sampath, H.; Ntambi, J.M. The fate and intermediary metabolism of stearic acid. Lipids 2005, 40, 1187–1191. [Google Scholar] [CrossRef]
  51. Yeon, S.H.; Lee, S.H.; Choi, B.H.; Lee, H.J.; Jang, G.W.; Lee, K.T.; Kim, K.H.; Lee, J.H.; Chung, H.Y. Genetic variation of FASN is associated with fatty acid composition of Hanwoo. Meat Sci. 2013, 94, 133–138. [Google Scholar] [CrossRef]
  52. Hoehne, A.; Nuernberg, G.; Kuehn, C.; Nuernberg, K. Relationships between intramuscular fat content, selected carcass traits, and fatty acid profile in bulls using a F 2-population. Meat Sci. 2012, 90, 629–635. [Google Scholar] [CrossRef]
  53. Lee, J.; Jin, M.; Lee, Y.; Ha, J.; Yeo, J.; Oh, D. Gene–gene interactions of fatty acid synthase (FASN) using multifactor-dimensionality reduction method in Korean cattle. Mol. Biol. Rep. 2014, 41, 2021–2027. [Google Scholar] [CrossRef]
  54. Abe, T.; Saburi, J.; Hasebe, H.; Nakagawa, T.; Misumi, S.; Nade, T.; Nakajima, H.; Shoji, N.; Kobayashi, M.; Kobayashi, E. Novel Mutations of the FASN Gene and Their Effect on Fatty Acid Composition in Japanese Black Beef. Biochem. Genet. 2009, 47, 397–411. [Google Scholar] [CrossRef]
  55. Li, C.; Aldai, N.; Vinsky, M.; Dugan, M.E.R.; McAllister, T.A. Association analyses of single nucleotide polymorphisms in bovine stearoyl-CoA desaturase and fatty acid synthase genes with fatty acid composition in commercial cross-bred beef steers. Anim. Genet. 2012, 43, 93–97. [Google Scholar] [CrossRef] [PubMed]
  56. Matsuhashi, T.; Maruyama, S.; Uemoto, Y.; Kobayashi, N.; Mannen, H.; Abe, T.; Sakaguchi, S.; Kobayashi, E. Effects of bovine fatty acid synthase, stearoyl-coenzyme A desaturase, sterol regulatory element-binding protein 1, and growth hormone gene polymorphisms on fatty acid composition and carcass traits in Japanese Black cattle. J. Anim. Sci. 2011, 89, 12–22. [Google Scholar] [CrossRef] [PubMed]
  57. Oh, D.; Lee, Y.; La, B.; Yeo, J.; Chung, E.; Kim, Y.; Lee, C. Fatty acid composition of beef is associated with exonic nucleotide variants of the gene encoding FASN. Mol. Biol. Rep. 2012, 39, 4083–4090. [Google Scholar] [CrossRef] [PubMed]
  58. Zhang, S.; Knight, T.J.; Reecy, J.M.; Beitz, D.C. DNA polymorphisms in bovine fatty acid synthase are associated with beef fatty acid composition. Anim. Genet. 2008, 39, 62–70. [Google Scholar] [CrossRef] [PubMed]
  59. Jo, B.-S.; Choi, S.S. Introns: The Functional Benefits of Introns in Genomes. Genom. Inform. 2015, 13, 112–118. [Google Scholar] [CrossRef] [PubMed]
  60. Maier, T.; Jenni, S.; Ban, N. Architecture of Mammalian Fatty Acid Synthase at 4.5 Å Resolution. Science 2006, 311, 1258–1262. [Google Scholar] [CrossRef]
  61. Li, B.; VanRaden, P.M.; Null, D.J.; O’Connell, J.R.; Cole, J.B. Major quantitative trait loci influencing milk production and conformation traits in Guernsey dairy cattle detected on Bos taurus autosome 19. J. Dairy Sci. 2021, 104, 550–560. [Google Scholar] [CrossRef]
  62. Maharani, D.; Jo, C.-R.; Jeon, J.-T.; Lee, J.-H. Quantitative Trait Loci and Candidate Genes Affecting Fatty Acid Composition in Cattle and Pig. Korean J. Food Sci. Anim. Resour. 2011, 31, 325–338. [Google Scholar] [CrossRef]
Figure 1. Single nucleotide polymorphisms on the FABP4 gene. (A) Clustering map of genetic variants; Agriculture 12 01171 i001 homozygotes similar to the reference sequence (Hereford), Agriculture 12 01171 i002 heterozygotes, and Agriculture 12 01171 i003 alternative homozygotes. (B) Correlations between detected SNP and carcass traits. * p < 0.05, ** p < 0.01, and *** p < 0.001.
Figure 1. Single nucleotide polymorphisms on the FABP4 gene. (A) Clustering map of genetic variants; Agriculture 12 01171 i001 homozygotes similar to the reference sequence (Hereford), Agriculture 12 01171 i002 heterozygotes, and Agriculture 12 01171 i003 alternative homozygotes. (B) Correlations between detected SNP and carcass traits. * p < 0.05, ** p < 0.01, and *** p < 0.001.
Agriculture 12 01171 g001
Figure 2. (A) Clustering map of genetic variants of the SCD SNP; Agriculture 12 01171 i001 homozygotes similar to the reference gene (Hereford), Agriculture 12 01171 i002 heterozygotes, and Agriculture 12 01171 i003 the alternative homozygotes. (B) Correlations between SNP and carcass traits. * p < 0.05, ** p < 0.01, and *** p < 0.001.
Figure 2. (A) Clustering map of genetic variants of the SCD SNP; Agriculture 12 01171 i001 homozygotes similar to the reference gene (Hereford), Agriculture 12 01171 i002 heterozygotes, and Agriculture 12 01171 i003 the alternative homozygotes. (B) Correlations between SNP and carcass traits. * p < 0.05, ** p < 0.01, and *** p < 0.001.
Agriculture 12 01171 g002
Figure 3. (A) Clustering map of genetic variants of the FASN SNP; Agriculture 12 01171 i001 homozygotes similar to the reference gene (Hereford), Agriculture 12 01171 i002 heterozygotes, and Agriculture 12 01171 i003 the alternative homozygotes. (B) Correlations between genes SNP and carcass traits. * p < 0.05, ** p < 0.01, and *** p < 0.001.
Figure 3. (A) Clustering map of genetic variants of the FASN SNP; Agriculture 12 01171 i001 homozygotes similar to the reference gene (Hereford), Agriculture 12 01171 i002 heterozygotes, and Agriculture 12 01171 i003 the alternative homozygotes. (B) Correlations between genes SNP and carcass traits. * p < 0.05, ** p < 0.01, and *** p < 0.001.
Agriculture 12 01171 g003
Figure 4. Multiple comparisons of hump height and loin eye muscle area (EMA) for FABP4 g.44677205A>G genotypic variants AA ( Agriculture 12 01171 i004), AG ( Agriculture 12 01171 i005), and GG ( Agriculture 12 01171 i006).
Figure 4. Multiple comparisons of hump height and loin eye muscle area (EMA) for FABP4 g.44677205A>G genotypic variants AA ( Agriculture 12 01171 i004), AG ( Agriculture 12 01171 i005), and GG ( Agriculture 12 01171 i006).
Agriculture 12 01171 g004
Figure 5. Multiple comparisons of P8 fat, marbling, backfat, and fat class of SCD g.21275851C>A genotypic variants CC ( Agriculture 12 01171 i004), CA ( Agriculture 12 01171 i005), and AA ( Agriculture 12 01171 i006).
Figure 5. Multiple comparisons of P8 fat, marbling, backfat, and fat class of SCD g.21275851C>A genotypic variants CC ( Agriculture 12 01171 i004), CA ( Agriculture 12 01171 i005), and AA ( Agriculture 12 01171 i006).
Agriculture 12 01171 g005
Table 1. Primer sequences for target gene amplification.
Table 1. Primer sequences for target gene amplification.
Gene 1FragmentPrimerSequence (5′ to 3′)Annealing Temp (°C)Product Size (bp)
SCD1ForwardGGAAGAAGACATCCGCCCTGAAAT605092
ReverseAGGAAGCGAGATTGGCACTGTATG60
SCD2ForwardGGAAGAAGACATCCGCCCTGAAAT6010,177
ReverseTGCCTCTGAGGGGATCTATTTGGT60
SCD3ForwardATGAGCCACACTGTGAACAAACCT602861
ReverseTTCTTTTTCTGGACAGGCAAGCCT60
FASN1ForwardTTGAGCTTCTGAGTATGATGGGAG687302
ReverseACCATCTATTATGCCTCCCTCAAC68
FASN2ForwardCTATAAGATCGGTGAGTCCTTGCA688648
ReverseTAGTATTATTCACAGCTCCCTGGC68
FABP4-ForwardGCTAAGACTGCCTGTATGTTCCCC603041
ReverseACCTAGAGAAATAGACAATCGCCC60
1 SCD, stearoyl-CoA desaturase; FASN, fatty acid synthase; FABP4, fatty acid binding protein 4.
Table 2. FABP4, SCD, and FASN gene polymorphisms, protein-coding sequence positions, and non-synonymous amino acid substitutions.
Table 2. FABP4, SCD, and FASN gene polymorphisms, protein-coding sequence positions, and non-synonymous amino acid substitutions.
Gene 1SNP (dbSNP ID) 2PCS Position 3Amino Acid Substitution
FABP4g.44677959T>C (rs110757796)220Isoleucine to Valine
SCDg.21272422C>T (rs41255693)878Alanine to Valine
FASNg.50784533C>G (rs481622676)2066Alanine to Glycine
g.50786496A>G3145Serine to Glycine
g.50788575T>C (rs41919993)4168Tyrosine to Histidine
g.50790973C>A (rs109149276)5572Leucine to Isoleucine
g.50794099T>C7277Isoleucine to Threonine
1 FABP4: fatty acid binding protein 4, SCD: stearoyl-CoA desaturase, FASN: fatty acid synthase. 2 SNP: single nucleotide polymorphism. Variant dbSNP ID are based on the Bovine Genome Variation Database (BGVD). SNP without dbSNP ID are not listed in BGVD. 3 PCS: protein-coding sequence.
Table 3. Least square means ± SD of carcass traits due to SNP in FABP4 g.44677205A>G, SCD g.21275851C>A, and FASN g.50784533C>G in northern Australian tropical crossbred beef cattle.
Table 3. Least square means ± SD of carcass traits due to SNP in FABP4 g.44677205A>G, SCD g.21275851C>A, and FASN g.50784533C>G in northern Australian tropical crossbred beef cattle.
Variable 1 p-Value 2
FABP4
g.44677205A>GTotal (n = 44)AA (n = 22)AG (n = 16)GG (n = 6)
HCW (kg)329.9 ± 26.15331.1 ± 20.06333.0 ± 32.81317.6 ± 27.260.16
P8 fat (mm)16.7 ± 5.8017.0 ± 5.9617.0 ± 6.4114.6 ± 3.270.66
Hump height (mm)112.5 ± 15.42109.3 ± 13.39111.8 ± 17.5125.8 ± 10.68<0.01
EMA (cm2)86.8 ± 7.2688.4 ± 8.7186.7 ± 4.9381.3 ± 3.880.05
Marbling (score)361.5 ± 79.77372.7 ± 83.05353.1 ± 81.79343.3 ± 67.130.86
Meat colour (score)2.3 ± 1.112.1 ± 1.182.3 ± 0.962.5 ± 1.380.70
Back fat (mm)14.7 ± 5.8015.0 ± 5.9615.0 ± 6.4112.6 ± 3.270.74
Fat class (score)3.8 ± 0.653.9 ± 0.753.8 ± 0.623.8 ± 0.410.98
MSA index51.11 ± 2.3251.3 ± 2.2851.2 ± 2.649.6 ± 0.420.39
SCD
g.21275851C>ATotal (n = 44)CC (n = 9)CA (n = 15)AA (n = 20)
HCW (kg)329.9 ± 26.15333.0 ± 38.21333.2 ± 23.32326.1 ± 22.440.75
P8 fat (mm)16.7 ± 5.8019.6 ± 6.7318.6 ± 5.0314.0 ± 4.880.01
Hump height (mm)112.5 ± 15.42105.5 ± 12.61116.6 ± 15.99112.5 ± 15.680.24
EMA (cm2)86.8 ± 7.2685.8 ± 5.4988.2 ± 8.6086.2 ± 7.070.82
Marbling (score)361.5 ± 79.77421.1 ± 75.9338.0 ± 83.85352.5 ± 67.740.02
Meat colour (score)2.3 ± 1.112.2 ± 1.202.3 ± 1.232.3 ± 1.030.98
Back fat (mm)14.7 ± 5.8017.5 ± 6.7316.6 ± 5.0512.0 ± 4.880.01
Fat class (score)3.8 ± 0.6554.2 ± 0.6674.0 ± 0.5943.6 ± 0.5980.02
MSA index51.1 ± 2.3252.3 ± 2.1850.9 ± 2.8250.7 ± 1.900.19
FASN
g.50784533C>GTotal (n = 44)CC (n = 30)CG (n = 12)GG (n = 2)
HCW (kg)329.9 ± 26.15327.2 ± 23.24329.5 ± 24.73377.5 ± 40.310.62
P8 fat (mm)16.7 ± 5.8015.5 ± 5.1018.4 ± 5.8425.0 ± 9.900.10
Hump height (mm)112.5 ± 15.42110.2 ± 14.97110.2 ± 15.60130.2 ± 14.140.92
EMA (cm2)86.8 ± 7.2686.5 ± 7.9087.4 ± 6.0588.5 ± 6.360.69
Marbling (score)361.5 ± 79.77380.2 ± 91.95350.2 ± 62.68345.2 ± 49.500.23
Meat colour (score)2.3 ± 1.112.5 ± 1.121.9 ± 1.152.0 ± 0.7100.31
Back fat (mm)14.7 ± 5.8013.5 ± 5.1016.4 ± 5.8423.0 ± 9.900.10
Fat class (score)3.8 ± 0.6503.8 ± 0.6483.9 ± 0.6694.5 ± 0.7070.38
MSA index51.1 ± 2.3250.9 ± 2.0151.4 ± 3.2451.8 ± 1.220.75
1 SCD, stearoyl-CoA desaturase; FASN, fatty acid synthase; FABP4, fatty acid binding protein 4; HCW, hot carcass weight; P8 fat; subcutaneous fat thickness at the rump site; EMA, loin eye muscle area; MSA index, Meat Standards Australia index. 2 ANOVA p-value.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mwangi, F.W.; Pewan, S.B.; Otto, J.R.; Adegboye, O.A.; Charmley, E.; Gardiner, C.P.; Malau-Aduli, B.S.; Kinobe, R.T.; Malau-Aduli, A.E.O. Single Nucleotide Polymorphisms in the Fatty Acid Binding Protein 4, Fatty Acid Synthase and Stearoyl-CoA Desaturase Genes Influence Carcass Characteristics of Tropical Crossbred Beef Steers. Agriculture 2022, 12, 1171. https://doi.org/10.3390/agriculture12081171

AMA Style

Mwangi FW, Pewan SB, Otto JR, Adegboye OA, Charmley E, Gardiner CP, Malau-Aduli BS, Kinobe RT, Malau-Aduli AEO. Single Nucleotide Polymorphisms in the Fatty Acid Binding Protein 4, Fatty Acid Synthase and Stearoyl-CoA Desaturase Genes Influence Carcass Characteristics of Tropical Crossbred Beef Steers. Agriculture. 2022; 12(8):1171. https://doi.org/10.3390/agriculture12081171

Chicago/Turabian Style

Mwangi, Felista W., Shedrach B. Pewan, John R. Otto, Oyelola A. Adegboye, Edward Charmley, Christopher P. Gardiner, Bunmi S. Malau-Aduli, Robert T. Kinobe, and Aduli E. O. Malau-Aduli. 2022. "Single Nucleotide Polymorphisms in the Fatty Acid Binding Protein 4, Fatty Acid Synthase and Stearoyl-CoA Desaturase Genes Influence Carcass Characteristics of Tropical Crossbred Beef Steers" Agriculture 12, no. 8: 1171. https://doi.org/10.3390/agriculture12081171

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop