|
|
||||||||
ANIMAL GENETICS |



* Departament de Ciéncia Animal i dels Aliments, Facultat de Veterinària, Universitat Autònoma de Barcelona, Bellaterra 08193, Spain;
and
Àrea de Producció Animal, Centre UdL-IRTA, Lleida 25198, Spain; and
and
Centre de Tecnologia de la Carn-IRTA, Monells 17121, Spain
| Abstract |
|---|
|
|
|---|
G181) and exon 5 (C458
G458). These two SNP are associated with Val (C)
Leu (G) and Ser (C)
Thr (G) conservative AA replacements at positions 61 and 153 of the DECR1 protein, respectively. Moreover, DECR1 genotyping in a representative sample of 184 pigs from the Large White, Piétrain, Iberian, Duroc, and Landrace breeds demonstrated the existence of disequilibrium linkage between these two SNP (Haplotype 1: C181C458; Haplotype 2: G181G458). An association analysis between DECR1 genotype and growth, carcass, and meat quality traits in a highly selected Landrace population (n = 470) revealed differences among genotypes for isocitrate dehydrogenase activity (highest posterior density [HPD] of 90%), longissimus thoracis pH (HPD of 95%), lightness (HPD of 90 to 95%), and redness (HPD of 95%). Because these associations were not consistently found in the three available genotype comparisons, we believe that exon 2 and 5 polymorphisms at the DECR1 gene might be in linkage disequilibrium with the true causal mutation influencing isocitrate dehydrogenase activity and muscle color and pH.
Key Words: Carcass Traits Fatty Acids Meat Quality Pig Quantitative Trait Loci Reductase
| Introduction |
|---|
|
|
|---|
2,4-Dienoyl CoA reductase 1 (DECR1) is a nuclear-encoded mitochondrial enzyme that participates in the ß-oxidation pathway by catalyzing the reduction of trans-2-cis-4-dienoyl-CoA to 3-enoyl-CoA (Kunau and Dommes, 1978
). The transcription unit of the human DECR1 gene includes 10 exons and nine introns of variable size that span 30 kb (Helander et al., 1997
). The deficiency of this enzyme in human has been associated with the presence of 2-trans-4-cis-decadienoylcarnitine, a metabolite produced by incomplete oxidation of the linoleic acid in urine and blood of the affected patients (Roe et al., 1990
).
The chromosomal location of the pig DECR1 gene, which coincides with the linoleic QTL previously described by Pérez-Enciso et al. (2000)
, and the crucial role of this enzyme in the ß-oxidation of PUFA, made evident the need for a more detailed characterization of the molecular features of this gene in the pig. The main objectives of our work were to identify polymorphisms in the DECR1 coding sequence and to investigate whether they are associated with phenotypic variation in carcass and meat quality traits.
| Materials and Methods |
|---|
|
|
|---|
Pigs were slaughtered when 180.0 d old (SD = 4.9 d) and weighing 104.6 kg (SD = 11.1 kg) live weight. Phenotypic records of the following traits were obtained: carcass weight, carcass length, fat thickness in the cervical region and at the last-rib level, and weight of the right and left hams, shoulders, cutlet, ribs, and bacon. Fat thickness was measured with a Fat-O-Meater (SFK Technology, Herlev, Denmark). The meat quality traits analyzed were pH and electrical conductivity determined at 45 min and 24 h after slaughter in the semimembranosus and the longissimus thoracis muscles. The pH was measured using a portable meter equipped with a xerolyte electrode (Crison, Barcelona, Spain), and electrical conductivity was measured using a Pork Quality Meter (Intek GmbH, Munich, Germany). Muscle color parameters in the CIELAB space, Lightness (L*), redness (a*), and yellowness (b*) (CIE, 1976
), were quantified in duplicate with a Minolta Chroma-Meter CR-200 (Konica Minolta, Madrid, Spain) at 24 h postmortem on the exposed cut surface of the muscle. Moreover, we analyzed the chemical composition of the muscle by measuring fat, CP, OM, and DM content from semimebranosus muscle samples (AOAC, 1990
).
Muscle samples for biochemical analyses were obtained at 24 h postmortem at the last-rib level. Samples for ELISA and enzyme activity analyses were taken from the longissimus thoracis core. They were frozen in liquid N and stored at 80°C until analysis. Samples for determination of heme pigment content were vacuum-packed and stored at 20°C until analysis.
The percentage of slow myosin heavy chain (MHC-I) in the muscle was determined with a specific MHC-I monoclonal antibody by using the ELISA technique (Picard et al., 1994
). The metabolic profile of the muscle was assessed by measuring the lactate dehydrogenase activity according to Ansay (1974)
and the isocitrate dehydrogenase (ICDH) activity according to Briand et al. (1981)
. These activities are expressed as micromoles of NADHmin1g of muscle1 (lactate dehydrogenase) and nanomoles of NADPHmin1g of muscle1 (ICDH). The concentration of heme pigment was determined according to the modified Hornsey (1956)
method. Results are given in micrograms of acid hematin per gram of fresh muscle.
Genomic DNA and RNA Extraction and cDNA Synthesis
Four hundred microliters of Tris-EDTA buffer (10 mM TrisHCl, pH = 8, 1 mM EDTA) was added to 0.4 mL of blood, and this mixture was centrifuged at 13,000 x g for 30 s. The supernatant was discarded and this washing step was repeated four or five times until a white pellet was obtained. Subsequently, cells were re-suspended in 0.4 mL of lysis buffer K (50 mM KCl, 10 mM TrisHCl, 0.5% Tween 20, and 0.1 mg/mL of proteinase K) and incubated for 5 h at 56°C. Genomic DNA was phenol-chloroform extracted and precipitated with 25 µL of 2 M NaCl and two volumes of ethanol. The genomic DNA pellet was centrifuged at 13,000 x g for 10 min, washed with 70% ethanol, and resuspended in 100 µL of Tris-EDTA.
Total RNA was extracted from ten Piétrain, Vietnamese, Large White, Iberian and Landrace pig liver samples and reverse transcribed to cDNA, as previously described (Amills et al., 2003
).
Amplification and Sequencing of the Pig DECR1 cDNA
We amplified 937 bp of the pig DECR1 cDNA by using two oligonucleotides, ENOIL-EXO2-5; 5'-AGT TTT TCA GTT ATG GGA CAA AAA-3', and DECR-3-CDNA; 5'-GAA CCT TTT GTC TTC CTG ATG AG-3'. The PCR mixture contained 1.5 mM MgCl2, 100 µM dNTP, 0.5 µM of each primer, 2 to 3 µL of the reverse transcription reaction, and 0.5 U Taq DNA polymerase (Ecogen S.R.L., Barcelona, Spain) in a final volume of 20 µL. The thermal profile comprised 35 cycles of 94°C for 1 min, 63°C for 2 min, and 72°C for 3 min. The amplified product was sequenced forward and reverse in 10 individuals with the BigDye Terminator v3.1 Cycle Sequencing kit (Applied Biosystems, Foster City, CA). Sequencing reactions were analyzed in a capillary electrophoresis device ABI PRISM 310 (Applied Biosystems, Foster City, CA). The primers used in the sequencing reactions were ENOIL-EXO2-5, DECR-3-CDNA, EN-OIL-EX3-3; 5'ATT AGG ATG TCC TGC AAC TTT GAT-3' and DECR-FW-EX5; 5'-GTG ATA AAC AAT GCA GCA GG-3'.
Genotyping of the G/C Polymorphisms at Exons 2 and 5
Primer sequences for amplifying the second exon of the DECR1 gene were ENOIL-EXO2-5 and ENOIL-EXO2-3, 5'-CAC TGA GCA CCT AGG CTG GA-3', whereas primers DECR-FW-EX5 and DECR-REV-EX5; 5'-CTT TCT GTG CTT TAA TTA GTT GC-3' were used for amplifying exon 5. Polymerase chain reactions contained 1.5 mM MgCl2, 100 µM dNTP, 0.5 µM of each primer, 30 (exon 2) or 60 ng (exon 5) of genomic DNA, and 0.5 (exon 2) or 0.75 U (exon 5) of Taq DNA polymerase (Ecogen S.R.L.) in a final 25- (exon 2) or 30-µL (exon 5) volume. The amplification of the second exon involved one denaturation step at 94° C for 1.5 min, 35 cycles at 94°C for 1.5 min, 58°C for 2 min, and 72°C for 2.5 min, and a final extension step of 72°C for 20 min. The thermal profile of the exon 5 PCR comprised 30 cycles of 94°C for 1 min, 61°C for 1 min, and 72°C for 1 min. Both polymorphisms were genotyped by primer extension analysis. The PCR products were purified with the ExoSAP-IT kit (Amersham Biosciences Europe GmbH, Cerdanyola del Vallés, Spain) and typed with the SnaPshot ddNTP Primer Extension kit (Applied Biosystems). The primers used in this typing procedure were SNAP2-DECR; 5'-CCA CCA AAT ACT TTT CAA GGA AAA-3' (exon 2), and SNAP5-DECR; 5'-CAT TAG GAG AGA GTC TTT CA-3'. The allelic frequencies of the exon 2 and 5 polymorphisms were calculated in a representative sample of 184 pigs from the Large White (n = 27), Piétrain (n = 28), Iberian (n = 22), Duroc (n = 31), and Landrace (n = 76) breeds.
Association Analyses with Carcass and Meat Quality Traits
The assumed model for the phenotypic data of each trait was as follows:
![]() |
where ß = the systematic effects (two sex effects, two Ryr1 configuration effects and three DECR1 genotype configurations), u = the vector of additive genetic effects, p = the litter effects, e = the residual vector, and X, Z1, Z2 are the incidence matrices that link phenotypic data with systematic, genetic, and permanent environmental effects, respectively. The likelihood of data is the following multivariate normal distribution:
![]() |
Prior distributions of ß were assumed flat between a range of possible values to ensure that the posterior distribution was proper. Prior distributions of the additive (u) and litter effects (p) are the following multivariate normal distributions:
![]() |
![]() |
where
and
p2 = the additive and litter variance, respectively. Prior distribution for
u2,
, and
were assumed flat between a range of possible values. Bayesian analyses were carried with the Gibbs Sampler algorithm (Geman and Geman, 1984
; Gelfand and Smith, 1990
; Tanner, 1993
) to obtain autocorrelated samples from the joint posterior density and subsequently from the marginal posterior densities of all the unknowns in the model. Specifics on distributions involved can also be found in previous studies (Wang et al., 1993
, 1994
). The posterior conditional distributions for the location parameters (sex, Ryr1, and DECR1 configuration effects) were univariate normal distributions, and the posterior conditional distributions of the variance components were inverted
2. The Gibbs sampler analysis was carried out for each analysis through a simple chain of 100,000 iterations, after discarding the first 5,000. The analysis of convergence was calculated using the algorithms of Raftery and Lewis (1992)
and García-Cortés et al. (1998)
. All iterations of the analysis were used to compute posterior means and SD invoking the ergodic property of the chain (Gilks et al., 1996
). By using this approach, all of the available information from the output of the Gibbs sampler could be considered.
The Bayesian approach considers the marginal posterior distribution of the differences between haplotypes assuming that several tests are carried out jointly; however, we did not perform any correction for multiple testing for different traits. Nevertheless, in this case, the posterior probability over (or below) zero is the probability of a difference bigger (or lower) than zero given the data, and it can be understood as a false discovery rate. For this reason, and under the assumption of Gaussian marginal posterior distributions, differences whose highest posterior density at 95% does not include the zero value have a false discovery rate lower than 2.5%, and differences whose highest posterior density at 90% does not include the zero value have a false discovery rate lower than 5%.
| Results and Discussion |
|---|
|
|
|---|
C polymorphism at exon 2 (position 181 of the coding sequence), previously described by Clop et al. (2002)
C polymorphism at exon 5 (position 458 of the coding sequence). The exon 2 and 5 polymorphisms are associated with Val (C)
Leu (G) and Ser (C)
Thr (G) conservative AA replacements at positions 61 and 153 of the DECR1 protein, respectively. We examined the allelic frequencies of both mutations. Our results indicate the existence of two segregating DECR1 haplotypes in the Landrace, Duroc, Large White, Piétrain, and Iberian pig breeds: Haplotype 1 (H1: C181C458) and Haplotype 2(H2: G181G458; Table 1
|
|
|
The associations between the polymorphism of the DECR1 gene and meat redness and pH might be explained by the fact that the two allelic variants we have detected are associated with a differential DECR1 enzymatic activity. In fact, the genomic location of this gene coincides with a QTL influencing the linoleic content, the fatty acid double-bond index, and the peroxidability index (Pérez Enciso et al., 2000
; Clop et al., 2003
). However, we do not favor this interpretation for two reasons. First, the two AA replacements identified at positions 61 (Val/Leu) and 153 (Ser/Thr) are conservative and, in principle, they are not expected to involve a dramatic change on the biochemical properties of DECR1. Second, the associations we have found are scarce and none of them are consistently found in the three available genotype comparisons. As a consequence, the more straightforward explanation would be that the DECR1 polymorphisms we have found are not the causal mutations of the associations. Exon 2 and 5 polymorphisms might be in linkage disequilibrium with the true causal mutation influencing ICDH and muscle color and pH, which may lie in another region of the DECR1 gene or even in a neighboring locus.
| Implications |
|---|
|
|
|---|
| Footnotes |
|---|
3 M. Amills and O. Vidal contributed equally to this work. ![]()
2 Correspondencephone: 34-93-5812876; fax: 34-93-5812106; e-mail: marcel.amills{at}uab.es.
Received for publication August 4, 2004. Accepted for publication November 19, 2004.
| Literature Cited |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. Lei, C. Luo, X. Peng, M. Fang, Q. Nie, D. Zhang, G. Yang, and X. Zhang Polymorphism of Growth-Correlated Genes Associated with Fatness and Muscle Fiber Traits in Chickens Poult. Sci., May 1, 2007; 86(5): 835 - 842. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |