Abstract
Introduction:
The microgravity environment astronauts experience during spaceflight can lead to an increased risk of oral diseases and possible changes in oral microecology. In this study, we aimed to assess changes in the microbial community of supragingival plaques to explore the effects of spaceflight microgravity environment on oral microecology.
Methods:
Sixteen healthy male volunteers were recruited, and supragingival plaque samples were collected under −6° head-down bed rest (HDBR) at five-time points: day 1 before HDBR; days 5, 10, and 15 of HDBR; and day 6 of recovery. Bacterial genomic DNA was sequenced using gene sequencing technology with 16S ribosomal ribonucleic acid V3–V4 hypervariable region amplification and the obtained data were analyzed bioinformatically.
Results:
Alpha diversity analysis showed a significant increase in species richness in supragingival plaque samples on day 15 of HDBR compared with that at pre-HDBR. Beta diversity analysis revealed that the community composition differed among the groups. Species distribution showed that, compared with those at pre-HDBR, the relative abundances of Corynebacterium and Aggregatibacter increased significantly during HDBR, while those of Veillonella, Streptococcus, and Lautropia decreased significantly. Moreover, compared with those at pre-HDBR, the relative abundance of Leptotrichia increased significantly on day 6 of recovery, whereas the relative abundances of Porphyromonas and Streptococcus decreased significantly. Network analysis showed that the interaction relationship between the dominant genera became simpler during HDBR, and the positive and negative correlations between them showed dynamic changes. Phylogenetic investigation of communities by reconstruction of unobserved states analysis showed that the amino acid metabolism function of plaque microorganisms was more enriched during HDBR.
Discussion:
In summary, in a 15-day simulated microgravity environment, the diversity, species distribution, interaction relationship, and metabolic function of the supragingival plaque microbial community changed, which suggests that microgravity may affect the oral microecosystem by changing the balance of supragingival plaque microbial communities and further leading to the occurrence and development of oral diseases.
1 Introduction
With the continuous exploration of space fields by humans, an increasing number of manned space missions have been successfully performed, and the construction of commercial spaceports is ongoing (Grenon et al., 2012). However, during space flights, the microgravity environment experienced by astronauts can affect human health and safety, including the musculoskeletal, hematopoietic, endocrine, and digestive systems (White and Averner, 2001; Graebe et al., 2004; Grimm et al., 2016).
Oral health is closely related to general health. The oral cavity is the gateway for pathogens to invade the human body and is a reservoir for infection-associated microorganisms that play a role in the spread of disease (Johnston and Dietlein, 1977). Oral dysbiosis can promote the development of oral diseases such as periodontal disease and caries, and may also affect systemic health by migration and colonization of pathogens, bacteremia, stimulating the body‘s inflammatory response and autoimmunity (Atarashi et al., 2017; Aarabi et al., 2018; Kwun et al., 2020). Therefore, oral health maintenance and monitoring of microbial communities are important for astronauts. According to a report, pain caused by pulpitis or gingivitis in astronauts during air flight may seriously affect their work efficiency (Brown et al., 1977). Astronauts experience increased dental calculus and gingival inflammation during space flights (Brown et al., 1976). In 1978, a Russian astronaut developed toothache symptoms during the last 2 weeks of a 96-day flight on Salyut 6; however, no emergency plans were made (Ball and Evans, 2001). A case of dental caries was reported in a medical incident at the Space Station MIR, which was treated with temporary fillings (Gontcharov et al., 2005). Rai (2009) confirmed through an 8-h−6° head-down-tilt bed rest experiment in 20 volunteers that a microgravity environment can cause toothache symptoms in healthy volunteers during tooth occlusion. Meanwhile, the results of space flight and simulation experiments showed that microgravity conditions could change the growth, morphology, virulence, antibiotic resistance, and secondary metabolism of microorganisms (Klaus and Howard, 2006; Wilson et al., 2007; Huang et al., 2018); change cellular interactions with the host (Foster et al., 2013); and cause negative effects on the immune system (Crucian et al., 2015). During life on the space station, astronauts’ normal microflora composition changes significantly, which results in their bodies becoming less resistant to the source of infection (Ilyin, 2005).
Microbial factors are key factors leading to oral infectious diseases such as periodontal disease and caries. Plaque biofilms are thought to initiate the development and progression of periodontal disease (Marsh, 2012). Plaque accumulation at the gingival margin may cause chronic inflammation in periodontal tissues. Disruption of the homeostatic relationship between the plaque bacteria and the host can initiate and promote periodontal disease progression (Zhang et al., 2020). A close relationship exists between the stimulatory activity of Toll-like receptors (TLR)2 and TLR4 in supragingival plaques and the periodontal status (Yoshioka et al., 2008). By counting the number of plaque microorganism clones over 2 years in 24 participants, it was found that the stability of plaque microbial composition may be a good predictor of periodontal health (Kumar et al., 2006). Furthermore, plaque microorganisms mediate the development of caries. Bacterial microorganisms in the community produce acids to erode hard dental tissues, resulting in their demineralization and, ultimately, the formation of cavities. The accumulation of supragingival plaque biofilms on teeth is usually the first manifestation of caries (Bowen, 2016). Generally, a dynamic balance exists between plaque microbes and the host, and changes in health status alter the composition of substances released by the bacteria and correspondingly alter host responses (Darveau et al., 1997).
In a previous study, our group found that the structure of salivary microbial communities in young men in the 15-day head-down bedridden state changed significantly, and differences were observed in species distribution (Sun et al., 2022). However, microbial structures can vary across different microecological sites within the oral cavity (Yang et al., 2021). Compared to saliva, plaque microbiota accumulates on non-detached teeth surfaces, which provides them with more time to develop into uniform and complex communities (Shi et al., 2018). Simultaneously, teeth have non-shedding surfaces suitable for microbial growth, which are key ecological determinants contributing to the persistence of supragingival plaque (Marsh, 2003). Various microorganisms in the plaque community interact to form complex biological networks that provide ideal sites for microbial growth and reproduction. At present, the research on the microbial community of human oral cavity in microgravity environment mainly focuses on saliva samples, while the research on supragingival plaque is very limited. Therefore, studying the microbial community structure of supragingival plaque is crucial.
The −6° head-down bed rest model is considered a well-established model for simulating the physiological effects of microgravity in astronaut spaceflight (Pavy-Le Traon et al., 2007; Hargens and Vico, 2016). In 2015, the International Academy of Astronautics (IAA) has organized a team to study and publish a standardized guideline of bed-rest studies in the spaceflight context (Sundblad et al., 2016), which clearly proposes the use of −6° head-tilt to simulate microgravity. This experiment was performed under the guidance of this guideline.
Early scientists used culture-based methods to collect supragingival plaque samples from astronauts before and after space flights and found that their microbial counts were altered (Brown et al., 1974, 1976; Johnston and Dietlein, 1977). However, due to the limitation of experimental techniques, the previously investigated microbial species are limited to a few culturable species. Meanwhile, their findings were limited to changes in the counts of specific microorganisms. These do not reflect the overall changes in the microbial community of astronauts’ supragingival plaque in the spaceflight microgravity environment, and do not report the co-occurring relationship between oral flora and predictions of microbial metabolic function. High-throughput sequencing technologies can provide new means to comprehensively understand the overall changes in the microbial community of supragingival plaques (Jagathrakshakan et al., 2015; Al-Hebshi et al., 2019). In addition, previous studies on microgravity have been limited to surveys before and after the flight or before and after simulation experiments and have not explored changes in plaque microflora at multiple time points during flight or experiments. Here, we used the −6° head-down bed rest method to simulate the microgravity environment and used 16S ribosomal ribonucleic acid (rRNA) gene sequencing technology based on the Illumina MiSeq platform to assess the changes in the microbial community of supragingival plaques at multiple time points and provide a basis for achieving precise microbiological prevention and treatment of space oral diseases in the future.
2 Materials and methods
2.1 Participants
The participants were openly recruited from the public for this study. The inclusion criteria were as follows: male; age 18–45 years; height 160–175 cm; weight 50–80 kg; physical health status. The exclusion criteria included dental caries, periodontal disease, other oral infectious diseases, smoking history, drug dependence history, mental and psychological disease history, and systemic disease history. After a physical examination and screening, 16 male participants were recruited, and none withdrew from the experiment. This study protocol conformed to the ethical principles of the Declaration of Helsinki and was approved by the Ethics Review Committee of the Space Science and Technology Institute (Shenzhen) (SISCJK202009001). All participants provided informed consent, signed an informed consent form, and received remuneration after the end of the experiment.
2.2 −6° head-down bed rest test
In this experiment, we used the −6° head-down bed rest (HDBR) method to simulate a microgravity environment in space. During the experiment, participants held the longitudinal axis of their body with a horizontal line at −6° head down. The experiment lasted for 28 days, including 7 days of adaptation pre-HDBR, 15 days of HDBR, and 6 days of recovery post-HDBR. Two participants were housed in the same room with beds separated by a movable curtain. Room temperature was maintained between 22°C and 26°C. The participants were provided oral hygiene education and training on the bass brushing method before the experiment. During HDBR, it was required that: (1) participants perform daily activities in bed, such as eating, bathing, and urinating; they could turn on their sides but not get up. (2) The participants woke up at 07:00 and retired at 22:30 every day. The schedule of the daily activities was strictly controlled. (3) Participants were provided a uniform diet according to nutritional standards by dietitians. No beverage, except pure water, was offered. (4) Specialists monitored the participants’ physical conditions and physiological changes to ensure that their physical indicators were within normal limits. This experiment was grouped according to 5-time points: day 1 before HDBR (AP group), during HDBR: day 5 of HDBR (BP group), day 10 of HDBR (CP group), day 15 of HDBR (DP group), and day 6 of recovery (EP group).
2.3 Samples collection
Supragingival plaque samples were collected at 9 a.m. Participants were asked to rinse their mouths and remove food debris. Following moisture removal using cotton rolls, supragingival plaque samples were collected from the buccal, lingual, and interproximal surfaces of six teeth, including the maxillary and mandibular first molars, the right maxillary central incisor, and the left mandibular central incisor, using sterile dental probes. The samples were placed in 1.5 mL Eppendorf tubes containing 1 mL of pre-chilled phosphate buffered saline, transferred on an ice bath, and stored in a−80°C freezer.
2.4 DNA extraction and purification
Target deoxyribonucleic acid (DNA) sequences in the samples were obtained using the cetyltrimethylammonium bromide (CTAB) method, and the purity and concentration of the DNA were detected using 1% agarose gel electrophoresis. An appropriate amount of DNA and a certain amount of sterile water were mixed in a centrifuge tube to dilute the sample to 1 ng/μl. Diluted genomic DNA was used as a template for PCR amplification using forward primer 341F and reverse primer 806R for the 16S rRNA V3–V4 region. PCR reactions used 15 μL Phusion® High-Fidelity PCR Master Mix (New England Biolabs, Ipswich, MA, USA), 2 μL forward and reverse primers and approximately 10 ng of template DNA. Thermal cycling consisted of initial denaturation at 98°C for 1 min, followed by 30 cycles of denaturation at 98°C for 10 s, annealing at 50°C for 30 s, and elongation at 72°C for 30 s. Finally, denaturation was performed at 72°C for 5 min. Samples were mixed equally with different tags and mixed thoroughly, and PCR amplification products were detected by 2% agarose gel electrophoresis. Target bands were tapped to obtain the recovered products. The products recovered by tapping were purified using the GeneJET Gel Extraction Kit (Qiagen, Hilden, Germany).
2.5 Sequencing and data analysis
Library construction was performed using a TruSeq® DNA polymerase chain reaction (PCR)-free sample preparation kit (Illumina, CA, USA). A Qubit fluorometer (Thermo Fisher, Carlsbad, CA, USA) and quantitative-PCR (Q-PCR) analysis were used to evaluate library quality. After the library was qualified, it was processed for sequencing using the Illumina NovaSeq platform (Nevogene, Beijing, China). Each data sample was split from the offboard data according to the barcode and PCR amplification primer sequences, and the sequencing reads of each sample were spliced, controlled, and filtered using FLASH (v1.2.7) after truncating the barcode and primer sequences to obtain clean tags (Magoč and Salzberg, 2011; Bokulich et al., 2013). The quality control process was completed according to the method described by Qiime (v1.9.1) (Caporaso et al., 2010; Poncheewin et al., 2019), and effective tags were obtained (Rognes et al., 2016). All effective tags were clustered using Uparse software (v7.0.1001) to cluster sequences into operational taxonomic units (OTUs) based on the similarity between sequences with 97% consistency (Haas et al., 2011). With SILVA based on the Mothur algorithm, the database performs species annotation analysis of OTUs sequences (setting a threshold of 0.8–1) to obtain the microbial community composition of each sample at different taxonomic levels (Edgar, 2013; Quast et al., 2013).1 Species accumulation curves were plotted using the R software (v2.15.3). Observed species, Chao1, ACE, Shannon, Simpson, and good coverage indices were calculated using the Qiime software based on species annotation. Fast multiple sequence alignment was performed using the MUSCLE software (v3.8.31) to obtain phylogenic relationships for representative sequences of all OTUs (Edgar, 2004). Principal coordinate analysis (PCoA) was performed by calculating the UniFrac distances using Qiime software, using phylogenetic relationships among OTUs (Lozupone et al., 2011). An unweighted pair-group method with an arithmetic mean (UPGMA) sample cluster tree was constructed. Multidimensional data were visualized by performing non-metric multidimensional scaling (NMDS) analysis based on the Bray-Curtis distance (Lozupone et al., 2007). The corresponding images were drawn using R software. Analysis of dissimilarities (ADONIS), analysis of similarities (ANOSIM), and analysis of molecular variance (AMOVA) were used to test the significance of the differences in community structure among the groups (Roewer et al., 1996; Chapman and Underwood, 1999; Anderson, 2001). Statistical analyses were performed using the SPSS software (v21.0.0.0). To determine whether species diversity and relative abundance differed among groups, we performed normality tests using the Shapiro–Wilk method, one-way analysis of variance (ANOVA), and Tukey multiple tests if a normal distribution was met; otherwise, the Friedman rank sum test and Dunn’s multiple tests were used. GraphPad Prism (v8.0.1) and R software were used for mapping. Network relationships among the genera were plotted using Spearman’s correlation coefficient and Graphviz software (v2.38.0). Significantly different metabolic pathways were identified by the phylogenetic investigation of communities by reconstruction of unobserved states (PICRUSt). Differences in metabolic function between groups were analyzed using a paired sample t-test. Statistical significance was set at p < 0.05.
3 Results
3.1 Sequencing data processing
The raw data and quality control results for all samples after being subjected to double-end sequencing on the Illumina NovaSeq platform are shown in Supplementary Table 1. Raw off-machine data yielded an average of 87,621 sequences per sample. After splicing and filtering low-quality and short-length sequences, chimeras were removed, and an average of 50,627 high-quality sequences with an average length of 418 bp were obtained. A total of 3,948,939 high-quality valid sequences were obtained for subsequent analysis.
3.2 OTUs clustering and species annotation
OTUs clustering analysis was performed on valid sequence data from supragingival plaque samples according to 97% concordance to analyze the composition and distribution of species in the supragingival plaque samples (Supplementary Figure 1). The clustering results showed that 4,941 OTUs were obtained, and 4,080 (82.57%) were annotated to the database. The OTUs at each level represented 2 bacterial boundaries, 29 phyla, 44 classes, 96 orders, 176 families, 379 genera, and 412 species. A Venn diagram showed that 1,436, 1843, 1787, 2,918, and 1833 OTUs were identified in the AP, BP, CP, DP, and EP groups, respectively, after 97% consistency treatment. Five groups of overlapping regions, i.e., 560 OTUs, were shared by the five groups (Figure 1A). Analyzing five groups of unique OTUs, we found 332, 343, 286, 986, and 521 unique OTUs in the AP, BP, CP, DP, and EP groups, respectively. The number of unique OTUs in the DP group was significantly higher than in the AP group. The results of the cumulative box plot of species indicated that the sample size in this experiment was sufficient to meet the analytical requirements for species richness (Figure 1B).
Figure 1
3.3 Alpha diversity analysis
Alpha diversity reflects species complexity within the sample community, and the results of each index are presented as box plots (Figure 2). The Chao1 index was used to estimate the number of species contained in the supragingival plaque samples, i.e., species richness. The results showed that the number of species in the samples showed an increasing trend with an increase in the-6° head-down bed rest experimental time compared to the AP group. The DP group had the largest number of species, which was significantly different from the AP, BP, and CP groups (p < 0.05). The number of species in the EP group did not differ significantly from those in the other four groups. The observed species index also showed that, at the same sequencing depth, the highest number of species was found in the DP group, which was significantly different from that in the AP group (p < 0.05). The Shannon index results showed no statistical difference in species diversity during HDBR compared with the AP group. Although species diversity in the EP group decreased significantly compared with the DP group (p < 0.05), it was not statistically different from that in the AP group. The above results showed that alpha diversity differed among the samples, with the DP group having the highest species richness.
Figure 2
3.4 Beta diversity analysis
Beta diversity analysis was performed to determine whether microgravity altered the microbial community structure. PCoA and NMDS analyses showed that the EP group was distinct from the other four groups in different directions, and the community structure was quite different (Figures 3A,B). The beta diversity index heatmap and UPGMA cluster analysis showed that the community structures of the AP and EP groups were quite different from those of the other four groups; however, the structures of the three groups during HDBR were similar (Figures 3C,D). To further verify the correctness of the above results, we also performed tests for differences between the groups (Table 1). These results showed that the community structures of the AP and EP groups were significantly different from those of the other four groups of samples under a simulated microgravity environment (p < 0.05).
Figure 3
Table 1
| Group | Anosim | Adonis | Amova | |||
|---|---|---|---|---|---|---|
| R | p | R2 | p | Fs | p | |
| AP-BP | 0.1012 | 0.016 | 0.06707 | 0.017 | 1.33637 | 0.094 |
| AP-CP | 0.1364 | 0.003 | 0.07504 | 0.003 | 1.80369 | 0.011 |
| AP-DP | 0.1304 | 0.005 | 0.07517 | 0.004 | 2.58057 | <0.001 |
| AP-EP | 0.1746 | 0.001 | 0.08985 | 0.006 | 3.68565 | <0.001 |
| BP-CP | 0.02656 | 0.219 | 0.05628 | 0.07 | 0.673817 | 0.905 |
| BP-DP | 0.00344 | 0.378 | 0.04035 | 0.204 | 1.04837 | 0.34 |
| BP-EP | 0.1966 | 0.001 | 0.09544 | 0.002 | 4.06971 | <0.001 |
| CP-DP | 0.08408 | 0.036 | 0.06574 | 0.012 | 1.14806 | 0.253 |
| CP-EP | 0.2694 | 0.002 | 0.11565 | 0.003 | 4.83831 | <0.001 |
| DP-EP | 0.1492 | 0.001 | 0.07624 | 0.004 | 4.29182 | <0.001 |
Dissimilarity analysis of supragingival plaque samples between groups.
Anosim analysis and Adonis analysis based on Bray-Curtis distance; Amova analysis based on Unifrac distance. In Anosim analysis, R > 0 indicates significant differences between groups, and R < 0 indicates that intra-group differences are greater than inter-group differences. In Adonis analysis, a larger value of R2 indicates a higher degree of explanation of differences by group. In Amova analysis, Fs is the F-test value. The reliability of the above statistical analyses were expressed as p values, and p < 0.05 indicated that the statistical differences were significant.
3.5 Species distribution
According to the species annotation results, 28 phyla were identified in the five groups of supragingival plaques at the phylum level. Columnar accumulation plots of the relative abundance of species drawn from the top 10 phyla with the maximum abundance were selected (Figure 4A). Among them, there were six dominant phyla (mean relative abundance >1%): Fusobacteria, Bacteroidetes, Proteobacteria, Actinobacteria, Firmicutes, and Spirochaetes according to their relative abundances from high to low, and the total abundance of the dominant phyla was 94.74%. Statistical analysis of differences in the relative abundance of dominant phyla revealed that compared to the AP group, Actinobacteria were significantly increased in the CP group (p < 0.05), whereas Firmicutes were significantly decreased in the CP (p < 0.01) and DP groups (p < 0.05); Actinobacteria were significantly decreased in the EP group compared with the CP group (p < 0.05) (Figure 4B).
Figure 4
In total, 379 genera were identified. Columnar accumulation plots of the relative abundance of the species were drawn for the top 10 most abundant genera (Figure 4C). There were 15 dominant genera with a total abundance of 80.10%. Statistical analysis of differences in the relative abundance of dominant genera revealed that compared with the AP group, Aggregatibacter were significantly increased (p < 0.01) and Lautropia were significantly decreased (p < 0.05) in the BP group; Corynebacterium (p < 0.05) and Aggregatibacter (p < 0.01) were significantly increased, Veillonella and Lautropia were significantly decreased (p < 0.05) in the CP group; Streptococcus were significantly decreased (p < 0.01) in the DP group; and Leptotrichia were significantly increased (p < 0.05), Porphyromonas (p < 0.05), and Streptococcus (p < 0.01) were significantly decreased in the EP group. Compared with the BP group, the EP group showed a significant increase in Leptotrichia (p < 0.05) and a significant decrease in Fusobacterium (p < 0.05) and Porphyromonas (p < 0.01). Compared to the CP group, the EP group showed a significant increase in Leptotrichia (p < 0.05) and Veillonella (p < 0.01) and a significant decrease in Corynebacterium, Porphyromonas, and Aggregatibacter (p < 0.01) (Figure 4D).
3.6 Genus-level species phylogenetic tree
Representative sequences of the top 100 dominant and abundant genera in the supragingival plaque microorganisms were analyzed, and the microbial communities in all samples were mainly composed of Firmicutes, Proteobacteria, Actinobacteria, Bacteroidetes, and Fusobacteria. Veillonella in Firmicutes; Neisseria in Proteobacteria; Corynebacterium in Actinobacteria; unidentified Prevotellaceae, Capnocytophaga, Prevotella, and Porphyromonas in Bacteroidetes; and Leptotrichia and Fusobacterium in Fusobacteria contributed more to the species composition (Figure 5).
Figure 5
3.7 Network analysis
Network relationship diagrams demonstrate the interaction relationships among microbial genera in different states. The interactions among microorganisms in the five groups of supragingival plaque samples were obtained by calculating the Spearman correlation coefficient for all samples and removing the filtering conditions, which were visualized using a network diagram (Supplementary Figure 2). The connections among the nodes represent a significant correlation among the genera. The AP group contained 64 nodes and 180 connectors, the BP group contained 82 nodes and 264 connectors, the CP group contained 70 nodes and 182 connectors, the DP group contained 98 nodes and 577 connectors, and the EP group contained 78 nodes and 266 connectors. It can be seen that the correlation of genera became progressively more complex in all samples during HDBR and did not return to pre-HDBR levels in the EP group. Furthermore, we analyzed the interactive relationships among the 15 dominant genera (Supplementary Figure 3). There were 17 pairs of interrelationships in the AP group, 13 in the BP group, 16 in the CP group, 5 in the DP group, and 16 in the EP group. With these data, we found that the interrelationship between the dominant genera became simpler during HDBR. Simultaneously, the positive and negative correlations between the dominant genera changed dynamically, and the specific correlation coefficient values are listed in Table 2. No correlation was observed between Porphyromonas, Prevotella, and an unidentified Spirochaetaceae in the AP group, whereas a positive correlation emerged in the DP group. Neisseria showed a positive correlation with Capnocytophaga in both the AP and EP groups; however, this was not observed during the experimental period. Fusobacterium and unidentified Prevotellaceae showed negative correlations in the AP group and positive correlations in the CP group. These results indicated that the correlations among the microbial communities changed under simulated microgravity conditions.
Table 2
| Dominant genera | R | ||||||||
|---|---|---|---|---|---|---|---|---|---|
| AP | BP | CP | DP | EP | |||||
| Leptotrichia | Fusobacterium | −0.67 | – | – | – | – | |||
| Neisseria | −0.66 | – | −0.67 | – | – | ||||
| unidentified_Prevotellaceae | 0.71 | – | – | – | – | ||||
| unidentified_Bacteria | 0.79 | – | – | – | – | ||||
| unidentified_Veillonellaceae | 0.73 | – | – | – | 0.66 | ||||
| Lautropia | −0.66 | – | – | – | – | ||||
| Fusobacterium | Corynebacterium | – | – | −0.69 | – | −0.60 | |||
| unidentified_Prevotellaceae | −0.63 | – | 0.64 | – | – | ||||
| Porphyromonas | – | – | – | – | 0.66 | ||||
| unidentified_Veillonellaceae | −0.61 | – | – | – | – | ||||
| Corynebacterium | Capnocytophaga | – | 0.75 | – | – | – | |||
| Prevotella | – | −0.70 | −0.66 | – | −0.65 | ||||
| Porphyromonas | – | – | −0.62 | – | – | ||||
| unidentified_Spirochaetaceae | – | – | −0.64 | – | – | ||||
| Neisseria | unidentified_Prevotellaceae | −0.84 | −0.78 | −0.66 | – | – | |||
| Capnocytophaga | 0.63 | – | – | – | 0.85 | ||||
| Prevotella | – | – | −0.71 | – | −0.76 | ||||
| Streptococcus | – | 0.69 | – | – | 0.62 | ||||
| Campylobacter | – | −0.60 | – | – | – | ||||
| unidentified_Spirochaetaceae | – | – | – | – | −0.68 | ||||
| unidentified_Veillonellaceae | −0.86 | – | – | – | – | ||||
| Lautropia | 0.72 | – | – | – | – | ||||
| unidentified_Prevotellaceae | Capnocytophaga | – | – | −0.66 | – | – | |||
| Prevotella | – | – | 0.81 | – | – | ||||
| unidentified_Bacteria | 0.68 | – | – | – | – | ||||
| Streptococcus | – | −0.83 | – | – | – | ||||
| Campylobacter | – | 0.80 | – | – | – | ||||
| unidentified_Spirochaetaceae | – | – | 0.65 | – | – | ||||
| unidentified_Veillonellaceae | 0.71 | 0.66 | 0.61 | 0.61 | – | ||||
| Lautropia | −0.65 | – | – | – | – | ||||
| Capnocytophaga | Prevotella | – | – | −0.61 | – | −0.75 | |||
| Aggregatibacter | – | – | 0.62 | – | – | ||||
| unidentified_Spirochaetaceae | – | – | −0.69 | – | −0.60 | ||||
| unidentified_Veillonellaceae | −0.71 | – | – | – | – | ||||
| Lautropia | – | – | – | – | 0.75 | ||||
| Prevotella | Porphyromonas | – | – | – | 0.75 | – | |||
| unidentified_Spirochaetaceae | 0.82 | – | 0.82 | 0.87 | 0.74 | ||||
| Lautropia | – | – | – | – | −0.65 | ||||
| Porphyromonas | unidentified_Spirochaetaceae | – | – | – | 0.76 | – | |||
| Veillonella | Campylobacter | – | 0.67 | – | – | – | |||
| unidentified_Bacteria | Aggregatibacter | – | −0.84 | – | – | −0.64 | |||
| unidentified_Veillonellaceae | 0.62 | 0.61 | – | – | – | ||||
| Streptococcus | Campylobacter | – | −0.73 | – | −0.65 | – | |||
| unidentified_Spirochaetaceae | – | −0.67 | – | – | – | ||||
| unidentified_Veillonellaceae | – | −0.77 | – | – | – | ||||
| Lautropia | 0.78 | – | – | – | – | ||||
| Campylobacter | unidentified_Veillonellaceae | – | 0.68 | – | – | 0.61 | |||
| unidentified_Veillonellaceae | Lautropia | −0.69 | – | – | – | – | |||
Correlations among the dominant genera in five groups.
Spearman’s correlation coefficient is expressed by R value, and the table only shows data with correlation coefficient |R| > 0.6 and p < 0.05 between two dominant genera.
3.8 PICRUSt functional prediction
The Kyoto Encyclopedia of Genes and Genomes (KEGG) database was used to analyze functional genes. Based on the annotation results of the database, the top 10 functional information of each grouping in terms of maximum abundance at each annotation level was selected to generate a functional relative abundance histogram (Figures 6A,B). The main functional gene categories were the same for all groups. According to the functional annotation and abundance information of the samples in the database, the top 35 abundant functions and their abundance information in each sample were selected to draw a heatmap, and clustering was performed based on the functional difference level (Figures 6C,D). Differences between the groups were analyzed based on the annotation results (Figure 6E).
Figure 6
The results of the first annotation hierarchy showed that the KEGG pathways for supragingival plaque microbes were mainly focused on metabolic and genetic information processing. Compared to the AP group, no significant difference was observed in gene function during HDBR; however, the metabolic function was more abundant, and the genetic information processing function decreased significantly in the EP group.
The second annotation hierarchy results showed that the top 10 major functional gene categories in relative abundance were the same for the five groups of samples, including membrane transport, carbohydrate metabolism, amino acid metabolism, replication and repair, and translation. Among them, the most abundant membrane transport function did not differ among the five groups; however, compared with the AP group, carbohydrate metabolism function was richer in the EP group, and amino acid metabolism function was richer in the BP, DP, and EP groups, although replication and repair function and translation function decreased in the EP group.
4 Discussion
Dysbiosis of the oral microbiota can lead to the development of several oral diseases. Therefore, exploring how changes in environmental factors can cause changes in the oral microbiota is crucial.
When clustering OTUs, we found a significant increase in the number of unique OTUs in the supragingival plaque microbial community on day 15 of HDBR compared to pre-HDBR. Meanwhile, Chao1 and the observed species indices showed a significant increase in species richness on day 15 of HDBR. This suggests that the number of species in the microbial community of the supragingival plaque increases under simulated microgravity. Ilyin (2005) collected oral microorganisms from astronauts at the Salyut and Mir orbital stations and found that the number of opportunistic pathogens increased during spaceflight, particularly the massive invasion of Actinomyces_naeslundii, Prevotella_melaninogenica, and Fusobacterium_nucleatum as well as bacterial species capable of maintaining the persistence of inflammation, but declined to preflight levels on day 14 after the return. Urbaniak et al. (2020) used 16S rRNA gene sequencing to analyze salivary microbes in 10 astronauts during missions of 2–9 months and observed that species richness increased significantly during space flight, but declined to preflight levels 6 months after returning to Earth. These studies are consistent with our experimental results, showing that the species richness of microorganisms increases significantly in microgravity environments. However, in this study, species richness did not return to pre-HDBR levels on day 6 of recovery, which may be because we did not observe species richness sufficiently long after HDBR. The Shannon index reflects the diversity of species, i.e., the combination of richness and evenness. The results showed no differential changes in microbial diversity during HDBR compared to pre-HDBR. Species diversity showed a significant decrease on day 6 of recovery compared to that during HDBR; however, considering no differences compared to pre-HDBR, we believe that species diversity returned to pre-HDBR levels on day 6 of recovery. Beta diversity analysis revealed significant differences in the structure of the supragingival plaque microbial community before and during HDBR and in the recovery state. Our group previously used 16S rRNA gene sequencing to examine saliva samples from volunteers under simulated microgravity conditions and found that the structure of the microbial community during HDBR differed from that before HDBR (Sun et al., 2022). However, Sun et al. (2022) did not observe differences in the abundance of salivary microbial communities during HDBR compared to before HDBR. This suggests that the microbial community changes in the free state saliva and adherent state supragingival plaque in the oral cavity are different in microgravity environments. This further demonstrates the necessity of exploring the changes in different oral microniches under simulated microgravity environments.
Species distribution was analyzed among the different groups, and we observed that the relative abundance of supragingival plaque microorganisms differed at different time points. The relative abundances of Corynebacterium and Aggregatibacter increased significantly during HDBR compared to pre-HDBR, and both declined to pre-HDBR levels on day 6 of recovery. One study found that Corynebacterium plays a critical role in plaque formation and underlies plaque structure and community interactions (Schoilew et al., 2019). Metagenomic sequence analysis showed that Corynebacterium is a key taxon of supragingival plaques, constituting a hedgehog-like structure and providing a structural backbone for the aggregation of other bacteria (e.g., Streptococcus and Porphyromonas), thereby forming a polymicrobial complex (Mark Welch et al., 2016; Fakhruddin et al., 2019). Chen et al. (2021) collected oropharyngeal samples from four healthy volunteers in a 180-day ground-simulated space-enclosed environment and found Corynebacterium as the dominant genus, whose relative abundance increased significantly throughout the experiment. In this study, a significant increase in the relative abundance of Corynebacterium during HDBR may have facilitated the aggregation of more pathogenic microorganisms. Aggregatibacter is gram-negative bacilli that is integral to healthy oral microbial communities; however, in conditions such as trauma or mucosal damage, it invades healthy tissues, causing infection. Members of this genus identified in our analysis were Aggregatibacter segnis and Aggregatibacter actinomycetemcomitans. A. actinomycetemcomitans is a well-known human periodontal pathogen (Fine et al., 2019). A. actinomycetemcomitans can produce various virulence factors, such as leukotoxins and cell expansion lethal toxins, destroy the local immunity of the host, and help periodontal pathogens evade host immunity, leading to the occurrence and development of periodontal disease (Shenker et al., 1999; Kachlany, 2010). Simultaneously, A. segnis has been suggested to have adhesiveness to the cellular structure and coexists with other bacteria, causing supragingival calculus formation (Baris et al., 2017).
In contrast, the relative abundances of Veillonella, Lautropia, and Streptococcus decreased significantly during HDBR compared to pre-HDBR. On day 6 of recovery, the relative abundances of Veillonella and Lautropia increased to pre-HDBR levels, whereas the relative abundance of Streptococcus remained different from that of pre-HDBR. Evidence exists that Veillonella lacks glycolytic capacity, relies solely on lactic acid and other organic acids as energy sources, produces lower pKa acids, and may increase the pH, suggesting that Veillonella may have beneficial effects on the oral environment (Gross et al., 2010). Quantitative data suggest a reduction in the frequency index of isolation of Veillonella gingival crevicular fluid of astronauts on day 1 after space flight (Volozhin et al., 2001). Similarly, Lautropia is a normal oral cavity colonizer and is the dominant genus in healthy oral environments (Ikeda et al., 2020; Fragkioudakis et al., 2021). Individuals with periodontal disease have been observed to exhibit a decrease in the content of Lautropia (Shaddox et al., 2012). Streptococcus species consists of numerous opportunistic pathogens associated with dental caries. Cheng et al. (2014) found that simulating microgravity conditions can change the structure of Streptococcus mutans biofilms and the distribution of extracellular polysaccharides and improve the acid resistance of S. mutans. Brown et al. (1974) detected a significant increase in S. mutans counts in the dental plaque of three astronauts during a 56-day simulated space laboratory mission. Because of limited experimental conditions, they only examined quantitative changes in species.
Furthermore, a significant increase was observed in the relative abundance of Leptotrichia and a significant decrease in the relative abundance of Porphyromonas on day 6 of recovery compared to pre-HDBR and during HDBR. Leptotrichia is highly glycolytic and ferments a variety of monosaccharides and disaccharides into lactic acid (Birkeland and Hofstad, 1985; Thompson and Pikis, 2012). This implies that Leptotrichia may be involved in the development of dental caries and may have cariogenic potential. As an opportunistic pathogen, Leptotrichia is associated with infections, particularly in immunocompromised hosts (Eribe and Olsen, 2008; Eribe and Olsen, 2017). Moreover, Brown et al. (1976) found that 4 days after the end of a 59-day space laboratory mission, the number of Leptotrichia in the plaque of major aircrew increased significantly compared with that before the flight, which is consistent with our findings. Porphyromonas are genetically a highly heterogeneous group, and some species in this group, such as P. pasteri and P. catoniae, can be used as potential oral health markers (Crielaard et al., 2011; Yasunaga et al., 2017). However, other species, such as P. gingivalis and P. endodontalis, are closely related to the occurrence and development of oral diseases (Hajishengallis et al., 2012). The primary colonizing habitat of P. gingivalis in a normal oral environment is the sulcus (Reinhardt et al., 1989). However, the flow of crevicular fluid in microgravity environments may lead to the translocation of microorganisms. Shi et al. (2016) found through a rotating cell culture system that simulated microgravity promotes the growth of P. gingivalis and alters its gene expression.
Network analysis showed that the interaction relationships between the 15 dominant genera in the supragingival plaque microbial community became simpler during HDBR than during pre-HDBR; however, new correlations emerged. We observed no correlation between Porphyromonas and Prevotella pre-HDBR or between Porphyromonas and unidentified Spirochaetaceae; however, a positive correlation was observed on day 15 of HDBR. The co-aggregation between gingivalis and Prevotella_oris has been proposed as a factor that promotes the advancement of periodontitis (Sato and Nakazawa, 2014). Simultaneously, P. intermedia exerted a synergistic effect with P. gingivalis W83 and exacerbated the virulence of P. gingivalis by regulating the expression of its virulence factors (Zhang et al., 2019). Moreover, a synergistic growth relationship was observed between P. gingivalis and T. denticola (Huang et al., 2011). P. gingivalis produces isobutyric acid, which promotes the growth of T. denticola (Grenier, 1992). T. denticola produces succinic acid found in the phospholipids and lipids on the cell envelope of P. gingivalis (Lev, 1979). Therefore, we believe that alterations in the interactions between microbes may disrupt the balance of the oral microbial community structure and increase the risk of opportunistic infections. The impact of the interactions between these microbes should not be underestimated.
We also explored the alterations in microbial gene function in supragingival plaques under simulated microgravity. We found that amino acid metabolism was more abundant in the supragingival plaque microbes during HDBR, and this alteration persisted until day 6 of recovery. Current studies suggest that a complex regulatory network exists for amino acid metabolism and that acidic, basic, aromatic, and sulfur-containing amino acid metabolites produced by oral microorganisms through amino acid metabolism may negatively affect the host and exert pro-inflammatory and cytotoxic effects (Barbour et al., 2022). Moreover, environmental factors can alter the metabolic properties of certain bacteria, such as Porphyromonas gingivalis, Prevotella intermedia, and Veillonella, allowing them to exert different pathogenicity (Li et al., 2022). Therefore, we believe that the metabolic balance between microorganisms and their hosts plays an important role in maintaining a stable state of the body. The gene function of supragingival plaque microorganisms changed under the influence of a simulated microgravity environment, and this change persisted during the recovery period, which may lead to disturbances in oral microecology.
Overall, this study combined high-throughput sequencing technology with simulated microgravity experiments to comprehensively analyze the gene sequences of supragingival plaque microorganisms. We chose multiple time points to investigate changes in the microbial community under simulated microgravity conditions. Moreover, we explored the functions of the microbial genome to provide new insights for future studies. In the future, we will combine microbiome changes with laboratory findings to conduct metabolomic studies of supragingival plaque microorganisms in a broader population, further establish disease metabolic regulatory networks, lay the foundation for establishing in vivo and ex vivo research targets for the effects of microgravity environments on oral supragingival plaque microorganisms, and deeply explore the relationship between oral dysbiosis and immune disorders and systemic changes in the microgravity environment, providing a theoretical basis for the prevention and treatment of oral aerospace diseases.
5 Conclusion
This experimental study showed that a simulated microgravity environment with a 15-day −6° head-down bed rest would have an effect on the supragingival plaque microbial community in young men. In this state, the richness of the plaque microbial communities was significantly altered, and the community composition differed between the groups. Simultaneously, the relative abundance and functional genes of the microorganisms changed. Some of these changes persisted until the end of the experiment. This suggests that the simulated microgravity environment may affect the oral ecosystem by changing the balance of the supragingival plaque microbial community, leading to infectious diseases. Astronauts need to focus on oral hygiene care during and after spaceflight, either through a balanced diet or supplementation with specific probiotics to maintain a balanced microbial community.
Statements
Data availability statement
The original contributions presented in the study are publicly available. This data can be found here: [https://www.ncbi.nlm.nih.gov/, PRJNA1003027].
Ethics statement
The studies involving humans were approved by the Ethics Review Committee of the Space Science and Technology Institute (Shenzhen). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
DZ: Formal analysis, Validation, Writing – original draft, Writing – review & editing. PQ: Conceptualization, Investigation, Writing – review & editing. QZ: Conceptualization, Investigation. HS: Formal analysis. BX: Methodology, Resources. BW: Methodology, Resources. CT: Conceptualization, Project administration, Supervision, Writing – review & editing.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Acknowledgments
We sincerely thank all volunteers who provided samples. We thank the Space Science and Technology Institute (Shenzhen) for providing the experimental conditions. We thank Novogene Company Limited (Beijing) for sequencing the data from this study. All members of our team are thanked for their efforts to successfully complete this study.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2024.1331023/full#supplementary-material
Footnotes
References
1
AarabiG.HeydeckeG.SeedorfU. (2018). Roles of Oral infections in the Pathomechanism of atherosclerosis. Int. J. Mol. Sci.19:1978. doi: 10.3390/ijms19071978
2
Al-HebshiN. N.BaraniyaD.ChenT.HillJ.PuriS.TellezM.et al. (2019). Metagenome sequencing-based strain-level and functional characterization of supragingival microbiome associated with dental caries in children. J. Oral Microbiol.11:1557986. doi: 10.1080/20002297.2018.1557986
3
AndersonM. J. (2001). A new method for non-parametric multivariate analysis of variance. Austral Ecol.26, 32–46. doi: 10.1111/j.1442-9993.2001.01070.pp.x
4
AtarashiK.SudaW.LuoC.KawaguchiT.MotooI.NarushimaS.et al. (2017). Ectopic colonization of oral bacteria in the intestine drives TH1 cell induction and inflammation. Science358, 359–365. doi: 10.1126/science.aan4526
5
BallJ. R.EvansC. H. (2001). Safe passage: astronaut care for exploration missions. Washington, DC: National Academies Press.
6
BarbourA.ElebyaryO.FineN.OveisiM.GlogauerM. (2022). Metabolites of the oral microbiome: important mediators of multikingdom interactions. FEMS Microbiol. Rev.46:fuab039. doi: 10.1093/femsre/fuab039
7
BarisO.DemirT.GulluceM. (2017). Investigation of in vitro mineral forming bacterial isolates from supragingival calculus. Niger. J. Clin. Pract.20, 1571–1575. doi: 10.4103/1119-3077.187316
8
BirkelandN. K.HofstadT. (1985). Oligosaccharides obtained by partial hydrolysis of lipopolysaccharides from Leptotrichia buccalis. Scand. J. Dent. Res.93, 432–435. doi: 10.1111/j.1600-0722.1985.tb01335.x
9
BokulichN. A.SubramanianS.FaithJ. J.GeversD.GordonJ. I.KnightR.et al. (2013). Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nat. Methods10, 57–59. doi: 10.1038/nmeth.2276
10
BowenW. H. (2016). Dental caries – not just holes in teeth! A perspective. Mol Oral Microbiol31, 228–233. doi: 10.1111/omi.12132
11
BrownL. R.FromeW. J.HandlerS. F.WheatcroftM. G.RiderL. J. (1977). Skylab oral health studies. Washington, DC: NASA.
12
BrownL. R.FrommeW. J.HandlerS. F.WheatcroftM. G.JohnstonD. A. (1976). Effect of Skylab missions on clinical and microbiologic aspects of oral health. J. Am. Dent. Assoc.93, 357–363. doi: 10.14219/jada.archive.1976.0502
13
BrownL. R.WheatcroftM. G.FromeW. J.RiderL. J. (1974). Effects of a simulated Skylab mission on the oral health of astronauts. J. Dent. Res.53, 1268–1275. doi: 10.1177/00220345740530053201
14
CaporasoJ. G.KuczynskiJ.StombaughJ.BittingerK.BushmanF. D.CostelloE. K.et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods7, 335–336. doi: 10.1038/nmeth.f.303
15
ChapmanM. G.UnderwoodA. J. (1999). Ecological patterns in multivariate assemblages: information and interpretation of negative values in ANOSIM tests. Mar. Ecol. Prog. Ser.180, 257–265. doi: 10.3354/meps180257
16
ChenY.XuC.ZhongC.LyuZ.LiuJ.ChenZ.et al. (2021). Temporal characteristics of the oropharyngeal and nasal microbiota structure in crewmembers stayed 180 days in the controlled ecological life support system. Front. Microbiol.11:617696. doi: 10.3389/fmicb.2020.617696
17
ChengX.XuX.ChenJ.ZhouX.ChengL.LiM.et al. (2014). Effects of simulated microgravity on Streptococcus mutans physiology and biofilm structure. FEMS Microbiol. Lett.359, 94–101. doi: 10.1111/1574-6968.12573
18
CrielaardW.ZauraE.SchullerA. A.HuseS. M.MontijnR. C.KeijserB. J. (2011). Exploring the oral microbiota of children at various developmental stages of their dentition in the relation to their oral health. BMC Med. Genet.4:22. doi: 10.1186/1755-8794-4-22
19
CrucianB.StoweR. P.MehtaS.QuiriarteH.PiersonD.SamsC. (2015). Alterations in adaptive immunity persist during long-duration spaceflight. NPJ Microgravity1:15013. doi: 10.1038/npjmgrav.2015.13
20
DarveauR. P.TannerA.PageR. C. (1997). The microbial challenge in periodontitis. Periodontol2000, 12–32. doi: 10.1111/j.1600-0757.1997.tb00190.x
21
EdgarR. C. (2004). MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res.32, 1792–1797. doi: 10.1093/nar/gkh340
22
EdgarR. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods10, 996–998. doi: 10.1038/nmeth.2604
23
EribeE. R.OlsenI. (2008). Leptotrichia species in human infections. Anaerobe14, 131–137. doi: 10.1016/j.anaerobe.2008.04.004
24
EribeE. R. K.OlsenI. (2017). Leptotrichia species in human infections II. J. Oral Microbiol.9:1368848. doi: 10.1080/20002297.2017.1368848
25
FakhruddinK. S.NgoH. C.SamaranayakeL. P. (2019). Cariogenic microbiome and microbiota of the early primary dentition: a contemporary overview. Oral Dis.25, 982–995. doi: 10.1111/odi.12932
26
FineD. H.PatilA. G.VelusamyS. K. (2019). Aggregatibacter actinomycetemcomitans (Aa) under the radar: myths and misunderstandings of aa and its role in aggressive periodontitis. Front. Immunol.10:728. doi: 10.3389/fimmu.2019.00728
27
FosterJ. S.KhodadadC. L.AhrendtS. R.ParrishM. L. (2013). Impact of simulated microgravity on the normal developmental time line of an animal-bacteria symbiosis. Sci. Rep.3:1340. doi: 10.1038/srep01340
28
FragkioudakisI.RiggioM. P.ApatzidouD. A. (2021). Understanding the microbial components of periodontal diseases and periodontal treatment-induced microbiological shifts. J. Med. Microbiol.70. doi: 10.1099/jmm.0.001247
29
GontcharovI. B.KovachevichI. V.PoolS. L.NavinkovO. L.BarrattM. R.BogomolovV. V.et al. (2005). In-flight medical incidents in the NASA-Mir program. Aviat. Space Environ. Med.76, 692–696. PMID:
30
GraebeA.SchuckE. L.LensingP.PutchaL.DerendorfH. (2004). Physiological, pharmacokinetic, and pharmacodynamic changes in space. J. Clin. Pharmacol.44, 837–853. doi: 10.1177/0091270004267193
31
GrenierD. (1992). Nutritional interactions between two suspected periodontopathogens, Treponema denticola and Porphyromonas gingivalis. Infect. Immun.60, 5298–5301. doi: 10.1128/iai.60.12.5298-5301.1992
32
GrenonS. M.SaaryJ.GrayG.VanderploegJ. M.Hughes-FulfordM. (2012). Can I take a space flight? Considerations for doctors. BMJ345:e8124. doi: 10.1136/bmj.e8124
33
GrimmD.GrosseJ.WehlandM.MannV.ReselandJ. E.SundaresanA.et al. (2016). The impact of microgravity on bone in humans. Bone87, 44–56. doi: 10.1016/j.bone.2015.12.057
34
GrossE. L.LeysE. J.GasparovichS. R.FirestoneN. D.SchwartzbaumJ. A.JaniesD. A.et al. (2010). Bacterial 16S sequence analysis of severe caries in young permanent teeth. J. Clin. Microbiol.48, 4121–4128. doi: 10.1128/JCM.01232-10
35
HaasB. J.GeversD.EarlA. M.FeldgardenM.WardD. V.GiannoukosG.et al. (2011). Chimeric 16S rRNA sequence formation and detection in sanger and 454-pyrosequenced PCR amplicons. Genome Res.21, 494–504. doi: 10.1101/gr.112730.110
36
HajishengallisG.DarveauR. P.CurtisM. A. (2012). The keystone-pathogen hypothesis. Nat. Rev. Microbiol.10, 717–725. doi: 10.1038/nrmicro2873
37
HargensA. R.VicoL. (2016). Long-duration bed rest as an analog to microgravity. J. Appl. Physiol.120, 891–903. doi: 10.1152/japplphysiol.00935.2015
38
HuangR.LiM.GregoryR. L. (2011). Bacterial interactions in dental biofilm. Virulence2, 435–444. doi: 10.4161/viru.2.5.16140
39
HuangB.LiD. G.HuangY.LiuC. T. (2018). Effects of spaceflight and simulated microgravity on microbial growth and secondary metabolism. Mil. Med. Res.5:18. doi: 10.1186/s40779-018-0162-9
40
IkedaE.ShibaT.IkedaY.SudaW.NakasatoA.TakeuchiY.et al. (2020). Japanese subgingival microbiota in health vs disease and their roles in predicted functions associated with periodontitis. Odontology108, 280–291. doi: 10.1007/s10266-019-00452-4
41
IlyinV. K. (2005). Microbiological status of cosmonauts during orbital spaceflights on Salyut and Mir orbital stations. Acta Astronaut.56, 839–850. doi: 10.1016/j.actaastro.2005.01.009
42
JagathrakshakanS. N.SethumadhavaR. J.MehtaD. T.RamanathanA. (2015). 16S rRNA gene-based metagenomic analysis identifies a novel bacterial co-prevalence pattern in dental caries. Eur. J. Dent.9, 127–132. doi: 10.4103/1305-7456.149661
43
JohnstonR. S.DietleinL. F. (1977). Biomedical results from Skylab. Washington, DC: NASA.
44
KachlanyS. C. (2010). Aggregatibacter actinomycetemcomitans leukotoxin: from threat to therapy. J. Dent. Res.89, 561–570. doi: 10.1177/0022034510363682
45
KlausD. M.HowardH. N. (2006). Antibiotic efficacy and microbial virulence during space flight. Trends Biotechnol.24, 131–136. doi: 10.1016/j.tibtech.2006.01.008
46
KumarP. S.LeysE. J.BrykJ. M.MartinezF. J.MoeschbergerM. L.GriffenA. L. (2006). Changes in periodontal health status are associated with bacterial community shifts as assessed by quantitative 16S cloning and sequencing. J. Clin. Microbiol.44, 3665–3673. doi: 10.1128/JCM.00317-06
47
KwunJ.-S.KangS.-H.LeeH.-J.ParkH.-K.LeeW.-J.YoonC.-H.et al. (2020). Comparison of thrombus, gut, and oral microbiomes in Korean patients with ST-elevation myocardial infarction: a case-control study. Exp. Mol. Med.52, 2069–2079. doi: 10.1038/s12276-020-00543-1
48
LevM. (1979). Sphingolipid biosynthesis and vitamin K metabolism in Bacteroides melaninogenicus. Am. J. Clin. Nutr.32, 179–186. doi: 10.1093/ajcn/32.1.179
49
LiY. C.ZhangS. W.PanY. P. (2022). Effects of oral bacterial metabolites on systemic diseases. Zhonghua Kou Qiang Yi Xue Za Zhi57, 648–653. doi: 10.3760/cma.j.cn112144-20220325-00132
50
LozuponeC. A.HamadyM.KelleyS. T.KnightR. (2007). Quantitative and qualitative beta diversity measures lead to different insights into factors that structure microbial communities. Appl. Environ. Microbiol.73, 1576–1585. doi: 10.1128/AEM.01996-06
51
LozuponeC.LladserM. E.KnightsD.StombaughJ.KnightR. (2011). UniFrac: an effective distance metric for microbial community comparison. ISME J.5, 169–172. doi: 10.1038/ismej.2010.133
52
MagočT.SalzbergS. L. (2011). FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics27, 2957–2963. doi: 10.1093/bioinformatics/btr507
53
Mark WelchJ. L.RossettiB. J.RiekenC. W.DewhirstF. E.BorisyG. G. (2016). Biogeography of a human oral microbiome at the micron scale. Proc. Natl. Acad. Sci. U. S. A.113, E791–E800. doi: 10.1073/pnas.1522149113
54
MarshP. D. (2003). Are dental diseases examples of ecological catastrophes?Microbiology149, 279–294. doi: 10.1099/mic.0.26082-0
55
MarshP. D. (2012). Contemporary perspective on plaque control. Br. Dent. J.212, 601–606. doi: 10.1038/sj.bdj.2012.524
56
Pavy-Le TraonA.HeerM.NariciM. V.RittwegerJ.VernikosJ. (2007). From space to earth: advances in human physiology from 20 years of bed rest studies (1986-2006). Eur. J. Appl. Physiol.101, 143–194. doi: 10.1007/s00421-007-0474-z
57
PoncheewinW.HermesG. D. A.van DamJ. C. J.KoehorstJ. J.SmidtH.SchaapP. J. (2019). NG-tax 2.0: a semantic framework for high-throughput amplicon analysis. Front. Genet.10:1366. doi: 10.3389/fgene.2019.01366
58
QuastC.PruesseE.YilmazP.GerkenJ.SchweerT.YarzaP.et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res.41, D590–D596. doi: 10.1093/nar/gks1219
59
RaiD. B. (2009). Human oral cavity in simulated microgravity: new prospects. Adv. Med. Dent. Sci.3, 35–39.
60
ReinhardtR. A.McDonaldT. L.BoltonR. W.DuBoisL. M.KaldahlW. B. (1989). IgG subclasses in gingival crevicular fluid from active versus stable periodontal sites. J. Periodontol.60, 44–50. doi: 10.1902/jop.1989.60.1.44
61
RoewerL.KayserM.DieltjesP.NagyM.BakkerE.KrawczakM.et al. (1996). Analysis of molecular variance (AMOVA) of Y-chromosome-specific microsatellites in two closely related human populations. Hum. Mol. Genet.5, 1029–1033. doi: 10.1093/hmg/5.7.1029
62
RognesT.FlouriT.NicholsB.QuinceC.MahéF. (2016). VSEARCH: a versatile open source tool for metagenomics. PeerJ4:e2584. doi: 10.7717/peerj.2584
63
SatoT.NakazawaF. (2014). Coaggregation between Prevotella oris and Porphyromonas gingivalis. J. Microbiol. Immunol. Infect.47, 182–186. doi: 10.1016/j.jmii.2012.09.005
64
SchoilewK.UeffingH.DalpkeA.WolffB.FreseC.WolffD.et al. (2019). Bacterial biofilm composition in healthy subjects with and without caries experience. J. Oral Microbiol.11:1633194. doi: 10.1080/20002297.2019.1633194
65
ShaddoxL. M.HuangH.LinT.HouW.HarrisonP. L.AukhilI.et al. (2012). Microbiological characterization in children with aggressive periodontitis. J. Dent. Res.91, 927–933. doi: 10.1177/0022034512456039
66
ShenkerB. J.McKayT.DatarS.MillerM.ChowhanR.DemuthD. (1999). Actinobacillus actinomycetemcomitans immunosuppressive protein is a member of the family of cytolethal distending toxins capable of causing a G2 arrest in human T cells. J. Immunol.162, 4773–4780. doi: 10.4049/jimmunol.162.8.4773
67
ShiX.LuY.TangC.NiuZ. (2016). Effect of simulated microgravity on Porphyromonas gingivalis strain ATCC33277 growth and gene expression. Int. J. Clin. Exp. Pathol.9, 3143–3151.
68
ShiW.TianJ.XuH.ZhouQ.QinM. (2018). Distinctions and associations between the microbiota of saliva and supragingival plaque of permanent and deciduous teeth. PLoS One13:e0200337. doi: 10.1371/journal.pone.0200337
69
SunH.ZhouQ.QiaoP.ZhuD.XinB.WuB.et al. (2022). Short-term head-down bed rest microgravity simulation alters salivary microbiome in young healthy men. Front. Microbiol.13:1056637. doi: 10.3389/fmicb.2022.1056637
70
SundbladP.OrlovO.AngererO.LarinaI.CromwellR. (2016). Standardization of bed rest studies in the spaceflight context. J. Appl. Physiol.121, 348–349. doi: 10.1152/japplphysiol.00089.2016
71
ThompsonJ.PikisA. (2012). Metabolism of sugars by genetically diverse species of oral Leptotrichia. Mol Oral Microbiol27, 34–44. doi: 10.1111/j.2041-1014.2011.00627.x
72
UrbaniakC.LorenziH.ThissenJ.JaingC.CrucianB.SamsC.et al. (2020). The influence of spaceflight on the astronaut salivary microbiome and the search for a microbiome biomarker for viral reactivation. Microbiome8:56. doi: 10.1186/s40168-020-00830-z
73
VolozhinA. I.TsarevV. N.MalnevaN. S.SashkinaT. I.SaldusovaI. V. (2001). Interaction peculiarities between microbial cenosis and local immunity of periodontium of humans under extreme conditions. Acta Astronaut.49, 53–57. doi: 10.1016/s0094-5765(00)00128-4
74
WhiteR. J.AvernerM. (2001). Humans in space. Nature409, 1115–1118. doi: 10.1038/35059243
75
WilsonJ. W.OttC. M.zu BentrupK. H.RamamurthyR.QuickL.PorwollikS.et al. (2007). Space flight alters bacterial gene expression and virulence and reveals a role for global regulator Hfq. Proc. Natl. Acad. Sci. U. S. A.104, 16299–16304. doi: 10.1073/pnas.0707155104
76
YangX.HeL.YanS.ChenX.QueG. (2021). The impact of caries status on supragingival plaque and salivary microbiome in children with mixed dentition: a cross-sectional survey. BMC Oral Health21:319. doi: 10.1186/s12903-021-01683-0
77
YasunagaH.TakeshitaT.ShibataY.FurutaM.ShimazakiY.AkifusaS.et al. (2017). Exploration of bacterial species associated with the salivary microbiome of individuals with a low susceptibility to dental caries. Clin. Oral Investig.21, 2399–2406. doi: 10.1007/s00784-016-2035-5
78
YoshiokaH.YoshimuraA.KanekoT.GolenbockD. T.HaraY. (2008). Analysis of the activity to induce toll-like receptor (TLR)2- and TLR4-mediated stimulation of supragingival plaque. J. Periodontol.79, 920–928. doi: 10.1902/jop.2008.070516
79
ZhangY.ShiW.SongY.WangJ. (2019). Metatranscriptomic analysis of an in vitro biofilm model reveals strain-specific interactions among multiple bacterial species. J. Oral Microbiol.11:1599670. doi: 10.1080/20002297.2019.1599670
80
ZhangS.YuN.ArceR. M. (2020). Periodontal inflammation: integrating genes and dysbiosis. Periodontol.82, 129–142. doi: 10.1111/prd.12267
Summary
Keywords
head-down bed rest, microgravity, supragingival plaque, microbial community, 16S rRNA high-throughput sequencing
Citation
Zhu D, Qiao P, Zhou Q, Sun H, Xin B, Wu B and Tang C (2024) Effect of 15 days −6° head-down bed rest on microbial communities of supragingival plaque in young men. Front. Microbiol. 15:1331023. doi: 10.3389/fmicb.2024.1331023
Received
31 October 2023
Accepted
04 January 2024
Published
24 January 2024
Volume
15 - 2024
Edited by
Zhenlin Han, University of Hawaii at Manoa, United States
Reviewed by
Rekha Arya, Sungkyunkwan University, Republic of Korea
Sahana Vasudevan, Institute for Stem Cell Science and Regenerative Medicine (inStem), India
Updates
Copyright
© 2024 Zhu, Qiao, Zhou, Sun, Xin, Wu and Tang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Chuhua Tang, tch306@126.com
†These authors have contributed equally to this work and share first authorship
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.