|
|
||||||||

* Laboratory of Animal Breeding and Genetics, Faculty of Agriculture, Kobe University 657-8501, Japan and
and
Experimental Farm, Kobe University, Kasai 675-2103, Japan
3 Correspondence:
Phone: +81-78-803-5803; fax: +81-78-803-5801; E-mail:
mannen{at}kobe-u.ac.jp).
| Abstract |
|---|
|
|
|---|
Key Words: Beef Cattle Carcass Quality Cytoplasmic Inheritance Mitochondrial DNA Nucleotide Sequence Ribosomal RNA
| Introduction |
|---|
|
|
|---|
A previous study determined the sequences of the displacement loop (D-loop) region of mtDNA in order to detect the effects of mtDNA variations on carcass traits of Japanese Black cattle (Mannen et al., 1998a). These results showed significant differences between the mitochondrial haplotypes for two carcass traits of Japanese Black cattle. The D-loop region does not code any gene products, so variations of this region were used as genetic markers for the remaining regions of mtDNA.
In this study, we determined the complete mtDNA sequences from eight Japanese Black cattle in order to identify the variation of the transcribed mtDNA product that may have an influence on carcass traits.
| Materials and Methods |
|---|
|
|
|---|
This study used DNA from 202 fattening Japanese Black steers. These were the same samples used in previous works, and the mitochondrial D-loop sequences from these samples were already determined (GenBank accession numbers U87633U87650, U87893U87905) (Mannen et al., 1998a,b). These mtDNA sequences showed 26 haplotypes and were classified into five major mitochondrial types (mtDNA D-loop types 1 to 5) (Mannen et al., 1998a). In order to determine complete mtDNA sequences, we used eight samples (JBC1 to 8), which provided two samples from each of four representative haplotypes (type 1 to 4). A complete sequence of bovine mtDNA sequence determined by Anderson et al. (1982) classified the additional mitochondrial D-loop type 5 (Mannen et al., 1998a).
Sequencing
Fifty-four primers were designed and used for PCR amplification and sequencing of mtDNA based on published bovine gene sequences (accession number V00654; Anderson et al., 1982). The PCR reactions were done with 50 ng of genomic DNA as a template in a volume of 50 µL of 1x reaction buffer; 25 µM of dNTPs; 0.5 µM of each primer; and 2.5 U of EX Taq polymerase (Takara Shuzo Co., Tokyo, Japan). Amplification of PCR products was carried out using a standard PCR program with 2-min denaturation at 94°C, 30 cycles for 1 min at 94°C, 1-min annealing at 65°C, 2-min extension at 72°C, and a final extension for 7 min at 72°C. After purification of the PCR product using a GENEMate DNA purification kit (ISC BioExpress, Kaysville, UT), standard double-strand DNA cycle sequencing was performed with approximately 200 to 500 ng of amplified product using SequiTherm EXCEL II DNA sequencing kits (Epicentre Technologies, Madison, WI) and was analyzed on a LI-COR automatic sequencer (model 4200L, Lincoln, NE). Variations in complete mitochondrial genome were defined by comparison with the bovine mtDNA reference sequence published by Anderson et al. (1982).
Detection of Polymorphisms
In order to detect mtDNA candidate substitutions for the 202 Japanese Black cattle, PCR-RFLP, mismatch PCR-RFLP (Haliassos et al., 1989a, b), and PCR-single-strand conformation polymorphism (SSCP) were used depending on the mtDNA substitution. The primer sequences used to amplify the mtDNA segment that contained substitution loci are shown in Table 1
.
|
Mismatch PCR-RFLP was performed to detect polymorphisms at bp 13,310 and 14,122. Primers BMT13334R and BMT14123R contained a purposeful mismatch sequence, so that when incorporated into the PCR products, they create a BbrPI restriction site with the allele from the reference sequence but not with a mutation allele, and a NlaIII restriction site with a mutation allele, but not with the reference allele, respectively. The PCR reactions were done with 50 ng of genomic DNA in a volume of 5 µL of 1x reaction buffer; 25 µM dNTPs; 1 µM of each primers; and 0.5 U of EX Taq polymerase. Cycling conditions were 1 min at 94°C, 30 cycles for 30 s at 94°C, 30 s at 53°C (bp 13,310) or 65°C (bp 14,122), 1 min at 72°C, and final extension for 7 min at 72°C. The PCR products were digested with BbrPI (bp 13,310) or NlaIII (bp 14,122), followed by electrophoresis in a 3% MetaPhor agarose gel (BioWhittaker Molecular Applications, Rockland, ME).
The PCR-SSCP was performed to detect the polymorphisms at bp 12,158, 12,908 and 14,140 using a modification of the method described previously (Mannen et al., 2000). Briefly, forward PCR primers (BMT12074F, BMT12849F, BMT14047F) were labeled with [
-32P] ATP (6,000 Ci/mmol) and T4 polynucleotide kinase (Takara Shuzo Co., Tokyo, Japan). The PCR reactions were done with 10 ng of genomic DNA in a volume of 5 µL of 1x reaction buffer; 25 µM of dNTPs; 1 µM of each primer; and 0.5 U of EX Taq polymerase. Cycling conditions were 1 min at 94°C, 20 cycles for 30 s at 94°C, 30 s at 66°C, 1 min at 72°C, and a final extension for 4 min at 72°C. The products were mixed with 245 µL of deionized formamide solution including 10 mM EDTA, 0.1% bromophenol blue, and 0.1% xylene cyanol. The samples were analyzed by 5% polyacrylamide gel electrophoresis under controlled gel temperature at 6°C (bp 12,158) or 12°C (bp 12,908 and 14,140).
Statistical Analysis
Six carcass traitscarcass weight, longissimus muscle area (LMA), rib thickness, subcutaneous fat thickness, yield estimate, and beef marbling score (BMS)of the animals were analyzed by a mixed linear model similar to that of a previous study (Mannen et al., 1998a). The model included herd year as a fixed effect, slaughter age and inbreeding coefficient as covariates, and both nuclear additive and maternal genetic effects and residuals as random effects. In addition, the absence or presence for each of the four substitutions at a protein-coding region was considered to be a fixed effect to detect the influence of each substitution on the traits. Pedigrees of the steers were traced back to individuals born in 1975, and the additive relationship matrix of steers and their ancestors included 1,722 animals. The variance components assumed were shown in a previous study (Mannen et al., 1998a). The general linear hypothesis tests were conducted for each substitution using best linear unbiased estimators (BLUE).
| Results |
|---|
|
|
|---|
|
In the seven candidate substitutions, the substitution at bp 2,232 was observed in JBC7 and JBC8 (Table 2
), the mitochondrial D-loop type 4 as described previously (Mannen et al., 1998a). The substitution at bp 12,158 was observed in JBC1, 2, 3, and 4, which are in both D-loop types 1 and 2. The substitution at bp 12,908 was observed in JBC 5 and 6, which are in the D-loop type 3. The substitution of bp 13,310 was observed in all samples used in this study, but not in reference sequence (EU1, Table 2
). The other three substitutions (bp 14,122, 14,140, and 14,565) were detected in only one sample from eight animals.
In order to assess the relationships between the candidate substitutions and the mitochondrial D-loop types defined, we subsequently performed polymorphic DNA analyses using PCR-RFLP, mismatch PCR-RFLP, and PCR-SSCP methods for 202 Japanese Black cattle whose mitochondrial D-loop sequences had been determined previously (Mannen et al., 1998a,b). These results were summarized in Table 3
.
|
In order to confirm the effects of these substitutions on carcass traits, differences of BLUE were examined for the substitutions at bp 2,232, 12,158, 12,908, and 14,565. Differences of BLUE with and without substitutions are shown in Table 4
. The significant differences were detected on LMA and BMS at bp 2,232 substitution. No significant differences were observed for the other traits and substitutions.
|
| Discussion |
|---|
|
|
|---|
In this study, the analysis of the complete mitochondrial sequence of eight animals identified seven haplotypes (Table 2
). The comparison of eight mtDNA haplotypes with a reference sequence showed 14 variable sites in the mtDNA encoding regions. Seven of the 14 substitutions were synonymous substitutions, which would not provide changes in mtDNA gene translation. In order to focus on the distribution of the remaining seven candidate substitutions in the Japanese cattle population, polymorphic DNA analyses were performed for the substitutions in 202 Japanese Black cattle whose D-loop sequence had been previously determined (Mannen et al., 1998a,b).
Of the seven mtDNA substitutions, one was in the 16S rRNA (bp 2,232). This was a consistent substitution throughout mtDNA D-loop type 4 and the substitution was observed only in this D-loop type (Table 3
). The cattle with D-loop type 4 mtDNA tended to have better meat quality than the cattle with other D-loop types (Mannen et al., 1998a). The substitution of bp 12,158 was observed in most of D-loop type 1 animals (99.2%), all of type 2 (100%) and a few of type 5 (7.3%) (Table 3
). The D-loop types 1 and 2 have been detected exclusively in Japanese Black cattle and are closely related (Mannen et al., 1998a,b). This substitution causes the amino acid replacement of Met to Thr in the ND5 protein at AA 17. The substitution of bp 12,908 was observed in all samples of D-loop type 3 (100%) and in one animal of D-loop type 5 (2.4%). This seems to be a common polymorphism in D-loop type 3 cattle, although the number of samples with the D-loop type 3 analyzed was small. This substitution causes the amino acid replacement of Thr to Met in the ND5 protein at AA 207. These amino acids differ by having polar uncharged side chains and a nonpolar side, respectively. These substitutions are likely candidate gene polymorphisms associated with the traits of meat quality.
The substitution of bp 13,310 was observed in all samples examined in this study, but was not observed in the reference sequence (Anderson et al., 1982). The substitution causes the amino acid replacement of Lys to Thr in the ND5 protein at AA 341. A search of gene databases determined that usual amino acids at this position are Thr, Ser, Asp, or Met for most vertebrates. While Met is the most frequent amino acid in primates, Thr is most common for other vertebrates at this amino acid position. Also, the amino acids have similar properties, even though Thr, Ser, and Asp are polar uncharged side-chain amino acids, and Met has nonpolar side chains. The replacement of Lys, which has basic side chains, may be a crucial AA substitution as it has different properties to the other amino acids. Further, no vertebrate species with Lys at this position of the ND5 protein has been identified apart from the reference sequence (data not shown). This shows a limitation of the reference sequence and the lack of the substitution in the reference may be a sequencing artifact or reading mistake. Our results suggest that an amino acid at position 341 in ND5 should be Thr in cattle.
The substitutions at bp 14,122, 14,140, and 14,565 were observed in only a few samples (Table 3
). These substitutions cause the amino acid replacements of Val to Met at bp 14,122 and Phe to Leu at bp 14,140 in the ND6 protein and Phe to Leu at bp 14,565 in the cytochrome b protein. In addition, these are all nonpolar side-chain amino acids, so that the protein character would not be changed by these replacements.
Our results also revealed intraspecific comparisons to assess variability in coding regions of complete mtDNA in cattle. Interestingly, the six substitutions causing amino acid replacement observed in this study were limited to regions from bp 12,158 to 14,565 that encode ND5, ND6, and cytochrome b proteins. In addition, the rate of amino acid replacement in this study (6/13, 46.2%) is higher than that of human (12/35, 34.3%), mouse (9/47, 19.1%), and sheep (2/31, 6.5%) (Hiendleder, 1998). This is useful information as it suggests that the genetic diversity of mtDNA in cattle is greater than for other species. The reasons for these differences among species remain unclear.
Differences in BLUE were examined for the substitutions at bp 2,232, 12,158, 12,908, and 14,565, and were significant at the bp 2,232 substitution in LMA and BMS (Table 4
). There were differences of 5.33 cm2 (P < 0.01) in LMA and 0.53 units (P < 0.05) in BMS. The BMS is the most important economic trait in Japanese beef markets. Given that BMS ranges from null (0) to very abundant (5) with the average of 1.60, a difference of 0.53 units is substantially large. This substitution corresponds with mtDNA D-loop type 4 mtDNA; cattle with this D-loop type had better meat quality than the cattle with other D-loop types (Mannen et al., 1998a). These results were consistent with the findings in the previous study. This substitution may alter the function of the mitochondrial ribosome and the rate of mitochondrial protein synthesis. This substitution can be regarded as a strong candidate in the mitochondrial genome for the expression of BMS.
In conclusion, we suggest one candidate gene polymorphism from mtDNA substitutions at bp 2,232 that is related to the carcass traits of LMA and BMS in Japanese Black cattle. This finding results from the effect of mtDNA sequence and not nuclear effects from dams and sires because additive nuclear genetics effects and nuclear maternal genetic effects were considered in the mixed linear model. The representative mtDNA haplotypes from European cattle breeds were also observed in Japanese Black cattle (Mannen et al., 1988b). This suggests that these polymorphisms might be used to detect genetic diversity for carcass traits in both Japanese and European breeds.
| Implications |
|---|
|
|
|---|
| Footnotes |
|---|
2 The authors gratefully acknowledge B. Tatham (Victorian Institute of Animal Science, Australia) for the preparation of and suggestions for this manuscript. ![]()
Received for publication January 16, 2002. Accepted for publication August 20, 2002.
| Literature Cited |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
A. I. Fernandez, E. Alves, A. Fernandez, E. de Pedro, M. A. Lopez-Garcia, C. Ovilo, M. C. Rodriguez, and L. Silio Mitochondrial genome polymorphisms associated with longissimus muscle composition in Iberian pigs J Anim Sci, June 1, 2008; 86(6): 1283 - 1290. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Guan, T. Geng, P. Silva, and E. J. Smith Mitochondrial DNA Sequence and Haplotype Variation Analysis in the Chicken (Gallus gallus) J. Hered., November 5, 2007; (2007) esm094v1. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |