Lipogenic Gene Single Nucleotide Polymorphic DNA Markers Associated with Intramuscular Fat, Fat Melting Point, and Health-Beneficial Omega-3 Long-Chain Polyunsaturated Fatty Acids in Australian Pasture-Based Bowen Genetics Forest Pastoral Angus, Hereford, and Wagyu Beef Cattle

This study used targeted sequencing aimed at identifying single nucleotide polymorphisms (SNP) in lipogenic genes and their associations with health-beneficial omega-3 long-chain polyunsaturated fatty acids (n-3 LC-PUFA), intramuscular fat (IMF), and fat melting point (FMP) of the M. longissimus dorsi muscle in Australian pasture-based Bowen Genetics Forest Pastoral Angus, Hereford, and Wagyu cattle. It was hypothesized that SNP encoding for the fatty acid-binding protein 4 (FABP4), stearoyl-CoA desaturase (SCD), and fatty acid synthase (FASN) genes will be significantly associated with health-beneficial n-3 LC-PUFA and the meat eating quality traits of IMF and FMP in an Australian pasture-based beef production system. Two SNP mutations, g.21267406 T>C and g.21271264 C>A, in the SCD gene were significantly (p < 0.05) associated with IMF, FMP, oleic acid (18:1n-9), linoleic acid (LA) 18:2n-6, alpha-linolenic acid (ALA) 18:3n-3, eicosapentaenoic acid (EPA) 20:5n-3, docosahexaenoic acid (DHA) 22:6-n-3, and docosapentaenoic acid (DPA) 22:5n-3. Significant positive correlations (p < 0.05) between FASN SNP g. 50787138 A>G and FMP, 18:1n-9, ALA, EPA, DHA, DPA, and total n-3 LC-PUFA were also detected. An SNP (g.44678794 G>A) in the FABP4 gene was associated with FMP. These results provide significant insights into the contributions of lipogenic genes to intramuscular fat deposition and the biosynthesis of health-beneficial n-3 LC-PUFA. The findings also unravel the potential use of lipogenic gene polymorphisms in marker-assisted selection to improve the content of health-promoting n-3 LC-PUFA and meat eating quality traits in Australian pasture-based Bowen Genetics Forest Pastoral Angus, Hereford, and Wagyu beef cattle.


Introduction
Consumers demand healthier, safer, and more consistent beef of high eating quality [1,2]. Currently, the Australian and global beef industry, supported by research bodies, has shifted its attention to finding ways of producing meat of consistent eating quality and health-promoting ≥C20 omega-3 long-chain polyunsaturated fatty acids (n-3 LC-PUFA) that play essential roles in human disease prevention [3,4]. Intramuscular fat content (IMF), fat melting point (FMP), and ≥C20 n-3 LC-PUFA composition are associated with consumer satisfaction and key eating quality attributes in the Australian, American, and Japanese beef

Animal Ethics
This project was approved by the James Cook University Animal Ethics Committee (Approval Number A2724), under the Animal Care and Protection Act for the Australian Code for the Care and Use of Animals for Scientific Purposes.

Animals and Management
The animals used in this study were Angus, Hereford, and Wagyu raised and maintained under the same on-farm management at the Bowen Genetics Forest Pastoral stud, Barraba, New South Wales, Australia. A total of 127 yearling Angus, Hereford and Wagyu cattle were utilized for this study. The animals were fed on ryegrass pasture, kept in the same herd, and maintained under the same management routine. This cohort of beef cattle was of similar average age (10-11 months), body condition score (3.0-3.5), and liveweight (450-455 kg) to minimize any phenotypic variation and potential confounding factors.

Muscle Biopsy Sampling Procedure
Based on the procedure first described by Malau-Aduli et al. [48], the animals were gently restrained in a crush, and the hair around the 12th and 13th ribs was shaved with clippers. An alcohol/chlorhexidine disinfectant solution was used to clean the shaved surgical site. Approximately 20 mL of a local anesthetic, lignocaine, was injected at the surgical site before 2 g of muscle biopsy was sampled. The skin was closed with non-absorbable braided suture material to close the incision firmly, and antibiotics were administered to prevent secondary infection. The muscle biopsy samples were immediately placed in a plastic bag on dry ice, flushed with nitrogen gas, transferred into a mobile refrigerator, and stored at −20 • C awaiting further analysis in the laboratory. Each animal was checked daily to ensure that there were no complications. After ten days, the sutures were removed, and the animals were placed in larger yards for further observation for 2-3 days before returning to the herd. The muscle biopsies were analyzed for IMF, FMP, and fatty acid composition.

Determination of Fatty Acid Profile, Intramuscular Fat, and Fat Melting Point
Procedures for the determination of fatty acid profile, IMF, and FMP had previously been described extensively [9,21,49], but for the purposes of clarity, they are repeated herein.

Determination of Intramuscular Fat
Briefly, the muscle sample was homogenized, and 1g transferred to a labelled 50 mL plastic tube containing 20 mL of chloroform: methanol (2:1) solvent and shaken vigorously for 5 min. A filter paper was used to collect the filtrate in another labelled 50 mL tube. Approximately 5 mL of 10% KCl was added to the filtrate to precipitate and separate the inorganic and lipid fractions into two distinct layers. The upper inorganic layer was removed and discarded, while the lower lipid layer was transferred into a clean, dry, pre-weighed and labelled ceramic crucible and evaporated in a laminar fume hood over a heating block. The crucible was cooled and further dried in a desiccator for 10-20 min before it was re-weighed. Samples were analyzed in duplicate to allow for replication and reproducibility. Intramuscular fat percentage was calculated as: [(Final crucible weight) − (Initial crucible weight)/(Initial sample weight)] × 100.

Determination of Fat Melting Point
The crucible containing the extracted IMF was placed in an oven at 100 • C for about 1-2 min to melt the fat. Using air suction, the melted fat was sucked into a thin capillary tube and placed in a refrigerator for about 10 min for the fat to solidify. The fat level in the capillary tube was marked with an indelible pen. The capillary tube was attached to a thermometer and vertically suspended in a beaker containing 80 mL of cold water, gradually heated over a heating block, and closely observed until the fat melted and "slipped" (rose above the mark) within the capillary tube. The temperature at which this Genes 2022, 13, 1411 4 of 21 slip occurred was recorded as the fat melting point. Samples were analyzed in duplicate to allow for replication and reproducibility.

Determination of Fatty Acid Composition
Total lipids in 1 g of un-homogenized muscle tissue samples were extracted overnight. The original phase was a single-phase overnight extraction utilizing CHCl3:MeOH:H2O (1:2:0.8 v/v). The second segment involved phase separation with the addition of CHCl3:saline Milli-Q H2O (1:1 v/v) followed by rotary evaporation of the lower chloroform phase at 40 • C to acquire total lipids. The extracted cumulative lipids were separated into lipid classes by thin-layer chromatography (TLC) using 100 mL of the lipid extract reconstituted in hexane. The extract was marked onto silica gel G plates (200 × 200 × 0.25 mm 3 ) using a micropipette. The TLC plate was developed in an acetone/petroleum ether (1:3.147 v/v) solvent system in a tank comprising a few crystals of butylated hydroxytoluene (BHT) to hinder oxidation. Triacylglycerols, cholesterol and free fatty acids migrated, while phospholipids remained at the origin of the plate. The phospholipids were scraped off the plate into clean screw-capped test tubes for transmethylation and eventual computation of the lipid conversion factor (LCF) of 0.912 based on g fatty acids/g total lipids (0.083 for phospholipids, 0.829 for triacylglycerols and 0% for cholesterol since cholesterol does not have any fatty acids). An aliquot from each total lipid extract was utilized for transmethylation with MeOH:CHCl3:HCl (10:1:1 v/v) for 2 h at 80 • C. Fatty acid methyl esters (FAME) were extracted thrice using hexane:CHCl3 (4:1 v/v). A known concentration of an internal standard (C19:0) was added in a 1500 µL vial encompassing the extracted FAME. The FAME was analyzed on a 7890B gas chromatograph (Agilent Technologies, Palo Alto, CA, USA) furnished with an Equi-tyTM -1 fused 15 m silica capillary column with 0.1 mm internal diameter and 0.1 µm film thickness (Supelco, Bellefonte, PA, USA), a flame ionization sensor, a split/splitless injector and an Agilent Technologies 7683 B Series autosampler. The gas chromatograph settings were splitless mode injection; carrier gas He; original oven temperature 120 • C and then increased to 270 • C at flow rates of 10 • C/min and to 310 • C at 5 • C/min. The Agilent Technologies ChemStation software (Palo Alto, CA, USA) was used to measure fatty acid peaks. The fatty acid identities were established using a Finnigan Thermoquest GCQTM GC/MS fitted with an on-column injector and Thermoquest Xcalibur software (Austin, TX, USA). Fatty acid percentages were calculated as follows: FA % = [(individual fatty acid area) × (100)]/(sum total area of fatty acids). Fatty acid contents were calculated as follows: FA mg/100 g = (Total lipid) × (LCF [0.912]) × ([%FA]/100) × 1000, where 0.912 was the resultant lipid conversion factor.

Blood Collection and Genomic DNA Isolation
Ten mL of blood samples were collected and stored at −80 • C until further analysis. Genomic DNA was extracted from 2 mL of blood using NucleoSpin Blood Kit (Macherey-Nagel GmbH & Co., Duren, Germany) following the procedures prescribed by the manufacturer. The yield and purity of the extracted genomic DNA was quantified using QuantiFluor ® dsDNA System (Promega, Madison, WI, USA).
2.6. Primer Design, Long-Range Amplification of Target Genes, PCR Clean-Up Products, Library Preparation, Sequencing, and Bioinformatics Analysis The procedures for primer design, PCR amplification, library preparation, SNP sequencing, bioinformatics, and data analysis have been described extensively by Pewan et al. [50]. The single coding sequence of the Hereford breed deposited in the Gen-Bank database for FASN, FABP4, and SCD genes with Accession Numbers NC_037346.1, NC_037353.1 and NC_037341.1, respectively, was utilized as a genotyping reference. All primer sequences used for amplifying target genes (Table 1) were designed using Geneious Prime software, synthesized at Integrated DNA Technologies Pte. Ltd., Singapore (itddna.com), and gel images of the resulting PCR products (Figures 1-3) confirmed successful amplification.

Primer Design, Long-Range Amplification of Target Genes, PCR Clean-Up Products, Library Preparation, Sequencing and Bioinformatics Analysis
The procedures for primer design, PCR amplification, library preparation, SNP sequencing, bioinformatics, and data analysis have been described extensively by Pewan et al. [50]. The single coding sequence of the Hereford breed deposited in the GenBank database for FASN, FABP4, and SCD genes with Accession Numbers NC_037346.1, NC_037353.1 and NC_037341.1, respectively, was utilized as a genotyping reference. All primer sequences used for amplifying target genes (Table 1) were designed using Geneious Prime software, synthesized at Integrated DNA Technologies Pte. Ltd., Singapore (itddna.com), and gel images of the resulting PCR products (Figures 1-3) confirmed successful amplification.

Long-Range PCR
Briefly, because of the different fragment lengths and DNA composition, it was necessary to use three different long-range PCR approaches to amplify the FASN, FABP4, and SCD genes. During optimization, all three approaches were tested for all three genes, but only the best performing combinations were utilized.

FASN Gene
TakaRa PrimeSTAR GXL Master Mix (TaKaRa Bio Inc., Kusatsu, Shiga, Japan) was used to perform FASN PCR amplification assay. PCR reaction assay was set up in a total volume of 50 µL containing 10 µL of 5× TakaRa PrimeSTAR GXL Buffer, 200 µM of TaKaRa dNTP Mixture, 1.25 units of TaKaRa PrimeSTAR GXL DNA Polymerase, 0.2 µM of each primer and 100 ng of DNA template. PCR was performed in a SimpliAmp™ Thermal Cycler (Thermofisher Scientific, Melbourne, Australia), in a 2-step protocol using the following conditions: 98 • C initial denaturation for 1 min (1 cycle); 98 • C denaturation for 10 s; 68 • C annealing/extension for 10 min for 30 cycles. PCR success was checked in a 0.8% agarose gel electrophoresis image, as depicted in Figure 1.

FABP4 and SCD Genes
Platinum™ SuperFi™ II PCR Master Mix (Thermofisher Scientific, Melbourne, Australia) and Hot Start II High-Fidelity PCR Master Mix (Thermofisher Scientific, Melbourne, Australia) were used to perform FABP4 and SCD PCR amplification assay under the same PCR conditions described below. The amplification reactions were performed in a total volume of 50 µL containing 25 µL of 2× Platinum™ SuperFi™ II PCR Master Mix or Phusion Hot Start II High-Fidelity PCR Master Mix (Thermofisher Scientific, Melbourne, Australia), 0.5 µM of each primer, and 100 ng of DNA template. PCR was performed in a SimpliAmp™ Thermal Cycler (Thermofisher Scientific, Melbourne, Australia), in a threestep protocol, using the following conditions: 98 • C initial denaturation in 1 min (1 cycle

Library Preparation, Quantification, Normalization, and Sequencing
Libraries were prepared using Nextera XT DNA Library Prep kit (Illumina, San Diego, CA, USA) following the manufacturer's protocols using the recommended input of 5 µL of 0.2 ng/µL gDNA per sample. This was followed by Sera-Mag™ SpeedBeads purification using 0.6× beads and two washes using 80% ethanol to select fragments >250 bp and remove unincorporated adapters. Each DNA library fragment size and concentration were determined using Agilent High Sensitivity D5000 reagents and ScreenTape on the Tape Station 4200 Instrument (Agilent Technologies, Santa Clara, CA, USA) according to the Agilent assay guide. Additionally, all individual libraries were quantified using QuantiFluor ® dsDNA System (Promega, Madison, WI, USA) to give an additional concentration estimate. The resultant size and concentration data from Tape Station and Quantifluor system were used to normalize each library to 4 nM by diluting with 10 mM Tris-HCl (pH 8.5) prior to pooling all samples. An equal volume of 5 µL was pooled and sequenced on an Illumina MiSeq benchtop sequencer (Illumina, Inc., San Diego, CA, USA), using a 500-cycle MiSeq Reagent Nano Kit v2 with a 10 pM input and 10% PhiX spike-in.

Bioinformatics and Next Generation Sequencing Data Analysis
Genomic data analysis was performed using the commercial bioinformatics program Geneious Prime software program 2022 v.1.1 (http://www.geneious.com, accessed on 15 January 2022) to analyze the sequences. The following reference sequences deposited in the NCBI database were used for comparative analysis: NC_037346.1, NC_037353.1, and NC_037341.1 for FASN, FABP4, and SCD genes, respectively. NGS sequenced data were retrieved from Illumina Dashboard-BaseSpace Sequence Hub (https://basespace.illumina. com/dashboard (accessed on 15 January 2022)) as paired read data in two separate forward and reverse read lists in FASTQ format. The retrieved raw reads were subjected to quality control measures. Reads were trimmed and adapters removed using the BBDuk trimmer in Geneious Prime 2022 v.1.1 with the default setting for paired-end reads. The Quality (Q) value of Phred score was set at 20 to improve sequenced data and increase the likelihood of calling true SNPs to 99%. Short reads with a minimum length of 20 bp were discarded, resulting in clean reads. Low coverage regions were excluded when calling SNPs using the Annotate and Predict→Find Low/High Coverage. The reads were mapped to reference in Geneious. The reference sequences were retrieved from the NCBI database (Genbank) of FASN, FABP4, and SCD of the Hereford bovine breed. The Sensitivity was set on the Medium Sensitivity/Fast and the Fine-Tuning (iterate up to 5 times) option was selected to improve the results by aligning reads to each other in addition to the reference sequence. Major allele frequencies from the next-generation sequencing data based on observed and expected genotypes were computed using the Hardy-Weinberg equilibrium principle.

Statistical Analyses
All statistical analyses of the associations between identified SNP of the SCD, FASN, and FABP4 genes with IMF, FMP, and fatty acid composition were performed using R software v.4.0.2 (R Foundation for Statistical Computing, Vienna, Austria). Descriptive summaries were presented as medians and inter-quartile ranges. Spearman correlations were used to estimate the relationships between SNP and meat quality traits (fatty acids, FMP, and IMF). Fatty acids, IMF, and FMP were analyzed as dependent variables using a non-parametric analysis of variance (Kruskal-Wallis test).
where γ ij is the quantitative phenotype (meat quality traits), µ is the overall mean, τ i is the SNP effect and ∈ ij is the random error. Differences between medians were compared using the Tukey-adjusted multiple comparisons test. The threshold for significance was set at p < 0.05.

Results
The candidate genes SCD, FASN, and FABP4 were targeted for this study because of their known physiological roles in intramuscular fat deposition and fatty acid composition. SNP identification in the sequences of SCD, FASN, and FABP4 lipogenic genes was made possible due to the successful primer design in Geneious Prime software, as depicted in Table 1 and corroborated by the gel images of the PCR products shown in Figures 1-3.

SCD, FASN, and FABP Variants and Genotypes
This study used the reference sequence extracted from the NCBI database for Hereford beef cattle, resulting in 12 identified SCD SNP with major allele frequencies ranging from 0.53 to 0.97, as shown in Table 2. A further identification of 18 SNP of the FASN gene with major allele frequencies ranging from 0.50 to 1.00 (Table 3) and five SNP in the FABP4 gene with major allele frequencies ranging from 0.50 to 0.98 is depicted in Table 4. It is evident from Table 2 that the Wagyu was heterozygous in SNP g.21267406 Genes 2022, 13, 1411 9 of 21 T>C, g.21267451 C>A and g.21267897 C>T, but homozygous in nine other SNP (g.21268157 G>A; g.21269020 A>G; g.21269495 G>C; g.21270740 A>G; g.21271264 C>A; g.21271726 C>G; g.21272247 A>G; g.21272423 C>T; g.21275852 C>A). Angus and Hereford were homozygous in most SCD SNP, except for SNP g.21268157 G>A ( Table 2). As depicted in Table 3, Wagyu was homozygous, while Angus was heterozygous in ten FASN SNP (g.50783867 C>A; g.50787010 A>G; g.50787138 A>G; g.50788320 C>T; g.50788576 T>C; g.50789409 T>C; g.50789883 T>C; g.50792348 C>T; g.50793358 A>G; g.50794100 T>C). The Hereford was heterozygous at five SNP loci (g.50787010 A>G; g.50789409 T>C; g.50789883 T>C; g.50792348 C>T; g.50794100 T>C). For the FABP4 gene, Angus was heterozygous in SNP g.44680268 C>G; whereas Hereford and Wagyu breeds were homozygous (Table 4).

Variable
FASN SNP g.50787138 A>G was associated with variations in FMP, oleic acid, ALA, EPA, DPA, EPA + DHA, EPA + DHA + DPA, and ALA + EPA + DHA + DPA (Figure 9). GG genotype had the lowest FMP (33 • C), but the highest n-3 LC-PUFA compared to AA and AG genotypes (Table 7). Significant differences existed between GG versus AG and GG versus AA genotypes for all variables (p < 0.05; Figure 9). There was no significant difference between AA and AG genotypes (Figure 9). In the FABP SNP g.44678794 G>A, only FMP differed between genotypes. AA genotype had the lowest FMP (35 • C), but both GG and AG were at 45 • C (Table 8). AA genotype differed significantly from AG (p = 0.056) and GG (p = 0.013) genotypes ( Figure 10). and AG genotypes (Table 7). Significant differences existed between GG versus AG and GG versus AA genotypes for all variables (p < 0.05; Figure 9). There was no significant difference between AA and AG genotypes (Figure 9). In the FABP SNP g.44678794 G>A, only FMP differed between genotypes. AA genotype had the lowest FMP (35 °C), but both GG and AG were at 45 °C (Table 8). AA genotype differed significantly from AG (p = 0.056) and GG (p = 0.013) genotypes ( Figure 10).

Discussion
The bioavailability of supplemented omega-3 fatty acids is significantly affected by biohydrogenation in ruminants [51] due to the conversion of unsaturated fatty acids to SFA moieties, which are absorbed in the small intestine and deposited in edible beef muscle tissues [52]. Genomic selection using breed-specific DNA markers may

Discussion
The bioavailability of supplemented omega-3 fatty acids is significantly affected by biohydrogenation in ruminants [51] due to the conversion of unsaturated fatty acids to SFA moieties, which are absorbed in the small intestine and deposited in edible beef muscle tissues [52]. Genomic selection using breed-specific DNA markers may considerably improve health-promoting n-3 LC-PUFA and meat eating quality traits in beef cattle and the accuracy of selection [53].

SCD Gene Polymorphism
The SCD gene plays a role in the biosynthesis of monounsaturated fatty acids (MUFA) from stearic and palmitic acids, by introducing a double bond between carbon atoms ∆ 9 and ∆ 10 [54]. SCD enzyme activity has been shown to be positively correlated with unsaturated fatty acid content and marbling in Angus [55], Hereford [31], Japanese Black [56], Jersey [57], and Italian Holstein [58] cattle. Most studies have focused on the association between SCD SNP and oleic acid, but in the present study, the finding that SCD SNP were significantly associated with health-beneficial n-3 LC-PUFA and meat quality traits (Tables 5 and 6, Figures 7 and 8) complements previous reports [31,44,54,59]. The fact that beef cattle breeds with heterozygous genotypes CT and AC in the SCD SNP g.21267406 T>C and g.21271264 C>A produced the highest oleic acid, ALA, EPA, DHA, DPA compared to CC, TT and AA genotypes (Tables 5 and 6), represents a genetic divergence between the three beef cattle breeds that provides insights into the role of the SCD gene in the metabolism of EPA, DHA, and DPA. This underscores the potential of SCD SNP in marker-assisted selection to improve the composition of health-promoting n-3 LC-PUFA and meat eating quality attributes in Australia's beef cattle grazing systems.

FASN Gene Polymorphism
Fatty acid synthase (FASN) catalyzes the de novo synthesis of palmitic and stearic acids [19,60]. Abe et al. [61] and Barton et al. [62] reported three FASN SNP loci (g.16024G>A, 16039T>C, and g.17924A>G), in exons 34 and 39, that caused an amino acid substitution of threonine to alanine in the thioesterase domain. These FASN polymorphisms were significantly associated with oleic, palmitic, and palmitoleic acids [19,62,63]. Yokota et al. [64] reported an additive effect of SNP g.16039T>C on several SFA and MUFA in the intramuscular fatty acids of Japanese Black cattle. The impact of FASN SNP g.17924A>G on fatty acid composition was also reported in Angus cattle [19], Hereford [63], and crossbred beef steers [65]. FASN SNP g.12870, g.13126, g.15532, g.16907, and g.17924 were reported to be associated with intramuscular fatty acids in Korean cattle. Results of the present study also demonstrated that polymorphisms in the FASN SNP g.50787138 A>G detected in intron 21 have positive associations with FMP, oleic acid, EPA, DHA, DPA, and total n-3 LC-PUFA, thus suggesting a synonymous mutation. Previous studies reported non-synonymous mutations where an amino acid substitution of threonine to alanine was observed [61,62]. In the present study, the finding that GG genotype had the lowest FMP and the highest oleic acid, ALA, EPA, DHA, DPA, and EPA+DHA+DPA compared to AG and AA genotypes (Table 7 and Figure 9) confirms that the g.50787138 A>G SNP can effectively be used as a selection marker due to its significant association with health-promoting n-3 LC-PUFA metabolism in pasture-based beef cattle.

FABP4 Gene Polymorphism
The fatty acid-binding protein-4 (FABP4) plays a role in the absorption, transport, and metabolism of long-chain fatty acids [66][67][68]. Previous studies reported an association between FABP4 polymorphism and carcass characteristics [35] and fat deposition [66,67]. In Japanese Black cattle, Hoashi et al. [69], reported that genotype 174V of the FABP4 gene was significantly associated with 16:1 and 18:2. FABP4 SNP g.2834C>G and g.3631G>A were associated with marbling [41]. Shin et al. [42] found a significant correlation of FABP4 SNP g.3691G>A with marbling and meat quality traits in Korean cattle. Other studies evaluated the associations of FABP4 SNP with carcass characteristics, 14:0, 16:1, 16:1n-7, LA, and oleic acid [62,67,69,70]. In the present study, FABP4 SNP g.44678794 G>A was significantly correlated with FMP, where the AA genotype had the lowest fat melting point of 35 • C (Table 8). This suggests that polymorphisms in intron 1 of the FABPA gene controlling fat melting point could be influential markers in specifically selecting for healthy polyunsaturated fatty acids in pasture-based Angus, Hereford, and Wagyu since genetic markers are often breed-specific and may not be extrapolated to other cattle breeds [46,62].

Conclusions
This study identified SCD SNP g.21267406 T>C and g.21271264 C>A in introns 2 and 4, FASN SNP g.50787138 A>G in intron 21 and FABP4 SNP g.44678794 G>A in intron 1 that were significantly associated with FMP, IMF, and health-beneficial n-3 LC-PUFA fatty acids. These results, taken together, provide an opportunity for SNP marker-assisted selection for improving meat eating quality and health-promoting n-3 LC-PUFA traits in Australian pasture-based Bowen Genetics Forest Pastoral beef cattle while they are young and still alive. As these data are actually laboratory tested, individualized, and customized to each cattle breed, as opposed to estimated breeding values based on herd averages, they increase the precision and accuracy of the selection program. They also enable Australian beef farmers to make early decisions in quantifying the genetic worth of their animals for meat eating quality and health-beneficial n-3 LC-PUFA based on muscle biopsy sampling and the development of a novel selection index for beef quality in their breeding program.