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


     


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary Data
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 Cagnazzo, M.
Right arrow Articles by Russo, V.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Cagnazzo, M.
Right arrow Articles by Russo, V.
J. Anim. Sci. 2006. 84:1-10
© 2006 American Society of Animal Science


ANIMAL GENETICS

Comparison of prenatal muscle tissue expression profiles of two pig breeds differing in muscle characteristics1

M. Cagnazzo*,{dagger},2, M. F. W. te Pas{dagger},2, J. Priem{dagger}, A. A. C. de Wit{dagger}, M. H. Pool{dagger}, R. Davoli*,3 and V. Russo*

* DIPROVAL University of Bologna, Sezione Allevamenti Zootecnici, Italy; and and {dagger} Wageningen University and Research Centre, Animal Sciences Group, ID-Lelystad, Division of Animal Resource Development, Animal Genomics Group, Lelystad, The Netherlands


    Abstract
 Top
 Abstract
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
The objective of this study was to compare purebred Duroc and Pietrain prenatal muscle tissue transcriptome expression levels at different stages of prenatal development to gain insight into the differences in muscle tissue development in these pig breeds. Commercial western pig breeds have been selected for muscle growth for the past 2 decades. Pig breeds differ for their muscle phenotypes (i.e., myofiber numbers and myofiber types). Duroc and Pietrain pig breeds are extremes; Duroc pigs have redder muscle fiber types with more intramuscular fat, and Pietrain pigs have faster-growing and whiter muscle fiber types. Pietrain pigs are more muscular than Duroc pigs, whereas Duroc pigs are fatter than Pietrain pigs. The genomic background underlying these breed-specific differences is poorly known. Myogenesis is a complex exclusive prenatal process involving proliferation and differentiation (i.e., fusion) of precursor cells called myoblasts. We investigated the difference in the prenatal muscle-specific transcriptome profiles of Duroc and Pietrain pigs using microarray technology. The microarray contained more than 500 genes affecting myogenesis, energy metabolism, muscle structural genes, and other genes from a porcine muscle cDNA library. The results indicated that the expression of the myogenesis-related genes was greater in early Duroc embryos than in early Pietrain embryos (14 to 49 d of gestation), whereas the opposite was found in late embryos (63 to 91 d of gestation). These findings suggest that the myogenesis process is more intense in early Duroc embryos than in Pietrain embryos but that myogenesis is more intense in late Pietrain fetuses than in Duroc fetuses. Transcriptomes of muscle structural genes followed that pattern. The energy metabolism genes were expressed at a higher level in prenatal Pietrain pigs than in prenatal Duroc pigs, except for d 35, when the opposite situation was found. Fatty acid metabolism genes were expressed at a higher level in early (14 to 49 d of gestation) Duroc embryos than in Pietrain embryos. Better understanding of the genomic regulation of tissue formation leads to improved knowledge of the genome under selection and may lead to directed breed-specific changes in the future.

Key Words: embryo • expression profile • microarray • muscle • pig breed • transcriptome


    INTRODUCTION
 Top
 Abstract
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Pig breeding during recent decades has focused on improving growth rate and muscularity (Merks, 2000Go). As a result, pig breeds differ in muscle traits such as muscularity, muscle fiber type, color, etc. Duroc pigs are slower growing with relatively high i.m. fat content (Sellier, 1998Go) and redder muscle fiber types. Pietrain pigs are faster-growing pigs with relatively low i.m. fat content (Jones, 1998Go; Sellier, 1998Go) and whiter muscle fiber types. In addition, overall fatness of Duroc pigs is greater than that of Pietrain pigs. These 2 breeds are considered to represent extremes of modern Western pig breeds.

Many genes regulating myogenesis are known. Often these genes affect the activity and/or expression of the muscle regulatory factors (MRF) and myocyte-specific enhancer binding factor2 (Mef2) gene families, which take central positions in the regulation of myogenesis. The MRF gene family consists of 4 genes regulating proliferation (myf-5 and MyoD), terminal differentiation (myogenin), and maintenance of muscle tissue (Olson, 1990Go; Weintraub et al., 1991Go). The Mef2 genes (Mef2-A,-B,-C,-D) are members of the MADS gene family, and Mef2-D is present in undifferentiated myoblasts and may participate in the earliest commitment events leading to myogenesis (Wagenknecht et al., 2003Go). The Mef2 proteins are expressed in many tissues, but Mef2 only activates transcription in developing muscle.

Using microarray technology, we studied expression of the many known genes simultaneously affecting myogenesis. Microarray technology can simultaneously measure the differential expression of a large number of genes in a given tissue and may identify the genes related to the different muscle phenotypes. The aim of this study was to compare purebred Duroc and Pietrain prenatal muscle tissue transcriptome expression levels at different stages of prenatal development to get insight into the differences in muscle tissue development in these pig breeds.


    MATERIALS AND METHODS
 Top
 Abstract
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Animals
Eighty-nine purebred Duroc and 87 Pietrain sows or gilts were bred with Duroc or Pietrain purebred boars, respectively. The sows were slaughtered around d 14, 21, 35, 49, 63, 77, and 91 after insemination, and the embryos/fetuses were collected. The litters of at least 14 sows per gestational age were collected. All handling of sows, embryos, and fetuses was done according to Dutch laws. The gestational ages of the litters varied for 4 d around the sampling ages. Fourteen-day-old embryos were collected by flushing the uterus. Thus, these samples represented all the embryos in a uterine horn. Twenty-one-day-old embryos were individually collected as whole embryos. Embryonic length was determined for each embryo. Fetuses of 35 to 91 d were weighed. For embryos with gestational age of 49 to 91 d, the gender was determined. Although the LM itself was not clearly visible at 35 d of development, the area where the muscle developed appeared and was extracted. Longissimus muscle tissue was dissected from the 49- to 91-d-old embryos. All collected materials were snap frozen in liquid N2 and stored at –80°C until RNA isolation. The gestational ages covered the whole period of myogenesis.

Microarray Construction
The microarray contained PCR fragments of 509 genes. A literature survey resulted in 200 genes with known effects on myogenesis, energy metabolism, and muscle structural genes. Some of the genes were members of larger gene families (See Supplemental Table; available online at http://jas.fass.org). A few other members of such gene families also were included in the microarray. The genes were partially cloned using a pool of RNA consisting of 2 samples per gestational age. Gene-specific primers were constructed based on sequence information from GenBank and used to amplify a PCR fragment on this cDNA pool. The fragments were cloned and partially sequenced to verify identity when doubt on identity existed. Furthermore, 309 clones of an adult pig muscle cDNA library (Davoli et al., 1999Go, 2002Go) were added. GenBank accession numbers, references, primers, and pathway information of all genes on the microarray are supplied in the Supplemental Table.

The clones were amplified in four 50-µL PCR reactions to obtain 200 µL total volume of PCR product for each clone, collected in a single well of a 96-well microtiter plate, and purified through Sephadex G100 (Sigma Aldrich, Zwijndrecht, The Netherlands). The purified products were precipitated and resuspended in 20 µL spotting buffer (1M phosphate buffer, pH 5.8; 50% dimethyl sulfoxide). After that, the length and quality of each fragment were checked on a 1% agarose gel. The fragments were spotted on glass slides in duplicate.

Microarray Hybridization
The Trizol-Phenol (Life Technologies, Breda, The Netherlands) method was used to isolate RNA from 2 embryos or fetuses per breed, except for the 14-d gestational age stage, for which Qiagen RNeasy kit (Qiagen, Cologne, Germany) was used. The RNA samples were pooled, and 2 µg of RNA for each stage and breed was labeled with Cy3 or Cy5 using the TSA labeling and amplification kit protocol (Perkin Elmer Life Sciences, Inc., Langen, Germany). Each microarray was hybridized with samples from a single gestational age with a sample from Duroc labeled with one dye and a sample from Pietrain labeled with the other dye. For each gestational age, hybridizations were performed in duplicate and in duplicate dye swap. A total of 7 gestational ages were investigated; thus, the number of microarray hybridizations was 28. Hybridization was done for 16 h at 65°C. After hybridization, the slides were rinsed according to the stringent washing protocol recommended by the manufacturer.

Microarray Analysis
Microarrays were scanned using the GeneTac2000 scanner (Genomics Solutions, Ann Arbor, MI). Each microarray was analyzed independently using the following steps. The first step was normalization of raw scanning data, which included background correction. The background signal was determined using blank spots. The background was normalized by including 33 blank spots per patch, which were evenly distributed. Blank spots are spots for which only spot buffer was used for printing and no DNA material. The blanks spots are therefore a better reference to account for nonspecific expression than local background surrounding expressed spots. Beyond that, we found that the number of lost spots attributable to negative values after background correction was less using blank spots for background correction. Normalization also was done 1) using all spots, and 2) per patch, all intensity-dependent using a locally weighted scatterplot smoothing (LOWESS) fit (Cleveland, 1974Go; Park et al., 2003Go), following the procedure described by Yang et al. (2002)Go and Pool et al. (2003)Go. Spots are represented by their M- and A-values: M = log2 (Cy5/Cy3); thus, the M-value indicates the differential expression between the Cy5 and Cy3 labeled RNA samples; and A = (log2 [Cy5 * Cy3])/2. Therefore, the A-value indicates a weighted mean expression level of the Cy5- and Cy3-labeled RNA samples (Yang et al., 2002Go). Additionally, the significance of the difference is indicated by the 2-sided P-value tested on a log-logistic distribution.

The second step in microarray analysis was that normalized spots with the difference Cy3 – Cy5 having P > 0.05 were discarded. The third step was that the remaining spots were analyzed for up- and/or downregulation of expression by comparing the differential expression in the breeds using the M-values, and for mean expression levels (i.e., A-values). The M-values were regarded significant if –1.58 < M > 1.58 (i.e., if the differential expression is larger by at least 3 times, as recommended by the manufacturer of the labeling kit; Perkin Elmer). Using the Spotfire Pro 7 software (BioASP, Amsterdam, The Netherlands), clustering was performed with the K-means clustering option after a principal component analysis using the M- and A-values. Clustering is based on A-values using Euclid distances. The software was arbitrarily set to divide the spots into 8 clusters. Finally, the last step was that the results were analyzed with biological interpretation of data using available data on the physiology of the genes (i.e., involved in energy metabolism, myogenesis, or muscle structural genes).

The P-values were calculated without accounting for multiple comparisons; thus, there is the need for protection against Type I errors. The Bonferroni test was not used because it is too conservative, especially for microarray data. We avoided the risk of Type I errors by considering multiple copies of the clones on the array. The false error rate of a single spot of a gene is 0.05 and is 0.05 x 0.05 = 0.0025 for 2 spots. If no effect (false positive), positive measurements are expected to occur only in 1 spot of a gene (with the chance to occur in 2 spots also being extremely small). To be included in the list of up- or downregulated genes, a gene should show 1 to 4 spots significant per microarray on at least 2 microarrays, of which at least 1 microarray should be a dye swap microarray, making a Type I error extremely unexpected. Therefore, finding multiple significant up-or downregulated copies of a gene on the array increases the overall chance of a gene to be differentially expressed on the array (i.e., the more copies are found significant, the smaller the risk of detecting falsely significant designated genes due to Type I errors). Furthermore, the chance of having a false positive upregulated copy would be as large as for a false positive downregulation. If the experiment had no power at all, one would expect as many false positive down-regulations as up-regulations. We have accounted for this by suggesting that a group of genes only has an effect when the number of either up- or down-regulated genes is more than one-half of the genes within that group. Therefore, we believe that we avoided the risk of Type I errors in this experiment.

Quantitative Real-Time PCR
To validate the results of the microarrays, 5 genes showing differential expression were selected and analyzed with real-time PCR using the Lightcycler (Roche Diagnostics, Almere, The Netherlands). Genes were chosen in different pathways: myogenesis-affecting (EPO-receptor, ß-catenin, and TGF ß2), energy metabolism (GAPDH), and muscle structural (collagen 3A1). Primers were designed on the cDNA sequence to amplify a 100 to 150 bp fragment, and probes containing fluorescein were designed according to the rules set by the manufacturer (Table 1Go). All reactions had an annealing temperature of 60°C, except TGF-ß2 (55°C), and a magnesium concentration of 3 mM, except ß-catenin (5 mM). For real-time PCR, each gestational age was represented by individual RNA samples isolated from 6 embryos or fetuses different from the animals and litters used for microarray analysis. Mean values are presented.


View this table:
[in this window]
[in a new window]
 
Table 1. Primers and probes used for real-time PCR verification of 5 genes showing differential expression in the microarray analysis
 
Quantitative results are usually presented normalized against an internal control, often house-keeping genes such as GAPDH or beta-actin. In our results, we show that these genes are also highly regulated during these prenatal differentiation events. It is clear that such genes cannot be used to normalize the expression of other genes (Radonic et al., 2004Go). In the microarray, normalization as a function of all genes on the microarray was used, which cannot be done with real-time PCR. Therefore, we also tried 18S rRNA; however, we showed that 18S rRNA also is strongly regulated during the 14-and 21-d gestational ages, as has been shown previously (Voronina, 2002Go; Radonic et al., 2004Go). Using 18S rRNA normalization, the 35- to 91-d period resulted in unrealistic differences between 21 and 35 d because of extremely low 18S rRNA levels (data not shown). This period is a very important stage in primary muscle differentiation. Therefore, we decided to use nonnormalized real-time PCR values. The mean value of 6 animals is presented. The results are presented as the log(Duroc/Pietrain); the sign of the results indicates the ratio and the direction of the differential expression between Duroc and Pietrain. It is important to note that no new data are produced with real-time PCR. These experiments were only used to verify the microarray data.


    RESULTS
 Top
 Abstract
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Microarray Analysis
Genes were grouped into 3 major groups: myogenesis, energy metabolism, and muscle structural genes. The first 2 groups have been subdivided into pathway-specific subgroups (Table 1Go). Results were analyzed for 1) up- or downregulation (i.e., the ratio between the expression level in Duroc and the expression level in Pietrain; M-value), and 2) for general expression level (A-value).

Duroc-Pietrain Expression Ratio
Data were analyzed for spots with P < 0.05 and for spots showing more than 3 times differential expression. Although the number of spots is substantially lower in the latter analysis, the results for each group are similar. The results are shown in Table 2Go. Column N shows the total number of genes on a microarray for each (sub)group, whereas column N-list indicates the total number of genes for each (sub)group of genes showing statistically significant differential expression on at least 1 gestational age. To be included, a gene should show 1 to 4 spots significant per microarray on at least 2 microarrays, of which at least 1 microarray should be a dye swap microarray. A (sub)group was considered significant when at least 50% of the genes was differentially expressed on at least 1 gestational age. Note that especially in the myogenesis group, genes may be included in 2 subgroups (e.g., differentiation-inhibiting and proliferation-stimulating effects).


View this table:
[in this window]
[in a new window]
 
Table 2. Numbers of genes showing differential expression within groups of genes related to energy metabolism, myogenesis, and muscle structure at 7 gestational ages covering the whole period of myogenesis
 
Energy Metabolism.
The results (Table 2Go) indicate a major difference between the energy metabolism in Duroc and Pietrain embryos and fetuses at all gestational ages. At all gestational ages except 35 d of gestation, the energy metabolism in the Pietrain was at a greater level than in the Duroc. At d 35 the situation was reversed.

The genes of the fatty acid metabolism had a different general expression profile. The results indicate that fatty acid metabolism was at a higher level in early Duroc embryos (d 14 to 49 of gestational age) compared with Pietrain embryos, whereas the reverse situation was found in older fetuses from d 63 of gestation and onward.

Finally, with the exception of the ATP metabolism subgroup, all subgroups had lower numbers of differentially expressed genes at 91 d of gestation. This finding may indicate that at birth the Duroc and Pietrain breeds differ less in energy metabolism.

Myogenesis.
Table 2Go indicates that, in general, myogenesis started earlier in Duroc than in Pietrain because the expression levels of all myoblast proliferation and differentiation affecting groups of genes were greater in Duroc embryos of 14 to 35 d than the expression levels in Pietrain. From 49 d and onwards, with the exception of the differentiation-inhibiting group of genes, Pietrain showed increased myogenesis. It should be noted that not all functional groups showed significant differential expression on all gestational ages.

Muscle Structural Genes.
The expression pattern of the muscle structural genes showed similarities with the myogenesis genes. In the early embryo until d 49, the expression of the genes was greater in Duroc than in Pietrain. Fetuses older than 49 d of gestation showed the opposite expression pattern.

Clustering Analysis
Clustering analysis was performed to investigate the absolute expression levels of the genes in both breeds in relation to function. A typical clustering analysis with principal component analysis (PCA) is shown in Figure 1Go, in which the PCA1 is the A-value and the PCA2 is the M-value. The principal component analysis and clustering analysis decreased the amount of individual data by grouping the genes according to their differential expression patterns.



View larger version (29K):
[in this window]
[in a new window]
 
Figure 1. A typical example of principal component K-means clustering analysis. The software was arbitrarily set to divide the spots into 8 clusters. Clustering was based on A-value (i.e., a weighted mean of absolute expression levels in the Cy3- and Cy5-labeled samples). Principal component analysis (PCA) is shown for PCA1 (A-value) and PCA2 (M-value).

 
Clusters are created based on the A-value of the spots, indicative for average expression levels in both breeds. The results as summarized in Table 3Go indicate that genes with similar biological function often cluster together, indicating pathway-regulated expression levels. Genes involved in myogenesis seemed to be an exception; however, myogenesis is a multiprocess pathway as indicated in the subgroups of Table 2Go. Although less clearly interpretable, there seems to be a coclustering of proliferation-stimulating and differentiation-inhibiting genes. Genes of other functional groups cluster more widespread in different clusters (data not shown).


View this table:
[in this window]
[in a new window]
 
Table 3. Cluster numbers of genes of similar biological functions at 7 gestational ages
 
Real-Time PCR Validation of Microarray Results
To validate the microarray results, we chose 5 genes, 3 myogenesis related (EPO-receptor, ß-catenin and TGF-ß2) and all belonging to both proliferation-stimulating and differentiation-inhibiting myogenesis sub-groups, an energy metabolism household gene (GAPDH) more specifically belonging to the glycolysis subgroup, and a muscle structural gene (collagen 3A1), a gene expressed early during myogenesis (i.e., from d 14 and onward) to represent the groups of genes.

Real-time PCR validated the results of the microar-rays (Figure 2Go). The results of the real-time PCR for each gene are shown in Figure 2Go as the log(Duroc/Pie-train). Presented in this way, the results directly indicate the direction of the differential expression: if the expression level in Duroc fetuses was greater than in Pietrain fetuses, the log of the ratio is positive (indicated with a bar to the right); the opposite differential expression is indicated by a bar to the left.



View larger version (21K):
[in this window]
[in a new window]
 
Figure 2. Real-time PCR of 5 genes representing myogenesis (EPO-receptor, ß-Catenin, and TGF-ß2), energy metabolism (GAPDH), and muscle structural (Col3A1) genes. The logs of the Duroc/Pietrain ratio of the expression are presented. If the expression in Duroc was greater than in Pietrain, the log of the product is >0. The reverse situation is indicated by a product that is <0.

 
The EPO-receptor quantitative real-time PCR showed that for all gestational ages, the expression in Duroc was greater than in Pietrain, except at d 49. Differences in expression levels were greatest at d 14. These results match best with the microarray results for the myogenesis differentiation-inhibiting subgroup. The results for ß-catenin indicate that the expression in Pietrain fetuses was greater than in Duroc fetuses, except at d 21. The difference in expression level was greater in early (14 and 21 d) embryos than in older fetuses. The TGF-ß2 results indicate that in young embryos, expression in Duroc was greater than in Pietrain, but the situation reversed from d 49 and onwards; embryos at 35 d of gestational age showed almost equal expression levels. The results for both ß-catenin and TGF-ß2 matched best with the microarray myogenesis proliferation-stimulating subgroup. Small differences are found between the microarray results and the results of the real-time PCR. The microarray showed no significant differences for 14, 49, and 63 d, whereas small differences were observed in the real-time PCR for both genes. This result may be related to the limits of the detection of differential expression of the microarray technology compared with real-time PCR technology.

The quantitative real-time PCR differential expression of Col3A1 indicated greater Duroc expression levels in 14- to 35-d-old embryos compared with Pietrain embryos, with the opposite being the case for older fetuses. The Col3A1 gene expression matched with the microarray muscle structural genes expression patterns, except for d 35 and 49. Because Col3A1 is expressed very early in prenatal development, the ratio-switch in expression from Duroc to Pietrain also may be early, so instead of after d 49, the ratio-switch in the expression of this gene occurs before d 49. It also should be noted that at d 14 of gestation, the gene is expressed in Duroc, but even real-time PCR was unable to detect expression in Pietrain.

Quantitative real-time PCR of GAPDH showed that the expression level was greatest in Pietrain fetuses, except for d 21. These results are similar to the microarray results, except that the changes in the expression in the microarray were at d 35, whereas the real-time PCR indicated the changes were at d 21.

In general, the results of the real-time PCR indicated that the microarray data were reliable for these 5 genes. Therefore, we extrapolate these findings to the reliability of all microarray data.


    DISCUSSION
 Top
 Abstract
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Genetic Regulation of Myogenesis
The formation of new skeletal muscle fibers or myogenesis is an exclusively prenatal process, determining muscle characteristics such as fiber numbers, which may be related to muscle strength and function (Rehfeldt et al., 2000Go). Myogenesis is under complex genetic regulation. The MRF and Mef2 gene families of transcription factors take central positions in the regulation of myogenesis. The MRF gene family consists of 4 genes regulating proliferation (myf-5 and MyoD), terminal differentiation (myogenin), and maintenance of muscle tissue (Olson, 1990Go; Weintraub et al., 1991Go). The Mef2 genes (Mef2A,-B,-C,-D) are members of the MADS gene family. Their proteins bind to a conserved A/T-rich sequence in the control regions of numerous muscle-specific genes. Although Mef2A and Mef2C proteins are specific to differentiated skeletal and cardiac muscle, Mef2D is present even in undifferentiated myoblasts, and it may participate in the earliest commitment events leading to myogenesis. (Wagenknecht et al., 2003Go). The Mef2 proteins are expressed in many tissues, but it is only in the developing cardiac, skeletal, and smooth muscle that Mef2 activates transcription. After birth, the expression of Mef2 is downregulated, but it is reinduced during skeletal muscle regeneration (Parker et al., 2003Go). In pigs, genetic variation in myogenin is related to muscle and general body growth (Soumillion et al., 1997Go; te Pas et al., 1999bGo), whereas no such relationship was reported for myf-5 (te Pas et al., 1999aGo). Furthermore, postnatal expression levels of myf-5, MyoD, and myogenin also were related to growth rate (te Pas et al., 2000Go). The activity and expression levels of the MRF genes are regulated by a network of genes together modulating skeletal muscle fiber development and differentiation and thereby regulating muscle and body growth potential. The MRF proteins act in concert with the Mef2 family in orchestrating differentiation.

Mammalian myogenesis takes place in 2 distinct waves: the primary and secondary waves of muscle fiber formation in the pig occur from d 30 to 60 and from d 54 to 90 of gestation, respectively (Wigmore and Stickland, 1983Go; Wigmore and Evans, 2002Go). Our samples are composed of both of these periods and the period preceding the primary muscle fiber formation when the myoblasts, precursor cells for myogenesis, are formed. It should be noted that especially the very young embryonic samples differed from the fetal samples that were older. The 14-d embryo samples were pools of all embryos in a uterine horn, and the 21-d embryos were composed of whole embryo samples. All other samples were LM specific. Therefore, comparison of the data of different gestational ages within the range of the young embryos should be treated with caution; however, we only present the comparison of the expression patterns between breeds for each specific gestational age. Because these data are derived from similar isolates, these comparisons are not affected by these sample differences. Nonetheless, it should be kept in mind that the nonmuscle tissue in these samples might affect the observed expression level, indicating expression that probably was not related to myogenesis. Because we selected genes involved in myogenesis, this is unlikely but cannot be excluded.

It is well known that pig-breed-specific differences in muscle fiber numbers exist (Dwyer and Stickland, 1991Go); however, the molecular background of this difference has never been investigated. Therefore, we compared the myogenesis of 2 pig breeds differing for muscle fiber characteristics. The high muscularity of the Pietrain pigs is mainly related to hypertrophy of muscle fibers leading to mostly white muscle fiber typing (Jones, 1998Go; Sellier, 1998Go). Contrary to this, Duroc pigs have smaller, redder fiber types (Sellier, 1998Go). The prenatal muscle tissue transcriptomes of both breeds were compared using microarray technology.

We performed a literature search to find as many genes as possible that have been related to myogenesis using data of all animal species and from in vitro studies. The pig homologues of these genes were cloned and placed on a microarray to investigate their expression profiles simultaneously. Because energy availability in the developing muscle tissue during myogenesis seemed to affect myogenesis, we included genes from several energy metabolism pathways on the microarray. Energy availability might be important to form mature muscle fibers that also contain enough energy for the period shortly after birth. If this fails, the perinatal muscles might be insufficient, leading to muscle failure such as splay leg and other leg weaknesses (Jorgensen and Vestergaard, 1990Go). Furthermore, a direct connection between energy metabolism and myogenesis has recently been suggested because energy content of satellite cell cultures has been related to hypertrophy of muscle fibers in vitro (Louis et al., 2004Go). Adding creatine to satellite cell cultures increases IGF-I level and expression of the 4 MRF genes, although only myogenin was increased more than 3-fold. Moreover, differential regulation of the expression of 6-phosphofructo-2-kinase — key regulatory enzymes in glycolysis and gluconeogenesis (Riera et al., 2003Go) and BMP9 — regulating the expression of phosphoenolpyruvate carboxykinase (PEPCK) and Akt kinase activity, a key enzyme in glycolysis and a key regulatory protein in differentiated myotubes, respectively (Chen et al., 2003Go), indicates a direct connection between energy metabolism and myogenesis.

Comparison of Myogenesis in Prenatal Duroc and Pietrain
Expression profile differences of myogenesis-related genes suggest that muscle fiber formation in Duroc was initiated earlier than in Pietrain because myogenesis-affecting genes were expressed to a greater level in young (14 to 35 d) prenatal Duroc pigs than in Pietrain pigs. The expression profile differences of muscle structural genes was similar, also suggesting that muscle fiber formation is taking place earlier in Duroc embryos compared with Pietrain embryos. Delayed muscle fiber formation has been related to increased muscle fiber numbers in selected quail lines (Coutinho et al., 1993Go). The delay may be associated with a prolonged proliferation period, not necessarily proliferation rate, which may be unaffected, giving rise to an increased numbers of myoblasts; this may result in increased primary muscle fiber numbers. Increased numbers of primary muscle fibers may give rise to increased secondary muscle fiber formation. Muscle fiber numbers have been related to muscle mass (Stickland and Goldspink, 1973Go; Dwyer et al., 1993Go); however, the increased muscle mass of Pietrain pigs has been related more to hypertrophy rather than to hyperplasia (Jones, 1998Go; Sellier, 1998Go). Alternatively, the delayed muscle fiber formation may result in the formation of a larger pool of satellite cells. The increased proliferation process during later developmental ages in Pietrain as compared with Duroc may enhance the number of satellite cells further, thereby increasing the hypertrophic capacity of Pietrain muscles further. Because hypertrophy is mainly a postnatal trait of muscle fibers, this seems to be related to the low expression level of energy metabolism genes relative to Duroc prenatal pigs. A high creatine level, suggestive of a high energy level, was reported to be associated with hypertrophic growth (Louis et al., 2004Go). Although speculative, this also may relate to the postnatal Pietrain leanness. Postnatal Pietrain pigs, which have highly hypertrophied muscles, may use a high energy metabolism to increase muscle hypertrophy, leaving less energy for fat deposition.

Clustering analysis indicated that the expressions of the genes involved in several myogenic processes are more or less similarly regulated. In particular, the inhibition of differentiation and stimulation of proliferation were related. This provides further molecular evidence that these 2 processes are often coregulated. Differentiating myoblasts leave cell division before terminal differentiation (Olson et al., 1991Go). Finally, expression levels of myogenesis affecting genes are often low compared with genes affecting energy metabolism. Considering that many of these genes are transcription factors involved in the regulation of the MRF genes, which also are transcription factors, and many energy metabolism genes are so-called household genes, this may not be surprising. Nonetheless, because clustering is based on A-values (i.e., mean expression level), this observed difference in expression level may interfere with clustering analysis. Because of this difference, genes of the energy pathways and the myogenesis pathways often do not cluster together, which may provide an alternative explanation for the observed clustering patterns. On the other hand, muscle structural gene expression level depends on the progression of myogenesis, so expression is low in early embryos and increases in time. Nevertheless, this expression pattern also does not interfere with pathway clustering.

Comparison of Energy Metabolism in Prenatal Duroc and Pietrain in Relation to Myogenesis
In general, the expression levels of the genes involved in energy metabolism were greater in prenatal Pietrain than in Duroc muscle tissue, except for d 35, where the opposite situation was found for all energy metabolism pathways. The period around d 35 of gestation is central to primary muscle fiber formation (Wigmore and Stickland, 1983Go; Wigmore and Evans, 2002Go). Using only Duroc prenatal developing muscle tissue for microarray analysis of daily comparison, we also observed that energy metabolism decreased during myogenic differentiation (te Pas et al., 2005Go). It seems clear that myogenic differentiation and energy metabolism are connected processes; however, a clear explanation for this observation is lacking. From the reported relationship between energy level and muscle fiber hypertrophy (Louis et al., 2004Go), one may suggest that muscle fiber hyperplasia is affected, but data are lacking.

The fatty acid metabolism genes show an opposite expression profile to energy metabolism genes. In early Duroc embryonic tissues, the expression level was greater than in Pietrain. Duroc pigs are known to have more adipose tissue than Pietrain pigs, which are considered the most lean pig breed. The inter- and intramuscular fat depots are especially larger in Duroc pigs (Sellier, 1998Go). During the isolation of the prenatal tissues, we observed that Duroc embryos from d 35 of gestation and onward showed a well-developed subcutaneous fat layer, whereas Pietrain did not. Thus, our results indicate that the fat depot differences are already initiated during early embryogenesis.

Finally, the differences between Duroc and Pietrain prenatal energy metabolism expression profiles seemed to be less in the 91-d-old fetuses than in younger fetuses, with the exception of ATP metabolism. This finding may indicate that at birth, general energy metabolism is more alike for these 2 pig breeds. On the other hand, ATP metabolism, especially, is a major determinant of energy metabolism, which by itself may have a large effect on energy availability in piglets after birth, with effects on survival rate and growth and body composition, including muscle fiber typing characteristics.


    Footnotes
 
1 The authors acknowledge the contributions and helpful discussions of the colleagues in the PorDictor EU project: K. Wimmers (University Bonn, Germany); K. C. Chang and N. de Costa (University Glasgow, UK); J. Merks and B. Harlizius (Institute for Pig Genetics (IPG); Beuningen, The Netherlands; and H. Henne (Züchtungszentrale Deutsches Hybridschwein Gmbh, Lueneburg, Germany). The prenatal pig samples were delivered from IPG and Züchtungszentrale Deutsches Hybridschwein Gmbh, Lueneburg, Germany. The Pietrain prenatal samples were isolated by K. Wimmers. This work was financially supported by an EU grant under Contract No. QLK5-CT-20000-01363 with additional finances from the director of ID-Lelystad. Back

2 These authors contributed equally to this work. Back

3 Corresponding author: roberta.davoli{at}unibo.it

Received for publication January 18, 2005. Accepted for publication August 25, 2005.


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


Chen, C., K. J. Grzegorzewski, S. Barash, Q. Zhao, H. Schneider, Q. Wang, M. Singh, L. Pukac, A. C. Bell, R. Duan, T. Coleman, A. Duttaroy, S. Cheng, J. Hirsch, L. Zhang, Y. Lazard, C. Fischer, M. C. Barber, Z.-D. Ma, Y.-Q. Zhang, P. Reavey, L. Hong, B. Teng, I. Sanyal, S. M. Ruben, O. Blondel, and C. E. Birse. 2003. An integrated functional genomics screening program reveals a role for BMP-9 in glucose homeostasis. Nat. Biotechnol. 21:294–301.[Medline]

Cleveland, W. S. 1974. Robust locally weighted regression and smoothing scatterplots. J. Am. Stat. Assoc. 74:829–836.

Coutinho, L. L., J. Morris, H. L. Marks, R. J. Buhr, and R. Ivarie. 1993. Delayed somite formation in a quail line exhibiting myofiber hyperplasia is accompanied by delayed expression of myogenic regulatory factors and myosin heavy chain. Development 117:563–569.[Abstract]

Davoli, R., L. Fontanesi, P. Zambonelli, D. Bigi, J. Gellin, M. Yerle, J. Milc, S. Braglia, V. Cenci, M. Cagnazzo, and V. Russo. 2002. Isolation of porcine expressed sequence tags for the construction of a first genomic transcript map of the skeletal muscle in pig. Anim Genet. 33:3–18. Erratum in Anim. Genet. 33: 168.[Medline]

Davoli, R., P. Zambonelli, D. Bigi, L. Fontanesi, and V. Russo. 1999. Analysis of expressed sequence tags of porcine skeletal muscle. Gene 233:181–188.[Medline]

Dwyer, C. M., J. M. Fletcher, and N. C. Stickland. 1993. Muscle cellularity and postnatal growth in the pig. J. Anim. Sci. 71:3339–3343.[Abstract]

Dwyer, C. M., and N. C. Stickland. 1991. Sources of variation in myofibre number within and between litters of pigs. Anim. Prod. 52:527–533.

Jones, G. F. 1998. Genetic aspects of domestication, Common breeds and their origin. Pages 17–50 in The Genetics of the Pig. M. F. Rothschild and A. Ruvinsky, ed. CAB Intl., Oxon, UK.

Jorgensen, B., and T. Vestergaard. 1990. Genetics of leg weaknesses in boars at the Danish pig breeding stations. Acta Agric. Scand. 40:59–69.

Louis, M., R. Van Beneden, M. Dehoux, J. P. Thissen, and M. Francaux. 2004. Creatine increases IGF-I and myogenic regulatory factor mRNA in C2C12 cells. FEBS Lett. 557:243–247.[Medline]

Merks, J. W. M. 2000. One century of genetic changes in pigs and the future needs. Pages 8–19 in The Challenge of Genetic Change in Animal Production. W. G. Hill, S. C. Bishop, B. McGuirk, J. C. McKay, G. Simm, and A. J. Webb, ed. Occasional Publication No. 27. Br. Soc. Anim. Sci., Edinburgh, UK.

Olson, E. N. 1990. Myo D family: A paradigm for development? Genes Dev. 4:1454–1461.[Free Full Text]

Olson, E. N., T. J. Brennan, T. Chakraborty, T.-C. Cheng, P. Cserjesi, D. Edmondson, G. James, and L. Li. 1991. Molecular control of myogenesis: Antagonism between growth and differentiation. Mol. Cell. Biochem. 104:7–13.[Medline]

Park, T., S.-G. Yi, S.-H. Kang, S. Y. Lee, Y.-S. Lee, and S. Simon. 2003. Evaluation of normalization methods for microarray data. BMC Bioinformatics 4:33.[Medline]

Parker, L., M. H. P. Seale, and M. Rudnicki. 2003. Looking back to the embryo: Defining transcriptional networks in adult myogenesis. Nat. Rev. Genet. 4:497–507.[Medline]

Pool, M. H., W. Kuurman, B. Hulsegge, L. L. G. Janss, J. M. J. Rebel, and S. van Hemel. 2003. Procedure for standardisation and normalisation of cDNA microarrays. Annu. Mtg. EAAP 2003. Rome, Italy. Abstr. G3.34.

Radonic, A., S. Thulke, I. M. Mackay, O. Landt, W. Siegert, and A. Nitsche. 2004. Guideline to reference gene selection for quantitative real-time PCR. Biochim. Biophys. Res. Comm. 313:856–862.[Medline]

Rehfeldt, C., I. Fiedler, G. Dietl, and K. Ender. 2000. Myogenesis and postnatal skeletal muscle cell growth as influenced by selection. Livest. Prod. Sci. 66:177–188.

Riera, L., M. Obach, A. Navasrro-Sabaté, J. Duran, J. C. Perales, F. Viñals, J. L. Rosa, F. Ventura, and R. Bartrons. 2003. Regulation of ubiquitous 6-phosphofructo-2-kinase by the ubiquitin-proteasome proteolytic pathway during myogenic C2C12 cell differentiation. FEBS Lett. 550:23–29.[Medline]

Sellier, P. 1998. Genetics of meat and carcass traits. Pages 463–510 in The Genetics of the Pig. M. F. Rothschild and A. Ruvinsky, ed. CAB Intl., Oxon, UK.

Soumillion, A., J. H. F. Erkens, J. A. Lenstra, G. Rettenberger, and M. F. W. te Pas. 1997. Genetic variation in the porcine myogenin gene locus. Mamm. Genome 8:564–568.[Medline]

Stickland, N. C., and G. Goldspink. 1973. A possible indicator muscle for the fibre content and growth characteristics of porcine muscle. Anim. Prod. 16:135–146.

te Pas, M. F. W., A. de Wit, J. Priem, M. Cagnazzo, R. Davoli, V. Russo, and M. H. Pool. 2005. Prenatal muscle transcriptome expression profiling in pigs related to myogenesis. J. Muscle Res. Cell Motil. (In press). Available: http://www.springerlink.com. Accessed Aug. 17, 2005.

te Pas, M. F. W., F. L. Harders, A. Soumillion, L. Born, W. Buist, and T. H. E. Meuwissen. 1999a. Genetic variation at the porcine myf-5 gene locus — Lack of association with meat production traits. Mamm. Genome 10:123–127.[Medline]

te Pas, M. F. W., A. Soumillion, F. L. Harders, F. J. Verburg, T. J. van den Bosch, P. Galesloot, and T. H. E. Meuwissen. 1999b. Influences of myogenin genotype on birth weight, growth rate, carcass weight, backfat thickness and lean weight of pigs. J. Anim. Sci. 77:2352–2356.[Abstract/Free Full Text]

te Pas, M. F. W., F. J. Verburg, C. L. M. Gerritsen, and K. H. de Greef. 2000. Messenger ribonucleic acid expression of the MyoD gene family in muscle tissue at slaughter in relation to Selection for porcine growth rate. J. Anim. Sci. 78:69–77.[Abstract/Free Full Text]

Voronina, A. S. 2002. Translation regulation in early development of eukaryotes. Mol. Biol. 36:773–784.

Wagenknecht, D., A. Stratil, H. Bartenschlager, M. van Poucke, L. J. Peelman, I. Majzlik, and H. Geldermann. 2003. Linkage and radiation hybrid mapping of the porcine MEF2D gene to chromosome 4q. Anim. Genet. 34:232–233.[Medline]

Weintraub, H., R. Davis, S. Tapscott, M. Thayer, M. Krause, R. Benezra, T. K. Blackwell, D. Turner, R. Rupp, S. Hollenberg, Y. Zhuang, and A. Lassar. 1991. The myoD gene family: Nodal point during specification of the muscle cell lineage. Science 251:761–766.[Abstract/Free Full Text]

Wigmore, P. M. C., and D. J. Evans. 2002. Molecular and cellular mechanisms involved in the generation of fiber diversity during myogenesis. Int. Rev. Cytol. 216:175–232.[Medline]

Wigmore, P. M. C., and N. C. Stickland. 1983. Muscle development in large and small pig fetuses. J. Anat. 137:235–245.

Yang, Y. H., S. Dudoit, P. Luu, P. Lin, D. M. Lin, V. Peng, J. Nagi, and T. Speed. 2002. Normalization of cDNA microarray data: A robust composite method addressing single and multiple slide systematic variation. Nucleic Acids. Res. 30, No. 4e15. Available: http://nar.oxfordjournals.org/cgi/content/abstract/30/4/e15. Accessed Feb. 24, 2003.


This article has been cited by other articles:


Home page
J ANIM SCIHome page
J. Tong, M. J. Zhu, K. R. Underwood, B. W. Hess, S. P. Ford, and M. Du
AMP-activated protein kinase and adipogenesis in sheep fetal skeletal muscle and 3T3-L1 cells
J Anim Sci, June 1, 2008; 86(6): 1296 - 1305.
[Abstract] [Full Text] [PDF]


Home page
Brief Funct Genomic ProteomicHome page
R. Davoli and S. Braglia
Molecular approaches in pig breeding to improve meat quality
Brief Funct Genomic Proteomic, January 21, 2008; (2008) elm036v1.
[Abstract] [Full Text] [PDF]


Home page
Physiol. GenomicsHome page
S. Ponsuksili, E. Murani, C. Walz, M. Schwerin, and K. Wimmers
Pre- and postnatal hepatic gene expression profiles of two pig breeds differing in body composition: insight into pathways of metabolic regulation
Physiol Genomics, May 11, 2007; 29(3): 267 - 279.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary Data
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 Cagnazzo, M.
Right arrow Articles by Russo, V.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Cagnazzo, M.
Right arrow Articles by Russo, V.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS