Skip to main content

Differential gene expression of human chondrocytes cultured under short-term altered gravity conditions during parabolic flight maneuvers

Abstract

Background

Chondrocytes are the main cellular component of articular cartilage. In healthy tissue, they are embedded in a strong but elastic extracelluar matrix providing resistance against mechanical forces and friction for the joints. Osteoarthritic cartilage, however, disrupted by heavy strain, has only very limited potential to heal. One future possibility to replace damaged cartilage might be the scaffold-free growth of chondrocytes in microgravity to form 3D aggregates.

Results

To prepare for this, we have conducted experiments during the 20th DLR parabolic flight campaign, where we fixed the cells after the first (1P) and the 31st parabola (31P). Furthermore, we subjected chondrocytes to isolated vibration and hypergravity conditions. Microarray and quantitative real time PCR analyses revealed that hypergravity regulated genes connected to cartilage integrity (BMP4, MMP3, MMP10, EDN1, WNT5A, BIRC3). Vibration was clearly detrimental to cartilage (upregulated inflammatory IL6 and IL8, downregulated growth factors EGF, VEGF, FGF17). The viability of the cells was not affected by the parabolic flight, but showed a significantly increased expression of anti-apoptotic genes after 31 parabolas. The IL-6 release of chondrocytes cultured under conditions of vibration was not changed, but hypergravity (1.8 g) induced a clear elevation of IL-6 protein in the supernatant compared with corresponding control samples.

Conclusion

Taken together, this study provided new insights into the growth behavior of chondrocytes under short-term microgravity.

Background

Joint friction at the extremities of long bones is reduced by articular cartilage. This kind of tissue is highly specialized, avascular, not innervated and consists mainly of a single cell type: the chondrocytes. The chondrocytes are tightly embedded in an extracellular matrix (ECM), which is composed of a network of collagens (predominantly collagen type II) and aggrecan. The collagen network contributes to the strength and mechanical resistance of cartilage tissue, whereas the aggrecan, a proteoglycan, comprising such molecules as chondroitin sulfate or keratin sulfate, is responsible for its viscoelasticity and flexibility due to its ability to absorb and retain considerable amounts of water [1-3]. Degenerative diseases of the cartilage like osteoarthritis are characterized by a progressive degradation of the ECM, caused by the increased secretion of matrix metalloproteinases (MMP) [4,5]. This process is triggered by pro-inflammatory cytokines, such as tumor necrosis factor-α (TNF-α) or interleukin-1β (IL-1β), which are transported into the cartilage via the synovial fluid [6,7]. The absence of vasculature and the extremely limited influx of chondrocyte progenitor cells [8] limit the tissues healing and restorative potential, so that in advanced stages of osteoarthritis a surgical replacement of the affected joint with a prosthesis is usually necessary.

Apart from the influences of age and use intensity, other factors have been found to mediate cartilage integrity. Most notably, microgravity (μg) has a strong impact on cartilage. It has been reported, that astronauts, after staying a longer time in Space, suffer from a reduction of cartilage mass [9] due to mechanical unloading. In addition, when cultured in Space, cartilage tissue showed a reduced aggrecan density and a less-organized collagen subtype 2 organization hinting towards an impaired resistance to mechanical stress and breaking [10,11]. It is therefore of high interest, to study the impact of microgravity on chondrocytes. On the one hand, this will help to understand the detrimental effects of prolonged stays in Space on the cartilage, but on the other it might also help to find ways to counteract this phenomenon in Space as well as to ameliorate cartilage problems caused by wear on Earth in the future.

Microgravity (μg) provides unique conditions for cell and tissue growth. It has been shown on various cell types, including chondrocytes, that cultivation under conditions of μg can induce the formation of 3D aggregates. These aggregates are especially interesting, as they do not require any potentially interfering scaffolding like those being generated under normal gravity conditions [12-16].

Cell cultivation in Space, however, is an extremely complicated and expensive venture. Therefore, pre-studies such as a parabolic flight, providing 31×22 s of short-term real μg, or simulated μg on ground-based facilities such as the rotating wall vessel bioreactor, clinostats or the random positioning machine (RPM) have been established [17-19]. Especially the parabolic flight is an attractive method to achieve real μg without going into Space. It should be taken into consideration, however, that every μg-phase is flanked by two 20 s-lasting hypergravity phases of 1.8 g. Moreover, during the flight vibrations occur and have to be taken into account for interpretation of the results [20].

This study aimed to investigate the influence of short-term real μg during a parabolic flight campaign on the gene expression profiles of cultivated chondrocytes. In addition, the effects of 1.8 g hypergravity as well as of vibration in an extent comparable to those during a parabolic flight were separately investigated.

Results

Influence of hypergravity on chondrocyte gene expression

Using the Microarray analysis technique, we found a total of 210 genes (top 30 are listed in Table 1, for a complete list see Additional file 1) which were significantly differentially expressed and showed a Fold-Change of >2 or < −2 under conditions of hypergravity of 1.8 g. An enrichment analysis for Gene Ontology (GO) Biological Processes (BP) terms (top 15 given in Table 2, for a complete list see Additional file 1) revealed that mainly development-related processes were detected. GO:0048598 (embryonic morphogenesis, p = 2.28×10−5), GO:0001501 (skeletal system development, p = 3.33×10−5), and GO:0048729 (tissue morphogenesis, p = 3.83×10−5) were the most prominent. However, many genes belonging to these BP are also implicated in cartilage development (GO:0051216, p = 5.59×10−4) or cell adhesion (GO:0007155, p = 5.96×10−4) and biological adhesion (GO:0022610, p = 6.05×10−4), such as the bone morphogenetic protein 4 (BMP4), fibroblast growth factor 9 (FGF9), hyaluronan and proteoglycan link protein 1 (HAPLN1), collagen type 2 α 1 (COL2A1), or laminin α 5 (LAMA5).

Table 1 Top 30 most differentially expressed genes under hypergravity as detected by microarray analysis
Table 2 Top 15 significantly enriched Gene Ontology biological processes under hypergravity

Influence of vibration on chondrocyte gene expression

Selected genes were analysed by quantitative real-time PCR. After 2 hours of vibration, gene expression of interleukin-6 (IL6, Figure 1B) and −8 (IL8, Figure 1D), were significantly elevated, while epidermal growth factor (EGF, Figure 1E), vascular endothelial growth factor D (VEGFD, Figure 1J), and fibroblast growth factor 17 (FGF17, Figure 1K) transcripts were significantly downregulated (p<0.05). No changes were observed for connective tissue growth factor (CTGF, Figure 1A), caveolin 2 (CAV2, Figure 1C), protein kinase, AMP-activated, α 1 (PRKAA, Figure 1F), vascular endothelial growth factor A (VEGFA, Figure 1G), protein kinase c alpha (PRKCA, Figure 1H), ezrin (VIL2, Figure 1I), and interleukin-15 (IL15, Figure 1L).

Figure 1
figure 1

Quantitative real-time PCR analysis of chondrocytes exposed to the Vibraplex device. The Vibraplex device provides vibration profiles corresponding to those occurring during parabolic flight and allows the isolated analysis of their effects on cultivated cells. In all diagrams, the x-axis represents the experiment conditions and the y-axis the relative gene expression in % of control. Results are given as mean value ± SD. Significant changes (p < 0.05) are indicated by brackets above the bars. The analyzed genes were A: CTGF; B: IL6; C: CAV2; D: IL8; E: EGF; F: PRKAA; G: VEGFA; H: PRKCA; I: VIL2; J: VEGFD; K: FGF17; L: IL15.

Influence of hypergravity and vibration on soluble factor release

Hypergravity induced a more than 2-fold increase in the release of IL-6 in the supernatant (Figure 2A). IL-8, EGF, VEGFD and FGF17 concentrations in the culture supernatant of hypergravity samples were below the detection limit.

Figure 2
figure 2

ELISA analysis of IL-6 released in the cell culture supernatant. A: IL-6 concentraion in the cell culture medium (pg/ml) from the hypergravity experiments. A significant increase in the secretion of IL-6 is detectable in cells cultured on the SAHC compared with static 1 g samples. B: IL-6 concentration in the cell culture medium (pg/ml) from the vibration experiments compared to the control sample. No significant changes could be detected. Significant changes (p<0.05) are indicated by brackets above the bars.

ELISA analysis revealved no significant change in the concentration of IL-6 protein in chondrocytes cultured under conditions of vibration (Figure 2B), whereas IL-8, EGF, VEGFD and FGF17 concentrations in the culture supernatant were below the detection limit.

Parabolic flight maneuvers induced expression changes in chondrocytes

The influence of the parabolic flight after 1P and 31 P was investigated. Shortly, the 1 g vs. first parabola (1P) vs. 31st Parabola (31P) set was subjected to an F-test. Resulting significant differential expressed probes (5% FDR) were clustered using k = 6. Individual expression characteristics of the 6 clusters are documented in Figure 3.

Figure 3
figure 3

Heatmap showing the clustering of differentially expressed transcripts for the parabolic flight microarray experiments. Analysis after one (1P) and 31 parabolas (31P); red: strong expression; blue: weak expression.

Influence of the parabolic flight

Gene array analysis

Cluster 1 (Table 3) of the parabolic flight set consists of 112 genes, which were upregulated after both 1P and 31P. The predominant biological processes (Additional file 1, for all following clusters) which were found to be significantly enriched were mostly transcription- and metabolism-related (GO:0045449 regulation of transcription, p = 6.08×10−8; GO:0051252 regulation of RNA metabolic process, p = 1.59×10−7; GO:0006355 regulation of transcription, DNA-dependent, p = 1.48×10−6; GO:0006357 regulation of transcription from RNA polymerase II promoter, p = 8.69×10−6; GO:0006350 transcription, p = 1.09×10−5; GO:0019220 ~ regulation of phosphate metabolic process, p = 6.90×10−5). Interestingly, there was also anti-apoptosis (GO:0006916, p = 1.42×10−4), represented by such genes as vascular endothelial growth factor A (VEGFA), neurogenic locus notch homolog protein 2 (NOTCH2), or high-mobility group protein B1 (HMGB1).

Table 3 Overview of the top 5 significantly enriched Gene Ontology biological processes in the 6 clusters of significantly differentially expressed genes during conditions of parabolic flight

51 genes were downregulated after 31P showed an intermediate state after 1P constitute cluster 2 (Table 3). Only few significantly enriched BP were identified and their composition resembled cluster 1 with transcriptional (GO:0045449 regulation of transcription, p = 6.35×10−4; GO:0006350 transcription, p = 1.44×10−4; GO:0006355 regulation of transcription, DNA-dependent, p = 1.90×10−2), metabolic (GO:0051252 regulation of RNA metabolic process, p = 2.15×10−2), and developmental (GO:0007498 mesoderm development, p = 9.48×10−3) processes.

Cluster 3 (Table 3), comprising of 59 genes which were downregulated after 31P, and also showed very few enriched BP, namely GO:0045449 (regulation of transcription, p = 8.41×10−3), and GO:0006350 (transcription, p = 1.00×10−2).

Cluster 4 (Table 3) represents 37 genes which were downregulated after 1P and showed intermediate effects after 31P. No significantly enriched BP was identified.

Cluster 5 (Table 3), consisting of 59 genes upregulated after 1P only, showed BP involved in transcription (GO:0006350 transcription, p = 2.75×10−5; GO:0045449 regulation of transcription, p = 3.01×10−5), but also regulation of tumor necrosis factor production (GO:0032680, p = 3.57×10−3) and rhythmic processes (GO:0048511 rhythmic process, p = 5.90×10-3; GO:0007623 circadian rhythm, p = 7.73×10−3).

73 genes were upregulated only after 31P were pooled into cluster 6 (Table 3). This cluster differs considerably from the other five, as a strong emphasis of (anti)-apoptotic BP was observed. Genes such as IL6, IL8, baculoviral IAP repeat-containing protein3 (BIRC3), induced myeloid leukemia cell differentiation protein (MCL1), or TNF receptor-associated factor 1 (TRAF1) lead to the enrichment of BP, such as GO:0006916 (anti-apoptosis, p = 1.66×10−5), GO:0043066 (negative regulation of apoptosis, p = 7.39×10−5), GO:0043069 (negative regulation of programmed cell death, p = 8.15×10−5, GO:0060548 (negative regulation of cell death, p = 8.31×10−5), or GO:0006915 (apoptosis, p = 1.15×10−4) (Figure 4).

Figure 4
figure 4

STRING analysis of the cluster 6 from the parabolic flight experiment. Chondrocytes were fixed after parabola 1 and 31 during a parabolic flight. In parallel, corresponding 1 g control samples were prepared. K-mean clustering of the resulting microarray data revealed 6 clusters of differentially expressed genes. Clusters 1–5 revealed mostly unspecific transcriptionally active genes, while cluster 6 showed a strong dominance by anti-apoptotic and cell-proliferative transcripts. Possible interactions of the corresponding proteins were visualized using the STRING software and genes involved in anti-apoptosis and cell proliferation were highlighted with white and black circles, respectively.

Quantitative real-time PCR

In addition to the microarray analysis, we also employed the quantitative real-time PCR technique to validate selected genes of interest. CCNA2 as well as IL8 were significantly upregulated only after 31P (Figure 5A, D). CD44 and TNFA were significantly upregulated only after 1P (Figure 5B, G). VCAM showed a significant downregulation after 31P (Figure 5E). No effects were observed for IL6, EDN1 and FGF9 (Figure 5C, F, H).

Figure 5
figure 5

Quantitative real-time PCR analysis of chondrocytes after exposure to parabolic flight. Chondrocytes were taken on a parabolic flight and fixated at two different timepoints. Subsequently, qPCR analyses were performed on these samples and the corresponding 1 g controls. In all diagrams, the x-axis represents the experiment conditions and the y-axis the relative gene expression in % of control. Results are given as mean value ± SD. Significant changes (p < 0.05) are indicated by brackets above the bars. The analyzed genes were A: CCNA2; B: CD44; C: IL6; D: IL8; E: VCAM1; F: EDN1; G: TNFA; H: FGF9. 1P: fixation after first parabola; 31P: fixation after 31st parabola.

Discussion

In this study we investigated the effects of short-term real μg, continuous hypergravity and vibration on human chondrocytes growing in monolayers. For these aims, we employed parabolic flight maneuvers as well as ground-based devices to expose the cells to isolated acceleration profiles and to vibrations as they occur in a combined manner during the flight conditions.

Short-term hypergravity affects chondrocytes

The microarray analysis of chondrocytes exposed for 2 h to 1.8 g revealed that only a very moderate amount of genes was affected in comparison to the parabolic flight effects. It is interesting to notice, that mainly biological processes were affected, which are involved in tissue morphogenesis or skeletal system development. This is a strong indication, that hypergravity directly affects cartilage development.

It has recently been shown, that mechanical load can induce vascular endothelial growth factor A (VEGF-A) expression [21]. VEGF-A belongs to a family of growth factors comprising VEGF-A, −B, −C, −D, −E, and placenta growth factor (PGF) [22,23] and is one of the key components to control angiogenesis, the development of new vessels from existing ones. VEGF is predominantly found in osteoarthritic cartilage/chondrocytes and contributes further to the disintegration of the cartilage by inducing matrix metalloproteinases, which are able to disrupt the extracellular matrix [24,25]. We have observed a similar tendency for VEGFA gene expression, which was significantly elevated by a factor of 1.65, but was omitted for the GO BP analysis due to our cut-off of a Fold-Change >2. Interestingly, MMP3 and MMP10 transcripts were downregulated under hypergravity. We speculate that this might be due to the fact that the chondrocytes form a monolayer in the culture flasks and are not embedded in an ECM-like matrix as in their physiological environment. This strong MMP downregulation might therefore be an attempt to build up a thick ECM, and a stronger stimulus than the VEGFA induction.

BMP4, on the other hand, coding for the bone morphogenetic protein 4 (BMP-4; Additional file 1) was found to be downregulated under hypergravity. BMP-4 stimulates the synthesis of collagen type 2 and aggrecan and thus, enhaces the production of articular cartilage [26,27]. EDN1 (endothelin 1) gene expression was observed to be enhanced (Additional file 1). It has been reported, that overexpression of endothelin 1 is associated with cartilage degeneration [28]. In contrast to this, the inhibitor of apoptosis, BIRC3 (Additional file 1) [29], was enhanced hinting towards an improved cell survival. Furthermore, the gene expression of wingless-type MMTV integration site family, member 5A (WNT5A) was decreased (Additional file 1). Wnt-5a has been shown to be able to induce cartilage degradation through upregulation of MMPs [30]. All in all, it is obvious that in our setup chondrocytes are sensitive to mechanical stress by hypergravity, but at the moment, no definite answer can be given about the nature of the effect.

Short-term vibration is detrimental to chondrocytes

In contrast to hypergravity, the effects of vibration on cultured chondrocytes were clearer. In our experimental setup, the vibrations which were transmitted into the culture flask also caused the culture medium to stir to a certain dregree, which, as we speculate, resulted in additional shear forces. Shear forces have been shown to have a negative effect on chondrocytes and cartilage [31,32]. It has been reported, that cartilage, that was treated in such a way or was degenerating, produced increased amounts of proinflammatory interleukins, such as IL-6 or IL-8 [33,34]. In our qPCR analysis we found a strong increase of both IL6 and IL8 (Figure 1B + D) gene expression, although no increase in IL-6 secretion (Figure 2), accompanied by decreases of EGF, VEGFD, and FGF17 (Figure 1E, J, K) gene expression. The presence of these factors has been described as beneficial for cartilage development [35,36]. Taken together, our results indicate that vibration drives chondrocytes towards an inflammatory, cartilage destabilzing state.

The influence of parabolic flight maneuvers

The microarray analysis showed, that after only 1P relatively unspecific effects on the cells were observed, mainly connected to transcription. This is an indication that the cells have perceived the change in gravity and that they were preparing their transcriptional apparatus for an altered gene expression as a reaction to this stimulus. After 31P, we observed an increase in the enrichment of anti-apoptotic genes. The qPCR analysis reflects the same tendency. Most of the investigated genes showed only transitional or no changes, such as CD44, IL6, EDN1, TNFA, and FGF9 (Figure 5B, C, F, G, H). This seems to hint toward a short μg “shock” that the cells are able to overcome very quickly. Only CCNA2 (Figure 5A), a cyclin involved in cell cycling and proliferaton [37] and the anti-apoptotic IL8 (Figure 5D) [38] are expressed in a manner that they exert a growth-promoting, cell-survival effect. It should be kept in mind, that these effects originate from only a short-term altered gravity (PFC) treatment and that longer exposure times have to be investigated in order to assess their significance. It is interesting to note that RPM exposure experiments resulted in increased expression of several genes responsible for cell motility, structure and integrity; control of cell growth, cell proliferation, cell differentiation and apoptosis [39] and that these results are also in very good accordance with earlier studies that also reported that chondrocytes are quite robust under μg stress [40].

Conclusions

We have shown that chondrocytes are very robust under conditions of parabolic flight maneuvers. They are able to adapt quickly to this new environment and actually profit from real μg by reducing their apoptotic rate. However, they are prone to damage/injury by hypergravity and especially by vibration/shear forces. All in all, these results are very promising and are a step further along the way to understand chondrocyte growth in μg, leading perhaps to new methods of scaffold-free preparation of cartilage grafts.

Materials and methods

Parabolic flight

The parabolic flight experiments were conducted aboard the Airbus A300 ZERO-G operated by Novespace and based in Bordeaux-Merignac, France [19,41-43]. Standard parabolic flights were performed, each with 31 parabolas in a row during the three to four hours flight. The flight manoeuver starts from the horizontal flight level followed by a 45° ascent for 20 s. During this time 1.5 g to 1.8 g are acting on the passengers and the experiments. Then the thrust is reduced and the aircraft follows the path of a parabola. The free fall (microgravity) phase starts and persists for 22 s. Afterwards, the engines are fully powered again and another phase of 1.8 g of 20 s terminates the parabola. Due to the aerodynamic forces and turbulences acting on the aircraft, the μg quality is in the range of about 10−2g.

Cells and cell culture medium

Commercially available human chondrocytes (Provitro®, Berlin, Germany) were cultured in Chondrocyte Growth Medium basal (CGM, Provitro®, Berlin, Germany) supplemented with 10% fetal calf serum (Provitro®, Berlin, Germany), 100 IU penicillin/mL and 100 μg streptomycin/mL (Provitro®, Berlin, Germany).

Cell culture procedure

Cells were grown as published recently [19,44]. Briefly, the cells were cultured in 24 T75 cell culture flasks (75 cm2; Sarstedt, Nümbrecht, Germany) until subconfluent monolayers were obtained. During this time, the cells were covered by 20 mL (T-75 flasks) CGM. One half of the flasks was used as ground control cells (1 g; n = 12), cultured and fixed further in the laboratory, the other half was taken on the parabolic flight (n = 12). During the parabolic flight, RNAlater (Applied Biosystems, Darmstadt, Germany) was injected via syringes containing the appropriate fixative. The syringes were connected to the T-75 flasks through a flexible tube and a 3-way-valve. One hour before each flight, the cell culture flasks were transported to the aircraft and placed into the 37°C preheated incubator on an experimental rack (Figure 6A).

Figure 6
figure 6

Overview of the facilities used to expose cells to parabolic flight, hypergravity, and vibration. A: The parabolic flight experimental rack; B: The Short Arm Human Centrifuge (SAHC) at the DLR, Cologne, Germany; C: Transportable incubator mounted on the SAHC; D: The Vibraplex device with T25 cell culture flasks mounted on it.

Cell fixation

The cells were fixed after the first parabola (1P) and after the 31st parabola (31P) using RNAlater (Applied Biosystems, Darmstadt, Germany) at a ratio of 4:1 (RNAlater:medium). After the flight, the fixative was discarded, the cells were briefly washed with PBS and covered with 10 ml of fresh RNAlater. Subsequently, the flasks were stored at 4°C and transported to the laboratory. For the quantitative real-time PCR, we collected n = 6 T75 cell culture flasks from both parabolic flight samples (1P and 31 P) and the 1 g control group, the remaining flasks were used for microarray analyses.

Hypergravity experiments

We performed experiments on the Short Arm Human Centrifuge (SAHC, DLR, Cologne, Germany) (Figure 6B), with cells cultured in T75 cell culture flasks (75 cm2; Sarstedt, Nümbrecht, Germany), growing in a monolayer. We installed two containers with floating mountings for the incubators on the SAHC (Figure 6C). In this configuration the T75 cell culture flasks were always exposed to a correct vertical gravity (acceleration) vector during centrifugation. By using the power supply on the SAHC the incubators were constantly heated to 37° C.

On the SAHC, we exposed the samples to a continuous hypergravity phase of 1.8 g of about 2 hours corresponding to the total time frame of 31 parabolas.

We designed a homogenous centrifuge profile with constant spin-up and spin-down times of each 34 seconds.

We collected n = 5 static 1 g controls and n = 5 1.8 g hyper-g samples for the Microarray analysis. The 1 g controls were cultivated in parallel in a neighboring identical incubator. Immediately after the run, the culture medium was discarded and replaced with 25 mL RNAlater solution. For the measurement of cytokines released in the supernatant, an additional run of the SAHC was performed to obtain n = 12 static 1 g samples and n = 12 1.8 hyper-g samples for the ELISA technique.

Vibration experiments

The detailed method was published earlier [20]. In short, the Vibraplex vibration platform (frequency range 0.2 Hz - 14 kHz) was used to create vibrations comparable to those occurring during parabolic flights (Figure 6D). Corresponding vibrations to the three phases of pull-up (1.8 g), free fall (μg), and pull-out (1.8 g) were recorded and analysed by Schmidt [45]. These data were then used for the simulation experiments with the Vibraplex. For quantitative real-time PCR analyses, we collected n = 5 samples of each of the two groups (1 g controls and cells subjected to a vibration profile corresponding to 31 parabolas of a parabolic flight). The 1 g controls without vibration were grown separately in a similar incubator.

For the measurement of cytokines released in the supernatant, three additional vibration experiments were performed to obtain n = 12 static 1 g samples and n = 12 vibrated samples for the ELISA technique.

RNA isolation and cDNA synthesis

After arrival in the laboratory, the RNAlater solution on the fixed cells was replaced by PBS (Invitrogen, Darmstadt, Germany). The cells were scraped off using cell scrapers (Sarstedt, Nümbrecht, Germany), transferred to 50 ml tubes, and pelleted by centrifugation (2500 g for 10 min at 4°C). An RNeasy Mini Kit (Qiagen, Hilden, Germany) was used according to the manufacturer’s instructions to isolate total RNA. RNA concentrations and quality were determined spectrophotometrically at 260 nm using an Ultrospec 2100 pro Spectrophotometer (Amersham Biosciences, Freiburg, Germany). The isolated RNA had an A260/280 ratio of >1.7. cDNA designated for the quantitative real-time PCR was then obtained with the First-Strand cDNA Synthesis Kit (Fermentas, St. Leon-Rot, Germany) using 1 μg of total RNA in a 20-μL reverse transcription reaction mixture.

Quantitative real-time PCR

Quantitative real-time PCR was used to determine the expression levels of the genes of interest. Primer Express software (Applied Biosystems) was applied to design appropriate primers with a Tm of ~60°C (Additional file 2). The primers were synthesized by TIB Molbiol (Berlin, Germany). All assays were run on a StepOnePlus Real-Time PCR System using the Power SYBR Green PCR Master Mix (both Applied Biosystems). The reaction volume was 25 μL including 1 μL of template cDNA and a final primer concentration of 500 nM. PCR conditions were as follows: 10 min at 95°C, 40 cycles of 30 s at 95°C and 1 min at 60°C, followed by a melting curve analysis step (temperature gradient from 60 to 95°C with +0.3°C/cycle).

If all amplicons showed one single Tm similar to the one predicted by Primer Express software, the PCR reactions were considered specific. Every sample was measured in triplicate, and relative quantification was effected by means of the comparative CT (ΔΔCT) method. 18S rRNA was used as a housekeeping gene to normalize the expression data.

ELISA

ELISAs of IL-6, IL-8, EGF, VEGFD (R&D Systems), and FGF17 (USCN Life Science Inc.) in the cell culture supernatant from vibration and hypergravity experiments have been performed according to the protocols supplied by the manufacturer.

Microarray analysis

Prior to the analysis, RNA integrity (RIN) was checked with the bioanalyzer. Only samples meeting the required quality were included in the analysis. The Illumina HumanWG-6_V2_0_R3 arrays have been normalized using the BeadStudio Gene Expression Module v3.3.7 and quantile normalization without background correction. After quantile normalization and exclusion of low or not expressed genes (minimum Illumina detection p-value > 0.05; performed in both analyses separately) the quality of arrays and the general expression profile has been checked by Principal Component Analysis (PCA) using Partek Genomic Suite 6.6, correlation as a dispersion matrix and normalized Eigenvector scaling. No obvious batch effect or outlier was found for the hypegravity analysis, while the outlying general expression profile of one sample from the parabolic flight experiment was removed before test statistic (for an overview: see Table 4).

Table 4 Microarray analyses: projects, samples and conditions

A parametric ANOVA comparing the conditions given in Table 5 was performed. The selection criteria for the significant differential expression are also given in Table 2. Differentiation of the expression profiles was performed using K-Mean clustering. The cluster analysis was done using Partek Genomic Suite 6.3 applying the Euclidean distance function on standardized log2 signal values. K was selected according to a local minimum of the Davies Bouldin K estimation procedure. Functional aspects of the differentially expressed probes were analyzed with g: Profiler using g: SCS threshold as significance criterion and the DAVID Bioinformatics Resources 6 [46,47]. Physical and functional interactions between proteins were determined using the String platform [48] at a low confidence score of 0.15.

Table 5 Microaray analyses: comparisons

Statistical analysis

All statistical analyses were performed using SPSS 16.0 software (SPSS, Inc., Chicago, IL, USA). We used either 1-way ANOVA or the Mann–Whitney U test. Differences were considered significant at the level of P < 0.05. All data are presented as means ± SD.

Abbreviations

BIRC3:

Baculoviral IAP repeat-containing protein3

BMP:

Bone morphogenetic protein

BP:

Biological processes

CASP:

Caspase

CAV:

Caveolin

CGM:

Chondrocyte growth medium

COL2A1:

Collagen type 2 α 1

CTGF:

Connective tissue growth factor

ECM:

Extracellular matrix

EDN1:

Endothelin 1

EGF:

Epidermal growth factor

ELISA:

Enzyme-linked immunosorbent assay

FCS:

Fetal calf serum

FDR:

False discovery rate

FGF:

Fibroblast growth factor

FKBP1A:

Peptidyl-prolyl cis-trans isomerase

FTC:

Follicular thyroid cancer

GO:

Gene Ontology

HAPLN1:

Hyaluronan and proteoglycan link protein 1

HMGB1:

High-mobility group protein B1

IL:

Interleukin

ITBG1:

Integrin beta 1

LAMA5:

Laminin α 5

μg :

Microgravity

MCL1:

Induced myeloid leukemia cell differentiation protein

MMP:

Matrix metalloproteinase

mRNA:

Messenger ribonucleic acid

NOTCH2:

Neurogenic locus notch homolog protein 2

NR4A2:

Nuclear receptor subfamily 4, group A, member 2

P:

Parabola

PCA:

Principal component analysis

PCR:

Polymerase chain reaction

PGF:

Placenta growth factor

PRKAA1:

Protein kinase, AMP-activated, α 1

PRKCA:

Protein kinase c α

PS:

Parabola set

SAHC:

Short arm human centrifuge

SD:

Standard deviation

SOCS3:

Suppressor of cytocine signaling 3

TBS:

Treated parabola set

TNF:

Tumor necrosis factor

TRAF1:

TNF receptor-associated factor 1

RPM:

Random positioning machine

TNFAIP3:

TNF α-induced protein 3

TNFRSF14:

TNF receptor superfamily member 14

VEGF:

Vascular endothelial growth factor

VIL2:

Ezrin

WNT5A:

Wingless-type MMTV integration site family, member 5A.

References

  1. Poole AR, Kojima T, Yasuda T, Mwale F, Kobayashi M, Laverty S. Composition and structure of articular cartilage: a template for tissue repair. Clin Orthop Relat Res. 2001;391(Suppl):S26–33.

    Article  PubMed  Google Scholar 

  2. Aigner T, Stöve J. Collagens–major component of the physiological cartilage matrix, major target of cartilage degeneration, major tool in cartilage repair. Adv Drug Deliv Rev. 2003;55:1569–93.

    Article  CAS  PubMed  Google Scholar 

  3. Wiberg C, Klatt AR, Wagener R, Paulsson M, Bateman JF, Heinegård D, et al. Complexes of matrilin-1 and biglycan or decorin connect collagen VI microfibrils to both collagen II and aggrecan. J Biol Chem. 2003;278:37698–704.

    Article  CAS  PubMed  Google Scholar 

  4. Tetlow LC, Adlam DJ, Woolley DE. Matrix metalloproteinase and proinflammatory cytokine production by chondrocytes of human osteoarthritic cartilage: associations with degenerative changes. Arthritis Rheum. 2001;44:585–94.

    Article  CAS  PubMed  Google Scholar 

  5. Clements KM, Flannelly JK, Tart J, Brockbank SM, Wardale J, Freeth J, et al. Matrix metalloproteinase 17 is necessary for cartilage aggrecan degradation in an inflammatory environment. Ann Rheum Dis. 2011;70:683–9.

    Article  CAS  PubMed  Google Scholar 

  6. Loeser RF. Aging and osteoarthritis: the role of chondrocyte senescence and aging changes in the cartilage matrix. Osteoarthr Cart. 2009;17:971–9.

    Article  CAS  Google Scholar 

  7. Young AA, Smith MM, Smith SM, Cake MA, Ghosh P, Read RA, et al. Regional assessment of articular cartilage gene expression and small proteoglycan metabolism in an animal model of osteoarthritis. Arthritis Res Ther. 2005;7:R852–61.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Alsalameh S, Amin R, Gemba T, Lotz M. Identification of mesenchymal progenitor cells in normal and osteoarthritic human articular cartilage. Arthritis Rheum. 2004;50:1522–32.

    Article  PubMed  Google Scholar 

  9. Zayzafoon M, Meyers VE, McDonald JM. Microgravity: the immune response and bone. Immunol Rev. 2005;208:267–80.

    Article  CAS  PubMed  Google Scholar 

  10. Greenwald E, Armstrong K, Leon P, Reinberger K. Growing cartilage in space-Understanding the physiological changes from microgravity. [http://chen2820.pbworks.com/w/page/11951455/Growing%20cartilage%20in%20space].

  11. Stamenković V, Keller G, Nesic D, Cogoli A, Grogan SP. Neocartilage formation in 1 g, simulated, and microgravity environments: implications for tissue engineering. Tissue Eng A. 2010;16:1729–36.

    Article  Google Scholar 

  12. Pietsch J, Ma X, Wehland M, Aleshcheva G, Schwarzwälder A, Segerer J, et al. Spheroid formation of human thyroid cancer cells in an automated culturing system during the Shenzhou-8 Space mission. Biomaterials. 2013;34:7694–705.

    Article  CAS  PubMed  Google Scholar 

  13. Grosse J, Wehland M, Pietsch J, Schulz H, Saar K, Hübner N, et al. Gravity-sensitive signaling drives 3-dimensional formation of multicellular thyroid cancer spheroids. FASEB J. 2012;26:5124–40.

    Article  CAS  PubMed  Google Scholar 

  14. Ulbrich C, Wehland M, Pietsch J, Aleshcheva G, Wise P, van Loon J, et al. The impact of simulated and real microgravity on bone cells and mesenchymal stem cells. Biomed Res Int. 2014;2014:928507.

    Article  PubMed Central  PubMed  Google Scholar 

  15. Ma X, Wehland M, Schulz H, Saar K, Hübner N, Infanger M, et al. Genomic approach to identify factors that drive the formation of three-dimensional structures by EAhy926 endothelial cells. PLoS One. 2013;8:e64402.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Ulbrich C, Westphal K, Pietsch J, Winkler HD, Leder A, Bauer J, et al. Characterization of human chondrocytes exposed to simulated microgravity. Cell Physiol Biochem. 2010;25:551–60.

    Article  CAS  PubMed  Google Scholar 

  17. Grimm D, Wehland M, Pietsch J, Aleshcheva G, Wise P, van Loon J, et al. Growing tissues in real and simulated microgravity: new methods for tissue Engineering. Tissue Eng B Rev. 2014;20:555–66.

    Article  Google Scholar 

  18. Herranz R, Anken R, Boonstra J, Braun M, Christianen PC, de Geest M, et al. Ground-based facilities for simulation of microgravity: organism-specific recommend-dations for their use, and recommended terminology. Astrobiology. 2013;13:1–17.

    Article  PubMed Central  PubMed  Google Scholar 

  19. Aleshcheva G, Wehland M, Sahana J, Bauer J, Corydon TJ, Hemmersbach R, et al. Moderate alterations of the cytoskeleton in human chondrocytes after short-term microgravity produced by parabolic flight maneuvers could be prevented by upregulation of BMP-2 and SOX-9. FASEB J.2015. Feb 13. pii: fj.14-268151. [Epub ahead of print].

  20. Wehland M, Ma X, Braun M, Hauslage J, Hemmersbach R, Bauer J, et al. The impact of altered gravity and vibration on endothelial cells during a parabolic flight. Cell Physiol Biochem. 2013;31:432–51.

    Article  CAS  PubMed  Google Scholar 

  21. Beckmann R, Houben A, Tohidnezhad M, Kweider N, Fragoulis A, Wruck CJ, et al. Mechanical forces induce changes in VEGF and VEGFR-1/sFlt-1 expression in human chondrocytes. Int J Mol Sci. 2014;15:15456–74.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Ferrara N. Role of vascular endothelial growth factor in regulation of physiological angiogenesis. Am J Physiol Cell Physiol. 2001;280:C1358–66.

    CAS  PubMed  Google Scholar 

  23. Ferrara N, Gerber HP, LeCouter J. The biology of VEGF and its receptors. Nat Med. 2003;9:669–76.

    Article  CAS  PubMed  Google Scholar 

  24. Pufe T, Lemke A, Kurz B, Petersen W, Tillmann B, Grodzinsky AJ, et al. Mechanical overload induces VEGF in cartilage discs via hypoxia-inducible factor. Am J Pathol. 2004;164:185–92.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Pufe T, Kurz B, Petersen W, Varoga D, Mentlein R, Kulow S, et al. The influence of biomechanical parameters on the expression of VEGF and endostatin in the bone and joint system. Ann Anat. 2005;187:461–72.

    Article  CAS  PubMed  Google Scholar 

  26. Steinert A, Weber M, Dimmler A, Julius C, Schütze N, Nöth U, et al. Chondrogenic differentiation of mesenchymal progenitor cells encapsulated in ultrahigh-viscosity alginate. J Orthop Res. 2003;21:1090–7.

    Article  CAS  PubMed  Google Scholar 

  27. Reddi AH. Bone morphogenetic proteins: from basic science to clinical applications. J Bone Joint Surg Am. 2001;83A Suppl 1:S1–6.

    Google Scholar 

  28. Yuan W, Zhao MD, Yuan FL, Che W, Duan PG, Liu Y, et al. Association of endothelin-1 expression and cartilaginous endplate degeneration in humans. PLoS One. 2013;8:e60062.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  29. Saleem M, Qadir MI, Perveen N, Ahmad B, Saleem U, Irshad T, et al. Inhibitors of apoptotic proteins: new targets for anticancer therapy. Chem Biol Drug Des. 2013;82:243–51.

    Article  CAS  PubMed  Google Scholar 

  30. Hosseini-Farahabadi S, Geetha-Loganathan P, Fu K, Nimmagadda S, Yang HJ, Richman JM. Dual functions for WNT5A during cartilage development and in disease. Matrix Biol. 2013;32:252–64.

    Article  CAS  PubMed  Google Scholar 

  31. Lane Smith R, Trindade MC, Ikenoue T, Mohtai M, Das P, Carter DR, et al. Effects of shear stress on articular chondrocyte metabolism. Biorheology. 2000;37:95–107.

    CAS  PubMed  Google Scholar 

  32. Smith RL, Carter DR, Schurman DJ. Pressure and shear differentially alter human articular chondrocyte metabolism: a review. Clin Orthop Relat Res. 2004;427(Suppl):S89–95.

    PubMed  Google Scholar 

  33. Zhang Y, Chee A, Shi P, Wang R, Moss I, Chen EY, et al. Allogeneic articular chondrocyte transplantation downregulates interleukin 8 gene expression in the degenerating rabbit intervertebral disk in vivo. Am J Phys Med Rehabil. 2014. Aug 14. [Epub ahead of print].

  34. Ryu JH, Yang S, Shin Y, Rhee J, Chun CH, Chun JS. Interleukin-6 plays an essential role in hypoxia-inducible factor 2α-induced experimental osteoarthritic cartilage destruction in mice. Arthritis Rheum. 2011;63:2732–43.

    Article  CAS  PubMed  Google Scholar 

  35. Civinini R, Nistri L, Martini C, Redl B, Ristori G, Innocenti M. Growth factors in the treatment of early osteoarthritis. Clin Cases Miner Bone Metab. 2013;10:26–9.

    PubMed Central  PubMed  Google Scholar 

  36. Krejci P, Krakow D, Mekikian PB, Wilcox WR. Fibroblast growth factors 1, 2, 17, and 19 are the predominant FGF ligands expressed in human fetal growth plate cartilage. Pediatr Res. 2007;61:267–72.

    Article  CAS  PubMed  Google Scholar 

  37. Blanchard JM. Cyclin A2 transcriptional regulation: modulation of cell cycle control at the G1/S transition by peripheral cues. Biochem Pharmacol. 2000;60:1179–84.

    Article  CAS  PubMed  Google Scholar 

  38. Abdollahi T, Robertson NM, Abdollahi A, Litwack G. Inhibition of TRAIL-induced apoptosis by IL-8 is mediated by the p38-MAPK pathway in OVCAR3 cells. Apoptosis. 2005;10:1383–93.

    Article  CAS  PubMed  Google Scholar 

  39. Aleshcheva G, Sahana J, Ma X, Hauslage J, Hemmersbach R, Egli M, et al. Changes in morphology, gene expression and protein content in chondrocytes cultured on a random positioning machine. PLoS One. 2013;8:e79057.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. Duke PJ, Daane EL, Montufar-Solis D. Studies of chondrogenesis in rotating systems. J Cell Biochem. 1993;51:274–82.

    Article  CAS  PubMed  Google Scholar 

  41. Chavagnac C, Gai F, Gharib T, Mora C. Astrium spaceplane for scientific missions. Acta Astronaut. 2013;92:172–7.

    Article  Google Scholar 

  42. Grosse J, Wehland M, Pietsch J, Ma X, Ulbrich C, Schulz H, et al. Short-term weightlessness produced by parabolic flight maneuvers altered gene expression patterns in human endothelial cells. FASEB J. 2012;26:639–55.

    Article  CAS  PubMed  Google Scholar 

  43. Ma X, Pietsch J, Wehland M, Schulz H, Saar K, Hübner N, et al. Differential gene expression profile and altered cytokine secretion of thyroid cancer cells in space. FASEB J. 2014;28:813–35.

    Article  CAS  PubMed  Google Scholar 

  44. Ulbrich C, Pietsch J, Grosse J, Wehland M, Schulz H, Saar K, et al. Differential gene regulation under altered gravity conditions in follicular thyroid cancer cells: relationship between the extracellular matrix and the cytoskeleton. Cell Physiol Biochem. 2011;28:185–98.

    Article  CAS  PubMed  Google Scholar 

  45. Schmidt W. Quickly changing acceleration forces (QCAFs) vibration analysis on the A300 ZERO-G. Microgravity Sci Technol. 2004;15:42–8.

    Article  PubMed  Google Scholar 

  46. Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID Bioinformatics Resources. Nat Protoc. 2009;4:44–57.

    Article  CAS  Google Scholar 

  47. Huang DW, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37:1–13.

    Article  PubMed Central  Google Scholar 

  48. Snel B, Lehmann G, Bork P, Huynen MA. STRING: a web-server to retrieve and display the repeatedly occurring neighbourhood of a gene. Nucleic Acids Res. 2000;28:3442–4.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The study was supported by the German Aerospace Center (DLR; BMWi grant 50WB1124; D.G.). Ganna Aleshcheva and Elisabeth Warnke are doctoral candidates of the Helmholtz Space Life Sciences Research School, German Aerospace Center Cologne, Germany.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Daniela Grimm.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

DG, MI and RH designed the study. MW, EW, MB, JP, XM, GA, TJC, and SR attended the parabolic flight mission and performed the flight experimens. MW and DG drafted the manuscript. MW, GA, TJC and DG supervised and supported the ground-based experiments. DG, RH and TF supervised the ground-based experiments in Cologne and supported the SAHC and Vibraplex experiments. KS, NH and HS performed the Gene Array analysis and the bioinformatics. XM performed the qPCR. All authors read and approved the final manuscript.

Markus Wehland and Ganna Aleshcheva contributed equally to this work.

Additional files

Additional file 1:

This file provides expanded versions of Tables 1 and 2 as well as the complete overview of the GO BP analysis of the microaray data.

Additional file 2:

Primers used for quantitative real-time PCR. All sequences are given in 5′-3′ direction.

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wehland, M., Aleshcheva, G., Schulz, H. et al. Differential gene expression of human chondrocytes cultured under short-term altered gravity conditions during parabolic flight maneuvers. Cell Commun Signal 13, 18 (2015). https://doi.org/10.1186/s12964-015-0095-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12964-015-0095-9

Keywords