Rapid coupling between gravitational forces and the transcriptome in human myelomonocytic U937 cells

The gravitational force has been constant throughout Earth’s evolutionary history. Since the cell nucleus is subjected to permanent forces induced by Earth’s gravity, we addressed the question, if gene expression homeostasis is constantly shaped by the gravitational force on Earth. We therefore investigated the transcriptome in force-free conditions of microgravity, determined the time frame of initial gravitational force-transduction to the transcriptome and assessed the role of cation channels. We combined a parabolic flight experiment campaign with a suborbital ballistic rocket experiment employing the human myelomonocytic cell line U937 and analyzed the whole gene transcription by microarray, using rigorous controls for exclusion of effects not related to gravitational force and cross-validation through two fully independent research campaigns. Experiments with the wide range ion channel inhibitor SKF-96365 in combination with whole transcriptome analysis were conducted to study the functional role of ion channels in the transduction of gravitational forces at an integrative level. We detected profound alterations in the transcriptome already after 20 s of microgravity or hypergravity. In microgravity, 99.43% of all initially altered transcripts adapted after 5 min. In hypergravity, 98.93% of all initially altered transcripts adapted after 75 s. Only 2.4% of all microgravity-regulated transcripts were sensitive to the cation channel inhibitor SKF-96365. Inter-platform comparison of differentially regulated transcripts revealed 57 annotated gravity-sensitive transcripts. We assume that gravitational forces are rapidly and constantly transduced into the nucleus as omnipresent condition for nuclear and chromatin structure as well as homeostasis of gene expression.

subsequently gene expression [12][13][14][15][16] . Small forces in the low piconewton range may finally trigger nuclear mechanotransduction 17 and force transduction into the chromatin can occur within seconds 12 . Thus, the nucleus is subjected to permanent small direct or indirect 18 forces induced by Earth's gravity, raising the fundamental question, if gene expression homeostasis is constantly shaped by the gravitational force on Earth.
Only experiments in microgravity allow to investigate gene expression under force-free conditions and therefore facilitate the elucidation of the role of Earth's gravity in gene expression homeostasis, while time-resolved studies would help to assess the adaptation potential in an altered gravitational environment. Microarray-based gene expression studies have been conducted previously with T cells or T cell lines in simulated microgravity [19][20][21][22] , in spaceflight experiments [23][24][25] , and with lymphatic tissue from space-flown animals 26, 27 . As a result, altered expression of microRNA in simulated microgravity conditions 21 correlated with the gene expression pattern of the transcription factor Rel 21 , which has been identified as microgravity-dependent gene expression regulator in a spaceflight experiment 24 . Previous gene expression studies in altered gravity have been mostly end-point measurements after time periods of hours or longer in microgravity and focused on the identification of particular gravity-responsive genes. Since gene expression responds very rapidly to altered gravity within or less than minutes 28,29 and force transduction into chromatin requires only seconds 12 , initial mechanisms can be studied in the minute range, where different microgravity platforms (parabolic flights, suborbital ballistic rockets) are available for multi-platform analysis at an integrative level 4,29 .
Therefore, we recently investigated the dynamics of gene expression response to different gravitational environments in human Jurkat T lymphocytic cells during parabolic flight and suborbital rocket experiments 4,29 , identified gravity-regulated genes, but also revealed an overall high stability of gene expression in microgravity 4 . Experiments with cells of the immune system not only address fundamental biological questions about the effects of gravity on cellular homeostasis, but also the important medical risk of exploration class long-term manned space missions requiring mitigation 30 . Thus, the immune system belongs to the most affected systems during spaceflight (reviewed in) [31][32][33] and sensitivity of cells of the human immune system to reduced gravity has been confirmed in numerous studies in real and simulated microgravity in T lymphocytes and cells of the monocyte-macrophage-system (MMS) [33][34][35][36][37][38][39] , but also indicated the existence of fast cellular adaptation 38 . In this study, we therefore focused on the first and initial transcriptome events in cells of the MMS. Due to the operational constraints of the conducted experiment missions, we used U937 human myelomonocytic cells, as established during other microgravity and space experiments 36,37 .
In addition to the hypothesis of direct force transduction into chromatin, force-sensitive ion channels have been discussed as trigger point of mechanotransduction 7,8,40,41 into complex cellular reactions such as gene expression. Macrophages harbor ATP-gated P2X channels, store-operated Orai channels and members of the transient receptor potential (TRP) cation (TRPC) channel family, which play important roles in inflammation and phagocytosis [42][43][44][45][46][47][48] . Because TRPC1 is activated by stretch [49][50][51] , it represents a candidate for transduction of gravitational forces in cells of the MMS. This hypothesis was corroborated by the finding that in the unicellular photosynthetic flagellate Euglena gracilis, knockdown of a putative TRP channel abolished gravitaxis 52 . In order to test a potential role of ion channels and in particular TRPC channels, we also conducted inhibitor experiments with SKF-96365, commonly used to characterize the potential functions of TRPC channels and blocking voltage-activated calcium channels 53,54 , sarco-and endoplasmic reticula Ca 2+ pumps 55,56 , voltage-gated sodium currents 57,58 , and ATP-sensitive and voltage-gated potassium channels 58,59 . Therefore, the combination of a wide range of ion channel inhibition with whole transcriptome analysis offers the possibility to study the functional role of ion channels in the transduction of gravitation forces at an integrative level.
Aim of this study was to investigate the transcriptome in force-free conditions of microgravity and in hypergravity, using different research platforms, gravity conditions and time points in order to identify potential molecular candidates of gravitational force-responsive gene regulation and to test the role of ion channels in gravitational force-dependent gene regulation. Our approach allowed the identification and validation of gravity-regulated gene expression through two fully independent large-scale research campaigns. Therefore, transcriptional changes identified after both campaigns are characterized by a high level of evidence due to independent sets of experiments in combination with independent research platforms.

Results
During the parabolic flight campaign (19 th DLR PFC), U937 cells were subjected to 20 s of hypergravity (1.8 g) and subsequently to 20 s of microgravity during the first parabola, and samples were obtained at the end of each flight phase. Control samples were obtained in-flight 5 min before the first parabola and on ground ( Table 1, Fig. 1). In case of the TEXUS-49 suborbital rocket flight, samples were acquired 75 s after lift-off after the hypergravity launch phase and before the microgravity phase and after 5 min of microgravity flight. Control samples were prepared on ground. Further sample sets in microgravity and on ground were obtained in presence of 25 µM SKF-96365, an ion channel inhibitor of TRPC channels, voltage-gated Ca 2+ channels and potassium channels. For both campaigns, RNA from at least four samples in each group was isolated, labeled and hybridized on a microarray chip.
Rapid and extensive transcriptome alterations after 20 s of altered gravity. During the first parabola in the parabolic flight experiment, the transcriptome of human U937 cells responded rapidly within 20 s of hypergravity with a total number of 17998 (17228 annotated) (BL-PFC hyp-g versus H/W 1 g GC) differentially regulated transcripts. After the subsequent first microgravity phase, 11810 (11361 annotated) transcripts (µg versus H/W 1 g GC) were differentially expressed, which is 2.7-fold the number of 4293 (4165 annotated) transcripts that were altered solely as a result of the flight conditions without altered gravity (comparison 1 g IF versus H/W 1 g GC, After exclusion of all transcripts already altered due to flight conditions (comparison 1 g IF versus H/W 1 g GC), we revealed 10345 hypergravity-sensitive annotated transcripts. In the microgravity-phase (µg versus BL-PFC hyp-g) we identified 1715 microgravity-sensitive annotated transcripts (Table 2, Table 3, and Fig. 2). Six transcripts were differentially expressed in the same direction in the microgravity phase as in the hypergravity phase and were therefore eliminated in order to exclude potential protracted effects, resulting in 1709 differentially regulated transcripts (Table 3 and Fig. 2). The expression fold change (FC) for the hypergravity-sensitive transcripts was in the range between +4.68 (up-regulation) and −7.8 (down-regulation) with average values of +1.97 and −1.79, respectively. The expression fold change (FC) for the microgravity-sensitive transcripts was in the range between +3.35 (up-regulation) and −2.66 (down-regulation) with average values of +1.52 and −1.58, respectively (Table 3).

Extensive transcriptome alterations after 75 s of hypergravity and 300 s of microgravity.
During the launch phase of the suborbital ballistic rocket experiment, the transcriptome of human U937 cells responded rapidly within 75 s covering the hypergravity phase with a total number of 10849 (10556 annotated) differentially regulated transcripts (comparison BL-TX hyp-g versus H/W 1 g GC, Table 4, Fig. 3). After the subsequent 300 s microgravity phase, 4783 (4668 annotated) transcripts were differentially expressed (µg vs BL-TX hyp-g, Table 4), which is 2.3-fold lower than the number of differentially expressed transcripts in hypergravity. 1291 (1247 annotated) transcripts were differentially expressed in microgravity when compared to the H/W 1 g ground control group (Table 4). One single transcript was differentially expressed in the same direction in hypergravity (BL-TX hyp-g vs H/W 1 g GC) and microgravity (µg vs BL-TX hyp-g) and was therefore eliminated in order to exclude potential protracted effects, resulting in 4667 microgravity-sensitive baseline (BL)-controlled annotated transcripts (Table 5 and Fig. 3a). The expression fold change (FC) was in the range between +8.85 (up-regulation) and -3.38 (down-regulation) with average values of +2.11 and −1.68 respectively for the hypergravity-sensitive transcripts and between +7.23 (up-regulation) and −3.26 (down-regulation) with average values of +1.85 and −1.71 respectively for microgravity-sensitive transcripts ( Table 5).
Role of SKF-96365-sensitive ion channels in gravity-regulated gene expression. Because forcesensitive ion channels have been discussed to induce transduction of mechanical forces into complex cellular reactions 7,8,40,41 , we used the wide-range ion channel inhibitor SKF-96365 to identify ion channel-dependent regulated gene expression in altered gravity. The number of differentially expressed transcripts with and without incubation with SKF-96365 are shown in Table 6 (comparisons H/W 1 g GC SKF versus H/W 1 g GC and µg SKF versus H/W 1 g GC SKF). SKF-96365 had distinct effects on gene expression: 4920 (4732 annotated) transcripts were differentially regulated in 1 g. In microgravity, SKF-96365 more than doubled the number of differentially regulated genes to 10959 (10413 annotated) transcripts (Table 6). Differentially regulated transcripts were now further verified by the appropriate control experiments (Table 7) and only transcripts were selected, which were (1) microgravity-sensitive (baseline controlled), (2) differentially expressed in the comparison µg SKF versus µg either in the same direction (enhancement of microgravity effect) or in the opposite direction (reversal of microgravity effect), and (3) that were not differentially expressed in 1 g (comparison H/W 1 g GC SKF-96365 versus H/W 1 g GC). In total, 13 annotated transcripts showed an enhanced sensitivity to microgravity and 99 annotated transcripts showed a decreased sensitivity to microgravity in the presence of SKF-96365 (Table 7, Supplementary  Table 1, and Fig. 3b). The expression fold change (FC) was in the range between +1.8 (up-regulation) and −1.43 (down-regulation) with average values of +1.5 and −1.43 respectively for annotated transcripts with an enhanced microgravity sensitivity and between +2.04 (up-regulation) and −4.17 (down-regulation) with average values of +1.5 and −1.89 respectively for annotated transcripts with an attenuated microgravity sensitivity (Table 7). In conclusion, we could attribute only 118 (112 annotated) out of 4782 (4667 annotated) microgravity-regulated transcripts, equivalent to 2.4%, to gravity-induced functional effects of SKF-96365-sensitive ion channels.
Gene ontology enrichment analyses. We performed Gene Ontology (GO) enrichment analyses for an integrative evaluation of the effects of altered gravity on cellular processes and functions. 20 s of microgravity influenced catabolic processes and cation channel transcripts, whereas 20 s of hypergravity had an impact on DNA replication and microtubule-based processes (Supplementary Figs S1 and S2). 75 s of hypergravity and 300 s of microgravity generated a higher number of affected cellular processes and functions with partially overlapping aspects ( Supplementary Figs S3 and S4). 300 s of microgravity influenced transcripts for intracellular transport, RNA and enzyme binding, mRNA processing, posttranscriptional regulation of gene expression, cell cycle and cell division ( Supplementary Fig. S3). 75 s spanning hypergravity impacts transcripts of metabolic processes, regulation of transcription, intracellular transport, cell cycle and cell division ( Supplementary Fig. S4).
No strong effects on apoptosis/necroptosis-associated pathways in altered gravity. In order to investigate proposed effects of altered gravity on apoptosis 60,61 , we screened our data sets for genes involved in     (Fig. 4a), longer periods of hypergravity resulted in primarily up-regulation of transcripts and longer periods of microgravity in down-regulation of transcripts (Fig. 4b).
In the case of the ion channel inhibitor SKF-96365 enhancing the effect of microgravity, transcripts were primarily up-regulated, in the case of SKF-96365 reversed the effect of microgravity, transcripts were almost equally distributed between up-and down-regulation (Fig. 4c).
Evidence for rapid counter-regulation of initially altered transcripts. We then analyzed all transcripts from the parabolic flight experiments that were differentially expressed in both micro-and hypergravity, and identified 1602 (1532 annotated) transcripts (Fig. 5). Surprisingly, all identified transcripts were regulated in the opposite direction in hyper-and microgravity: Transcripts that were up-regulated in hypergravity, were down-regulated in microgravity and vice versa. None of the identified transcripts was regulated in the same direction in hypergravity and in microgravity (Fig. 5). We continued this analysis with the data from the suborbital ballistic rocket experiment and identified 3574 (3514 annotated) transcripts ( Fig. 6) that were differentially expressed in hypergravity and in microgravity and revealed the same observation: All transcripts, except for one, were expressed also in the opposite direction in hyper-and microgravity, respectively (Fig. 6). Due to the nature of the flight profiles, the microgravity phase always followed the hypergravity phase in both flight experiments.   We therefore investigated the hypothesis that reversely regulated transcripts are caused by counter-regulatory mechanisms. For this reason, we performed a cross-platform comparison of the hypergravity and microgravity phases, separately and with regard to their different duration. We identified 2192 (2150 annotated) transcripts that were altered in both platforms after 20 s and 75 s of hypergravity, respectively. We found that 2079 (2039 annotated) of these transcripts were regulated in the opposite direction, while 113 (111 annotated) transcripts   were differentially expressed in the same direction (Fig. 7). The cross-platform comparison of the microgravity phases revealed 110 (106 annotated) transcripts that were altered in both platforms after 20 s and 300 s of microgravity, respectively. We found that 98 (97 annotated) transcripts were regulated in the opposite direction, while 12 (9 annotated) transcripts were differentially expressed in the same direction ( Fig. 8). That means that more than 92% of microgravity-sensitive annotated transcripts that were differentially expressed on both platforms were regulated in opposite directions and are therefore adapting to the new gravitational environment by counter-regulation.

Identification of gravity-regulated genes.
Inter-platform comparisons of all transcripts which were differently regulated in hyper-and microgravity revealed 58 (57 annotated) transcripts. All of these transcripts were reversely regulated between the two gravity conditions within one platform. Additionally, all transcripts but one were reversely regulated in the same gravity condition comparing both platforms (Fig. 9, Table 8). Interestingly, one of these 57 transcripts was ATP6V1E1, representing a vacuolar H + -ATPase (V-ATPase), already identified previously as a gravity regulated gene in human Jurkat T cells 29 .
Almost complete adaptation of initially differentially altered transcripts. The time course of differentially regulated transcripts in hypergravity is summarized in Fig. 10 and in microgravity in Fig. 11  altered transcripts adapted after 5 min. In hypergravity, 98.93% of all initially altered transcripts adapted after 75 s. Therefore, we identified two pools of altered transcripts: A first one, which responded after seconds and adapted within 5 min latest, and a second one, which appeared at least after 5 min of altered gravity.

Discussion
Fast transcriptome response to altered gravity and potential gravity-sensitive cellular structures. One of the first theoretical studies about the physical background of cellular microgravity concluded that cells with a diameter of 10 μm and more would experience gravity 63 . However, because the weight of single normal-sized cells of 10 µm is too small compared with other cellular forces for direct "sensing" of the gravity vector 18 , gravity may thus be "sensed" indirectly at least as the magnitude of gravity-induced weight forces surrounding tissues and fluids 18 . Since decades it has been asked for mechanism by which the gravitational force becomes transmitted into a biological process 64 , but was rarely addressed in experimental approaches. In a very recent review about perception of gravity in eukaryotes 65 , gravity perception by cytoskeletal processes and mechanosensitive ion channels have been discussed, but the presented theories about gravitational force transmission in mammalian cells has not been validated experimentally so far. Therefore, we used a multi-platform approach to determine the time period between gravitational force change and transcriptome response, and also investigated the role of ion channels in gravitational force-induced transcriptome regulation. We were able to detect profound alterations in the transcriptome as early as 20 s after the onset of altered gravitational force. Therefore, the preceding hypothetical transduction processes have to be fast enough for inducing transcriptome changes within 20 s of force alteration. We therefore compiled the transduction velocities of potential transduction processes: DNA decondensation was shown to occur within the time frame of minutes to hours 66 , transcription factor binding takes a couple of seconds 67 , and the elongation rate of RNA polymerase II was found to range between 0.37 kb/ min 68 24 . In summary, mechanosensation at the plasma membrane leads to downstream nucleocytoplasmic shuttling of various transcription regulators 17 . Thus, force transduction along activation of transcription factors is very well described and a potential reason for fast transcriptome response to altered gravity, but requires 10 min or more into the nucleus 86,87 , probably not fast enough to explain rapid transcriptome alterations after 20 s. In this context, experiments with the TRPC channel inhibitor SKF-96365 revealed that only a very low share of all detected gravity-induced transcriptome alterations can be attributed to this type of ion channels.
Transcriptional or posttranscriptional regulation processes after rapid mechanical transduction. Fast transcriptome alterations could be initiated also at the posttranscriptional level: RNA binding proteins (RBPs) bind mature mRNA and exert important regulatory effects on mRNA stability and translation. RBPs can either function to promote mRNA degradation as in the case of AUF1 or TTP shortening the mRNA half-life, or they promote mRNA stability like HuR extending the mRNA half-life. Furthermore, RBPs have been shown to work cooperatively with miRNAs to regulate mRNA turnover. For example, HuR is able to recruit let-7/RISC to inhibit c-MYC mRNA translation 88 and it competes with miR-494 and miR-548c-3p for the binding of nucleolin and TOP2A mRNA, respectively 89,90 . RBPs were found to be associated to the cytoskeleton 91,92 and play critical roles in the localization of mRNA 93 . Apart from that, it has been observed that binding of the extracellular matrix to integrins on the cell surface and mechanical tension induce movement of mRNA to focal adhesions 94 . These findings suggest that mechanical stimuli act on RBP function and thus on mRNA regulation. mRNAs that are translationally repressed, have been found to accumulate in cellular structures called P-or GW-bodies where they are degraded 95,96 . Alternatively, mRNA can be stored as a response to stress conditions in cytoplasmic aggregates named stress granules (SGs) 97 . In P-bodies proteins are enriched that play an important role in RNA deadenylation, decapping and degradation. Interestingly, we identified the 5′-3′ exoribonuclease xrn1, representing a key player in the cellular mRNA degradation process, being differentially regulated when exposed to altered gravity. Furthermore, several other genes involved in RNA splicing, transcriptional regulation and RNA processing were found highly sensitive to gravitational changes. miRNAs, which were recently found differentially expressed after 20 s of microgravity in Jurkat T cells 29 , appear to repress translation and promote decay by recruiting P-body components to individual mRNAs 95 . SGs and P-bodies share protein components and can be found in close vicinity in  99 . The cellular pause of protein synthesis triggered by environmental stress conditions (heat, UV irradiation, osmotic pressure etc.) is important for the cell to focus on repair processes, e.g. due to DNA damage. The translational stress that follows the environmental stress is highly specific for different types of genes. While some selected groups of RNAs are stabilized, others are destabilized and degraded 100 . In the hypothetical model presented by Anderson and Kedersha 97,101 it is shown that the sorting of mRNAs into SGs or P-bodies is reversible, meaning that the remaining RNAs can return to the polysome fraction after the cell reaches the equilibrium again. This process of re-distribution occurs rapidly and shuttling times in and out of SGs have been measured with a half-life between 2 and 8 seconds 98 . Therefore, significant SG-bodies-caused transcriptome alterations are theoretically possible after 20 s, if the preceding signal cascade to the cytoplasm requires only seconds.
In conclusion, while biochemical force-transduction processes appear unlikely to explain the rapid transcriptome response, all hypothetical mechanisms such as transcriptional regulation, RBPs or P-body or SG-based regulations are only fast enough if followed by rapid mechanical transduction either into the nucleus or into the cytoplasm lasting in the range of microseconds 71 . To explain how cells can not only sense but also spatially transfer physical forces over relevant distances, the model of "tensegrity" architecture of the cytoskeleton implies that cytoskeletal components are interconnected with a certain level of tensional pre-stress, established by a force balance between extracellular adhesions, contractile microfilaments and microtubules 6,8,102 . This architecture not only provides shape stability to the cell, but it also implicates that tensional forces are propagated across the whole tensegrity structure which includes intracellular structures [102][103][104] . Propagation of stress along cytoskeletal filaments takes only 2 µs to cover 50 µm, whereas diffusion-based biochemical signaling through the cytosol needs seconds to span the same distance 71 .

Rapid force transduction chains to the chromatin via the cytoskeletal-nucleoskeletal network.
Although the majority of mechanotransduction research has focused on the perception of mechanical forces at and across the cell membrane to induce signaling pathways originating in the cytoplasm 17,71 , many experimental evidences demonstrated a very close connection between cytoskeletal force transduction and chromosome organization and gene expression: Alterations in cell geometry resulted in cytoskeletal reorganization, leading to nuclear morphology remodeling, affecting orientation, 3D radial position, compaction, and intermingling of chromosome territories 16 and chromatin condensation 105 , accompanied by differential gene expression patterns 106 . Force transmission through the actin cytoskeleton to the nuclear envelope resulted in direct displacements of Cajal body-associated protein complexes 107 and through the LINC complex to chromatin deformation and force-induced expression up-regulation of specific genes 14 . Therefore, mechanical forces can regulate gene expression independent of molecular relays by opening or closing chromatin configurations 72 . Mechanical force transduction into the chromatin seems to occur within seconds, subsequently resulting in transcription alteration proportional to the magnitude of chromatin stretching 12 and was observed to induce epigenetic changes in promoter regulation 108 . Nuclear lamins, emerin, LINC complex proteins, heterochromatin protein-1 (HP1) and BAF have been recently proposed to be involved in force transduction to the chromatin, resulting in stretching and increased accessibility of the transcriptional machinery 12 . BAF binds directly to linker histone H1.1 and core histone H3 in vitro and in vivo 109 and is associated with poly(ADP-ribose) polymerase-1 (PARP-1) 110 , a central modulator of chromatin structure and transcription 111 . Interestingly, PARP-1 has been demonstrated to be involved in hyperstretch-induced mechanotransduction in bronchial epithelial cells 112 . PARP-1 and PARP-2 also interact physically and functionally with HP1 113 . Importantly, HP1 binds to histone H3 methyl lysine 9 and SuvH39, a histone H3 lysine methyltransferase, where methylation of histone H3 lysine 9 by SUV39H1 creates a binding site for HP1 proteins 114 bridging H3K9me3 in condensed chromatin 115 . Because silencing H3K9 methyltransferase SUV39H1 completely blocked the force-induced H3K9 methylation 116 , and H3K9me3 has been shown to be mechanical strain-driven 13 , Suv39H and H3K9me3 are potential components of the gravity-induced force transmission cascade to the chromatin.
Surprisingly, more than 98% of all initially altered transcripts adapted rapidly to the new gravitational environment after 75 s of hypergravity or 300 s of microgravity, respectively. Thus, the same gravity condition obviously resulted in qualitatively different directions of transcript regulation as a function of time. Gene ontology enrichment analyses ( Supplementary Figs 1-4) revealed regulatory effects of altered gravity on RNA binding, regulation of transcription and posttranscriptional regulation of gene expression. In this context, the transcription factor AKNA, the histone-binding transcriptional repressor L3MBTL2, and the RNA binding protein LOC643980 belong to the group of non-adapting microgravity-sensitive transcripts (Supplementary Table S3). Although molecular mechanisms of rapid transcriptome adaptations to altered gravity are unknown, they appear probably at different steps of transcriptional processing. Since the macrophageal oxidative burst reaction adapted rapidly to microgravity in only 42 s 39 , we assume rapid adaptation mechanisms are occurring not only inside the nucleus, but also in the cell membrane. It could be therefore possible, that the cytoskeletal-nucleoskeletal network is reacting to gravitational force alterations as a whole interconnected system.
No pro-apoptotic/necroptotic effects in altered gravity. In order to exclude potential pro-apoptotic/ necroptotic effects during the entire cell preparation process for the experiment mission and to investigate potential pro-apoptotic effects of altered gravity as described earlier 60,61 , we analyzed apoptotic/necroptotic genes identified from the KEGG database and from relevant literature and revealed only slight transcript expression differences for them. We compared these values to previously published ones of apoptotic gene expression in human monocytes displaying fold changes of >+5 for caspase 6 117 , in human chronic myelogenous leukemia K562 cells with fold changes of +1.779 and +1.633 for caspase 6 and caspase 8 respectively after chemical stimulation 118 ,  119 . In comparison to this data we observed a rather mild cellular reaction in altered gravity (Supplementary Tables S3 and S5). Additionally, all microgravity-altered pro-apoptotic/necroptotic transcripts adapted after 5 min, as well as all hypergravity-altered transcripts, except of FADD (BC000334). However, FADD is down-regulated after 20 s and 75 s of hypergravity, indicating a rather protective effect of the altered gravity condition. In this context, transient up-regulation of apoptotic markers with a subsequent reversibility and cellular recovery has been described previously [120][121][122] .
The experiments of this study were performed in different very low gravity environments (10 −2 -10 −3 g for parabolic flight experiments and 10 −4 g for TEXUS experiments). Currently, it is unknown, in which extent different levels of very low gravity are transduced into a cellular response. Whereas a 2D clinostat study with 1F6 melanoma cells reported differences in guanylyl cyclase A mRNA expression in the range between 0.012-0.036 g 123 , the response of the oxidative burst reaction in NR8383 macrophages did not differ between the range of 10 −2 -10 −3 g for parabolic flight experiments 35 and the <10 −5 g for the ISS experiment 39 . Additionally, indications for a gravitational threshold between 0.3 g and 0.5 g were found 39 . Therefore, all microgravity levels are probably below this threshold, but current knowledge about biological effects of gravitational changes in very low gravity environments is limited. A recent study re-analyzed transcriptome data from two space experiments (ISS, human T cells and Progress 40 P, HUVEC) and two ground-based studies with "simulated" microgravity (humans T cells and PBMC) and identified the reduction of second messenger molecules generation as a statistically significant pathway in all datasets, which has been discussed as potential effect of microgravity-induced low shear stress 124 . Due to the rapid response and adaptation of the transcriptome (Figs 10 and 11), the described alterations after longer microgravity exposure time could re-represent a stable systemic regulatory status, a steady-state, and therefore the last and final phase of microgravity-induced transcriptome alterations. Of course, the rapid cellular response to altered gravity 35,39 , detected in cells of the monocyte-macrophage system in this study, might not be present in other cell types, particularly if caused by cell-type specific chromosome territory architectures 16,105,106 .  In conclusion, we detected an immediate transcriptome response after 20 seconds of altered gravity, which in turn adapted rapidly afterwards. Theoretical explanations for rapid transcriptome regulation require rapid transduction into the nucleus, which is likewise theoretically possible and has been demonstrated for other mechanical forces than gravity. Thus, we assume that gravitational forces are rapidly and constantly transduced into the nucleus, acting as an omnipresent condition in the nuclear and chromatin structure and subsequently resulting in homeostasis of gene expression. Due to the proposed fundamental role of gravitational forces for transcriptome homeostasis, the existence of rapid adaptation mechanisms is not surprising. Gravity could potentially represent an essential environmental condition for the genomic functionality of life on Earth.

Material and Methods
Cell culture. The human myelomonocytic cell line U937 (ATCC CRL1593.2) was used as a model cell line to analyze the differential gene expression under altered gravity conditions in the human monocyte-macrophage system. U937 cells were cultured in RPMI-1640 (Biochrom/Merck Millipore, Germany), supplemented with 10% fetal bovine serum (FBS Superior; Biochrom/Merck Millipore, Germany), 1% glutamine (200 mM; Gibco/Life Technologies, Germany) and 1% penicillin/streptomycin (10,000 U/ml and 10,000 µg/ml respectively; Gibco/Life Technologies, Germany). Cells were cultured with a density of 0.2 × 10 6 cells/ml and medium exchange was performed every 48 hours. Cells were centrifuged at 300 g for 5 min at room temperature, the supernatant was discarded, and the cell pellet was resuspended in fresh medium. An aliquot was taken, diluted with trypan blue solution and the number of vital cells was counted. Cells were reseeded at a concentration of 0.2 × 10 6 cells/ml in fresh medium. Parabolic flight experiment platform. As described previously 29 , parabolic flights are an ideal platform to study initial and primary effects in mammalian cells and the associated rapid responsive molecular alterations excluding influences and interferences of secondary signal cascades. Parabolic flights offer a sequence of consecutive gravity conditions including 1 g, 1.8 g, and microgravity (µg) with a quality of 10 −2 to 10 −3 g. We designed and constructed an experimental system which allows cell culture experiments during parabolic flights on board the Airbus A300 ZERO-G (reg. no. F-BUAD), that has been used already for different parabolic flight experiments 28,29,36 . Primary importance was placed on realizing the direct safety technique during the development activity. The experimental structure (Fig. 1c) consists of three experiment racks (storage rack for cell culture containers before the experiments at 36.5 °C, cooling rack for storage of cell culture containers after cell lysis at 4 °C, and a working rack for handling and execution of the experiments). The modular system is able to accommodate up to 54 cell culture containers (double containment) for each flight and allows storage of cell cultures until the start of the experiment, injection of a fluid (culture medium) at any defined time during the parabolic maneuver, and automatic injection of a second fluid (lysis buffer) after 20 s at the end of a defined gravity phase. Appropriate in-flight controls were obtained during the 1 g flight phase directly before the first parabola. Injection of all fluids operates automatically and is pre-programmed, while exchange of cell culture containers and supervision of the experiment was performed manually. During the 19 th DLR parabolic flight campaign (PFC), we investigated the gene expression in human U937 cells in microgravity and hypergravity (1.8 g) compared to in-flight 1 g. Experiments were only conducted during the first parabola to assure that detected differential gene expression levels were a result of the effect of gravitational change and not an accumulated long-term effect.
Preparation and execution of the parabolic flight experiments. During the 19 th DLR PFC, 1 × 10 7 U937 cells in 10 ml medium (RPMI 1640 supplemented with 1% penicillin/streptomycin, amphotericin (Gibco/ Life Technologies, Germany), 1% glutamine and 2% FBS (i.e. serum starved) were filled into 200 ml Nutrimix bags (B. Braun Melsungen, Germany) and transported from the home laboratory to the pre-flight preparation laboratories at the NOVESPACE premises in Bordeaux, France. After arrival, cells were de-starved with 0.8 ml FBS per sample and stored at 36.5 °C overnight and used for the flight experiment on the following morning. 36.5 °C were chosen instead of 37 °C to rule out any thermic activation of the cells caused by regulatory oscillation of the storage rack. For the flight day, the Nutrimix bags were placed in a solid plastic housing to create a double containment that prevents spillage of fluids in the aircraft in case of leakage of the hardware system. Rapid lysis of U937 cells in the respective gravity phase was achieved by fast injection of 5 volumes of RLT buffer (Qiagen, Germany) and mixing by inverting the samples three times immediately at the appropriate time point (1 g in-flight samples 5 min before the first parabola, 1.8 g and microgravity samples during the first parabola). After landing, 1 g ground controls were performed immediately using the same hardware inside the aircraft. Post-flight,  all samples were directly transported to the on-site laboratory where total RNA was purified. In total, 28 samples were obtained from two parabolic flight days: 6x H/W 1 g GC, 8x 1 g IF, 6x BL-PFC hyp-g, 8x µg (see Table 1).
RNA isolation after the parabolic flight. RNA was isolated as described previously 29 . After landing of the aircraft and transport of the samples to the laboratory facilities, the protective plastic housings were disassembled, Figure 10. Time course of differential gene expression in hypergravity. Human myelomonocytic U937 cells were exposed to 20 s and 75 s of hypergravity. The numbers of annotated transcripts grouped according to their regulation after the two exposure times are shown. Continuous response means that annotated transcript are either up-or down-regulated at both time points. Adaption is either disappearance of the hypergravity induced effect or regulation in the opposite direction after 75 s. Late response is a regulation that only appears after 75 s. Most transcripts did not respond to hypergravity at both time points.
the Nutrimix bags were gently agitated and the lysed cell solution was filled into a T75 straight neck cell culture flask. The cell solution was mixed for 10 s by vortexing and sheared by passing four times through a Ø 0.8 × 120 mm needle (B. Braun Melsungen, Germany) fitted to a sterile 50 ml syringe. 50 ml of absolute ethanol were added and precipitates were resuspended by vigorous shaking. A Qiavac 24 plus vacuum system (Qiagen, Figure 11. Time course of differential gene expression in microgravity. Human myelomonocytic U937 cells were exposed to 20 s and 300 s of microgravity. The numbers of annotated transcripts grouped according to their regulation after the two exposure times are shown. Continuous response means that annotated transcript are either up-or down-regulated at both time points. Adaption is either disappearance of the microgravity induced effect or regulation in the opposite direction after 300 s. Late response is a regulation that only appears after 300 s. Most transcripts did not respond to microgravity at both time points. Germany) was prepared by placing 24 valves and sterile connective pieces on the Qiavac 24 plus vacuum manifold and an RNA maxi column (Qiagen, Germany) was attached to each connective piece. The system was set to a vacuum level of -200 mbar and the RNA maxi columns were loaded with the lysed cell suspensions. Subsequently, the valves were closed and the RNA maxi columns were centrifuged at 3220 g for 3 min at room temperature. The, 15 ml of buffer RW1 (Qiagen, Germany) were carefully applied to the column to wash the membrane-bound RNA. After centrifugation at 3220 g for 7 min at room temperature, the flow through was discarded and additional two washing steps were performed with 10 ml RPE buffer (Qiagen, Germany) followed by centrifugation at 3220 g for 3 min and 10 min at room temperature, respectively. The column-bound RNA was eluted by application of 600 µl of pre-warmed RNase-free water (Qiagen, Germany), incubation for 1 min at room temperature and centrifugation for 4 min at 3220 g again at room temperature. The elution step was repeated with the first eluate, the column was centrifuged for 7 min at 3220 g, and the purified RNA was stored in a sterile 1 ml cryotube on dry ice. Finally, the extracted RNA was transported on dry ice and stored at −80 °C until the processing of the RNA for the microarray analysis.   experimental unit was immediately recovered and returned to the launch site within 1.5 h after lift-off by helicopter. The general experimental composition consists of multiple sets of three syringes, filled with cell suspension (human U937 cells), cell culture medium with or without SKF-96365, and lysis solution (Trizol LS). All three syringes were connected by a T-piece, while small plugs at the outlet ports prevented premature contact of the fluids. The syringe systems were housed in a temperature-controlled, vacuum-resistant container (Fig. 1d). The temperature-controlled syringe systems were placed at microgravity positions inside the payload structure. Before launch and during flight, syringes were activated by a pneumatic system at pre-set time points. Several pre-flight tests and development tests were conducted: Biocompatibility tests, chemical stability tests, culture medium optimization with regard to buffer systems and supplements, sterilization tests, viability tests, cell lysis tests (different lysis compounds and concentrations). The entire mission procedure was standardized and tested several times. Margins and possible holding times were determined. The experimental setup consisted of the baseline group (lysis after hypergravity phase and before onset of microgravity), in-flight microgravity group (lysis after 5 min of microgravity and before reentry into the Earth's atmosphere), and 1 g ground control reference inside the experimental hardware. Cells, medium and lysis fluid (Trizol LS) syringes were prepared directly before the launch. All procedures started 7 hours before launch. The experimental containers were integrated into the payload structure by a "late access" port between 1:15 h and 0:45 h before launch. Sample temperature was maintained at 36.5 °C ± 0.5 °C until lysis. On landing and payload recovery, the experimental containers were immediately removed and returned to the ESRANGE laboratory for further processing. The cell suspension was transferred from the syringes into sterile plastic reaction tubes, and cells were homogenized with subsequent isolation of RNA. The purified RNA was stored and transported on dry ice or in liquid nitrogen and analyzed afterwards by means of genome-wide expression arrays.

Experimental preparation and integration for TEXUS-49. U937 cells were cultured in the ESRANGE
laboratories on site. Cells were cultivated with a density of 0.2 × 10 6 cells/ml, and the medium was exchanged every 48 hours (see above). During the countdown phase, cells were visually inspected, harvested, the vital cell number was counted, and cells were pooled to a concentration of 5 × 10 7 cells/ml. 0.5 ml of cells (i.e. 25 million cells) were filled in sterile 3 ml plastic syringes shortly before the handover to the launch team. Additionally, a second set of syringes was filled with 0.3 ml of cell culture medium with or without 0.3 ml SKF-96365 (25 µM; Sigma Aldrich, Germany), and a third set with 1 ml Trizol LS (Life Technologies, Germany) per sample unit. The three syringes with small plugs at the outlet ports were mounted on a sterilized plastic T-block with a connecting tubing system. 26 experiment units were prepared and stored at 36.5 °C ± 0.5 °C until integration into the payload of the rocket or until manual execution of the ground controls, respectively. These experimental units were finally integrated into the automatically operated experiment system. During the experimental run, firstly 0.3 ml of cell culture medium with or without SKF96365 and secondly 1 ml of Trizol LS were injected to the cell suspension at defined time points to lyse the cells and preserve the current status of differential gene expression. Directly before the µg phase, a set of samples was lysed at the time point of 75 s after launch (baseline, BL), representing the effect of hypergravity, spin and vibrations during the launch and rocket engine burn. Two further sets of samples (with and without SKF-96365) were fixed at 375 s after launch, shortly before the end of the µg phase. Additionally, 1 g ground controls were kept on ground in the incubator analogously to the µg sample group. In total, 26 samples were obtained after the TEXUS-49 rocket flight: 6x H/W 1 g GC, 5x BL-TX hyp-g, 7x µg, 4x µg SKF, 4x H/W 1 g GC SKF (see Table 1).

RNA isolation after TEXUS-49 landing.
The sample processing has been described previously 29 . Directly after landing, localization and recovery of the payload by helicopter, the experiment modules were dismantled and handed over for processing. The sample containing syringes were connected to a sterile 20G needle (B. Braun Melsungen, Germany), the 1.8 ml of cell suspension were sheared three times and distributed equally in two 2.0 ml reaction tubes. 0.1 ml of chloroform (Sigma-Aldrich, Germany) were added, the homogenate was vortexed for 15 s and incubated for 5 min at room temperature before a 15 min centrifugation step at 11000 g and 4 °C. The upper phase from both 2.0 ml tubes was transferred into one 15 ml tube and 4 ml of RLT buffer (Qiagen, Germany), as well as 3 ml of absolute ethanol were added and the suspension was mixed. 4 ml of this solution were pipetted on an RNA Midi column (Qiagen, Germany) and centrifuged for 30 s at 3000 g and room temperature. The flow through was discarded and the residual 4 ml of RNA solution were loaded on the column. All samples were centrifuged for 5 min at 3000 g at room temperature. Then, the columns were washed twice with 2.5 ml of RPE buffer (Qiagen, Germany) and centrifuged firstly for 2 min and additionally for 5 min at 3000 g at room temperature. The RNA was eluted by the addition of pre-warmed 250 µl RNase-free water (Qiagen, Germany) to the column, incubation for 1 min at room temperature, and centrifugation for 3 min at 3000 g and room temperature. The flow through was loaded again onto the column, incubated for 1 min at room temperature, and centrifuged for 5 min at 3000 g and room temperature. The isolated RNA was transferred into sterile 1 ml cryotubes and stored and transported at −80 °C. After arrival in the home laboratory, samples were stored at −80 °C until the processing of the RNA for the microarray analysis.
RNA sample processing and microarray data analysis. The RNA quantity and quality of the samples of the 19 th DLR parabolic flight campaign and the TEXUS-49 sounding rocket mission was analyzed using a Nanodrop 1000 (Thermo Scientific). All RNA samples were of high quality with 260/280 nm ratios between 1.9 and 2.1. The RNA integrity number (RIN) was measured using an Agilent 2100 Bioanalyzer (Agilent Technologies, USA) and only RNA samples with a RIN of >8.7 were used for the following microarray analysis. 400 ng total RNA were cy3-labeled with the Low RNA Input Linear Amplification Kit, PLUS, One-Color" (Agilent Technologies) and hybridized for 17.5 h to a NimbleGen expression microarray (12 × 135,000 features) SCIEntIfIC REPORTS | (2018) 8:13267 | DOI:10.1038/s41598-018-31596-y applying the "Gene Expression Hybridization Kit" (Agilent Technologies, USA). Microarrays were washed and scanned with the Micro Array Scanner G2505B (Agilent Technologies, USA). The image files of the scanner were analyzed with the NimbleScan Software 2.6 using the Robust Multi-Array Analysis (RMA) with default parameters. RMA represents a probe-level summarization method that identifies probes that are outliers in the overall behavior of the expression measured for a given gene. Differential expression of transcripts was determined based on the normalized microarray data using Excel 2013 and expression fold changes (FCs) of all transcripts on the microarray were calculated. FCs were used for comparisons of experimental groups within one experimental platform. For this, averages of the linear expression values were determined for each experimental group. The ratio was calculated by dividing the average value of one experimental group by the average value of the experimental group to which it should be compared. If the ratio is ≥1 it is equal to the FC, if the ratio is <1 the FC was determined by building the negative reciprocal of the ratio. Furthermore, t-tests were performed for all comparisons. FCs ≤ −1.3 or ≥1.3 with p-value < 0.05 were regarded to represent a significantly differential expression.

Intra-platform and inter-platform comparisons.
To further validate the relation between the different gravitational conditions and the detected differential expressions, intersections were made between the experimental comparison and respective control-comparison whenever available. Transcripts which were differentially expressed in the experimental and in the control-comparison were consequently excluded from the pools of gravisensitive transcripts (Tables 3, 5, 7, 8 and Venn diagrams). To analyze if transcripts respond preferentially to hyper-or to microgravity or to both conditions, intersections were made between the pools of hypergravity-and microgravity-sensitive transcripts. Furthermore, to evaluate the gravisensitivity of transcripts over time, intersections were made between the pools of gravisensitive transcripts of both experimental platforms.

Gene annotation enrichment analysis (GOrilla analysis). The gene annotation enrichment analysis
was carried out by using DAVID 6.8. Enriched GO terms visualized in ranked lists of genes were generated using GOrilla 125 . For illustration purposes only end points of processes and functions are summarized. The p-value threshold was set to <0.05.

Data Availability
The datasets generated and analyzed during the current study can be accessed in the GEO (Gene Expression Omnibus) repository (www.ncbi.nlm.nih.gov/projects/geo) under accession no. GSE101309.