|
Advertisement | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 282, Issue 45, 33086-33097, November 9, 2007
Determination of Dual Effects of Parathyroid Hormone on Skeletal Gene Expression in Vivo by Microarray and Network Analysis*
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ABSTRACT |
|---|
|
|
|---|
1,000), but differentially, in both regimes. PTH-(1-31) regulated a similar number of genes in the intermittent regimen but fewer in the continuous regimen, consistent with its less potent catabolic effect. PTH-(3-34) regulated very few genes in both regimes, which suggests the protein kinase C pathway plays a limited role in mediating the dual effects of PTH, whereas the cAMP-dependent protein kinase A pathway appears to predominate. In the intermittent treatment, many genes encoding signaling mediators, transcription factors, cytokines, and proteases/protease inhibitors are regulated rapidly and cyclically with each PTH injection; genes associated with skeletal development show a slowly accruing pattern of expression. With continuous treatment, some genes are regulated from 6 h, and the mRNA levels are sustained with a longer infusion, whereas others show a kinetic decrease and then increase later. Significant up-regulation of genes stimulating osteoclastogenesis in the anabolic regime suggests a provocative and paradoxical theme for the anabolic effect of PTH that a full anabolic response requires a transient up-regulation of genes classically associated with a resorptive response. Ingenuity pathway analysis was performed on the microarray data. A novel signaling network was established that is differentially regulated in the two PTH treatment regimes. Key regulators are suggested to be AREG, CCL2, WNT4, and cAMP-responsive element modulator. | INTRODUCTION |
|---|
|
|
|---|
Bone is a highly specialized form of connective tissue and dynamic organ in all higher vertebrates; it undergoes continuous regeneration. Bone remodeling occurs at discrete sites within the skeleton and proceeds in an orderly fashion, with coupling of osteoclastic resorption and osteoblastic formation (3). This physiological process is coordinated and tightly regulated by local and endocrine factors to ensure that there is a balance between bone formation and bone resorption (4). PTH is known to regulate bone metabolism directly by its receptor on the osteoblast plasma membrane; binding of PTH to its receptor causes release of cytokines and growth factors from the osteoblast. As a result, it indirectly affects osteoclasts and other bone cells in a paracrine fashion (5, 6). Upon binding to its receptor, early signaling events associated with transient increases in PTH levels are sufficient to induce the anabolic response. However, late events associated with continuously elevated PTH lead to the catabolic response and bone disease. Therefore, understanding the molecular events of these two administration regimens, intermittent versus continuous treatment, is essential for determining the signal transduction pathways that mediate PTH-induced parathyroid bone disease versus bone formation. A detailed characterization of the genome-wide dynamic regulation of gene expression by PTH is a basis to investigate the molecular mechanism of its dual actions, anabolic versus catabolic, in response to intermittent or continuous treatment of this hormone in bone.
There have been many reports describing genes that are regulated by PTH in osteoblastic cells; however, at the time this project was launched, there were no established data for the entire differentially expressed gene profiles from PTH administration in anabolic versus catabolic animal models. We think the description of the genome-wide dynamic regulation of gene expression by PTH is critical to investigate the molecular mechanism of its dual actions in bone. The goal of this study was to utilize well established in vivo PTH anabolic versus catabolic treatment regimens to identify genes that mediate bone anabolic versus catabolic effects of the hormone.
Furthermore, the involvement of the PKC and the PKA pathway in mediating the anabolic effects of PTH is controversial, and much less is known about the involvement of the PKC pathway in mediating the bone-resorbing effects of PTH. Based on previous findings, unlike the full activity PTH peptide PTH-(1-34), either amino-terminally or carboxyl-terminally truncated PTH analogs selectively activate cAMP and intracellular calcium signaling pathways or PKC and the intracellular calcium system. When PTH-(3-34) binds to PTH1R, it activates the PKC/calcium pathways but does not elicit a cAMP response (7), whereas PTH-(1-31) activates cAMP/PKA pathway but has no effect on PKC signaling (8). Comparison of gene expression profiles associated with the different PTH peptides in both the intermittent and continuous administration regimes will provide important information to understand the relationship between specific PTH-activated signaling pathways and the dual effects of the hormone on bone metabolism.
DNA microarray technology is a relatively new technology that allows researchers to monitor the expression level of thousands of genes simultaneously in a single experiment. Previously we have successfully applied this technique to identify the PTH-regulated gene profile in osteoblastic cells in vitro (9). However, in vitro data are not sufficient to comprehensively understand the gene profile regulated by PTH, especially when comparing the anabolic and catabolic effects of the hormone, because the in vitro system cannot represent the dual effects of PTH. With the complete rat genome sequence on Affymetrix GeneChips, we can now monitor the gene expression profile change induced by PTH treatment for all genes in rats, the best known physiological model to recapitulate the dual effects of PTH. The skeletal anabolic response to PTH in rats is a universal phenomenon regardless of age, strain, and gender (10-12), whereas different strains of mice require different dosages of PTH to show its anabolic responses (13-16).
More excitingly, a structured network knowledge-based approach to analyze genome-wide transcriptional responses in the context of known functional interrelationships among proteins, small molecules, and phenotypes has been established. This approach has been successfully used to analyze changes in blood leukocyte gene expression patterns in human subjects receiving an inflammatory stimulus (17) and other systems (18, 19). In this study, we applied this knowledge-based approach analysis, and for the first time, we demonstrate the interactions and network of signaling elicited by anabolic and catabolic treatment regimens of PTH peptides.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
Animals
Sprague-Dawley rats (3-month-old female, weight about 250 g) were purchased from Taconic Farms, Inc. (Germantown, NY). For anabolic models, PTH-(1-34), -(1-31), and -(3-34) peptides dissolved in vehicle (0.9% saline with 1% rat heat-in-activated serum) in a final volume of 200 µl were injected subcutaneously daily for the indicated times at 8 µg/100 g/day to 3-month-old female Sprague-Dawley rats for up to 14 days (10-12, 20, 22; at least four animals per group). To permit measurements of bone apposition rates and the calculation of dynamic parameters in 14-day-treated rats, animals received two subcutaneous injections of calcein, at 9 and 2 days before being killed (2 mg/100 g).
For catabolic models, PTH-(1-34), -(1-31), and -(3-34) dissolved in vehicle (0.9% saline with 1% rat heat-inactivated serum) in a final volume of 200 µl were continuously infused into animals at a nominal pumping rate of 1 µl/h by Alzet osmotic pumps (DURECT Corp.) implanted subcutaneously; continuous infusion of PTH-(1-34), -(1-31), and -(3-34) peptides was at 4 µg/100 g/day for the indicated times to 3-month-old female Sprague-Dawley rats (12, 22; at least four animals per group). The animal protocols were approved by the Robert Wood Johnson Medical School Institutional Animal Care and Use Committee.
Tissue Collection and Processing
At the time of implantation of osmotic pumps into rats with anesthesia, 1 ml of blood was collected by a glass capillary tube being inserted at the medial canthus to gain access to the rat orbital venous plexus, the area that offers a readily accessible site for blood collection. Animals were euthanized using CO2 right after infusion ceased or at the indicated time after injection. Long bones, both pairs of femora and tibiae, were harvested. After the animals were unconscious, a cardiac puncture was performed to collect blood (4 ml) to measure serum calcium, phosphate, exogenous hPTH-(1-34) levels, and markers of bone turnover. For histomorphometry and peripheral quantitative computed tomography, femora or tibiae were preserved in 70% ethanol until they were processed for further fixation and plastic embedding in methyl-methacrylate resin, cutting sections, and Goldner's Trichrome staining for microscopic analysis of the various histomorphometric parameters.
Blood was left at room temperature for at least 30 min before centrifuging at 200 x g for 10 min to separate sera. Total calcium was determined by a kit from Sigma. Serum phosphate levels were determined using an automatic hemogas analyzer by the Department of Pathology and Laboratory Medicine at UMDNJ-Robert Wood Johnson Medical School.
RNA Isolation
Femora were cut, and all connective tissue, including periosteum, completely removed. The distal epiphysis, including the growth plate, was removed, and a subjacent 3-mm wide band of the metaphyseal trabecular primary spongiosa was dissected (20) and kept in RNAlater (Qiagen) at 4 °C until mRNA analyses were performed. Total RNA was extracted by homogenization in TRI Reagent® (Sigma) using a PT 10-35 Polytron homogenizer (Kinematica GA, Switzerland) as recommended by the manufacturer. Isolated RNA was quantitated using spectrophotometry by measuring the absorbance at 260 nm, and the 260/280 nm ratio was calculated, and RNA was kept only when this ratio was 1.9-2.0 to ensure the absence of protein contamination and limited RNA degradation. The integrity of sample RNAs (collected at 1 h after the last injection of 14-day intermittent treatment or right after 14-day infusion ceased) used for microarray analysis was confirmed by capillary electrophoresis resolving 18 S and 28 S ribosomal RNA on an Agilent Technologies 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA).
Enzyme-linked Immunosorbent Assay (ELISA)
The exogenous hPTH-(1-34) levels in sera of rats treated with 14-day infusion were examined by ELISA kit (Immutopics International, San Clemente, CA). This kit specifically detects bioactive hPTH-(1-34). RatLaps ELISA was performed using sera collected from rats in the 14-day infusion groups to measure collagen type I C-telopeptide degradation products (CTX-I) using a specific monoclonal antibody in a competitive ELISA form (Nordic Bioscience, Diagnostics A/S Herlev, Denmark) according to the manufacturer's instructions.
Histology
Longitudinal undemineralized sections (5 µm thick) were cut from methyl methacrylate plastic-embedded blocks of frontal sections of tibiae or femora, using a Riechert-Jung microtome. These were mounted on gelatin-coated 25 x 75 frosted slides and stained with Goldner's Trichrome stain for the static measurements, and additional sections were cut at 10 µm and left unstained for dynamic (fluorescent) measurements of calcein injections.
A region of interest was selected that was exactly 1 mm distal to the growth plate, and extending 2 mm downward (thereby avoiding the primary spongiosa) through the metaphysis of the tibia or femora. The region of interest is the same size in every animal, regardless of the shape of the section. Standard bone histomorphometry was performed by the methods of Parfitt et al. (21). Image analysis software (SPOT Advanced; Diagnostic Instruments, Inc.) was used for tibial histological analysis. Femoral histological analysis was performed by the Department of Pathology at the University of Alabama, Birmingham.
Microarray Analysis
Microarray analyses were performed utilizing 15 µg of total RNA. It should be noted that our analyses are focused on the cancellous and endosteal cortical bone compartments because the periosteum was removed. Affymetrix Rat Genome RAE 230A chips, including probe sets for 5,399 known rat genes and 10,467 ESTs, were used for all intermittent treatment regimen experiments (four replicates in vehicle and PTH-(1-34)-treated samples, three replicates in PTH-(1-31)-treated and PTH-(3-34)-treated samples); RAE 230 2.0. chips (contains over 31,000 probe sets representing more than 30,200 transcripts and variants, including 5,399 known rat genes, the same genes as 230A, and 25,643 ESTs, an upgrade from 230A, and the standard chip available at the time of the continuous experiment) were used for all 14-day continuous treatment regimen experiments (four replicates with all PTH peptides and vehicle-treated samples). Gene expression values were determined using GCOS software (version 1.4 from Affymetrix) with global scaling normalization method (target value 150). The gene expression data after GCOS were further analyzed using GeneSpring software from Agilent or ArrayTools software from the Biometric Research Branch of the NCI, National Institutes of Health.
Detecting Significantly Changed Genes—We used the following two methods to determine significantly regulated genes by PTH peptide/administration categories, 1) For the ANOVA test, we conducted ANOVA analysis using GeneSpring software from Agilent with a parametric test, variances not assumed equal (Welch ANOVA). Then we performed multiple testing corrections using the Benjamini and Hochberg False Discovery Rate. 2) For class comparison, we used ArrayTools developed by NCI Biometric Research Group to conduct class comparison to determine significantly regulated genes between different PTH peptide/administration groups.
Systematic Analysis of Differentially Expressed Genes in Different Peptide Treatments—Because of the fact that the two sets of experiments were conducted using different Affymetrix chipsets, we used the following steps to determine the significantly regulated genes in each experiment. First we identified the significantly regulated gene probe sets compared with control in each treatment regime, intermittently or continuously, using the class comparison method provided in ArrayTools at a parametric p value cutoff of 0.001. Then we imported the significantly regulated genes into an Access data base. We queried the data base to obtain the commonly regulated genes in both treatment regimens and genes that are uniquely affected by a single treatment regime.
Real Time RT-PCR
Two-step real time RT-PCR was performed to confirm the differential expression results obtained by microarray experiments and the expression levels at indicated time points during intermittent and continuous treatments as indicated. As described earlier (9), total RNA was isolated using TRI Reagent® (Sigma) followed by purification with the RNeasy kit (Qiagen). A TaqMan® reverse transcription kit (Applied Biosystems) was used to reverse transcribe mRNA into cDNA. Following this, PCR was performed on Opticon (MJ Research) using a SYBR® Green PCR core kit (Applied Biosystems). The primers used for the RT-PCR are summarized in supplemental Table 5.
-Actin was used as an internal control.
| RESULTS |
|---|
|
|
|---|
Histomorphometric analyses of the proximal femoral metaphyses confirmed the expected morphological changes that occur with PTH-(1-34) intermittent treatment (22). Substantial increases in bone formation rate, as measured by the labeled surface to bone surface ratio, and in osteoblast parameters were confirmed in intermittent treatment groups relative to control animals (Table 1). Consistently, 14-day intermittent administration of PTH-(1-34) increased tibial total and trabecular bone mineral density level (Table 3). PTH-(1-31) also showed anabolic effects by significantly increasing osteoblast parameters as well as bone formation rate.
|
|
We measured the hPTH-(1-34) level in continuously treated animals three times (at the beginning, the middle, and at the end of infusion) to ensure the performance of the osmotic pumps (Fig. 1C) and found it was elevated at both times. We confirmed the catabolic effect of PTH-(1-34) in the continuous regime by detecting an increase in degraded type I collagen levels (Fig. 1D) and a significant decrease in total bone mineral density level in tibiae (Table 3). Significant increases in osteoclast parameters in tibiae were found in the groups infused with PTH-(1-34) and PTH-(1-31) (Table 2). Although a significant decrease in bone volume only occurred in the group infused with PTH-(1-34), significant increases in bone resorption, indicated by the elevation of the carboxyl-terminal product of type I collagen, were found in both PTH-(1-34) and -(1-31) but not in PTH-(3-34)-infused groups (Fig. 1D).
|
From the bone sections, we can easily observe the strong similarity in effects of PTH-(1-34) and PTH-(1-31) in both treatment regimes as follows: increasing trabecular bone volume both in tibiae and femora with intermittent injection (Table 1; Fig. 2, A and B; Table 3); strong peritrabecular fibrosis after infusion (Fig. 2C). PTH-(3-34), unlike the other two PTH peptides, failed to induce significant bone responses in either treatment protocol (Fig. 2). Overall, these results suggest that the PKA pathway is the predominant pathway for PTH action in bone, whereas the PKC pathway has a limited role.
|
Similarly, we identified 3088 genes showing significant gene expression change with p values less than 0.05 between any different PTH peptide and vehicle-treated groups in the continuous treatment experiment. We further trimmed the gene list by filtering for genes with p values less than 0.01 and only included genes with Presence call in at least 3 of 16 chip results. There are 754 genes that passed the filtering criteria to conduct hierarchical clustering with continuously treated samples (Fig. 3B). Our data indicate that continuous treatment with PTH-(3-34) (CP334) induced very limited gene expression changes that cluster together with the controls. CP134 induced the most significant gene expression changes. CP131 also induced some gene expression changes and clusters together with CP134. However, the four CP131 samples were spread on either side of the CP134 samples indicating the variation of gene expression values in CP131-treated animals is relatively high.
Systematic Analysis of Differentially Expressed Genes in Different Peptide Treatments—Because of the fact that the two sets of experiments were conducted using different Affymetrix chipsets, we adopted the following steps to determine the significantly regulated genes in each experiment:
Using the class comparison method provided in ArrayTools, we determined PTH-(1-34) intermittent treatment induced 1154 unique probe sets at a parametric p value cutoff of at least 0.001 versus controls. Among them, 345 are known or predicted genes (according to NCBI RefSeq record) with a fold change of at least 1.5-fold uniquely regulated by intermittent treatment with PTH-(1-34) (supplemental Tables 1 and 6); similarly, PTH-(1-34) continuous treatment induced 1342 significantly changed probe sets at a parametric p value cutoff of less than 0.001. There are 462 known or predicted genes uniquely regulated by continuous treatment with PTH-(1-34) by at least 1.5-fold (supplemental Tables 2 and 7). The number of commonly regulated probe sets with both intermittent and continuous PTH-(1-34) treatments is 675. For the commonly regulated probe sets, there are 188 known or predicted genes with at least 1.5-fold change in both treatment regimens (supplemental Tables 3 and 8). Interestingly, all the commonly regulated genes share the same direction of gene expression value change compared with control in both experiments, with the exception of two genes as follows: sfrp2 (secreted frizzled related protein 2) and mepe (matrix extracellular phosphoglycoprotein). The former is proposed to be a secreted Wnt antagonist that directly interacts with Wnt ligands to inhibit signaling (23-25). MEPE was first identified in tumor conditioned medium and osteosarcoma cell lines. Multiple experimental observations support that MEPE is correlated positively to serum Pi and has an important role in negative regulation of mineralization in the skeleton (26, 27).
|
|
We conclude that CP134 is more potent in inducing gene expression changes in the rat compared with CP131 and CP334 as indicated by the number of genes regulated. A notable significant difference between CP134 and CP131 is their effect on the ctnnb1 gene, which encodes
-catenin, a WNT pathway member that is known to stimulate cell proliferation and repress apoptosis. It is decreased (5-fold) by CP131 treatment and unchanged by the other peptides. This drop in transcription level may indicate CP131 represses cell proliferation.
|
|
Kinetic Expression of PTH-regulated Genes—We performed microarray with RNA isolated from animals treated for 14 days to obtain the whole profiles of gene expression in response to anabolic or catabolic effects of PTH treatment regimens. However, the concern may be raised that the gene profiles may contain PTH secondary response genes, which are regulated through some PTH-induced primary response genes or could be the consequence but not the principal contributors of the anabolic or catabolic effects of PTH. To test if the genes are early response genes to either continuous or intermittent PTH administration, we also collected RNA from animals treated with single PTH-(1-34) injection or 6-h infusion (Fig. 5, C and D).
Furthermore, we collected RNA from rats receiving different durations of both intermittent (1, 7, and 14 days) and continuous (6 h and 3, 6, and 14 days) administration of PTH-(1-34) and vehicle. For the intermittent regimen, RNA were also collected at three different time points (1, 2, and 24 h after each injection), which enables us to obtain the kinetic expression pattern of PTH-regulated genes. We found interesting kinetic expression patterns of different genes that we examined by real time RT-PCR. At least two categories of gene expression pattern can be found in both PTH treatment regimens according to their kinetic expression changes. In the intermittent treatment, genes that encode signaling mediators (socs3), transcription factors (c-fos), cytokines (gro1), and proteases/protease inhibitors (spli and tfpi2) are regulated cyclically by each PTH injection; whereas genes associated with skeletal development, such as those indicative of extracellular matrix production and mineralization, show a different pattern of expression. For instance the gene encoding osteocalcin is not up-regulated until after several PTH injections, but then its mRNA level is sustained to the next PTH challenge. Because the stimulation is compared with the vehicle control at each time point instead of the basal for each injection, we are unable to tell if a gene such as osteocalcin has returned to basal at the time of PTH injection. However, we compared all the control vehicle-injected samples at the three time points, 1, 2, and 24 h, and no significant differences in the gene expression levels were found, at least not in the genes we have checked.
In the continuous regimen, similarly extracellular matrix genes (postn and alp) are not stimulated in the 6-h- and 3-day-treated samples but are at 6- and 14-day infusion. The up-regulation of these extracellular matrix genes in the infusion regimen is consistent with an increase in formation activity as reflected by the substantial increases in bone formation rate, as well as bone resorption, and increased bone turnover found in the histomorphometric analyses. The mRNA expression levels of other genes that are involved in different signaling pathways were kinetically regulated by continuous PTH administration. The anti-apoptotic gene bcl2 was up-regulated at the beginning of the infusion and then decreased significantly by 4-fold at 6 days and 6-fold at 14 days. pthr1 was substantially down-regulated after a 3-day infusion but bounced back and was up-regulated 2-fold at 6 days and 14 days of infusion. These results suggest the outcome of the catabolic effects of long term continuous treatment with PTH-(1-34) initiates with a wave of regulation of expression of these genes.
RANKL (tnfsf11) is a well known PTH-regulated gene, which plays a critical role during osteoclastogenesis. It was previously thought to be stimulated by continuous PTH to promote osteoclastic bone resorption. However, to our surprise, we observe a 5-fold increase in the anabolic (intermittent) regime but no significant increase in the catabolic regime by microarray. Consistent with our finding, RANKL is not regulated by continuous PTH treatment in the study from Onyia et al. (12). To confirm this result, we performed real time RT-PCR at different time points for RANKL. A strong and cyclic stimulation of RANKL mRNA is associated with intermittent injections of PTH-(1-34), which increases with increasing injections (Fig. 6A); a moderate but sustained stimulation is associated with continuous infusion of the hormone (Fig. 6D). The level of osteoprotegrin (OPG) and the ratio of RANKL/OPG were also examined because the latter is thought to be more relevant to bone resorption (Fig. 6, B, C, E, and F). With the intermittent PTH-(1-34) regime, the increases in OPG (Fig. 6B) lagged considerably behind the increases in RANKL resulting in a transient stimulation of the RANKL/OPG ratio (Fig. 6C). This may provide the molecular basis for a cyclic sharp up-regulation of bone resorption in the anabolic regime compared with a moderate but constantly elevated level of RANKL/OPG in the catabolic regime (Fig. 6F). The important point for anabolism may be the transient stimulation of osteoclast activity followed by the increase in OPG suppressing osteoclast activation. In contrast, in the continuous treatment, OPG mRNA levels decreased as time went on maintaining an elevated RANKL/OPG ratio.
Ingenuity Pathway Analysis—Our Ingenuity pathway analysis (IPA) (see Supplemental Material) indicates IP131 and IP134 induce similar if not identical effects. This conclusion is in agreement with our physiological results showing that they induce similar biological outcomes in the rats. In the network of the signaling pathways regulated by PTH we generated from IPA analysis (see Supplemental Material), potential interactions of three genes are proposed, all of which showed very high levels of gene expression change and are uniquely regulated by intermittent administration of PTH-(1-34) as follows: ccl2 (encodes monocyte chemoattractant protein-1), 142-fold; areg (amphiregulin), 20-fold; and crem (cAMP-responsive element modulator), 20-fold. Based on the IPA network graph, the overexpression of ccl2 and areg appears to be the most outstanding PTH-regulated genes unique to IP134 and IP131, i.e. the anabolic effect of PTH. The network analysis also revealed that differential regulation of Wnt signaling exists in the different PTH treatment regimens. In our accompanying article (42), we proposed a model to explain the importance of the transient stimulation of CCL2 for the dual effects of PTH. We have confirmed the stimulation of these genes in vivo by real time RT-PCR as follows: ccl2, areg, crem, wnt4, and sfrp4 are dramatically changed with intermittent PTH-(1-34) and PTH-(1-31) treatment in vivo (Fig. 7).
|
|
| DISCUSSION |
|---|
|
|
|---|
A recent report on the in vivo study of regulated gene expression profile of PTH studied the anabolic versus catabolic effects of PTH (12). This investigation detailed for the first time the broad profiling of the gene and pathway changes following treatment of intermittent versus continuous PTH-(1-34) in intact female adult rats. They utilized microarray analysis to examine gene expression either 24 h after the last injection or at cessation of the continuous infusion. These results extended previous observations of gene expression changes and revealed the in vivo regulation of bmp3 and multiple neuronal genes by PTH treatment. More than 50% (96 of 174 genes) of the continuous PTH-regulated genes reported in their study can be found in our results. Considering that we used different rat chips in the microarray experiments, and we treated animals for 14 days instead of 6 days, the overlap between the two studies of effects of continuous PTH administration in regulation of gene expression indicates very high reproducibility from one study to the other. Because we also performed microarray with 6-day PTH-(1-34) infused rats (data not shown), we were able to compare our 6-day continuous PTH-(1-34)-regulated gene list with the results from Onyia et al. (12). The overlap is the same as we see with the 14-day results indicating the differences between our results and theirs are largely due to the differences in chips and the duration of treatments.
However, with the intermittent PTH regimen, Onyia et al. (12) performed the microarray with RNA isolated 24 h after the last injection. As can be seen in our experiment and many other reports on PTH-regulated genes, it is notable that most of the regulated genes respond to PTH injection very promptly (16, 20). Hence, effects in many cases occur within 1 h after PTH injection and decline to basal levels within 24 h until the next PTH challenge due, we think, to feedback effects in PKA signaling. By examining the gene expression profile 24 h after the last PTH treatment, Onyia et al. (12) risked missing optimal time points to detect the responsive gene profile in the intermittent PTH administration regimen. Compared with this study, they found a very limited number of genes (41 genes) were moderately regulated (the highest fold change is 3.4-fold increase of Cdh2, which is up-regulated 2.6-fold at 1 h in our study) by intermittent PTH treatment. At least 24 of these 41 genes have also been found regulated at 1 h after the last injection according to our results, most of them are matrix-associated proteins, adhesion proteins, or chaperones. They only found three signal transduction genes regulated, protein-tyrosine phosphatase receptor (1.7-fold), olfactory receptor 50 (1.6-fold), and angiotensinogen (-1.6-fold). The first gene was also detected in our data with a 1.6-fold increase indicating this receptor is constantly up-regulated by intermittent administration of PTH-(1-34). In conclusion, our study is the first to examine the optimal time point after PTH administration and thus achieves the most informative and thorough gene expression profile regulated by PTH in vivo. Moreover, we not only performed PTH treatment intermittently and continuously, which enabled us to compare the gene expression changes associated with the paradoxical effects of PTH treatments, but also applied three different PTH peptides simultaneously and compared the bone responses and gene expression profiles among the different PTH peptides to be able to determine the role of each PTH-stimulated signaling pathway (PKA versus PKC) in anabolic and catabolic effects of PTH.
Although most of the known PTH-regulated genes are PKA pathway-dependent according to in vitro experiments, the role of PKA and PKC pathways and their relation/interaction during the dual effects of PTH in vivo has never been addressed. To elucidate the signaling pathways for the anabolic effects of PTH on bone, in this study, we used amino-terminally and carboxyl-terminally truncated PTH analogs, PTH-(3-34), which binds to the PTH1R and activates the PKC/calcium pathways but does not elicit a cAMP response, and PTH-(1-31), which activates the cAMP/PKA pathway and calcium signaling but has no effect on PKC activation, based on previous findings (8, 34). The comparison of the different PTH peptide effects in both intermittent and continuous administration regimes in vivo provides important information to understand the relationship between the PTH-activated signaling pathway and the dual effects of the hormone on bone metabolism. Despite the fact that PTH-(3-34) may have a lower binding affinity to the PTH1R than PTH-(1-34) and PTH-(1-31), it still displays substantial binding affinity (34, 35), and it has been shown that PTH-(3-34) can specifically activate calcium/PKC systems in vivo and in vitro at the same concentration as PTH-(1-34) (36). From our data, it seems that the PKC system, i.e. PTH-(3-34), is not able to contribute to the catabolic effects of PTH because no alterations in genes and bone resorption markers are observed in response to this peptide with continuous infusion. Nevertheless, PTH-(1-31) is not as effective as PTH-(1-34) in regulation of gene expression with the continuous treatment regimen. This indicates that the PKC system is likely to be involved in the enhancement of the catabolic effects by interacting with the PKA system. Interestingly, activation of the PKA system is sufficient for the anabolic effects because similar bone phenotypes and very similar changes in gene expression profiling are seen with intermittent administration of either PTH-(1-34) or PTH-(1-31). Intermittent administration of PTH-(3-34) regulated many fewer genes compared with the other two PTH peptides and did not have significant anabolic effects on bone. Thus, we conclude that activation of the PKC system in the intermittent treatment regimen of PTH-(1-34) has a minor role in its anabolic effects. More importantly, PTH-(1-31) demonstrates similar (if not identical) anabolic effects to PTH-(1-34) but is less effective in inducing catabolic effects. This may suggest that PTH-(1-31) could serve as a better anabolic agent than PTH-(1-34) therapeutically.
Besides the complexity of different PTH peptides initiating different signaling pathways and contributing differentially to anabolic and catabolic effects of the hormone, the kinetic profile of PTH-responsive genes is another important factor causing the paradoxical physiological responses to intermittent and continuous PTH administration. We have shown that some extracellular matrix genes changed only after treatment for 1 week, whereas some cytokines, growth factors, and transcription factors respond to PTH administration very rapidly and dramatically but are not sustained when the hormone levels decline. Continuous exposure to the hormone may cause some of the same gene expression changes, but the changes tend to be moderate (usually about 3-fold) and sustained. The different regulation patterns of the same gene provide a clue to understanding the molecular basis of the anabolic and catabolic effects of PTH. It is known that both intermittent and continuous treatment of PTH increase bone remodeling, whereas the anabolic effect is associated with the intermittent protocol, and the catabolic effect is associated with the continuous protocol.
In this study, we demonstrate that a key regulator of osteoclastogenesis and activity, RANKL, is increased with both treatment regimes. It is interesting that the increase in RANKL in the intermittent treatment is very dramatic but transient but increases with increasing injections, whereas in the continuous treatment, its level is moderately up-regulated but sustained throughout the treatment. We think this is because the coupling of osteoclastic bone resorption and osteoblastic bone formation requires transient activated resorption ahead of osteoblasts lining the remodeling sites and building new bone. One of the reasons the anabolic effect is achieved may be because transient up-regulation of RANKL causes short term increased bone resorption that is succeeded by increased bone formation; on the contrary, the catabolic effect produces sustained elevated osteoclastic bone resorption from the sustained up-regulation of RANKL from continuous infusion of PTH. It should be noted that there is a potential contribution of osteocytes to these gene expression data.
Overall, the microarray expression data support a provocative and paradoxical theme, i.e. for PTH therapy to be fully anabolic in nature, there needs to be a transient and cyclic up-regulation of genes classically associated with resorption such as RANKL and interleukin-6. This theme has elements similar to the process of an acute inflammatory response supported by the microarray and kinetic data (apart from RANKL and interleukin-6, CXCL1 and CCL2/MCP-1 are highly increased, see accompanying article Ref. 42). This is in distinct contrast to the catabolic regime where expression of these genes is much lower but, importantly, is sustained.
In summary, we have conducted a comprehensive analysis of the genes that are differentially regulated in vivo by intermittent and continuous treatment with different PTH peptides using a gene array approach. These data support earlier histomorphometric observations in our hands and other publications (37-39). The profiling of gene expression highlights important similarities and differences between intermittent and continuous treatment regimens of PTH. PTH-(1-34) has strong physiological effects and differentially regulates a large number of genes in both treatment regimens. PTH-(1-31) reproduces the effects of PTH-(1-34) in the anabolic regimen but has a lesser catabolic effect in the continuous regimen and thus regulated a lesser number of genes. PTH-(3-34) regulates a very limited number of genes in the intermittent regimen and none in the continuous regimen. We conclude that activation of the PKC system in the intermittent treatment regimen of PTH-(1-34) has a minor role in its anabolic effects but may interact with PKA signaling to enhance the catabolic effect of PTH-(1-34) even though sole activation of the PKC system fails to elicit the catabolic effects.
A large number of genes such as slpi, tfpi2, socs3, and gro1, that were previously not considered to be expressed in bone or to be regulated by either PTH treatment regimen were identified and validated. The kinetic expression patterns of some of the PTH-regulated genes in both treatment regimens were also described. Recently, follicle-stimulating hormone (FSH) was reported to directly regulate bone resorption independent of estrogen (40). Our microarray data found the FSH-responsive gene (fshprh1), which is stimulated by FSH (41), was down-regulated by PTH-(1-34) in both treatment regimes suggesting the interaction of FSH and PTH signaling in bone. The function of these newly identified genes in bone and their precise role, if any, in the anabolic and catabolic bone state deserve further investigation.
Finally, by using the knowledge base of biological networks provided by IPA, we set up a novel signaling network that is differentially regulated in the different PTH treatment regimes and suggest an association of this network with the dual effects of PTH in vivo. This information will be useful for further distinguishing the critical mediators of modeling and remodeling processes and for identification of key targets for therapeutic intervention in disorders of bone metabolism.
| FOOTNOTES |
|---|
The on-line version of this article (available at http://www.jbc.org) contains supplemental Materials, Results and Discussion, Tables 1-8, Fig. 1, and additional references. ![]()
1 To whom correspondence should be addressed: Dept. of Physiology and Biophysics, UMDNJ-Robert Wood Johnson Medical School, 675 Hoes Ln., Research Tower 561, Piscataway, NJ 08854. Tel.: 732-235-4552; Fax: 732-235-3977; E-mail: partrinc{at}umdnj.edu.
2 The abbreviations used are: PTH, parathyroid hormone; h, human; RT, reverse transcription; ELISA, enzyme-linked immunosorbent assay; RANKL, receptor activator of NF-
B ligand; PKA, cAMP-dependent protein kinase; PKC, protein kinase C; PMA, phorbol 12-myristate 13-acetate; OPG, osteoprotegrin; IPA, Ingenuity pathway analysis; FSH, follicle-stimulating hormone. ![]()
| ACKNOWLEDGMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
D. Onan, E. H. Allan, J. M. W. Quinn, J. H. Gooi, S. Pompolo, N. A. Sims, M. T. Gillespie, and T. J. Martin The Chemokine Cxcl1 Is a Novel Target Gene of Parathyroid Hormone (PTH)/PTH-Related Protein in Committed Osteoblasts Endocrinology, May 1, 2009; 150(5): 2244 - 2253. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. M. Kent, L. J. C. Smyth, J. Plumb, C. L. Clayton, S. M. Fox, D. W. Ray, S. N. Farrow, and D. Singh Inhibition of Lipopolysaccharide-Stimulated Chronic Obstructive Pulmonary Disease Macrophage Inflammatory Gene Expression by Dexamethasone and the p38 Mitogen-Activated Protein Kinase Inhibitor N-cyano-N'-(2-{[8-(2,6-difluorophenyl)-4-(4-fluoro-2-methylphenyl)-7-oxo-7,8-dihydropyrido[2,3-d] pyrimidin-2-yl]amino}ethyl)guanidine (SB706504) J. Pharmacol. Exp. Ther., February 1, 2009; 328(2): 458 - 468. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Ogita, M. T. Rached, E. Dworakowski, John. P. Bilezikian, and S. Kousteni Differentiation and Proliferation of Periosteal Osteoblast Progenitors Are Differentially Regulated by Estrogens and Intermittent Parathyroid Hormone Administration Endocrinology, November 1, 2008; 149(11): 5713 - 5723. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. B. Lowry, S. Lotinun, A. A. Leontovich, M. Zhang, A. Maran, K. L. Shogren, B. K. Palama, K. Marley, U. T. Iwaniec, and R. T. Turner Osteitis Fibrosa Is Mediated by Platelet-Derived Growth Factor-A Via a Phosphoinositide 3-Kinase-Dependent Signaling Pathway in a Rat Model for Chronic Hyperparathyroidism Endocrinology, November 1, 2008; 149(11): 5735 - 5746. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. M. Teplyuk, M. Galindo, V. I. Teplyuk, J. Pratap, D. W. Young, D. Lapointe, A. Javed, J. L. Stein, J. B. Lian, G. S. Stein, et al. Runx2 Regulates G Protein-coupled Signaling Pathways to Control Growth of Osteoblast Progenitors J. Biol. Chem., October 10, 2008; 283(41): 27585 - 27597. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. Alam, Q. Sun, L. Liu, D. L. Koller, Y. Liu, H. J. Edenberg, M. J. Econs, T. Foroud, and C. H. Turner Genomic expression analysis of rat chromosome 4 for skeletal traits at femoral neck Physiol Genomics, October 8, 2008; 35(2): 191 - 196. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Li, L. Qin, M. Bergenstock, L. M. Bevelock, D. V. Novack, and N. C. Partridge Parathyroid Hormone Stimulates Osteoblastic Expression of MCP-1 to Recruit and Increase the Fusion of Pre/Osteoclasts J. Biol. Chem., November 9, 2007; 282(45): 33098 - 33106. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| All ASBMB Journals | Molecular and Cellular Proteomics |
| Journal of Lipid Research | ASBMB Today |