- Research article
- Open Access
Alternative splicing in osteoclasts and Paget’s disease of bone
BMC Medical Genetics volume 15, Article number: 98 (2014)
Mutations in the SQSTM1/p62 gene have been reported in Paget’s disease of bone (PDB), but they are not sufficient to induce the pagetic osteoclast (OC) phenotype. We hypothesized that specific RNA isoforms of OC-related genes may contribute to the overactivity of pagetic OCs, along with other genetic predisposing factors.
Alternative splicing (AS) events were studied using a PCR-based screening strategy in OC cultures from 29 patients with PDB and 26 healthy donors (HD), all genotyped for the p62P392L mutation. Primer pairs targeting 5223 characterized AS events were used to analyze relative isoform ratios on pooled cDNA from samples of the four groups (PDB, PDBP392L, HD, HDP392L). Of the 1056 active AS events detected in the screening analysis, 192 were re-analyzed on non-amplified cDNA from each subject of the whole cohort.
This analysis led to the identification of six AS events significantly associated with PDB, but none with p62P392L. The corresponding genes included LGALS8, RHOT1, CASC4, USP4, TBC1D25, and PIDD. In addition, RHOT1 and LGALS8 genes were upregulated in pagetic OCs, as were CASC4 and RHOT1 genes in the presence of p62P392L. Finally, we showed that the proteins encoded by LGALS8, RHOT1, USP4, TBC1D25, and PIDD were expressed in human OCs.
This study allowed the identification of hitherto unknown players in OC biology, and our findings of a differential AS in pagetic OCs may generate new concepts in the pathogenesis of PDB.
Osteoclasts (OCs) are multinucleated cells, whose phenotype includes the expression of the calcitonin receptor (CTR) and of the Receptor Activator of NF-κB (RANK), and the singular ability to resorb bone ,. OC formation, activity and survival are regulated by two major signaling pathways, which are activated by M-CSF and RANKL (Receptor Activator of NK-κB Ligand).
Paget’s disease of bone (PDB) is characterized by focal and disorganized increases in bone turnover, and OCs have been identified as the cells primarily affected in PDB . Pagetic OCs are both larger and more numerous than those in healthy individuals. They are overactive and hypersensitive to osteoclastogenic factors , and are resistant to apoptosis . The SQSTM1 gene encodes the ubiquitin-binding protein sequestosome 1, also known as p62. Mutations in the SQSTM1 gene have been identified in a high proportion of PDB patients , the p62P392L substitution being the most frequent . In PDB, p62P392L contributes at least in part to the induction of an activated stage in OCs by stimulating signaling pathways that can lead to NF-κB activation ,. In vivo, knock-in mice expressing p62P394L (the murine equivalent of human p62P392L) develop focal osteolytic lesions, some of which resemble PDB lesions , although in another study, co-expression of the measles virus nucleocapsid gene in the OC lineage was required . Therefore SQSTM1 mutations may not be sufficient to induce the pagetic OC phenotype, and environmental factors such as viral infection may contribute . Two recently published genome-wide association studies (GWAS) in PDB patients free of SQSTM1 mutations identified seven significant genetic variants for susceptibility to PDB located at the 1p13 (CSF1), 7q33 (CNOT4, NUP205, SLC13A4), 8q22 (TM7SF4), 10p13 (OPTN), 14q32 (RIN3), 15q24 (PML, GOLGA6A), and 18q21 (TNFRSF11A) loci ,. This implies that genes other than SQSTM1 may also contribute to the pathogenesis of PDB, although so far only disease-causing mutations in SQSTM1 have been identified. Besides gene expression modulation or mutations, other mechanisms regulate phenotypic variations in cells through post-transcriptional processes such as alternative splicing (AS) that leads to a vast array of transcripts with diverse functions. We hypothesized that specific RNA isoforms of OC-related genes may contribute to the pagetic OC phenotype. Our objective was therefore to look for PDB-specific AS events in these cells, and to investigate the impact of the p62P392L mutation.
Opti Eagle’s minimum essential media (Opti-MEM), penicillin, streptomycin, fungizone, glutamine, and fetal calf serum (FCS) were purchased from Wisent (Montreal, QC). Ficoll–Paque was purchased from Amersham Biosciences (Montreal, QC). Human recombinant (hr) M-CSF, and hrGM-CSF were purchased from R&D (R&D Systems, Minneapolis, MN); soluble hrRANKL was produced in our laboratory. Rabbit polyclonal antibodies against human PIDD (#ab78389), Galectin 8 (#ab41649), and RHOT1 (#ab83779) were purchased from Abcam (Cambridge, MA); rabbit polyclonal antibodies against human TBC1D25 (OATL1) (#HPA029197), and USP4 (#U0635) were purchased from Sigma-Aldrich (St. Louis, MO), and fluorescent Alexa antibodies, Di Aminido Phenyl lndol (DAPI), and siRNA from Invitrogen (Burlington, ON).
Clinical investigation and phenotype classification
Blood samples from each of the healthy donors (HD) and PDB patients were provided by authors JPB and LM. This research has been approved by the Ethics Committees of the CHU de Québec. All participants signed an informed-consent document before entering the study. Phenotypic expression of PDB had been previously established (JPB) in all participants based on: total serum alkaline phosphatase, total body bone scan, skull and enlarged view pelvis X-rays, and if needed, additional X-rays to confirm PDB, as described . Exons 7 and 8 of the gene encoding SQSTM1/p62 had been previously sequenced in every donor and patient . The p62P392L mutation was the only SQSTM1 mutation identified in individuals participating in this study. Fifty-five participants, 29 PDB (13 females, 16 males) and 26 HD (16 females, 10 males), were divided into four groups: healthy donors exempt from any known mutation in the SQSTM1/p62 gene (HDwt); healthy donors carrying the p62P392L mutation (HDP392L), these donors came from families with PDB history (carrying the p62P392L mutation), but had no phenotypic expression of PDB at the time of study; PBD patients exempt from any known mutation in the SQSTM1/p62 gene (PDBwt); PDB patients carrying the p62P392L mutation (PDBP392L).
Peripheral blood mononuclear cells (PBMCs) were isolated from heparinized blood by density-gradient centrifugation, and suspended in Opti-MEM with the antibiotics, glutamine, and 2% FCS. After incubating overnight, the cells were gently washed to remove any non-adherent cells. The selected PBMCs were cultured in Opti-MEM supplemented with GM-CSF (100 pg/ml) for the first 3 days, and then for a further 3 weeks in the same medium supplemented with M-CSF (25 ng/ml) and RANKL (75 ng/ml). The medium was changed every 2–3 days. Alternatively, blood was harvested from human umbilical cord at delivery after obtaining informed consent from parturient women, as approved by our institution’s review board (University of Sherbrooke). Cord blood monocytes (CBMs) were isolated and processed as described above. These culture conditions generate multinucleated cells (MNC) that express OC markers and have the ability to resorb bone ,. Human embryonic kidney-293 T (293 T) cells were cultured in DMEM containing 10% FCS, glutamine and antibiotics.
Total cellular RNA was extracted from fully matured OCs, reverse transcribed, and analyzed by high-throughput PCR amplification at the Université de Sherbrooke RNomics Platform as previously described ,. Five (5) ng of total RNA was used for each PCR experiment. For the detection screening we used RNA pre-amplified using a linear isothermal RNA amplification (Transplex Whole Transcriptome Amplification Kit, Sigma, Markham, ON) following the manufacturer’s protocol. AS events were characterized by end-point PCR. Primers were designed to flank the AS events, such that following amplification and analysis by microcapillary electrophoresis on Caliper LC-90 instruments (Caliper Life Sciences, Hopkinton, MA), the relative ratio of the isoforms can be deduced . AS events which were amenable to characterization by high-thoughput PCR, that is, whose isoform sizes differ by between 10 and 450 nucleotides at a particular event, were selected from the RefSeq database . Automated querying of this highly curated database has identified a genome-wide selection of 5223 AS events which fit the selection criteria from the full set of over 20 K human gene entries .
In the first step of the screen, termed the “detection screen”, primer pairs targeting the 5223 AS events were used to perform PCR and microcapillary analysis on pre-amplified RNA in representative samples of the four groups, each sample being a pooled RNA from five patients (60 ng of total RNA per individual needed prior to pre-amplification). The purpose of this screen was to detect AS events in the sample pools for which both isoforms were detectable to a minimum of 10% of the total signal in at least one of the samples. Such AS events are designated as “active”. If only one form of the event is detected in all samples, or if no amplification of the expected products occurs, these AS events were not considered for subsequent validation.
Next, during what we term the “validation screen”, we analyzed the active AS events in PCR experiments on non-amplified RNA in our cohort of 55 individuals, using the same primers, however in this screen, the relative amounts of long and short isoforms issued from the microcapillary analysis of the amplification reaction was considered. The splicing profile from each PCR reaction was expressed as a percent splicing index (PSI or Ψ) which is the ratio of the long isoform concentration to the total (long + short) isoform concentrations, expressed as a percentage. The Ψ values of the active AS events were compared in the four groups. The use of multiple samples reduces the impact of individual aberrant AS, artefacts related to the pre-amplification procedure, and to confirm the PDB-specific events.
Gene expression study
For the quantitative real-time RT-PCR analysis, total RNA was extracted from fully matured OCs using the RNeasy plus kit (Qiagen, Mississauga, Canada), including a DNase digestion step. At least 50 μg of RNA was harvested from each cell sample. After being assayed and checked for quality, 1 μg of total RNA was reverse transcribed using the high-efficiency QuantiTect Primer Assay kit (Qiagen), and sent for quantitative PCR amplification (RNomics Platform, University of Sherbrooke). Human primers of candidate genes were generated and validated, and the real-time PCR reaction was conducted with 200 nM of primers, 10 ng of cDNA and 5 μl of Power SYBR-green Master Mix in a total volume of 10 μl. Amplification and detection of the candidate genes and of three reference genes (Splicing factor 3a, subunit 1 (SF3A1), Proteasome 26S subunit ATPase 4 (PSMC4), Pumilio homolog 1 (PUM1)) were conducted with a Realplex 2 Master Cycler (Eppendorf, Mississauga, Canada).
We investigated the relative expression of genes for which AS events were significantly associated with PDB in OC cultures. All the primers were designed based on sequences reported in the Aceview database . Three primer pairs were designed for each gene, in order to choose the best primer pair after validation. All primers pairs were validated by qPCR amplification of a total XPressRef Universal human mRNA (Cedarlane, Burlington, ON). Amplification products were submitted to capillary electrophoresis in order to validate the size of the expected product. In addition, a melting curve was generated and only those primers giving one product were selected. If all the three primer pairs gave such results, the pair giving the lowest Ct was selected. Finally, the primer pair efficiency was obtained from a standard curve experiment where a series of dilution of the same sample was correlated to the Ct values.
All samples were run in triplicate for the target and reference genes. Relative expression levels were normalized with respect to a set of reference primer pairs (i.e. the geometric mean of the three reference genes). Relative expression quantification analysis relied on the qBase method . This method constitutes an improvement over the classical delta-delta-Ct method, which was extended to take into account multiple stably expressed reference genes for improved normalization.
Seventy-two to 96 hours before lysing cells, 293 T cells were transfected either with a specific siRNA (50 nM), or with a negative control siRNA (scrambled siRNA) in a solution containing Lipofectamine 2000 Transfection Reagent in serum-free Opti-MEM. One specific siRNA was used to inhibit each gene (PIDD, TBC1D25, RHOT1, LGALS8), two siRNAs against USP4 were used in combination (Table 1). The cells were incubated at 37°C for 2 hours. The subsequent downregulation of all targeted proteins was assessed by western-blotting.
Cells were cultured as described above. Proteins were extracted and lysates were quantified by modified Bradford assay (Bio Rad, Mississauga, ON) and run through SDS-PAGE. Western-blots were performed by incubating a primary antibody against Galectin 8, RHOT1, USP4, TBC1D25, PIDD, or β-actin where appropriate overnight at 4°C. HRP-conjugated secondary antibodies were used to achieve detection with a chemiluminescent system.
At the end of the OC cultures, cells were washed quickly with cold PBS, fixed with 1% para-formaldehyde, and then permeabilized in PBS-Triton X-100 (0.2%). Specific antibodies directed against USP4, TBC1D25, and PIDD were incubated in Antibody diluent (DAKO) overnight at 4°C. Alexa-488 (green) anti-rabbit antibodies were incubated for 2 h at room temperature, along with Alexa Fluor 633 phalloidin. DAPI counter-staining was performed, so that OCs containing more than three nuclei could be visualized. Sequences of pictures were taken with appropriate filters to show all three colors, and superimposed with Simple PCI software.
We studied four groups (PDBwt, PDBP392L, HDwt, HDP392L) in six intergroup analyses (all PDB vs all HD; all PDB patients vs HDwt; PDBwt vs HDwt, PDBwt vs HDP392L, PDBwt vs PDBP392L; and HDwt vs HDP392L). Student’s t-tests were applied to the Ψ values for each AS event in each group and Ψ value changes with values of p < 0.05 were considered to be significant, with a correction for multiple testing using the false discovery rate or q-value . For genes in which we confirmed the presence of specific AS events of PDB OCs, we looked for association with the P392L mutation using Fisher exact tests. For gene expression analysis, the fold change was calculated as the ratio of the mean relative expression for each group, and results were expressed as the base 2 logarithm of the fold change. Statistical significance was defined as p < 0.05 using a Student’s t-test.
Alternative splicing in Paget’s disease of bone
To perform our study of AS in PDB, we used a validated model for the in-vitro study of human OCs involving peripheral blood mononuclear cells (PBMCs) . We studied OC cultures from a cohort of patients (PDB) and healthy donors (HD), all genotyped for the p62P392L mutation, and representing four distinct groups: HDwt (n = 16, mean age 61 years; 39–84), HDP392L (n = 10, mean age 51.3 years; 33–70; of which three were under 52 years old), PDBwt (n = 17, mean age 70.7 years; 45–86), PDBP392L (n = 12, mean age 75 years; 59–85). In the first step of the AS screen, termed the “detection screen”, primer pairs targeting 5223 characterized AS events were used to amplify and evaluate the relative isoform ratios on pooled cDNA from our four groups (5 patients in each group). 1056 of the targeted AS events showed evidence of regulation, and were thus classified as “active”. Three genes of the 5223 tested AS events (VDR, OPTN, and PML) had been previously associated with PDB ,, of which OPTN and PML displayed active AS events.
From this active list (n = 1056), 192 AS events (in 164 genes) were selected based on genes that have a confirmed or potential function in the OC phenotype, function or signaling pathways, or those whose involvement in PDB has been suggested ,, (Additional file 1: Table S1). The 192 active AS events were then tested in a “validation screen”. Five individuals from each group were analyzed separately using the computed Ψ values (see Methods) to compare profiles between groups. This screen yielded 21 AS events that were either significantly associated with PDB or with the p62P392L mutation (t-test, p < 0.05, Additional file 2: Table S2). We then confirmed this association by analyzing these 21 active AS events by PCR on cDNA from the remaining non-tested individuals (12 PDBwt, 7 PBDP392L, 11 HDwt, 5 HDP392L), and Ψ values were compared (Table 2).
This led to the identification of six AS events that were significantly associated with PDB, corresponding to: LGALS8 (galectin-8), RHOT1 (mitochondrial Rho GTPase1), CASC4 (cancer susceptibility candidate 4), USP4 (ubiquitin specific protease 4), TBC1D25 (TBC1 domain family, member 25), and PIDD (p53-induced death domain protein) (Table 3). No significant association of AS profile with the p62P392L mutation was found.
Expression of genes with PDB-related AS events in OCs
The OC gene expression profiles of the six genes LGALS8, RHOT1, CASC4, USP4, TBC1D25, and PIDD, were analyzed in the four groups (12 PDBwt, 10 PBDP392L, 12 HDwt, 10 HDP392L), and the results in PBD patients relative to healthy donors, in presence of the p62P392L mutation or not, are shown in Figure 1. In the comparison between the gene expression (mean +/− SEM) in all PDB patients regardless of whether the p62P392L mutation was present or not versus non-mutated healthy donors, there was a significant up-regulation in the RHOT1 gene (1.06 ± 0.04 in PBD vs 0.85 ± 0.03 in HDwt, p =0.002), and to a lesser extent in the LGALS8 gene (1.07 ± 0.03 vs 0.94 ± 0.03, p =0.015) in PDB patients. When comparing the relative gene expression between OCs from PDB patients carrying p62P392L and PDB patients with no mutation in the SQSTM1 gene, we observed a significant increase in the expression of the CASC4 gene in presence of the mutation (1.34 ± 0.09 in PDBP392L vs 1.01 ± 0.04 in PDBwt, p =0.002), and in the expression of RHOT1 (1.17 ± 0.06 vs 0.97 ± 0.04, p =0.021). Finally, the relative expression of two genes was significantly upregulated in OC cultures from HD carriers of p62P392L compared to healthy donors with no mutation in the SQSTM1 gene, with an upregulation of CASC4 (1.39 ± 0.40 in HDP392L vs 1.05 ± 0.11 in HDwt, p = 0.015), and RHOT1 genes (1.02 ± 0.06 in HDP392L vs 0.85 ± 0.06 in HDwt, p = 0.026) (Figure 1).
Protein expression of PIDD, TBC1D25, RHOT1, LGALS8 and USP4 in mature OCs
A major interest of this study is that it enabled us to identify AS events from genes not previously known to be associated with OC biology. We therefore studied the expression of these genes in human OCs at the protein level, focusing on five genes (PIDD, TBC1D25, RHOT1, LGALS8 and USP4), as they were previously reported for their involvement in NF-κB signaling, apoptosis, autophagy, and ubiquitination in other systems, which are collectively recognized as pathways essential to OC biology . As we were unable to find functioning antibodies against CASC4, it was excluded from further study. To evaluate the expression of the encoded proteins, we used OCs derived from cord blood monocytes (CBMs), and investigated the protein-level expression by immunofluorescence and western blot using specific antibodies. Our objective was to confirm the expression of the proteins encoded by the selected genes in human OCs, for which the CBM-derived OCs represent a reliable model ,,.
Western blotting showed that all the proteins encoded by the PIDD, TBC1D25, RHOT1, USP4, and LGALS8 genes were indeed expressed in human OC cultures (Figure 2). To identify and validate some of the isoforms of the candidate proteins observed in OCs, we also analyzed the effects of siRNA-mediated depletion of these genes (Table 1) at the protein level. 293 T cells were used for this as OCs cannot be transfected so efficiently. The antibodies used were all directed against a peptide shared by all the described isoforms for each of the proteins, except for anti-TBC1D25 antibodies, which are directed against a peptide within isoform 1 but containing only 25% of the isoform 2 sequence. PIDD isoform 4 (70 kDa) was the main isoform expressed by OCs, and its expression was confirmed by siRNA knockdown. Two other bands of higher MW (130 kDa and 85 kDa) were also detected in OCs and may correspond to the other isoforms (1/2 and 3, respectively), but were not clearly down-regulated by siRNA. The expression of the long isoform 1 of TBC1D25 was specifically detected (80 kDa). Two isoforms of RHOT1 were primarily detected in OCs (55 and 65 kDa), but we were not able to differentiate between specific isoforms as the observed bands migrated faster than predicted. Two proteins of high MW (95 and 110 kDa) were detected in the USP4 analysis, which may correspond to the two known isoforms. Only one Galectin 8 isoform of 42 kDa was specifically detected (Figure 2).
The same protein profile was observed in OC cultures from non-mutated PDB or HD, and we were not able to quantify the expression of the short and long isoforms (as defined in Table 3), as the long and short isoforms of four candidate proteins (PIDD, RHOT1, USP4, and galectin-8) showed similar migration rates. The relative expression of TBC1D25 long isoform 1 (the short isoform being non-detectable) was significantly lower in PDB OCs compared to HD OCs (p =0.009) (Figure 3).
In order to compensate for effects potentially attributable to non-osteoclastic cells included in the populations used for western blotting, we included single-cell analysis using immunofluorescence to evaluate the expression of the PIDD, TBC1D25, and USP4 encoded proteins in human OCs. We found that mature multinucleated OCs expressed PIDD, predominantly in the nucleus. The expression of TBC1D25 and USP4 was cytoplasmic in OCs (Figure 4).
AS plays a central role in protein diversity and post-transcriptional gene regulation. The role of AS in the production of oncogenes and tumor suppressors is of major interest, as variants of these genes are often found specifically in tumors from diverse tissues ,. In non-neoplastic diseases, the functional importance of AS is best demonstrated by its capacity to cause diseases such as muscular dystrophies and premature-aging disorders . In spite of its imposing biological relevance, little is known about the effects of AS isoforms in bone diseases. Splicing site SQSTM1 mutations in PDB have been reported at least twice ,. Alternative splicing mutations have also been reported in other bone diseases, such as in TCIRG1-linked autosomal recessive osteopetrosis , however the present study is the first to investigate AS in OCs in PDB. A major strength of our approach lies in using human OCs derived from circulating OC precursors from PBD patients and healthy controls. Culture of OCs in vitro may help to identify the genetic contribution that may cause intrinsic alterations in OCs, by comparison with appropriate normal donors. We hypothesized that specific RNA splicing of OC-related genes may contribute to the OC phenotype and the clinical expression of PDB. Using a screening strategy that has been applied to cancer research in our institution , we found AS events significantly associated with PDB, but not with p62P392L, in six genes: LGALS8, RHOT1, CASC4, USP4, TBC1D25, and PIDD. While gene selection biases may have been present in both the detection and validation screens, as, for economic reasons, not all known AS events were screened, the identification of the six candidate genes from this small subset underscores potentially significant role of AS in the development and progression of PDB. Upregulation of LGALS8 and RHOT1 genes was observed in pagetic OCs, and the p62P392L mutation was associated with an upregulation of RHOT1 and CASC4. The encoded proteins were all detected in immunoblotting analyses in OC cultures, and three of them (PIDD, TBC1D25 and USP4) were expressed in human mature OCs. While not all isoforms have been clearly discriminated by western blot analysis, thus limiting our capacity to examine isoform ratios between PBD and HD at the protein level, subtle changes in the splice-isoform ratios might be associated with disease susceptibility . In addition to providing evidence of an AS modulation in PDB, our findings also enabled us to identify AS events from genes not previously known to be associated with OC biology.
In OCs, the p62 scaffolding protein is one of the functional links reported between RANKL and TNFR-associated factor 6 (TRAF6)-mediated NF-κB activation ,. Protein p62 promotes the binding of CYLD to TRAF6 , a de-ubiquitinating enzyme (DUB) that negatively regulates NF-κB activity by reducing TRAF6 auto-ubiquitination . The p62P392L mutation abolished CYLD interaction and enhanced OC formation and activity . USP4 (Ubiquitin-specific peptidase 4) is another DUB, of which two isoforms are produced by AS in humans . Its functional properties include the de-ubiquitination of TAK1 , an important kinase in the RANKL-induced signaling complex . In HEK293 cells, USP4 has also been shown to associate with TRAF6 and to inhibit TRAF6 ubiquitination, and consequently NF-κB activation . USP4 could therefore play a role in the p62/TRAF6-related signaling in pagetic OCs.
Pagetic OCs are overactive and resistant to apoptosis, a phenotype partly related to p62P392L. We found a differential splicing in RNA of genes involved in apoptosis (CASC4, RHOT1, PIDD). The increased expression of CASC4 has been associated with HER-2/neu proto-oncogene overexpression, a membrane tyrosine kinase linked to cell proliferation and survival pathways . RHOT1, also named MIRO-1, is a Rho GTPase involved in mitochondrial transport. Overexpression of its constitutively active mutant resulted in an increased apoptotic rate . LGALS8 encodes Galectin-8, a beta-galactoside-binding lectin which has been implicated in cell spreading ,, and in the modulation of neutrophil functions . Galectin-8 also appears as a potent pro-apoptotic agent in T cells and inflammatory cells ,, and is a strong inducer of ERK activation ,. In addition to USP4 that may modulate p53 degradation, PIDD, also named LRDD (leucine-rich repeats and death domain containing), is an effector of p53-induced apoptosis. PIDD is able to induce NF-κB activation , but is also able to form a caspase- 2-activating platform to cause cell death . PIDD isoform 1 is the only isoform able to activate caspase-2, whereas isoforms −1, −2, −3 activate NF-κB . Isoform 4 is regulated by p53 and when overexpressed independently promotes apoptosis . While isoform 4 did not show significant differences between PDB and non-PDB samples, and isoform 2 did not pass the detection screen selection criteria, we found a slight decrease in isoform 1 relative to isoform 3 in PDB samples, implying a more anti-apoptotic behaviour of PIDD in PDB than non-PDB samples. None of these genes (PIDD, RHOT1, LGALS8) has ever been associated with PDB, and all represent good candidates for further investigation of the OC survival in PDB.
A pathognomonic feature found in PDB patients is the presence of inclusion bodies in OCs within affected bone; because these inclusion bodies resemble the p62-aggregates observed in diseases involving defective autophagy, dysregulation of the autophagy process may well be part of the pathogenesis of PDB . A dysregulation of the autophagy process has been suggested in OCs from knock-in mice expressing p62P394L, and an impaired autophagosomal degradation of p62 has been observed in HEK 293 cell lines transfected with PDB-related p62 variants . OATL1/TBC1D25, a Rab-GTPase-activating protein (Rab-GAP), of which two spliced isoforms exist, is a newly-discovered binding partner of Atg8/LC3, and has been shown to be involved in a late stage of autophagosome maturation in MEFs and COS-7 cells . Our results have identified TBC1D25 as a potential candidate that may be worth investigating in PDB.
A limitation of our study is that expression of the identified genes and their RNA isoforms occurred in cells that were analyzed in long-term in-vitro cultures. These conditions may have been subject to stress-induced aberrant gene expression, such as cancer susceptibility CASC4, or genes related to p53 which is induced by DNA damage or stress stimuli such as PIDD, and USP4. Yet, the proteins encoded by these genes were found to be expressed in human OCs, and their splicing profiles were significantly modified in OCs from PBD patients compared to controls. In addition, except for the PIDD isoforms, the functional impact of the different splice variants of the six genes reported here remain elusive.
Our screening strategy led to the identification of AS events that were significantly associated with pagetic OCs, and our findings allowed the identification of hitherto unknown players in OC biology. Further studies are warranted to identify the exact function of the LGALS8, RHOT1, CASC4, USP4, TBC1D25, and PIDD genes in human OCs, to investigate the regulation and functional impact of the identified specific isoforms in OCs using isoform-specific inhibitors, and their respective contribution to the OC phenotype in PDB. Our results may identify future avenues for research, as AS modulation and its transposition to OC biology represent a new area of research with wide-ranging potential interest for explaining how these cells are regulated. This new approach may also generate new concepts in the pathogenesis of Paget’s disease.
Boyle WJ, Simonet WS, Lacey DL: Osteoclast differentiation and activation. Nature. 2003, 423 (6937): 337-342. 10.1038/nature01658.
Quinn JM, Morfis M, Lam MH, Elliott J, Kartsogiannis V, Williams ED, Gillespie MT, Martin TJ, Sexton PM: Calcitonin receptor antibodies in the identification of osteoclasts. Bone. 1999, 25 (1): 1-8. 10.1016/S8756-3282(99)00094-0.
Roodman GD, Windle JJ: Paget disease of bone. J Clin Invest. 2005, 115 (2): 200-208. 10.1172/JCI24281.
Kurihara N, Hiruma Y, Yamana K, Michou L, Rousseau C, Morissette J, Galson DL, Teramachi J, Zhou H, Dempster DW, Windle JJ, Brown JP, Roodman GD: Contributions of the measles virus nucleocapsid gene and the SQSTM1/p62(P392L) mutation to Paget’s disease. Cell Metab. 2011, 13 (1): 23-34. 10.1016/j.cmet.2010.12.002.
Chamoux E, Couture J, Bisson M, Morissette J, Brown JP, Roux S: The p62 P392L mutation linked to Paget’s disease induces activation of human osteoclasts. Mol Endocrinol. 2009, 23 (10): 1668-1680. 10.1210/me.2009-0066.
Chung PY, Van Hul W: Paget’s disease of bone: evidence for complex pathogenetic interactions. Semin Arthritis Rheum. 2012, 41 (5): 616-641. 10.1016/j.semarthrit.2011.07.005.
Morissette J, Laurin N, Brown JP: Sequestosome 1: mutation frequencies, haplotypes, and phenotypes in familial Paget’s disease of bone. J Bone Miner Res. 2006, 21 (Suppl 2): P38-P44. 10.1359/jbmr.06s207.
Chamoux E, McManus S, Laberge G, Bisson M, Roux S: Involvement of kinase PKC-zeta in the p62/p62(P392L)-driven activation of NF-kappaB in human osteoclasts. Biochim Biophys Acta. 2013, 1832 (3): 475-484. 10.1016/j.bbadis.2012.12.008.
Daroszewska A, Hof RJ V ’t, Rojas JA, Layfield R, Landao-Basonga E, Rose L, Rose K, Ralston SH: A point mutation in the ubiquitin-associated domain of SQSMT1 is sufficient to cause a Paget’s disease-like disorder in mice. Hum Mol Genet. 2011, 20 (14): 2734-2744. 10.1093/hmg/ddr172.
Albagha OM, Visconti MR, Alonso N, Langston AL, Cundy T, Dargie R, Dunlop MG, Fraser WD, Hooper MJ, Isaia G, Nicholson GC, del Pino MJ, Gonzalez-Sarmiento R, di Stefano M, Tenesa A, Walsh JP, Ralston SH: Genome-wide association study identifies variants at CSF1, OPTN and TNFRSF11A as genetic risk factors for Paget’s disease of bone. Nat Genet. 2010, 42 (6): 520-524. 10.1038/ng.562.
Albagha OM, Wani SE, Visconti MR, Alonso N, Goodman K, Brandi ML, Cundy T, Chung PY, Dargie R, Devogelaer JP, Falchetti A, Fraser WD, Gennari L, Gianfrancesco F, Hooper MJ, Van Hul W, Isaia G, Nicholson GC, Nuti R, Papapoulos S, Montes Jdel P, Ratajczak T, Rea SL, Rendina D, Gonzalez-Sarmiento R, Di Stefano M, Ward LC, Walsh JP, Ralston SH: Genome-wide association identifies three new susceptibility loci for Paget’s disease of bone. Nat Genet. 2011, 43 (7): 685-689. 10.1038/ng.845.
Laurin N, Brown JP, Lemainque A, Duchesne A, Huot D, Lacourciere Y, Drapeau G, Verreault J, Raymond V, Morissette J: Paget disease of bone: mapping of two loci at 5q35-qter and 5q31. Am J Hum Genet. 2001, 69 (3): 528-543. 10.1086/322975.
Laurin N, Brown JP, Morissette J, Raymond V: Recurrent mutation of the gene encoding sequestosome 1 (SQSTM1/p62) in Paget disease of bone. Am J Hum Genet. 2002, 70 (6): 1582-1588. 10.1086/340731.
Roux S, Lambert-Comeau P, Saint-Pierre C, Lepine M, Sawan B, Parent JL: Death receptors, Fas and TRAIL receptors, are involved in human osteoclast apoptosis. Biochem Biophys Res Commun. 2005, 333 (1): 42-50. 10.1016/j.bbrc.2005.05.092.
Klinck R, Bramard A, Inkel L, Dufresne-Martin G, Gervais-Bird J, Madden R, Paquet ER, Koh C, Venables JP, Prinos P, Jilaveanu-Pelmus M, Wellinger R, Rancourt C, Chabot B, Abou Elela S: Multiple alternative splicing markers for ovarian cancer. Cancer Res. 2008, 68 (3): 657-663. 10.1158/0008-5472.CAN-07-2580.
Venables JP, Klinck R, Koh C, Gervais-Bird J, Bramard A, Inkel L, Durand M, Couture S, Froehlich U, Lapointe E, Lucier JF, Thibault P, Rancourt C, Tremblay K, Prinos P, Chabot B, Elela SA: Cancer-associated regulation of alternative splicing. Nat Struct Mol Biol. 2009, 16 (6): 670-676. 10.1038/nsmb.1608.
Thierry-Mieg D, Thierry-Mieg J: AceView: a comprehensive cDNA-supported gene and transcripts annotation. Genome Biol. 2006, 7 (Suppl 1): S12 11-S12 14. 10.1186/gb-2006-7-s1-s12.
Hellemans J, Mortier G, De Paepe A, Speleman F, Vandesompele J: qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol. 2007, 8 (2): R19-10.1186/gb-2007-8-2-r19.
Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proc Natl Acad Sci U S A. 2003, 100 (16): 9440-9445. 10.1073/pnas.1530509100.
Michou L, Chamoux E, Couture J, Morissette J, Brown JP, Roux S: Gene expression profile in osteoclasts from patients with Paget’s disease of bone. Bone. 2010, 46 (3): 598-603. 10.1016/j.bone.2009.11.012.
McManus S, Roux S: The adaptor protein p62/SQSTM1 in osteoclast signaling pathways. J Mol Signal. 2012, 7: 1-10.1186/1750-2187-7-1.
Hodge JM, Kirkland MA, Aitken CJ, Waugh CM, Myers DE, Lopez CM, Adams BE, Nicholson GC: Osteoclastic potential of human CFU-GM: biphasic effect of GM-CSF. J Bone Miner Res. 2004, 19 (2): 190-199. 10.1359/JBMR.0301232.
Christofk HR, Vander Heiden MG, Harris MH, Ramanathan A, Gerszten RE, Wei R, Fleming MD, Schreiber SL, Cantley LC: The M2 splice isoform of pyruvate kinase is important for cancer metabolism and tumour growth. Nature. 2008, 452 (7184): 230-233. 10.1038/nature06734.
Venables JP: Aberrant and alternative splicing in cancer. Cancer Res. 2004, 64 (21): 7647-7654. 10.1158/0008-5472.CAN-04-1910.
Tazi J, Bakkour N, Stamm S: Alternative splicing and disease. Biochim Biophys Acta. 2009, 1792 (1): 14-26. 10.1016/j.bbadis.2008.09.017.
Beyens G, Wuyts W, Cleiren E, de Freitas F, Tiegs R, Van Hul W: Identification and molecular characterization of a novel splice-site mutation (G1205C) in the SQSTM1 gene causing Paget’s disease of bone in an extended American family. Calcif Tissue Int. 2006, 79 (5): 281-288. 10.1007/s00223-006-0122-3.
Cavey JR, Ralston SH, Hocking LJ, Sheppard PW, Ciani B, Searle MS, Layfield R: Loss of ubiquitin-binding associated with Paget’s disease of bone p62 (SQSTM1) mutations. J Bone Miner Res. 2005, 20 (4): 619-624. 10.1359/JBMR.041205.
Susani L, Pangrazio A, Sobacchi C, Taranta A, Mortier G, Savarirayan R, Villa A, Orchard P, Vezzoni P, Albertini A, Frattini A, Pagani F: TCIRG1-dependent recessive osteopetrosis: mutation analysis, functional identification of the splicing defects, and in vitro rescue by U1 snRNA. Hum Mutat. 2004, 24 (3): 225-235. 10.1002/humu.20076.
Bertram L, Hiltunen M, Parkinson M, Ingelsson M, Lange C, Ramasamy K, Mullin K, Menon R, Sampson AJ, Hsiao MY, Elliott KJ, Velicelebi G, Moscarillo T, Hyman BT, Wagner SL, Becker KD, Blacker D, Tanzi RE: Family-based association between Alzheimer’s disease and variants in UBQLN1. N Engl J Med. 2005, 352 (9): 884-894. 10.1056/NEJMoa042765.
Moscat J, Diaz-Meco MT, Wooten MW: Signal integration and diversification through the p62 scaffold protein. Trends Biochem Sci. 2007, 32 (2): 95-100. 10.1016/j.tibs.2006.12.002.
Jin W, Chang M, Paul EM, Babu G, Lee AJ, Reiley W, Wright A, Zhang M, You J, Sun SC: Deubiquitinating enzyme CYLD negatively regulates RANK signaling and osteoclastogenesis in mice. J Clin Invest. 2008, 118 (5): 1858-1866. 10.1172/JCI34257.
Kovalenko A, Chable-Bessia C, Cantarella G, Israel A, Wallach D, Courtois G: The tumour suppressor CYLD negatively regulates NF-kappaB signalling by deubiquitination. Nature. 2003, 424 (6950): 801-805. 10.1038/nature01802.
Sundaram K, Shanmugarajan S, Rao DS, Reddy SV: Mutant p62P392L stimulation of osteoclast differentiation in Paget’s disease of bone. Endocrinology. 2011, 152 (11): 4180-4189. 10.1210/en.2011-1225.
Frederick A, Rolfe M, Chiu MI: The human UNP locus at 3p21.31 encodes two tissue-selective, cytoplasmic isoforms with deubiquitinating activity that have reduced expression in small cell lung carcinoma cell lines. Oncogene. 1998, 16 (2): 153-165. 10.1038/sj.onc.1201537.
Fan YH, Yu Y, Mao RF, Tan XJ, Xu GF, Zhang H, Lu XB, Fu SB, Yang J: USP4 targets TAK1 to downregulate TNFalpha-induced NF-kappaB activation. Cell Death Differ. 2011, 18 (10): 1547-1560. 10.1038/cdd.2011.11.
Walsh MC, Kim GK, Maurizio PL, Molnar EE, Choi Y: TRAF6 autoubiquitination-independent activation of the NFkappaB and MAPK pathways in response to IL-1 and RANKL. PLoS One. 2008, 3 (12): e4064-10.1371/journal.pone.0004064.
Xiao N, Li H, Luo J, Wang R, Chen H, Chen J, Wang P: Ubiquitin-specific protease 4 (USP4) targets TRAF2 and TRAF6 for deubiquitination and inhibits TNFalpha-induced cancer cell migration. Biochem J. 2012, 441 (3): 979-986. 10.1042/BJ20111358.
Gutierrez C, Schiff R: HER2: biology, detection, and clinical implications. Arch Pathol Lab Med. 2011, 135 (1): 55-62.
Fransson A, Ruusala A, Aspenstrom P: Atypical Rho GTPases have roles in mitochondrial homeostasis and apoptosis. J Biol Chem. 2003, 278 (8): 6495-6502. 10.1074/jbc.M208609200.
Carcamo C, Pardo E, Oyanadel C, Bravo-Zehnder M, Bull P, Caceres M, Martinez J, Massardo L, Jacobelli S, Gonzalez A, Soza A: Galectin-8 binds specific beta1 integrins and induces polarized spreading highlighted by asymmetric lamellipodia in Jurkat T cells. Exp Cell Res. 2006, 312 (4): 374-386. 10.1016/j.yexcr.2005.10.025.
Levy Y, Ronen D, Bershadsky AD, Zick Y: Sustained induction of ERK, protein kinase B, and p70 S6 kinase regulates cell spreading and formation of F-actin microspikes upon ligation of integrins by galectin-8, a mammalian lectin. J Biol Chem. 2003, 278 (16): 14533-14542. 10.1074/jbc.M207380200.
Nishi N, Shoji H, Seki M, Itoh A, Miyanaka H, Yuube K, Hirashima M, Nakamura T: Galectin-8 modulates neutrophil function via interaction with integrin alphaM. Glycobiology. 2003, 13 (11): 755-763. 10.1093/glycob/cwg102.
Eshkar Sebban L, Ronen D, Levartovsky D, Elkayam O, Caspi D, Aamar S, Amital H, Rubinow A, Golan I, Naor D, Zick Y: The involvement of CD44 and its novel ligand galectin-8 in apoptotic regulation of autoimmune inflammation. J Immunol. 2007, 179 (2): 1225-1235. 10.4049/jimmunol.179.2.1225.
Norambuena A, Metz C, Vicuna L, Silva A, Pardo E, Oyanadel C, Massardo L, Gonzalez A, Soza A: Galectin-8 induces apoptosis in Jurkat T cells by phosphatidic acid-mediated ERK1/2 activation supported by protein kinase A down-regulation. J Biol Chem. 2009, 284 (19): 12670-12679. 10.1074/jbc.M808949200.
Janssens S, Tinel A, Lippens S, Tschopp J: PIDD mediates NF-kappaB activation in response to DNA damage. Cell. 2005, 123 (6): 1079-1092. 10.1016/j.cell.2005.09.036.
Tinel A, Janssens S, Lippens S, Cuenin S, Logette E, Jaccard B, Quadroni M, Tschopp J: Autoproteolysis of PIDD marks the bifurcation between pro-death caspase-2 and pro-survival NF-kappaB pathway. EMBO J. 2007, 26 (1): 197-208. 10.1038/sj.emboj.7601473.
Cuenin S, Tinel A, Janssens S, Tschopp J: p53-induced protein with a death domain (PIDD) isoforms differentially activate nuclear factor-kappaB and caspase-2 in response to genotoxic stress. Oncogene. 2008, 27 (3): 387-396. 10.1038/sj.onc.1210635.
Huang L, Han D, Yang X, Qin B, Ji G, Yu L: PIDD4, a novel PIDD isoform without the LRR domain, can independently induce cell apoptosis in cytoplasm. Biochem Biophys Res Commun. 2011, 407 (1): 86-91. 10.1016/j.bbrc.2011.02.114.
Helfrich MH, Hocking LJ: Genetics and aetiology of Pagetic disorders of bone. Arch Biochem Biophys. 2008, 473 (2): 172-182. 10.1016/j.abb.2008.02.045.
Hocking LJ, Whitehouse C, Helfrich MH: Autophagy: a new player in skeletal maintenance?. J Bone Miner Res. 2012, 27 (7): 1439-1447. 10.1002/jbmr.1668.
Itoh T, Kanno E, Uemura T, Waguri S, Fukuda M: OATL1, a novel autophagosome-resident Rab33B-GAP, regulates autophagosomal maturation. J Cell Biol. 2011, 192 (5): 839-853. 10.1083/jcb.201008107.
We acknowledge the contribution of Danielle Poulin, Research Nurse in the recruitment of PDB patients and healthy donors, and of Edith R. Gagnon, in sequencing the p62P392L mutation. This study was supported by grants from the Canadian Institutes of Health Research (CIHR) (SR) (reference MOP-102707), and SR and LM were supported by the FRQ-S (Fonds de la Recherche du Québec- Santé).
All authors confirm that they have no competing interest.
RK and SR devised the study and analyzed the data. RK conducted the splicing and gene expression studies. MB, GL and SM carried out cell cultures and immunoblots. GL carried out immunofluorescence studies. LM and JPB provided blood samples from each of the healthy donors and PDB patients. SR drafted the paper. All authors reviewed the data, revised the manuscript and approved the final manuscript.
About this article
Cite this article
Klinck, R., Laberge, G., Bisson, M. et al. Alternative splicing in osteoclasts and Paget’s disease of bone. BMC Med Genet 15, 98 (2014). https://doi.org/10.1186/s12881-014-0098-1
- Alternative splicing
- Paget’s disease of bone