Fat quality is determined by the composition of fatty acids. Genetic relationships between this composition and single nucleotide polymorphisms (SNPs) in the stearoyl-CoA desaturase1 (SCD1) gene were examined using 513 Korean native cattle. Single and epistatic effects of 7 SNP genetic variations were investigated, and the multifactor dimensionality reduction (MDR) method was used to investigate gene interactions in terms of oleic acid (C18:1), mono-unsaturated fatty acids (MUFAs) and marbling score (MS). The g.6850+77 A>G and g.14047 C>T SNP interactions were identified as the statistically optimal combination (C18:1, MUFAs and MS permutation p-values were 0.000, 0.000 and 0.001 respectively) of two-way gene interactions. The interaction effects of g.6850+77 A>G, g.10213 T>C and g.14047 C>T reflected the highest training-balanced accuracy (63.76%, 64.70% and 61.85% respectively) and was better than the individual effects for C18:1, MUFAs and MS. In addition, the superior genotype groups were AATTCC, AGTTCC, GGTCCC, AGTCCT, GGCCCT and AGCCTT. These results suggest that the selected SNP combination of the SCD1 gene and superior genotype groups can provide useful inferences for the improvement of the fatty acid composition in Korean native cattle.