Microgravity influences circadian clock oscillation in human keratinocytes

Highlights • Changes of gravitational forces affect oscillations of circadian clock genes.• The effects of microgravity on the circadian clock persist during the recovery process from microgravity.• Linking of mechanotransduction to circadian rhythms reveals changes in physiology.


Introduction
Mechanical stimuli are crucial to many biological processes that take place in cells, tissue and organs. Almost any cell is able to respond to mechanical cues that are registered by specific sensor molecules and commuted into biochemical signals that regulate cell shape and morphology, but also participate to the execution of numerous programs, including cell proliferation, differentiation and migration. The most common and typical forces applied to living cells in vivo are shear flow, tensile stretch, compressive forces, and a variety of mechanical stimuli sometimes different for magnitude, frequency and duration.
The molecular sensors responsible to trace forces applied to living cells, are numerous and include extracellular-matrix proteins, receptors and transmembrane proteins, components of the cytoskeleton, voltage or ligand sensitive ion channels [1][2][3]. These structures, organized to contain or function themselves as ''mech anoreceptors" respond to mechanical stress by generating signals that become finely integrated with the signaling network with the intent to override, synergize or simply tune biological responses.
In this context, gravity and sudden variations of gravitational forces have significant impact on living cells, with reflection on tissue and organ physiology. This is evident in astronauts that, exposed to micro-gravitational status for a long period of time, suffer of loss bone, immunosuppression [4], cardiovascular problems [5], changes in metabolism and body temperature activity [6,7], sleep defects [8,9], blood pressure alteration, hormone and neurotransmitter impaired secretion [10]. Some of these effects (e.g. loss bone or heart rate changes) can be variably ascribed to the altered equilibrium of signaling events, in cells anymore exposed to compressional mechanical forces [11], whereas others might be related to circadian desynchronization with the environment [7].
As a yet unexplored hypothesis, we wanted to verify whether microgravity would impact on circadian gene regulation at individual cell compartment. To test this, we exposed HaCaT human keratinocytes to microgravity in the Random Positioning Abbreviations: FAK, focal adhesion kinases; GPCR, G-protein-coupled receptors; hg, hypergravity; mAbs, monoclonal antibodies; RPM, Random Positioning Machine Machine (RPM) device for different periods of time, monitoring variation in the expression of the core clock genes Bmal1 and Rev-erba. Selection of this cellular model was based on previous studies demonstrating the presence in human keratinocytes and melanoma cells of functional clock genes subjected to circadian oscillation for gene expression [12][13][14]. Using this model we recently described the expression and nucleolar localization of a novel splicing variant of Period 2 gene (Per2S) [15], suggesting a critical role for the nucleolus in the regulation of circadian rhythms. On the other hand, keratinocytes are also known to respond to mechanical stretches by increasing, in particular conditions, proliferation and migration rates [16].
In a simple experimental setting we found that microgravity amplifies Bmal1/Rev-erba-sustained circadian oscillations, indicating the existence of integrated gene networks connecting the machinery of mechanotransduction to the cellular timing system of the circadian clock.

Cells and treatments
Human spontaneously transformed keratinocyte HaCaT cells were cultured in Dulbecco's Modified Eagles Medium (DMEM) supplemented with 10% Fetal Bovine Serum (FBS) plus antibiotics. Cells were seeded in Opticell units (Biocrystal Ltd, Westerville, OH, USA) and allowed to attach on both the two inner membrane surfaces of the chamber by repetitive turning upside down every 15 min during 4 h. All the experiments were performed after 24 h from cell plating.
Microgravity conditions were obtained by positioning the Opticell chambers containing HaCaT cells in the Random Positioning Machine (RPM) device (Dutch Space, Leiden, Netherlands) inside a humidified incubator (5% CO 2 at 37°C), setting the angular velocity of rotation at 90°/sec as maximum and 30°/sec as minimum, in a random mode [17]. Under this experimental condition, the cells were exposed to simulated microgravity conditions ranging from 0 to 0.01g, and experiments performed for different periods of time (48 h or 60 h at 37°C).

RNA extraction and cDNA synthesis
Total RNAs were extracted from HaCaT cells using the Quick-RNA MiniPrep (Zymo Research, Irvine, CA, USA) and concentrations were determined using NaNoDrop spectrophotometer (Thermo scientific, Huntsville, TX, USA). All the cDNAs were synthesized by the reverse transcription of 1 lg of total RNA using the iScript cDNA synthesis kit (Bio-Rad Laboratories, Richmond, CA, USA). Reactions were performed according to manufacturer's instructions.
For each primer pair, we performed no-template control and no-reverse-transcriptase control (RT negative) assays, which produced negligible signals. Real-time PCR were carried out in 96-well plate using iQ SYBR Green Supermix (Bio-Rad) and the iCycler Real-Time Detection System (iQ5; Bio-Rad). The relative quantification of gene expression was assessed by comparative threshold cycle method (2-DDCT ), normalizing the gene expression to ribosomal 18S RNA and expressing the data as fold increase. All the experiments have been done in triplicates and are reported as mean ± SD.

Immunofluorescence
HaCaT cells exposed to gravitational forces and pharmacological synchronization of circadian genes were fixed with 4% paraformaldehyde in PBS for 30 min at 25°C. Fixed cells were permeabilized using 0.1 M glycine for 20 min at 25°C followed by additional 5 min of incubation with 0.1% Triton X-100 at 25°C. Per- System (Zeiss, Oberkochen, Germany) connected with an Axiovert 200 inverted microscope (Zeiss) for the analysis of fluorescent signals; all the images were processed using the Axiovision software (Zeiss). The three dimensional reconstruction shown in each figure is a selection of three of the total number of serial optical sections evaluated under the microscope. Quantitative analysis of the percentage of Ki67 positive cells was assessed evaluating 400 cells from each sample, randomly observed in 5 optical fields from three different experiments. Results have been expressed as mean values ± standard errors (SE). p value were calculated using one-way analysis of variance -ANOVA -followed by Bonferroni post hoc test, and significant level has been defined as p < 0.05.

Simulated microgravity modifies Bmal1 and Rev-erba circadian oscillations
The effects of microgravity can be monitored analyzing the oscillations of the molecular clock genes taken at different endpoint times in cells previously synchronized by pharmacological treatment. In our experiments, we treated HaCaT keratinocytes with dexamethasone 1 lM for 1 h. This treatment allows for the synchronization of the cell clock machinery [12] and precedes the exposure of cells to the variation of gravitational forces. We choose four different conditions for the analysis: exposure of unsynchronized and Dex-synchronized cells to normal gravitational force (1g), and exposure of unsynchronized and Dex-synchronized cells to micro-g environment (lg) (Fig. 1A). In this experimental setting we evaluated the oscillations of two circadian clock genes: Bmal1 and its repressor Rev-erba.
Quantitative analysis for mRNA expression of Bmal1 in Dex-treated synchronized cells exposed to normal gravity indicates the existence of an oscillatory period of 24 h, with peak of expression in the first 6 h, and return to the basal levels within the 24 h. Importantly, changes in the mechanical forces applied to Dex-synchronized cells subjected to microgravity had a strong impact on the Bmal1 kinetics. Its expression peaks in the first 6 h (3.9-fold increase), remains constantly high at 12th hour (3.2-fold increase), tends to normalize only at the approach of the 24th hour (1.1-fold increase) and then recycles maintaining the same of oscillatory phase, amplitude and period kinetics for the next 24 h. Dex-treated cells exposed to microgravity display a significant enhancement of amplitude in comparison to Dex-synchronized cells kept in 1g (Fig. 1B), (p < 0.05, one-way analysis of variance -ANOVA -followed by Bonferroni post hoc test). As expected, the expression of Bmal1 in unsynchronized cells exposed to either 1g or microgravity did not show any oscillatory period during the 48 h of observation.
The rhythmic expression of Bmal1 is coordinately controlled by the transcriptional repressor Rev-erba directly binding to the ROR elements in the Bmal1 promoter. This tightly regulated mechanism accounts for an anti-phasic expression of Rev-erba to Bmal1. Therefore, we chose to evaluate the expression of the clock gene Rev-erba mRNA to asses the specificity of Bmal1 oscillations observed in cells exposed to microgravity. As expected, to the highest peaks of Rev-erba gene expression corresponds a minimal synthesis of Bmal1 mRNA, which occurs at the onset of the 24th hour. However, the amplified oscillations of Rev-erba observed in cells exposed to microgravity were not significantly different from Dex-synchronized cells kept a 1g, indicating that Rev-erba may not be subjected to variation of gravitational forces. Alternatively, exposure to microgravity produces limited but still efficient lM Dex for gene clock synchronization were exposed to gravitational forces of 1g or lg. Total RNA samples were collected every 6 h for consecutive 48 h and subjected to qRT-PCR analysis for Bmal1 and Rev-erba mRNA expression. (B) Bmal1 mRNA transcripts show no oscillation in unsynchronized HaCaT cells exposed to either lg or 1g. Both Dex-synchronized cells cultures at 1g (Dex + 1g) or in lg (Dex + lg) show Bmal1 circadian oscillations, with higher Bmal1 transcription in cells exposed to lg (p < 0.05, one-way ANOVA). Rev-erba mRNA show no oscillation in unsynchronized HaCaT cells exposed to either lg or 1g. Both Dex-synchronized cells cultures at 1g (Dex + 1g) or in lg (Dex + lg) show minimal increase of Rev-erba circadian oscillations, with higher transcripts in cells exposed to lg. GAPDH mRNA expression is not subjected to circadian oscillations. All the experiments were done in triplicates and were evaluated for statistical significance (*p < 0.05) using the one-way ANOVA followed by Bonferroni post hoc test. Results are expressed as mean fold increase ± standard deviation (SD). (C-F) Unsynchronized and Dex-synchronized cells exposed to microgravity for 24 h were stained fluctuation of Rev-erba proteins that might be sufficient to negatively regulate transcription of Bmal1. As an internal control for all the experimental conditions, we monitored glyceraldehyde-3phosphate dehydrogenase (GAPDH) mRNA expression, which is not subjected to circadian rhythm oscillations (Fig. 1B).
Since the exposure of cells to decreased gravitational forces is known to induce alterations of the cytoskeleton dynamics [20,21] with changes of the adhesion properties, proliferation rate and apoptosis [17,[22][23][24], we wanted to verify these potential effects on cells exposed for 24 h to microgravity. We choose cytokeratin and b1 integrin as markers for the evaluation of cytoskeleton dynamics, Ki67 to detect cycling cells, and the cleaved form of cytokeratin 18 for the evaluation of the early events of apoptosis. We could not appreciate any morphological change or variation in the organization of the cytokeratin filaments and b1 integrin adhesion molecules in cells cultured in microgravity environment (Fig. 1C and D). Also, comparative evaluation of proliferative and apoptotic rates in cells either exposed to microgravity or normal gravity did not show differences in the experimental condition tested ( Fig. 1E and F, Table 1).
These results demonstrate that microgravitational forces applied to cells specifically amplify the synchronized molecular clock in the absence of evident physiological changes of human keratinocytes and that the experimental conditions utilized for our experiments do not represent a major source of cell stress.

The cycle periods and amplitudes of Bmal1 and Rev-erba are enhanced during the recovery phase
It has been recently shown that lg influences gene expression in human keratinocytes and that these effects are evident for a period time after recovery to 1g [25]. We reasoned that the oscillations of the clock genes expression Rev-erba and Bmal1 would also be present after the exposure to microgravity, for a period of time during recovery at 1g. To verify this hypothesis, HaCaT cells were exposed to lg or 1g for consecutive 60 h and only then synchronized with Dex in 1g conditions to evaluate the expression of Rev-erba and Bmal1 measured by qRT-PCR at 6 h, 12 h, 24 h, 30 h, 36 h, 42 h and 48 h post synchronization ( Fig. 2A).
We evaluated the circadian oscillations of Rev-erba and Bmal1 during the recovery phase from microgravity for consecutive 48 h. The results indicate that Rev-erba and Bmal1 oscillations are variably affected for intensity of expression, time of peaking and period length during recovery. In particular, Bmal1 oscillations in Dex-synchronized cells previously exposed to microgravity were significantly and constantly higher as compared to cells continuously cultured at 1g conditions (Fig. 2B). Furthermore, we noticed a significant anticipation of Bmal1 expression with peak at 12 h from synchronization and a lengthened cycle period, more pronounced in Dex-treated cells in recovery (p < 0.05, one-way ANOVA followed by Bonferroni post hoc test). By contrast, the analysis of Bmal1 expression revealed no oscillations in asynchronous cells pre-cultured at 1g or lg conditions. The parallel evaluation of Rev-erba indicated a general anti-phasic expression as compared to Bmal1, with nadir at the 12 h from synchronization corresponding to the zenith of Bmal1, with a period of about 30 h for synchronized cells, either cultured at 1g or in the recovery phase and with amplitudes and peaks of the same intensity. Again, GAPDH mRNA did not show any circadian oscillations (Fig. 2B).
Cell integrity and functional behaviour were anyhow assessed at the end of the 60 h of culture in lg, through the evaluation of cytokeratin organization, integrin dynamics, number of apoptotic or proliferative cells (Fig. 2C and F, Table 1).

The effects of microgravity on clock genes synchronization do not depend on the pharmacologic treatment
Genes of the peripheral clock can be synchronized by numerous pharmacologic treatments, each one interfering with identified signaling pathway [26]. To verify the specificity of our treatment for synchronization of the clock machinery and to ascertain the absence of interference with other pathways, we duplicate the experiment of microgravity preconditioning exposing HaCaT cells to 1 h of 1 lM forskolin after 60 h of lg. Bmal1 mRNA was quantified by qRT-PCR at 12 h, corresponding to the zenith of expression (Fig. 3A). As expected, the two synchronization protocols (Dex and Forsk) showed overlapping results with peaks of Bmal1 expression equally high in Dex-and Forsk-treated cells (Fig. 3B and C) compared to cells continuously cultured at 1g conditions.

Discussion
Numerous and substantiated evidence indicates that light is the primary environmental cue that entrains circadian clocks. However, gravity and change of gravitational forces affect circadian patterns with significant changes of general behavior, hormone synthesis, body temperature and also metabolism [6][7][8]10,[27][28][29]. These observations highlight the importance of the biological responses that cells display to counteract mechanical cues and, among these, changes of gravitational forces applied to cells.
It has been previously shown that rats exposed to hypergravity (hg) from mid-gestation though early lactation manifest severe alterations of multiple organs regulating metabolic activity of mammary, liver and adipose tissues with dramatic effects on maternal behavior and pup survival. These changes of the homeorhetic responses leading also to an unsuccessful lactation, are sustained by defective hormones incretion and alteration of circadian gene expression [30]. However, there is no direct evidence that microgravity and sudden changes of gravitational forces may influence the circadian oscillations of clock genes in cells not subjected to hormonal or metabolic control.
Through mechanoreceptors, epithelial cells initiate bidirectional inside-outside mechanotransduction, coupling ECM molecules and cell-cell adhesion structures, e-cadherin/a-catenin, integrins-focal adhesion kinases (FAK) complexes and G-proteincoupled receptors (GPCR) to intracellular signaling pathways through the recruitment and activation of PKC, PI3K, MAP and Rho kinases. In skin, convergence of these signals seems to be mediated by mechanosensitive transcription factors such as b-catenin, YAP and AP-1 [1]. Skin is a suitable model to study these effects: in fact, circadian rhythm in HaCaT cells is sustained by high-amplitude oscillation of the clock genes, Per2, Per3, Bmal1 and Rev-Erba [12]. However, increased amplitudes of circadian genes expression by microgravity have not been reported before. We report that microgravitational forces applied to HaCaT cells specifically amplify the molecular clock gene Bmal1, whose amplitude, time of peaking and period, parallel the expression of the Rev-erba clock genes. We also describe amplitude changes and significant lengthened of Bmal1 and Rev-erba cycle periods in the recovery phase.
In fact, quantitative analysis of Bmal1 mRNA in Dex-treated synchronized HaCaT cells exposed to microgravity indicates a general increase of amplitudes with peak of a 3.9-fold increase at 6 h from synchronization. Furthermore, we observed similar effects on Bmal1 oscillation in cells exposed to lg for 60 h, then suddenly exposed to normal gravitational force. In this case, Bmal1 mRNA expression is 5.5-fold higher after 6 h from synchronization as compared to HaCaT cells that, maintained at the inertial exposure of 1g constantly show a sinusoidal Bmal1 oscillation of 24 h. By contrast, an increased length period and higher peaks of Bmal1 expression characterize the recovery phase. As expected, to the highest peaks of Rev-erba gene expression correspond a minimal synthesis of Bmal1 mRNA, indicating a general maintenance of the transcriptional control of Rev-erba over Bmal1 along the circadian cycle. Interestingly, 6 h is the time constantly observed and necessary for HaCaT cells to produce the maximal Bmal1 peak of expression in the transition from 1g to microgravity, while 12 h is the time necessary to observe Bmal1 peak in cells readjusting to 1g. However, while the mRNA expression levels of Rev-erba faithfully recapitulate the anti-phasic expression to Bmal1, lack of fluctuation of Rev-erba seems to indicate the existence of more complex mechanisms that govern their expression in cells exposed to variation of gravitational forces, and the molecular links sustaining the elongation of cycle periods are presently unknown. These effects can be explained in light of the existence of mechanotransduction machinery able to translate variation of mechanical forces applied to cells by gravitational changes into transcriptional regulated clock gene expression.
Taken together, our results demonstrate that microgravity and sudden changes of gravitational status specifically influence the molecular clock, without affecting the physiology of human keratinocytes, at least in the window of observation of 48 h. Importantly, the effects of lg on clock oscillation are detectable also during the recovery period in 1g conditions. In this case, enhanced Bmal1 transcription would be the consequence of additional mechanical stress on cells shifting from microgravity to 1g. This ''double mechanical hit" in cells first exposed to microgravity, then back to 1g, would not only explain persistence of high Bmal1 amplitudes during recovery period, but even justify its further increase found in our experimental setting (3.9 in microgravity; 5.5-fold from microgravity to 1g), with a cycle period that appears cultured either in 1g (1g + Dex) or in microgravity (lg + Dex). The kinetics of Bmal1 mRNA transcripts indicates that cells exposed to lg display atypical oscillation recovery as compared to cells recovered from 1g, different for intensity, period and time of peaking (*p < 0.05, one-way ANOVA), with a 5.5-fold increase of Bmal1 mRNA. Results are expressed as mean values ± standard deviation (SD). The Rev-erba kinetics in synchronized cells either cultured at 1g or in the recovery phase is anti-phasic to Bmal1, with a period of 30 h and with amplitudes and peaks of equal intensity. GAPDH expression is used as an internal control. All the experiments were done in triplicates and evaluated for statistically significance using the one-way ANOVA test. Results are expressed as mean fold increase ± standard deviation (SD). (C-F) Cells were stained for immunofluorescence after 60 h of lg or 1g. Immunofluorescence was performed using specific antibodies and nuclei visualized by DAPI staining. Bar: 10 lm. to be significantly lengthened (36 h period for Bmal1, 30 h for Reverba during the recovery phase). Translated to the experience of astronauts returning from space, these findings indicate the early phase of the recovery period as a critical moment for rehabilitation from microgravity, and would suggest the need of rapid and harmonic re-entrainment with our zeitgeber (see Table 1).
We did not observe appreciable morphological changes as measured by reorganization of cytokeratins and b1 integrins, rather than functional modification expressed by growth or cell death rates. This is not surprising, since only in a few cases are strength and amplitude of mechanical signals sufficient to induce adaptive changes in cell morphology, or modify the patterns of proliferation and differentiation. In general, a single type of mechanical load is unlikely to be sufficient to generate competent mechanotransductive signals to specify for defined phenotypes. Nevertheless, these signals are generated and utilized to refine transductional signals and implement transcription. To date, no mechanosensors have been identified that specifically transduce variation of gravitational forces and the pathways involved in the transcriptional control of clock genes. It will be important to pursue studies that will shed light on mechanisms of skin physiology and the interrelationships between circadian clock, mechanotransduction and cancer.