|
|
||||||||
Department of Animal Science, Cornell University, Ithaca, NY 14853
2 Correspondence:
314 Morrison Hall (phone: 607-255-2851; fax: 607-255-9829; E-mail:
rgm9{at}cornell.edu).
| Abstract |
|---|
|
|
|---|
Key Words: Gene Expression Muscles Sheep
| Introduction |
|---|
|
|
|---|
Testosterone could exert an effect on muscle growth through the IGF-I axes. Locally produced IGF-I is an important growth regulator acting in an autocrine and paracrine manner (Weimann and Kiess, 1990), but different muscles may possess different IGF-I sensitivities (Boge et al., 1995) and(or) IGF-I synthesis rates (Pfaffl et al., 1998b; Thissen et al., 1994) and, therefore, exhibit different growth rates. Also, testosterone action is mediated by the androgen receptor (AR), which transduces the steroid signal within cells. Different receptor densities, as postulated by Sauerwein and Meyer (1989), may account for relatively higher sensitivity to testicular steroids in the neck muscles. Testosterone could also regulate muscle-specific genes. One of these is myostatin, a gene shown to be a negative regulator of muscle mass.
The objective of this study was to elucidate the possible role of IGF-I, AR, and myostatin genes in differential growth. The approach was to quantify IGF-I, AR, and myostatin mRNA expression in splenius (a sexually dimorphic muscle) and semitendinosus (not a sexually dimorphic muscle) muscles in rams and to assess whether the levels of expression of these genes were different.
| Materials and Methods |
|---|
|
|
|---|
Ribonucleic Acid Isolation
Total RNA was isolated using the guanidine thiocyanate, acid phenol:chloroform procedure of Chomczynski and Sacchi (ToTALLY RNA Total RNA Isolation Kit, Ambion, Austin, TX). Using a hand-held Tissue Tearor (Biospec Products, Bartlesville, OK), 140 mg of muscle tissue and 40 mg of liver tissue were homogenized on ice in denaturation solution. After isopropanol precipitation, RNA pellets extracted from liver were washed with LiCl to remove glycogen. The RNA pellet was resuspended in 40 to 50 µL (muscles) and 80 to 100 µL (liver) of 1 mM sodium citrate, pH 6.4 ± 0.2 (The RNA Storage Solution, Ambion). The concentration of RNA in the final preparations was calculated from the OD260. The integrity of RNA was verified by denaturing agarose gel electrophoresis and ethidium bromide stain.
Testosterone Assay
Testosterone concentrations in plasma were determined in the Cornell Veterinary Medicine Diagnostic Laboratory using a radioimmunoassay method (Total Testosterone Coat-A-Count Kit, Diagnostic Products, Los Angeles, CA).
Ribonuclease Protection Assay
Preparation of the AR, Myostatin, GAPDH-Labeled Antisense and Unlabeled Sense Riboprobes.
From 500 ng of total ovine muscle RNA, cDNA was generated using a gene-specific primer (Table 1
) in the reverse transcription (Display THERMO-RT Kit, Display Systems Biotech, Boston, MA). Total RNA, 1x RT buffer, 0.5 mM dNTP (each of dATP, dGTP, dCTP, dTTP), 1 µM reverse primer, and Display Terminator Mix were incubated for 40 min at 42°C followed by 10 min at 65°C. Specific regions were amplified from 1 µL of reverse trancription reaction in a 10-µL PCR amplification reaction (Table 1
) using the Idaho Technologies Rapid Cycler (Idaho Falls, ID).
|
|
F' competent cells (Stratagene, La Jolla, CA). Positive transformants were selected on the basis of ampicillin resistance and a selection of these were grown in LB broth containing ampicillin (50 µg/µL). The ovine cDNA structures for AR, myostatin, and GAPDH and the orientations in the plasmid were confirmed by sequencing. The myostatin and GAPDH PCR products were found to be in reverse sequence in the vector. The plasmids were named pCR(myostatin) and pCR(GAPDH) to designate cDNA templates for ovine myostatin, GAPDH, and pCR(ARs) and pCR(ARas) to designate cDNA templates for ovine AR sense and AR antisense (because the pCR 2.1 had only the T7 promoter, two plasmids were required, one for each orientation), respectively. Large-scale plasmid DNA purification was carried out using a Qiagen Plasmid Kit (Qiagen, Valencia, CA). The vectors pCR(myostatin), pCR(ARs), and pCR(ARas) were linearized by BamHI, BamHI and HincII digestion, respectively, and the cloned ovine cDNA were transcribed with T7 RNA polymerase at 37°C for 1 h using the reagents supplied with the MaxiScript Kit (Ambion) to produce sense myostatin, sense AR, and antisense AR, respectively. The vectors pCR(myostatin) and pCR(GAPDH) were linearized by DraI and EcoO109I digestion, respectively, and the cloned ovine cDNA were transcribed with SP6 RNA polymerase (Ambion) at 37C for 1 h to produce antisense myostatin and antisense GAPDH. At the end of the reaction, 1 µL of RNase-free DNase I (10 U) was added and the reaction was incubated for 15 min at 37°C. Subsequently, all sense and antisense RNA transcripts were gel-purified and the concentration was calculated from the OD260 after gel elution. The antisense RNA transcripts were labeled with Psoralen-Biotin using the BrightStar Psoralen-Biotin Nonisotopic Labeling Kit (Ambion). Due to the abundance of GAPDH, the probe concentration was adjusted to generate a signal in the same range with the myostatin and AR signal. This adjustment was accomplished by mixing biotin-labeled and unlabeled GAPDH probe at a ratio of 1:25.
Ribonuclease Protection Assay.
Forty micrograms of splenius and semitendinosus muscle total RNA were coprecipitated with 1 ng of myostatin, AR, and GAPDH riboprobes and hybridization was performed at 42°C overnight using the protocol and reagents supplied in the RPA III Kit (Ambion) as described in the standard procedure. The next day, the mixture was treated with an RNase cocktail (0.04 U of RNase A and 1.5 U of RNase T1) to degrade single-stranded RNA. The protected fragments were precipitated, the pellets were washed once with 75% ethanol, dried, resuspended in 10 µL of gel-loading buffer, denatured at 95°C for 5 min, loaded onto a 6% polyacrylamide/8M urea denaturing gel and run at 200V for approximately 1 h. The RNA was then electrophoretically transferred to a positively charged nylon membrane (Bright Star-Plus, Ambion) using a mini-gel blotting apparatus (BioRad, Hercules, CA) and 1x TBE as transfer buffer. The transfer was carried out for 1 h at 100 mA, and the membrane was UV-irradiated to crosslink the RNA to the nylon. The BrightStar BioDetection Kit (Ambion) was used for chemiluminescent detection of the protected fragments. At the end of the procedure, the membrane was exposed to x-ray film (Kodak) overnight at room temperature and for shorter periods (between 15 and 30 min) the following day.
Quantitative Analysis of Myostatin and AR mRNA.
Known amounts (60, 40, 27, 18, 12, and 8 pg) of in vitro synthesized myostatin and AR sense RNA were hybridized with an excess of labeled antisense probe (1 ng) in RPA to construct standard curves. The values were plotted as the log of the amount of sense RNA vs the intensity of the protected band in the standard curve. The absolute amount of the protected RNA species in the muscle sample RNA was determined by comparing the intensity of the band to the standard curve, after correcting for the size difference and for the RNA loading difference. The protected fragments in the experimental samples differed in size from the protected fragments in the standard curve; therefore, the difference in the incorporation of biotin had to be taken into account. The following equation was used to correct for size difference:
![]() |
To correct for variations in RNA loading, the band densities (corrected for size differences) were multiplied by the ratio between the strongest GAPDH band (considered to correspond to the 40 µg RNA loaded into the gel) and the GAPDH band from the respective lane.
Reverse Transcription-Polymerase Chain Reaction
Construction of a Heterologous Competitor.
The relationship between the initial amount A of target mRNA present in the tissue and the amount Yn of DNA produced after n PCR cycles can be expressed as Yn = A (1 + E)n, where E is the amplification efficiency of one reaction cycle (Chelly et al., 1988). Small variations in reaction efficiency, therefore, translate into large differences in the amount of RT-PCR product generated after n cycles. These limitations in quantitative analyses can be overcome by parallel co-amplification of the native mRNA with known amounts of an internal standard cRNA. The amplification efficiency should affect both templates similarly. Several designs have been used in quantitative RT-PCR to obtain an internal standard cRNA that has amplification efficiency identical to that of the native mRNA template and is easily distinguishable from it. One common practice is the use of a PCR MIMIC, a heterologous internal standard that "mimics" the primer binding and amplification characteristics of the target (Gaudette and Crain, 1991; Siebert and Larrick, 1993), obtained with the addition of the specific primers to gene products such as ubiquitously expressed "housekeeping" genes such as GAPDH (Murphy et al., 1990; Gaudette and Crain, 1991).
To quantify ovine muscle and liver IGF-I mRNA abundance by RT-PCR, we constructed a heterologous competitor (MIMIC) based on the ovine GAPDH sequence obtained by RT-PCR using 20 nt primers homologous to ovine GAPDH sequence to which an additional 20 nt were appended on the 5' ends. The additional 20 nucleotides corresponded to the IGF-I forward and reverse primers, designed based on the ovine IGF-I mRNA sequence (GenBank accession number M31734). These composite primers (Table 1
) amplified a 325-bp product, including the 20 bp for IGF-I at both ends. Reverse transcription and amplification techniques were applied to produce first-strand cDNA. Using the Display THERMO-RT Kit (Display Systems Biotech, Boston, MA), 500 µg of total ovine muscle RNA was amplified with an IGFGAPDH reverse primer in the reverse transcription. One microliter of cDNA from the reverse transcription reaction was used in the PCR amplification reaction using the Idaho Technologies Rapid Cycler (Idaho Falls, ID) and the two composite primers (IGFGAPDH forward primer and IGFGAPDH reverse primer). The 325-bp GAPDH MIMIC PCR product was cloned into pCR II using a TA cloning kit (Invitrogen, Carlsbad, CA). The plasmid was subsequently transfected into INV
F' competent cells (Stratagene, La Jolla, CA). Positive transformants were selected on the basis of ampicillin resistance and, from all positive colonies, 10 were grown in LB broth containing ampicillin (50 µg/µL). The orientation in the plasmid was checked by digestion with AluI restriction enzyme, which has 27 restriction sites in the plasmid and only one in the GAPDH MIMIC. Large-scale plasmid DNA purification was carried out using the Qiagen Plasmid Kit (Valencia, CA). The vector was linearized by EcoRV digestion, and the cloned ovine cDNA was transcribed with SP6 RNA polymerase at 37°C for 1 h using the reagents supplied with the MaxiScript kit (Ambion) to create a 422-bp template containing the 325-bp MIMIC and 97 bp of the plasmid. At the end of the reaction, 1 µL of RNase-free DNase I (10 U) was added and the reaction was incubated for 15 min at 37°C. The concentration of RNA MIMIC was determined using spectrophotometry.
Establishment and Validation of the Quantitative RT-PCR.
To ensure a parallel start in all individual reactions and to increase specificity, yield, and precision of the PCR, Zombie-Taq (Clontech, Palo Alto, CA) DNA polymerase was used. In the Zombie-Taq DNA polymerase, the polymerase activity is inhibited by the Taq antibody until the first denaturation step. To quantify the native IGF-I mRNA in different tissues, a preliminary estimation of the IGF-I cRNA start-molecule concentration range to be used for individual tissues was required. Initially, we titrated against a broad range of dilutions to obtain a rough estimate of the amount of RNA present. This was performed by six titration steps ranging from 2,500 to 0.025 cRNA start-attomoles in a 1:10 dilution scheme added to 0.5 µg of total tissue RNA. We then performed a finer titration over a narrower range of MIMIC (24.1 to 0.1 attomoles for the muscle and 217.2 to 0.9 attomoles for liver in a 1:3 dilution scheme).
Competitive RT-PCR.
For each of the three tissues, six RT reactions were performed using 596 ng of total RNA and fixed concentrations of the MIMIC IGF-I cRNA (217.17, 72.39, 24.13, 8.04, 2.68, and 0.89 attomoles for liver and 24.13, 8.04, 2.68, 0.89, 0.30, and 0.1 attomoles for muscle). First-strand cDNA of target and MIMIC were generated using the Omniscript RT Kit (Qiagen, Valencia, CA) in an Amplitron (Barnstead/Thermolyne, Dubuque, IA) thermocycler. The RNA, 1x RT buffer, 0.5 mM of each dNTP, 1 µM reverse primer, and 2 U Omniscript Reverse Transcriptase were incubated for 60 min at 37°C. At the end of the reaction, the Omniscript Reverse Transcriptase was inactivated by heating the reaction mixture to 93°C for 5 min followed by rapid cooling on ice. The following PCR was performed in a Rapid Cycler (Idaho Technologies) with 1 µL cDNA; Idaho Technologies buffer (50 mM Tris, pH 8.3; 0.25 mM crystalline BSA; 2 mM MgCl2; and sucrose), 0.4 U Taq DNA polymerase (Gibco BRL), 0.2 mM dNTPs (Promega), and 0.5 µM each primer. We used 30 cycles (0 s at 94°C, 0 s at 60°C, 15 s at 75°C) followed by an additional 30 s at 75°C for complete amplification of all PCR products. For the Rapid Cycler, 0 s means that the denaturing and annealing temperatures were reached without setting a hold time, because both denaturing and annealing occur almost instantaneously after the sample has reached the appropriate temperature. In each experiment, a negative control, with water replacing template, was included.
Quantification of PCR Products.
The 10-µL PCR product was electrophoresed on a 4% NuSieve agarose gel (East Rutherford, NJ), which was then stained in a solution of 0.1 µg/mL ethidium bromide (Figure 2a
). The gel was scanned using a gel documentation and analysis system (Alpha Innotech, San Leandro, CA) under UV light. The area of optical densities of the target and the MIMIC peaks were quantified using the Alpha Innotech software. Because the size of the MIMIC cDNA was smaller than the size of the target cDNA, the decrease in the incorporation of ethidium bromide had to be taken into account. The optical density associated with the 325-bp competitor (MIMIC) product was corrected by a multiplication factor corresponding to the base-pair ratio of both fragments as follows:
|
![]() |
The regression of log10 of ratio of the corrected MIMIC band densities to the target band densities (Y variable) against log10 of concentration of the MIMIC (X variable) was performed. As shown in Figure 2b
, the association was linear. When the ratio of MIMIC to target was 1 (band densities were equal) the amounts of the target and the MIMIC were equal. To determine the amount of target IGF-I in the sample, we calculated the value of X for a Y value of 0 (log10 1 = 0).
Statistical Analyses
Our main objective was to compare the gene expression in two specific muscles. Because the two muscles were sampled from each individual, the gene expression of the two muscles represented paired observations with respect to the individual. A new variable was created as the difference between the amount of mRNA expressed in the two muscles for the three genes, Di = (splenius - semitendinosus)i, where i = 1, 2, ..., 11 (animals in the study).
To determine whether age and weight of the animals affected the difference in gene expression in the semitendinosus and splenius muscles, the following regression model was used:
![]() |
where Di = difference in attomoles/µg RNA between splenius and semitendinosus muscle for the ith animal, Ai = age of the ith animal, Wi = weight of the ith animal, and (AW)i = age x weight of the ith animal.
Prior to a formal statistical analysis, we plotted the data to assess the type and strength of the relationship. Significant explanatory variables were chosen by a step-down procedure. First, the crossproduct between the two explanatory variables was evaluated and, if not significant, was removed from the model. Second, each explanatory variable was evaluated and, if not significant, was removed from the model.
To test whether there was a difference in gene expression between the two muscles, a hypothesis test was performed. The null hypothesis was Ho: d = 0 and the alternative hypothesis was Ha: d > 0 for IGF-I and AR and Ha: d < 0 for myostatin.
The assumptions of normality and equal variance were checked with a normal probability plot and a plot of the absolute values of residuals against the predicted values.
| Results and Discussion |
|---|
|
|
|---|
Testosterone
Plasma testosterone concentrations ranged from 0.16 ng/mL to 4.80 ng/mL. The concentrations of testosterone were not correlated with the expression of the three genes investigated. All testosterone concentrations were within the physiological range for rams at this age (Arnold et al., 1997).
Insulin-like Growth Factor I mRNA Expression
Insulin-like growth factor I is a ubiquitous peptide known to stimulate myoblast survival, proliferation, and differentiation (Florini et al., 1995; Stewart and Rotwein, 1996) and it is a potent anabolic agent in regulation of skeletal muscle protein metabolism (Harper et al., 1987; Roeder et al., 1988; Frost et al., 1997).
The IGF-I gene is expressed in many tissues, but liver, and to a lesser extent bone, are the primary sources of circulating IGF-I (Jones and Clemmons, 1995). Several studies have shown that IGF-I can be produced by skeletal muscle, although it is only found at low levels in normal, mature animals (Isgaard, 1992; McGuire et al., 1992). Knock-out mice with a liver-specific IGF-I deletion (Yakar et al., 1999) showed a 75% reduction in IGF-I circulating levels but had normal postnatal and pubertal growth, suggesting that autocrine/paracrine IGF-I production is sufficient for normal growth and development. In all extracellular tissues, IGF-I is bound to a family of six high-affinity IGF-binding proteins (IGFBP), designated IGFBP-1 to -6 (Jones and Clemmons, 1995). The IGFBP function not only as carrier proteins for IGF-I in circulation, protecting it from degradation and transporting it to specific tissues, but also as modulators of IGF-I action.
Insulin-like growth factor-I mRNA expression was determined using a competitive RT-PCR, and the mean difference between the splenius and semitendinosus muscles is shown in Table 2
. The variation in IGF-I mRNA levels was related neither to age nor to weight of the animals within the ranges covered in the study. The splenius muscle had a higher concentration (P < 0.05) of IGF-I mRNA than the semitendinosus muscle (Table 2
). The mean difference between the two muscles was statistically significant (P < 0.05). This result is consistent with a study of IGF-I gene expression in cattle (Pfaffl et al., 1998a) in which IGF-I mRNA expression was measured in two muscles selected because of their overproportional (splenius) and underproportional (gastrocnemius) growth response to testicular steroids. In bulls, higher IGF-I mRNA concentration was found in the splenius than in the gastrocnemius muscle. Our data support the hypothesis that local differences in IGF-I expression might be one of the mediators of the differential growth of these individual muscles in intact males. However, the concomitant regulation of all other factors of the IGF system, such as the binding proteins and receptors, needs to be considered along with the observed differential regulation of local IGF-I mRNA.
|
Androgen Receptor mRNA Expression
The AR is a member of the superfamily of nuclear transcription factors that mediate the action of steroid hormones. When activated by ligand binding, these transcription factors bind to specific DNA sequences on target genes, referred to as hormone response elements, and regulate the transcriptional activity of those genes. It is now recognized that unliganded androgen receptors are primarily localized to the cell nucleus (Prins, 2000), and androgens diffuse into the nucleus, where they bind to available receptors.
Testosterone stimulates muscle growth by affecting the rate of protein synthesis, protein breakdown, and the net gain or loss of muscle protein (Rooyackers and Nair, 1997). When testosterone increases muscle protein synthesis, intramuscular mRNA concentrations of IGF-I are increased and concentrations of the inhibitory IGF binding protein 4 are decreased (Urban et al., 1995), and locally produced IGF-I subsequently stimulates muscle protein synthesis and decreases protein breakdown. The response to testosterone differs among muscle groups, and this differential response may be explained by the variation of AR number among skeletal muscles (Krieg et al., 1977). Thus, sexual dimorphism can be explained partly by higher androgen sensitivities in muscles with pronounced growth under androgen stimulation (Sauerwein and Meyer, 1989). A recent study in cattle showed that androgen receptor mRNA concentrations in muscles with different fiber type compositions and growth stimuli is positively related to the individual growth patterns (Brandstetter et al., 2000).
Androgen receptor gene expression in the two muscles showing differential growth patterns was measured by a ribonuclease protection assay (Figure 3
) using GAPDH to normalize AR data for differences in RNA loading. Age and weight did not contribute to the variation in AR mRNA levels. Androgen receptor mRNA levels per unit of total RNA measured in splenius were 37% greater than those measured in semitendinosus (P < 0.05, Table 2
). The increased responsiveness of neck muscles to testosterone could be correlated with the AR level, because changes in the number of receptors are a point of modulation of sensitivity to hormones. Divergent densities of AR in individual muscles have been reported previously in cattle (Sauerwein and Meyer, 1989), and data on the developmental regulation of AR mRNA expression in three bovine muscles that differed in muscle fiber composition, metabolic activity, and growth pattern showed a relationship between AR mRNA concentration and differential growth (Brandstetter et al., 2000). The higher AR mRNA concentrations in splenius muscle shown in this experiment are in agreement with these studies and help to explain the pronounced muscle growth in the neck of maturing rams.
|
The results from myostatin null mice and from double-muscled cattle breeds indicate that lack of myostatin activity during embryonic development allows growth of muscle tissue. There is also increasing evidence that myostatin plays a role in postnatal muscle growth. Myostatin levels are increased in the serum of HIV-infected patients with wasting (Gonzalez-Cadavid, 1998), elderly men and women with sarcopenia (Yarashesky, 1999), and in skeletal muscle in a rat model of aging (Mallidis, 1999) and during hindlimb unloading (Baldwin, 1996).
Recently, myostatin was shown to function by controlling the proliferation of muscle precursor cells (Thomas et al., 2000). Myostatin promotes an increase in p21 expression, a cyclin-dependent kinase inhibitor, and a decrease in Cdk2 protein and activity, thus resulting in an accumulation of hypophosporylated Rb protein. This, in turn, leads to the arrest of myoblasts in the G1-phase of the cell cycle. Thus, myoblast number and, hence, fiber number, following differentiation, is regulated (limited), suggesting that the increased number of myofibers in cattle and mice with increased muscle mass is a result of deregulated myoblast proliferation caused by the absence of functional myostatin (Thomas et al., 2000).
The RPA measured the presence of myostatin mRNA in both skeletal muscles (Figure 3
). No statistically significant difference in myostatin mRNA concentration was found between the splenius and semitendinosus muscles and age and weight of the animals had no effect. We were expecting myostatin gene expression to be lower in the splenius muscle compared to semitendinosus to help explain the increased muscle mass in splenius. Our study did not implicate myostatin in the sexual dimorphism of muscle growth. Nevertheless, myostatin could still play a role in this process, and a study comparing the myostatin gene expression in rams and wethers or in rams before and after puberty might address this issue.
| Implications |
|---|
|
|
|---|
| Footnotes |
|---|
Received for publication October 3, 2001. Accepted for publication February 21, 2002.
| Literature Cited |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
R. G. Mateescu and M. L. Thonney Effect of testosterone on insulin-like growth factor-I, androgen receptor, and myostatin gene expression in splenius and semitendinosus muscles in sheep J Anim Sci, April 1, 2005; 83(4): 803 - 809. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. Sinha-Hikim, W. E. Taylor, N. F. Gonzalez-Cadavid, W. Zheng, and S. Bhasin Androgen Receptor in Human Skeletal Muscle and Cultured Muscle Satellite Cells: Up-Regulation by Androgen Treatment J. Clin. Endocrinol. Metab., October 1, 2004; 89(10): 5245 - 5255. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |