Spaceflight missions affect the behavior of microbes that are inevitably introduced into space environments and may impact astronauts’ health. Current studies have mainly focused on the biological characteristics and molecular mechanisms of microbes after short-term or long-term spaceflight, but few have compared the impact of various lengths of spaceflight missions on the characteristics of microbes. Researchers generally agree that microgravity (MG) is the most critical factor influencing microbial physiology in space capsules during flight missions. This study compared the growth behavior and transcriptome profile of
Keywords
- simulated microgravity
- length of time
- phenotype
- transcriptomics
The field of space life science has become a new research frontier following the advent of space exploration and the development of space technology (Belobrajdic et al. 2021; Macaulay et al. 2021). During space exploration, microbes are inevitably introduced into space habitats by astronauts and may grow on surfaces in the space capsule (Singh et al. 2018; Checinska Sielaff et al. 2019). Due to extreme environmental factors, such as microgravity (MG), high vacuum, and intense radiation, microbes in space undergo a series of changes in gene expression at the transcriptional level and eventually perform global physiological alterations, including metabolism, growth rate, motility, virulence, and biofilm formation (Taylor 2015; Senatore et al. 2018). For instance, strains of
MG is defined as a state of weightlessness that occurs due to decreased physical force exerted by gravity, where gravity ranges from 10–3 to 10–6
Moreover, several studies have shown that results obtained from HARV experiments are similar to those obtained in a true MG. For example,
Fig. 1
Schematic diagram of mechanical principle of HARV. MG is simulated by rotating the samples around a plane perpendicular to the gravity vector.

Fig. 2
Growth curve of PML (red) and PMS (blue). The growth curves of PML and PMS were quantified by measuring the OD600 values every 2 h for 24 h.

OD600 value of PML and PMS.
Time (h) | PML | PMS |
---|---|---|
0 | 0.548 ± 0.005 | 0.547 ± 0.006 |
2 | 0.581 ± 0.008 | 0.591 ± 0.005 |
4 | 0.674 ± 0.010 | 0.726 ± 0.005 |
6 | 0.863 ± 0.026 | 0.955 ± 0.012 |
8 | 1.206 ± 0.018 | 1.411 ± 0.026 |
10 | 1.560 ± 0.022 | 1.636 ± 0.013 |
12 | 1.709 ± 0.019 | 1.806 ± 0.027 |
14 | 1.776 ± 0.016 | 1.872 ± 0.024 |
16 | 1.826 ± 0.021 | 1.938 ± 0.025 |
18 | 1.876 ± 0.018 | 1.983 ± 0.015 |
20 | 1.927 ± 0.018 | 1.955 ± 0.014 |
22 | 1.900 ± 0.019 | 1.957 ± 0.028 |
24 | 1.889 ± 0.018 | 1.916 ± 0.014 |
to determine the quantity of the biofilm. The results showed that PML decreased the OD570 value in contrast to PMS (0.096 ± 0.006 vs. 0.118 ± 0.011,
Fig. 3
Biofilm formation ability of PML (red) and PMS (blue). Biofilm formation ability of PML and PMS was examined by measuring the OD570 values in microtiter plates.

Fig. 4
Swarming motility of PML (red) and PMS (blue). Swarming motility of PML and PMS was determined by measuring the translucent zone diameters in swarming agar.

Fig. 5
The Venn diagram of common and specific genes in PML and PMS.

DEGs were determined using the criteria of abs log2 (FoldChange) ≥ 1 and adjusted
Fig. 6
Scatter plot of DEGs. The x-axis and y-axis represent the logarithm of the gene expression for PML and PMS, respectively. Red spots and blue spots represent upregulated and downregulated genes, respectively.

Statistics of differentially expressed genes (DEGs).
Gene ID | Expression (up/down) | Log2 FoldChang (PML/PMS) | Adjusted |
Gene products |
---|---|---|---|---|
up | 1.308 | < 0.001 | cysteine sulfinate desulfinase | |
up | 1.156 | < 0.001 | cysteine desulfuration protein | |
up | 1.323 | 0.004 | choline trimethylamine-lyase | |
up | 1.686 | 0.002 | alkaline phosphatase | |
up | 1.077 | 0.003 | phage shock protein A | |
up | 1.107 | 0.002 | phage shock protein B | |
up | 1.423 | 0.049 | phosphate transport system ATP-binding protein | |
up | 1.060 | 0.023 | large subunit ribosomal protein | |
up | 1.439 | 0.002 | phosphate transport system ATP-binding protein | |
up | 1.377 | < 0.001 | phosphate transport system ATP-binding protein | |
up | 1.228 | 0.010 | phosphate transport system ATP-binding protein | |
up | 1.575 | 0.002 | phosphate transport system ATP-binding protein | |
up | 1.148 | < 0.001 | fimbria/pilus periplasmic chaperon | |
up | 1.031 | 0.022 | cysteine desulfurase | |
up | 1.064 | 0.009 | sensor histidine kinase | |
up | 2.720 | < 0.001 | small multidrug resistance pump | |
up | 1.424 | 0.017 | putative phosphotransacetylase | |
up | 1.473 | 0.023 | amino acid transport system substrate-binding protein | |
up | 1.520 | < 0.001 | Cd2+/Zn2+-exporting ATPase | |
down | –1.009 | < 0.001 | arginine transport system ATP-binding protein | |
down | –1.290 | 0.022 | cytochrome o ubiquinol oxidase subunit I | |
down | –1.011 | 0.007 | D-erythrose 4-phosphate dehydrogenase | |
down | –1.330 | 0.016 | 3-hydroxyacyl-CoA dehydrogenase | |
down | –1.045 | 0.047 | acyl-CoA dehydrogenase | |
down | –1.277 | < 0.001 | ferrous iron transport protein A | |
down | –1.465 | 0.005 | glycerol-3-phosphate dehydrogenase | |
down | –1.348 | 0.023 | glycerol-3-phosphate dehydrogenase | |
down | –1.221 | 0.003 | glycerol-3-phosphate transporter | |
down | –1.080 | 0.018 | acetolactate synthase I/III small subunit | |
down | –1.178 | 0.003 | D-methionine transport system permease protein | |
down | –1.021 | < 0.001 | peptide-methionine oxide reductase | |
down | –1.181 | 0.046 | thiosulfate reductase | |
down | –1.099 | < 0.001 | PTS system, glucose-specific IIB component | |
down | –1.105 | < 0.001 | protein NrfC | |
down | –1.210 | 0.048 | cytochrome ubiquinol oxidase subunit III | |
down | –1.283 | 0.036 | cytochrome ubiquinol oxidase subunit IV | |
down | –1.255 | 0.028 | – | |
down | –1.008 | 0.048 | toxin CptA | |
down | –2.137 | < 0.001 | major pilin subunit PapA | |
down | –1.001 | < 0.001 | phosphorelay signal transduction system | |
down | –1.190 | < 0.001 | TetR/AcrR family transcriptional regulator | |
down | –1.029 | < 0.001 | putative transport protein | |
down | –1.315 | 0.017 | Mat/Ecp fimbriae adhesion | |
down | –1.430 | 0.012 | Mat/Ecp fimbriae outer membrane usher protein | |
down | –1.349 | 0.037 | Mat/Ecp fimbriae periplasmic chaperone | |
down | –2.151 | 0.048 | Mat/Ecp fimbriae adhesion | |
down | –1.003 | 0.027 | cation/acetate symporter | |
down | –1.250 | < 0.001 | glycerol uptake facilitator protein | |
down | –1.018 | 0.020 | minor fimbrial subunit | |
down | – 1.060 | 0.023 | putative oxidoreductase |
Fourteen categories, including 50 DEGs (identical DEGs may fall into different categories), were identified between PML and PMS according to the KEGG pathway classification (Fig. 7). Compared with PMS, PML is characterized by regulating genes related to cellular community, membrane transport, energy metabolism, carbohydrate metabolism, lipid metabolism, and amino acid metabolism.
Fig. 7
KEGG pathway analysis of DEGs. The x-axis and y-axis represent the numbers of DEGs and KEGG pathway category, respectively.

Notably, some DEGs were found to belong to diverse functional categories including stress response (
It is the first study to analyze the growth behavior and transcriptome profile of
The effects of MG on the growth rate and final cell density of microorganisms are controversial. Previous studies have shown that microorganisms have an increased growth rate, with a growth curve characterized by a shorter lag phase and a more prolonged exponential phase under MG or SMG. For example,
Biofilms produced by microbes in the water system and device surface of the spacecraft cabin represent a significant phenotype concerning virulence and antibiotic resistance and induce corrosion of metals and structural materials (Novikova et al. 2006; Liu 2017). Therefore, it is necessary to analyze the biofilm formation ability of microbes under real and simulated MG. A representative investigation demonstrated that the biofilm of
Motility is a critical microbial behavior that plays a crucial role in nutrient absorption, tissue localization and invasion, biofilm formation, virulence, and survival (Acres et al. 2021). Alterations in microbial motility often originate from changes in the gene expression of the flagella-related system (Morgenstein et al. 2010). Although few studies have focused explicitly on the effects of spaceflight and SMG on the phenotype of microbial motility, genetic and transcriptomic analyses have provided a general understanding of gene expression that can predict microbial motility according to the genes involved in flagella function (Milojevic and Weckwerth 2020).
This study presents the only data available to compare long-term versus short-term SMG exposure effects on
China’s space station will be completed in the early 2020s, and spaceflight missions will require astronauts to remain on the space station for long periods. With the reality of long-duration missions lasting years, avoiding infectious diseases induced by microbes and maintaining the safety of the space station are of great significance in the future. Due to the challenges of sending microbial samples into the space environment, most research still needs to be done in laboratories using ground-based bioreactors under SMG. Considering that one of the new challenges facing astronauts is reaching increasingly long-duration spaceflight targets, studying the influences of long-term versus short-time SMG exposure on the growth behavior of
Fig. 1

Fig. 2

Fig. 3

Fig. 4

Fig. 5

Fig. 6

Fig. 7

Statistics of differentially expressed genes (DEGs).
Gene ID | Expression (up/down) | Log2 FoldChang (PML/PMS) | Adjusted |
Gene products |
---|---|---|---|---|
up | 1.308 | < 0.001 | cysteine sulfinate desulfinase | |
up | 1.156 | < 0.001 | cysteine desulfuration protein | |
up | 1.323 | 0.004 | choline trimethylamine-lyase | |
up | 1.686 | 0.002 | alkaline phosphatase | |
up | 1.077 | 0.003 | phage shock protein A | |
up | 1.107 | 0.002 | phage shock protein B | |
up | 1.423 | 0.049 | phosphate transport system ATP-binding protein | |
up | 1.060 | 0.023 | large subunit ribosomal protein | |
up | 1.439 | 0.002 | phosphate transport system ATP-binding protein | |
up | 1.377 | < 0.001 | phosphate transport system ATP-binding protein | |
up | 1.228 | 0.010 | phosphate transport system ATP-binding protein | |
up | 1.575 | 0.002 | phosphate transport system ATP-binding protein | |
up | 1.148 | < 0.001 | fimbria/pilus periplasmic chaperon | |
up | 1.031 | 0.022 | cysteine desulfurase | |
up | 1.064 | 0.009 | sensor histidine kinase | |
up | 2.720 | < 0.001 | small multidrug resistance pump | |
up | 1.424 | 0.017 | putative phosphotransacetylase | |
up | 1.473 | 0.023 | amino acid transport system substrate-binding protein | |
up | 1.520 | < 0.001 | Cd2+/Zn2+-exporting ATPase | |
down | –1.009 | < 0.001 | arginine transport system ATP-binding protein | |
down | –1.290 | 0.022 | cytochrome o ubiquinol oxidase subunit I | |
down | –1.011 | 0.007 | D-erythrose 4-phosphate dehydrogenase | |
down | –1.330 | 0.016 | 3-hydroxyacyl-CoA dehydrogenase | |
down | –1.045 | 0.047 | acyl-CoA dehydrogenase | |
down | –1.277 | < 0.001 | ferrous iron transport protein A | |
down | –1.465 | 0.005 | glycerol-3-phosphate dehydrogenase | |
down | –1.348 | 0.023 | glycerol-3-phosphate dehydrogenase | |
down | –1.221 | 0.003 | glycerol-3-phosphate transporter | |
down | –1.080 | 0.018 | acetolactate synthase I/III small subunit | |
down | –1.178 | 0.003 | D-methionine transport system permease protein | |
down | –1.021 | < 0.001 | peptide-methionine oxide reductase | |
down | –1.181 | 0.046 | thiosulfate reductase | |
down | –1.099 | < 0.001 | PTS system, glucose-specific IIB component | |
down | –1.105 | < 0.001 | protein NrfC | |
down | –1.210 | 0.048 | cytochrome ubiquinol oxidase subunit III | |
down | –1.283 | 0.036 | cytochrome ubiquinol oxidase subunit IV | |
down | –1.255 | 0.028 | – | |
down | –1.008 | 0.048 | toxin CptA | |
down | –2.137 | < 0.001 | major pilin subunit PapA | |
down | –1.001 | < 0.001 | phosphorelay signal transduction system | |
down | –1.190 | < 0.001 | TetR/AcrR family transcriptional regulator | |
down | –1.029 | < 0.001 | putative transport protein | |
down | –1.315 | 0.017 | Mat/Ecp fimbriae adhesion | |
down | –1.430 | 0.012 | Mat/Ecp fimbriae outer membrane usher protein | |
down | –1.349 | 0.037 | Mat/Ecp fimbriae periplasmic chaperone | |
down | –2.151 | 0.048 | Mat/Ecp fimbriae adhesion | |
down | –1.003 | 0.027 | cation/acetate symporter | |
down | –1.250 | < 0.001 | glycerol uptake facilitator protein | |
down | –1.018 | 0.020 | minor fimbrial subunit | |
down | – 1.060 | 0.023 | putative oxidoreductase |
OD600 value of PML and PMS.
Time (h) | PML | PMS |
---|---|---|
0 | 0.548 ± 0.005 | 0.547 ± 0.006 |
2 | 0.581 ± 0.008 | 0.591 ± 0.005 |
4 | 0.674 ± 0.010 | 0.726 ± 0.005 |
6 | 0.863 ± 0.026 | 0.955 ± 0.012 |
8 | 1.206 ± 0.018 | 1.411 ± 0.026 |
10 | 1.560 ± 0.022 | 1.636 ± 0.013 |
12 | 1.709 ± 0.019 | 1.806 ± 0.027 |
14 | 1.776 ± 0.016 | 1.872 ± 0.024 |
16 | 1.826 ± 0.021 | 1.938 ± 0.025 |
18 | 1.876 ± 0.018 | 1.983 ± 0.015 |
20 | 1.927 ± 0.018 | 1.955 ± 0.014 |
22 | 1.900 ± 0.019 | 1.957 ± 0.028 |
24 | 1.889 ± 0.018 | 1.916 ± 0.014 |