J. Anim Sci.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J. Anim Sci. 2008. 86:2821-2832. doi:10.2527/jas.2007-0676
© 2008 American Society of Animal Science

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
jas.2007-0676v1
86/11/2821    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Abe, T.
Right arrow Articles by Kobayashi, E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Abe, T.
Right arrow Articles by Kobayashi, E.

ANIMAL GENETICS

Bovine quantitative trait loci analysis for growth, carcass, and meat quality traits in an F2 population from a cross between Japanese Black and Limousin

T. Abe*,1, J. Saburi*, H. Hasebe*, T. Nakagawa*, T. Kawamura{dagger},2, K. Saito{ddagger},3, T. Nade*,4, S. Misumi*, T. Okumura*,5, K. Kuchida§, T. Hayashi#, S. Nakane*,6, T. Mitsuhasi*,7, K. Nirasawa*,7, Y. Sugimoto|| and E. Kobayashi*

* National Livestock Breeding Center, Nishigo, Fukushima 961-8511, Japan; and {dagger} Tokachi Branch of National Livestock Breeding Center, Otofuke, Hokkaido 080-0572, Japan; and {ddagger} Ohu Branch of National Livestock Breeding Center, Shichinohe, Aomori 039-2567, Japan; and § Obihiro University of Agriculture and Veterinary Medicine, Obihiro, Hokkaido 080-8555, Japan; and # National Institute of Agrobiological Sciences, Tsukuba, Ibaraki 305-8602,Japan; and || Shirakawa Institute of Animal Genetics, Nishigo, Fukushima 961-8061, Japan


    Abstract
 Top
 Abstract
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 LITERATURE CITED
 
A genome-wide scan for QTL affecting economically important traits in beef production was performed using an F2 resource family from a Japanese Black x Limousin cross, where 186 F2 animals were measured for growth, carcass, and meat-quality traits. All family members were genotyped for 313 informative microsatellite markers that spanned 2,382 cM of bovine autosomes. The centromeric region of BTA2 contained significant QTL (i.e., exceeding the genome-wide 5% threshold) for 5 carcass grading traits [LM area, beef marbling standards (BMS) number, luster, quality grade, and firmness), 8 computer image analysis (CIA) traits [LM lean area, ratio of fat area (RFA) to LM area, LM area, RFA to musculus (M.) trapezius area, M. trapezius lean area, M. semispinalis lean area, RFA to M. semispinalis area, and RFA to M. semispinalis capitis area], and 5 meat quality traits (contents of CP, crude fat, moisture, C16:1, and C18:2 of LM). A significant QTL for withers height was detected at 80.3 cM on BTA5. We detected significant QTL for the C14:0 content in backfat and C14:0 and C14:1 content in intermuscular fat around the 62.3 to 71.0 cM region on BTA19 and for C14:0, C14:1, C18:1, and C16:0 content and ratio of total unsaturated fatty acid content to total SFA content in intramuscular fat at 2 different regions on BTA19 (41.1 cM for C14:1 and 62.3 cM for the other 4 traits). Overall, we identified 9 significant QTL regions controlling 27 traits with genome-wide significance of 5%; of these, 22 traits exceeded the 1% genome-wide threshold. Some of the QTL affecting meat quality traits detected in this study might be the same QTL as previously reported. The QTL we identified need to be validated in commercial Japanese Black cattle populations.

Key Words: bovine • F2 family • meat quality • quantitative trait loci


    INTRODUCTION
 Top
 Abstract
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 LITERATURE CITED
 
Many studies have successfully detected QTL for economically important traits of beef cattle such as growth, carcass, and meat quality traits by using crossbred experimental populations (Keele et al., 1999Go; Stone et al., 1999Go; Casas et al., 2000Go, 2003Go; MacNeil and Grosz, 2002Go; Kim et al., 2003Go). Alexander et al. (2007aGo, b)Go recently reported the results of QTL analysis of a population for which Japanese Black and Limousin cattle were the parents.

Since 1994, we have generated an F2 resource population derived from crosses between Japanese Black sires and Limousin dams to map loci affecting economically important traits. A unique characteristic of the Japanese Black breed is the high fat content in the meat (so-called highly marbled beef), which is an important criterion for beef quality in the present Japanese market. May et al. (1993)Go described the difference in fatty acid compositions of the intramuscular fat of Wagyu crossbred and Angus steers, and Kuber et al. (2004)Go reported that Wagyu steaks had lower Warner-Bratzler shear force values than did Limousin steaks (note that most cattle known as Wagyu are Japanese Black breeds). Given those findings, Japanese Black cattle may have other economically favorable traits, in addition to marbling, compared with other breeds. In contrast to Japanese Black cattle, Limousin cattle produce leaner meat and have a larger body size. We chose these 2 breeds, the phenotypes of which differ dramatically, to construct an experimental F2 resource family for bovine QTL analysis.

In this report, we describe QTL underlying the difference in growth, carcass, and meat quality traits between Japanese Black and Limousin cattle. We incorporated physicochemical property traits of the F2 beef, including the fatty acid composition of backfat, intermuscular fat, and LM i.m. fat. We also identified QTL for computer image analysis (CIA) traits.


    MATERIALS AND METHODS
 Top
 Abstract
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 LITERATURE CITED
 
Animal care and use was according to the protocol approved by the National Livestock Breeding Center Animal Care and Use Committee.

Generation and Feeding of F2 Population

F2 Population. An F2 resource population was generated at the Tokachi and Ohu branches of the National Livestock Breeding Center in Japan. The animals used as parents were 2 Japanese Black (JB) sires (JB-A and JB-B) and 2 Limousin (L) dams (L-A and L-B). The F1 animals were obtained by crossing JB-A with L-A (family A) and JB-B with L-B (family B). Family A consisted of 2 F1 males and 17 F1 females, and family B consisted of 2 F1 males and 15 F1 females. To avoid obtaining progeny homozygous for latent recessive hereditary disease loci that may be present in the 2 JB sires, F2 animals were obtained by crossing F1 males and their nonsibling F1 females (between family A and B) using embryo transfer techniques. Both F1 families were produced and raised at the Tokachi branch, and fertilized eggs were collected in a frozen state. The majority of the frozen eggs were then sent to the Ohu branch and used to produce F2 animals. We produced 37 F2 animals at the Tokachi branch from July 1999 to January 2000. The remaining 149 F2 animals were produced at the Ohu branch during 4 periods: October to December 1999 (18 cattle); April to June 2000 (44 cattle); January to March 2001 (52 cattle); and July to October 2001 (35 cattle). The recipient dams of the F2 population consisted of 6 breeds, and ages were distributed from 2 to 13 yr (Table 1Go).


View this table:
[in this window]
[in a new window]

 
Table 1. Composition of the recipient dams of the F2 population
 
Feeding Conditions. A total of 186 F2 animals were weaned at 56 d of age. Calves were raised by artificial suckling. In the period from weaning to immediately before the fattening stage (rearing period), animals were fed with mixed feed (Snow Brand Seed Co. Ltd., Sapporo, Japan), with free intake of hay, water, and mineral salts. At 8 mo of age, F2 animals were moved to fattening stalls and began to receive the mixed feed for the fattening stage. The mixed feed comprised 30.2% barley corn, 39.4% dent corn, 15.1% wheat bran, 3.4% rice bran, 4.0% wheat flour, 5.7% soybean waste, 1.0% NaCl, and 1.0% monocalcium phosphate. The digestible CP of this diet was 11.4% and the total digestible nutrients were 83.1%. All of these percentages were calculated on a DM basis. The amount of mixed feed provided during the fattening stage was calculated considering the energy requirement given in the Japanese Feeding Standard for beef cattle (JLIA, 1995Go). To accurately control feed intake, every F2 animal was tagged with individual recognition equipment for an automatic feeding system. The allowed proportions of hay and mixed feed for the fattening stage were 25:75 for animals <14 mo old, 15:85 for animals 15 to 20 mo old, and 10:90 for animals 21 to 24 mo old. Hay, water, and mineral salts were fed without restriction.

Phenotype Measurements

Growth and Carcass Traits. The F2 animals were weighed at birth, 4 wk of age, and on the first and last day of the fattening period. Withers height (from the ground to the peak between shoulder blades), hip cross height (from the ground to the intersecting point of hip points line and the median line), body length (from lower edge of the scapula to the end of the ischial tuberosity), chest girth and width (close behind the fore foot), hip length (from the hip cross to back end of the ischial bone), hip point width (between the points of the hip), rump width, and pin bone (back end of the ischial bone) width were measured 1 d before slaughter. The F2 animals were slaughtered at the age of 24 mo (731.62 ± 5.01 d of age). At the meat plant, HCW of the F2 animals were obtained, and carcass quality was evaluated about 48 h after slaughter by certified graders belonging to the Japan Meat Grading Association (Tokyo, Japan). Graded traits were LM area, rib thickness (the length from the abdominal lining to the external side of latissimus dorsi at around the midpoint of entire rib bone of the cross section), backfat thickness (the length from the external side of latissimus dorsi to carcass surface on the vertical line from the lower end of iliocostalis to carcass surface), beef marbling standards (BMS; No. 1 to 12), beef color standards (No. 1 to 7), beef fat standards (No. 1 to 7), luster, firmness, and texture. All graded traits were measured at the sixth rib bone side of the cross section between the sixth and seventh rib bones.

Meat Quality Traits. Physicochemical property traits of the F2 beef were measured. The rib roast blocks of the seventh to eighth rib bone were sampled in all F2 animals. The LM was excised from the block and minced for analysis of moisture, crude fat, and CP content as described by Okumura et al. (2007)Go, where approximately 50 g of LM was excised and put in a plastic bag, and then incubated for 1 h in a constant-temperature bath at 70°C, and reweighed to calculate the cooking loss value. This incubated muscle was then cut thinly to yield pieces that were 1 cm x 1 cm x 5 cm cuboids) and used to measure the Warner-Bratzler shear force (Salter, Kent, UK). Meat color was measured as described by Sato et al. (2003)Go. In addition, we determined the fatty acid content of 3 parts of the rib loin block: backfat (on M. trapezius), intermuscular fat (between M. rhomboideus and LM), and intramuscular fat (of LM). Fat extractions were done as described by Folch et al. (1957)Go, and extracted fat was saponified with potassium hydrate-ethanol solution and methylesterified with boron trifluoride-methanol complex. Processed fat was analyzed by gas chromatography (6890A, Agilent Technologies Japan Ltd., Tokyo, Japan) under the following conditions: the temperature of the inlet was 150°C, the oven was warmed from 150 to 220°C, and the temperature of the detector sensor was 220°C. We used helium gas as a carrier, a capillary column (TC-70, 0.25 mm i.d. x 60 m, df (the phase thickness of the inside of the capillary column) = 0.25 µm; GL Science, Tokyo, Japan), and flame-ionization detector for detection.

Computer Image Analysis Traits. Digital images of the carcass cross section were taken between the sixth and seventh ribs using photographic equipment developed by Kuchida et al. (2001a)Go. This equipment comprised 2 parts: a dome with 570 white light-emitting diodes and a digital camera (2 megapixels, FinePix2900Z, Fuji Film, Tokyo, Japan) with a wide conversion lens (WL-FX29, Fuji Film). The distance between the camera and the surface of the carcass was fixed, and the lens was always parallel to the carcass cross section. As a result, area and length could be measured with high accuracy using the equipment. Obtained digital images were then analyzed using software developed by Kuchida et al. (2000)Go. The total muscle area, lean area, and fat area of LM, M. trapezius, M. semispinalis, and M. semispinalis capitus were calculated by this software. Here, the total muscle area represents the internal area of the muscle outline form. Therefore, the lean and fat areas are summed to give total muscle area. The ratio of the length of minor and major axes of LM was also calculated.

Genotyping

We extracted DNA from blood using automatic extraction equipment (NA1000, Kurabo, Osaka, Japan), and the final DNA concentration was adjusted to 20 ng/µL. A genome screen was conducted with microsatellite markers (Kappes et al., 1997Go; Ihara et al., 2004Go). Polymerase chain reaction amplification was performed in a volume of 15 µL containing 20 ng of genomic DNA, 1.67 mM MgCl2, 6.25 pmol of each primer, 0.2 mM de-oxynucleotides, and 0.375 U of Taq DNA polymerase (ABgene, Epsom, UK). The annealing temperatures of each marker in thermocycling steps were optimized by referencing those recommended by Ihara et al. (2004)Go. Amplifications were performed under the following conditions: 5 min at 94°C, 30 cycles of 30 s at 94°C, 30 s at annealing temperature, 30 s at 72°C, and a final extension of 7 min at 72°C. After PCR amplification, reaction products were fractionated on an ABI377 DNA sequencer (Applied Biosystems, Foster City, CA), and fragment analysis was performed with GeneScan and Genotyper software (Applied Biosystems).

Linkage Analysis

Linkage maps for the 29 bovine autosomes were constructed by using CRI-MAP (Green et al., 1990Go), and the constructed map was used for the whole-genome QTL scan. The information content of markers was calculated by the method described by Knott et al. (1998)Go.

A QTL analysis for each trait was performed by the method developed by Haley et al. (1994)Go. The statistical model is based on linear regression of phenotypes on the probabilities of QTL genotypes at a given location. We assumed that the grandparental breeds, Limousin and Japanese Black, were fixed for alternative alleles at a QTL. Two alleles at a putative QTL at a given location were denoted by Q and q. There are 3 possible genotypes, QQ, Qq, and qq, for a QTL at the given location on an autosome. The probabilities of the QTL genotypes [denoted as Prob(QQ), Prob(Qq), and Prob(qq)] were calculated from the observed genotypes of markers linked to the QTL. The calculation was done as described by Haley et al. (1994)Go. In analyses of actual data, some fixed effects other than QTL effects including sex-associated differences, breeds, and ages of the recipient cows, seasons, and locations were taken into account.

Let the effects of genotypes QQ, Qq, and qq be denoted by a, d, and –a, respectively. We assumed that the phenotypic value of a trait is written for the ith individual in F2 as follows:


Formula

where bj is the jth element of the vector of fixed effects, which includes overall mean, sex effect, breeds of the recipients (6 breeds), ages of the recipients (10 levels), and combinations of seasons and locations (5 levels); xij is the (i,j)th element of the design matrix associating bj to yi; cai is the coefficient for the additive component for individual i at the given location that is calculated from the probabilities 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 the sum of squares,


Formula

is minimized, where n is the number of individuals of F2. Denoting least squares estimators of bj, a, and d by the terms Formulaj,

Formula

, and

Formula

, the minimal sum of squares is obtained as


Formula

Under the null model corresponding to no QTL, where a = d = 0 is assumed, the minimal sum of squares is denoted by S0. Detection of a significant QTL is declared based on the ratio involving S1 and S0. In this report, we used the F-ratio, [(S0 S1)/2]/[S1/(n – 20)], as a statistic for detecting QTL, where it should be noted that degrees of freedom of 20 is assigned to the fixed effects. Significance thresholds were obtained by a permutation test with 1,000 repetitions for each trait.

Correlation coefficients among the 27 F2 phenotypes for which significant QTL were detected were calculated by PROC CORR (SAS Inst. Inc., Cary, NC).


    RESULTS AND DISCUSSION
 Top
 Abstract
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 LITERATURE CITED
 
Phenotype Measurement and Marker Selection

The 76 traits measured are summarized in Table 2Go. One trait of particular interest was the BMS number, because in the Japanese market, the value of a beef carcass is heavily dependent on this grading score. Although BMS numbers are categorized as 1 to 12 in Japanese grading systems, the highest number in our F2 population was 7. Furthermore, the distribution of the BMS numbers was extremely biased and was skewed toward lower scores (Figure 1Go).


View this table:
[in this window]
[in a new window]

 
Table 2. Performance, growth, traits of carcass grade, meat quality, computer image analysis (CIA), and fatty acid composition of F2 animals from an intercross of F1 animals derived from 2 Japanese Black sires and 2 Limousin dams
 

Figure 1
View larger version (14K):
[in this window]
[in a new window]

 
Figure 1. Bar chart of the beef marbling standards (BMS) numbers of 186 F2 animals. The x-axis indicates the BMS numbers (1 to 12), the y-axis indicates the percentage, and the numeral on each bar represents the number of F2 animals with that particular BMS.

 
The 4 parents of our F2 family were genotyped with bovine autosomal microsatellite markers (Kappes et al., 1997Go; Ihara et al., 2004Go) to select informative markers. Contrary to our expectations, the allele types were quite similar between the 2 parental breeds for many markers, so we could not help eliminating a large part of them (data not shown). Overall, we genotyped these 4 parents for a total of 1,755 markers. We selected the markers in consideration of the marker distances based on the published bovine linkage map (Kappes et al., 1997Go; Ihara et al., 2004Go) and the number of characteristic alleles that could distinguish the origin among 4 parental animals or 2 parental breeds.

Linkage Analysis and QTL Mapping

From the marker linkage analysis, 313 markers were mapped to 29 bovine autosomal chromosomes over 2,382 cM, and the average distance between markers was 8.4 cM. With this linkage map, we detected QTL on BTA2, 5, and 19 (Table 3Go). We identified 9 QTL for 27 traits at the 5% genome-wide threshold level; QTL for 22 traits were significant at the 1% genome-wide level. Details of the significant QTL are presented in Figures 2Go, 3Go, and 4Go. In the QTL analysis, we took into account the fact that our F2 population was produced at 2 different stations over different time periods. We also took into account the effects of recipient dams of the F2 population, including their breed and age, as they might have an effect on the performance of the offspring (Table 1Go).


View this table:
[in this window]
[in a new window]

 
Table 3. Summary of QTL location, genome-wide probability, and additive and dominance effects
 

Figure 2
View larger version (34K):
[in this window]
[in a new window]

 
Figure 2. Plot of the F-ratios from multilocus least squares analysis (Haley et al., 1994Go) of carcass grade and physicochemical property traits on BTA2. The x-axis indicates the relative position on the linkage map; the left-hand y-axis represents the F-ratio; and the right-hand y-axis (dotted curve) indicates information content (IC). Triangles on the x-axis indicate marker positions. Markers were MNS-2, DIK621, ILSTS026, DIK1081, DIK1140, BM4440, RM041, TGLA226, DIK1109, MM8D3, INRA135, IDVGA-37, and IDVGA-2. The horizontal lines indicate genome-wide threshold values for 5% level (dotted line) and 1% level (solid line). (A) QTL profile of carcass grade traits: {square} = LM area; {blacksquare} = beef marbling standards; {circ} = luster; {circ} = carcass grade; and + = firmness. (B) QTL profile of physicochemical property traits: {square} = CP content; {blacksquare} = crude fat content; –= moisture content; {circ} = C18:2 content; + = C16:1 content of intramuscular fat. (C) QTL profile of lean area: {square} = LM; {blacksquare} = M. trapezius; {circ} = M. semispinalis. (D) QTL profile of muscle area: {square} = LM. (E) QTL profile of ratio of fat area (RFA) to muscle area: {square} = LM; {blacksquare} = M. trapezius; {circ} = M. semispinalis capitus; – = M. semispinalis. There is no public information for marker DIK621. The primer sequences of marker DIK621 were forward primer = TCATGGCCATCATACATCAAG, reverse primer = CCCCTTTCCAAACCCATAAT.

 

Figure 3
View larger version (19K):
[in this window]
[in a new window]

 
Figure 3. Plot of the F-ratios from multilocus least squares analysis (Haley et al., 1994Go) of withers height ({square}) on BTA5. The x-axis indicates the relative position on the linkage map; the left-hand y-axis represents the F-ratio; and the right-hand y-axis (dotted curve) indicates information content (IC). Triangles on the x-axis indicate marker positions. Markers were BMS1095, BMS610, BP1, RM103, BMS1898, MS2106, CA084, ETH10, MNS-44, BMS1248, BM315, DIK2206, DIK2287, DIK2122, BM733, DIK2035, and BMS597. The horizontal line indicates threshold values for genome-wide 1% level.

 

Figure 4
View larger version (21K):
[in this window]
[in a new window]

 
Figure 4. Plot of the F-ratios from multilocus least squares analysis (Haley et al., 1994Go) for fatty acid composition traits on BTA19. The x-axis indicates the relative position on the linkage map; the left-hand y-axis represents the F-ratio; and the right-hand y-axis (dotted curve) indicates information content (IC). Triangles on the x-axis indicate marker positions. Markers were DIK2452, X82261, INRABERN148, URB044, CSSME070, BMS2389, BM17132, IOBT34, NLB-CMK39, and DIK688. The horizontal lines indicate threshold values for genome-wide 5% level (dotted line) and genome-wide 1% level (solid line). (A) QTL profile of fatty acid composition in backfat: {square} = C14:0 content. (B) QTL profile of fatty acid composition in intermuscular fat: {square} = C14:0 content; {blacksquare} = C14:1 content. (C) QTL profile on fatty acid composition of LM intramuscular fat: {square} = C14:0 content; {blacksquare} = C14:1 content; {circ} = C18:1 content; – = C16:0 content; and + = the ratio between total unsaturated fatty acid and total SFA.

 
We mapped significant QTL for 5 carcass grading traits (Figure 2Go, panel A), 8 CIA traits (Figure 2Go, panels C to E), and 5 meat quality traits (Figure 2Go, panel B) to the same centromeric region of BTA2. Among the CIA traits, LM lean area showed the greatest F-ratio of all QTL identified in this study (Figure 2Go, panel C; Table 3Go). In addition, QTL for M. semispinalis and M. trapezius lean area were detected at 4.7 and 2.0 cM, respectively, on BTA2 (Figure 2Go, panel C). Animals that inherited the Limousin alleles had larger lean area in those 3 muscles than did those that inherited the Japanese Black alleles. A QTL for LM area was detected at the same position of 1.0 cM on BTA2 (Figure 2Go, panel D); animals that inherited the Limousin allele had a greater muscle area. Because a positive correlation between muscle area and lean area in our F2 family occurred in every muscle measured by CIA (data not shown), animals that inherited the Limousin alleles at this QTL had greater muscle area with leaner meat.

In contrast, the QTL for RFA to LM, M. semispinalis, M. semispinalis capitus, and M. trapezius area were detected at 4.7 cM on BTA2 (Figure 2Go, panel E), as was the QTL for BMS number (Figure 2Go, panel A). Furthermore, QTL for crude fat content of LM and the C16:1 content of intramuscular fat were detected at 5.7 and 0 cM, respectively, on BTA2 (Figure 2Go, panel B). In all these cases, animals that inherited the Japanese Black alleles had greater values. Kuchida et al. (2000Go; 2001bGo) described significant relationships between crude fat content and RFA (r = 0.98), and between RFA and BMS number (r = 0.93), respectively. They suggested that crude fat content and RFA are useful data for evaluating marbling objectively. This suggestion is consistent with the results of our study.

All of our F2 animals showed BMS numbers of 2 to 5, except for one animal that had a 6 and one that had a 7 (Figure 1Go). This indicates the lower intramuscular fat content of the F2 population than of purebred Japanese Black cattle. In fact, Okumura et al. (2007)Go noted that the crude fat content (%) of Japanese Black cattle slaughtered at 24 mo of age (identical to the age at slaughter for our F2 population) was 37.0 ± 4.4, whereas that in our F2 population was 17.7 ± 5.9 (Table 2Go). Summarizing these results, we suggest that the Limousin alleles, which produced larger lean and muscle area, had a more extreme effect than those of the Japanese Black alleles, which produced beef with greater BMS, RFA, crude fat, and C16:1 content. Using a population that inherited the muscle hypertrophy locus, Casas et al. (1998)Go detected a QTL for ribeye area, marbling, and fat thickness. Their results were similar to ours not only in the type of phenotypes affected, but also in the QTL regions reported; the multiple QTL were in the same chromosomal region. Furthermore, using a Wagyu x Limousin crossbred F2 population, Alexander et al. (2007a)Go detected a QTL for LM area in the centromeric region of BTA2.

Recently, Sellick et al. (2007)Go reported the effect of the F94L mutation of the myostatin gene. They treated this gene as a positional candidate of the QTL for meat percentage, eye muscle area, and silverside (meat block composed of M. gluteobiceps and M. semitendinosus) percentage detected in the 0 to 15 cM region of BTA2. They analyzed a population derived from a Jersey x Limousin cross and explained that the F94L mutation of the myostatin gene originated from the Limousin breed and significantly increased these traits. Their results seem applicable to our findings. Myostatin is one of the strong candidate genes for QTL of BTA2 detected in our F2 population.

Interestingly, a QTL for C18:2 content of LM was detected at 2 cM on BTA2 (Figure 2Go, panel B). Animals that inherited the Limousin alleles had greater content of C18:2. This fatty acid is a constituent of CLA, which has recently been studied for its favorable effect on human health, especially for reducing human cancer cell growth (De La Torre et al., 2006Go). The C18:2 fatty acid in beef cannot be synthesized in the bovine body, but originates from feedstuffs. Therefore, it seems strange that the genetic effect was observed on this trait. On the other hand, several studies have reported the difference in C18:2 content between muscle and adipose of cattle (Hristov et al., 2005Go; Noci et al., 2005Go). According to these studies, C18:2 is more abundant in muscle than in adipose. In our F2 population, animals that inherited the Limousin allele had leaner meat, as described above. Furthermore, there was a weak but positive correlation between C18:2 content and both CP content and LM lean area (r = 0.11 and 0.12, respectively). Conversely, there was a negative correlation between C18:2 and both crude fat content and RFA to LM area (r = –0.12 and –0.13, respectively). These results may relate to the fact that the QTL for C18:2 content was detected at the centromeric end of BTA2.

There are several reports on candidate gene analysis of BMS and subcutaneous fat depth (SFD) traits using similar Wagyu x Limousin F2 populations. Jiang et al. (2005)Go detected genetic variation in the mitochondrial transcription factor A (TFAM) gene and determined its significant effect on both BMS and SFD. Michal et al. (2006)Go analyzed the bovine fatty acid binding protein 4 (FABP4) gene as a candidate, and found a significant relationship between detected SNP and these 2 traits. Wibowo et al. (2007)Go reported a significant effect of mutations detected in the corticotrophin-releasing hormone (CRH) gene on BMS and SFD. The first TFAM gene is located on BTA28, and the latter 2 genes (FABP4 and CRH) are both located on BTA14. We detected no significant (or suggestive) QTL for BMS or backfat thickness on those chromosomes. One possible explanation for this observation was the difference in parental individuals between those 2 F2 populations. In our F2 population, the mutations detected in these 3 genes might be fixed in the 2 breeds. In addition, the difference in the measuring procedure of BMS and SFD between the United States and Japan might have caused the different results.

We detected a QTL for withers height at 80.3 cM on BTA5 (Figure 3Go). Animals that inherited the Limousin alleles were taller than those that inherited the Japanese Black alleles. Quantitative trait loci for birth weight were detected in this region (Casas et al., 2003Go; Kim et al., 2003Go) by use of a Bos indicus x Bos taurus crossbred family. Mizoshita et al. (2004)Go detected a QTL for carcass yield on BTA5 in a half-sib population of purebred Japanese Black cattle, but the position was different from our QTL. Li et al. (2004)Go detected a QTL for preweaning ADG and ADG on feed in the 73.5- to 77.6-cM region on BTA5 using a crossbred population developed from several bovine breeds. Those investigators considered IGF-I to be a positional candidate and included information regarding IGF-I polymorphisms in their analysis. Although we measured 14 growth-associated traits in total, including birth weight and ADG during the fattening period, the genome-wise significant QTL was detected only for withers height.

On BTA19, we detected QTL for fatty acid composition (Figure 4Go, panels A to C). In addition, QTL for C14:0 content were detected at 62.3 cM for backfat (Figure 4Go, panel A); QTL for C14:0 and C14:1 content were detected at 62.3 and 71.0 cM for intermuscular fat (Figure 4Go, panel B) and at 62.3 and 41.1 cM for intramuscular fat (Figure 4Go, panel C). For each of these 3 loci, animals that inherited the Limousin alleles showed increased C14:0 and C14:1 content. We also detected QTL for C16:0 and C18:1 content and the ratio of total unsaturated fatty acid content to total SFA content of intramuscular fat at 62.3 cM on BTA 19 (Figure 4Go, panel C). Individuals that inherited the Japanese Black allele at this QTL demonstrated reduced C16:0 content, but increased C18:1 content and the ratio of total unsaturated fatty acid content to total SFA content. The degree of fatty acid composition in the intramuscular fat is an important factor for the eating quality of beef. Generally, the melting points of unsaturated fatty acids are less than those of SFAs, so beef with more unsaturated fatty acid in the intramuscular fat has superior eating quality and good texture. In addition, Mandell et al. (1998)Go suggested that C18:1 content has a favorable effect on beef flavor. In contrast, Fernandez and West (2005)Go stated that C12:0, C14:0, and C16:0 are considered to be associated with hypercholesterolemia, because they increase the concentration of low-density lipoprotein in human plasma, and Bláha et al. (2000)Go suggested that SFA concentrations and coronary atherosclerosis are related. Considering these points, we suggest that our findings here may facilitate the production of beef that is both pleasant to eat and healthier for human consumption. Furthermore, Viitala et al. (2003)Go detected a QTL for milk fat percentage at 67 cM on BTA19. Subsequently, Roy et al. (2006)Go studied the bovine fatty acid synthase (FASN) gene as a candidate gene for the QTL and found various SNP that had significant effects on milk fat percentage. Morris et al. (2007)Go detected QTL for fatty acid composition in both adipose tissue and milk fat in the 60 to 80 cM region on BTA19; the locations of those QTL overlap those that we detected. Morris et al. (2007)Go also analyzed FASN as a candidate gene for this QTL and found that the SNP haplotype had a significant effect on fatty acid composition; FASN may also be a strong candidate gene for controlling fatty acid composition in our F2 family. On the other hand, Alexander et al. (2007b)Go analyzed the fatty acid composition of the LM of their Wagyu x Limousin F2 population, and carried out QTL mapping on this trait, but did not detect significant QTL on BTA19. A possible reason for this result is the difference in parental individuals of the 2 Wagyu x Limousin reference populations. Whereas their F2 population originated from 8 Wagyu bulls and 108 Limousin females, our family was constructed from only 2 Japanese Black sires and 2 Limousin females. Our 2 Japanese Black sires were considered excellent individuals in 1995, so several favorable genes (for meat qualities) might be fixed in the 2 sires. It could be said that the structure of our F2 population was more suitable for detecting effective QTL. Interestingly, Alexander et al. (2007b)Go detected significant QTL for fatty acid composition on the centromeric region of BTA2, where we also detected the QTL on C16:1 and C18:2. The latent factor with effects on beef fatty acid composition may also be located in this region.

We observed several pairs of traits that showed highly positive or negative correlations among the 27 QTL-detected traits in this study. For example, BMS number had a strong positive relationship with RFA and LM muscle area (r = 0.85), but also showed a negative correlation with CP content (r = –0.78). Notably, QTL of those 3 traits were located at the same position, 4.7 cM on BTA2. This suggests that markers targeted for one trait may improve performance of the other trait. The opposite result might occur for other trait combinations. We may have to pay attention to this matter when consideration is given to marker-assisted selection.

Overall, the findings we report here provide fundamental information on the transmission of bovine quantitative traits. Because the QTL we detected may represent only breed-associated differences between Japanese Black and Limousin cattle, we need to confirm these QTL effects in a purebred Japanese Black population to obtain information useful in breeding Wagyu cattle.


    Footnotes
 
2 Current address: National Livestock Breeding Center, Nishigo, Fukushima 961-8511, Japan. Back

3 Current address: Agriculture, Forestry and Fisheries Research Council, Ministry of Agriculture, Forestry and Fisheries, Kasumigaseki, Tokyo 100-8950, Japan. Back

4 Current address: Nippon Veterinary and Life Science University, Musashino, Tokyo 180-8602, Japan. Back

5 Current address: Tokachi Branch of National Livestock Breeding Center, Otofuke, Hokkaido 080-0572, Japan. Back

6 Current address: Kinki Regional Agricultural Administration Office, Ministry of Agriculture, Forestry and Fisheries, Kamigyo, Kyoto 602-8054, Japan. Back

7 Current address: National Institute of Livestock and Grassland Sciences, Tsukuba, Ibaraki 305-0901, Japan. Back

1 Corresponding author: t0abe{at}nlbc.go.jp

Received for publication October 23, 2007. Accepted for publication June 27, 2008.


    LITERATURE CITED
 Top
 Abstract
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 LITERATURE CITED
 


Alexander, L. J., T. W. Geary, W. M. Snelling, and M. D. MacNeil. 2007a. Quantitative trait loci with additive effects on growth and carcass traits in a Wagyu-Limousin F(2) population. Anim. Genet. 38:413–416.[CrossRef][Medline]

Alexander, L. J., M. D. MacNeil, T. W. Geary, W. M. Snelling, D. C. Rule, and J. A. Scanga. 2007b. Quantitative trait loci with additive effects on palatability and fatty acid composition of meat in a Wagyu-Limousin F2 population. Anim. Genet. 38:506–513.[CrossRef][Medline]

Bláha, V., D. Solichová, D. Cernohorsky, M. Brátová, P. Vyroubal, and Z. Zadák. 2000. Bioanalysis of PUFA metabolism and lipid peroxidation in coronary atherosclerosis. J. Pharm. Biomed. Anal. 22:563–572.[CrossRef][Medline]

Casas, E., J. W. Keele, S. D. Shackelford, M. Koohmaraie, T. S. Sonstegard, T. P. Smith, S. M. Kappes, and R. T. Stone. 1998. Association of the muscle hypertrophy locus with carcass traits in beef cattle. J. Anim. Sci. 76:468–473.[Abstract/Free Full Text]

Casas, E., S. D. Shackelford, J. W. Keele, M. Koohmaraie, T. P. Smith, and R. T. Stone. 2003. Detection of quantitative trait loci for growth and carcass composition in cattle. J. Anim. Sci. 81:2976–2983.[Abstract/Free Full Text]

Casas, E., S. D. Shackelford, J. W. Keele, R. T. Stone, S. M. Kappes, and M. Koohmaraie. 2000. Quantitative trait loci affecting growth and carcass composition of cattle segregating alternate forms of myostatin. J. Anim. Sci. 78:560–569.[Abstract/Free Full Text]

De La Torre, A., E. Debiton, P. Juanéda, D. Durand, J. M. Chardigny, C. Barthomeuf, D. Bauchart, and D. Gruffat. 2006. Beef conjugated linoleic acid isomers reduce human cancer cell growth even when associated with other beef fatty acids. Br. J. Nutr. 95:346–352.[CrossRef][Medline]

Fernandez, M. L., and K. L. West. 2005. Mechanisms by which dietary fatty acids modulate plasma lipids. J. Nutr. 135:2075–2078.[Abstract/Free Full Text]

Folch, J., M. Lees, and G. H. S. Stanley. 1957. A simple method for the isolation and purification of total lipids from animal tissues. J. Biol. Chem. 226:497–509.[Free Full Text]

Green, P., K. Falls, and S. Crooks. 1990. Documentation for CRI-MAP, Version 2.4. Washington University School of Medicine, St. Louis, MO.

Haley, C. S., S. A. Knott, and M. Elsen. 1994. Mapping quantitative trait loci in crosses between outbred lines using least squares. Genetics 136:1195–1207.[Abstract]

Hristov, A. N., L. R. Kennington, M. A. McGuire, and C. W. Hunt. 2005. Effect of diets containing linoleic acid- or oleic acid-rich oils on ruminal fermentation and nutrient digestibility, and performance and fatty acid composition of adipose and muscle tissues of finishing cattle. J. Anim. Sci. 83:1312–1321.[Abstract/Free Full Text]

Ihara, N., A. Takasuga, K. Mizoshita, H. Takeda, M. Sugimoto, Y. Mizoguchi, T. Hirano, T. Itoh, T. Watanabe, K. M. Reed, W. M. Snelling, S. M. Kappes, C. W. Beattie, G. L. Bennett, and Y. Sugimoto. 2004. A comprehensive genetic map of the cattle genome based on 3802 microsatellites. Genome Res. 14:1987–1998.[Abstract/Free Full Text]

Jiang, Z., T. Kunej, J. J. Michal, C. T. Gaskins, J. J. Reeves, J. R. Busboom, P. Dovc, and R. W. Wright Jr. 2005. Significant associations of the mitochondrial transcription factor A (TFAM) promoter polymorphisms with marbling and subcutaneous fat depth in Wagyu x Limousin F2 crosses. Biochem. Biophys. Res. Commun. 334:516–523.[CrossRef][Medline]

JLIA. 1995. Japanese Feeding Standard for Beef Cattle. Japan Livestock Industry Association, Tokyo.

Kappes, S. M., J. W. Keele, R. T. Stone, R. A. McGraw, T. S. Sonstegard, T. P. Smith, N. L. Lopez-Corrales, and C. W. Beattie. 1997. A second-generation linkage map of the bovine genome. Genome Res. 7:235–249.[Abstract/Free Full Text]

Keele, J. W., S. D. Shackelford, S. M. Kappes, M. Koohmaraie, and R. T. Stone. 1999. A region on bovine chromosome 15 influences beef longissimus tenderness in steers. J. Anim. Sci. 77:1364–1371.[Abstract/Free Full Text]

Kim, J. J., F. Farnir, J. Savell, and J. F. Taylor. 2003. Detection of quantitative trait loci for growth and beef carcass fatness traits in a cross between Bos taurus (Angus) and Bos indicus (Brahman) cattle. J. Anim. Sci. 81:1933–1942.[Abstract/Free Full Text]

Knott, S. A., L. Marklund, C. S. Haley, K. Andersson, W. Davies, H. Ellegren, M. Fredholm, I. Hansson, B. Hoyheim, K. Lundström, M. Moller, and L. Andersson. 1998. Multiple marker mapping of quantitative trait loci in a cross between outbred wild boar and Large White pigs. Genetics 149:1069–1080.[Abstract/Free Full Text]

Kuber, P. S., J. R. Busboom, E. Huff-Lonergan, S. K. Duckett, P. S. Mir, Z. Mir, R. J. McCormick, M. V. Dodson, C. T. Gaskins, J. D. Cronrath, D. J. Marks, and J. J. Reeves. 2004. Effects of biological type and dietary fat treatment on factors associated with tenderness: I. Measurements on beef longissimus muscle. J. Anim. Sci. 82:770–778.[Abstract/Free Full Text]

Kuchida, K., M. Hasegawa, M. Suzuki, and S. Miyoshi. 2001a. Prediction of beef color standard number from digital image obtained by using photographing equipment for the cross section of carcass. Anim. Sci. J. (Jpn.) 72:321–328.

Kuchida, K., S. Kono, K. Konishi, L. D. Van Vleck, M. Suzuki, and S. Miyoshi. 2000. Prediction of crude fat content of longissimus muscle of beef using the ratio of fat area calculated from computer image analysis: Comparison of regression equations for prediction using different input devices at different stations. J. Anim. Sci. 78:799–803.[Abstract/Free Full Text]

Kuchida, K., M. Suzuki, and S. Miyoshi. 2001b. Development of photographing equipment for the cross section of carcass and prediction of BMS number by using obtained image from that equipment. Anim. Sci. J. (Jpn.) 72:224–231.

Li, C., J. Basarab, W. M. Snelling, B. Benkel, B. Murdoch, C. Hansen, and S. S. Moore. 2004. Assessment of positional candidate genes myf5 and igf1 for growth on bovine chromosome 5 in commercial lines of Bos taurus. J. Anim. Sci. 82:1–7.[Abstract/Free Full Text]

MacNeil, M. D., and M. D. Grosz. 2002. Genome-wide scans for QTL affecting carcass traits in Hereford x composite double back-cross populations. J. Anim. Sci. 80:2316–2324.[Abstract/Free Full Text]

Mandell, I. B., J. G. Buchanan-Smith, and C. P. Campbell. 1998. Effects of forage vs grain feeding on carcass characteristics, fatty acid composition, and beef quality in Limousin-cross steers when time on feed is controlled. J. Anim. Sci. 76:2619–2630.[Abstract/Free Full Text]

May, S. G., C. A. Strudicant, D. K. Lunt, R. K. Miller, and S. B. Smith. 1993. Comparison of sensory characteristics and fatty acid composition between Wagyu crossbred and Angus steers. Meat Sci. 35:289–298.[CrossRef]

Michal, J. J., Z. W. Zhang, C. T. Gaskins, and Z. Jiang. 2006. The bovine fatty acid binding protein 4 gene is significantly associated with marbling and subcutaneous fat depth in Wagyu x Limousin F2 crosses. Anim. Genet. 37:400–402.[CrossRef][Medline]

Mizoshita, K., T. Watanabe, H. Hayashi, C. Kubota, H. Yamakuchi, J. Todoroki, and Y. Sugimoto. 2004. Quantitative trait loci analysis for growth and carcass traits in a half-sib family of purebred Japanese Black (Wagyu) cattle. J. Anim. Sci. 82:3415–3420.[Abstract/Free Full Text]

Morris, C. A., N. G. Cullen, B. C. Glass, D. L. Hyndman, T. R. Manley, S. M. Hickey, J. C. McEwan, W. S. Pitchford, C. D. Bottema, and M. A. Lee. 2007. Fatty acid synthase effects on bovine adipose fat and milk fat. Mamm. Genome 18:64–74.[CrossRef][Medline]

Noci, F., F. J. Monahan, P. French, and A. P. Moloney. 2005. The fatty acid composition of muscle fat and subcutaneous adipose tissue of pasture-fed beef heifers: Influence of the duration of grazing. J. Anim. Sci. 83:1167–1178.[Abstract/Free Full Text]

Okumura, T., K. Saito, H. Sakuma, T. Nade, S. Nakayama, K. Fujita, and T. Kawamura. 2007. Intramuscular fat deposition in principal muscles from twenty-four to thirty months of age using identical twins of Japanese Black steers. J. Anim. Sci. 85:1902–1907.[Abstract/Free Full Text]

Roy, R., L. Ordovas, P. Zaragoza, A. Romero, C. Moreno, J. Altarriba, and C. Rodellar. 2006. Association of polymorphisms in the bovine FASN gene with milk-fat content. Anim. Genet. 37:215–218.[CrossRef][Medline]

Sato, S., Y. Oyamada, K. Atsuji, T. Nade, S. Sato, E. Kobayashi, T. Mitsuhashi, K. Nirasawa, A. Komatsuda, Y. Saito, S. Terai, T. Hayashi, and Y. Sugimoto. 2003. Quantitative trait loci analysis for growth and carcass traits in a Meishan x Duroc F2 resource population. J. Anim. Sci. 81:2938–2949.[Abstract/Free Full Text]

Sellick, G. S., W. S. Pitchford, C. A. Morris, N. G. Cullen, A. M. Crawford, H. W. Raadsma, and C. D. Bottema. 2007. Effect of myostatin F94L on carcass yield in cattle. Anim Genet. 38:440–446.[CrossRef][Medline]

Stone, R. T., J. W. Keele, S. D. Shackelford, S. M. Kappes, and M. Koohmaraie. 1999. A primary screen of the bovine genome for quantitative trait loci affecting carcass and growth traits. J. Anim. Sci. 77:1379–1384.[Abstract/Free Full Text]

Viitala, S. M., N. F. Schulman, D. J. de Koning, K. Elo, R. Kinos, A. Virta, J. Virta, A. Mäki-Tanila, and J. H. Vilkki. 2003. Quantitative trait loci affecting milk production traits in Finnish Ayr-shire dairy cattle. J. Dairy Sci. 86:1828–1836.[Abstract/Free Full Text]

Wibowo, T. A., J. J. Michal, and Z. Jiang. 2007. The corticotrophin releasing hormone is a promising candidate gene for marbling and subcutaneous fat depth in beef cattle. Genome 50:939–945.[Medline]


This article has been cited by other articles:


Home page
J ANIM SCIHome page
M. D. Garcia, R. M. Thallman, T. L. Wheeler, S. D. Shackelford, and E. Casas
Effect of bovine respiratory disease and overall pathogenic disease incidence on carcass traits
J Anim Sci, February 1, 2010; 88(2): 491 - 496.
[Abstract] [Full Text] [PDF]


Home page
J DAIRY SCIHome page
A. Schennink, W. M. Stoop, M. H. P. W. Visker, J. J. van der Poel, H. Bovenhuis, and J. A. M. van Arendonk
Short communication: Genome-wide scan for bovine milk-fat composition. II. Quantitative trait loci for long-chain fatty acids
J Dairy Sci, September 1, 2009; 92(9): 4676 - 4682.
[Abstract] [Full Text] [PDF]


Home page
J ANIM SCIHome page
L. J. Alexander, L. A. Kuehn, T. P. L. Smith, L. K. Matukumalli, B. Mote, J. E. Koltes, J. Reecy, T. W. Geary, D. C. Rule, and M. D. MacNeil
A Limousin specific myostatin allele affects longissimus muscle area and fatty acid profiles in a Wagyu-Limousin F2 population
J Anim Sci, May 1, 2009; 87(5): 1576 - 1581.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
jas.2007-0676v1
86/11/2821    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Abe, T.
Right arrow Articles by Kobayashi, E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Abe, T.
Right arrow Articles by Kobayashi, E.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS