Application of linkage disequilibrium mapping methods to detect QTL for carcass quality on chromosome 6 using a high density SNP map in Hanwoo.
In the last decade, many advanced technologies have been developed in genomics such as next-generation sequencing and high-throughput genotyping platforms, which provide high-density SNP arrays as a state-of-the-art tool for genetic and genomic analyses in domestic animals (Fan et al., 2010). Due to the developments of high density SNP maps and high-throughput genotyping techniques, marker density is no longer a limiting factor in QTL finemapping studies. However, genetic analysis using large amount of SNPs, e.g. whole genome association (WGA) study, requires efficient statistical methods (Druet et al., 2008).
Linkage disequilibrium (LD) mapping methods have been practiced for fine QTL mapping studies in livestock species. The efficiency of the LD mapping depends on sample size, SNP map density, magnitude of QTL effect, effective population size, population history, and statistical models (Zhao et al., 2007). To evaluate mapping power and mapping efficiency between LD mapping models , Grapes et al. (2004), Grapes et al. (2006), and Zhao et al. (2007) compared several types of fine LD mapping models with single marker regression, regression based on haplotypes using more than two markers, and identity-by-descent (IBD) mapping. They made a conclusion that the single marker regression analysis had power to detect QTL and mapping precision not worse than the haplotype-based regression or the IBD mapping method. However, Hayes et al. (2007) reported that, when haplotypes with more than four markers were used, the LD mapping method using haplotypes had greater mapping accuracy. Calus et al. (2008) also reported that the IBD approach gave greater accuracy for breeding value estimation than either the single maker regression or the haplotype-based regression method.
Korean consumers prefer a beef of the Korean native cattle, Hanwoo, because the breed has good genetic characteristics in meat palatability and chewiness (Cho and Ko, 1998). In a previous study, we reported three QTL affecting backfat thickness, average daily gain, and final weight on BTA6 in a Hanwoo population by applying a linkage mapping method (Lee et al., 2010b). Also, Lee et al. (2008; 2010c) reported that the genotypes of the 12273 165 SNP within the CCDC158 gene are significantly associated with carcass cold weight.
The aim of this study was to further refine the QTL region on BTA6, and to search for candidate genes affecting carcass qualities in Hanwoo by using a high density SNP map and different LD mapping methods.
MATERIALS AND METHODS
Animals, phenotypes and molecular data
The steers (N = 473) with phenotype and molecular data were chosen among the progeny of candidate bulls for progeny testing in the Hanwoo Improvement Center of National Agriculture Cooperative Federation in Seosan, Chungnam province, Korea. The data set comprised 45 sires and their 473 steers that were born between spring of 2005 and fall of 2007. The number of steers for each of the 45 paternal sire families ranged from 6 to 13 with the average eight steers.
The steers were weaned at 5 or 6 months of age, and each group of 10 steers were raised in a pen. The feeding program was divided into early, middle, and late stage, each with six months of interval. In the early and middle stages, the steers were fed with concentrates with the amount of 1.8% of the body weight and ad labium in the late stage. The concentrates were composed of 15%, 13%, and 11% of crude protein, and 71%, 72%, and 73% totally digestible nutrients (TDN), in the respective feeding stages. Roughages with 4.5% crude protein and 37.5 TDN were offered ad labium with other additives such as vitamin and minerals. All steers were slaughtered at an approximately 24 months of age. Traits measured were carcass weight after slaughter (CWT), backfat thickness (BFT), longissimus dorsi muscle area (LMA), and marbling score (Marb). The Marb score was numbered as 1 through 9 according to the Korean Beef Marbling Standard (1 = trace, 9 = very abundant). Details about measurement of the traits and management practices were described in Lee et al. (2010a). Table 1 shows the summary statistics for the observed carcass quality traits.
A dense marker map covering BTA6 was constructed with the 2,535 SNPs that were embedded in the Illumina Bovine SNP50K BeadChip (Matukumalli et al., 2009). Details on the SNP genotyping experiments were described in Lee et al. (2010a). Among the 2,535 SNPs, those were removed based on the following criteria: i) monomorphic SNPs or when the minor allele frequency was smaller than 0.05, ii) proportion of individuals with genotype completeness was smaller than 90%, iii) markers with significant departure from H-W equilibrium (p<0.001).
To test the association between a SNP and QTL for carcass quality, the following four models were used.
The fixed model fitting a single marker (Fix_Sig),
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII.]
The random model fitting a single marker (Ran_Sig),
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII.]
The random model fitting haplotypes using two adjacent markers (Ran_Hap),
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII.]
The random model fitting haplotypes under the hidden states (Ran_HS),
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII.]
Where [y.sub.i] is the phenotypic record of animal i, m is the average phenotypic performance, sirei is the random sire effect for animal i, [c.sub.ji] is the value of the [j.sup.th] covariate or fixed effect for the animal i, [[beta].sub.j] is an estimate of the [j.sub.th] fixed effect or covariate, [g.sub.ij]_fixed is the fixed effect of SNP genotype values (e.g. -1 for BB, 0 for BA, and 1 for AA) for animal i, [g.sub.ij]_random is the random effect of the SNP alleles (A vs. B) for animal i, [haplotype.sub.ijk] is a random effect for a paternal (k = 1) or maternal (k = 2) haplotype of the two adjacent loci of animal i, [HS.sub.ij] is a random effect for the K defined hidden states for animal i, and ei is a random residual for animal i. Construction of haplotype and reconstruction of haplotype hidden states was based on the approach of Druet and Georges (2010). Haplotypes were inferred from the fastPHASE probability model, which is a family rule (Mendelian segregation and linkage) based algorithm. Observed haplotypes were modeled as mosaics of K hidden states (or ancestral haplotypes), with K held constant throughout the genome (Druet et al., 2008).
Two fixed effects were fitted in the models across all of the fixed and random models; year and season of birth (5 levels) and region where the steers were born (40 levels), and slaughter age was fitted as a covariate. For the fixed model in which a single maker was fitted, phenotypes were preadjusted using SAS GLM procedure of SAS v9.1 (SAS Inst., Inc., Cary, NC), and the residuals of the phenotypes were used in the Fix_Sig model. For the random models, effect of sire was fitted as a random term. The models were tested at each SNP map position on BTA6. The lack of fit (LOF) tests was applied to obtain test statistics in the Fix_Sig model (Lee et al., 2010a),
[LOF = [RSS.sub.Ho] - [RSS.sub.Ho]/[RSS.sub.Ho]
where the RSS is the residual sum of squares, and Ho and Ha indicate the models without and with fitting the SNP, respectively (Neter et al., 1990). Threshold values for the LOF tests were obtained from the F distribution.
For the three random models, likelihood ratio tests (LRT) were applied by comparing the maximum likelihoods between the full model with the testing SNP or haplotypes and the reduced model without the QTL effect;
LRT = -2([LogLikehood.sub.reduced_ model] - [LogLikehoodsub.full_ model])
The LRT test statistics approximately followed chisquared distributions with 1 degree of freedom (Olsen et al., 2004). For significance threshold, 0.1% comparison-wise p value was applied for QTL detection.
RESULTS AND DISCUSSION
Among the 2,535 SNPs covering BTA6 in the Illumina bovine 50K SNP chip, 1,855 SNPs (73.2%) were available for the association tests. This set of markers covered 122.54 megabase (Mb) of the BTA6 with the 66.04 ([+ or -]65.04) kb average distance between adjacent markers. A total of 12 QTL were detected, for which four, two, five, and one QTL were detected for CWT, BFT, LMA, and MS, respectively (Table 2). Mapped QTL can be accessed in the Ensembl (Hubbard et al., 2009). Three QTL for CWT were detected in the proximal region of BTA6 (30 to 45 Mb), on which several genes including FAM190A and ABCG2 were located (Figure 1). The former encodes actin-binding LIM protein family, member 2. A QTL for BFT was detected in the proximal region of the BTA6 (Table 2; Figure 1). In the region, several genes including SYNPO2 and CAMK2D were located, the latter of which encodes calcium/calmodulin-dependent protein kinase type II delta chain. Four QTL for LMA were detected in the distal region of the test chromosome (100 to 120 Mb), in which one gene, SCD5, that encodes stearyl-CoA desaturase 5, was located. One QTL was detected for Marb in the proximal region of the BTA6 (43Mb), around which several genes such as SLIT2 and STIM2 were located (Figure 1).
[FIGURE 1 OMITTED]
The four different LD mapping methods enabled the detection of the four, one, three, and four QTL in the Fix_Sig, Ran_Sig, Ran_Hap, and Ran_HS models, respectively (Table 2 and Figure 1). The QTL was detected only in one of the four mapping models, i.e. no QTL was detected in more than two mapping models (Table 2). This may be partly due to the small sample size (N = 473), such that significant QTL in one mapping model can be undetected in another models. Another reason of the inconsistency of QTL detection between the models may be small magnitude of QTL effects. Previously, we scanned BTA6 to find QTL for the carcass quality traits using an interval mapping approach in a Hanwoo population, but only one QTL was detected for BFT at the 5% chromosome-wise level (Lee et al., 2010b). Therefore, our results suggest that the application of multiple LD mapping approaches is beneficial in QTL detection, when QTL effect is of no great magnitude in an experiment with a small sample size. However, care must be taken in determining QTL definition, due to the possibility of false positive QTL.
Some QTL were located in close distances, e.g. the two SNPs (BTA-75768-no-rs and BFGL-NGS-114855) for CWT at 35.3 Mb and 36.1 Mb, and the two SNPs (ARSBFGL-NGS-32948 and ARS-BFGL-NGS-103739) for LMA at 113.6 Mb, respectively (Table 2 and Figure 1). Some of the SNPs may be not causal mutation, i.e. the SNPs were so close to the causal SNP as to form high LDs with the causal mutation for the trait.
Threshold values for the statistical analysis were based on the point-wise level, i.e. one test basis, and chromosomewise threshold values were not applied in this study. Because evidences of the QTL for the carcass traits on BTA6 were not strong, i.e. the nominal P (-[log.sub.10]P) values were greater than 0.00001 (<4.0) (Table 2 and Figure 1), no QTL would be detected when multiple tests were taken into account in this study.
Several QTL studies have been conducted to find any association with carcass quality traits in Hanwoo cattle (Kim et al., 2004; Cheong et al., 2008; Chung et al., 2008). Most reports were based on candidate gene studies, in which genes with physiological functions that were related with the trait of interest were pre-selected and tested. However, there was no report on candidate genes that were located on BTA6 (Cheong et al., 2007; Shin et al., 2007; Cho et al., 2008). Previously, we scanned BTA6 using an interval mapping method in Hanwoo paternal halfsib families (Lee et al., 2010b). Lee et al. (2010b) reported that one QTL for BFT were detected at the 5% chromosomewise level at 0 cM, which was close to the marker, IL90 (0 cM). Another QTL for LMA was detected at the 10% chromosome-wise level at 67 cM, close to BM4621. In the QTL regions similar to Lee et al. (2010b), we detected QTL for BFT (13.2 Mb) and LMA (55.3 Mb) (Table 2 and Figure 1). In our previous study (Lee et al., 2010a), no significant SNP was found for the carcass quality traits on BTA6 in a Hanwoo population. This may be partly due to the small sample size used in the previous study (N = 289). However, we applied several LD mapping methods with a larger sample (N = 473), and detected twelve QTL for the four carcass quality traits (p<0.001).
In conclusion, our results suggest that the use of multiple LD mapping approaches are beneficial in increasing power to detect QTL given a limited sample size and magnitude of QTL effect. The detected QTL need to be verified in other Hanwoo populations for commercial application via marker-assisted selection.
This research was supported by the Technology Development Program for Agriculture and Forestry, Ministry of Agriculture, Forestry and Fisheries, Republic of Korea, 2010.
Received January 26, 2011; Accepted February 30, 2011
Calus, M. P., T. H. Meuwissen, A. P. de Roos and R. F. Veerkamp. 2008. Accuracy of genomic selection using different methods to define haplotypes. Genetics 178(1):553-561.
Cheong, H. S., D. H. Yoon, L. H. Kim, B. L. Park, H. W. Lee, C. S. Han, E. M. Kim, H. Cho, E. R. Chung, I. Cheong and H. D. Shin. 2007. Titin-cap (TCAP) polymorphisms associated with marbling score of beef. Meat Sci. 77:257-263.
Cheong, H. S., D. H. Yoon, B. L. Park, L. H. Kim, J. S. Bae, S. Namgoong, H. W. Lee, C. S. Han, J. O. Kim, I. C. Cheong and H. D. Shin. 2008. A single nucleotide polymorphism in CAPN1 associated with marbling score in Korean cattle. BMC Genet. 9:33-39.
Cho, B. D. and Y. D. Ko. 1998. Hanwoo meat. In: Domestic animal industry in Korea (Ed. J. K. Jung). World Association of Animal Science in Korea. pp. 3-25.
Cho, S., T. S. Park, D. H. Yoon, H. S. Cheong, S. Namgoong, B. L. Park, H. W. Lee, C. S. Han, E. M. Kim, I. C. Cheong, H. Kim and H. D. Shin. 2008. Identification of genetic polymorphisms in FABP3 and FABP4 and putative association with back fat thickness in Korean native cattle. BMB Rep. 41(1):29-34.
Chung, E. R., S. C. Shin, K. H. Shin and K. Y. Chung. 2008. SNP discovery in the leptin promoter gene and association with meat quality and carcass traits in Korean Cattle. Asian-Aust. J. Anim. Sci. 21(12):1689-1695.
Druet, T., S. Fritz, M. Boussaha, S. Ben-Jemaa, F. Guillaume, D. Derbala, D. Zelenika, D. Lechner, C. Charon, D. Boichard, I. G. Gut, A. Eggen and M. Gautier. 2008. Fine mapping of quantitative trait loci affecting female fertility in dairy cattle on BTA03 using a dense single-nucleotide polymorphism map. Genetics 178(4): 2227-2235.
Druet, T. and M. Georges. 2010. A hidden markov model combining linkage and linkage disequilibrium information for haplotype reconstruction and quantitative trait locus fine mapping. Genetics 184(3):789-798.
Fan, B., Z. Q. Du, D. M. Gorbach and M. F. Rothschild. 2010. Development and application of high-density SNP arrays in genomic studies of domestic animals. Asian-Aust. J. Anim. Sci. 23(7): 833-847.
Grapes, L., J. C. Dekkers, M. F. Rothschild and R. L. Fernando. 2004. Comparing linkage disequilibrium-based methods for fine mapping quantitative trait loci. Genetics 166(3):1561 1570.
Grapes, L., M. Z. Firat, J. C. Dekkers, M. F. Rothschild and R. L. Fernando. 2006. Optimal haplotype structure for linkage disequilibrium-based fine mapping of quantitative trait loci using identity by descent. Genetics 172(3):1955-1965.
Hayes, B. J., A. J. Chamberlain, H. McPartlan, I. Macleod, L. Sethuraman and M. E. Goddard. 2007. Accuracy of markerassisted selection with single markers and marker haplotypes in cattle. Genet. Res. 89(4):215-220.
Hubbard, T. J., B. L. Aken, S. Ayling, B. Ballester, K. Beal, E. Bragin, S. Brent, Y. Chen, P. Clapham, L. Clarke, G. Coates, S. Fairley, S. Fitzgerald, J. Fernandez-Banet, L. Gordon, S. Graf, S. Haider, M. Hammond, R. Holland, K. Howe, A. Jenkinson, N. Johnson, A. Kahari, D. Keefe, S. Keenan, R. Kinsella, F. Kokocinski, E. Kulesha, D. Lawson, I. Longden, K. Megy, P. Meidl, B. Overduin, A. Parker, B. Pritchard, D. Rios, M. Schuster, G. Slater, D. Smedley, W. Spooner, G. Spudich, S. Trevanion, A. Vilella, J. Vogel, S. White, S. Wilder, A. Zadissa, E. Birney, F. Cunningham, V. Curwen, R. Durbin, X. M. Fernandez-Suarez, J. Herrero, A. Kasprzyk, G. Proctor, J. Smith, S. Searle and P. Flicek. 2009. Ensembl 2009. Nucleic Acids Res. 37(Database issue): D690-697.
Kim, N. K., Y. W. Seo, G. H. Kim, J. H. Joh, O. H. Kim, E. R.
Chung and C. S. Lee. 2004. A previously unreported DraI polymorphism within the regulatory region of the bovine growth hormone gene and its association with growth traits in Korean Hanwoo cattle. Anim. Genet. 35:152-154.
Lee, Y. M., C. M. Han, Y. Li, J. J. Lee, L. H. Kim, J. H. Kim, D. I. Kim, S. S. Lee, B. L. Park, H. D. Shin, K. S. Kim, N. S. Kim and J. J. Kim. 2010a. A whole genome association study to detect single nucleotide polymorphisms for carcass traits in hanwoo populations. Asian-Aust. J. Anim. Sci. 23(4):417-424.
Lee, Y. M., Y. S. Lee, C. M. Han, J. H. Lee, J. S. Yeo and J. J. Kim. 2010b. Detection of quantitative trait loci for growth and carcass traits on BTA6 in a Hanwoo population. Asian-Aust. J. Anim. Sci. 23(3):287-291.
Lee, Y.-S., D. Y. Oh, J.-J. Kim, J.-H. Lee, H.-S. Park and J.-S. Yeo. 2010c. A single nucleotide polymorphism in LOC534614 as an unknown gene associated with body weight and cold carcass weight in Hanwoo (Korean cattle). Asian-Aust. J. Anim. Sci. 23(12): 1543-1551.
Lee, Y. -S., J. H. Lee, J. Y. Lee, J. J. Kim, H. S. Park and J. S. Yeo. 2008. Identification of candidate SNP (single nucleotide polymorphism) for growth and carcass traits related to QTL on chromosome 6 in Hanwoo (Korean Cattle). Asian-Aust. J. Anim. Sci. 21(12):1703-1709.
Matukumalli, L. K., C. T. Lawley, R. D. Schnabel, J. F. Taylor, M. F. Allan, M. P. Heaton, J. O'Connell, S. S. Moore, T. P. Smith, T. S. Sonstegard and C. P. Van Tassell. 2009. Development and characterization of a high density SNP genotyping assay for cattle. PLoS One, 4(4):e5350.
Neter, J., W. Wasserman and M. H. Kutner. 1990. Applied linear statistical models. 3rd ed. Irwin. Boston.
Olsen, H. G., S. Lien, M. Svendsen, H. Nilsen, A. Roseth, M. Aasland Opsal and T. H. Meuwissen. 2004. Fine mapping of milk production QTL on BTA6 by combined linkage and linkage disequilibrium analysis. J. Dairy Sci. 87(3):690-698.
Shin, S. C. and E. R. Chung. 2007. Association of SNP marker in the thyroglobulin gene with carcass and meat quality traits in Korean cattle. Asian-Aust. J. Anim. Sci. 20(2):172-177.
Zhao, H. H., R. L. Fernando and J. C. Dekkers. 2007. Power and precision of alternate methods for linkage disequilibrium mapping of quantitative trait loci. Genetics 175(4):1975-1986.
* Corresponding Author : Jong-Joo Kim. Tel: +82-53-810-3027, Fax: +82-53-810-4769, E-mail: firstname.lastname@example.org
Y. Li (a), J.-H. Lee (1,a), Y.-M. Lee and J.-J. Kim*
School of Biotechnology, Yeungnam University, Gyeongsan, Gyeongbuk, 712-749, Korea
(1) Gyeongbuk Livestock Research Institute, Yeungju, 750-780, Korea.
(a) Both contributed equally.
Table 1. Summary statistics of observations on carcass traits for the 427 Hanwoo steers Trait Average Std. Dev. (a) Minimum Carcass weight (kg) 356 40 158 Backfat thickness (cm) 10.2 0.5 0.3 Longissimus dorsi muscle 79.1 9.3 22 area ([cm.sup.2]) Marbling score (1-9) 3.4 1.8 1 Trait Maximum CV (b) Carcass weight (kg) 481 11.1 Backfat thickness (cm) 60 48.8 Longissimus dorsi muscle 109 11.8 area ([cm.sup.2]) Marbling score (1-9) 9 52.2 (a) Standard deviation. (b) Coefficisent of variation (%). Table 2. Identities, positions of the SNPs associated with carcass traits with statistical significance at the point-wise 0.001 level on BTA6 Trait Model (b) Position(Mb) -[log.sub.- SNP Marker (a) (c) 10]P (e) Carcass weight BTA-75768-no-rs Ran_HS 35.25 3.06 BFGL-NGS-114855 Ran_HS 36.09 3.15 BTB-00406718 Ran_Hap 40.84 3.01 Hapmap27076-BTC-037281 40.89 ARS-BFGL-NGS-41730 Fix_Sig 122.51 3.29 Backfat thickness BTB-00244322 Ran_Hap 13.16 3.79 Hapmap40587-BTA-76947 13.23 Hapmap27519-BTC-047206 Ran_Sig 54.55 3.03 Longissimus dorsi muscle area BTB-00843812 Fix_Sig 55.29 3.62 BTB-01096127 Fix_Sig 101.41 3.02 ARS-BFGL-NGS-69648 Ran_Hap 110.49 4.63 BTB-00280906 110.53 ARS-BFGL-NGS-32948 Ran_HS 113.58 3.02 ARS-BFGL-NGS-103739 Ran_HS 113.62 3.02 Marbling score Hapmap59861-rs29027897 Fix_Sig 42.94 3.12 Trait Gene close to Marker SNP Marker (a) Carcass weight BTA-75768-no-rs FAM190A (35.76Mb-36.10Mb) BFGL-NGS-114855 FAM190A (35.76Mb-36.10Mb) BTB-00406718 SLIT2 (41.19Mb-41.35Mb) Hapmap27076-BTC-037281 ARS-BFGL-NGS-41730 ABLIM2 (122.44Mb-122.54Mb) Backfat thickness BTB-00244322 KCC2D (12.89Mb-13.32Mb) Hapmap40587-BTA-76947 Hapmap27519-BTC-047206 PCDH7 (51.52Mb-51.53Mb) Longissimus dorsi muscle area BTB-00843812 LOC781512 (57.54Mb-57.79Mb) BTB-01096127 SCD5 (101.23Mb-101.41Mb) ARS-BFGL-NGS-69648 CLNK (109.94Mb-110.19Mb) BTB-00280906 ARS-BFGL-NGS-32948 BOD1L (113.44Mb-113.46Mb) ARS-BFGL-NGS-103739 BOD1L (113.44Mb-113.46Mb) Marbling score Hapmap59861-rs29027897 GPR125 (43.20Mb-43.29 Mb) (a,c) SNP marker annotations and their positions were based on the bovine reference genome (btau4.0). (b) Declared LD mapping model type: Fix_Sig = The fixed model fitting a single marker; Ran_Sig = The random model fitting a single marker; Ran_Hap = The random model fitting haplotypes using two adjacent markers; Ran_HS = The random model fitting haplotypes under the hidden states. (e) Negative logarithm of the comparison-wise p-value of the test statistic against the null hypothesis of no SNP effect at the SNP position.
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||quantitative trait loci, single-nucleotide polymorphism|
|Author:||Li, Y.; Lee, J.-H.; Lee, Y.-M.; Kim, J.-J.|
|Publication:||Asian - Australasian Journal of Animal Sciences|
|Date:||Apr 1, 2011|
|Previous Article:||Factor analysis of biometric traits of Kankrej cows to explain body conformation.|
|Next Article:||Single nucleotide polymorphisms of the GnRHR gene associated with reproductive traits of Japanese flounder (Paralichthys olivaceus).|