|
|
||||||||





* National Livestock Breeding Center, Nishigo, Fukushima 961-8511, Japan;
and
Ibaraki Branch of National Livestock Breeding Center, Fujigaya, Ibaraki 308-0112, Japan;
and
National Institute of Agrobiological Sciences, Tsukuba, Ibaraki 305-0901, Japan; and
and
Shirakawa Institute of Animal Genetics, Nishigo, Fukushima 961-8061, Japan
| Abstract |
|---|
|
|
|---|
Key Words: Growth Linkage Analysis Meat Quality Pigs Quantitative Trait Loci
| Introduction |
|---|
|
|
|---|
We generated an F2 resource population from a cross between a Meishan female and a Duroc male in order to map loci affecting economically important traits. Duroc boars are extensively used in terminal crosses as sires for the production of pork meat with improved meat quality in Japan because the Duroc breed is characterized by high meat quality and large muscle mass. Meishan, characterized by large litter size, differ from commercial breeds in many respects, including genetic characters. These two breeds exhibit large differences not only in growth performance, carcass composition, and meat quality, but also in reproductive traits. Therefore, we selected Meishan and Duroc breeds as parents for a resource population to detect QTL affecting growth, carcass and reproductive traits in pig and to obtain information and tools for breeding programs, such as marker-assisted introgression. We report porcine QTL results for 48 traits, including carcass traits and meat quality.
| Materials and Methods |
|---|
|
|
|---|
Coat Color.
Coat color was recorded from six angled photographs taken of each animal, and scored as follows: Class 1 score = gray coat color; Class 2 score = brown coat color; Class 3 score = black; Class a score = spotting; Class b score = no spotting. Although coat color is a qualitative trait, we treated it as a quantitative trait for analysis.
Data Collection.
Pigs were weighed at birth and at 21, and 30 d of age. From these pigs, 292 were weighed at 60 d of age. Testicular weights were recorded from F2 males that were castrated at 2 mo of age. Among the castrated males, 165 were raised until they obtained a body weight of 90 kg. Average daily gain was calculated between 30-kg and 90-kg body weights. These males were slaughtered at a live weight of approximately 90 kg. At slaughter, these males were recorded for body weight, size, and length, including the circumference of the chest and cannon bone, the height at the withers, and the chest depth and width. Carcasses were scalded and dehaired, and chilled overnight. Then carcass measurements were recorded for weight; lengths I, II, and III; thickness; and vertebra number. Carcass lengths I, II, and III refer to the lengths from the first cervical to the pubic bone, from the first rib to the pubic bone, and from the first rib to last lumbar vertebra, respectively. Depth of backfat over the midline was recorded at the first rib (shoulder), the thinnest depth at ribs (back), and the first lumbar vertebra (loin) with a ruler. Depth of belly over the midline was recorded at the last rib (forebelly), at the diaphragm (middle belly), and at the last lumbar vertebra (hind belly). One side of the carcass was split between the 4th and 5th ribs, last rib, and the first lumbar vertebra (at the middle); the longissimus muscle was traced on acetate paper; and the area was determined using computerized morphometric planimetry. The remaining side was then weighed and cut into the major wholesale cuts. The weight of each wholesale cut was recorded. Intramuscular fat (IMF) content was measured as described by Gerbens et al. (1999)
. After minced meat was dried by heating to 102°C for 24 h in a drying oven, moisture was calculated from weights taken before and after drying and was expressed as a percentage. Objective measurements of meat color were taken with a color meter, model ZE 2000 (Nippon Denshoku Industries Co., Tokyo, Japan); L values measure light reflectance, a values represent the degree of redness, and b values represent the degree of yellowness. Cooking loss and shear force were measured as described by Hovenier et al. (1992)
.
Table 1
summarizes the data of 48 traits, including ADG, vertebra number, backfat thickness (BFT), and IMF. Because the 165 F2 males were produced in three independent groups and in different parities, possible phenotypic differences derived from groups and parities were corrected for using the mixed-model least squares, maximum-likelihood computer program described by Harvey (1977)
.
|
Linkage Analysis
Linkage maps were constructed using CRI-MAP (Green et al., 1990
) for 18 autosomes and the sex chromosome. The sex-averaged map was used for the whole-genome scan of QTL. The information content was calculated using the method described by Knott et al. (1998)
. When a given location belonged to a pseudoautosomal region of a sex chromosome, analyses were performed using the same method as for autosomes. For testing QTL in sex-specific regions of the sex chromosome, the analysis model was modified as described by Knott et al. (1998)
. We assumed that the Y chromosome contained no loci. A test statistic for detecting QTL and the threshold was constructed as for autosomes.
A QTL analysis for each trait was performed using the method developed by Haley et al. (1994)
. The statistical model in the method is based on a linear regression of phenotypes on the probabilities of QTL genotypes at a given location. We assumed that the grandparental breeds were fixed for alternative alleles at a QTL. Two alleles at a putative QTL at a given location were denoted by Q and q. Probabilities of QTL genotypes, denoted by Prob(QQ), Prob(Qq), and Prob(qq), were calculated from the observed genotypes of markers linked to the QTL. The calculation was described by Haley et al. (1994)
. In the analyses of real data, sex difference was taken into account.
There are three possible genotypes, QQ, Qq, and qq, for a QTL at the given location on an autosome. Let the effects of genotypes QQ, Qq, and qq be denoted by a, d, and -a. We assume that phenotypic value of a trait is written for the ith individual in F2 as follows:
![]() |
where µ is the mean; si is the indicator of the sex of individual i, which is equal to 1 or -1 for male or female, respectively; h is a sex effect; cai is the coefficient for the additive component for individual i at the given location that is calculated from the probability of QTL genotypes, and equal to Prob(QQ) - Prob(qq); cdi is the coefficient for the dominance component for individual i at the given location, which is equal to Prob(Qq); and ei is the residual error. Model parameters, µ, h, a, and d, are estimated by a least squares method. That is, estimators of the parameters are obtained such that a sum of squares,
![]() |
is minimized, where n is the number of individuals of F2. Denoting least squares estimators of µ, h, a, and d, by the terms
,
, â, and
, the minimum sum of squares is obtained as
![]() |
Under the null model corresponding to no QTL, where a = d = 0 is assumed, the minimum sum of squares is denoted by S0. Significance detection of QTL is declared based on the ratio of S1 and S0. In this report, we used F-ratio, [(S0 - S1)/2]/[S1/(n - 2)], as a statistic for detecting QTL. Significance thresholds were obtained with 10,000 repeated permutation test cycles for each trait.
In chromosomal regions exceeding the suggestive and significance level on some QTL, P-values of the nearby markers were computed by SAS-GLM (SAS Inst. Inc., Cary, NC). When two parents have more than two marker alleles, an analysis of the individual marker genotypes will efficiently detect the associations of the different alleles in one or both breeds with specific traits. Phenotypic difference among genotypes of nearby markers was estimated with SAS-GLM. The number of marker genotypes was different depending on the genotypes of the parents. This analysis of variance supported results of QTL analysis and estimated within-breed allele effects rather than between-breed allele effects.
Linkage Map
We searched for informative microsatellites from the USDA-MARC linkage map (Rohrer et al., 1996
) and selected 180 microsatellites at approximately 20-cM intervals in the whole pig genome. The 180 microsatellites were genotyped across parents, 28 F1, and 865 F2 animals. A sex-averaged linkage map was constructed with multipoint analysis. The order of the markers was the same as in the USDA-MARC linkage map (Rohrer et al., 1996
). The average interval was 12.6 cM, and the total length was 2,263.6 cM, which was comparable to that of the USDA-MARC linkage map (Rohrer et al., 1996
). Two large gaps existed between SW274 and SW72 on SSC 3 (34.6 cM) and SW949 and SW980 on SSC X (40.7 cM). The average information content was 0.80.
As this F2 population consisted of only full-sib families, each genotype was checked for segregation distortion. To confirm the presence of a lethal allele, we examined whether the segregation ratio at each marker fit with mendelian ratios 1:2:1, Duroc homozygous, heterozygous, and Meishan homozygous, respectively, by chi-squared test.
| Results and Discussion |
|---|
|
|
|---|
Mapping of Coat Color Traits
The F2 animals exhibited coat colors of black (n = 373), brown (n = 383), and gray (n = 109). Several brown and gray pigs had black spotting (spotting, n = 289; no spotting, n = 576). The regions associated with coat color were mapped on SSC 1 (107.7 cM) at an F-ratio of 28.2 and SSC 6 (1 cM) at an F-ratio of 122.0. The spotting region was also mapped on SSC 6 (2 cM) with an F ratio of 111.7. The locus on SSC 1 was not previously reported. On SSC 6, the genotypes of marker S0035 at the coat color and the spotting loci were tightly linked to black and brown coat color and spotting. Because the melanocortin-1 receptor (MC1R) gene was mapped on the SSC 6 telomeric region (Mariani et al., 1996
) and the mutations were associated with pig coat color and spotting (Kijas et al., 1998
, 2001
), partial MC1R was amplified by PCR using genomic DNA from F1 and F2 pigs, and sequenced. We confirmed a previous report that Duroc and Meishan pigs are homozygous for MC1R*4 and heterozygous for MC1R*2/3 (Kijas et al., 1998
), respectively, suggesting that the MC1R gene is responsible for coat color and spotting on SSC 6.
QTL Mapping Results
The QTL mapping results are summarized in Table 2
. We identified 38 QTL for 28 traits at the 5% genomewise level. Of the 38 QTL, 24 QTL for 17 traits were significant at the 1% genomewise level. The significant QTL are presented in Figures 1
through 7
in detail. Some of the detected QTL replicated previous findings.
|
|
|
|
|
|
For meat quality, we located a significant QTL for IMF on SSC 7, where a potential QTL for muscle moisture was detected (Figure 5
). In addition, we located both IMF and muscle moisture QTL at the same regions of SSC 9 (near SW983) and SSC 13 (near S0289; Figure 6
). On the SSC 9 QTL, the F-ratios were 10.75 for IMF and 5.18 for muscle moisture. On the SSC 13 QTL, F-ratios were 6.20 for IMF and 11.04 for muscle moisture. Paszek et al. (2001)
mapped both QTL on the same region of SSC 6. Whether IMF is directly associated with muscle moisture is not known. Uncoupling proteins 2 and 3 (UCP2 and UCP3), which are mitochondrial membrane transporters involved in thermogenesis, were mapped to SSC 9p21-p24 (Werner et al., 1999
). Because UCP2 is associated with hyperinsulinemia and obesity in mouse and human (Fleury et al., 1997
), UCP2 and UCP3 might affect lipid metabolism related to porcine IMF phenotypes.
|
|
|
Our results identified the presence of several previously identified QTL, such as for testicular weight on SSC X (Ford et al., 2001
), growth on SSC 4 (Knott et al., 1998
; Paszek et al., 1999
), vertebra number on SSC 1 (Wada et al., 2000
), and BFT on SSC 7 (Rohrer and Keele, 1998a
; de Koning et al., 1999
; Wada et al., 2000
; Malek et al., 2001b
) and on SSC X (Knott et al., 1998
; Rohrer and Keele, 1998a
). On the other hand, we detected IMF QTL on SSC 7, SSC 9, and possibly SSC 13, although IMF and marbling have also been mapped on SSC 1 (Malek et al., 2001a
) and SSC 2, SSC 4, and SSC 5 (de Koning et al., 1999
).
| Implications |
|---|
|
|
|---|
| Footnotes |
|---|
3 Present address: Ibaraki Prefectural Livestock Research Center, Yasato, Ibaraki, 315-0132, Japan. ![]()
1 Correspondence: Odakura (phone: 81-248-25-6165; fax: 81-248-25-3990; E-mail: s0sato{at}nlbc.go.jp).
Received for publication February 27, 2003. Accepted for publication August 20, 2003.
| Literature Cited |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
Y. Uemoto, Y. Nagamine, E. Kobayashi, S. Sato, T. Tayama, Y. Suda, T. Shibata, and K. Suzuki Quantitative trait loci analysis on Sus scrofa chromosome 7 for meat production, meat quality, and carcass traits within a Duroc purebred population J Anim Sci, November 1, 2008; 86(11): 2833 - 2839. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Abe, J. Saburi, H. Hasebe, T. Nakagawa, T. Kawamura, K. Saito, T. Nade, S. Misumi, T. Okumura, K. Kuchida, et al. Bovine quantitative trait loci analysis for growth, carcass, and meat quality traits in an F2 population from a cross between Japanese Black and Limousin J Anim Sci, November 1, 2008; 86(11): 2821 - 2832. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. B. Edwards, C. W. Ernst, R. J. Tempelman, G. J. M. Rosa, N. E. Raney, M. D. Hoge, and R. O. Bates Quantitative trait loci mapping in an F2 Duroc x Pietrain resource population: I. Growth traits J Anim Sci, February 1, 2008; 86(2): 241 - 253. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. B. Edwards, C. W. Ernst, N. E. Raney, M. E. Doumit, M. D. Hoge, and R. O. Bates Quantitative trait locus mapping in an F2 Duroc x Pietrain resource population: II. Carcass and meat quality traits J Anim Sci, February 1, 2008; 86(2): 254 - 266. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. A. Kuehn, G. A. Rohrer, D. J. Nonneman, R. M. Thallman, and K. A. Leymaster Detection of single nucleotide polymorphisms associated with ultrasonic backfat depth in a segregating Meishan x White Composite population J Anim Sci, May 1, 2007; 85(5): 1111 - 1119. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Sato, K. Atsuji, N. Saito, M. Okitsu, S. Sato, A. Komatsuda, T. Mitsuhashi, K. Nirasawa, T. Hayashi, Y. Sugimoto, et al. Identification of quantitative trait loci affecting corpora lutea and number of teats in a Meishan x Duroc F2 resource population J Anim Sci, November 1, 2006; 84(11): 2895 - 2901. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Kim, M. Kong, Y. J. Nam, and C. Lee A Quantitative Trait Locus for Oleic Fatty Acid Content on Sus scrofa Chromosome 7 J. Hered., September 1, 2006; 97(5): 535 - 537. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. M. Stearns, J. E. Beever, B. R. Southey, M. Ellis, F. K. McKeith, and S. L. Rodriguez-Zas Evaluation of approaches to detect quantitative trait loci for growth, carcass, and meat quality on swine chromosomes 2, 6, 13, and 18. II. Multivariate and principal component analyses J Anim Sci, November 1, 2005; 83(11): 2471 - 2481. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Mikawa, T. Hayashi, M. Nii, S. Shimanuki, T. Morozumi, and T. Awata Two quantitative trait loci on Sus scrofa chromosomes 1 and 7 affecting the number of vertebrae J Anim Sci, October 1, 2005; 83(10): 2247 - 2254. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. M. Stearns, J. E. Beever, B. R. Southey, M. Ellis, F. K. McKeith, and S. L. Rodriguez-Zas Evaluation of approaches to detect quantitative trait loci for growth, carcass, and meat quality on swine chromosomes 2, 6, 13, and 18. I. Univariate outbred F2 and sib-pair analyses J Anim Sci, July 1, 2005; 83(7): 1481 - 1493. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Nii, T. Hayashi, S. Mikawa, F. Tani, A. Niki, N. Mori, Y. Uchida, N. Fujishima-Kanaya, M. Komatsu, and T. Awata Quantitative trait loci mapping for meat quality and muscle fiber traits in a Japanese wild boar x Large White intercross J Anim Sci, February 1, 2005; 83(2): 308 - 315. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Nonneman, G. A. Rohrer, T. H. Wise, D. D. Lunstra, and J. J. Ford A Variant of Porcine Thyroxine-Binding Globulin Has Reduced Affinity for Thyroxine and Is Associated with Testis Size Biol Reprod, January 1, 2005; 72(1): 214 - 220. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |