Gene expression profiling of early intervertebral disc degeneration reveals a down-regulation of canonical Wnt signaling and caveolin-1 expression: implications for development of regenerative strategies
- Lucas A Smolders1,
- Björn P Meij1,
- David Onis1,
- Frank M Riemers1,
- Niklas Bergknut1, 2,
- Richard Wubbolts3,
- Guy CM Grinwis4,
- Martin Houweling3,
- Marian JA Groot Koerkamp5,
- Dik van Leenen5,
- Frank CP Holstege5,
- Herman AW Hazewinkel1,
- Laura B Creemers6,
- Louis C Penning1 and
- Marianna A Tryfonidou1Email author
© Smolders et al.; licensee BioMed Central Ltd. 2012
Received: 28 August 2012
Accepted: 10 January 2013
Published: 29 January 2013
Early degeneration of the intervertebral disc (IVD) involves a change in cellular differentiation from notochordal cells (NCs) in the nucleus pulposus (NP) to chondrocyte-like cells (CLCs). The purpose of this study was to investigate the gene expression profiles involved in this process using NP tissue from non-chondrodystrophic and chondrodystrophic dogs, a species with naturally occurring IVD degeneration.
Dual channel DNA microarrays were used to compare 1) healthy NP tissue containing only NCs (NC-rich), 2) NP tissue with a mixed population of NCs and CLCs (Mixed), and 3) NP tissue containing solely CLCs (CLC-rich) in both non-chondrodystrophic and chondrodystrophic dogs. Based on previous reports and the findings of the microarray analyses, canonical Wnt signaling was further evaluated using qPCR of relevant Wnt target genes. We hypothesized that caveolin-1, a regulator of Wnt signaling that showed significant changes in gene expression in the microarray analyses, played a significant role in early IVD degeneration. Caveolin-1 expression was investigated in IVD tissue sections and in cultured NCs. To investigate the significance of Caveolin-1 in IVD health and degeneration, the NP of 3-month-old Caveolin-1 knock-out mice was histopathologically evaluated and compared with the NP of wild-type mice of the same age.
Early IVD degeneration involved significant changes in numerous pathways, including Wnt/β-catenin signaling. With regard to Wnt/β-catenin signaling, axin2 gene expression was significantly higher in chondrodystrophic dogs compared with non-chondrodystrophic dogs. IVD degeneration involved significant down-regulation of axin2 gene expression. IVD degeneration involved significant down-regulation in Caveolin-1 gene and protein expression. NCs showed abundant caveolin-1 expression in vivo and in vitro, whereas CLCs did not. The NP of wild-type mice was rich in viable NCs, whereas the NP of Caveolin-1 knock-out mice contained chondroid-like matrix with mainly apoptotic, small, rounded cells.
Early IVD degeneration involves down-regulation of canonical Wnt signaling and Caveolin-1 expression, which appears to be essential to the physiology and preservation of NCs. Therefore, Caveolin-1 may be regarded an exciting target for developing strategies for IVD regeneration.
Degeneration of the intervertebral disc (IVD) is a major cause of low back pain, a major health problem in the Western world . Low back pain resulting from IVD degeneration may be treated medically in combination with physiotherapy. Surgical therapies include decompression with partial discectomy (removal of the diseased IVD tissue), spinal fusion of the affected segment, or partial or total artificial IVD replacement [2, 3]. Although these surgical therapies are generally successful, they are suboptimal since they are not curative and are associated with various complications: decompression/partial discectomy results in spinal instability, spinal fusion can result in adjacent segment degeneration, and IVD replacements/prostheses are associated with failure of the surgical implants [4–9]. Therefore, within the field of regenerative medicine the focus has been on strategies concentrating on biological repair of the degenerating disc using adult stem or progenitor cells, growth factors, and/or gene therapy . However, the biomolecular events involved in IVD degeneration remain largely unexplored [11–13].
Like humans, dogs suffer from spontaneous IVD degeneration that involves similar macroscopic (for example, dehydration of the NP (nucleus pulposus), decrease in disc height), histopathological (for example, chondrocyte proliferation, disorganization of the annulus fibrosus), and biochemical changes (for example, decrease in proteoglycan content, increase in matrix metalloproteinase (MMP) activity) . In humans and dogs, the juvenile, healthy NP of the IVD consists of a dense population of notochordal cells (NCs) embedded in a modest amount of extracellular matrix [15, 16]. The NC has restorative capacity upon other cells, such as chondrocyte-like cells (CLCs) and mesenchymal stem cells, with significant regenerative potential, and thus is an interesting focus for regenerative strategies [13, 17–21]. In humans and dogs, aging and early degeneration of the IVD involves chondroid metaplasia of the NP, which is characterized by the replacement of NCs by CLCs [11, 16, 22]. With regard to this cellular phenomenon, the dog is a unique species, because two subspecies can be distinguished, namely, chondrodystrophic and non-chondrodystrophic dog breeds [22, 23]. Chondrodystrophic breeds are characterized by a disturbed endochondral ossification, resulting in disproportionally short limbs relative to the length of the spine. In these breeds, replacement of NCs by CLCs occurs before 1 year of age, with a concurrent onset of IVD degeneration at all spinal levels . In contrast, non-chondrodystrophic breeds have normal growth of the long bones, and in these dogs the NC remains the predominant cell type of the NP until middle or old age. In non-chondrodystrophic dogs, IVD degeneration generally occurs at older ages compared with chondrodystrophic dogs (> 4 to 5 years of age), and mainly at selected locations (caudal cervical and lumbosacral spine), probably due to a high workload at these spinal levels [22, 24–29].
Therefore, these two dog types reflect a naturally occurring animal model for IVD degeneration, representing differential maintenance of the NC with differential causative factors [14, 30]. Hence, the dog can be considered a unique model for studying the (patho)physiology of the NC and associated early IVD degeneration.
Although the process of early IVD degeneration has been described histopathologically [14, 16, 22], biomolecular signaling pathways involved in the transition from the NP rich in NCs to the NP rich in CLCs (that is, early IVD degeneration) require further investigation. The aim of this study was to investigate the biomolecular signaling profiles associated with NC maintenance and replacement of NCs by CLCs in non-chondrodystrophic and chondrodystrophic dogs, to identify possible targets for IVD regeneration. In the present study, apart from expected biomolecular signaling pathways involved in early IVD degeneration, including Wnt/β-catenin signaling, new pathways were identified. In particular, Caveolin-1, a regulator of Wnt/β-catenin signaling, was found to be a crucial factor in the maintenance of NC health and physiology, and in the initiation of IVD degeneration, being significantly different between the two canine subspecies. These results indicate that Caveolin-1 is an exciting target for further studies.
Materials and methods
All materials used in this study were collected from animals euthanized in other, unrelated experiments approved by the Ethics Committee on Animal Experimentation (DEC) of Utrecht University. Canine IVD tissue was collected from dogs euthanized in studies investigating osteoarthritis (Experiment numbers DEC 2007.III.08.110 and DEC 2009.III.05.037; euthanasia performed by way of an intravenous overdose of pentobarbital) [31–33] and cardiovascular disease (Experiment number DEC 2007.II.01.029; euthanasia performed under general anesthesia by way of fibrillation and subsequent excision of the heart) [, other, unpublished data].
Murine IVD tissue was collected from mice euthanized for studies investigating the role of Caveolin-1 in liver regeneration (Experiment number DEC 2008.III.01.001; euthanasia performed by way of exsanguination under isoflurane anesthesia; work not yet published). In these unrelated experiments, the animals were sacrificed for the collection of tissue other than IVD. All experimental procedures were performed strictly according the guidelines set by the Ethics Committee of Utrecht University. The Ethics Committee of Utrecht University approved post-mortem harvesting of the IVD tissue employed in the present study.
Cervical (vertebrae C2 to T1) and thoracolumbar (vertebrae T10 to S1) spines were collected from five mongrel dogs (non-chondrodystrophic group, age range 13 to 60 months, body weight range 26.6 to 32.1 kg) and six Beagle dogs (chondrodystrophic group, age range 25 to 36 months, weight range, 13.6 to 16.0 kg). Spines were resected en bloc and cut into single functional spinal units (FSUs; endplate-IVD-endplate). Each FSU was cut in the sagittal plane into two equal parts. The NP was collected from one part, snap frozen in liquid nitrogen, and stored at -70°C. The other part was fixed in 4% neutral buffered formaldehyde, decalcified in 10% ethylenediaminetetraacetic acid (EDTA) for 3 months at room temperature, and embedded in paraffin.
Histopathological grading of IVD samples
To further assess the cellular phenomenon involving the transition from NC-rich to CLC-rich NP, and thus the applied histopathological grading, the gene expression of the notochordal markers Brachyury [12, 37] and Cytokeratin 8 [38, 39] was investigated in all groups (qPCR analysis, see below).
DNA microarray analysis
Four NPs from the three histopathological groups for each breed type were randomly selected and used to isolate total RNA, using the RNeasy Fibrous Tissue Mini Kit (Qiagen, Venlo, the Netherlands) according to the manufacturer's instructions. After on-column DNase-I treatment (Qiagen RNAse-free DNase kit, Venlo, the Netherlands), RNA was quantified spectrophotometrically using Nanodrop ND-1000 (Isogen Life Science, De Meern, the Netherlands) and RNA integrity was determined using a Bioanalyzer 2100 (Agilent Technologies, Amstelveen, the Netherlands). RNA integrity numbers of the samples varied from 4.5 to 7.4 (mean ± SD 6.2 ± 0.7), indicating that these samples could be used to perform valid gene expression analysis .
A two-color DNA microarray with a reference experiment design was performed on 44 k Canine Gene Expression Microarrays V1 (G2519F, Agilent Technologies, Amstelveen, the Netherlands) . Double round RNA amplification and labeling were performed as described before  on an automated system (Caliper Life Sciences, 's-Hertogenbosch, the Netherlands) with 10 to 50 ng total RNA input from each sample. The common RNA reference pool consisted of a multitude of canine organs, including liver, spleen, kidney, lung, heart, intestine, and bone. Hybridizations were performed on an HS4800PRO system with QuadChambers (Tecan Benelux B.V.B.A., Giessen, the Netherlands), using 300 to 500 ng labeled cRNA per channel as described previously . From each group two samples were labeled with Cy5 and hybridized against the common reference cRNA (Cy3) on dual channel arrays, and two samples were hybridized in dye swap. Hybridized slides were scanned on an Agilent scanner (G2565BA) at 100% laser power, 100% photomultiplier tube (PMT). After automated data extraction using Imagene 8.0 (BioDiscovery, El Segundo, CA, United States of America), print-tip loess normalization was performed  on mean spot intensities. Dye bias was corrected based on a within-set estimate . Patterns of gene expression were compared between the three stages of IVD degeneration (NC-rich, mixed, CLC-rich NP) within each breed type and between the two breed types. Data were analyzed using microarray analysis of variance (MAANOVA) . In a fixed-effect analysis, sample, array, and dye effects were modeled. P-values were determined with a permutation F2-test, in which residuals were shuffled 5,000 times globally. Genes with a P-value < 0.05 after Benjamini-Hochberg determination of false discovery rate (FDR) were considered significantly changed; a change cutoff value of 1.3-fold was applied.
Differentially expressed genes were converted to their human homologues, and the following array comparisons were included in functional pathway analysis using the GeneGo MetaCore platform : 1) NC-rich group vs. mixed group, mixed group vs. CLC-rich group, and NC-rich group vs. CLC-rich group in both chondrodystrophic and non-chondrodystrophic dogs; 2) non-chondrodystrophic vs. chondrodystrophic breeds: NC-rich group, mixed group, and CLC-rich group. Pathways showing significant changes in gene expression were selected and analyzed further as described below.
Quantitative PCR (qPCR)
Six NPs from the three histopathological groups for each breed type were analyzed (cDNA samples used for microarray and additional, similarly processed samples). Unlabeled microarray cDNA from the first round of amplification was used. For the additional samples, cRNA was synthesized by extracting RNA as described for the microarray, followed by one round of amplification using in vitro transcription , assuring similar pre-PCR treatment of both cRNA sources.
For all samples, cDNA was synthesized using the iScript™ cDNA Synthesis Kit (Biorad, Veenendaal, the Netherlands). qPCR was performed using a MyIQ thermal cycler, IQ SYBRGreen SuperMix (BioRad, Veenendaal, the Netherlands) and dog-specific primers (Eurogentec, Maastricht, the Netherlands) (see Additional file 1). Each sample was treated as an individual sample, and samples from different histopathological grades and subspecies were analyzed separately (in tehnical duplicates). Primers were designed for the notochordal marker genes Brachyury and Cytokeratin 8, and based on the microarray analyses primers were designed for the following Wnt target genes: Wnt3a, Wnt7b, Wnt inhibitory factor 1 (Wif1), Frizzled 1 (Fzd1), Low density lipoprotein receptor-related protein 5 (Lrp5), Dickkopf homolog 3 (Dkk3), Integrin-linked kinase (Ilk), Caveolin-1 (Cav1) and Axin2 [49–57]. Conditions for the qPCR experiments were carefully validated as described previously . The amplification efficiency was between 90% and 110%. Relative expression was calculated by the efficiency-corrected delta-delta Ct (ΔΔCt) method  using a set of five reference genes (see Additional file 1).
Using the ΔΔCt as a parameter value, linear mixed models [59, 60] were designed to statistically analyze the obtained data for each individual target gene. P-values were calculated to analyze differences in relative gene expression between groups (chondrodystrophic and non-chondrodystrophic dogs) and degeneration stages (NC-rich, mixed, CLC-rich). The Benjamini-Hochberg correction was used to correct for multiple comparisons . P < 0.05 was considered statistically significant. For a complete description of the statistical analyses, see Additional file 2.
β-Catenin protein expression
To validate differences in Wnt/β-catenin signaling between non-chondrodystrophic and chondrodystrophic dogs, immunohistochemistry of β-catenin was performed as described before , and analyzed focussing on differences between non-chondrodystrophic and chondrodystrophic dogs. Integrated density for β-catenin staining of the NP was corrected for the surface area of the NP, positively stained for β-catenin, in order to correct for the difference in disc size between the two breed groups. Furthermore, total protein was extracted from lumbar CLC-rich NPs from non-chondrodystrophic and chondrodystrophic dogs using radioimmunoprecipitation assay (RIPA) buffer; protein concentration was determined with a Lowry assay. Aliquots of protein were subjected to 10% SDS-PAGE (15 μg/lane). Protein (5 μg) from a human insulinoma CM cell line  served as positive control. The proteins were electroblotted onto a Hybond-C nitrocellulose membrane (Amersham Biosciences, RPN203C, Bath, United Kingdom). After blocking with 4% non-fat dry milk in Tris-buffered saline containing 0.1% Tween20 (TBST0.1), the membrane was incubated overnight at 4°C with β-catenin antibody (Ab6302, Abcam, Cambridge, United Kingdom, 1:500 in 4% BSA in TBST0.1). After three 5-minute washes in TBST0.1, the membrane was incubated for 90 minutes with anti-rabbit horseradish peroxidase (HRP)-conjugated secondary antibody (R&D, Minneapolis, MN, United States of America; HAF008, 1:20.000 in TBST0.1). Protein expression was detected using an enhanced chemiluminescence substrate (ECL Advance, Amersham RPN2135, Bath, United Kingdom) in a ChemiDoc XRS System (Bio-Rad Laboratories, Veenendaal, the Netherlands). Control experiments were included by omitting the primary antibodies. After completing the western blot for β-catenin, the membranes were washed in TBST0.1 and re-used to determine α-tubulin protein expression. For this purpose, all steps as described above were performed except for the antibody incubations: the membrane was incubated with primary antibody for α-tubulin (Sigma T6199, 1:750, Zwijndrecht, the Netherlands) for 2 hours at room temperature, followed by secondary anti-mouse HRP-conjugated antibody (R&D HAF007,1:20000, Minneapolis, MN, United States of America). All experiments were performed in triplicate.
Paraffin-embedded IVD (as described above) sections from five non-chondrodystrophic and five chondrodystrophic dogs for each group (NC-rich, mixed, CLC-rich; n = 10 per group) were subjected to antigen retrieval in 10 mM citrate buffer (pH 6.0), followed by blocking of endogenous peroxidase activity . Nonspecific background staining was minimized by pre-incubation with blocking buffer containing 10% normal goat serum (Sigma-Aldrich, Zwijndrecht, the Netherlands)/0.1% Tween-20 (Tween-20, Boom BV, Meppel, the Netherlands) in PBS for 30 minutes, and an overnight incubation at 4°C with the primary antibody monoclonal mouse anti-Caveolin-1 (Transduction Laboratories, Breda, the Netherlands mAb2297, 2.5 μg/ml, 1:100 in PBS with 0.1% Tween-20). After sections were washed in PBS buffer/0.025% Triton X, Caveolin-1 was visualized with the goat anti-mouse Envision System and the liquid diaminobenzidine (DAB) chromogen system (Dako, Heverlee, Belgium) and counterstained with hematoxylin (Hematoxylin QS, Vector Laboratories Inc., Peterborough, United Kingdom). In negative control sections, the primary antibody was omitted or replaced with its respective serum. All sections were stained in the same session.
Detailed overview images of each stained slide were made using a Colorview IIIU digital camera (Olympus, Zeist, the Netherlands) mounted to a BX-40 microscope (Olympus, Zeist, the Netherlands). The total NP surface of each sample was measured by defining the perimeter of the NP, excluding the transition zone. A custom-made color range selection optimized for Caveolin-1 specific staining was used to calculate the proportion of the NP surface area that stained for Caveolin-1, and the mean gray value (staining intensity) for Caveolin-1 staining in each sample, as described previously .
For each parameter, linear mixed models were used to calculate P-values to analyze differences between groups (chondrodystrophic and non-chondrodystrophic) and degeneration stages (NC-rich, mixed, CLC-rich) (for a complete description of the statistical analyses: Additional file 2). The Benjamini-Hochberg correction was used to correct for multiple comparisons . P < 0.05 was considered statistically significant.
To investigate the relationship between Caveolin-1 and canonical Wnt signaling in the NP, their colocalization was investigated by way of simultaneous immunofluorescence analysis of both proteins. Paraffin-embedded IVD slides were used for immunofluorescent labeling of Caveolin-1 as described above, except that the slides were incubated overnight with a combination of primary antibodies monoclonal mouse anti-Caveolin-1 (Transduction Laboratories, Breda, the Netherlands, mAb2297, 2.5 μg/ml, 1:100 in PBS with 0.1% Tween-20) and polyclonal rabbit anti-β-catenin (Abcam, Cambridge, United Kingdom, ab6302, 1:50 in PBS with 0.1% Tween-20). The secondary antibodies used were 1:100 goat anti-mouse antibody conjugated to Alexa488 (2.5 μg/ml; Invitrogen, Breda, the Netherlands) and 1:100 goat anti-rabbit antibody conjugated to Alexa568 (2.5 μg/ml; Invitrogen, Breda, the Netherlands). Topro-3 iodide (2 μg/ml; Invitrogen, Breda, the Netherlands, T3605) was used to stain the nucleus.
To outline the proximity of the fluorescently-marked Caveolin-1 and β-catenin proteins, profile intensity plots were generated in LAS-AF imaging software (Leica microsystems, Wetzlar, Germany). Straight lines were drawn across representative cell bodies and intensity profiles were extracted from the channels visualizing and measuring the fluorescence of Caveolin-1, β-catenin, and Topro-3. The data were subsequently exported to Microsoft Excel (Microsoft Corporation, Amsterdam, the Netherlands) and plotted.
Caveolin-1 expression in NCs in vitro
On the basis of the previous analyses, Caveolin-1 was investigated in NCs in vitro. NCs were isolated from the NPs of cervical (C2 to T1) and lumbar (L1 to S1) IVDs from six, young-adult, mongrel dogs (non-chondrodystrophic, aged 16 to 18 months, and weighing 16 to 24 kg). The NCs were cultured in their original cluster-like conformation as described previously  on coverslips in 6-well plates (Falcon Multiwell Primaria, Becton Dickinson, Breda, the Netherlands) in penicillin/streptomycin (P/S)-FCS, 10%)-DMEM-F12 under normoxic conditions (5% CO2) at 37°C for 10 days.
For RNA isolation, non-adherent cells at days 0 and 2 (NCs first adhered on day 4) were collected by centrifuging the medium at 1500 RPM at 4°C for 1 minute; on days 4, 6, 8, and 10, the medium was removed, the wells were washed with RNase-free Hank's solution, and the adherent cells were lysed and used for analysis. For all time points, total cellular RNA was isolated using the RNeasy Mini Kit (Qiagen, Venlo, the Netherlands) according to the manufacturer's instructions. The relative gene expression of Caveolin-1 was analyzed as described above.
Non-adhered cells on days 0 and 2 were collected from the culture medium and mounted on positively charged slides (Klinipath, Duiven, the Netherlands) using a Shandon Cytospin 4 (Therma Scientific, Landsmeer, the Netherlands). Cells were fixed on days 0, 2, 4, 6, 8, and 10, and were used for immunofluorescent labeling of Caveolin-1 as described above, except that the secondary antibody used was 1:100 donkey anti-mouse antibody conjugated to Alexa488, (2.5 μg/ml; Invitrogen, Breda, the Netherlands). Topro-3 iodide (2 μg/ml; Invitrogen, Breda, the Netherlands, T3605) was used to stain the nucleus. The cells were mounted with Fluorsave (Calbiochem, Darmstadt, Germany). Images at 5 random locations in each sample were acquired by a sequential recording procedure on a multiphoton imaging station (MP2100, Zeiss, Herfordshire, United Kingdom). Immunofluorescent images of the cells were analyzed with CellProfiler 2.0 software package (Massachusetts Institute of Technology, Massachusetts, United States of America), with cell nuclei detected on the Topro-3 images (nuclear staining) and the Caveolin-1 expression signal on the Caveolin-1 images. For each image, the total number of cells and the total Caveolin-1 signal were measured, and the estimated mean intensity of Caveolin-1 protein staining per cell was calculated by dividing the total Caveolin-1 signal, after subtraction of background, by the total number of cells.
For both ΔCT for Caveolin-1 gene expression, and Caveolin-1 protein expression per cell, linear mixed models were designed to calculate P-values to analyze differences in the Caveolin-1 gene and protein expression between time points in culture. The Benjamini-Hochberg correction was used to correct for multiple comparisons . P < 0.05 was considered statistically significant (see Additional file 2 for a complete description of the statistical analyses).
NP in Caveolin-1 knock-out (KO) mice
To further assess the role of Caveolin-1 in NC preservation and early IVD degeneration, spines were collected from 3-month-old Caveolin-1 KO-mice (Cav tm1Mls , JAX®, the Jackson Laboratory, Maine, United States of America) and wild-type mice (strain B6129SF2, JAX®). Spines were fixed in 4% neutral buffered formalin, decalcified (7 days in EDTA at 4°C), and embedded in paraffin. Mid-sagittal sections (4 μm) were mounted on Microscope KP+ slides (Klinipath B.V., Duiven, the Netherlands) and stained with H&E and alcian blue/picosirius red, the latter highlighting proteoglycan (blue) and collagen content (red) . Multiple sections of the NP of multiple cervical and lumbar IVDs (n = 4) from each mouse were histopathologically evaluated.
Microarray: changes in gene expression
In NPs from non-chondrodystrophic dogs, the NC-rich, mixed, and CLC-rich groups consisted of 100.0% NCs, 45.4% to 87.0% NCs, and 100.0% CLCs, respectively; in NPs from chondrodystrophic dogs, these proportions were 93.9% to 100.0% NCs, 39.7% to 89.9% NCs, and 100.0% CLCs, respectively. Therefore, the applied histopathological classification resulted in the NC-rich, mixed, and CLC-rich groups being comparable between chondrodystrophic and non-chondrodystrophic dogs, allowing reliable investigation of early IVD degeneration within both types of breed, and comparison between the two types of breed.
Top 25 up- and downregulated genes for the microarray comparison of notochordal cell (NC)-rich nucleus pulposus (NP) (reference) vs. chondrocyte-like cell-(CLC)-rich NP in chondrodystrophic dogs
Chondrodystrophic dogs: NC-rich NP vs. CLC-rich NP
Total upregulated genes: 1.178
Total downregulated genes: 1,741
GO term: biological Process
Protein modification process
Cellular iron ion homeostasis
Negative regulation of canonical Wnt receptor signaling pathway
Peptide cross-linking via chondroitin 4-sulfate glycosaminoglycan
Cartilage oligomeric matrix protein
Negative regulation of bone mineralization
Collagen fibril organization
Cellular response to erythropoietin
RAN binding protein 3-like
Retinol binding protein 4, plasma
Protein complex assembly
Cysteine dioxygenase, type I
Response to glucagon stimulus
Adenylate cyclase 2
Activation of adenylate cyclase activity by G-protein signaling pathway
Microfibrillar associated protein 5
S100 calcium binding protein A12
Phosphotyrosine interaction domain containing 1
Cell wall macromolecule catabolic process
SPARC-like 1 (hevin)
Serpin peptidase inhibitor, clade G (C1 inhibitor), member 1
Regulation of proteolysis
Sphingomyelin phosphodiesterase, acid-like 3A
Sphingomyelin catabolic process
Golgi to plasma membrane CFTR protein transport
tRNA-yW synthesizing protein 3 homolog
A kinase (PRKA) anchor protein 12
Phospholipase C-like 1
Lipid metabolic process
Myosin, light chain 9, regulatory
Regulation of muscle contraction
Mitochondrial ribosomal protein S27
Ectonucleotide pyrophosphatase/phosphodiesterase 2
Regulation of cell migration
Cell differentiation involved in embryonic placenta development
Carbohydrate metabolic process
Interspecies interaction between organisms
Nucleosome assembly protein 1-like 1
RAB20, member RAS oncogene family
Small GTPase mediated signal transduction
Positive regulation of protein binding
Potassium voltage-gated channel, delayed-rectifier, subfamily S, member 3
Positive regulation of phosphorylation
Sorbin and SH3 domain containing 2
Phosphatidylcholine transfer protein
Cholesterol metabolic process
Kv channel interacting protein 1
Carbonic anhydrase II
Thy-1 cell surface antigen
RAB38, member RAS oncogene family
GTP catabolic process
Sema domain, immunoglobulin domain (Ig), short basic domain, secreted, (semaphorin) 3C
Neural tube development
Metacore pathway map analysis could not be performed on the gene regulation results from non-chondrodystrophic dogs, because relatively too few genes were down- or upregulated in this breed group (see Additional file 4).
Wnt/β-catenin signaling was analyzed further because it is involved in both the regeneration and the degeneration of various tissues . The expression of the Wnt7b (Wnt ligand), Wif1 (inhibits by binding to Wnt ligands), Ilk (inhibits glycogen synthase kinase 3-β), and Lrp5 (Wnt co-receptor) genes was significantly changed and these Wnt/β-catenin target genes were analyzed further by qPCR, as were additional targets involved in canonical Wnt signaling: Wnt3a (Wnt ligand), Fzd1 (Wnt receptor), Dkk3 (negative regulator of Wnt), and Axin2 (Wnt read-out) (see Additional file 3).
Quantitative PCR of the canonical Wnt signaling pathway and β-catenin protein expression
The relative gene expression of axin2, which is a highly reliable read-out for the activity of Wnt/β-catenin signaling [69–71], was significantly lower in the CLC-rich group than in the NC-rich and mixed groups in both non-chondrodystrophic and chondrodystrophic dogs (Figure 2). In chondrodystrophic dogs, this decrease in Axin2 gene expression may be related to the gene expression of Wif1 (inhibits Wnt ligands), which was significantly upregulated in the CLC-rich group compared with the NC-rich and mixed groups. In non-chondrodystrophic dogs, no significant changes in Wif1 gene expression were found. However, gene expression of the Wnt ligand Wnt7b, which activates canonical Wnt signaling through interactions with Fzd and LRP5 , was significantly higher in the CLC-rich group compared with the mixed group in both non-chondrodystrophic and chondrodystrophic dogs.
Compared with non-chondrodystrophic dogs, Axin2 gene expression was significantly higher in chondrodystrophic dogs in all histopathological groups, indicating an overall higher Wnt signaling activity in chondrodystrophic dogs. Accordingly, the integrated density of β-catenin corrected for the NP surface area positively stained, was significantly higher in the CLC-rich NP of chondrodystrophic dogs compared with non-chondrodystrophic dogs (see Figure S1A and B in Additional file 6, and Additional file 7). Due to the abudance of matrix protein in the native CLC-rich NP tissue, western blot analysis of active β-catenin expression was cumbersome and quantification of the data was not reliable. However, the chondrodystrophic CLC-rich NP appeared to have less degraded β-catenin compared to non-chondroystrophic dogs (see Figure S1C in Additional file 6, and Additional file 7). These findings are sustained by investigation of the gene expression of the Wnt receptor Fzd1, which was significantly higher in chondrodystrophic dogs than in non-chondrodystrophic dogs for all three histopathological stages. Gene expression of Wnt7b was significantly higher in non-chondrodystrophic dogs compared with chondrodystrophic dogs in all histopathological stages.
The relative gene expression of Lrp5, Ilk, and Dkk3 remained unchanged in both non-chondrodystrophic and chondrodystrophic dogs. Gene expression of the Wnt3a was undetectable in all groups in both breed types.
The microarray analyses showed significant changes in Caveolin-1, -2 and -3. Early IVD degeneration involved significant downregulation of Caveolin-1 and -2, and significant upregulation of Caveolin-3 (see Additional files 3 and 4). Given the role of Caveolin-1 in the regulation of canonical Wnt signaling [55, 56, 73] and the reported upregulation of Caveolin-1 in degenerated human IVDs , its gene and protein expression were further investigated by way of qPCR and immunohistochemistry.
Caveolin-1 protein was predominantly located in the cell membranes of NCs, and occasionally in their cytoplasm; Caveolin-1 protein was seldom observed in CLCs. The proportion of the NP surface area that stained for Caveolin-1 was significantly lower in the CLC-rich NP than in the NC-rich and mixed NP in both non-chondrodystrophic and chondrodystrophic dogs, indicating decreased Caveolin-1 protein expression (see Additional file 8). No significant differences were found in the mean gray value (staining intensity) between the NC-rich, mixed, and CLC-rich groups in both non-chondrodystrophic and chondrodystrophic dogs. The intensity of Caveolin-1 staining in all three histopathological groups was significantly higher in non-chondrodystrophic than in chondrodystrophic dogs.
To further assess the relationship between Caveolin-1 and canonical Wnt signaling in the NP, co-immunofluorescence of Caveolin-1 and β-catenin was performed. Profile intensity plots (Figure 3) showed clear signal peaks of β-catenin expression within the cell nucleus of NCs, indicative of active canonical Wnt signaling. Also, clear protein expression of Caveolin-1 and β-catenin was colocalized at the cell membrane of NCs, which is supportive of interaction of these proteins within the NC-rich NP. In CLCs, Caveolin-1 and β-catenin protein expression was rarely observed, which is indicative of less active canonical Wnt signaling and no colocalization and potential protein interaction.
To verify the possible role of Caveolin-1 in NCs, its expression and intracellular distribution in primary NCs was investigated
The physiological role Caveolin-1 in the preservation of NCs was further investigated in Caveolin-1 KO mice
Gene expression profiling of chondroid metaplasia of the NP, which is seen in physiological aging as well in early (subclinical) degeneration of the IVD, revealed novel genes and pathways involved in the preservation of the NC in the healthy IVD. The results obtained indicate that IVD degeneration in dogs involved significant downregulation of Wnt/β-catenin signaling and Caveolin-1; the latter seems to be crucial for NC maintenance and an important factor in IVD degeneration.
Material was obtained from non-chondrodystrophic and chondrodystrophic dogs, a unique animal model in view of investigating the NC and degeneration/regeneration of the IVD: in non-chondrodystrophic dogs, the NC remains the predominant cell type of the NP during the majority of life, with IVD degeneration only occurring at old age and in selected IVDs; conversely, in chondrodystrophic dogs the NC is lost early in life, with concurrent degeneration of all IVDs. Therefore, these two dog types reflect a naturally occurring animal model representing differential maintenance of the NC and associated differences in maintenance of optimal matrix health of the IVD [14, 30].
In the present study the main focus was the transition of a healthy NC-rich NP into an aging/degenerating CLC-rich NP. However, recent studies have indicated that both the healthy and degenerating NP, apart from NCs and CLCs, contain stem/progenitor cells (approximately 1% of the NP cell population) [75–77]. These NP progenitor cells express so-called stemness genes, and have been shown to functionally differentiate into the adipogenic, chondrogenic, and neurogenic lineage . As these progenitor cells are present within the NP, it is very well possible that these cells had a profound effect upon the signaling of NCs and CLCs and on the studied biomolecular signaling events, particularly with respect to the different phenotypes of the non-chondrodystrophic and chondrodystrophic NP investigated. Additional studies are warranted on this specific matter.
The relative expression of the notochordal markers Brachyury and Cytokeratin 8 [12, 37–39] revealed marked differences between chondrodystrophic and non-chondrodystrophic dogs. In the initial stages of early IVD degeneration (classified as the mixed group), Brachyury and Cytokeratin 8 gene expression remained constant, indicating the preservation of the NC phenotype in the NP at this stage in both chondrodystrophic and non-chondrodystrophic dogs. However, Brachyury and Cytokeratin 8 gene expression was significantly decreased in the CLC-rich NP from chondrodystrophic, but not non-chondrodystrophic dogs. Decreased expression of NC markers is also seen in degenerated human IVDs  and is indicative of a significant loss of NCs from the NP [12, 37–39]. The maintenance of high levels of Brachyury and Cytokeratin 8 gene expression in the non-chondrodystrophic, CLC-rich NP indicates that the NP cells in non-chondrodystrophic dogs undergo significant morphological changes, but can preserve characteristics of NCs on the basis of their gene expression pattern. These findings are consistent with the observation that there are niches of NCs in degenerated, human CLC-rich NP [11, 12, 38, 78]. Therefore, based on histopathological assessment and the expression of notochordal marker genes, and in accordance with previous studies [17, 19, 79, 80], loss of the notochordal phenotype would appear to be associated with the accelerated IVD degeneration seen in chondrodystrophic dogs, whereas in non-chondrodystrophic dogs the NCs or CLCs retain notochordal characteristics, accounting for the relatively low prevalence of IVD degeneration observed in this breed type [22, 81].
Canonical Wnt signaling activity is decreased in early IVD degeneration, with a clear differences between non-chondrodystrophic and chondrodystrophic dogs
An increase in canonical Wnt signaling activity results in increased gene expression of Axin2, which is considered a reliable read-out for the activity of canonical Wnt signaling [69–71]. Axin2 gene expression was significantly downregulated in the CLC-rich NP compared with the NC-rich and mixed NP in both non-chondrodystrophic and chondrodystrophic dogs, indicating that early (subclinical) IVD degeneration involves a significant reduction in canonical Wnt signaling activity. These results seem contradictory to our previous findings , showing that the CLC-rich NP had a higher Axin2 gene expression than the NC-rich NP. However, in that particular study CLC-NP tissue from chondrodystrophic dogs was compared with NC-rich NP tissue from non-chondrodystrophic dogs. Indeed, chondrodystrophic NPs exhibit significantly higher levels of Axin2 gene expression than non-chondrodystrophic dogs in all histopathological groups investigated in this study (Figure 2), which explains the apparently contradictory results between the previous report  and the present one. The significantly higher Axin2 gene expression found in chondrodystrophic dogs was further sustained by the significantly higher signal intensity for β-catenin-protein expression within the CLC-rich NP tissue as shown by immunohistochemistry, and most probably less degraded β-catenin as shown by the western blot analysis (see Additional files 6 and 7).
The gene expression of Wnt7b, which activates canonical Wnt signaling through interacting with the Wnt receptors LRP5 and Fzd1 , was increased during NP chondrification and was consistently higher in all histopathological groups in non-chondrodystrophic dogs than in chondrodystrophic dogs. Increased Wnt7b gene expression has also been reported in the cartilage of patients with osteoarthritis (OA) and rheumatoid arthritis (RA) . In that respect, the significantly higher gene expression of Wnt7b may reflect a response to increase canonical Wnt signaling activity and an attempt to preserve NP health. Interestingly, canonical Wnt signaling activity (Axin2 and Fzd1 expression) was markedly higher in chondrodystrophic dogs in all histopathological stages, as was also reported earlier in canine NP cells . In chondrodystrophic dogs, the decrease in Wnt signaling activity during NP chondrification may be explained by the increased expression of Wif1 and Frizzled-related protein (Frzb) (Wnt inhibitors, Figure 2 see alos Additional file 3), and the reduced expression of r-spondin-3 (rspo3; Wnt activator, see Additional file 3). It is tempting to hypothesize about the potential role of Wnt signaling in the transition from the healthy, NC-rich NP to the aged or degenerated CLC-rich NP. As Wnt signaling regulates notochord fate and stem cell renewal and apoptosis [83, 84], decreased Wnt signaling may result in increased apoptosis and decreased self-renewal of NCs or NP stem cells, and ultimately in chondroid metaplasia of the NP. The higher canonical Wnt signaling seen in chondrodystrophic dogs may reflect an ineffective attempt to preserve the notochordal phenotype of NP cells or to regulate stem cell renewal and apoptosis. Conversely, given the involvement of Wnt signaling in tissue degeneration [85, 86], higher Wnt signaling in chondrodystrophic dogs might reflect a diminished capacity to limit Wnt signaling, resulting in accelerated extracellular matrix breakdown, as is observed in patients with OA or RA [82, 85]. This line of thought may also be in accordance with recent findings showing increased β-catenin protein expression in degenerated human IVD tissue as compared with healthy controls . Altogether, these data illustrate the dual role of Wnt signaling, which requires further elucidation with respect to the transition from the NC-rich to the CLC-rich NP, also with respect to the role played by progenitor cells within the NP in various stages of degeneration [75–77].
Loss of NC phenotype involves significant downregulation of Caveolin-1 expression
Caveolin-1 is required for notochord development , and caveolins regulate canonical Wnt signaling by recruiting β-catenin to caveolae membranes, thereby inhibiting Wnt/β-catenin signaling and reinforcing cell-cell adhesion mechanisms , and by internalizing LRP6 (Wnt receptor), thereby activating canonical Wnt signaling . Furthermore, Caveolin-1 stimulates canonical Wnt signaling by activating integrin-linked kinase which inhibits glycogen synthase kinase 3-β, a key enzyme in Wnt/β-catenin signaling that phosphorylates β-catenin leading to the subsequent degradation of this molecule . Also, Caveolin-1 stimulates canonical Wnt signaling through the accumulation of β-catenin to caveolae membranes, thereby preventing degradation of β-catenin by glycogen synthase kinase 3-β . The performed co-immunofluorescence analysis of Caveolin-1 and β-catenin (the effector protein of canonical Wnt signaling) showed protein expression peaks of both proteins localized at the cell membrane, which may indicate an interaction between these proteins within the NC. However, in CLCs no such protein expression/colocalization was observed. It is tempting to hypothesize that within the NC, caveolae function to regulate β-catenin signaling by preventing its degradation and to reinforce cell-cell clusters, which is a morphological characteristic of NCs . In line with the decrease in canonical Wnt signaling in early IVD degeneration, Caveolin-1 gene expression was significantly downregulated in both breed types. Caveolin-1 protein expression was observed almost exclusively in NCs and decreased significantly with chondrification/degeneration of the NP (Figure 3). However, other authors have reported increased Caveolin-1 expression to be associated with senescence of NP cells [74, 89] and chondrocytes, with subsequent IVD degeneration and OA , respectively. This discrepancy is probably because the other studies investigated advanced stages of degeneration involving CLCs only, whereas we investigated early degeneration, involving the transition from the NC-rich to CLC-rich NP. Apart from cell senescence, Caveolin is known to play a significant role in various cellular processes, including stem cell regulation and proliferation . Therefore, Caveolin-1 may act differently according to the triggering signals and cellular context, such as the stage of IVD degeneration and in different cell types . The relationship between Wnt/β-catenin signaling and Caveolin-1 in the NP in different stages of degeneration and in NCs compared with CLCs requires further investigation. Interestingly, the expression of Caveolin-1 protein was consistently lower in chondrodystrophic NP than in non-chondrodystrophic NP regardless of histopathological stage, indicating a direct relationship between the absence of Caveolin-1 and the accelerated loss of NCs from the NP in chondrodystrophic dogs.
Caveolin-1 appears to fulfill essential functions in the NC cytoskeleton
In an attempt to understand the role of caveolin-1 in NC physiology, we studied its distribution in NCs in culture. Caveolin-1 protein was detected in the cell membrane and in intracellular membranes (Figure 4), as reported previously . As expected, Caveolin-1 gene and protein expression increased significantly when the cells adhered to the culture plate (days 4 and 6). Caveolin interacts with actin filaments of the cytoskeleton  and its increased expression may be involved in the formation of the NC cytoskeleton and adherence of NCs to the culture plate. In accordance with these observations, microarray analysis revealed significant changes in cytoskeletal components, supporting a role for Caveolin-1 in the NC cytoskeleton. NCs are known to contain a dense network of intracellular Actin, which may be involved in the homeostasis of the intracellular vesicles and cluster formation of these cells . An important function of Caveolin-1 may be to regulate the interaction of caveolae with the actin cytoskeleton, thereby controling whether caveolae are at the cell surface or traveling to interior sites, or regulating the homeostasis of intracellular vesicles and intercellular clusters .
Absence of Caveolin-1 coincides with decreased NC preservation and early IVD degeneration
The essential role of Caveolin-1 in NC physiology was corroborated by the IVD phenotype of the Caveolin-1 KO mice. Unlike NP from wild-type mice, NP from Caveolin-1 KO mice showed relatively few healthy NC clusters; most NP cells lacked the morphological characteristics of NCs and showed signs of apoptosis, and the NP contained an abundance of intercellular chondroid matrix, similar to the CLC-rich NP (Figure 5). These changes are similar to the histopathological changes observed in the transition from an NC-rich to CLC-rich NP (Figure 1) [14, 22]. Therefore, these findings suggest that Caveolin-1 is essential for NC maintenance, and that decreased Caveolin-1 expression is an important factor in NC physiology and IVD degeneration. To further investigate the role of Caveolin-1 and its relationship with Wnt/β-catenin, future studies need to focus on Wnt/β-catenin signaling in Caveolin-1 KO mice.
Apart from the involvement in IVD degeneration, Caveolin proteins are involved in stem cell regulation and proliferation, as well as in the pathogenesis of cancers, pulmonary hypertension, cardiomyopathy, diabetes, and muscular dystrophy . There has been an increasing interest in the application of Caveolin-mimetic peptides for the treatment of both cancer and pulmonary hypertension . The findings of this study suggest that Caveolin-1 is crucial for NC maintenance and IVD health, and this protein may be regarded an exciting target for developing ways to regenerate the IVD. For example, the degenerated IVD may be treated by locally applying Caveolin-1-mimetic peptide or by promoting the expression of Caveolin-1 in the cells of the degenerated IVD, thereby promoting regeneration of the degenerated tissue. Also, it should be taken into consideration that Caveolin-1 may exert different actions dependent on the cell context and the stage of degeneration , and further in vitro mechanistic studies are required to test this concept.
In the present study the main focus was to investigate NCs and CLCs. However, the process of early IVD degeneration may also involve significant changes in NP progenitor/stem cells, and these cells may have a significant influence on the biomolecular signaling events within the NP.
A potential limitation of this study is that relatively few genes were differentially expressed when comparing the NC-rich, mixed, and CLC-rich NPs in the non-chondrodystrophic dogs, whereas in chondrodystrophic dogs relatively many genes were up- or downregulated. This might be because of the genetic heterogeneity of the non-chondrodystrophic sample (mongrels of different age and variable size), whereas the chondrodystrophic sample consisted of Beagle dogs of the same age and standardized size. For the chondrodystrophic breed group (Beagles), an appropriate sample size for the NC-rich, mixed, and CLC-rich groups could be obtained using young Beagles (25 to 36 months), as IVD degeneration occurs relatively early in life. In contrast, to obtain an appropriate sample size for the non-chondrodystrophic group, IVD material was obtained from young and older dogs, since most IVDs in these non-chondrodystrophic dogs remain rich in NCs and IVD degeneration mainly occurs at older age. Therefore, IVDs graded as mixed and CLC-rich could only be found in relatively older non-chondrodystrophic dogs, explaining the relatively large variation in age.
Apart from the distinct IVD phenotype, Caveolin-1 KO mice also exhibit a distinct bone phenotype . Since the endplates are proposed to play a role in the pathophysiology of IVD degeneration, endplate changes may also influence IVD and NC physiology.
Early (subclinical) degeneration of the IVD, which is characterized by changes in the NP cell population, involves significant changes in the expression of genes involved in canonical Wnt signaling, ultimately leading to downregulation of this pathway. The expression of Caveolin-1, which regulates canonical Wnt signaling, is decreased in the CLC-rich NP and appears to be essential to NC physiology and preservation. In view of the high resemblance between humans and dogs regarding the biochemistry and molecular biology of IVD degeneration, it is concluded that Caveolin-1 may play an important role in IVD aging/degeneration in humans as well. Caveolin-1 may serve as an interesting target for developing novel treatment strategies for IVD degeneration.
bovine serum albumin
Dickkopf homolog 3
Dubecco's modified Eagle's serum
fetal calf serum
false discovery rate
Frizzled related protein
functional spine unit
Gene Expression Omnibus
hematoxylin and eosin
Low density lipoprotein receptor-related protein 5
microarray analysis of variance
quantitative polymerase chain reaction
Wnt inhibitory factor 1.
The authors would like to thank Agria International (Sweden) for funding, Ms Jane Sykes for language corrections, Mr. Hans Vernooij (Utrecht University) for guidance in statistical analysis, and the Multimedia Department of the Faculty of Veterinary Medicine for technical assistance.
- Luoma K, Riihimaki H, Luukkonen R, Raininko R, Viikari-Juntura E, Lamminen A: Low back pain in relation to lumbar disc degeneration. Spine (Phila Pa 1976). 2000, 25: 487-492. 10.1097/00007632-200002150-00016.View ArticleGoogle Scholar
- Bono CM, Lee CK: Critical analysis of trends in fusion for degenerative disc disease over the past 20 years: influence of technique on fusion rate and clinical outcome. Spine (Phila Pa 1976). 2004, 29: 455-463. 10.1097/01.BRS.0000090825.94611.28. discussion Z455View ArticleGoogle Scholar
- Gibson JN, Waddell G: Surgical interventions for lumbar disc prolapse: updated Cochrane Review. Spine (Phila Pa 1976). 2007, 32: 1735-1747. 10.1097/BRS.0b013e3180bc2431.View ArticleGoogle Scholar
- Brinckmann P, Grootenboer H: Change of disc height, radial disc bulge, and intradiscal pressure from discectomy. An in vitro investigation on human lumbar discs. Spine (Phila Pa 1976). 1991, 16: 641-646. 10.1097/00007632-199106000-00008.View ArticleGoogle Scholar
- Eck JC, Humphreys SC, Hodges SD: Adjacent-segment degeneration after lumbar fusion: a review of clinical, biomechanical, and radiologic studies. Am J Orthop. 1999, 28: 336-340.PubMedGoogle Scholar
- Frei H, Oxland TR, Rathonyi GC, Nolte LP: The effect of nucleotomy on lumbar spine mechanics in compression and shear loading. Spine (Phila Pa 1976). 2001, 26: 2080-2089. 10.1097/00007632-200110010-00007.View ArticleGoogle Scholar
- Goel VK, Nishiyama K, Weinstein JN, Liu YK: Mechanical properties of lumbar spinal motion segments as affected by partial disc removal. Spine (Phila Pa 1976). 1986, 11: 1008-1012. 10.1097/00007632-198612000-00007.View ArticleGoogle Scholar
- Lee CK: Accelerated degeneration of the segment adjacent to a lumbar fusion. Spine (Phila Pa 1976). 1988, 13: 375-377. 10.1097/00007632-198803000-00029.View ArticleGoogle Scholar
- MacDougall J, Perra J, Pinto M: Incidence of adjacent segment degeneration at ten years after lumbar spine fusion. Spine J. 2003, 3: 675-View ArticleGoogle Scholar
- Tow BP, Hsu WK, Wang JC: Disc regeneration: a glimpse of the future. Clin Neurosurg. 2007, 54: 122-128.PubMedGoogle Scholar
- Erwin WM: The enigma that is the nucleus pulposus cell: the search goes on. Arthritis Res Ther. 2010, 12: 118-10.1186/ar3001.PubMed CentralView ArticlePubMedGoogle Scholar
- Risbud MV, Schaer TP, Shapiro IM: Toward an understanding of the role of notochordal cells in the adult intervertebral disc: from discord to accord. Dev Dyn. 2010, 239: 2141-2148. 10.1002/dvdy.22350.PubMed CentralView ArticlePubMedGoogle Scholar
- Risbud MV, Shapiro IM: Notochordal cells in the adult intervertebral disc: new perspective on an old question. Crit Rev Eukaryot Gene Expr. 2011, 21: 29-41. 10.1615/CritRevEukarGeneExpr.v21.i1.30.PubMed CentralView ArticlePubMedGoogle Scholar
- Bergknut N, Rutges JP, Kranenburg HJ, Smolders LA, Hagman R, Smidt HJ, Lagerstedt AS, Voorhout G, Hazewinkel HAW, Grinwis GC, Creemers LB, Meij BP, Dhert WJ: The dog as an animal model for intervertebral disc degeneration?. Spine (Phila Pa 1976). 2012, 37: 351-358. 10.1097/BRS.0b013e31821e5665.View ArticleGoogle Scholar
- Hunter CJ, Matyas JR, Duncan NA: The functional significance of cell clusters in the notochordal nucleus pulposus: survival and signaling in the canine intervertebral disc. Spine (Phila Pa 1976). 2004, 29: 1099-1104. 10.1097/00007632-200405150-00010.View ArticleGoogle Scholar
- Trout JJ, Buckwalter JA, Moore KC, Landas SK: Ultrastructure of the human intervertebral disc. I. Changes in notochordal cells with age. Tissue Cell. 1982, 14: 359-369. 10.1016/0040-8166(82)90033-7.View ArticlePubMedGoogle Scholar
- Erwin WM, Inman RD: Notochord cells regulate intervertebral disc chondrocyte proteoglycan production and cell proliferation. Spine (Phila Pa 1976). 2006, 31: 1094-1099. 10.1097/01.brs.0000216593.97157.dd.View ArticleGoogle Scholar
- Purmessur D, Schek RM, Abbott RD, Ballif BA, Godburn KE, Iatridis JC: Notochordal conditioned media from tissue increases proteoglycan accumulation and promotes a healthy nucleus pulposus phenotype in human mesenchymal stem cells. Arthritis Res Ther. 2011, 13: R81-10.1186/ar3344.PubMed CentralView ArticlePubMedGoogle Scholar
- Erwin WM, Islam D, Inman RD, Fehlings MG, Tsui FW: Notochordal cells protect nucleus pulposus cells from degradation and apoptosis: implications for the mechanisms of intervertebral disc degeneration. Arthritis Res Ther. 2011, 13: R215-10.1186/ar3548.PubMed CentralView ArticlePubMedGoogle Scholar
- Erwin WM, Las Heras F, Islam D, Fehlings MG, Inman RD: The regenerative capacity of the notochordal cell: tissue constructs generated in vitro under hypoxic conditions. J Neurosurg Spine. 2009, 10: 513-521. 10.3171/2009.2.SPINE08578.View ArticlePubMedGoogle Scholar
- Erwin WM, Inman RD: Notochord cells regulate intervertebral disc chondrocyte proteoglycan production and cell proliferation. Spine (Phila Pa 1976). 2006, 31: 1094-1099. 10.1097/01.brs.0000216593.97157.dd.View ArticleGoogle Scholar
- Hansen HJ: A pathologic-anatomical study on disc degeneration in dog, with special reference to the so-called enchondrosis intervertebralis. Acta Orthop Scand Suppl. 1952, 11: 1-117.View ArticlePubMedGoogle Scholar
- Braund KG, Ghosh P, Taylor TK, Larsen LH: Morphological studies of the canine intervertebral disc. The assignment of the beagle to the achondroplastic classification. Res Vet Sci. 1975, 19: 167-172.PubMedGoogle Scholar
- da Costa RC, Parent JM, Partlow G, Dobson H, Holmberg DL, Lamarre J: Morphologic and morphometric magnetic resonance imaging features of Doberman Pinschers with and without clinical signs of cervical spondylomyelopathy. Am J Vet Res. 2006, 67: 1601-1612. 10.2460/ajvr.67.9.1601.View ArticlePubMedGoogle Scholar
- De Decker S, Gielen IM, Duchateau L, Van Soens I, Bavegems V, Bosmans T, van Bree HJ, Van Ham LM: Low-field magnetic resonance imaging findings of the caudal portion of the cervical region in clinically normal Doberman Pinschers and Foxhounds. Am J Vet Res. 2010, 71: 428-434. 10.2460/ajvr.71.4.428.View ArticlePubMedGoogle Scholar
- Rossi F, Seiler G, Busato A, Wacker C, Lang J: Magnetic resonance imaging of articular process joint geometry and intervertebral disk degeneration in the caudal lumbar spine (L5-S1) of dogs with clinical signs of cauda equina compression. Vet Radiol Ultrasound. 2004, 45: 381-387. 10.1111/j.1740-8261.2004.04083.x.View ArticlePubMedGoogle Scholar
- Seiler GS, Häni H, Busato AR, Lang J: Facet joint geometry and intervertebral disc degeneration in the L5-S1 region of the vertebral column in German Shepherd Dogs. Am J Vet Res. 2002, 63: 86-90. 10.2460/AJVR.2002.63.86.View ArticlePubMedGoogle Scholar
- Johnson JA, da Costa RC, Bhattacharya S, Goel V, Allen MJ: Kinematic motion patterns of the cranial and caudal canine cervical spine. Vet Surg. 2011, 40: 720-727. 10.1111/j.1532-950X.2011.00853.x.View ArticlePubMedGoogle Scholar
- Breit S, Kunzel W: Shape and orientation of articular facets of cervical vertebrae (C3-C7) in dogs denoting axial rotational ability: an osteological study. Eur J Morphol. 2002, 40: 43-51.PubMedGoogle Scholar
- Erwin WM: The Notochord, Notochordal cell and CTGF/CCN-2: ongoing activity from development through maturation. J Cell Commun Signal. 2008, 2: 59-65. 10.1007/s12079-008-0031-5.PubMed CentralView ArticlePubMedGoogle Scholar
- Vos P, Interma F, van El B, DeGroot J, Bijlsma JW, Lafeber F, Mastbergen S: Does loading influence the severity of cartilage degeneration in the canine Groove-model of OA?. J Orthop Res. 2009, 27: 1332-1338. 10.1002/jor.20897.View ArticlePubMedGoogle Scholar
- Vos PA, Degroot J, Barten-van Rijbroek AD, Zuurmond AM, Bijlsma JW, Mastbergen SC, Lafeber FP: Elevation of cartilage AGEs does not accelerate initiation of canine experimental osteoarthritis upon mild surgical damage. J Orthop Res. 2012, 30: 1398-1404. 10.1002/jor.22092.View ArticlePubMedGoogle Scholar
- Van Meegeren MER, Roosendaal G, Barten-van Rijbroek AD, Schutgens REG, Lafeber FPJG, Mastbergen S: Coagulation aggravates blood-induced joint damage. Arthritis & Rheumatism. 2012, 64: 3231-3239. 10.1002/art.34552.View ArticleGoogle Scholar
- van Horssen P, Siebes M, Hoefer I, Spaan JA, van den Wijngaard JP: Improved detection of fluorescently labeled microspheres and vessel architecture with an imaging cryomicrotome. Med Biol Eng Comput. 2010, 48: 735-744. 10.1007/s11517-010-0652-8.PubMed CentralView ArticlePubMedGoogle Scholar
- Smolders LA, Meij BP, Riemers FM, Licht R, Wubbolts R, Heuvel D, Grinwis GC, Vernooij HC, Hazewinkel HA, Penning LC, Tryfonidou MA: Canonical Wnt signaling in the notochordal cell is upregulated in early intervertebral disk degeneration. J Orthop Res. 2012, 30: 950-957. 10.1002/jor.22000.View ArticlePubMedGoogle Scholar
- Hunter CJ, Matyas JR, Duncan NA: The three-dimensional architecture of the notochordal nucleus pulposus: novel observations on cell structures in the canine intervertebral disc. J Anat. 2003, 202: 279-291. 10.1046/j.1469-7580.2003.00162.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Herrmann BG, Kispert A: The T genes in embryogenesis. Trends Genet. 1994, 10: 280-286. 10.1016/0168-9525(90)90011-T.View ArticlePubMedGoogle Scholar
- Gilson A, Dreger M, Urban JP: Differential expression level of cytokeratin 8 in cells of the bovine nucleus pulposus complicates the search for specific intervertebral disc cell markers. Arthritis Res Ther. 2010, 12: R24-10.1186/ar2931.PubMed CentralView ArticlePubMedGoogle Scholar
- Gotz W, Kasper M, Fischer G, Herken R: Intermediate filament typing of the human embryonic and fetal notochord. Cell Tissue Res. 1995, 280: 455-462. 10.1007/BF00307819.View ArticlePubMedGoogle Scholar
- Opitz L, Salinas-Riester G, Grade M, Jung K, Jo P, Emons G, Ghadimi BM, Beissbarth T, Gaedcke J: Impact of RNA degradation on gene expression profiling. BMC Med Genomics. 2010, 3: 36-10.1186/1755-8794-3-36.PubMed CentralView ArticlePubMedGoogle Scholar
- Holstege Lab, Microarray Facility. [http://microarrays.holstegelab.nl/modx/index.php]
- Roepman P, de Koning E, van Leenen D, de Weger RA, Kummer JA, Slootweg PJ, Holstege FC: Dissection of a metastatic gene expression signature into distinct components. Genome Biol. 2006, 7: R117-10.1186/gb-2006-7-12-r117.PubMed CentralView ArticlePubMedGoogle Scholar
- van de Peppel J, Kemmeren P, van Bakel H, Radonjic M, van Leenen D, Holstege FC: Monitoring global messenger RNA changes in externally controlled microarray experiments. EMBO Rep. 2003, 4: 387-393. 10.1038/sj.embor.embor798.PubMed CentralView ArticlePubMedGoogle Scholar
- Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J, Speed TP: Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Res. 2002, 30: e15-10.1093/nar/30.4.e15.PubMed CentralView ArticlePubMedGoogle Scholar
- Margaritis T, Lijnzaad P, van Leenen D, Bouwmeester D, Kemmeren P, van Hooff SR, Holstege FC: Adaptable gene-specific dye bias correction for two-channel DNA microarrays. Mol Syst Biol. 2009, 5: 266-PubMed CentralView ArticlePubMedGoogle Scholar
- Wu H, Kerr MK, Cui X, Churchill GA: MAANOVA: a software package for the analysis of spotted cDNA microarray experiments. The analysis of gene expression data: methods and software. Edited by: Parmigiani G, Garett ES, Orozarru RA, Zeger SL. 2003, New York: Springer-Verlag, i-xlviii.Google Scholar
- Ekins S, Nikolsky Y, Bugrim A, Kirillov E, Nikolskaya T: Pathway mapping tools for analysis of high content data. Methods Mol Biol. 2007, 356: 319-350.PubMedGoogle Scholar
- Roepman P, Wessels LF, Kettelarij N, Kemmeren P, Miles AJ, Lijnzaad P, Tilanus MG, Koole R, Hordijk GJ, van der Vliet PC, Reinders MJ, Slootweg PJ, Holstege FC: An expression profile for diagnosis of lymph node metastases from primary head and neck squamous cell carcinomas. Nat Genet. 2005, 37: 182-186. 10.1038/ng1502.View ArticlePubMedGoogle Scholar
- Hsieh JC, Kodjabachian L, Rebbert ML, Rattner A, Smallwood PM, Samos CH, Nusse R, Dawid IB, Nathans J: A new secreted protein that binds to Wnt proteins and inhibits their activities. Nature. 1999, 398: 431-436. 10.1038/18899.View ArticlePubMedGoogle Scholar
- Cadigan KM, Nusse R: Wnt signaling: a common theme in animal development. Genes Dev. 1997, 11: 3286-3305. 10.1101/gad.11.24.3286.View ArticlePubMedGoogle Scholar
- Mao J, Wang J, Liu B, Pan W, Farr GH, Flynn C, Yuan H, Takada S, Kimelman D, Li L, Wu D: Low-density lipoprotein receptor-related protein-5 binds to Axin and regulates the canonical Wnt signaling pathway. Mol Cell. 2001, 7: 801-809. 10.1016/S1097-2765(01)00224-6.View ArticlePubMedGoogle Scholar
- McDonald PC, Fielding AB, Dedhar S: Integrin-linked kinase-essential roles in physiology and cancer biology. J Cell Sci. 2008, 121: 3121-3132. 10.1242/jcs.017996.View ArticlePubMedGoogle Scholar
- Zhurinsky J, Shtutman M, Ben-Ze'ev A: Plakoglobin and beta-catenin: protein interactions, regulation and biological roles. J Cell Sci. 2000, 113 (Pt 18): 3127-3139.PubMedGoogle Scholar
- Jho EH, Zhang T, Domon C, Joo CK, Freund JN, Costantini F: Wnt/beta-catenin/Tcf signaling induces the transcription of Axin2, a negative regulator of the signaling pathway. Mol Cell Biol. 2002, 22: 1172-1183. 10.1128/MCB.22.4.1172-1183.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Galbiati F, Volonte D, Brown AM, Weinstein DE, Ben-Ze'ev A, Pestell RG, Lisanti MP: Caveolin-1 expression inhibits Wnt/beta-catenin/Lef-1 signaling by recruiting beta-catenin to caveolae membrane domains. J Biol Chem. 2000, 275: 23368-23377. 10.1074/jbc.M002020200.View ArticlePubMedGoogle Scholar
- Yamamoto H, Komekado H, Kikuchi A: Caveolin is necessary for Wnt-3a-dependent internalization of LRP6 and accumulation of beta-catenin. Dev Cell. 2006, 11: 213-223. 10.1016/j.devcel.2006.07.003.View ArticlePubMedGoogle Scholar
- Nusse R: Developmental biology. Making head or tail of Dickkopf. Nature. 2001, 411: 255-256. 10.1038/35077199.View ArticlePubMedGoogle Scholar
- Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29: e45-10.1093/nar/29.9.e45.PubMed CentralView ArticlePubMedGoogle Scholar
- Pinheiro J, Bates D, DebRoy S, Sarkar D, the R Core team: nlme: linear and nonlinear mixed effects models. R package version 3. 2009, 1-96.Google Scholar
- R Development Core Team: R: A language and environment for statistical computing. R Foundation for Statistical Computing. 2010, Vienna: AustriaGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Stat Soc B. 1995, 57: 289-300.Google Scholar
- Baroni MG, Cavallo MG, Mark M, Monetini L, Stoehrer B, Pozzilli P: Beta-cell gene expression and functional characterisation of the human insulinoma cell line CM. J Endocrinol. 1999, 161: 59-68. 10.1677/joe.0.1610059.View ArticlePubMedGoogle Scholar
- Gruber HE, Ingram J, Hanley EN: An improved staining method for intervertebral disc tissue. Biotech Histochem. 2002, 77: 81-83.View ArticlePubMedGoogle Scholar
- Stosiek P, Kasper M, Karsten U: Expression of cytokeratin and vimentin in nucleus pulposus cells. Differentiation. 1988, 39: 78-81. 10.1111/j.1432-0436.1988.tb00083.x.View ArticlePubMedGoogle Scholar
- Salisbury JR, Isaacson PG: Demonstration of cytokeratins and an epithelial membrane antigen in chordomas and human fetal notochord. Am J Surg Pathol. 1985, 9: 791-797. 10.1097/00000478-198511000-00002.View ArticlePubMedGoogle Scholar
- Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30: 207-210. 10.1093/nar/30.1.207.PubMed CentralView ArticlePubMedGoogle Scholar
- Gene expression Omnibus (GEO) data for LS001A: Aging/Early Degeneration Of The Intervertebral Disc Involves Downregulation Of Canonical Wnt Signaling and Caveolin-1 expression. 666, [http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=phonhywsqeosgra&acc=GSE35717]
- Corr M: Wnt-beta-catenin signaling in the pathogenesis of osteoarthritis. Nat Clin Pract Rheumatol. 2008, 4: 550-556.View ArticlePubMedGoogle Scholar
- Yan D, Wiesmann M, Rohan M, Chan V, Jefferson AB, Guo L, Sakamoto D, Caothien RH, Fuller JH, Reinhard C, Garcia PD, Randazzo FM, Escobedo J, Fanti WJ, Williams LT: Elevated expression of axin2 and hnkd mRNA provides evidence that Wnt/beta -catenin signaling is activated in human colon tumors. Proc Natl Acad Sci USA. 2001, 98: 14973-14978. 10.1073/pnas.261574498.PubMed CentralView ArticlePubMedGoogle Scholar
- Lustig B, Jerchow B, Sachs M, Weiler S, Pietsch T, Karsten U, van de Wetering M, Clevers H, Schlag PM, Birchmeier W, Behrens J: Negative feedback loop of Wnt signaling through upregulation of conductin/axin2 in colorectal and liver tumors. Mol Cell Biol. 2002, 22: 1184-1193. 10.1128/MCB.22.4.1184-1193.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Jho EH, Zhang T, Domon C, Joo CK, Freund JN, Costantini F: Wnt/beta-catenin/Tcf signaling induces the transcription of Axin2, a negative regulator of the signaling pathway. Mol Cell Biol. 2002, 22: 1172-1183. 10.1128/MCB.22.4.1172-1183.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang Z, Shu W, Lu MM, Morrisey EE: Wnt7b activates canonical signaling in epithelial and vascular smooth muscle cells through interactions with Fzd1, Fzd10, and LRP5. Mol Cell Biol. 2005, 25: 5022-5030. 10.1128/MCB.25.12.5022-5030.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Chun J, Hyun S, Kwon T, Lee EJ, Hong SK, Kang SS: The subcellular localization control of integrin linked kinase 1 through its protein-protein interaction with caveolin-1. Cell Signal. 2005, 17: 751-760. 10.1016/j.cellsig.2004.10.016.View ArticlePubMedGoogle Scholar
- Heathfield SK, Le Maitre CL, Hoyland JA: Caveolin-1 expression and stress-induced premature senescence in human intervertebral disc degeneration. Arthritis Res Ther. 2008, 10: R87-10.1186/ar2468.PubMed CentralView ArticlePubMedGoogle Scholar
- Erwin WM, Islam D, Eftekarpour E, Inman RD, Karim MZ, Fehlings MG: Intervertebral disc-derived stem cells: implications for regenerative medicine and neural repair. Spine (Phila Pa 1976). 2013, 38: 211-216. 10.1097/BRS.0b013e318266a80d.View ArticleGoogle Scholar
- Blanco JF, Graciani IF, Sanchez-Guijo FM, Muntion S, Hernandez-Campo P, Santamaria C, Carrancio S, Barbado MV, Cruz G, Gutierrez-Cosio S, Herrero C, San Miguel JF, Briñon JG, del Cañizo MC: Isolation and characterization of mesenchymal stromal cells from human degenerated nucleus pulposus: comparison with bone marrow mesenchymal stromal cells from the same subjects. Spine (Phila Pa 1976). 2010, 35: 2259-2265. 10.1097/BRS.0b013e3181cb8828.View ArticleGoogle Scholar
- Risbud MV, Guttapalli A, Tsai TT, Lee JY, Danielson KG, Vaccaro AR, Albert TJ, Gazit Z, Gazit D, Shapiro IM: Evidence for skeletal progenitor cells in the degenerate human intervertebral disc. Spine (Phila Pa 1976). 2007, 32: 2537-2544. 10.1097/BRS.0b013e318158dea6.View ArticleGoogle Scholar
- Minogue BM, Richardson SM, Zeef LA, Freemont AJ, Hoyland JA: Transcriptional profiling of bovine intervertebral disc cells: implications for identification of normal and degenerate human intervertebral disc cell phenotypes. Arthritis Res Ther. 2010, 12: R22-10.1186/ar2929.PubMed CentralView ArticlePubMedGoogle Scholar
- Aguiar DJ, Johnson SL, Oegema TR: Notochordal cells interact with nucleus pulposus cells: regulation of proteoglycan synthesis. Exp Cell Res. 1999, 246: 129-137. 10.1006/excr.1998.4287.View ArticlePubMedGoogle Scholar
- Erwin WM, Ashman K, O'Donnel P, Inman RD: Nucleus pulposus notochord cells secrete connective tissue growth factor and up-regulate proteoglycan expression by intervertebral disc chondrocytes. Arthritis Rheum. 2006, 54: 3859-3867. 10.1002/art.22258.View ArticlePubMedGoogle Scholar
- Cappello R, Bird JL, Pfeiffer D, Bayliss MT, Dudhia J: Notochordal cell produce and assemble extracellular matrix in a distinct manner, which may be responsible for the maintenance of healthy nucleus pulposus. Spine (Phila Pa 1976). 2006, 31: 873-882. 10.1097/01.brs.0000209302.00820.fd.View ArticleGoogle Scholar
- Nakamura Y, Nawata M, Wakitani S: Expression profiles and functional analyses of Wnt-related genes in human joint disorders. Am J Pathol. 2005, 167: 97-105. 10.1016/S0002-9440(10)62957-4.PubMed CentralView ArticlePubMedGoogle Scholar
- Ukita K, Hirahara S, Oshima N, Imuta Y, Yoshimoto A, Jang CW, Oginuma M, Saga Y, Behringer RR, Kondoh H, Sasaki H: Wnt signaling maintains the notochord fate for progenitor cells and supports the posterior extension of the notochord. Mech Dev. 2009, 126: 791-803. 10.1016/j.mod.2009.08.003.PubMed CentralView ArticlePubMedGoogle Scholar
- Cho HH, Kim YJ, Kim SJ, Kim JH, Bae YC, Ba B, Jung JS: Endogenous Wnt signaling promotes proliferation and suppresses osteogenic differentiation in human adipose derived stromal cells. Tissue Eng. 2006, 12: 111-121. 10.1089/ten.2006.12.111.View ArticlePubMedGoogle Scholar
- Corr M: Wnt-beta-catenin signaling in the pathogenesis of osteoarthritis. Nat Clin Pract Rheumatol. 2008, 4: 550-556.View ArticlePubMedGoogle Scholar
- Sen M: Wnt signalling in rheumatoid arthritis. Rheumatology (Oxford). 2005, 44: 708-713. 10.1093/rheumatology/keh553.View ArticleGoogle Scholar
- Wang M, Tang D, Shu B, Wang B, Jin H, Hao S, Dresser KA, Shen J, Im HJ, Sampson ER, Rubery PT, Zuscik MJ, Schwarz EM, O'Keefe RJ, Wang Y, Chen D: Conditional activation of beta-catenin signaling in mice leads to severe defects in intervertebral disc tissue. Arthritis Rheum. 2012, 64: 2611-2623. 10.1002/art.34469.PubMed CentralView ArticlePubMedGoogle Scholar
- Nixon SJ, Carter A, Wegner J, Ferguson C, Floetenmeyer M, Riches J, Key B, Westerfield M, Parton RG: Caveolin-1 is required for lateral line neuromast and notochord development. J Cell Sci. 2007, 120: 2151-2161. 10.1242/jcs.003830.View ArticlePubMedGoogle Scholar
- Heathfield SK, Hoyland JA: Interleukin-1 induces caveolin-1 expression in human nucleus pulposus cells [abstract]. J Bone Joint Surg Br. 2011, 93-B: 484-10.1302/0301-620X.93B4.25899.View ArticleGoogle Scholar
- Dai SM, Shan ZZ, Nakamura H, Masuko-Hongo K, Kato T, Nishioka K, Yudoh K: Catabolic stress induces features of chondrocyte senescence through overexpression of caveolin 1: possible involvement of caveolin 1-induced down-regulation of articular chondrocytes in the pathogenesis of osteoarthritis. Arthritis Rheum. 2006, 54: 818-831. 10.1002/art.21639.View ArticlePubMedGoogle Scholar
- Mercier I, Jasmin JF, Pavlides S, Minetti C, Flomenberg N, Pestell RG, Frank PG, Sotgia F, Lisanti MP: Clinical and translational implications of the caveolin gene family: lessons from mouse models and human genetic disorders. Lab Invest. 2009, 89: 614-623. 10.1038/labinvest.2009.23.PubMed CentralView ArticlePubMedGoogle Scholar
- Boscher C, Nabi IR: Caveolin-1: role in cell signaling. Adv Exp Med Biol. 2012, 729: 29-50. 10.1007/978-1-4614-1222-9_3.View ArticlePubMedGoogle Scholar
- Li WP, Liu P, Pilcher BK, Anderson RG: Cell-specific targeting of caveolin-1 to caveolae, secretory vesicles, cytoplasm or mitochondria. J Cell Sci. 2001, 114: 1397-1408.PubMedGoogle Scholar
- Anderson RG: Caveolae: where incoming and outgoing messengers meet. Proc Natl Acad Sci USA. 1993, 90: 10909-10913. 10.1073/pnas.90.23.10909.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu P, Rudick M, Anderson RG: Multiple functions of caveolin-1. J Biol Chem. 2002, 277: 41295-41298. 10.1074/jbc.R200020200.View ArticlePubMedGoogle Scholar
- Rubin J, Schwartz Z, Boyan BD, Fan X, Case N, Sen B, Drab M, Smith D, Aleman M, Wong KL, Yao H, Jo H, Gross TS: Caveolin-1 knockout mice have increased bone size and stiffness. J Bone Miner Res. 2007, 22: 1408-1418. 10.1359/jbmr.070601.View ArticlePubMedGoogle Scholar
- Mudunuri U, Che A, Yi M, Stephens RM: bioDBnet: the biological database network. Bioinformatics. 2009, 25: 555-556. 10.1093/bioinformatics/btn654.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.