Meta-Analysis-Assisted Detection of Gravity-Sensitive Genes in Human Vascular Endothelial Cells

Gravity affects the function and maintenance of organs, such as bones, muscles, and the heart. Several studies have used DNA microarrays to identify genes with altered expressions in response to gravity. However, it is technically challenging to combine the results from various microarray datasets because of their different data structures. We hypothesized that it is possible to identify common changes in gene expression from the DNA microarray datasets obtained under various conditions and methods. In this study, we grouped homologous genes to perform a meta-analysis of multiple vascular endothelial cell and skeletal muscle datasets. According to the t-distributed stochastic neighbor embedding (t-SNE) analysis, the changes in the gene expression pattern in vascular endothelial cells formed specific clusters. We also identified candidate genes in endothelial cells that responded to gravity. Further, we exposed human umbilical vein endothelial cells (HUVEC) to simulated microgravity (SMG) using a clinostat and measured the expression levels of the candidate genes. Gene expression analysis using qRT-PCR revealed that the expression level of the prostaglandin (PG) transporter gene SLCO2A1 decreased in response to microgravity, consistent with the meta-analysis of microarray datasets. Furthermore, the direction of gravity affected the expression level of SLCO2A1, buttressing the finding that its expression was affected by gravity. These results suggest that a meta-analysis of DNA microarray datasets may help identify new target genes previously overlooked in individual microarray analyses.


INTRODUCTION
Long-term space missions, such as those to Mars, have become more feasible. During space missions, astronauts experience microgravity, which is the condition of being weightless. During spaceflight, the astronauts' cardiovascular function undergoes tremendous changes, such as blood redistribution, increased compliance with lower extremity veins, heart contraction, and orthostatic intolerance (Delp et al., 2016;Tanaka et al., 2017). Thus, it is crucial to study the impact of spaceflight on various organs, such as the cardiovascular system. A recent study on twin astronauts demonstrated that a 1-year spaceflight induced various changes, such as carotid artery distension and increased intimamedia thickness (Garrett-Bakelman et al., 2019). This study suggests that retinal edema formation in space may be promoted by angiogenesis. Therefore, it is vital to elucidate the effect of microgravity on vascular endothelial cells for the coming space age.
Microgravity affects the morphology and physiology of the cardiovascular system. While the microgravity-induced changes are dramatic, they are thought to be mostly secondary to headward fluid shifts rather than direct effects of microgravity on the cardiovascular system (Norsk, 2020). Another factor that microgravity exerts on the orthostatic intolerance is a decrease in constrictor responses of mesenteric artery and vein (Behnke et al., 2013).
On the other hand, the endothelium is one of the most sensitive tissues to changes in gravity (Versari et al., 2013;Maier et al., 2015). Microgravity directly acts on vascular endothelial cells and affects processes, such as angiogenesis (Cotrupi et al., 2005). For example, microgravity causes the down-regulation of actin (Carlsson et al., 2003) and cytoskeletal disorganization (Villa et al., 2005) and decreases the metabolism of the endothelial cells (Kapitonova et al., 2012).
Several models, such as human subjects in bed rest, the suspended hind limbs of a rodent, clinorotated cultured cells, and humans, animals, and cultured cells in spaceflight (Fukuda and Iida, 2000;Perhonen et al., 2001;Globus and Morey-Holton, 2016), are used to evaluate the effect of microgravity on organs. On the other hand, DNA microarray is often performed using RNA extracted from various tissues to analyze the gene expressions affected by microgravity.
However, the studies examining the effect of microgravity on gene expression tend to yield inconsistent results. Although many investigations generate gene expression data, only a few have cross-sectionally examined data, especially in space biology and medicine. Two factors make such comparisons of gene expression among multiple datasets difficult. First, the studies used different microarrays to analyze different gene sets. Second, the nomenclatures of the same genes differ in the studies. We hypothesized that it was possible to identify common, essential microgravity-induced changes in gene expression from the results obtained using varying materials and methods. This study aimed to uncover the overlooked effects of microgravity on gene expression from accumulated valuable database resources by using different statistical methods. In this study, we applied a method developed by Kristiansson et al. (2013) to categorize genes into homologous groups and perform a meta-analysis to resolve the problem of genes having different names in various DNA arrays. Since microgravity causes not only cardiovascular deconditioning but also muscle atrophy, we analyzed the changes in gene expression in endothelial cells and skeletal muscle obtained from microarray datasets in multiple simulated or inflight microgravity experiments.
It would also be fascinating and important to identify the similarity in the impact of microgravity on the gene expression in different species. However, due to the variety of materials and methods used, it is difficult to perform a consistent analysis of the genes that respond to microgravity. To our knowledge, the similarity between the results from gene expression analysis obtained with different materials and methods has never been shown systematically in space biology. Here, we demonstrated such similarity by applying the t-distributed stochastic neighbor embedding (t-SNE) analysis (van der Maaten and Hinton, 2008), a method extensively used for visualizing high-dimensional data from flow cytometry to DNA array datasets.
Furthermore, after the gravity-sensitive genes were identified from the meta-analysis, their change in expression in human umbilical vein endothelial cells (HUVECs) was induced by simulated microgravity (SMG) using a clinostat and tested with qRT-PCR.
We found that the levels of LOX, SLCO2A1, and TXNIP expression were significantly altered in response to SMG. Also, we found that the direction of gravity affected the expression of SLCO2A1. Instead of being specific to the condition of microgravity, the expression of SLCO2A1 may be altered by gravity from a direction different from the arrangement of its cytoskeleton. These results suggest that the meta-analysis of accumulated microarray data is useful for identifying candidate genes with altered expression in response to gravity changes. This strategy could potentially contribute to the advancement of space biology and medicine.

Microarray Analysis
The gene expression data for endothelial cells were obtained from two studies, Gene Express Omnibus (GSE43582) (Versari et al., 2013) and Ma et al. (2013) 1 (Table 1), which investigated the effects of microgravity. In the GSE43582 study, HUVECs were cultured in a plastic cylindrical chamber. After the launch, the cells were exposed to microgravity for 7 days in space. The control group was exposed to 1 G on the ground for 7 days. The cells were preserved in RNAlater solution for later analysis. In the study conducted by Ma et al. (2013), EA.hy926 cells, a permanent cell line of the HUVEC, were cultured in a T-25 flask. The cells were exposed to SMG using the random positioning machine for 7 days. The adherent cells were scraped off for RNA isolation. Ma et al. (2013) made two comparisons. They compared the 1 G and SMG adherent groups, and between the 1 G and SMG threedimensional (3D) aggregate groups on day 7. Then, we performed a meta-analysis to test the null hypothesis that a gene of interest was not differentially expressed among all genes in the databases using a method described by Kristiansson et al. (2013). The crossexperimental p-values were corrected for multiple testing using the false discovery rate (FDR) determined via the Benjamini-Hochberg procedure. For the muscle samples, gene expression data from 12 experiments investigating the effects of microgravity were obtained from Gene Express Omnibus (Table 2). We performed the same meta-analysis and p-value correction for endothelial cells.

t-SNE Analysis
A t-SNE analysis was performed using R's tsne package 2 to group the multiple DNA array datasets. The perplexity parameter in the t-SNE analysis indicates the effective number of neighbors, and the result of the t-SNE plot may differ depending on its setting. The grouping by t-SNE analysis is shown to the same regardless of the value of perplexity by performing the t-SNE analysis with varying perplexity parameter 1, 2, or 3. In this study, the t-SNE analysis was performed at least 10 times at each perplexity parameter to confirm that the results showed a constant pattern. Furthermore, cluster analysis was performed using cluster 3.0 to 2 https://cran.r-project.org/web/packages/tsne/ group multiple DNA array datasets by a method different from the above t-SNE analysis (de Hoon et al., 2004).

Exposure to SMG
The HUVECs were seeded in the flasks and cultured for 2 days at 37 • C at 5% CO 2 to ensure the cells attach to the bottom of the flasks. Then, the flasks were split into two groups, SMG and 1 G control. The SMG flasks were clinorotated for 7 days using a clinostat (CL-1000, Zeromo; Kitagawa Iron Works Co., Ltd., Hiroshima, Japan) to achieve 10 −3 G by rotating samples on two axes to evenly distribute gravity in the incubator. During clinorotation, air bubbles were removed by filling the flasks with the culture medium to prevent shear stress. The 1 G control flasks were placed next to the clinostat and kept static for 7 days.

Cell Culture Under Different Gravitational Configurations
The HUVECs were seeded in the flasks simultaneously and cultured for 2 days to confirm the cells' attachment. Then, the flasks were divided into the "up, " "side, " and "down" configures, where the cells were on a horizontal surface, on a vertical surface, and hung on a horizontal surface, respectively. The "up" configuration was the control condition. This experiment was performed independently from the clinostat experiment in Section "Exposure to SMG."

Phase-Contrast Images of HUVECs
For the SMG experiments, the cells' condition was observed by acquiring daily phase-contrast images of the same areas in the culture flask using a TE2000 inverted microscope. For the experiment involving different gravitational configurations, the same areas in the culture flask were imaged on the day before exposure to the gravitational conditions and after 7 days of exposure to the conditions to observe the morphology and growth of the cells.

qRT-PCR
Total RNA was harvested from cells using a High Pure RNA Isolation Kit (Roche Diagnostics Indianapolis, IN, United States), according to the manufacturer's protocols. The RNA concentration was measured using a spectrometer (Nanodrop ND-1000, Thermo Fisher Scientific, United States). RNA samples with an OD 260/280 ratio of more than 1.80 were used for subsequent procedures. The samples were stored at −80 • C. For qRT-PCR, RNA samples were first heated at 70 • C for 5 min to destroy the secondary structures. Then, the 1 ng of RNA was reverse-transcribed into cDNA using a Verso cDNA Synthesis Kit (Thermo Scientific, Lithuania, United States). The cDNA was used for qRT-PCR to quantify the expression of 20 genes ( Table 3) using forward and reverse primers from previous publications or designed using the Primer3 (Version 0.4.0 3 ) software. qRT-PCR was performed using StepOnePlus thermal cycler (Thermo Fisher Scientific, United States). The following thermal cycling conditions were used: 95 • C for 10 min, followed by 40 cycles of 95 • C for 15 s and 60 • C for 60 s. The PCR primers are listed in Table 3. SYBR Green reagents (Life Technologies) were used for the relative quantification of RNA levels according to the manufacturer's instructions. The expression levels of target genes were normalized to that of 18S rRNA. Each measurement was obtained by averaging the measurements of three technical replicates. 3 http://bioinfo.ut.ee/primer3-0.4.0/

Statistical Analysis
Statistical analysis was performed using Prism (version 6.07, GraphPad Software, La Jolla, CA, United States). The data were expressed as mean ± standard error of the mean. The significance of the differences between the experimental groups was determined using the Student's t-test or Kruskal-Wallis test with Dunn's Multiple Comparison. A p-value of 0.05 or less was considered statistically significant.

Statistical Meta-Analysis of Gene Expression in Endothelial Cells in Response to Microgravity Conditions
We analyzed the gene expression data from microgravity exposure experiments performed in spaceflight or with clinorotation in the laboratory ( Table 1) to identify genes in endothelial cells with altered expression in response to microgravity. A total of 20,358 genes were retrieved from the spaceflight experiment, and 903 genes were retrieved from the clinorotation experiment. At first, GEO2R (Smyth, 2004) was used to identify the genes differentially expressed in the ground control and flight samples. The p-values obtained from this analysis were used with the p-values from the comparison between the 1 G and clinorotated adherent groups provided by Ma et al. (2013) for the subsequent meta-analysis. Applying the method described by Kristiansson et al. (2013), 108 genes with a p-value of 0.001 or less were identified, of which 46 genes had a FDR of 0.05 or less ( Figure 1A). We randomly chose 10 genes from the 46 genes for subsequent gene expression assays using qPCR. Of the 10 genes, CCN5, encoding WISP2 (WNT1 Inducible Signaling Pathway Protein 2), had the lowest p-value of 2.13 × 10 −8 , also the point with a high −log 10 FDR value in the volcano plot (Figures 1B,C). LOX (encoding lysyl oxidase) and MT1X (encoding metallothionein 1X) showed higher fold change (FC) values (logFC) at 4.11 and 1.71, respectively, in clinorotated EA.hy926 cells, indicating that microgravity facilitated the changes in their expression ( Figure 1C). However, the logFC values of LOX and MT1X were −0.35 and −0.45, respectively, in the space-flown HUVECs, indicating the decreased expression of these genes. The logFC value of SLCO2A1 (encoding solute carrier organic anion transporter family member 2A1) was −2.10 in the clinorotated EA.hy926 cells and −0.23 in the space-flown HUVECs, with an FDR of 5.18 × 10 −2 .
Then, we performed another meta-analysis using the same p-values from the space-flown HUVECs and those from another microgravity condition, the 3D aggregates in clinorotated EA.hy926 cells. As a result, 79 genes with a p-value of 0.001 or less were identified, of which 64 genes had a FDR of 0.05 or less (Figure 2A). Again, LOX showed a significant p-value of 2.95 × 10 −7 and a high logFC value of 3.44 in the clinorotated EA.hy926 cells, with the logFC of −0.35, in opposite direction in the space-flown HUVECs (Figures 2A-C). In contrast, SLCO2A1 showed a significant logFC value of −2.66 and an FDR of

Primer
Forward sequence Reverse sequence Reference Accession

Statistical Meta-Analysis of Gene Expression in Skeletal Muscle in Response to Microgravity Conditions
Next, we analyzed the gene expression data from the experiments using multiple microgravity loads on skeletal muscle ( Table 2). The p-values were obtained using GEO2R from each dataset. For example, the differential expression of genes across pre and postbed rest experiments in the soleus muscle was analyzed using the GSE14798 dataset. A meta-analysis of 12 datasets identified 6,328 genes with a FDR of 0.05 or less, of which 10 genes exhibited changes in expression in the same direction (both increasing or decreasing) in most of the datasets with valid FC values (Figure 3). These genes were selected for the subsequent gene expression assays using qPCR. Then, we performed the t-SNE analysis to gain insight into the similarity between the DNA array datasets. We analyzed 3,290 genes with valid FC values in all the datasets listed in Table 2. The individual datasets were marked A-L. The cost function in the t-SNE analysis reflected the difference between the distance of each gene in multidimensions, 3,290 dimensions in this case, and that in 2D. After performing the t-SNE analysis, the cost value converged between 1,000 and 5,000 steps ( Figure 4A). In step 10,000, the cost value converged to 0.480, 0.293, and 0.259 when the perplexity parameter was 1, 2, and 3, respectively. When the perplexity was 1, the array datasets were divided into two groups, i.e., A/B/C/D/E/I/J and F/G/H/K/L ( Figure 4B). A similar tendency was observed when the perplexity was 2 and 3 ( Figures. 4C,D). Furthermore, a similar pattern was observed in the results of the cluster analysis ( Figure 4I).
It would also be interesting to see how the microgravityinduced changes in the gene expression pattern in endothelial cells differ from those in muscle tissue. We investigated this by attempting to perform the t-SNE analysis by combining the 2 array datasets obtained from the endothelial cells ("X" and "Y" in Table 1) with the 12 array datasets obtained from skeletal muscle ("A-L" in Table 2). However, only 52 genes with valid FC values were identified in the combined datasets with 14 arrays. The removal of datasets A and B from the combined datasets resulted in the identification of 177 genes with valid FC values. We conducted the t-SNE analysis using these 177 genes. As a result, the cost value converged within 700 steps ( Figure 4E). In step 10,000, the cost value converged to 0.488, 0.321, and 0.206 when the perplexity was 1, 2, and 3, respectively. When the perplexity was 1, the array datasets were divided into three groups: C/F/H/I/J/X/Y, D/E, and G/K/L ( Figure 4F). A similar pattern was observed when the perplexity was 2 and 3 ( Figures 4G,H); the endothelium-derived datasets X and Y were relatively close to each other and within the same group.

Morphological Analysis of HUVECs Under SMG Using a Clinostat
The shape of the cells in the SMG was compared to those in the 1 G control groups after exposure to gravitational conditions. The groups exhibited similar cell shapes. Also, the difference in the cell numbers between the groups was not statistically significant, with 622.33 ± 264.66 in the SMG group and 563 ± 101.60 in the 1 G control group (p > 0.05) (Figures 5B-D).

Comparison of mRNA Expression Between Control and SMG Conditions by qRT-PCR
After exposing the HUVECs to different gravity conditions, we measured the mRNA level of the 10 genes selected from the DNA array meta-analysis of endothelial cells (Figure 6A). LOX expression was higher in the SMG group than in the 1 G group at 1.91 ± 0.27 and 1.20 ± 0.16, respectively (p < 0.05). Conversely, SLCO2A1 (SMG, 0.43 ± 0.12; 1G, 0.94 ± 0.05; p < 0.01) and TXNIP expression (SMG, 0.43 ± 0.10; 1G, 0.96 ± 0.06; p < 0.01) was lower under the SMG condition than the 1 G condition.
Next, we measured the expression of the 10 genes selected from the DNA array meta-analysis of muscle tissue (Figure 6B). In contrast to the result from genes selected from the meta-analysis of endothelial cell DNA arrays, there were no statistically significant differences in gene expression between the SMG and 1 G groups.

Morphological Analysis of HUVECs Under Different Gravitational Configurations
We tested the effects of the direction of gravity on the morphology and growth of endothelial cells by incubating HUVECs in flasks at the up, side, or down configuration ( Figure 7A). In the down configuration, the cells were less dense than the cells in the up and side configurations (Figures 7B-D). Moreover, the cells in the down configuration appeared "brushed" and "branched" (Figure 7D). The number of cells was the smallest in the down configuration at 209.4 ± 112.6 and largest in the side configuration at 638.6 ± 349.4 ( Figure 7E). The difference in cell numbers between the down and side configurations was statistically significant (p < 0.05).

Comparison of mRNA Expression Among the Up, Side, and Down Configurations by qRT-PCR
In order to confirm whether the changes in gene expression observed in the SMG could be observed also in the different flask configurations, gene expression analysis was performed using the gene set used in Figure 5A. After exposing HUVECs to different gravity configurations, we measured the expression of the 10 genes selected from the DNA array meta-analysis of endothelial cells (Figure 8). Kruskal-Wallis test revealed that gravity configuration affected the expression of SLCO2A1, with a statistically significant difference (Dunn's post hoc test, p < 0.05) in expression between the side and down configurations.

Meta-Analysis of Microarray Data
Researchers have used different types of species, organs, and conditions of gravity to evaluate microgravity's effect. However, FIGURE 3 | A meta-analysis of gene expression in muscle samples in response to microgravity conditions. The 10 genes exhibiting changes in expression in the same direction (all increasing or decreasing) in most of the datasets were listed. vastus l, vastus lateralis; gcnemius, gastrocnemius; EDL, extensor digitorum longus; longis, longissimus dorsi.  Table 2). (A) The convergence of cost values, reflecting the difference between the distance of each gene in multidimensions and that in 2D, over time at the perplexity of 1, 2, and 3. The t-SNE plots at perplexity = 1 (B), 2 (C), and 3 (D) are shown. (E-H) The t-SNE analysis of DNA array datasets obtained from endothelial cells (Table 1) and skeletal muscle ( Table 2). (E) The convergence of cost values over time at the perplexity of 1, 2, and 3. The t-SNE plots at perplexity = 1 (F), 2 (G), and 3 (H) are shown. Dashed lines are drawn to guide the grouping of datasets. (I) The cluster analysis of DNA array datasets obtained from skeletal muscle ( Table 2).  The gene sets were obtained from the meta-analyses of datasets from muscle samples. The relative quantity of each gene was measured using qRT-PCR and normalized to the level of 18S rRNA. For each condition n = 3. The data were expressed as the mean ± S.E.M. *p < 0.05 and **p < 0.01.  it is challenging to extract meaningful knowledge from these results. In this study, we performed a meta-analysis of gene expression using DNA microarray datasets obtained from experiments performed on diverse samples and species using various methods to mimic or achieve microgravity.
Interestingly, the t-SNE analysis revealed that gene expression patterns were grouped in DNA array datasets. For example, the datasets (G, H, K, and L) obtained from the skeletal muscle of space-flown mice formed a different cluster from datasets (A, B, C, D, E, and J) from skeletal muscle under SMG, except that one dataset from the space-flown sample (I) was in the "simulated microgravity" cluster and another dataset from SMG (F) was in the "spaceflight" cluster ( Figures 4B-D). Meanwhile, the datasets from hind limb unloading (D and E) were located far from the "spaceflight" cluster. This tendency was also observed when the number of analyzed genes was reduced to 177 (Figures 4F-H). However, it should also be noted that the difference in experimental conditions between the datasets G/H/K/L and the datasets D/E was recognized not only in the mode of microgravity but also in the duration of microgravity and the method of sample storage. Nevertheless, according to the t-SNE analysis, the gene expression pattern in datasets D and E differed from G, K, and L. This fact was also supported by the result of cluster analysis ( Figure 4I).
In addition, the DNA array datasets obtained from endothelial cells under actual (X) and simulated (Y) microgravity are located close to each other (Figures 4F-H), suggesting a remarkable similarity in gene expression patterns, given that the datasets come from independent research groups. This result provides a valid basis for estimating gene expression change under microgravity in space with SMG datasets in endothelial cells. These results suggest that when considering changes in gene expression, mixing different types of cells in a meta-analysis can cause problems with the accuracy of the results.

Changes in Gene Expression in Response to Gravity
Next, we investigated the genes whose expression may be affected by microgravity. The study results demonstrated that SLCO2A1 expression was decreased in both space-flown HUVECs (FC value = −0.23) and clinorotated EA.hy926 cells (FC value = −2.66) (Figure 2A). Furthermore, the qRT-PCR analysis revealed that SLCO2A1 expression was decreased in response to a 7-day exposure to SMG using a clinostat, in parallel to the results in space-flown HUVECs and clinorotated EA.hy926 cells (Figure 6). Importantly, the decreased expression of SLCO2A1 was not found in previous studies using space-flown HUVECs (Versari et al., 2013) and clinorotated EA.hy926 cells (Ma et al., 2013). In other words, the decreased expression of SLCO2A has been identified via a meta-analysis of multiple sets of previously published data.
Other genes responsive to microgravity have also been uncovered, such as LOX, encoding LOX with a copper-dependent amine oxidase activity and catalytic activity in the cross-linking of collagen and elastin (Vallet and Ricard-Blum, 2019), was affected. LOX expression was increased, consistent with the results in EA.hy926 cells (Ma et al., 2013). However, the expression of LOX is reduced in space-flown HUVECs (Versari et al., 2013). Such discrepancy may be due to the differences in cell types and the conditions of microgravity, spaceflight vs. clinorotation. Another gene whose expression is altered by gravity changes is LRP4, which is involved in forming and maintaining neuromuscular junctions (Barik et al., 2014). Interestingly, the expression of LRP4 was increased in this metaanalysis (Figure 3). It is tempting to speculate that one of the mechanisms of muscle atrophy under microgravity is the altered expression of LRP4; this hypothesis should be carefully studied in the future.
On the other hand, we analyzed the genes with identified sensitivity to gravity in the context of skeletal muscle but observed no changes in their expression in HUVECs exposed to SMG ( Figure 6B). This result suggests that modulation of gene expression is organ-specific. This is also suggested because the two DNA array datasets from endothelial cells are located outside of the muscle cluster in the t-SNE analysis.
Furthermore, the direction of gravity was found to affect the morphology, proliferation, and gene expression of HUVECs in this study. The decrease in the cell number in the "down" configuration is likely due to the change in the cell adhesion in response to gravity. However, the cells are unlikely to simply fall off, as the adhesion of extracellular matrix proteins to the culture flask is relatively stable. Cells are known to respond to gravity. In animals, the direction of gravity affects the proliferation of osteoblasts (Kacena et al., 2002).
Interestingly, in our study, osteoblast proliferation is inhibited in the inverted culture, similar to the HUVECs cultured in the down configuration. On the other hand, plant cells respond to gravity's direction, thus modulating the direction of root growth (Herranz and Medina, 2014). In plants, it is thought that the sedimentation of heavy statolith increases the tension of the actin cytoskeleton and subsequently activates mechanosensitive ion channels (Tatsumi et al., 2014). Suppose heavy nuclei and mitochondria also contribute to an animal cell's gravity sensing mechanism similar to plant statolith. In that case, this mechanism may explain the change in gene expression depending on the direction of gravity observed in the current study.
Intriguingly, SLCO2A1's expression was modulated in response to the direction of gravity in this study. In line with the results from previous studies using SMG, gravity appears to affect the expression of SLCO2A1, which encodes a prostaglandin (PG) transporter. Furthermore, it is reported that shear stress regulates the expression of SLCO2A1 in vascular endothelial cells (Topper et al., 1998). This fact supports the hypothesis that the expression of SLCO2A1 in vascular endothelial cells is modulated by mechanical stimuli, such as shear stress and gravity. The protein encoded by SLOC2A1 is involved in mediating the uptake and clearance of PGs, including PGE1, PGE2, PGD2, PGF2α, and TXB2 (Kanai et al., 1995;Lu et al., 1996). SLCO2A1 protein is expressed in vascular endothelial cells of various tissues, such as the heart, lung, kidney (Topper et al., 1998), brain (Dalvi et al., 2015;Nakamura et al., 2018), intestine (Dalvi et al., 2015;Yamaguchi et al., 2018), and skin (Syeda et al., 2012;Liu et al., 2015). PGE2 is one of the factors that promote angiogenesis (Zhang and Daaka, 2011). Interestingly, SLCO2A1, a transporter of the PGE2, is involved in the angiogenesis during wound healing (Liu et al., 2015) and the development of colon cancer (Nakanishi et al., 2017). Angiogenesis in microgravity is a subject that needs further research. Vyas et al. (2020) reported that the patterning of vasculature was decreased in the retinas of astronauts in spaceflight-associated neuro-ocular syndrome (SANS). Although it is tempting to speculate that the altered expression of SLCO2A1 in microgravity might be involved in the development of SANS, this hypothesis should be carefully tested in the future.

CONCLUSION
Using the meta-analysis of DNA microarray datasets obtained from the vascular endothelial cells, we have identified gene candidates whose expression levels change under microgravity. Accordingly to the t-SNE analysis, among the mixed datasets from endothelial cells and skeletal muscle, the gene expression patterns in response to gravity in endothelial cells form specific clusters. Experiments using endothelial cells under SMG and different directions of gravity confirmed that the expression of one of the candidate genes, SLCO2A1 encoding the PG transporter, is indeed altered. These results suggest that meta-analysis of multiple microarray datasets is useful for identifying novel target genes.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/ repositories and accession number(s) can be found in the article/Supplementary Material.

AUTHOR CONTRIBUTIONS
KT and KN contributed to the conception and design of the study. YiL and MW performed the experiments. MW, YiL, YuL, CW, and KT wrote sections of the manuscript. All authors contributed to manuscript revision, and all authors read and approved the submitted version.

FUNDING
This study was supported by Grant-in-Aid for Scientific Research on Innovative Areas (No. 15H05936) and for Challenging Exploratory Research (20K21896).