Identification of Stearoyl-CoA Desaturase (SCD) Gene Interactions in Korean Native Cattle Based on the Multifactor-dimensionality Reduction Method.

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.

INTRODUCTION

Fat quality in beef market is determined by its fatty acid composition and fat tissue containing mono-unsaturated fatty acids (MUFAs), which have a lower melting point and thus contribute to the flavor and tenderness of beef (Miyazaki et al., 1979; Melton et al., 1982; Matsumoto et al., 2007). Oleic acid (C18:1) is one of the major MUFAs in beef fat and thought to be the source of the aroma in cooked beef (Melton et al., 1982; Mandell et al., 1998). For high-quality beef, it is important to improve its fatty acid composition, particularly its C18:1. In addition, the fatty acid composition is a heritable trait with heritability ranging from 0.31 to 0.73 (Inoue et al., 2008).

The stearoyl-CoA desaturase (SCD) gene, an enzyme in the chemical conversion of saturated fatty acids to unsaturated fatty acids, has been reported as a main factor for change of stearic acid (C18:0) and palmitic acid (C16:0) to C18:1 on unsaturated fatty acid of w-9; and many studies have attempted to increase unsaturated fatty acids for the improved flavor of meat and health benefits (Kuchel et al., 2004; Scollan et al., 2006; Alexander et al., 2007; Erkkila et al., 2008; Webb and O’Neill, 2008; Zhang et al., 2008). Recently, genetic relationships between the fatty acid composition of beef and multiple nucleotide sequence variants have been examined in encoding the SCD1 gene (Taniguchi et al., 2004). Three variants c.702 Aunsaturated fatty acids and MUFAs in the adipose tissue of Wagyu, Canadian Holstein, Jersey breeds, Fleckvieh bulls and Korean native cattle (p<0.05) (Kgwatalala et al., 2007; Milanesi et al., 2008; Tsuji, 2008; Barton et al., 2010; Ohsaki et al., 2009; Oh et al., 2011). SNPs in the UTR region of the SCD1 gene are significantly related to carcass and beef traits without an evaluation of fatty acids in Hanwoo (Shin et al., 2006). Based on the above discussion, the SCD gene, which is closely related to C18:1 as an important factor influencing the beef flavor, is considered to improve beef quality. This indicates that it is important not only to reveal the independent effect of each factor but also to analyze the combined effects. However, previous studies have paid little attention to single-factor effects and have not always considered independent and combined effects.

The multifactor dimensionality reduction (MDR) method has been employed (Ritchie et al., 2001; Chung et al., 2006) as a technique for efficiently detecting multiple genes and their interaction effects. This method is designed to address high-dimensional data and uncover complex relationships without relying on specified models fitting interactions between multiple genes (Bastione et al., 2004). With the increased availability of SNP information from an increasing number of SNPs in livestock genome databases (www.animalgenome.org) as well as human and other mammalian genomes, functional SNPs represent a better option for the simultaneous identification of several SNPs based on high-throughput tools such as DNA chips (Barendse et al., 2007).

In this study, the MDR method was used to detect superior SNP combinations influencing C18:1, MUFAs and marbling score (MS) in Korean native cattle. Given that the detection of genetic information on 15 SNPs in the intron, exon5 and 3′UTR regions of the SCD1 gene are significantly associated with fatty acid composition and MS (Oh et al., 2011), the interaction effects of SNPs in the SCD1 gene on C18:1, MUFAs and MS were examined. Here multiple SNPs were used by employing the MDR method to test the main and interaction effects of multiple SNPs on the meat quality of Korean native cattle, and a number of SNP factors believed to affect C18:1, MUFAs and MS were considered. Finally, the superior genotype groups were explored based on interactions between candidate genes.

MATERIALS AND METHODS

Animals and the phenotype

A total of 513 Korean native cattle (called Hanwoo in Korea) were bred in Gyeongsangbuk-Do, Korea. The cattle were fed according to the feeding program of each farmer (n = 14) and were the progeny of 18 sires. The pedigree records of these cattle were collected by the Charmpoom Hanwoo Research Institution. In general, they were weaned at 6 months of age, castrated at 6 months of age, fed growth-stage feed for 18 months, and given a high-concentration diet in the last 6 months. All were slaughtered at 941±72 d of age, and the marbling score was measured 24 h after slaughter. First, the carcass was dissected at the last rib and the first lumber vertebra according to the Animal Product Grading System of Korea. The MS was measured in the left carcass cut across the vertebra between the last thoracic vertebra and the first lumbar vertebra. The marbling degree was scored from 1 to 9 (mean = 5.43) such that the higher the score, the more abundant the intramuscular fat. Genomic DNA was extracted from the longissimus muscle by using the LaboPass Tissue Mini kit (Cosmo Genetech, Seoul, Korea).

Total lipids were extracted from approximately 500 mg of the longissimus muscle with chloroform/methanol (2/1, v/v) (Folch et al., 1957) and then methylated with sodium methylate (O’Keefe et al., 1968). They were filtered through a filter paper in a water bath (40°C). The filtrate was mixed with distilled water, from which a layer of methanol and water was removed. After the removal of chloroform and lipid layers by using nitrogen gas, the sample was treated with BF3-methanol (14%) and then subjected to transmethylation at 65°C. Fatty acid contents were analyzed using gas-chromatography (PerkinElmer, Inc., Waltham, MA, USA).

SNP genotyping

Fifteen SNPs of the SCD1 gene were selected based on the reference SNP of the bovine SCD1 gene (GenBank No. AY241932) in dbSNP. Genotypes at the 15 SNPs were subjected to a preliminary analysis. Primers for their amplification and extension were designed for single-base extension (SBE) (Vreeland et al., 2002) by using forward, reverse, and extension primer sequences (Supplementary Table 1). Primer extension reactions were performed using the SNaPshot ddNTP Primer Extension Kit (Applied Biosystems, Foster City, CA, USA). For the purification of the reaction of the primer extension, mixtures containing exonuclease 1 and shrimp alkaline phosphatase were added to the reaction mixtures. Samples were cultured at 37°C for 1 h and then inactivated at 72°C for 15 min. The PCR product was mixed with the Genescan 120 LIZ standard and HiDi formamide (Applied Biosystems), followed by denaturation at 95°C for 5 min. Electrophoresis was performed using the ABI PRISM 3130XL Genetic Analyzer, and then the electrophoresis product was assayed using GeneMapper v.4.0 (Applied Biosystems).

Supplementary Table 1

The information of primer sequence about SNPs of SCD (Accession number : NM_AY241932) gene of BTA 26 in Korean native cattle

SNP	Region		Sequences	Temperature	Product size
g.6850+77 A>G	Intron	F1	ACTCTGGCCTGAGGAAAACA	60	348
		R2	CCGGAACTTAACCACAAGGA
		E3	GTGGCGGAATGGAGGCGATC
g.8243-72 A>G		F	AGGCATCAGAATGGAACTGC	60	401
		R	AGACACCCTCTCAGGGGAAT
		E	AGCTGTGGGCTCTGGAAGAC
g.8646+128 A>G		F	AGCTGGTGATGTTCCAGAGG	60	390
		R	TTTGTGCCTCTCCTCGCTAT
		E	AAATTACAATTTAAAACCCT
g.10153 A>G	Exon 5	F	CAGAAAATTTCCTTGCCCATT	60	552
		R	TGTTGCTTAACTTTCAAGGGTTT
		E	ATCCTGCCCACACTCGTGCC
g.10213 T>C		F	CAGAAAATTTCCTTGCCCATT	60	552
		R	TGTTGCTTAACTTTCAAGGGTTT
		E	TTTGCCACCTTATTCCGTTA
g.10329 C>T		F	CAGAAAATTTCCTTGCCCATT	60	552
		R	TGTTGCTTAACTTTCAAGGGTTT
		E	TCTGGTTTCCCTGGGAGCTG
g.13565 G>T	3′UTR	F	TTTTGTCCCTGTCACCTTCC	60	425
		R	AGGGAGCATACTGGGTGATG
		E	TCCAAAGCCTAGACAACCCT
g.13757 C>A		F	TTTTGTCCCTGTCACCTTCC	60	425
		R	AGGGAGCATACTGGGTGATG
		E	ACCAAACTGAAAAAACAGGT
g.14047 C>T		F	GTCTGCCCATCACCCAGTAT	60	397
		R	TCACCACCATACCACCAATG
		E	GAATTGTTTTCTTTCTTTCA
g.14578 A>G		F	AGAGGGGGCTTGAAGAAGAG	60	469
		R	GCTTGCCTGTCCAGAAAAAG
		E	GTTGTTGGGTGTCCTGATAT
g.15001 A>G		F	TGACAGAATGGTCAGGGTCA	60	328
		R	TGACTCAGGGTAGGCAGCTT
		E	ACTTGGCCTCCAGAGTCTCC
g.15101 C>T		F	TGACAGAATGGTCAGGGTCA	60	328
		R	TGACTCAGGGTAGGCAGCTT
		E	GCATTCGCACTCATGAGGCA
g.15202 G>A		F	TGACAGAATGGTCAGGGTCA	60	328
		R	TGACTCAGGGTAGGCAGCTT
		E	CAGAACCAGAGAGCATGCTG
g.15491 G>A		F	GGGAAGCACAAGAAGCAAAG	60	320
		R	CAGTGGCAGAGCATTTGGTA
		E	CAACACAGTCTAGGGCAGAC
g.16524 C>T		F	GGTAGTTTTAATTCCACAGGTCAT	58	332
		R	TCCAAAATTATTCTTAGGTCACAA
		E	ATGAAAAGTTCAAGAAAGCT

Forward sequence.

Reverse sequence.

Extension sequence.

Linear discriminant and statistical analyses

For 45 extreme cattle, the significance of SNPs allele frequencies between the two groups (high and low) was analyzed using the chi-square test. Those SNP alleles for which there was a significant frequency difference were included in a linear discriminant analysis to develop a function classifying animals into the two groups (Tsuji et al., 2008). Only the influential fragments for discrimination (F >2.5) were included in the function.

The Hardy-Weinberg equilibrium (HWE) was determined for each locus by comparing expected frequency with observed genotype frequencies by using the chi-squares test for the Korean cattle population. The relationships between phenotypes (C18:1, MUFAs and MS) and individual SNPs for the 513 samples were analyzed using the following mixed analysis of covariance (ANCOVA) model based on SPSS v19.0 (SPSS Inc., Chicago, IL, USA).

where y is the phenotype, μ is the overall mean, S is the random effect of sire i with an assumption of independent and identical normal distribution, P is the fixed effect of calving place j (14 classes), G is the fixed effect of genotype k, β is the regression coefficient, age is a covariate for age in days at slaughter, e is the random residual assumed to have independent and identical normal distribution.

MDR analysis

The MDR method, a non-parametric and model-free technique, was initially implemented in case-control studies (Hahn et al., 2003). For its application to continuous data, the classification and regression tree (CART) algorithm has been developed and combined with the MDR method (Paolo, 2003). The expanded MDR method involves classification into two groups using a regression tree, i.e., high and low groups for phenotypes, and the group’s impurity can be defined as:

Where is the mean for the node (group) and y is the observation of the jth individual in a total of N individuals in the group (g = high or low). Each individual is assigned to a cell (e.g., if there are two SNPs with three genotypes per SNP, then there are 9 possible cells), and then each cell can be defined as a high or low group. The expected numbers of high and low groups are:

Respectively, where n is the number of total cells and I is an indicator function in which

Then the average squared error (ASE) is:

where is the sum of squared errors for the high group, is that for low group, and y is the phenotype of the j individual in the i cell.

In this study, the procedures for the expanded MDR method were as follows:

Step 1: Data were randomly divided into 10 equal parts (one testing set and nine training sets for cross-validation purpose).

Step 2: A set of n SNPs was selected from the pool of all SNPs.

Step 3: Based on the observed level of each n, steers were partitioned into classes, referred to as cells. If n = 2, then a set of two SNPs was selected, and, because a SNP had three genotypes, there were 32 = 9 possible cells. Phenotypic means were calculated within each cell.

Step 4: The impurity function in the CART algorithm used the variance impurity such that the group with the higher average value was labeled as high and the rest was labeled as low.

Step 5: The expanded MDR model with the smallest ASE was chosen among all two-factor combinations.

Step 6: For the evaluation of the predictive ability of the model, the predicted ASE (P_ASE) was estimated using a 10-fold cross-validation method. These six steps were repeated for each possible combination of given n. The model with the minimum predicted ASE was selected as the best model, but for the selected model, statistical significance was not determined by the predicted ASE. Therefore, permutation tests were performed to determine the empirical significance thresholds by applying the same MDR method (Good, 1994). Before the MDR implementation, phenotypes were adjusted for contemporary and, sire/steer’s age effects by using residuals obtained after fitting the general linear model without SNP effects.

RESULTS AND DISCUSSION

The SCD is an endoplasmic reticulum (ER) enzyme that catalyzes the biosynthesis of MUFAs from saturated fatty acids that are either synthesized de novo or derived from the diet (Smith et al., 2006). Variations in SCD enzyme activity in mammals are likely to affect a variety of key physiological variables, including cellular differentiation, insulin sensitivity, metabolic rates, adiposity, atherosclerosis, cancer, and obesity (Paton and Ntambi, 2009). Two types of SCD gene isoforms have been characterized in cattle, namely SCD1, expressed mainly in adipose tissue, and SCD5, expressed mainly in the brain (Lengi and Corl, 2007). The expression of SCD1 in bovine adipose tissue is regulated by numerous factors. The SCD1 mRNA level in muscle tissue has been reported to increase after weaning until 12 mon of age (Lee et al., 2005) but also during the late fattening stages (Kwon et al., 2009), and its development was dependent on the breed and diet (Chung et al., 2007). Breed differences in the SCD1 mRNA level have been found Japanese Black and Holstein steers (Taniguchi et al., 2004). Based on these reports, we collected 15 SNPs suggested by dbSNP of NCBI to identify those SNPs related to fatty acids in the SCD1 gene (Oh et al., 2011). These SNPs were located at 3 SNPs in the intron region, 3 SNPs in the exon 5 region and 9 SNPs in the 3′UTR region. By SBE method for verification of 15 SNPs in 45 extreme steers for the high and low MUFAs group of Korean native cattle, 15 polymorphic SNPs were identified (Table 1). And they were excluded in the following genetic association analysis. The other loci had their minor allele frequency ranging from 0.233 to 0.500, and their genotypes did not deviate from the Hardy-Weinberg equilibrium (p>0.05, Table 1).

Table 1

Genotype and frequency of 15 polymorphic SNPs of the SCD gene in Korean native cattle

SNP	Region	Genotype (No. of animals)				H1	MAF2	HWE3
		Frequency
g.6850+77 A>G	Intron	AA (2)	GA (17)	GG (26)	N4 (45)	0.358	0.244	0.708
		0.044	0.378	0.578	1
g.8243−72 A>G		AA (3)	GA (22)	GG (20)	N (45)	0.429	0.233	0.345
		0.067	0.489	0.444	1
g.8646+128 A>G		AA (9)	GA (21)	GG (15)	N (45)	0.491	0.311	0.738
		0.200	0.467	0.333	1
g.10153 A>G	Exon 5	AA (10)	GA (19)	GG (16)	N (45)	0.491	0.433	0.345
		0.222	0.422	0.356	1
g.10213 T>C		CC (13)	CT (22)	TT (10)	N (45)	0.498	0.448	0.905
		0.289	0.489	0.222	1
g.10329 C>T		CC (10)	CT (20)	TT (15)	N (45)	0.494	0.467	0.502
		0.222	0.444	0.333	1
g.13565 G>T	3′UTR	GG (8)	GT (21)	TT (16)	N (45)	0.484	0.444	0.808
		0.178	0.467	0.356	1
g.13757 C>A		AA (12)	GA (22)	GG (11)	N (45)	0.499	0.489	0.884
		0.267	0.489	0.244	1
g.14047 C>T		CC (10)	CT (23)	TT (12)	N (45)	0.499	0.478	0.871
		0.222	0.511	0.267	1
g.14578 A>G		AA (11)	GA (20)	GG (14)	N (45)	0.498	0.467	0.472
		0.244	0.444	0.311	1
g.15001 A>G		AA (10)	GA (20)	GG (15)	N (45)	0.494	0.444	0.502
		0.222	0.444	0.333	1
g.15101 C>T		CC (10)	CT (23)	TT (12)	N (45)	0.499	0.478	0.871
		0.222	0.511	0.267	1
g.15202 G>A		AA (11)	GA (23)	GG (11)	N (45)	0.500	0.500	0.881
		0.244	0.511	0.244	1
g.15491 G>A		AA (13)	GA (23)	GG (9)	N (45)	0.496	0.456	0.839
		0.289	0.511	0.200	1
g.16524 C>T		CC (8)	CT (26)	TT (11)	N (45)	0.498	0.467	0.280
		0.178	0.578	0.244	1

Heterozygosity.

Minor allele frequency.

Hardy–Weinberg equilibrium.

Total number.

Table 2 shows the efficiency of discrimination for the 15 polymorphic SNPs in the high and low MUFAs groups (n = 45). For the classification of 45 extreme animals into these two groups, the discriminant coefficients were calculated for the 15 polymorphic SNPs by using information on SNP genotypes the high- and low-MUFAs groups. The discriminant coefficients for 7 of these 15 SNPs (g.6850+77 A>G, g.8646+128 A>G, g.10153 A>G, g.10213 T>C, g.10329 C>T, g. 14047 C>T and g.14578 A>G) were −1.786, −0.500, −1.267, −0.424, −0.954, −0.966, and −0.944, respectively. These 7 SNPs showed significant relationships with MUFAs (p<0.05, F>2.5). This is equivalent to a discriminative error of 2.1% and a high correlation ratio of 0.91. These results indicate that the discriminant function successfully separated extreme cattle (Tsuji et al., 2008). Therefore, these 7 SNPs were selected for large scale genotyping in the Korean native cattle (n = 513).

Table 2

Efficiency of discrimination for 15 ploymorphic SNPs in high and low MUFAs groups of Korean native cattle (n = 45)

SNP	High MUFAs group			Low MUFAs group			Discriminant coefficient	F-value	p-value
	observed genotypeNo. of animals (Frequency)			observed genotypeNo. of animals (Frequency)
g.6850+77 A>G	AA	GA	GG	AA	GA	GG	−1.786	5.041	0.011
	2(1.00)	12(0.71)	8(0.31)	0(0.00)	5(0.29)	18(0.69)
g.8243−72 A>G	AA	GA	GG	AA	GA	GG	−0.305	0.241	0.787
	2(0.67)	11(0.5)	9(0.45)	1(0.33)	11(0.5)	11(0.55)
g.8646+128 A>G	AA	GA	GG	AA	GA	GG	−0.500	3.041	0.058
	4(0.44)	14(0.67)	4(0.27)	5(0.56)	7(0.33)	11(0.73)
g.10153 A>G	AA	GA	GG	AA	GA	GG	−1.267	5.358	0.008
	7(0.70)	12(0.63)	3(0.19)	3(0.30)	7(0.37)	13(0.81)
g.10213 T>C	TT	CT	CC	TT	CT	TT	−0.424	3.907	0.028
	4(0.40)	15(0.68)	3(0.23)	6(0.60)	7(0.32)	10(0.77)
g.10329 C>T	CC	CT	TT	CC	CT	TT	−0.954	4.255	0.021
	6(0.60)	13(0.65)	3(0.20)	4(0.40)	7(0.35)	12(0.80)
g.13565 G>T	GG	GT	TT	GG	GT	TT	0.185	0.248	0.781
	3(0.38)	11(0.52)	8(0.50)	5(0.62)	10(0.48)	8(0.50)
g.13757 C>A	CC	CA	AA	CC	CA	AA	−0.811	1.919	0.159
	8(0.73)	10(0.46)	4(0.33)	3(0.27)	12(0.54)	8(0.67)
g.14047 C>T	CC	CT	TT	CC	CT	TT	−0.966	3.739	0.032
	6(0.60)	14(0.61)	2(0.17)	4(0.40)	9(0.39)	10(0.83)
g.14578 A>G	AA	GA	GG	AA	GA	GG	−0.944	3.337	0.045
	7(0.64)	12(0.60)	3(0.21)	4(0.36)	8(0.40)	11(0.79)
g.15001 A>G	AA	GA	GG	AA	GA	GG	−0.069	0.021	0.979
	5(0.50)	10(0.50)	7(0.47)	5(0.50)	10(0.50)	8(0.53)
g.15101 C>T	CC	CT	TT	CC	CT	TT	0.179	0.199	0.821
	4(0.40)	12(0.52)	6(0.50)	6(0.60)	11(0.48)	6(0.50)
g.15202 G>A	GG	GA	AA	GG	GA	AA	0.000	0.095	0.909
	5(0.46)	12(0.52)	5(0.46)	6(0.54)	11(0.48)	6(0.54)
g.15491 G>A	GG	GA	AA	GG	GA	AA	0.185	0.098	0.907
	4(0.44)	11(0.48)	7(0.54)	5(0.56)	12(0.52)	6(0.46)
g.16524 C>T	CC	CT	TT	CC	CT	TT	−0.302	0.343	0.712
	5(0.63)	12(0.46)	5(0.46)	3(0.37)	14(0.54)	6(0.54)

The MDR method, a novel technique, was applied to examine the interactions between the 7 polymorphic SNP genotypes in complex traits (C18:1, MUFAs and MS). These were applied to the MDR method, and superior two- and three-way SNP combinations were used for each trait. Here p<0.05 was used to determine statistical significance for the critical value in the t-test and permutation test. For the investigation of two- and three-way SNP gene interactions, the top 5 combinations were listed for each trait (Tables 3 and 6). The ranking was determined based on the training-balanced accuracy of results for each trait.

Table 3

Significant SNP combinations of top 5 two-way gene interactions from 3 traits based on MDR method in Korean native cattle

Rank	SNP4	C18:11		SNP	MUFAs2		SNP	MS3
		Training-balanced accuracy	Testing-balanced accuracy		Training-balanced accuracy	Testing-balanced accuracy		Training-balanced accuracy	Testing-balanced accuracy
1	g.10213 T>C, g.14047 C>T	0.6016	0.5742	g.10329 C>T, g.14047 C>T	0.6140	0.5585		g.10153 A>G, g.10213 T>C	0.5678	0.4839
2	g.14047 C>T, g.14578 A>G	0.5997	0.5625	g.14047 C>T, g.14578 A>G	0.6131	0.5644	g.6850+77 A>G, g.10329 C>T	0.5642	0.4664
3	g.10329 C>T, g.14047 C>T	0.5386	0.5640	g.6850+77 A>G, g.14047 C>T	0.6109	0.5387	g.6850+77 A>G, g.14578 A>G	0.5634	0.4675
4	g.10153 A>G, g.14047 C>T	0.5968	0.5683	g.10213 T>C, g.14047 C>T	0.6056	0.5821	g.6850+77 A>G, g.14047 C>T	0.5608	0.4827
5	g.6850+77A>G, g.14047 C>T	0.5961	0.5464	g.8646+128 A>G, g.14047 C>T	0.6013	0.5762	g.8646+128 A>G, g.10213 T>C	0.5575	0.4927

Oleic acid.

Monounsaturated fatty acids.

Marbling score.

SNP combination.

Table 6

Significant SNP combinations of top 5 three-way gene interactions from 3 traits based on the MDR method in Korean native cattle

Rank	SNP4	C18:11		SNP	MUFAs2		SNP	MS3
		Training-balanced accuracy	Testing-balanced accuracy		Training-balanced accuracy	Testing-balanced accuracy		Training-balanced accuracy	Testing-balanced accuracy
1	g.6850+77 A>G, g.10213 T>C, g.14047 C>T	0.6376	0.5851	g.6850+77 A>G, g.10213 T>C, g.14047 C>T	0.6470	0.5789	g.6850+77 A>G, g.10213 T>C, g.14047 C>T	0.6185	0.5900
2	g.8646+128 A>G, g.10153 A>G, g.14047 C>T	0.6351	0.5860	g.8646+128 A>G, g.10153 A>G, g.14047 C>T	0.6432	0.5758	g.6850+77 A>G, g.8646+128 A>G, g.10213 T>C	0.6092	0.5722
3	g.6850+77 A>G, g.10153 A>G, g.14047 C>T	0.6314	0.5927	g.6850+77 A>G, g.10329 C>T, g.14047 C>T	0.6401	0.5619	g.6850+77 A>G, g.10213 T>C, g.13757 C>A	0.6072	0.5608
4	g.6850+77 A>G, g.10329 C>T, g.14047 C>T	0.6287	0.5923	g.8646+128 A>G, g.10329 C>T, g.14047 C>T	0.6339	0.6117	g.6850+77 A>G, g.10213 T>C, g.14578 A>G	0.6057	0.5737
5	g.6850+77 A>G, g.14047 C>T, g.14578 A>G	0.6233	0.5608	g.13757 C>A, g.14047 C>T, g.14578 A>G	0.6310	0.6263	g.6850+77 A>G, g.10153 A>G, g.10329 C>T	0.5979	0.5586

Oleic acid.

Monounsaturated fatty acids.

Marbling score.

SNP combination.

Table 3 shows that the g.10213 T>C and g.14047 C>T combination had the highest score for both training-balanced accuracy (60.16%) and testing-balanced accuracy (57.42%) in the two-way interaction model of C18:1. The g.10329 C>T and g.14047 C>T combination represented the best two-way interaction model of MUFAs. In case of MS, the g.10153 A>G and g.10213 T>C combination represented the highest training-balanced accuracy (56.78%).

The best interaction models of C18:1, MUFAs and MS were not the same, and therefore 4 superior SNP combinations (g.10213 T>C and g.14047 C>T; g.14047 C>T and g.14578 A>G; g.10329 C>T and g.14047 C>T; and g.6850+77 A>G and g.14047 C>T) were calculated and compared (Table 4). Then, the g.6850+77 A>G and g.14047 C>T SNP interaction was selected as the statistically best combination (C18:1, MUFAs and MS permutation t-values were 7.091, 7.017 and 3.293 respectively) of two-way interactions.

Table 4

The selected of 4 combinations based on two-way gene interactions for 3 traits in Korean native cattle

SNP1	Traits2	Genotype	N3	LSmean±SE	t-value	Permutation p4
g.10213 T>C, g.14047 C>T	C18:1	CCTT, CCTC, CTTC, CTCC	315	44.92±0.14	6.921***	0.000
		Others	198	43.32±0.18
	MUFA	CCTT, CCTC, CTTC, CTCC	315	54.17±0.16	6.691***	0.000
		Others	198	52.44±0.20
	MS	CCTT, CCTC, CTTT, CTTC, CTCC, TTTT, TTTC	397	5.56±0.10	2.879**	0.005
		Others	116	4.97±0.17
g.14047 C>T, g.14578 A>G	C18:1	AACC, AGCC, AGCT	294	44.97±0.15	6.941***	0.000
		Others	219	43.40±0.17
	MUFA	AACC, AACT, AGCC, AGCT	300	54.24±0.16	6.952***	0.000
		Others	213	52.47±0.19
	MS	AACC, AGCC, AGCT, AGTT, GGCC	326	5.59±0.11	2.560*	0.014
		Others	187	5.14±0.14
g.10329 C>T, g.14047 C>T	C18:1	CCCC, CCCT, CTCT, CTTT	318	44.91±0.14	6.888***	0.000
		Others	195	43.31±0.18
	MUFA	CCCC, CCCT, CTCT, CTTT	318	54.18±0.16	6.828***	0.000
		Others	195	52.41±0.20
	MS	CCCC, CCCT, CTCT, CTTT, TTCC, TTCT	373	5.60±0.10	3.290**	0.002
		Others	140	4.97±0.16
g.6850+77 A>G, g.14047 C>T	C18:1	CCAA, CCAG, CCGG, CTAG	214	45.24±0.18	7.091***	0.000
		Others	299	43.63±0.14
	MUFA	CCAA, CCAG, CTAA, CTAG	181	54.70±0.23	7.017***	0.000
		Others	332	52.86±0.15
	MS	CCAA, CCAG, CCGG, CTAG, TTAA, TTAG	214	5.76±0.13	3.293**	0.001
		Others	299	5.19±0.11

SNP combination.

C18:1 = Oleic acid, MUFA = Monounsaturated fatty acids, MS = Marbling score.

Number of animals.

Permutation p value.

p<0.05.

p<0.01.

p<0.001.

The single and combination effects of g.6850+77 A>G and g.14047 C>T were examined (Table 5). As shown in the table, the combination effects of g.6850+77 A>G and g.14047 C>T on C18:1, MUFAs and MS traits were better than their individual effects. In particular, g.6850+77 A>G and g.14047 C>T SNPs were not significant (with permutation p values of 0.060 and 0.192 respectively) in MS. However, the g.6850+77 A>G and g.14047 C>T combination was highly significant with p-value 0.001. In terms of superior genotypes CCAA, CCAG, CCGG and CTAG were obtained from C18:1 (45.24±0.18); CCAA, CCAG, CTAA and CTAG from MUFAs (54.70±0.23) and CCAA, CCAG, CCGG and CTAG from MS (5.76±0.13). Therefore, the CCAA, CCAG and CTAG genotype of the g.6850+77 A>G and g.14047 C>T combination had the greatest interaction effects on all 3 traits. Because all these combinations were significant, the best combination was determined based on the t-value.

Table 5

Single and combinations effect of 2 SNPs for 3 traits in Korean native cattle

Traits1	SNP	Genotype	N2	LSmean±SE	t-value	p3
C18:1	g.6850+77 A>G	AA, AG	220	45.03±0.20	5.380***	0.000
		GG	293	43.75±0.14
	g.14047 C>T	CC, CT	332	44.86±0.14	6.688***	0.000
		TT	181	43.28±0.19
	Combination	CCAA, CCAG, CCGG, CTAG	214	45.24±0.18	7.091***	0.000
		Others	299	43.63±0.14
MUFAs	g.6850+77 A>G	AA, AG	220	54.28±0.22	5.136***	0.000
		GG	293	52.92±0.15
	g.14047 C>T	CC, CT	332	54.09±0.16	6.268***	0.000
		TT	181	52.43±0.21
	Combination	CCAA, CCAG, CTAA, CTAG	181	54.70±0.23	7.017***	0.000
		Others	332	52.86±0.15
MS	g.6850+77 A>G	AA, AG	220	5.29±0.13	1.802	0.060
		GG	293	5.60±0.11
	g.14047 C>T	CC, CT	332	5.24±0.11	1.348	0.192
		TT	181	5.51±0.15
	Combination	CCAA, CCAG, CCGG, CTAG, TTAA, TTAG	214	5.76±0.13	3.293**	0.001
		Others	299	5.19±0.11

C18:1 = Oleic acid, MUFA = Monounsaturated fatty acids, MS = Marbling score.

Number of animals.

Permutation p value.

p<0.05.

p<0.01.

p<0.001.

For the three-way interaction model, the g.6850+77 A>G, g.10213 T>C and g.14047 C>T combination showed the highest training-balanced accuracy (63.76%, 64.70% and 61.85% respectively) (Table 6). Noteworthy is that g.6850+77 A>G and g.14047 C>T were also observed in the two-way interaction model. As shown in Table 7, the individual and interaction effects of g.6850+77 A>G, g.10213 T>C and g.14047 C>T were examined and compared for C18:1, MUFAs and MS. The superior genotype groups were 45.45±0.13 for C18:1, 54.73±0.22 for MUFAs and 5.95±0.12 for MS and significant (all permutation p-values were 0.000) (Table 7). By contrast, the average for the ht1*ht2 haplotype group (g.10153 A>G, g.10213 T>C and g.10329 C>T) was 44.92 for C18:1, as described in Oh et al. (2011); However, in this study, the average for genotype of g.6850+77 A>G (AA, AG), g.10213 T>C (TT, TC) and g.14047 C>T (CC, CT) SNP combinations was 45.45, better than 0.53 for the ht1*ht2 haplotype group. The combinations effect of SNPs was greater than the single effect for C18:1, MUFAs and MS (Table 7). As in the case of the two-way interaction effect, the individual effect of g.6850+77 A>G and g.14047 C>T on the MS were not significant. In addition, g.6850+77 A>G, g.10213 T>C and g.14047 C>T had a significant three-way interaction effect, and AATTCC, AGTTCC, GGTCCC, AGTCCT, GGCCCT and AGCCTT belonged to the superior group for all 3 traits.

Table 7

Single and combinations effect of 3 SNPs for 3 traits in Korean native cattle

Traits1	SNP	Genotype	N2	LSmean±SE	t-value	p3
C18:1	g.6850+77 A>G	AA, AG	220	45.03±0.20	5.380***	0.000
		GG	293	43.75±0.14
	g.10213 T>C	TT	112	44.91±0.27	2.743**	0.013
		TC, CC	401	44.13±0.13
	g.14047 C>T	CC, CT	332	44.86±0.14	6.688***	0.000
		TT	181	43.28±0.19
	Combination	AATTCC, AGTTCC, AATCCC, AGTCCC, GGTCCC, AGTCCT, GGCCCT, AGCCTT	189	45.45±0.13	7.484***	0.000
		Others	324	43.63±0.20
MUFAs	g.6850+77 A>G	AA, AG	220	54.28±0.22	5.136***	0.000
		GG	293	52.92±0.15
	g.10213 T>C	TT, TC	393	53.78±0.15	3.796***	0.003
		CC	120	52.62±0.27
	g.14047 C>T	CC, CT	332	54.09±0.16	6.268***	0.000
		TT	181	52.43±0.21
	Combination	AATTCC, AGTTCC, AATCCC, AGTCCC, GGTCCC, GGTTCT, AATCCT, AGTCCT, GGCCCT, AGCCTT	204	54.73±0.22	7.768***	0.000
		Others	309	52.70±0.15
MS	g.6850+77 A>G	AA, AG	220	5.60±0.13	1.802	0.060
		GG	293	5.29±0.11
	g.10213 T>C	TT, TC	393	5.56±0.10	2.828**	0.006
		CC	120	4.99±0.17
	g.14047 C>T	CC, CT	332	5.51±0.11	1.348	0.192
		TT	181	5.27±0.15
	Combination	AATTCC, AGTTCC, GGTCCC, AATTCT, GGTTCT, AGTCCT, AACCCT, AGCCCT, GGCCCT, AATTTT, GGTTTT, AATCTT, GGTCTT, AACCTT, AGCCTT	237	5.95±0.12	5.884***	0.000
		Others	276	4.97±0.11

C18:1 = Oleic acid, MUFA = Monounsaturated fatty acids, MS = Marbling score.

Number of animals.

Permutation p value.

p<0.05.

p<0.01.

p<0.001.

These results indicate that g.6850+77 A>G (AA, AG), g.10213 T>C (TT, TC) and g.14047 C>T (CC, CT) provided the best SNP combination for C18:1, MUFAs and MS in Korea native cattle. The MDR method is a fast and reliable technique for multiple SNP analyses of potential applications in marker-assisted selection. In addition, this technique, a non-parametric statistical method for detecting gene interactions, is useful for small samples and performs an exhaustive search of all n-locus models, collapsing multi-locus genotypes into high- and low-risk classes (Ritchie et al., 2001). This suggests that, after further development, the MDR method may be used for the genetic assessment of quantitative traits. Finally, these genotypes of individual SNPs and their combinations may be important genetic markers for improving beef quality, which is one of the most important goals in the Korean cattle industry.