Issue cover
Research Article | 26 Feb 2026

Identification of novel fatty acid-binding protein 4 polymorphisms and association of g.5002C>T with longissimus dorsi thickness in Indonesian Bali cattle (Bos javanicus)

Dairoh Dairoh1 ORCID , Mokhamad Fakhrul Ulum2 ORCID , Sutikno Sutikno1 ORCID , Ahmad Furqon1 ORCID , Isyana Khaerunnisa1 ORCID , Cece Sumantri3 ORCID , and Jakaria Jakaria3 ORCID Show more
VETERINARY WORLD | Article No. 22 | pg no. 760-770 | Vol. 19, Issue 2 | DOI: 10.14202/vetworld.2026.760-770
Citations:

Cite this Article

  • APA
  • MLA
  • Chicago
  • Vancouver
  • Harvard

                            
                        

ABSTRACT

Background and Aim: The fatty acid-binding protein 4 (FABP4) gene is a key regulator of lipid metabolism and has been linked to carcass and meat quality traits in several cattle breeds. However, genetic variation in FABP4 and its phenotypic relevance in Bali cattle (Bos javanicus), an important indigenous Indonesian breed, remains poorly characterized. This study aimed to identify novel single-nucleotide polymorphisms (SNPs) in the FABP4 gene of Bali cattle and to evaluate their associations with in vivo ultrasound-measured carcass traits and fatty acid composition.

Materials and Methods: Blood samples were collected from 95 Bali cattle bulls sourced from a single population. Genomic DNA was extracted, and a 721-bp fragment of FABP4 spanning intron 2 to exon 4 was amplified and sequenced using the Sanger method. SNP detection and sequence alignment were performed using standard bioinformatics tools. Carcass and meat quality traits, including longissimus dorsi thickness (LDT), backfat thickness (BFT), marbling score (MS), and intramuscular fat (IMF), were measured in live animals using ultrasound imaging. Fatty acid composition was determined from post-mortem muscle samples using gas chromatography–flame ionization detection. Associations between FABP4 genotypes and phenotypic traits were analyzed using a general linear model, with statistical significance set at p < 0.05.

Results: Four novel FABP4 SNPs were identified: g.4631T>C, g.4724T>C, and g.4769G>A located in intron 3, and g.5002C>T located in exon 4. The g.5002C>T variant was a nonsynonymous substitution resulting in a valine-to-alanine amino acid change. Among the identified polymorphisms, only g.5002C>T showed a significant association with LDT (p = 0.0296). Animals carrying the CT genotype exhibited greater LDT compared with CC and TT genotypes. No significant associations were observed between any FABP4 SNPs and BFT, MS, IMF, or fatty acid composition.

Conclusion: This study provides the first evidence of FABP4 genetic variation in Bali cattle and demonstrates that the g.5002C>T polymorphism is significantly associated with LDT. The CT genotype at this locus may serve as a potential genetic marker for improving muscle development in Bali cattle, subject to validation in larger and independent populations.

Keywords: Bali cattle, Bos javanicus, FABP4 gene, fatty acid composition, meat quality traits, single-nucleotide polymorphism, ultrasound imaging.

INTRODUCTION

Meat quality and carcass traits are economically important in beef production because they directly affect processing efficiency, market value, and consumer acceptance [1]. Key attributes, including muscle development, tenderness, intramuscular fat (IMF), and fatty acid (FA) composition [2], play central roles in determining meat palatability and nutritional quality. Improvements in these traits can substantially enhance the market appeal of meat products by meeting consumer expectations for superior taste and texture while also providing health benefits through improved FA profiles [3]. Inoue et al. [4] further reported that additional meat quality traits, particularly tenderness, contribute to unit meat price alongside the marbling score (MS).

In Indonesia, Bali cattle (Bos javanicus) represent an important indigenous breed and a valuable animal genetic resource domesticated from Banteng (Bibos banteng) [5]. Bali cattle are well known for their adaptability to tropical environments, tolerance to limited feed resources, and significant contribution to national beef production. Compared with commercial breeds, B. javanicus exhibits superior meat quality, with carcass yields ranging from 54.07% to 55.61% [6]. These cattle are characterized by desirable meat texture and are highly preferred in the marketplace. The meat of B. javanicus is reported to be tender [7], and previous studies have demonstrated favorable meat quality with a healthy FA composition, including omega-3 and omega-6 fatty acids [8].

Meat quality and carcass traits are complex phenotypes regulated by multiple genes. The application of DNA markers in marker-assisted selection (MAS) has markedly advanced livestock breeding [9]. Developments in molecular genetics have enabled the identification of candidate genes associated with economically important traits, thereby improving breeding efficiency through MAS [10, 11]. This approach allows researchers to detect specific genetic variations linked to desirable phenotypes. However, integrative genotype–phenotype analyses combining in vivo ultrasound measurements and FA profiling have not previously been applied to B. javanicus.

The fatty acid-binding protein 4 (FABP4) gene plays a crucial role in lipid metabolism by promoting lipid hydrolysis and the intracellular transport of FAs among tissues [12]. In cattle, FABP4 is located on chromosome 14 and consists of four exons and three introns, encoding a protein that binds and transports FAs within cells [13]. FABP4 facilitates FA uptake, transport, and storage in adipocytes and regulates fat deposition through interactions with lipid signaling pathways, including PPARγ [14], thereby influencing fat accumulation and meat quality [15]. Recent studies have reported associations between single-nucleotide polymorphisms (SNPs) in FABP4 and marbling traits in Caracu × Nellore cows, as well as in Caracu, Canchim, and Angus sires [16], and with FA composition in Angus, Brahman, Creole, Hereford, Holstein, Limousine, and Nellore cattle [17, 18]. Nishimura et al. [19] demonstrated a significant association between IMF deposition and meat tenderness, underscoring the importance of FABP4 in determining meat quality. Through its effects on lipid hydrolysis and intracellular FA transport, FABP4 supports adipocyte development and differentiation [20], promoting IMF accumulation that enhances tenderness by altering muscle fiber structure and increasing juiciness and flavor [2123]. Furthermore, Bayraktar et al. [24] evaluated the expression of FABP3 and FABP4 in a native Indonesian cattle breed, highlighting their important roles in FA metabolism and local breed adaptation.

Despite growing evidence that FABP4 polymorphisms are associated with carcass traits, IMF deposition, and FA composition in several taurine and crossbred cattle populations, comparable information for Bali cattle remains very limited. In particular, the genetic diversity of FABP4 in Bali cattle and its relationship with carcass and meat quality traits assessed in vivo has not been clearly established. Most existing studies have examined either molecular variation or phenotypic traits in isolation, while integrated genotype–phenotype approaches combining in vivo ultrasonography and FA profiling are still lacking in Bali cattle. This gap constrains the identification of robust genetic markers that could be implemented in MAS programs to improve meat quality in this indigenous breed.

Therefore, this study aimed to identify and characterize SNPs in the FABP4 gene of Bali cattle and to evaluate their associations with in vivo ultrasound-measured carcass traits, including longissimus dorsi thickness (LDT), backfat thickness (BFT), MS, IMF, and FA composition. By integrating molecular genetic information with phenotypic measurements, the study sought to assess the potential of FABP4 polymorphisms as candidate markers for MAS and to provide foundational data to support genetic improvement and conservation strategies in Bali cattle.

MATERIALS AND METHODS

Ethical approval

The Animal Ethics Committee of the Banjarmasin City Food Security, Agriculture, and Fisheries Service approved the experimental procedures (approval number: 520/624/DKP3/X11/2021). All sampling procedures complied with established animal welfare guidelines. Gentle restraint, appropriately sized needles, and rapid venipuncture performed by trained personnel were used to minimize pain and stress during blood collection. The puncture site was cleaned before and after sampling to reduce discomfort and promote healing. Slaughter was conducted at a licensed slaughterhouse in accordance with Indonesian regulations for humane handling and slaughter of cattle.

Study period and location

The experiment was conducted from December 2022 to June 2023. A total of 95 Bali cattle bulls from a single population were included, and all available animals were used as representative samples. The animals weighed 250–350 kg and were 18–36 months of age. Bali cattle were sourced from Kupang, East Nusa Tenggara Province, Indonesia, and transported by ship to the Basirih slaughterhouse in Banjarmasin. Upon arrival, cattle were maintained under an intensive management system for 2 weeks before slaughter in South Banjarmasin, South Kalimantan. During this period, forage was provided at 10% of body weight, concentrate at 2%, and water ad libitum. Animals were restrained using a cattle chute to minimize movement during blood collection by a trained veterinarian.

Blood collection and DNA extraction

Blood samples were collected from the jugular vein using a Venoject with a sterile disposable 21-gauge needle (Becton, Dickinson and Company, Franklin Lakes, NJ, USA). Approximately 5 mL of blood was transferred into vacuum tubes containing 1.5 mL ethylenediaminetetraacetic acid (EDTA). Samples were immediately placed on ice and stored at 4°C until DNA extraction. Genomic DNA was extracted using the Geneaid Genomic DNA Mini Kit (Geneaid Biotech Ltd., Taipei, Taiwan) according to the manufacturer’s instructions. Extracted DNA was stored at −20°C. DNA purity was assessed spectrophotometrically, with acceptable A260/280 ratios of 1.8–2.2 and concentrations ≥20 ng/μL.

Polymerase chain reaction (PCR) amplification of FABP4

Primer sequences were designed using Primer3 and BLAST primer tools. The forward (5′-CCC TCC ATC ATT GTA ATC ACT-3′) and reverse (5′-GGA CAA CGT ATC CAG CAG AAA-3′) primers for the FABP4 gene amplified a fragment spanning intron 2 to exon 4. The expected PCR product length was 721 base pairs (bp). Amplification was performed using an AB System thermocycler (Applied Biosystems, Foster City, CA, USA) under the following conditions: pre-denaturation at 95°C for 1 min, followed by 35 cycles of denaturation at 95 °C for 15 s, annealing at 57°C for 15 s, extension at 72°C for 10 s, and a final extension at 72°C for 3 min.

Agarose gel electrophoresis and DNA sequencing

PCR products were separated on a 1% agarose gel stained with Florosafe (1st BASE, Singapore, Singapore) and visualized using a UV transilluminator (Bio-Rad, Hercules, CA, USA). Sequencing was performed at 1st BASE Laboratory Services (Selangor, Malaysia) using an ABI PRISM system with the BigDye Terminator kit v3.1 (Applied Biosystems). Sequencing quality was manually examined using FinchTV, BioEdit, and Molecular Evolutionary Genetic Analysis software. The first and last 30 bp were excluded because of low peak quality, and only SNPs with clear chromatogram peaks were retained. Sequence alignment was conducted against the Ensembl Bos taurus RefSeq transcript ENSBTAG00000037526, and SNP positions were determined using the Bos taurus ARS-UCD2.0: CM008181.2 reference genome assembly.

Ultrasound imaging measurements

In vivo carcass and meat traits were measured using a 7.1-MHz ultrasound transducer applied to the left rib area. A portable ultrasound device (SIUI CTS-800, SIUI Medical Systems, Shantou, China) was used to measure LDT, BFT, MS, and IMF at the 12th–13th rib position (Figure 1A) using transverse and longitudinal scans [25, 26]. IMF sonograms were obtained 3 days before slaughter. Hair at the scanning site was shaved, and ultrasound gel was applied to ensure adequate transmission of sound waves. Images were displayed in real time, saved, and analyzed using ImageJ software (National Institutes of Health, Bethesda, MD, USA). ImageJ was used to calibrate scale units to millimeters. Measurements of BFT, LDT (Figure 1B), and IMF (Figure 1C) were recorded at three points per animal. IMF values were subjected to simple linear regression to estimate the IMF percentage. Marbling was scored according to AUS-MEAT and MSA standards, with scores ranging from 0 to 9 (AUS-MEAT, 2018).

Figure 1

Figure 1. Ultrasound image analysis using ImageJ software in live Bali cattle. (A) Ribeye ultrasound image obtained between the 12th and 13th ribs. (B) Longitudinal ultrasound image used to measure backfat thickness and longissimus dorsi muscle thickness. (C) A transverse ultrasound image was used to estimate the percentage of intramuscular fat.

FA methyl ester preparation

Approximately 250 g of tenderloin muscle was collected from the right side of each carcass about 1 h post-mortem and stored at −20°C. From this sample, 40 g was used for FA analysis. Lipids were extracted using a chloroform–methanol solution, followed by transesterification to produce fatty acid methyl esters (FAMEs). FAMEs were extracted, centrifuged, dried in hexane, resuspended in chloroform to remove impurities, and purified by solid-phase extraction.

Gas chromatography–flame ionization detector conditions

FAMEs (1 μL) were analyzed using gas chromatography equipped with a flame ionization detector. Separation was achieved on a cyanopropyl methyl sil capillary column (60 m × 0.25 mm, 0.25 μm film thickness) with nitrogen as the carrier gas at 30 mL/min. Injector and detector temperatures were set at 220°C and 240°C, respectively. The oven temperature program was 125°C for 5 min, then increased to 225°C at 10, 5, and 3°C/min, with holding times of 5, 10, and 7 min, respectively. The split ratio was 1:80, the injection volume was 1 μL, and the linear velocity was 23.6 cm/s. FAMEs were quantified as relative percentages using a mixed FAME standard for retention-time and peak-area comparisons. Quality control was ensured by standard injections before sample analysis. Fatty acids were grouped into saturated fatty acids (SFA), monounsaturated fatty acids (MUFA), and polyunsaturated fatty acids (PUFA) [27]. Total unsaturated fatty acids (UFA) were calculated as the sum of MUFA and PUFA. Analyses were performed in duplicate, and instruments were calibrated according to the manufacturer’s guidelines.

Population genetics and SNP diversity analysis

Allele and genotype frequencies, observed and expected heterozygosity, and Hardy–Weinberg equilibrium were calculated using PopGen 1.32 software. Polymorphic information content (PIC) was calculated using the formula:

where pi and pj represent allele frequencies at a given locus.

Statistical analysis

Associations between FABP4 genotypes and carcass traits, meat quality traits, and FA composition were analyzed using a general linear model (GLM) in SAS version 9.4 (SAS Institute Inc., Cary, NC, USA). Duncan’s multiple-range test was used for post hoc comparisons. Statistical significance was set at p < 0.05. The GLM was expressed as:

Yij = μ + Gi + eij,

where Yij is the phenotypic observation, μ is the overall mean, Gi is the genotype effect, and eij is the random error. Carcass and meat traits, including FA composition, were corrected to 36 months of age, body weight, and similar environmental conditions using the equation:

Xi(corrected) = (X̅standard / X̅observation) × Xi(observation)

where Xi(corrected) is the corrected value, X̅standard is the standard group mean, X̅observation is the observation group mean, and Xi(observation) is the observed value.

RESULTS

FABP4 polymorphism

Genetic variation in FABP4 was identified by SNP analysis of a 721 bp DNA sequencing product. SNPs were detected and characterized by examining double peaks in sequencing chromatograms using FinchTV software (Figure 2). Four novel SNPs were identified in FABP4 from the sequencing map. Three SNPs (g.4631T>C, g.4724T>C, and g.4769G>A) were located in intron 3, whereas one SNP (g.5002C>T) was located in exon 4 (Table 1). Among these variants, g.5002C>T was a nonsynonymous substitution resulting in an amino acid change from valine to alanine, whereas the remaining three SNPs were synonymous and did not alter the amino acid sequence. Analysis of allele and genotype frequencies indicated that the AA genotype was more frequently observed for g.4724T>C, g.4769G>A, and g.5002C>T, whereas the BB genotype predominated for g.4631T>C (Table 2).

Figure 2

Figure 2. Partial sequencing maps of the FABP4 gene in Bali cattle.

Table 1. SNP information of the FABP4 gene

GeneSNPLocationVariation typedbSNPAmino acids
FABP4g.4631T>CIntron 3TransitionNovel
FABP4g.4724T>CIntron 3TransitionNovel
FABP4g.4769G>AIntron 3TransitionNovel
FABP4g.5002C>TExon 4TransitionNovelVal/Ala

Ala = Alanine, SNP = Single-nucleotide polymorphism, Val = Valine.

Table 2. Allelic and genotypic frequencies and diversity parameters of FABP4 SNPs.

GeneSNPNAAABBBABHoHeχ² testPIC
FABP4g.4631T>C950.000.080.920.040.960.0840.0810.160 ns0.077
FABP4g.4724T>C950.930.070.000.960.040.0740.0710.119 ns0.068
FABP4g.4769G>A950.780.210.010.880.120.2110.2060.051 ns0.184
FABP4g.5002C>T950.840.140.020.910.090.1370.1642.721 ns0.150

AA = Reference genotype (wild-type), AB = Heterozygous genotype, BB = Mutant genotype, FABP4 = Fatty acid-binding protein 4 gene, He = Expected heterozygosity, Ho = Observed heterozygosity, ns = Not significant (χ² test < 3.84 at α = 0.05, df = 2), PIC = Polymorphic information content, SNP = Single-nucleotide polymorphism.

Genetic association of FABP4 variants with meat quality traits

The associations between FABP4 SNPs and carcass and meat quality traits are summarized in Table 3. The SNP g.5002C>T showed a significant association with LDT (p < 0.05), suggesting a potential role in muscle development. In contrast, no significant associations were detected between FABP4 SNPs and BFT, MS, or IMF in Bali cattle. Notably, animals carrying the CT genotype at g.5002C>T exhibited significantly greater LDT compared with those carrying the CC and TT genotypes.

Table 3. Association of FABP4 SNPs with carcass ultrasound and meat characteristics in live Bali cattle.

SNPGenotype (N)LDT (mm)BFT (mm)MSIMF (%)
g.4631T>CTC (8)47.24 ± 6.111.84 ± 0.301.36 ± 0.192.22 ± 0.47
CC (83)47.16 ± 5.771.88 ± 0.311.51 ± 0.562.58 ± 1.39
g.4724T>CTT (84)47.23 ± 5.761.89 ± 0.311.51 ± 0.562.57 ± 1.39
TC (7)46.51 ± 6.211.79 ± 0.301.39 ± 0.192.28 ± 0.48
g.4769G>AGG (70)46.71 ± 5.121.89 ± 0.321.51 ± 0.572.57 ± 1.41
GA (20)49.36 ± 6.991.84 ± 0.281.42 ± 0.362.35 ± 0.91
AA (1)35.17 ± nc1.54 ± nc2.61 ± nc5.30 ± nc
g.5002C>TCC (77)46.61b ± 5.731.86 ± 0.311.51 ± 0.572.57 ± 1.43
CT (12)51.18a ± 4.422.00 ± 0.231.39 ± 0.232.28 ± 0.56
TT (2)44.67ab ± 7.401.77 ± 0.571.76 ± 0.413.20 ± 1.03

BFT = Backfat thickness, FABP4 = Fatty acid-binding protein 4 gene, IMF = Intramuscular fat, LDT = Longissimus dorsi thickness, MS = Marbling score, nc = Not counted. Means in the same column with different superscripts differ significantly (p < 0.05).

Genetic association of FABP4 variants with fatty acid composition

The relationship between FABP4 polymorphisms and FA composition in Bali cattle is presented in Table 4. None of the identified FABP4 SNPs showed a significant association with FA composition. Although FABP4 polymorphisms were hypothesized to influence FA profiles, the results indicate that these variants did not contribute significantly to FA variation in Bali cattle. This suggests that other genetic factors or environmental influences may play a more dominant role in determining FA composition. Further studies are required to elucidate the interactions between FABP4, other candidate genes, and non-genetic factors involved in the regulation of meat quality traits in Bali cattle.

Table 4. Association of FABP4 SNPs with FA composition in Bali cattle.

Fatty acid (%)g.4631T>C TC (1)g.4631T>C CC (43)g.4724T>C TT (43)g.4724T>C TC (1)g.4769G>A GG (35)g.4769G>A GA (9)g.5002C>T CC (40)g.5002C>T CT (4)
Fat content3.57 ± nc3.18 ± 1.233.18 ± 1.233.57 ± nc3.17 ± 1.223.26 ± 1.273.26 ± 1.222.43 ± 0.92
C8:00.00 ± nc0.06 ± 0.200.06 ± 0.200.00 ± nc0.06 ± 0.210.04 ± 0.110.06 ± 0.200.07 ± 0.14
C12:00.18 ± nc0.07 ± 0.020.07 ± 0.020.18 ± nc0.07 ± 0.030.07 ± 0.010.07 ± 0.030.06 ± 0.02
C13:00.04 ± nc0.03 ± 0.020.03 ± 0.020.04 ± nc0.03 ± 0.020.03 ± 0.010.03 ± 0.020.03 ± 0.03
C14:03.35 ± nc2.17 ± 0.532.17 ± 0.533.35 ± nc2.15 ± 0.552.37 ± 0.562.21 ± 0.582.00 ± 0.27
C14:10.09 ± nc0.28 ± 0.370.28 ± 0.370.09 ± nc0.28 ± 0.380.26 ± 0.310.29 ± 0.380.18 ± 0.29
C15:00.81 ± nc0.59 ± 0.360.59 ± 0.360.81 ± nc0.61 ± 0.370.55 ± 0.260.60 ± 0.370.52 ± 0.07
C16:032.95 ± nc20.82 ± 3.6120.82 ± 3.6132.95 ± nc21.10 ± 4.0921.05 ± 3.8921.16 ± 4.1720.43 ± 1.69
C16:12.00 ± nc1.30 ± 0.351.30 ± 0.352.00 ± nc1.31 ± 0.381.33 ± 0.311.33 ± 0.371.19 ± 0.12
C17:02.13 ± nc1.99 ± 0.771.99 ± 0.772.13 ± nc2.01 ± 0.831.92 ± 0.422.00 ± 0.761.96 ± 0.84
C17:10.44 ± nc0.26 ± 0.140.26 ± 0.140.44 ± nc0.27 ± 0.140.26 ± 0.160.27 ± 0.140.22 ± 0.17
C18:036.55 ± nc31.91 ± 3.6431.91 ± 3.6436.55 ± nc32.00 ± 4.0132.09 ± 2.0132.11 ± 3.6531.05 ± 4.32
C18:1n9c0.10 ± nc12.45 ± 4.5112.45 ± 4.510.10 ± nc11.64 ± 5.0014.21 ± 3.6512.06 ± 5.0413.25 ± 1.70
C18:1n9t1.72 ± nc3.03 ± 1.253.03 ± 1.251.72 ± nc2.86 ± 1.253.56 ± 1.173.07 ± 1.282.27 ± 0.71
C18:2n6c3.88 ± nc1.75 ± 0.761.75 ± 0.763.88 ± nc1.83 ± 0.911.67 ± 0.231.75 ± 0.782.27 ± 1.18
C18:3n30.00 ± nc0.23 ± 0.270.23 ± 0.270.00 ± nc0.21 ± 0.270.28 ± 0.260.21 ± 0.260.41 ± 0.29
SFA77.19 ± nc58.56 ± 5.9958.56 ± 5.9977.19 ± nc58.98 ± 7.0559.01 ± 4.4259.16 ± 6.6057.16 ± 6.66
MUFA10.40 ± nc20.31 ± 4.3220.31 ± 4.3210.40 ± nc19.47 ± 4.5022.48 ± 4.0220.04 ± 4.6820.61 ± 2.90
PUFA5.14 ± nc17.65 ± 4.7517.65 ± 4.755.14 ± nc16.71 ± 5.1219.91 ± 4.1117.38 ± 5.2717.27 ± 2.30
UFA5.26 ± nc2.66 ± 1.092.66 ± 1.095.26 ± nc2.76 ± 1.262.57 ± 0.512.66 ± 1.123.34 ± 1.34
Total FA87.60 ± nc78.66 ± 7.0378.66 ± 7.0387.60 ± nc78.33 ± 7.2080.96 ± 6.5179.01 ± 7.1077.44 ± 7.67

FABP4 = Fatty acid-binding protein 4 gene, FA = Fatty acid, MUFA = Monounsaturated fatty acid, nc = Not counted, PUFA = Polyunsaturated fatty acid, SFA = Saturated fatty acid, UFA = Unsaturated fatty acid. Means in the same row with different superscripts differ significantly (p < 0.05).

DISCUSSION

Genetic diversity and population structure

Genetic diversity is a key determinant of adaptability and long-term sustainability in livestock populations. Based on allele and genotype frequency data, the AA genotype was more frequently observed for SNPs g.4724T>C, g.4769G>A, and g.5002C>T, whereas the BB genotype predominated for g.4631T>C (Table 2). Adequate genetic variation is essential for breeding programs aimed at improving performance and ensuring food security [28, 29]. Preservation of genetic diversity reduces the risks associated with inbreeding and genetic uniformity, thereby maintaining resilience, fitness, and protection against inbreeding depression [30].

All SNPs identified in this study were polymorphic, with major allele frequencies 1% [31]. The observed heterozygosity (Ho) was 50% for all SNPs, indicating relatively low heterozygosity in the Bali cattle population. This pattern may reflect inbreeding, genetic drift, and a small effective population size, which can increase homozygosity and allele fixation [32]. Reduced genetic variation may negatively affect overall health and adaptability, increasing susceptibility to diseases and environmental challenges [30]. All SNPs conformed to Hardy–Weinberg equilibrium, suggesting stable allele and genotype frequencies that were not significantly influenced by selection, migration, mutation, or genetic drift [33]. This equilibrium indicates population stability and limited disturbance from breeding practices or environmental pressures [34]. In addition, all SNPs in FABP4 exhibited low PIC values, consistent with the limited genetic variability commonly reported in indigenous breeds with narrow genetic bases [35, 36].

Association of FABP4 polymorphisms with meat quality traits

The FABP4 gene has been widely investigated for its role in regulating meat quality traits in cattle, particularly those related to lipid metabolism and muscle development. Genes involved in these pathways influence key traits such as IMF deposition and muscle thickness. In the present study, SNP g.5002C>T showed a significant association with LDT (p < 0.05), indicating its potential involvement in muscle development. This nonsynonymous substitution results in an amino acid change that may alter FABP4 protein function, supporting its candidacy as a genetic marker for selective breeding in Bali cattle.

Previous studies have reported significant associations between FABP4 SNPs and MS and BFT, particularly for variants located in exon 3 [37]. Gao et al. [38] observed associations between FABP4 diversity and LDT as well as ultrasonographic IMF in Qinchuan cattle. Associations with hump height and eye muscle area have also been reported in northern Australian crossbred cattle [39], while Shin et al. [40] demonstrated relationships between FABP4 variants and MS in Hanwoo cattle. In contrast, the lack of association between SNPs and MS, IMF, and BFT in the present study may be attributed to unbalanced genotype distributions. For example, only one AA genotype was observed for SNP g.4769G>A, precluding robust association analysis. Differences in population structure, sample size, and genetic background across studies can lead to variable allele and genotype frequencies and influence statistical outcomes [41].

Association of FABP4 polymorphisms with fatty acid composition

The absence of significant associations between FABP4 polymorphisms and FA composition in Bali cattle contrasts with findings reported by Bayraktar et al. [24], who demonstrated strong associations between FABP4 variants and milk composition and FA profiles in Native Southern Yellow cattle. These observations highlight tissue- and breed-specific functions of the FABP gene family, suggesting that the influence of FABP4 on lipid metabolism may be more pronounced in milk than in muscle or adipose tissues of Bali cattle.

Other studies have reported significant relationships between FABP4 variants and specific FAs. Hoashi et al. [42] identified associations between FABP4 polymorphisms and palmitoleic acid (C16:1) and linoleic acid (C18:2n6c) in Qinchuan cattle. Palmitoleic acid has been linked to improved metabolic health, including reduced insulin resistance [43, 44], while linoleic acid is known for its cardiovascular benefits [45, 46]. Yin et al. [47] identified six FABP4 SNPs associated with fat content and MS in Yellow Yanbian cattle, and Oh et al. [48] reported associations with oleic acid (C18:1n9c) and monounsaturated FA content. The lack of association observed in Bali cattle may reflect the polygenic nature of FA composition, limited genetic diversity within the population, and strong environmental influences such as diet, management system, and physiological status, which can mask genetic effects [49, 50].

Implications for marker-assisted selection

The SNP g.5002C>T in FABP4 represents a promising candidate for MAS in Bali cattle. However, validation in larger and geographically diverse populations under varying environmental conditions is required before practical implementation. Following validation, incorporation of this SNP into national breeding programs could facilitate the identification of animals with favorable genotypes, improve meat quality efficiency, and accelerate genetic improvement in Bali cattle [51].

CONCLUSION

This study identified four novel SNPs in the FABP4 gene of Bali cattle, three located in intron 3 and one nonsynonymous variant (g.5002C>T) in exon 4. Among the detected variants, only g.5002C>T showed a significant association with LDT (p < 0.05), with the CT genotype exhibiting superior muscle thickness compared with CC and TT genotypes. No significant associations were observed between FABP4 SNPs and BFT, MS, IMF, or FA composition, indicating that the genetic effect of FABP4 in Bali cattle is trait-specific and primarily related to muscle development rather than fat deposition or FA profiles.

The significant association between FABP4 g.5002C>T and LDT highlights its potential utility as a genetic marker for improving muscle-related carcass traits in Bali cattle. Incorporation of this SNP into MAS schemes could support more efficient selection strategies aimed at enhancing meat yield without adversely affecting other quality traits. This is particularly relevant for breeding programs focused on indigenous cattle improvement under tropical production systems.

A major strength of this work lies in the integration of molecular genetics with in vivo phenotyping, combining SNP discovery, ultrasound-based carcass measurements, and FA profiling. This comprehensive genotype–phenotype approach provides a robust framework for evaluating candidate genes in indigenous cattle populations. Additionally, the identification of previously unreported FABP4 variants expands the genomic resource base for Bali cattle.

Despite these strengths, the study was conducted on a single population with a relatively limited sample size, which may have constrained the detection of associations for traits with low genetic variance or unbalanced genotype frequencies. The low heterozygosity observed for several SNPs may also have reduced statistical power, particularly for traits such as MS, IMF, and FA composition that are strongly influenced by multiple genes and environmental factors.

Future studies should validate the effect of FABP4 g.5002C>T in larger and geographically diverse Bali cattle populations and across different management and feeding systems. Integration of this marker with additional candidate genes and genomic tools could further enhance MAS accuracy. Functional studies examining the biological impact of the g.5002C>T amino acid substitution on FABP4 activity would also help clarify the mechanistic basis of its association with LDT.

In conclusion, this study provides the first evidence that FABP4 genetic variation contributes to muscle development in Bali cattle, with g.5002C>T emerging as a promising marker associated with LDT. While FABP4 polymorphisms were not linked to fat-related traits or FA composition, the findings underscore the importance of breed-specific genetic effects and support the targeted use of FABP4 in future genetic improvement and conservation programs for Bali cattle.

DATA AVAILABILITY

All the generated data are included in the manuscript.

AUTHORS’ CONTRIBUTIONS

DD and JJ: Conceptualization of the study, investigation, methodology, data analysis, laboratory work, interpretation of data, and drafted, reviewed, and edited the manuscript. DD, MFU, and JJ: Collected the animal sample and performed ultrasound imaging. SS, AF, and IK: Validation, investigation, and reviewed and edit the manuscript. CS: Data analysis and revised the manuscript. All authors have read and approved the final version of the manuscript.

COMPETING INTERESTS

The authors declare that they have no competing interests.

PUBLISHER’S NOTE

Veterinary World remains neutral with regard to jurisdictional claims in the published institutional affiliations.

ACKNOWLEDGMENTS

This research was funded by PMDSU program (contract number: 3677/IT3/L1/PT/01.03/P/B/2022) years 2022 from the Ministry of Research and Technology or granted by the Ministry of Education and Culture Republic of Indonesia (by Dairoh); Postdoctoral Fellowship at the Research Center for Applied Zoology, National Research and Innovation Agency (BRIN), Indonesia No: B-3464/III.5/SI.06/6/2025 (By Dairoh), and Penelitian Fundamental Reguler (PFR) from the Ministry of Education, Culture, Research, and Technology No: 006/C3/DT.05.00/PL/2025 (by Jakaria). The authors would like to express their gratitude to the field personnel for their support and to the National Research and Innovation Agency’s Talent Management for providing the postdoctoral program.

REFERENCES

  1. Desmukh SS, Elayadeth-Meethal M, Abraham BL, Asaf M, Lali FA. Sustainable livestock production:Prospects for innovation. Hyderabad: MANAGE; 2023. p. 111. [Google Scholar]
  2. Vázquez-Mosquera JM, Fernandez-Novo A, de Mercado E, Vázquez-Gómez M, Gardon JC, Pesántez-Pacheco JL. Beef nutritional characteristics, fat profile and blood metabolic markers from purebred Wagyu, crossbred Wagyu and crossbred European steers raised on a fattening farm in Spain. Animals 2023;13(5):846. [Google Scholar]
  3. Baik M, Lee J, Kim SY, Ranaweera KKTN. Factors affecting beef quality and nutrigenomics of intramuscular adipose tissue deposition. Anim Biosci 2023;36(2):350-363. [Google Scholar]
  4. Inoue K, Shoji N, Hinda T, Oyama K. Genetic relationships between meat quality traits and fatty acid composition in Japanese Black cattle. Anim Sci J 2017;88(1):11-18. [Google Scholar]
  5. Martojo H. Indigenous Bali cattle is most suitable for sustainable small farming in Indonesia. Reprod Domest Anim 2012;47(1):10-14. [Google Scholar]
  6. Tahuk PK, Oktavianus R, Nahak R, Gerson FB. The effect of complete feed to carcass characteristics and meat quality of male Bali cattle fattened in West Timor, Indonesia. Vet World 2020;13(11):2515-2527. [Google Scholar]
  7. Dairoh, Ulum MF, Jakaria, Sumantri C. Calpain 1 gene expression in liver tissue and the association of novel calpain 1 single-nucleotide polymorphisms with meat quality in Bali cattle. Arch Anim Breed 2024;67((3)):311-321. [Google Scholar]
  8. Dairoh, Sutikno, Ishak AB, Priyanto R, Sumantri C, Ulum MF, Jakaria. Fatty acid profiling of Bali and Wagyu cattle using principal component analysis. Bull Anim Sci 2024;48(1):64-69. [Google Scholar]
  9. Sharma P, Doultani S, Hadiya KK, George LB, Highland HN. Overview of marker-assisted selection in animal breeding. J Adv Biol Biotechnol 2024;27(5):303-318. [Google Scholar]
  10. Meuwissen T, Hayes B, Goddard M. Genomic selection:A paradigm shift in animal breeding. Anim Front 2016;6(1):6-14. [Google Scholar]
  11. Ribeca C, Bonfatti V, Cecchinato A, Albera A, Gallo L, Carnier P. Effect of polymorphisms in candidate genes on carcass and meat quality traits in double muscled Piemontese cattle. Meat Sci 2014;96(3):1376-1383. [Google Scholar]
  12. Josephrajan A, Hertzel AV, Bohm EK, McBurney MW, Imai SI, Mashek DG. Unconventional secretion of adipocyte fatty acid binding protein 4 is mediated by autophagic proteins in a sirtuin-1-dependent manner. Diabetes 2019;68((9)):1767-1777. [Google Scholar]
  13. Zhang Y, Zhang J, Ren Y, Lu R, Yang L, Nie G. Tracing the evolution of fatty acid-binding proteins in organisms with heterogeneous fat distribution. FEBS Open Bio 2020;10(5):861-872. [Google Scholar]
  14. Furuhashi M. Fatty acid-binding protein 4 in cardiovascular and metabolic diseases. J Atheroscler Thromb 2019;26(3):216-232. [Google Scholar]
  15. Xiang J, Li H, Guo Z, Li T, Yamada T, Li X. Effect of FABP4 gene polymorphisms on fatty acid composition, chemical composition, and carcass traits in Sonid sheep. Animals 2025;15(2):226. [Google Scholar]
  16. Blecha IM, Siqueira F, Ferreira AB, Feijó GL, Torres Junior RA, Medeiros SR. Identification and evaluation of polymorphisms in FABP3 and FABP4 in beef cattle. GMR 2015;14(4):16353-16363. [Google Scholar]
  17. Dawood M, Kramer LM, Shabbir MI, Reecy JM. Genome-wide association study for fatty acid composition in American Angus cattle. Animals 2021;11(18):2424. [Google Scholar]
  18. Goszczynski DE, Papaleo-Mazzucco J, Ripoli MV, Villarreal EL, Rogberg-Muñoz A, Mezzadra CA. Genetic variation in FABP4 and evaluation of its effects on beef cattle fat content. Anim Biotechnol 2017;28(3):211-219. [Google Scholar]
  19. Nishimura T. The role of intramuscular connective tissue in meat texture. Anim Sci J 2010;81(1):21-27. [Google Scholar]
  20. Mallick R, Basak S, Das RK, Banerjee A, Paul S, Pathak S, Duttaroy AK. Fatty acids and their proteins in adipose tissue inflammation. Cell Biochem Biophys 2024;82:35-51. [Google Scholar]
  21. Hocquette JF, Gondret F, Baéza E, Médale F, Jurie C, Pethick DW. Intramuscular fat content in meat-producing animals:Development, genetic and nutritional control, and identification of putative markers. Animal 2010;4(2):303-319. [Google Scholar]
  22. Listrat A, Lebret B, Louveau I, Astruc T, Bonnet M, Lefaucheur L. How muscle structure and composition influence meat and flesh quality. Sci World J 2016;2016:1-14. [Google Scholar]
  23. Chen J, Chen F, Lin X, Wang Y, He J, Zhao Y. Effect of excessive or restrictive energy on growth performance, meat quality, and intramuscular fat deposition in finishing Ningxiang pigs. Animals 2021;11(1):27. [Google Scholar]
  24. Bayraktar M, Goncu S, Ergul A, Karaman R, Ozcan BD, Ergul S. Association between FABP3 and FABP4 genes with changes in milk composition and fatty acid profiles in native Southern Yellow cattle. Vet Sci 2025;12(19):893. [Google Scholar]
  25. Guidelines for uniform beef improvement program. North Carolina: North Carolina State University; 2016. p. 35-46. [Google Scholar]
  26. Jakaria, Khasanah H, Priyanto R, Baihaqi M, Ulum MF. Prediction of meat quality in Bali cattle using ultrasound imaging. Indones Trop Anim Agric 2017;42(2):59-65. [Google Scholar]
  27. Dairoh, Ulum MF, Jakaria, Ishak ABL, Sumantri C. A novel SNPs of the SREBF1 and SCARB1 genes and the association with fatty acid profile in Bali cattle. Trop Anim Sci J 2023;46(4):428-438. [Google Scholar]
  28. Blackburn HD, Wilson CS, Krehbiel B. Conservation and utilization of livestock genetic diversity in the United States of America through gene banking. Diversity 2019;11(12):244. [Google Scholar]
  29. Martyniuk E. Policy effects on the sustainability of animal breeding. Sustainability 2021;13(14):7787. [Google Scholar]
  30. Hoban S, Arche FI, Bertol LD, Bragg JG, Breed MF, Bruford MW. Global genetic diversity status and trends:Towards essential biodiversity variables for genetic composition. Biol Rev Camb Philos Soc 2022;97(4):1511-1538. [Google Scholar]
  31. Allendorf FW, Luikart G, Aitken SN. Conservation and the genetics of populations. Hoboken: Wiley-Blackwell; 2013. p. 602. [Google Scholar]
  32. Kardos M, Taylor HR, Ellegren H, Luikart G, Allendorf FW. Genomics advances the study of inbreeding depression in the wild. Evol Appl 2016;9(10):1205-1218. [Google Scholar]
  33. Abramovs N, Brass A, Tassabehji M. Hardy-Weinberg equilibrium in the large-scale genomic sequencing era. Front Genet 2020;11:210. [Google Scholar]
  34. Boes J, Boettcher P, Honkatukia M. Innovations in cryoconservation of animal genetic resources. Rome: FAO; 2023. p. 41-47. [Google Scholar]
  35. Harish A, Lopes Pinto FA, Eriksson S, Johansson AM. Genetic diversity and recent ancestry based on whole-genome sequencing of endangered Swedish cattle breeds. BMC Genomics 2024;25(1):89. [Google Scholar]
  36. Strucken EM, Gebrehiwot NZ, Swaminathan M, Joshi S, Al Kalaldeh M, Gibson JP. Genetic diversity and effective population sizes of thirteen Indian cattle breeds. Genet Sel Evol 2021;53(1):47. [Google Scholar]
  37. Ardicli S, Samli H, Alpay F, Dincel D, Soyudal B, Balci F. Association of SNP in FABP4 gene with carcass characteristics and meat quality in Holstein bulls. Ann Anim Sci 2017;17(1):117-130. [Google Scholar]
  38. Gao YY, Cheng G, Cheng ZX, Bao C, Yamada T, Cao GF. Association of variants in FABP4, FASN, SCD, SREBP1 and TCAP genes with intramuscular fat and carcass traits in Qinchuan cattle. Meat Sci 2022;192:108882. [Google Scholar]
  39. Mwangi FW, Pewan SB, Otto JR, Adegboy O, Charmley E, Gardiner CP. SNPs in FABP4, FASN and SCD genes influence carcass traits of tropical crossbred beef steers. Agriculture 2022;12(8):1171. [Google Scholar]
  40. Shin SC, Heo JP, Chung ER. Genetic variants of FABP4 gene associated with marbling score and meat quality in Hanwoo cattle. Mol Biol Rep 2012;39(5):5323-5330. [Google Scholar]
  41. Serdar CC, Cihan M, Yücel D, Serdar MA. Sample size, power and effect size revisited. Biochem Med 2021;31((1)):010502. [Google Scholar]
  42. Hoashi S, Hinenoya T, Tanaka A, Ohsaki H, Sasazaki S, Taniguchi M. Association between fatty acid composition and FABP4 and LXR-alpha genotypes in Japanese Black cattle. BMC Genet 2008;9(1):84. [Google Scholar]
  43. Bergman BC, Howard D, Schauer IE, Maahs DM, Snell-Bergeon JK, Clement TW. Importance of palmitoleic acid to insulin sensitivity. J Clin Endocrinol Metab 2013;98(1):E40-E50. [Google Scholar]
  44. Bermúdez MA, Pereira L, Fraile C, Valerio L, Balboa MA, Balsinde J. Roles of palmitoleic acid and positional isomers in inflammation and metabolic diseases. Cells 2022;11(14):2146. [Google Scholar]
  45. Anzhany D, Toharmat T, Despal, Lozicki A. Fatty acid biohydrogenation, fermentation, and digestibility of ration containing napier and king grass with different harvest ages and altitudes:in vitro study. Trop Anim Sci J 2023;47(1):68-78. [Google Scholar]
  46. Mercola J, D'Adamo CR. Linoleic acid and chronic disease associations. Nutrients 2023;15(14):3129. [Google Scholar]
  47. Yin BZ, Fang JC, Zhang JS, Zhang LM, Xu C, Xu H. SNPs in FABP4 correlate with meat quality and lipid metabolism in Yanbian yellow cattle. PLOS One 2020;15(6):1-14. [Google Scholar]
  48. Oh D, Lee Y, La B, Yeo J. SNP identification for fatty acid composition associated with FABP4 in Korean cattle. Asian-Australas J Anim Sci 2012;25(7):913-920. [Google Scholar]
  49. Kostusiak P, Slósarz J, Puppel K, Gołębiewski M, Grodkowski G. Polymorphism of genes and impact on beef quality. Curr Issues Mol Biol 2023;46(6):4749-4762. [Google Scholar]
  50. Sakowski T, Grodkowski G, Gołębiewski M, Puppel K, Kostusiak P, Slósarz J. Genetic and environmental determinants of beef quality. Front Vet Sci 2022;9:819605. [Google Scholar]
  51. Neumann GB, Korkuć P, Arends D, Wolf MJ, May K, Reißmann M. Design and performance of bovine 200k SNP chip for endangered cattle. BMC Genomics 2021;22(1):905. [Google Scholar]