Open access peer-reviewed chapter

Gene Expression Profile of HDF in SMG Partially Overlaps with That in the NASA Twins Study

By Jade Q. Clement

Submitted: May 8th 2019Reviewed: August 1st 2019Published: September 26th 2019

DOI: 10.5772/intechopen.88957

Downloaded: 257


Microgravity research is an important field in biomedical sciences not only due to our interest in exploring and living in space, but also because of the insights it gives on earthbound health conditions. Using a human dermal fibroblast (HDF) cell line cultured in simulated microgravity (SMG) in combination with high throughput cDNA microarrays and quantitative Northern analysis, 271 differentially regulated genes were identified and 72% of these genes were also reported in the high throughput gene expression data of the recent National Aeronautics and Space Administration (NASA) Twins Study. The identification of the large number of overlapping microgravity sensitive genes between the skin fibroblast in microgravity and astronaut’s peripheral blood mononuclear cells (PBMCs) indicated that microgravity alone, without space radiation, was able to elicit an adaptive response involving a set of about 200 genes. Further analysis of the overlapping genes with the same direction of regulation (86 genes) and opposite direction of regulation (108 genes) revealed important pathways and cellular processes in the microgravity adaptation responses.


  • gene expression profiling
  • DNA microarray
  • northern blotting
  • rotating wall vessel (RWV)
  • rotary cell culture system (RCCS)
  • fibroblast
  • microgravity

1. Introduction

Humans have been traveling to space since 1961. During the close to 60-year period, hundreds of astronauts and cosmonauts have experienced microgravity as well as radiation in space. Exposure to microgravity environment has been shown to have potentially negative effects on human health. Some of the factors include cardiovascular deconditioning [1, 2], decline in immune response [3, 4], bone deterioration [5, 6, 7], and muscular atrophy [8, 9, 10].

Human immune response dysregulation has been shown during and after space flight [11]. Although acute onset life-threatening incidents directly caused by microgravity or spaceflight have not been reported, many conditions of health concern and symptoms related to immune dysfunction during orbital spaceflight have been described [12, 13]. Among the 70 tabulated clinical symptoms and medical conditions pertinent to immune dysfunction on board the International Space Station (ISS), skin rash and hypersensitivity account for 23 events; skin infections, 6 events; cold sores (caused by herpes simplex virus 1 infection), 6 events. About 50% of the tabulated incidents are of or closely related to skin symptoms or abnormal skin conditions during the long duration space flights [12]. Skin is the essential outer cover that protects the internal tissues and organs from potential physical, chemical, and biological assaults of the environment. In addition to immune cells and keratinocytes, dermal fibroblasts also play an important immunomodulation role such as in antimicrobial defense [14]. Epidermal keratinocytes can sense the presence of pathogen invasion and other environmental stimuli such as the presence of UV light and foreign chemicals and produce cytokines, chemokines, and growth factors in response. Communication between keratinocytes and dermal fibroblasts through cytokines is fundamental in skin immunity. A recent report shows that in a 3D skin model-based study, just keratinocytes and fibroblasts alone embedded in a collagen matrix are able to activate CD4+ T cells in response to microbial invasion [15]. The dermal fibroblasts play an essential role in antimicrobial response by integrating signals among cells in the skin.

To date, much microgravity research work has been done in ground-based research using microgravity analogs. Due to cost and limits to the technology much less has been done directly in the space environment. The recent NASA Twins Study is a tremendously important study because it is the first study that uses an integrated approach to study human adaptation to a space environment by documenting the molecular, physiological and cognitive effects during long term spaceflight [16]. The study also highlights the need for further study on important aspects such as vascular changes and immunological stress associated with the weightlessness of space flight [16].

Sudden gravity change has altered gene expressions from many cell types [17, 18, 19]. High-throughput gene expression analysis have great potential for application to research involving changes in environmental conditions [20]. Various high throughput studies such as cDNA microarrays and transcriptome RNA sequencing have been increasingly used to assess the mRNA levels in microgravity research [16, 17]. This is an effective approach because the control of mRNA abundance of genes is efficiently adapted by cells through controlling transcription (especially transcription initiation), nuclear pre-mRNA processing, mRNA transport, mRNA stability, etc. The cellular abundance of mRNAs is critical to gene function and protein production which are intriguingly fine-tuned by non-coding regulatory RNAs such as miRNAs. There have been many gene expression studies done on various cell lines grown both in space and using ground-based microgravity analogs. Many of these studies have yielded valuable data, but correlation of gene expression data between studies has been relatively low [17, 18].

To further understand the cellular and molecular mechanisms by which space flight alters skin immune defense activities such as those analyzed from the ISS [12], the effects of microgravity on various human skin cell lines need to be studied to identify the genes whose functions are altered by microgravity. In our previous study, we found expression changes in certain genes (such as HLA-G and IL-1β among many others) in response to simulated microgravity [21]. The current report is on the gene expression profile of HDF in response to SMG. Interestingly, a substantial overlap in gene expression profiles between the HDF under SMG and that from the human blood cells of the NASA Twins Study, especially the peripheral blood mononuclear cells (PBMCs) from inflight in the ISS. The comparative analysis yielded 194 differently expressed genes in both studies, of which 86 genes were regulated in the same direction (trend) while 108 genes were regulated oppositely. The significance of these findings was discussed.


2. Materials and methods

2.1 Simulated microgravity and cell culture

The HDF cell line, AG 1522, was generously provided by Dr. Honglu Wu of NASA. The HDF cell line displayed regular monolayer spindle shaped growth in conventional 2-D cell culture flasks. When HDF cells were subjected to SMG treatment, in a 3-D culture environment, they formed spherical aggregates. Ground based simulated microgravity was achieved using the 50 ml high aspect ratio vessels (HARVs) or rotating wall vessels (RWVs) of a rotary cell culture system (RCCS-4D) bioreactors from Synthecon, Inc. Cell viability and cell concentration were determined by Vi-Cell 1.01 cell counter of Beckman Coulter. For the three parallel experiments of 5 day modeled microgravity exposures, a density of 2.0 × 106 cells/ml with viability of 95.5% of the HDF AG1522 cells were cultured in RWVs at 20 rpm rotary setting to achieve the constant free-fall experience for cell aggregates. At the end of the five-day microgravity exposure period, the content of the bioreactor vessels was poured out into a 50 ml sterile centrifuge tube to collect cell pellet and 5 ml of the cell suspension from the bioreactor vessels were transferred to T75 flasks for morphological observation. Non-exposed stationary normal gravity control AG1522 cells were cultured in tissue culture flasks with vented caps (TPP Techno Plastic) in the same incubator at 37°C, 5% CO2.

2.2 Total RNA isolation and DNA microarray hybridization

HDF cells cultured in three SMG bioreactor vessels and control flasks were removed at the end of the five-day SMG exposure period, washed with phosphate buffered saline three times and lysed in Guanidinium Isothiocyanate Buffer. The cell lysates were stored at −80°C prior to ultracentrifugation for total RNA isolation [22, 23]. Total cellular RNA was labeled using the Agilent Low RNA Input Fluorescent Linear Amplification Kit following manufacturer’s protocols [24]. The fluorescently labeled cRNA probes were further purified and hybridized to Agilent 22 K Human Microarray V2 according to the specified procedures within the kit.

2.3 Microarray scanning, feature extraction and functional grouping

The microarrays were scanned using a ScanArray microarray scanner (Perkin-Elmer). The images generated from the scanning were imported into GenePix 6.0 (Molecular Devices, Sunnyvale, CA) for alignment and initial quantitation. The Gene Pix Results (GPR) files were then uploaded to CARMAweb [25] for normalization and statistical analysis. Background was subtracted and then each array was normalized using loess normalization within the array. A paired moderated T-Test was applied with Benjamini-Hochberg correction to control the false discovery rate [26]. Cut-offs were set at a P-value ≤0.01 and fold change of ≥1.5.

2.4 Northern blotting and quantitative gene expression analysis

Some of the significantly regulated microgravity sensitive genes identified from the DNA microarray analysis were further verified using Northern blotting. Briefly, 10 ug total RNA was loaded per lane on a 1% formaldehyde agarose gel for electrophoresis separation of RNA species. RNA Ladders from Fermentas Life Sciences were used as RNA size markers. The gel-separated RNAs were capillary transferred onto a nylon membrane which was subjected to a hybridization procedure using chemiluminescent (Pierce Biotechnology) labeled cDNA probe fragments. The blot was sequentially hybridized and striped and hybridized again with 12 different probe fragments. The cDNA probe fragments were generated from reverse transcription polymerase chain reaction (RT-PCR) using total cellular RNA as the templates. RT-PCRs were carried out using the Reverse Transcription System of Promega and the BioLine Red Polymerase PCR kit. Northern blot quantitation and association with the cDNA array results were done as described previously [21, 23].

2.5 Comparative analysis with the data from the NASA twins study

The NASA Twin gene expression data was identified from Suppemental Table 2 in the NASA Twins Study [16]. Prior to further bioinformatics analysis, the high throughput gene expression data from transcriptome RNA sequencing analysis for the inflight, first half and in-flight, second half was extracted from the NASA Twin study [16]. This data was then compared to the gene expression data from the current study with HDF in SMG. Mainly the transcriptome RNA sequencing data of the PBMC RNA from NASA Twins Study was compared with the HDF SMG data here, since it offered the most abundant overlapping genes. The PBMC data was extracted from the Excel spreadsheets and combined with the HDF data for further comparative analysis using similar previously published method [17].

2.6 Pathway and gene ontology analysis

To determine the Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathways [27, 28] and Gene Ontology, the genes that were determined to be differentially regulated were uploaded to the Database for Annotation, Visualization and Integrated Discovery (DAVID) [29]. The pathway and gene ontology information were used to build the tables. The process was done for both the HDF SMG data alone, and for its comparative analysis with the data from the NASA Twins Study. Prior to further bioinformatics analysis, the gene expression data for the inflight, first half and in-flight, second half was extracted from Supplemental Table 2 of the NASA Twins study [16]. This data was then compared to the gene expression data from the current study with HDF in SMG.

2.7 Constructing heat map and Venn diagram

The heat map was generated using Genesis 1.8.1 [30]. The Venn diagram was generated using the Excel plugin Array File Maker 4.0 [31].

3. Results

3.1 Gene expression profiling and the identification of microgravity sensitive genes from SMG treated HDF

HDF cells from each of the three SMG bioreactors and normal gravity controls were removed after 5-day SMG exposure for RNA extraction and microarray experiments. After normalization and statistical analysis (student ttest), the gene expression data were used to identify the initial set of significantly differentially regulated genes at the statistically significance level of P ≤ 0.01 and cut off point of ≥1.5 fold up or down regulation. The volcano plot ( Figure 1 ) shows the overall profile of the gene expression data from the three sets of parallel SMG experiments. Each dot on the volcano plot represents a gene selected in the initial set of differentially expressed ones if the dot is to the left (down) or right (up) of the red lines. The genes found outside of the vertical parallel red lines correspond to the genes differentially expressed by 1.5 fold or greater in the SMG. A total of 271 genes were identified from these three SMG experiments on HDF to be the initial set of significant genes that were used for further analysis as follows. The two most notable differentially regulated genes by the SMG evident on the plot were the matrix metallopeptidase 1 (MMP1) and connective tissue growth factor (CTGF) genes ( Figure 1 ).

Figure 1.

Volcano plot scattering the average M values (x axis) against the corrected p values (y axis) created using CARMAweb tool. Dots to the left of the red line and above the dashed line have a p-value of 0.01 or less and are down-regulated by at least 1.5 fold. The dots to the right of the red line and above the dashed lines have a p-value of 0.01 or less and are up-regulated by at least 1.5 fold.

The significantly (≥ 1.5 fold change, P ≤ 0.01) differentially expressed genes under simulated microgravity were further analyzed using Heatmap and Venn diagram to visually display the directions and centralization of gene expression. A high level of consistency of both the Heatmap and Venn diagram results were found among the microarray data from the three SMG bioreactors ( Figure 2 ). The Heatmap indicated that the expression levels of the three replicates were similar, with very minor variations in magnitude of expression. The Venn diagram further showed that there was a complete match among the three replicate microarrays from the three SMG bioreactor RNA samples. Thus, there was a very high level of consistence among data from the three SMG bioreactors. Among the 271 microgravity sensitive genes that differentially regulated by 1.5 fold or greater with a P value of ≤0.01, 129 were down-regulated and 142 were upregulated ( Figure 2B ).

Figure 2.

Heatmap and Venn diagram comparing the expression levels among the HDF cells cultured in three separate bioreactors. (A) The Heatmap shows the expression levels were very similar among the three replicate experiments. Green represents the down-regulated genes and red represents the up-regulated genes. (B) Venn diagram comparing the centralization level of the differentially regulated genes among each of the three replicates. The red represents the up-regulated genes and the green represents the down-regulated genes.

3.2 KEGG pathways of the microgravity sensitive genes from HDF

The identified 271 SMG sensitive genes were subjected to further bioinformatics analysis using the DAVID v6.8, which uses a modified Fisher Exact P-value for gene enrichment analysis and statistically determines the over-representation of functional gene categories in a gene list. P values equal to or smaller than 0.05 are considered significantly enriched [29]. Through the DAVID analysis, 16 statistically significant (p ≤ 0.05) KEGG Pathways from the SMG gene list were identified and genes in all the 23 KEGG pathways had >2 fold enrichment ( Table 1 ).

PathwaysP valueGene nameFE
hsa03010: Ribosome1.25 × 10−16Ribosomal Protein L7(RPL7), L9 (RPL9), L10A (RPL10A), L11 (RPL11), L12 (RPL12), L18A (RPL18), L27 (RPL27), L27A (RPL27A), L31 (RPL30), L32 (RPL32), L34 (RPL34), L35 (RPL35), L36 (RPL36), L39 (RPL39), S2 (RPS2), S3 (RPS3), S3A (RPS3A), S7 (RPS7), S10 (RPS10), S15A (RPS15A), S17 (RPS17), S18 (RPS18), S19 (RPS19), S29 (RPS29), S23 (RPS23), S24 (RPS24),Mitochondrial Ribosomal ProteinS6 (MRPS6)7.99
hsa04978: Mineral absorption9.38 × 10−7ATPase Na+/K+ Transporting Subunit Alpha 1 (ATP1A1), ATPase Na+/K+ Transporting Subunit Beta 3 (ATP1B3), Metallothionein 1A (MT1A), 1B (MT1B), 1G MT1G), 1H (MT1H), 1X (MT1X), 2A (MT2A), Ferritin Heavy Chain 1(FTH1), Ferritin Light Chain (FTL)9.14
hsa04510: Focal adhesion1.61 × 10−4Actin Beta (ACTB), Actin Gamma 1(ACTG1), Caveolin 1 (CAV1), Collagen Type 3 Alpha 1 Chain (COL3A1), Type I Alpha 1 Chain (COL1A1), Type 1 Alpha 2 Chain (COL1A2), Type 6 Alpha 3 Chain (COL6A3) Actinin Alpha 1, Myosin Light Chain 9 (MYL9), 12A (MYL12A), 12B (MYL12B), Calpain 2 (CAPN2),Integrin Beta 1(ITGB1), Filamin A (FLNA),, Vascular endothelial growth factor B (VEGFB),Fibronectin 1(FN1)3.12
hsa05205: Proteoglycans in cancer1.30 × 10−3ACTB, ACTG1,Cathepsin L (CTSL), Human leukocyte antigen A (HLA-A), B (HLA-B), C (HLA-C), G (HLA-G), Hepatocyte Growth Factor (HGS),Tubulin Beta 6 Class V (TUBB6), Tubulin alpha 1 B (TUBA1B),ITGB1,2.82
hsa04612: Antigen processing and presentation2.60 × 10−3Beta-2-Microglobulin (B2M), CTSL, HLA-A, HLA-B, HLA-C, HLA-G, Protein Family A (Hsp70) Member 1A Shock Protein (HSPA1A), Protein Family A (Hsp70) Member 8A (HSPA8)4.23
hsa05012: Parkinson’s disease2.61 × 10−3Peptidylprolyl Isomerase F (PPIF), NADH–ubiquinone oxidoreductase chain 1 (ND1), NADH–ubiquinone oxidoreductase chain 3 (ND3)ATP synthase F1 subunit epsilon (ATP5E), ATP Synthase Peripheral Stalk Subunit OSCP (ATP50),NADH:Ubiquinone Oxidoreductase Subunit B4 (NDUFB4),Ubiquinol-Cytochrome C Reductase, Complex III Subunit XI (UQCR11),Cyclooxygenase 1 (COX1) and 2 (COX2), Ubiquitin C-Terminal Hydrolase L1 (UCHL1), Ubiquitin B (UBB)3.12
hsa05416: Viral myocarditis2.67 × 10−3ACTB, ACTG1, CAV1,HLA-A, HLA-C, HLA-B, HLA-G4.94
hsa04145: Phagosome3.88 × 10−3ACTB, ACTG1,CTSL, HLA-A, HGS,TUBB6,HLA-C, HLA-B,TUBA1B,ITGB1, HLA-G2.95
hsa05130: Pathogenic Escherichia coliinfection8.15 × 10−3ACTB, ACTG1, ARPC2, TUBB6, TUBA1B,ITGB14.73
hsa04141: Protein processing in endoplasmic reticulum8.86 × 10−3Heat Shock Protein 90 Beta Family Member 1 (HSP90B1), Valosin Containing Protein (VCP), Defender Against Cell Death 1 (DAD1), DnaJ homolog subfamily A member 1 (DNAJA1), Protein Disulfide Isomerase Family A Member 4 (PDIA4) and 6 (PDIA6), HSPA1A, Heat Shock Protein Family A (Hsp70) Member 5 (HSPA5),CAPN2,HSPA8, Signal Sequence Receptor Subunit 1 (SSR1)2.62
hsa00190: Oxidative phosphorylation1.69 × 10−2ND1,ATP5E,NDUFB4,UQCR11,COX2, COX1, ND3,ATP5O, PPA12.72
hsa04530: Tight junction2.03 × 10−2ACTB, ACTG1, ACTN1, MYL12B, MYL12A,member RAS oncogene family(RAB13),MYL93.24
hsa04512: ECM-receptor interaction2.03 × 10−2CD44,COL3A1, COL6A3, COL1A2, COL1A1,ITGB1,FN13.24
hsa04974: Protein digestion and absorption2.13 × 10−2HSP90B1, VCP, DAD1, DNAJA1, PDIA6, HSPA1A, HSPA5, PDIA4,CAPN2,HSPA8, SSR13.2
hsa04670: Leukocyte transendothelial migration2.34 × 10−2ACTB, ACTG1, ACTN1, MYL12B, MYL12A,ITGB1,MYL9, THY12.8
hsa04260: Cardiac muscle contraction3.75 × 10−2UQCR11,ATP1B3, COX2, COX1, ATP1A1,tropomyosin 2 (beta)(TPM2)3.22
hsa04611: Platelet activation4.17 × 10−2ACTB, ACTG1, COL3A1, COL1A2, MYL12B, COL1A1, MYL12A,ITGB12.48
hsa05100: Bacterial invasion of epithelial cells4.32 × 10−2ACTB, ACTG1, CAV1, ARPC2,ITGB1,FN13.09
hsa05332: Graft-versus-host disease4.71 × 10−2HLA-A, HLA-C, HLA-B, HLA-G4.88
hsa05010: Alzheimer’s disease5.59 × 10−2tumor necrosis factor receptor superfamily member 1A(TNFRSF1A), ATP5E,NDUFB4,UQCR11,COX2, COX1,ATP5O, CAPN2, Calmodulin 2 (CALM2)2.16
hsa05330: Allograft rejection6.25 × 10−2HLA-A, HLA-C, HLA-B, HLA-G4.35
hsa05131: Shigellosis7.30 × 10−2ACTB, ACTG1,CD44,ARPC2,ITGB13.14
hsa04940:Type I diabetes mellitus8.45 × 10−2HLA-A, HLA-C, HLA-B, HLA-G3.83
hsa04918:Thyroid hormone synthesis9.43 × 10−2HSP90B1, ATP1B3, ATP1A1, HSPA5, PDIA42.87

Table 1.

KEGG pathway analysis of the 271 differential regulated genes from the HDF cells exposed to 5 days SMG.

Red indicates that the genes were up-regulated in the data and green indicates the genes were down-regulated in the data. The significance of gene enrichment is evaluated chiefly by the p value with a modified Fisher’s exact test. Fold enrichment (FE) is a measure of the magnitude of gene enrichment. P value ≤0.05 and FE ≥1.5 were considered significant and interesting.

The Ribosome KEGG Pathway included 27 of the SMG sensitive genes; 26 of these genes were down-regulated between 1.5 and 2 fold; only MRPS6 was upregulated (4.5 fold). Interestingly, a number of ribosomal protein genes have also been found to be down-regulated in other studies including our previous study on keratinocytes [21] as well as the NASA Twins Study [16]. The mineral absorption pathway had 10 genes which were all significantly up-regulated in SMG. ATP1A1 and ATP1B3 were both up-regulated by over 2.8 fold. The ferritin genes FTH1 and FTL were up-regulated by 6.6 and 1.6 fold, respectively. Ferritin genes have also been found to be up-regulated as a result of exposure to simulated microgravity in other studies [32]. The metallothionein genes were up-regulated from 7.5 to close to 12 fold. In a previous study on human keratinocytes in SMG, metallothioneins were also up-regulated [21]. Metallothionein isoforms MT1 and MT2 have been identified as gravity sensitive genes [17, 18, 21, 33, 34, 35]. Focal adhesion pathway genes were found to be generally down-regulated by more than 2 fold in the present study. The only exceptions were ITGB1 and VEGFB which were up-regulated by 1.5 and 1.8 fold, respectively. There were 14 genes in the proteoglycans in cancer category, with 9 genes being up-regulated and 4 down-regulated; these included several major gravity sensitive genes identified previously, CAV1, FN1, DCN, and CD44 [17]. In the antigen processing and presentation pathway, the 8 genes represented were all up-regulated. These include the HLA genes, the closely related B2M, as well as CTSL, HSPA1A and HSPA8. In a previous study, HLA-G is also up-regulated in the SMG treated human keratinocytes [21]. In normal gravity environment, HLA-G has been shown to have direct inhibitory effect on T, APC, and NK cell functions and induces suppressor T-cells [36]. It has been found to be present in neurological disorders [37]. HLA-G is considered a stress inducible gene [38]. HLA-G also plays a role in tumor-driven immune escape mechanism of cancer cells during the later phase in host and tumor cell interactions. When HLA-G is expressed, it can result in an immune suppressive function [38, 39, 40]. HLA-G expression has been found to correlate with low frequency of rejection in some forms of organ transplants [39]. Heat shock proteins HSPA1A and HSPA8 have been identified as being gravity sensitive in several studies. In some of these studies, HSPA1A has been shown to be up-regulated [21, 32, 41]. They were also shown to be down-regulated [42].

3.3 Validation of microarray results through northern blotting analysis

The quality of a gene list from a high throughput study is essential for a successful functional analysis in DAVID [29]. The high throughput microarray data of HDF in SMG presented above was further validated by performing Northern blot analysis.

Northern blotting measures the abundance as well as the size of the RNA of interest [21]. In agreement with the microarray data, the Northern results showed that the mRNA levels of β-actin, CTGF, and VIM were down regulated, while that of FTH1, HLA-A, HLA-B (data not shown), HLA-G, MMP1, and MT2A were upregulated ( Figure 3 ). In addition, Northern analysis showed that the expression level of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was consistent in the presence and absence of microgravity ( Figure 3I ).

Figure 3.

Microarray data validation using northern blot analysis. 10 μg of total RNA of the fibroblast cells from each of the three RWV bioreactors (3-D spheres) were loaded in lanes 2–4. Lanes 1 and 5 are control RNA samples from the cells grown at normal gravity. The northern blot was probed and striped repeatedly using cDNA probe fragments from genes indicated at the bottom of each panel (A–I). RNA size markers in nucleotides were labeled along the left side of each panel.

3.4 Comparison of HDF in SMG gene expression data with that from the NASA twins study

The NASA Twins study is very substantial in undertaking with extensive data sets [14]. It was very appealing to compare the current data of HDF in SMG with the recently published gene expression data from the NASA Twins study. Since a large amount of data is available, it is more manageable to first focus on the gene expression data from the samples taken from inflight first half (up to 6 months) and the second half (up to 1 year), generated from the inflight astronaut’s PBMCs. Such comparison of the two gene expression data sets gave 194 overlapping differentially expressed genes, or about 72%, of the genes found to be differentially regulated in the current study with HDF in SMG were also differentially regulated in the selected data set from the NASA study. Of these 194 genes, 86 had a similar expression pattern (regulated in the same direction) to the HDF in SMG data and 108 had the opposite expression pattern. A Heatmap was generated as a way of better visualizing the similarities between the data points ( Figure 4 ).

Figure 4.

Heatmap comparison of HDF in SMG gene expression data with that from the NASA twins study Heatmap is used to compare the HDF SMG data (column A) with the inflight, first half of the twins study ribodepleted (column B), PolyA+ (column C), multivariant (column D), and the inflight, second half of the twins study ribodepleted (column E), PolyA+ (column F), and multivariant (column G). Green were down-regulated genes, red represented up-regulated genes, and gray indicates the presence of no corresponding data.

3.5 KEGG pathways of the overlapping genes

As a way of comparing and better understanding the potential relationships between the HDF SMG gene expression data and the NASA Twin gene expression data from Table 2 [16], we uploaded the gene list of the genes that showed same direction expression patterns ( Table 2 ) and the genes that showed opposite expression patterns ( Table 3 ) to DAVID in order to generate KEGG Pathway information.

PathwaysP ValueGene nameFE
hsa03010:Ribosome3.93 × 10−27RPL35, RPL27A, RPL36, RPS15A, RPS2, RPL39, RPS3, RPL32, RPL7, RPL31, RPS3A, RPL9, RPL34, RPL11, RPL10A, RPL12, RPS23, RPS24, RPL27, RPS7, RPS18, RPS19, RPL18A, RPS17, RPS1022.18
hsa04530:Tight junction5.32 × 10−3ACTB, ACTG1, MYL12B, MYL12A, MYL96.936
hsa05130:Pathogenic Escherichia coli infection8.07 × 10−3ACTB, ACTG1, ARPC2, TUBA1B9.465
hsa04670:Leukocyte transendothelial migration1.40 × 10−2ACTB, ACTG1, MYL12B, MYL12A, MYL95.247
hsa04510:Focal adhesion2.55 × 10−2ACTB, ACTG1, MYL12B, MYL12A, CAPN2, MYL93.515
hsa04810:Regulation of actin cytoskeleton2.74 × 10−2ACTB, ACTG1, ARPC2, MYL12B, MYL12A, MYL93.448
hsa04921:Oxytocin signaling pathway3.32 × 10−2ACTB, MYL6, ACTG1, CALM2, MYL94.023
hsa04611:Platelet activation8.88 × 10−2ACTB, ACTG1, MYL12B, MYL12A3.713
hsa05131:Shigellosis9.53 × 10−2ACTB, ACTG1, ARPC25.657

Table 2.

KEGG pathways generated from comparison of genes differentially regulated in the same direction between the HDF and NASA twins data.

The gene symbols in green color indicated down-regulation. The significance of gene enrichment is evaluated chiefly by the p value with a modified Fisher’s exact test. FE is a measure of the magnitude of gene enrichment. P value ≤0.05 and FE ≥1.5 were considered significant and interesting, respectively.

PathwaysP valueGene nameFE
hsa04141:Protein processing in endoplasmic reticulum2.93 × 10−5HSP90B1, VCP, DAD1, DNAJA1, PDIA6, HSPA1A, HSPA5, PDIA4, HSPA8, SSR15.81
hsa04512:ECM-receptor interaction1.70 × 10−4CD44,COL3A1, COL6A3, COL1A2, COL1A1, ITGB1, FN17.91
hsa04612:Antigen processing and presentation7.62 × 10−4HLA-A, HLA-C, HSPA1A, HLA-B, HSPA8, B2M7.76
hsa04978:Mineral absorption7.86 × 10−4ATP1B3, ATP1A1, MT1X, FTH1, FTL11.2
hsa04974:Protein digestion and absorption1.48 × 10−3ATP1B3,COL3A1, COL6A3, COL1A2,ATP1A1,COL1A16.7
hsa04918:Thyroid hormone synthesis4.41 × 10−3HSP90B1, ATP1B3, ATP1A1, HSPA5, PDIA49.1
hsa04510:Focal adhesion1.37 × 10−2CAV1, COL3A1, COL6A3, COL1A2, COL1A1, ITGB1, FN17.02
hsa05134:Legionellosis1.48 × 10−2ARF1, VCP, HSPA1A, HSPA83.44
hsa04151:PI3K-Akt signaling pathway1.67 × 10−2HSP90B1,COL3A1, COL6A3, COL1A2,YWHAQ ,COL1A1, ITGB1, FN1,DDIT43.34
hsa05416:Viral myocarditis1.71 × 10−2CAV1,HLA-A, HLA-C, HLA-B6.9
hsa04144:Endocytosis2.74 × 10−2CAV1,ARF1, HLA-A, HLA-C, HSPA1A, HLA-B, HSPA82.56
hsa05169:Epstein–Barr virus infection2.92 × 10−2CD44,VIM,HLA-A, HLA-C, HLA-B4.64
hsa05332:Graft-versus-host disease3.96 × 10−2HLA-A, HLA-C, HLA-B4.06
hsa05205:Proteoglycans in cancer4.23 × 10−2CAV1,CD44, CD63, DDX5,ITGB1, FN12.85
hsa05330:Allograft rejection4.87 × 10−2HLA-A, HLA-C, HLA-B8.93
hsa04514:Cell adhesion molecules (CAMs)5.37 × 10−2HLA-A, HLA-C, HLA-B,ITGB17.97
hsa04145:Phagosome6.32 × 10−2HLA-A, HLA-C, HLA-B,ITGB13.46
hsa04940:Type I diabetes mellitus6.60 × 10−2HLA-A, HLA-C, HLA-B3.28
hsa03050:Proteasome7.16 × 10−2PSMB4, PSMB7, PSMA67.02
hsa05320:Autoimmune thyroid disease9.55 × 10−2HLA-A, HLA-C, HLA-B6.7
hsa05145:Toxoplasmosis9.79 × 10−2PPIF, HSPA1A, ITGB1, HSPA85.67

Table 3.

KEGG pathways generated from a comparison of genes differentially regulated in the opposite direction between the HDF and NASA twins data.

The green indicates genes that were down-regulated and the red indicates genes that were up regulated in the HDF data set. The significance of gene enrichment is evaluated chiefly by the p value with a modified Fisher’s exact test. FE is a measure of the magnitude of gene enrichment. P value ≤0.05 and FE ≥1.5 were considered significant and interesting, respectively.

After processing the list of genes with the same expression patterns through DAVID, a total of 7 KEGG pathways were statistically significant (P ≤ 0.05) Interestingly, the genes enriched in all 9 KEGG pathways had >3.5 fold enrichment ( Table 2 ). All the genes represented in the KEGG pathways were down-regulated in both sets of data. The Ribosome pathway had the greatest number of genes represented. The down-regulation of the ribosomal protein genes is consistent with our previous study in keratinocytes in SMG [21].

KEGG pathway analysis of the list of genes with opposite direction of expression regulation produced a total of 15 KEGG Pathways that were statistically significant (P ≤ 0.05). Interestingly, the genes in all the 21 KEGG pathways had >2.5 fold enrichment ( Table 3 ).

The protein processing in endoplasmic reticulum (ER) pathway had the greatest number of genes represented. In the current study of HDF in SMG, the genes represented in the ER pathway were all up-regulated whereas they were down-regulated in the NASA Twin Study data [16]. Heat shock proteins HSPA1A and HSPA8 have been identified as being gravity sensitive in several studies. In some of these studies, HSPA1A has been shown to be up-regulated [21, 32, 41] and in one study in addition to the NASA Twin Study was shown to be down-regulated [42]. Collagens also seem to be mixed in differential regulation in microgravity. In our current study, COL1, 2, and 6 were down-regulated but were up-regulated in the Twin study. In some studies COL1 has been shown to be up-regulated [43] while in others it has been shown to be down-regulated [33, 44]. FN1, which we identified as a putative “space gene” [17, 18] has also been shown to have variations in expression patterns. It has been shown to be down-regulated in several studies [34, 44, 45] and up-regulated in others [16, 43, 46].

4. Discussion

In the current study, the gene expression profile of HDF grown in 5 day SMG was first displayed and validated ( Figures 1 3 , Table 1 ). The high throughput cDNA microarray data of HDF in 5 day SMG was then used to compare with the high throughput RNA sequencing data from an astronaut’s PBMCs during a long term inflight ISS ( Figure 4 , Tables 2 and 3 ). Amazingly, about 72% of the 271 microgravity sensitive genes of HDF in SMG, were also differentially regulated in the NASA Twins 6- and 12-month inflight data. These 194 overlapping genes were identified as putative microgravity sensitive space genes, because the human dermal fibroblast cell line was only exposed to SMG, no radiation nor other space related environmental factor was involved. However, other factors such as cell type difference and space radiation, may also influence the expression of microgravity sensitive genes. Indeed, among these microgravity sensitive genes, 86 genes showed the same expression pattern in both simulated and real microgravity conditions while the other 108 genes displayed opposite direction of regulation.

It is remarkable that as many as 86 genes were found to have the same directions of expression regulation between very different settings of studies. It is most likely that these genes were the main players in cellular response to microgravity environment. When the 86 microgravity sensitive genes with the same expression regulation trends were subjected to KEGG pathway analysis, they were represented in seven significant pathways where they were all downregulated ( Table 2 ). Most notably, both sets of global gene expression data showed the down regulation of 25 ribosomal protein genes. The genes in pathogenic Escherichia coli infection pathway and leukocyte trans-endothelial migration pathway were all down-regulated in microgravity ( Table 2 ), which may contribute to the decreased immune resistant to microbial infection in spaceflight. In addition, genes in the cytoskeleton network (ACTB, ACTG1, ARPC2, MYL12B, MYL12A, MYL9), focal adhesion (ACTB, ACTG1, MYL12B, MYL12A, CAPN2, MYL9), as well as tight junction (ACTB, ACTG1, MYL12B, MYL12A, MYL9) pathways were also downregulated in both sets of data. In combination with the downregulation of extracellular matrix proteins ( Table 1 ) such as COL1A1, COL1A2, COL3A1 in HDF, the data indicated an overall decrease in bone matrix and skeletal muscle synthesis and increased catabolism (e.g. MMP1 increased sharply). Furthermore, genes in the oxytocin signaling pathway (ACTB, MYL6, ACTG1, CALM2, MYL9), which is involved in smooth muscle contraction and stress management, were also down in their expression levels in both microgravity data sets. Malfunction of this pathway has been implicated in depression, autism, and schizophrenia [47]. Overall, the data in Table 2 gave a strong mechanistic connection to the main symptoms, such as skin problems and immunological stress, vascular changes, muscle atrophy and bone density alteration that were associated with the weightlessness of space flight. Indeed, the altered expression of these 86 microgravity sensitive genes affected many fundamental molecular functions (data not shown), including structural constituent of ribosome, RNA binding, protein binding, metal ion binding, structural constituent of the cytoskeleton, cadherin binding involved in cell-cell adhesion, extracellular matrix binding, etc. Many biological processes in these cells (data not shown), such as SRP-dependent co-translational protein targeting to membrane, translation initiation and elongation, mRNA stability, muscle contraction, regulation of cell shape, among others, were also significantly impacted. These genes with common trend of expression regulation in microgravity would most likely expand the list of putative major space genes and microgravity sensitive pathways [17, 18]. A substantial amount of information was derived from the current work which may necessitate more detailed analysis and discussion in future communications.

The number of overlapping microgravity sensitive genes was substantial considering the many differences between the two study settings. The result from this comparative analysis further validated the effectiveness of the bioreactors for SMG cell culture. The identification of 86 genes ( Table 2 ) with the same direction of regulation in two different study settings is very substantial and unique. These genes should be considered best candidates for major microgravity sensitive genes because one of the two studies, the current study, only involves simulated microgravity, while the other study, the NASA Twins Study involves true spaceflight environment with microgravity and space radiation. This comparative analysis here enabled the differentiation of the microgravity effect alone on the differentially expressed genes from the human astronaut spaceflight gene expression data. The identification of the overlapping significant genes regulated in the opposite direction rendered important insight into human gene activity changes in very different study systems. The HDF cell line in SMG versus the human astronaut in ISS, adjusted their expressions toward adaptation to the simple SMG as well as the true space environment of both space microgravity and space radiation. The 108 microgravity sensitive genes with opposite directions of expression regulation could also be of major significance in the microgravity adaptation process. Compared to the single cell line in SMG alone for the HDF cells, the cell samples from the inflight astronaut was exposed to various other factors such as space radiation, in addition to the microgravity of the ISS. The more complex space environment may require the significant genes to modify their expression toward adaptation. Expression patterns in this group of genes could provide insight into our understanding regarding the interplay among different cellular gene functions in human adaptation to microgravity and space radiation ( Table 3 ).

Our previous studies on gene expression in HEK cells grown in SMG has 43 genes overlapping with the HDF data in the current communication; of which, 23 genes were regulated in the same direction and 20 were regulated in the opposite direction [21]. These two different types of cells require different culture conditions and perform different roles in the skin. It is understandable that they have their characteristic expression profiles in response to the simulated microgravity environment. However, the number of overlapping genes were also substantial. In a previous review paper comparing various microarray based gene expression studies on microgravity effects, an initial list of 129 genes were identified as putative microgravity sensitive genes or major space genes [17]. In the current study, 12 out of the 194 genes that were significantly differentially regulated in both the HDF cells and the PBMCs, are also in the group of the putative major spaces genes [17], with 4 genes regulated in the same direction (MMP1, GPNMB, RPL10A, and ANXA2) and 8 genes (CAV1, CD44, CD59, CYR61, FN1, HSPA1A, MT1X, and PDIA4) changed their expression in the opposite directions in response to microgravity. Continued microgravity research in space and the readily controlled simulated microgravity bioreactors would provide valuable information toward the identification of major gravity sensitive genes, or simply, the major space genes. With more data available, the molecular and cellular mechanisms underlying the microgravity response could be better understood. Elucidation the molecular mechanism of human space adaptation response is an important aspect toward safer space experience and human health in general. It is evident that continued microgravity research is beneficial to healthy living in space and on earth.

5. Conclusion

The identification of 271 genes of HDF significantly differentially regulated by SMG provided a set of data for more detailed mechanistic studies; 72% of these microgravity sensitive genes were also reported in the high throughput gene expression data in the recent NASA Twins’ Study.

The identification of the large number of overlapping genes between the HDF in SMG and astronaut’s PBMCs indicates microgravity alone, without space radiation, was able to elicit an adaptive response involving a set of about 200 genes.



This work was partially supported by NASA JSC Grant No. NCC9-165.


ACTBactin beta
ACTG1actin gamma 1
ANXA2annexin A2
ATP1A1ATPase Na+/K+ transporting subunit alpha 1
ATP1B3ATPase Na+/K+ transporting subunit beta 3
ATP50ATP synthase peripheral stalk subunit OSCP
ATP5EATP synthase F1 subunit epsilon
CALM2calmodulin 2
CAPN2calpain 2
CARMAwebcomprehensive R based microarray analysis web frontend
CTGFconnective tissue growth factor
CTSLcathepsin L
CYR61cysteine rich angiogenic inducer 61
DADdefender against cell death 1
DAVIDdatabase for annotation, visualization and integrated discovery
DNAJA1DnaJ homolog subfamily A member 1
FLNAfilamin A
FTH1ferritin heavy chain 1
FTLferritin light chain
GAPDHglyceraldehyde-3-phosphate dehydrogenase
GOgene ontology
GPNMBglycoprotein NMB
HARVshigh aspect ratio vessels
HDFhuman dermal fibroblast
HLAhuman leukocyte antigen
HSPheat shock protein
ISSinternational space station
ITGB1integrin beta 1
KEGGkyoto encyclopedia of genes and genomes
MMP1matrix metallopeptidase 1
MRPS6mitochondrial ribosomal protein 6
MYLmyosin light chain
NDNADH–ubiquinone oxidoreductase chain
NDUFB4NADH:ubiquinone oxidoreductase subunit B4
PBMCsperipheral blood mononuclear cells
PDIAProtein disulfide isomerase family A member
PPIFpeptidylprolyl isomerase F
RAB13member RAS oncogene family
RCCSrotary cell culture system
RPribosomal protein
RWVrotating wall vessel
SMGsimulated microgravity
SSR1signal sequence receptor subunit 1
TNFRSF1Atumor necrosis factor receptor superfamily member 1A
TPM2tropomyosin beta chain
TUBA1Btubulin alpha 1 B
TUBB6tubulin beta 6 class V
UBBubiquitin B
UCHL1ubiquitin C-terminal hydrolase L1
UQCR11ubiquinol-cytochrome C reductase, complex III subunit XI
VCPvalosin containing protein
VEGFBvascular endothelial growth factor B

© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Jade Q. Clement (September 26th 2019). Gene Expression Profile of HDF in SMG Partially Overlaps with That in the NASA Twins Study, Gene Expression and Phenotypic Traits, Yuan-Chuan Chen and Shiu-Jau Chen, IntechOpen, DOI: 10.5772/intechopen.88957. Available from:

chapter statistics

257total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Environmental Factors Affecting the Expression of Bilateral-Symmetrical Traits in Plants

By Sergey Baranov, Igor Vinokurov and Lubov Fedorova

Related Book

First chapter

Introductory Chapter: Gene Editing Technologies and Applications

By Yuan-Chuan Chen

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us