- Research article
- Open Access
Transcriptional profiling of PRKG2-null growth plate identifies putative down-stream targets of PRKG2
BMC Research Notesvolume 8, Article number: 177 (2015)
Kinase activity of cGMP-dependent, type II, protein kinase (PRKG2) is required for the proliferative to hypertrophic transition of growth plate chondrocytes during endochondral ossification. Loss of PRKG2 function in rodent and bovine models results in dwarfism. The objective of this study was to identify pathways regulated or impacted by PRKG2 loss of function that may be responsible for disproportionate dwarfism at the molecular level.
Microarray technology was used to compare growth plate cartilage gene expression in dwarf versus unaffected Angus cattle to identify putative downstream targets of PRGK2.
Pathway enrichment of 1284 transcripts (nominal p < 0.05) was used to identify candidate pathways consistent with the molecular phenotype of disproportionate dwarfism. Analysis with the DAVID pathway suite identified differentially expressed genes that clustered in the MHC, cytochrome B, WNT, and Muc1 pathways. A second analysis with pathway studio software identified differentially expressed genes in a host of pathways (e.g. CREB1, P21, CTNNB1, EGFR, EP300, JUN, P53, RHOA, and SRC). As a proof of concept, we validated the differential expression of five genes regulated by P53, including CEBPA, BRCA1, BUB1, CD58, and VDR by real-time PCR (p < 0.05).
Known and novel targets of PRKG2 were identified as enriched pathways in this study. This study indicates that loss of PRKG2 function results in differential expression of P53 regulated genes as well as additional pathways consistent with increased proliferation and apoptosis in the growth plate due to achondroplastic dwarfism.
PRKG2 is a key regulatory kinase in the proper temporal and spatial development of growth plate cartilage. Loss of PRKG2 function results in a disorganized mixture of hypertrophic and proliferative chondrocytes instead of the highly organized columnar array observed in properly formed growth plates . Mouse knockouts and natural mutations in PRKG2 in rats and cattle exhibit achondroplastic dwarfism [1-3]. A deletion of PRKG2 on human chromosome 4q21 was recently associated with growth restriction and mental retardation . Thus, the functional role of PRKG2 in growth plate development is highly conserved across species. Based on the impact of PRKG2 on growth plate elongation across species, it is likely that many of the pathways regulated downstream of PRKG2 are also highly conserved.
Several targets of PRKG2 have been identified previously. PRKG2 signaling regulates growth plate chondrocyte hypertrophy and hyperplasia via SRY, sex determining region Y, −box 9 (SOX9), glycogen synthase 3 beta (GSK3b), and additional unknown factors [2,5,6]. PRKG2 phosphorylation is required for nuclear translocation of SOX9. Phosphorylated SOX9 modulates collagen expression from a proliferative (collagen 2, COL2) to a hypertrophic (collagen 10, COL10) program . It is currently unclear whether PRKG2 directly or indirectly regulates SOX9 signaling . Phosphorylation of GSK3b at Serine 9 by a protein complex comprised of PRKG2 and Axin1 causes increased beta-catenin (CTNNB1) activity and increased hypertrophy of growth plate chondrocytes . Haploinsufficiency of GSK3b was capable of 30-40% rescue of skeletal growth in PRKG2 −/− mice .
Several regulators of PRKG2 have also been identified. Notable regulators of PRKG2 include C-type naturetic peptide (CNP), cyclic guanosine monophosphate (cGMP) and nitric oxide (NO). CNP regulates multiple pathways required for proper endochondral ossification. Various mouse crosses have demonstrated that transgenic modulation of CNP can rescue mice from achondroplasia caused by disruption of fibroblast growth factor receptor 3 (FGFR3), mitogen-activated protein kinase (MAPK), or PRKG2 signaling [7-9]. However, no direct studies have been performed to determine if FGFR3 and MAPK are regulated by PRKG2. Additionally, cGMP and NO act as upstream triggers for PRKG2 signaling. These small molecules have been tied to calcium signaling, apoptosis, and proliferation processes [10,11] as well as regulation of MAPK, nuclear factor of kappa light polypeptide gene enhancer in B-cells (NF-κB), mitogen-activated protein kinase 8 (JNK/ MAPK8), vascular endothelial growth factor (VEGF), and cAMP responsive element binding protein 1 (CREB) signaling .
To date, no studies have dissected the role of PRKG2 regulation of global gene expression in the growth plate. Our objective was to characterize transcriptional changes in pathways downstream of PRKG2 in hopes of discovering additional regulators responsible for the switch from proliferative to hypertrophic growth plate development. A secondary objective is to understand the changes in transcriptional networks in the growth plate due to achondroplastic dwarfism. Transcriptional profiling of Angus PRKG2 R678X/R678X (dwarf) versus PRKG2 R678X/+ (unaffected) cattle was used as a means to identify downstream targets of PRKG2.
Analysis of differential expression
No differentially expressed (DE) genes were detected across genotypes after correcting for multiple testing. Since the statistical power to detect differences in gene expression was limited, pathway analyses were used as a filter to prioritize DE genes. The rationale was that transcripts within a DE signaling pathway or biological process were likely to be over-represented even when all genes were not DE. A nominal p < 0.05 significance level was used for each photomechanical transfer unit (PMT) level to declare “significance” as a first filter for the data. We focused on pathway analysis of the PMT 70 DE gene list, because it appeared the most valid based on the p-value histogram. The PMT 80 and PMT 90 gene lists were used in combined analyses only to corroborate results from the PMT 70 DE gene list.
Summary of nominally DE genes
Using a significance level of p < 0.05, 1284 genes were detected as DE at scanning PMT 70, 860 genes at PMT 80, and 628 genes at PMT 90. Overlap between all three PMT levels included 180 genes. A summary of the significant probes is provided in Additional file 1.
DAVID enrichment results are consistent with known PRKG2 biology
The gene ontology (GO) classification of genes represented at PMT 70 showed the greatest over-representation of genes involved in apoptosis, cell cycle regulation, cell death, protein modification, and phosphorylation (p < 0.01, Family-wise false discovery rate (FWFDR) < 0.10) using the database for annotation, visualization and integrated discovery (DAVID) tool (Table 1). Terms over-represented from the Panther molecular function database included: non-receptor serine threonine protein kinase, nucleic acid binding, transferase, protein kinase, non-motor actin and guanyl-nucleotide exchange factor (p < 0.01, FWFDR < 0.10) (Table 2). Functionally enriched gene clusters with statistically significant enrichment score (ES) included: proliferation (ES = 2.9), actin cytoskeletal development (ES = 2.25), and differentiation and apoptosis (ES = 2.17). The Angiotensin II mediated activation of JNK Pathway via protein tyrosine kinase 2 beta (PYK2) dependent signaling pathway was over-represented in the Biocarta pathway database (p < 0.05). Over-represented KEGG pathways included: MAPK signaling, Adherens junction, and Glycan structures - biosynthesis 2 (p < 0.05).
Enriched biological processes from the combined analysis, including genes from the PMT70, PMT 80 and PMT 90 gene lists, were very similar to those found from the PMT 70 gene list (Table 3). Notable biological processes over-represented in this analysis include: apoptosis, cell cycle regulation, and cell death (p < 0.005, FWFDR < 0.10). Notably, wingless-type MMTV integration site family (WNT) signaling was over-represented in the combined analysis (p < 0.05) (Table 4). Functional gene clustering results were also similar to those observed in the PMT 70 data, including clusters enriched for proliferation (ES = 3.69) and differentiation and apoptosis (ES = 3.69) (Additional file 2). Note, an ES ≥ 1.3 is equivalent to p < 0.05. .
Pathway studio identifies hub genes important to apoptosis and cell cycle
We used Pathway Studio software to search for genes that group into pathways based on published results in PubMed. Many pathways were differentially represented between the two PRKG2 genotypes. For this reason, we focused our attention on pathways relevant to the known biological function of PRKG2. A sample of highly relevant hub genes indentified by Pathway studio included: cAMP responsive element binding protein 1 (CREB1), cyclin-dependent kinase inhibitor 1A (CDKN1A/ P21/ Cip1), catenin (cadherin-associated protein), beta 1, 88 kDa (CTNNB1), epidermal growth factor receptor (EGFR), E1A binding protein p300 (EP300), jun proto-oncogene (JUN), tumor protein P53 (P53), ras homolog family member A (RHOA), and SRC proto-oncogene, non-receptor tyrosine kinase (SRC). Each of these genes had many differentially regulated genes (p < 0.05) within their pathway across PRKG2 genotypes using the common regulators analysis (Additional file 3). The direct regulators analysis also suggested that P21, CTNNB1, EGFR, and JUN pathways were differentially regulated across PRKG2 genotypes. Additional pathways/hub genes with differentially expressed target genes across genotypes are listed in Additional file 3. A combined analysis of all differentially expressed genes found across PMT levels identified additional target genes in the P53 pathway (n = 476 genes downstream of P53). Figure 1 shows the P53 target genes detected from the PMT 70 gene list by Pathway Studio.
Real-time PCR results confirmed differential expression of P53 target genes
We tested eight genes to validate their status as differentially expressed, including seven genes from the P53 pathway. The primer sequences used in real-time qPCR for these genes are provided in Additional file 4. Five of these seven P53 target genes were confirmed as DE (p < 0.05). Three of these genes (breast cancer 1, early onset (BRCA1), CCAAT/ enhancer binding protein (C/EBP), alpha (CEBPA), and CD58 molecule (CD58)) were detected as DE using the ddCT method while two genes, BUB1 mitotic checkpoint serine/threonine kinase (BUB1) and vitamin D (1,25- dihydroxyvitamin D3) receptor (VDR), were detected as DE by difference in natural log starting copy number (p < 0.05). Figure 2 presents the coordinated change in expression of these five genes.
DAVID results indicated that specific transcripts involved in apoptosis, cellular hypertrophy, DNA damage repair, and proliferation are differentially enriched as a result of bovine achondroplasia
The enrichment results from DAVID identified global changes in gene expression related to cell development, proliferation, and death. The most significantly enriched GO biological process for the PMT 70 data was cell development, likely representing genes related to hypertrophy and terminal differentiation. Furthermore, cell proliferation and cell cycle genes were also identified; likely due to the increased proliferative phenotype of PRKG2 null (PRKG2 R678X/R678X) chondrocytes. Last, and not unexpected, were changes related to cell death and response to stress. These changes in gene expression are likely related to apoptosis in the growth plate due to uncoordinated development of the growth plate, which results in the inappropriate mixing of hypertrophic and proliferative growth zones. A combined analysis of differentially expressed gene lists across all of the PMT levels resulted in similar, more statistically significant, enrichment of transcripts for biological processes related to PRKG2 function. Enriched biological processes from combined PMT data include: apoptosis, cell hypertrophy/differentiation, cell proliferation and regulation of cell cycle, and processes related to proper organization of the growth plate (cell motility, locomotory behavior and established cell localization). WNT and calcium signaling pathways were over-represented in a KEGG pathway enrichment analysis using the combined gene lists. Additionally, changes were observed in glycan structure processes, likely due to the disorganization of the growth plate in PRKG2 null individuals. A key result from the enrichment analysis across all PMT was the detection of over-represented WNT pathway transcripts, which was recently shown to regulate PRKG2 .
The functional gene-clustering tool in DAVID identified modules of genes enriched for actin cytoskeletal and cell motility, apoptosis, kinase/kinase function, proliferation/ cell cycle and metal ion binding from the PMT 70 gene list (ES > 2.9). Combining all of the differentially expressed gene lists across PMT levels resulted in similar results with the addition of DNA damage repair (ES = 1.66) and embryonic appendage/limb development gene clusters (ES = 1.62). The biological functions enriched in these gene clusters are consistent with increased proliferative and altered morphology observed in PRKG2 null (dwarf) vs. wild-type (unaffected) individuals. In summary, DAVID enrichment results are consistent with known biological functions of PRKG2 and phenotypes observed in null animals.
Pathway studio identified two known PRKG2 targets in CTNNB1 and SRC as well as a host of other possible down-stream targets involved in apoptosis and proliferation
Many over-represented pathways detected by Pathway studio appear to be involved primarily in the control of proliferation, apoptosis or a combination of both processes [14-22]. Gene pathways differentially expressed between genotypes include: P21, CREB1, EP300 (p300), EGFR, JUN, SRC, RHOA, and CTNNB1. P21 appears to be largely inhibited or down-regulated in PRKG2 wild-type individuals compared to PRKG2 null individuals. JUN, CREB1, EP300, RHOA, and SRC target genes appear to be up-regulated in PRKG2-null individuals. Beta-catenin (CTNNB1) was previously identified as a downstream target of PRKG2 and an important regulator of bone growth in the WNT pathway . Additionally, it has been suggested that other regulators identified in this study (i.e. EP300 and EGFR) interact with CTNNB1 to regulate cellular differentiation and proliferation respectively [15,21]. SRC is regulated by PRKG2 mediated protein tyrosine phosphatase, non-receptor type 6 (SHP-1) phosphorylation in mechanically stimulated osteoblasts . In turn, SRC acts as a mitogen by promoting cell growth through RAS and the Raf-MEK-ERK pathways. This study also finds a RAS homologue pathway, RHOA, to be differentially represented between the PRKG2 null and unaffected individuals.
P53 is known to regulate proliferation and apoptosis in the growth plate
The P53 pathway is a logical candidate for PRKG2 regulation because P53 is a key regulator of proliferation, apoptosis and response to DNA damage by regulating target gene expression [24-28]. P53 regulates critical checkpoints in the cell cycle. Increased P53 is known to increase P21 expression, which is essential for progression from G1 to S phase of cell cycle . Notably, P21 is upregulated in a variety of achondroplasia phenotypes in humans [29,30]. The location and time where apoptosis occurs in the growth plate is critical for proper long bone growth. To allow proper proliferation of chondrocytes, P53 is down regulated in the resting zone to prevent apoptosis . At the same time, apoptosis is required for deposition of new bone at the epiphyseal-metaphyseal junction to allow elongation of long bones [32-34]. P53 induced apoptosis may be important to facilitate this process in long bone growth. Consistent with this hypothesis, data indicates that P53 −/− mouse chondrocyte progenitor cells have impaired terminal differentiation and enhanced proliferation as well as progression into osteosarcoma [24,35]. Interestingly, PRKG2 is located underneath a quantitative trait locus on mouse chromosome 5 for alpha-radiation induced osteosarcoma . The connection between PRKG2 and P53 may be GSK3. GSK3 has been shown to regulate P53 via MDM2 proto-oncogene (MDM2) phosphorylation, leading to down regulation of P53 .
To validate the changes in P53 target gene’s expression, we confirmed five of seven P53 target genes as differentially expressed by real-time PCR (p < 0.05). Gene expression profiles indicated down-regulation of P53 in PRKG2 null individuals, consistent with increased cellular proliferation. The P53 regulated transcripts BRCA1, and VDR were confirmed as up-regulated in PRKG2 null individuals, while BUB1, CD58, CEBPA were down-regulated. These results indicate that PRKG2 may act to turn on or turn off P53 either directly or indirectly. Our microarray experiment indicated that MDM2 was significantly (p < 0.001) up regulated in PRKG2 null individuals, but we were unable to find more than a trend towards up-regulation of MDM2 in PRKG2 nulls (p < 0.20). MDM2 acts as part of a feedback loop to regulate P53, first promoting and then inhibiting P53 expression. No protein or functional data was generated in this experiment to confirm these results.
Since MDM2 phosphorylation down regulates P53 signaling, we asked if MDM2 could be phosphorylated by PRKG2, acting as a direct target of PRKG2 regulation that mediates P53 repression. Using the GPS software , we were able to predict several PRKG2 phosphorylation sites on both MDM2 and P53. The maximum log of odds (LOD) score for MDM2 phosphorylation by PRKG2 was 7.0, where LOD = 1.69 is equivalent to p < 0.05. The most significant LOD score for P53 phosphorylation by PRKG2 was 4.67. In comparison, SOX9, a known PRKG2 substrate, had a PRKG2 phosphorylation site predicted with a LOD = 9.33.
We identified differentially represented pathways in growth plate cartilage due to dwarfism caused by loss of PRKG2 function. Biological processes differentially regulated due to PRKG2 loss of function based on ontology enrichment analysis include: apoptosis, cellular differentiation, organization and proliferation. Hub genes with differentially expressed targets included: P21, CREB1, EP300 (P300), EGFR, JUN, SRC, RHOA, CTNNB1, and P53. Both SRC and CTNNB1 were previously identified as downstream targets of PRKG2. We validated five of eight DE genes from a microarray experiment by real-time PCR, including five P53 regulated genes, demonstrating that pathway enrichment can identify truly DE genes when no other DE genes are detected. Differential expression of P53 target genes between PRKG2 null and wild-type individuals may indicate another way in which PRKG2 regulates growth plate chondrocyte proliferation and apoptosis. Alternatively, P53 expression across PRKG2 genotypes may vary as a result of cellular response to stress caused by DNA damage and cellular death in cells that are unable to terminally differentiate in a reactive fashion to loss of PRKG2 function. Regardless of the cause of P53 differential expression, the change in P53 target gene expression is important given the role of P53 in growth plate development and regulation of apoptosis. These results demonstrate that important differences in gene expression in response to or as a result of PRKG2 loss of function were identified by pathway enrichment analysis. Identification of differentially enriched pathways due to loss of PRKG2 function may be useful in suggesting potential treatments for achondroplasia and other bone diseases.
Dwarf (PRKG2 R678X/R678X) and unaffected (PRKG2 R678X/+) calves were produced by mating a mother, daughter (carrier, dwarf) pair to a single carrier bull. Standard embryo transfer methods were used to allow implantation into surrogate mothers. More than twenty embryos were produced; however, only six survived to 60 days post-implantation. Six calves, three from each female, were born full term from surrogate mothers in two calving groups. Calves within the two groups were born during the same week and were reared and euthanized together. Two dwarves and one unaffected individual were born in each calving group. Each animal was genotyped to verify the dwarf or unaffected phenotype . All calf rearing, experimental procedures, tissue collection and euthanization were carried out under the guidelines of the Iowa State University Institutional Animal Care and Use Committee under approved protocol number 6-04-5697-B-1.
Calves were euthanized at approximately 210 days of age during a period of rapid long bone growth . Whole growth plate was collected from the tibia bone of all the animals and preserved using RNAlater™ (Ambion, Austin, TX). Each tibia was sliced horizontally at the growth plate to allow cartilage pieces to be removed using a sterile razor blade. Extra care was taken to avoid harvest of bone marrow and osteoblasts by leaving a thin layer of growth plate on the bone epiphysis and diaphysis. Samples were archived at - 80°C until processed for RNA isolation.
RNA was isolated from growth plate cartilage using TRIzol™ (Invitrogen, Carlsbad, CA). A mass of 0.15 to 0.20 g of growth plate cartilage was placed into a sterile Petri dish with one milliliter (ml) of TRIzol™. Next, each sample was diced into a slurry using a sterile scalpel. Then, samples were homogenized using a rotor star probe in a total volume of 6 ml TRIzol™. Sample yield was measured by spectrophotometer. These initial RNA yield values were used to estimate how much purified RNA could be eluted by column purification (Qiagen, Valencia, CA). Column purified RNA was stored at - 80°C for future use.
An Oligonucleotide 70 mer microarray created by the bovine microarray consortium was used to assay the expression of approximately 24,000 transcripts . Slides were UV cross-linked prior to each hybridization. Synthesis of cDNA was carried out using Superscript II reverse transcriptase (Invitrogen, Carlsbad, CA) along with the Genisphere 3DNA 350 HS kit components (Genisphere, Hatfield, PA). Unaffected animals, PRKG2 R678X/+, were compared within calving group to each of the two dwarf, PRKG2 R678X/R678X , calves, resulting in two microarrays per calving group. Cy 3 and Cy 5 dyes were used for differential display of cDNA binding between genotype treatments to the array. Dye swapping was used for genotype comparisons to account for differential dye incorporation. A total of four microarrays were hybridized to compare PRKG2 R678X/R678X and PRKG2 R678X/+ animals using the Genisphere kit. Microarrays were hybridized using GE Healthcare's Lucidea SlidePro hybridization oven (General Electric Healthcare, Piscataway, NJ) at the Iowa State University DNA facility according to manufacturer’s protocol. All microarray data was submitted to NCBI GEO under accession number GSE31627.
Hybridized slides were scanned for visualization of Cy 3 and Cy 5 signal using a Pro Scan Array HT scanner at the microarray facility within the Iowa State University, Plant Sciences Institute. Slides were scanned at three different laser PMT intensities, 70, 80 and 90 PMT, to account for high and low expression levels of various transcripts  as different genes will appear as expressed or not expressed at the different PMT levels. Images were analyzed to obtain raw signal intensities and localized background signals with Imagene 5.1 software. Text files were produced for each slide channel for statistical analysis.
Image files were normalized by channel to correct for background, differential dye incorporation using Lowess normalization , and heterogeneous variance by median centering and natural log transformation using R scripts kindly provided by Dr. Dan Nettleton, Iowa State University. Data were analyzed using PROC MIXED of SAS accounting for fixed effects of the dye, genotype and calving group. The slide was included as a random effect. To account for random effects having an estimated variance of zero, we used the Kenworth-Rodgers correction. Multiple testing corrections were made by calculating q-values to determine false discovery rate (FDR) levels for each gene list .
Pathway and gene function analysis
Genes on the bovine microarray consortium 70-mer array were annotated where annotation or BLAST analysis could not identify the gene corresponding to each probe. We used the SWISS-PROT protein IDs to perform analysis of gene function and search for over-represented gene pathways using DAVID [13,43] and Pathway Studio (Ariande Genomics, Rockville, MD). Gene expression analysis with DAVID identified over-represented gene categories, genes of specific molecular function, functional gene clusters and pathways. Pathways identified by DAVID were mined from the KEGG  and Biocarta pathway analysis tools. The default settings for EASE were used for all of these analyses. Pathway studio was used to find specific gene pathways that were highly represented within a gene list. Pathway studio uses a wave propagation algorithm to text mine databases (i.e. curated literature data from PubMed) for information about groups of entities used in analysis and tries to relate these groups back to the gene list. We used the direct regulators (relationship between genes on a gene list) and common regulators (identifies upstream regulators that have targets within a gene list) analyses. Gene pathways were selected as differentially represented based on the number of connections originating from a common regulator that differed between the two genotypes.
The cDNA from growth plate cartilage was analyzed for differential expression by replicating each sample three times for the target and housekeeping gene in two separate plate replicates (i.e. six total replicates of both the target and housekeeping gene) using primers listed in Additional file 4 for the following seven P53 regulated genes: BRCA1, BUB1 CEBPA, CD58, MDM2, PRL and VDR. In addition, the SOX9 responsive gene, phosphate regulating endopeptidase homolog, X-linked (PHEX), was assayed across genotypes by real-time PCR. We used beta-actin as our housekeeping gene. Primer template sequences were retrieved using the 3’ human exonic sequence as a BLAST template to find the bovine template sequence. Genes analyzed by real-time PCR were selected based on the Pathway Studio analysis. Eight target genes and one house keeping gene (Beta-actin) were analyzed to determine the validity of microarray gene expression for the P53 pathway. Primers, PMT levels, and thermocycling protocols are listed in Additional file 4. Standard curves for each target and reference gene were made and evaluated for linearity in amplification. All real-time PCR data were analyzed using MyiQ™ software to check for primer dimmer, linear amplification (BioRad; Hercules, CA).
Statistical analysis of real-time PCR
Real-time PCR data was analyzed using the delta-delta cycle threshold (ddCt) method  in Proc GLM of SAS (SAS Institute; Cary, NC). The model used to analyze the expression data was
where y = ddCt (i.e. the target gene – housekeeping gene (beta-actin) for the dwarf subtracted from the target gene– housekeeping gene for an unaffected individual), p = real-time plate replicate, g = genotype, gr = common calving and rearing group, and e = error. When genes were not detected as differentially expressed using the ddCT method, we analyzed the difference in log starting copy number where the copy number of each target gene was compared against a standard curve for the same transcript. The model used to analyze the natural log (ln) starting copy number data was
where y = ln starting value of the target gene mRNA adjusted for ln starting value of housekeeping gene mRNA (ln starting number target – ln starting number beta-actin), p = real-time plate replicate, g = R678X PRKG2 genotype, gr = common calving and rearing group, and e = error.
Basic local alignment search tool
Breast cancer 1, early onset
BUB1 Mitotic checkpoint serine/threonine kinase
Cyclin-dependent kinase inhibitor 1A (p21, Cip1)
Complimentary deoxyribonucleic acid
CCAAT/ Enhancer binding protein (C/EBP), alpha
Cyclic guanosine monophosphate
C-type naturetic peptide
cAMP Responsive element binding protein 1
Catenin (cadherin-associated protein), beta 1, 88 kDa
Database for annotation, visualization and integrated discovery
Epidermal growth factor recptor
E1A Binding protein p300
Fibroblast growth factor receptor 3
Family-wise false discovery rate
Glycogen synthase 3 beta
Mitogen-activated protein kinase 8 (MAPK8)
Kyoto encyclopedia of genes and genomes
Mitogen-activated protein kinase
MDM2 proto-oncogene, E3 ubiquitin protein ligase
Major histocompatability complex
Mucin 1, cell surface associated
Tumor protein P53
Nuclear factor of kappa light polypeptide gene enhancer in B-cells
Polymerase chain reaction
Phosphate regulating endopeptidase homolog, X-linked
Photomechanical transfer unit
cGMP-dependent, type II, protein kinase
Protein tyrosine kinase 2 beta
Ras homolog family member A
Protein tyrosine phosphatase, non-receptor type 6 (PTPN6)
SRY, sex determining region Y, −box 9
SRC proto-oncogene, non-receptor tyrosine kinase
Vitamin D (1,25- dihydroxyvitamin D3) receptor
Vascular endothelial growth factor
Wingless-type MMTV integration site family
Pfeifer A, Aszodi A, Seidler U, Ruth P, Hofmann F, Fassler R. Intestinal secretory defects and dwarfism in mice lacking cGMP-dependent protein kinase II. Science. 1996;274(5295):2082–6.
Chikuda H, Kugimiya F, Hoshi K, Ikeda T, Ogasawara T, Shimoaka T, et al. Cyclic GMP-dependent protein kinase II is a molecular switch from proliferation to hypertrophic differentiation of chondrocytes. Genes Dev. 2004;18(19):2418–29.
Koltes JE, Mishra BP, Kumar D, Kataria RS, Totir LR, Fernando RL, et al. A nonsense mutation in cGMP-dependent type II protein kinase (PRKG2) causes dwarfism in American Angus cattle. Proc Natl Acad Sci U S A. 2009;106(46):19250–5.
Bonnet C, Andrieux J, Beri-Dexheimer M, Leheup B, Boute O, Manouvrier S, et al. Microdeletion at chromosome 4q21 defines a new emerging syndrome with marked growth restriction, mental retardation and absent or severely delayed speech. J Med Genet. 2010;47(6):377–84.
Kawasaki Y, Kugimiya F, Chikuda H, Kamekura S, Ikeda T, Kawamura N, et al. Phosphorylation of GSK-3beta by cGMP-dependent protein kinase II promotes hypertrophic differentiation of murine chondrocytes. J Clin Invest. 2008;118(7):2506–15.
Kugimiya F, Chikuda H, Kamekura S, Ikeda T, Hoshi K, Ogasawara T, et al. Involvement of cyclic guanosine monophosphate-dependent protein kinase II in chondrocyte hypertrophy during endochondral ossification. Modern Rheumatol. 2005;15(6):391–6.
Agoston H, Khan S, James CG, Gillespie JR, Serra R, Stanton LA, et al. C-type natriuretic peptide regulates endochondral bone growth through p38 MAP kinase-dependent and -independent pathways. BMC Dev Biol. 2007;7:18.
Miyazawa T, Ogawa Y, Chusho H, Yasoda A, Tamura N, Komatsu Y, et al. Cyclic GMP-dependent protein kinase II plays a critical role in C-type natriuretic peptide-mediated endochondral ossification. Endocrinology. 2002;143(9):3604–10.
Yasoda A, Komatsu Y, Chusho H, Miyazawa T, Ozasa A, Miura M, et al. Overexpression of CNP in chondrocytes rescues achondroplasia through a MAPK-dependent pathway. Nat Med. 2004;10(1):80–6.
Gudi T, Hong GK, Vaandrager AB, Lohmann SM, Pilz RB. Nitric oxide and cGMP regulate gene expression in neuronal and glial cells by activating type II cGMP-dependent protein kinase. Faseb J. 1999;13(15):2143–52.
Gudi T, Huvar I, Meinecke M, Lohmann SM, Boss GR, Pilz RB. Regulation of gene expression by cGMP-dependent protein kinase. Transactivation of the c-fos promoter. J Biol Chem. 1996;271(9):4597–600.
Pilz RB, Broderick KE. Role of cyclic GMP in gene regulation. Front Biosci. 2005;10:1239–68.
da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.
Behrens A, Sibilia M, Wagner EF. Amino-terminal phosphorylation of c-Jun regulates stress-induced apoptosis and cellular proliferation. Nat Genet. 1999;21(3):326–9.
Hu T, Li C. Convergence between Wnt-beta-catenin and EGFR signaling in cancer. Mol Cancer. 2010;9:236.
Kitagawa M, Lee SH, McCormick F. Skp2 suppresses p53-dependent apoptosis by inhibiting p300. Mol Cell. 2008;29(2):217–31.
Konishi J, Yi F, Chen X, Vo H, Carbone DP, Dang TP. Notch3 cooperates with the EGFR pathway to modulate apoptosis through the induction of bim. Oncogene. 2010;29(4):589–96.
Mantamadiotis T, Lemberger T, Bleckmann SC, Kern H, Kretz O, Martin Villalba A, et al. Disruption of CREB function in brain leads to neurodegeneration. Nat Genet. 2002;31(1):47–54.
Pedraza LG, Stewart RA, Li DM, Xu T. Drosophila Src-family kinases function with Csk to regulate cell proliferation and apoptosis. Oncogene. 2004;23(27):4754–62.
Rosemann M, Kuosaite V, Kremer M, Favor J, Quintanilla-Martinez L, Atkinson MJ. Multilocus inheritance determines predisposition to alpha-radiation induced bone tumourigenesis in mice. Int J Cancer. 2006;118(9):2132–8.
Teo JL, Kahn M. The Wnt signaling pathway in cellular proliferation and differentiation: a tale of two coactivators. Adv Drug Deliv Rev. 2010;62(12):1149–55.
Wang G, Beier F. Rac1/Cdc42 and RhoA GTPases antagonistically regulate chondrocyte proliferation, hypertrophy, and apoptosis. J Bone Miner Res. 2005;20(6):1022–31.
Rangaswami H, Schwappacher R, Marathe N, Zhuang S, Casteel DE, Haas B, et al. Cyclic GMP and protein kinase G control a Src-containing mechanosome in osteoblasts. Sci Signal. 2010;3(153):ra91.
Jacks T, Remington L, Williams BO, Schmitt EM, Halachmi S, Bronson RT, et al. Tumor spectrum analysis in p53-mutant mice. Curr Biol. 1994;4(1):1–7.
Lowe SW, Cepero E, Evan G. Intrinsic tumour suppression. Nature. 2004;432(7015):307–15.
Massague J. G1 cell-cycle control and cancer. Nature. 2004;432(7015):298–306.
Olashaw N, Bagui TK, Pledger WJ. Cell cycle control: a complex issue. Cell Cycle. 2004;3(3):263–4.
Zaman F, Menendez-Benito V, Eriksson E, Chagin AS, Takigawa M, Fadeel B, et al. Proteasome inhibition up-regulates p53 and apoptosis-inducing factor in chondrocytes causing severe growth retardation in mice. Cancer Res. 2007;67(20):10078–86.
Sahni M, Ambrosetti DC, Mansukhani A, Gertner R, Levy D, Basilico C. FGF signaling inhibits chondrocyte proliferation and regulates bone development through the STAT-1 pathway. Genes Dev. 1999;13(11):1361–6.
Su WC, Kitagawa M, Xue N, Xie B, Garofalo S, Cho J, et al. Activation of Stat1 by mutant fibroblast growth-factor receptor in thanatophoric dysplasia type II dwarfism. Nature. 1997;386(6622):288–92.
Hurst-Kennedy J, Zhong M, Gupta V, Boyan BD, Schwartz Z. 24R,25-Dihydroxyvitamin D3, lysophosphatidic acid, and p53: a signaling axis in the inhibition of phosphate-induced chondrocyte apoptosis. J Steroid Biochem Mol Biol. 2010;122(4):264–71.
Mansfield K, Pucci B, Adams CS, Shapiro IM. Induction of apoptosis in skeletal tissues: phosphate-mediated chick chondrocyte apoptosis is calcium dependent. Calcif Tissue Int. 2003;73(2):161–72.
Teixeira CC, Mansfield K, Hertkorn C, Ischiropoulos H, Shapiro IM. Phosphate-induced chondrocyte apoptosis is linked to nitric oxide generation. Am J Physiol Cell Physiol. 2001;281(3):C833–9.
Zhong M, Carney DH, Ryaby JT, Schwartz Z, Boyan BD. Inhibition of phosphate-induced apoptosis in resting zone chondrocytes by thrombin peptide 508. Cells Tissues Organs. 2009;189(1–4):56–9.
Tataria M, Quarto N, Longaker MT, Sylvester KG. Absence of the p53 tumor suppressor gene promotes osteogenesis in mesenchymal stem cells. J Pediatr Surg. 2006;41(4):624–32. discussion 624–632.
Kulikov R, Boehme KA, Blattner C. Glycogen synthase kinase 3-dependent phosphorylation of Mdm2 regulates p53 abundance. Mol Cell Biol. 2005;25(16):7170–80.
Xue Y, Zhou F, Zhu M, Ahmed K, Chen G, Yao X. GPS: a comprehensive www server for phosphorylation sites prediction. Nucleic Acids Res. 2005;33(Web Server issue):W184–7.
Mori R, Kodaka T, Soeta S, Sato J, Kakino J, Hamato S, et al. Preliminary study of histological comparison on the growth patterns of long-bone cortex in young calf, pig, and sheep. J Vet Med Sci. 2005;67(12):1223–9.
Garbe JR, Elsik CG, Antoniou E, Reecy JM, Clark KJ, Venkatraman A, et al. Development and application of bovine and porcine oligonucleotide arrays with protein-based annotation. J Biomed Biotechnol. 2010;2010:453638.
Skibbe DS, Wang X, Zhao X, Borsuk LA, Nettleton D, Schnable PS. Scanning microarrays at multiple intensities enhances discovery of differentially expressed genes. Bioinformatics. 2006;22(15):1863–70.
Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J, et al. Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Res. 2002;30(4):e15.
Storey JD, Tibshirani R. Statistical significance for genomewide studies. Proc Natl Acad Sci U S A. 2003;100(16):9440–5.
Dennis Jr G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, et al. DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003;4(5):3.
KEGG: Kyoto Encyclopedia of Genes and Genomes. [http://www.genome.jp/kegg/pathway.html]
Livak KJST. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. Methods. 2001;25(4):402–8.
The authors acknowledge M. Ruble, R. Berryman, and the ISU farm crews for animal care & handling; C. Tuggle for use of Pathway Studio software; the ISU DNA facility & the Schnable lab for microarray hybridization chamber & scanner use. DK and RSK thankfully acknowledge the financial support received as overseas fellowship from DBT, Govt. of India. Additionally, we thank M. Mayes, J.R. Tait, J. Frerichs, M. Schneider, B. Mote for assistance with tissue collection; D. Nettleton, D. Wang, M. Rothschild for statistical advice; Y. Cheng, R. Wagner for technical assistance. This journal paper of the Iowa Agriculture and Home Economics Experiment Station, Ames, Iowa (Project No. NRSP-8) was supported by Hatch Act and State of Iowa funds.
JMR and JEK hold a patent for a genetic test to identify R678X PRKG2 affected and carrier individuals. The authors declare they have no competing interest.
JEK performed the microarray experiments, statistical analysis, enrichment analyses, real-time PCR and conceived and designed the experiments. DK and RSK developed the methodology and assisted with RNA isolation from growth plate. VC provided professional assistance in sample collection and identification of proper tissues during dissection of samples. JMR provided conceptional advice and guidance and oversaw all work in the project. All authors assisted in writing and editing of the paper, with JEK taking the lead. All authors read and approved the final manuscript.
List of annotations for all putative DE genes and enriched pathways. Additional file 1 contains all probe and annotation information for the cDNA microarray used in this study. Specifically, the data is split into worksheets by the specific analysis and annotation information. Each work sheet contains probe IDs differentially expressed at p<0.05 at one of three different PMT measures. SwissProt IDs are the human or mouse IDs for the corresponding bovine transcript. Gene names from DAVID are the corresponding gene names from Swissprot as Swissprot IDs were used to run the DAVID analysis.
Gene Ontology (GO) enrichment results. Additional file 2 contains all GO enrichment results for the DE genes in the PMT70 and combined datasets.
Pathway studio enrichment results divided by up and down regulated genes in PRKG2 null (dwarf) animals. Additional file 3 contains all Pathway studio pathway enrichment results. It is divided into worksheets by the two analysis methods detecting direct or common regulators and the entities which provide descriptions of the curated regulators. In addition, the data is divided by which genes were up regulated in the wild-type (PRKG2 +/-) or the dwarf (PRKG2 -/-) genotypes.