Results
To investigate the difference in bacterial community of S-ECC among the original dental plaque, and cultured bacteria, we recruited nine S-ECC patients, half of the samples were used for culture and the other half for 16S rRNA gene high-throughput sequencing directly without culture. Alpha diversity indices of ACE, Chao1, Shannon and Simpson showed clinical dental plaque samples had greater alpha diversity than either all their cultured counterparts ( Figure 1A ) or the cultured bacteria under normoxia and anaerobic conditions ( Figure 1B ), although the Venn diagram showed that the OTU numbers of cultured bacteria at species level under anaerobic and normoxia conditions were higher than that in clinical samples ( Figure 2A ). Core flora of clinical S-ECC samples had 14 common features, including two unclassified Streptococcus , one unclassified Granulicatella , two unclassified Veillonella , one unclassified Fusobacterium , one unclassified Aggregatibacter , one unclassified Kingella , one unclassified Lachnoanaerobaculum , one unclassified Lactobacillales , one unclassified Abiotrophia , one unclassified Rothia , one Rothia aeria , and one Capnocytophaga endodontalis ( Fig. S2 ).

To visualize the distribution of samples in coordinate system, NMDS analysis showed that anaerobic cultured bacteria and normoxia cultured bacteria samples had similar culture features. However, both had a markedly different bacterial structure in comparison with clinical sample group ( Figure 2B ).

To display the composition of dominant bacteria, the stacked column plot showed the percentage of each high-abundance bacteria. In clinical samples, except Veillonella , Leptotrichia and Streptococcus , other bacteria occupied the most of bacterial community. No matter under anaerobic and normoxia conditions, the percentages of other bacteria, Veillonella and Leptotrichia were clearly decreased. Meanwhile, there were observed increases in Desemzia , Actietobacter , Lactobacillus , Limosilactobacillus , Ligilactobacillus , Vagococcus and Streptococcus ( Figure 3A ). The results indicated that there were more bacterial categories in clinical samples than cultures by media. NEMM and SHI media only could cultivate part of clinical bacteria in S-ECC plaques. There is still a lot of room for improvement in the culture medium.

To detect the difference among the groups, ANOVA was used to display the relative difference of bacteria at the genus level. The results showed that both anaerobic and normoxia cultures were beneficial for the growth of Streptococcus , Bacillus (including other unclassified Bacillus ) and Bifidobacterium . The relative higher abundance of most other bacteria, including Leptotrichia , Veillonella , Corynebacterium , Actinomyces , Fusobacterium , Neisseria , Campocytophaga , Prevotella , Selenomonas , Campylobacter and Aggregatibacter in clinical samples presented marked lower percentage in both anaerobic and normoxia cultures. It suggested that most of the above-mentioned bacteria categories could not be successfully enriched by NEMM and SHI media ( Figure 3B ).

To assess the predicted function of the bacterial community in clinical samples and the cultured bacteria samples, microbial phenotype prediction was evaluated by BugBase followed by the Mann–Whitney U test. The results showed that there were decreases in Gram-positive bacteria, facultatively anaerobic bacteria and stress tolerance ability, but increases in Gram-negative bacteria, anaerobic respiration, biofilms forming and potential pathogenicity in clinical samples compared with cultured samples under aerobic and anaerobic conditions ( Figure 4 ). The results indicated that uncultured Gram-negative, anaerobic, potentially pathogenic and biofilm-forming bacteria may contribute to S-ECC development.

All cultured microbiota were examined in terms of both species proportion and abundance, with a specific focus on the enrichment of species. Of the organisms presented in the inoculum, all 32 phyla, 87 classes, 260 orders, 501 families, 1,107 genera and 1,502 bacterial species were captured. We identified 3,719, 5,605, 3,705, and 3,460 OTUs in the aNEMM, nNEMM, nSHI and aSHI groups, respectively. Among them, 608 OTUs were uniform, occupying 4.9% of all the OTUs detected, indicating that a steady facultative anaerobe microorganism composition was found in all four fluid media. The other OTUs were not shared in all groups. Among these variable microbiomes, 2,061 OTUs were unique to the aNEMM group, while 1,848 OTUs could be detected only in aSHI. Moreover, 4,010 OTUs could be found only in nNEMM groups, while 2,096 OTUs could be found only in nSHI groups.

To compare NEMM with SHI, which are suitable for cultivating specific bacterial types, we further analyzed 16S rRNA gene high-throughput sequencing data. Alpha diversity indices of ACE and Chao1 showed no significant difference between NEMM and SHI under the normoxia condition. However, under anaerobic culture conditions, NEMM had greater alpha diversity indices of ACE and Chao1 than SHI ( Figure 5A ). NMDS analysis showed that regardless of whether anaerobic or normoxia, cultured bacteria samples had similar features ( Figure 5B ). Venn diagram showed that the OTU numbers of NEMM cultured bacteria under anaerobic and normoxia conditions were higher than that in SHI cultured bacteria under the same conditions ( Figure 5C ).

Differential predominant bacteria between NEMM and SHI were analyzed in the relative abundance of bacteria at the genus level by ANOVA. Of the 19 most predominant bacterial genera found in NEMM and SHI media, Ligilactobacillus , Limosilactobacillus , Lactobacillus , and Desemzia showed the highest abundance levels. Based on the relative abundance of cultured bacteria, NEMM was suitable for cultivating bacteria including Lactobacillus, Desemzia , Veillonella , Psychrobacterium and Leptotrichia , while SHI was good for culturing Ligilactobacillus , Limosilactobacillus , Bacillus , Neisseria , Leuconostoc , Cryptobacterium , Aggregatibacter and Kingella ( Figure 6 ).

Bugbase was used to predict the function of bacterial community in NEMM and SHI under anaerobic and normoxia culture conditions. The results showed that the microbiota developed in NEMM were predicted stronger anaerobic ability, whether under anaerobic conditions or normoxic conditions. Under normoxic culture conditions, there were no significant differences in Gram-positive bacteria among the four groups (aNEMM, aSHI, nNEMM and nSHI). However, under anaerobic conditions, NEMM supported the growth of more Gram-negative bacteria and less Gram-positive bacteria. Dental plaque samples cultured in NEMM were more likely to be stress tolerant and less likely to be potential pathogens and form biofilms. Moreover, it had weaker facultatively anaerobic ability when under anaerobic conditions ( Figure 7 ).

To assess the effects of the two media on the predicted gene categories (COGs), we compared the predicted COGs between aNEMM and aSHI, nNEMM and nSHI ( Figure 8A ). NEMM showed more ability in carbohydrate transport and metabolism than SHI medium under both anaerobic conditions and normoxia conditions. Additionally, those in NEMM displayed an increased relative abundance of protein functions linked to defense metabolism, inorganic ion transport and metabolism under anaerobic conditions. Furthermore, pathways of membrane transport and cellular community-prokaryotes were associated with aNEMM group ( Figure 8B ).

The predominant bacteria were largely consistent among the four groups (aNEMM, aSHI, nNEMM, nSHI), but different relative abundances could be observed ( Fig. S3 ). More genera, species and uncultured bacteria were captured in aNEMM groups at day 1 ( Table S2 ) under anaerobical condition compared with those in aSHI Day1 group. For normoxia conditions, there was a significant difference of genera and uncultured bacteria between NEMM and SHI media ( Table 1 ). Table S3 showed that genera, species and uncultured bacteria detected in nNEMM were higher compared to nSHI at each time point.

Notes:

The nonparametric Wilcoxon test was performed to evaluate differences.

SD standard deviation

standard deviation

Principal component analysis of the development of the oral microbiome shows distinct shifts in the composition of the microbiome. We could see the changes in bacterial community structure over time in Figs. 9A and 9B . The proportions of Streptococcus were reduced along with time in the four groups. In this study, the NEMM as a growth and colonization substrate resulted in a dynamic microbial succession, characterized by an initial proportional reduction in Veillonella , corresponding to an increase of Lactobacillus and unclassified_Bacilli proportions under anaerobic conditions. Bifidobacterium reached the highest percentage on the 15th day and then decreased in aNEMM group. In the aSHI group, the proportions of Lactobacillus and unclassified_Bacilli increased gradually before the 15th day. While Ligilactobacillus increased rapidly from day one to day three in both aSHI and nSHI groups. In nNEMM group, Ligilactobacillus , unclassified_Bacilli , Lactobacillus , Lacticaseibacillus and Lactiplantibacillus showed the time-dependent accumulation, while Acinetobacter exhibited the opposite tendency.

The taxa that most likely explains the differences between NEMM and SHI media from different time points are defined by LEfSe. Figures 10A and 10B showed cladograms representing the potential biomarkers of different groups. At the genus level, Veillonella was significantly enrichment in the aNEMM cultured Day1 group (aNEMM-Day1), while Streptococcus , Gardiobacterium and Granulicatella exhibited relatively higher abundance in the aSHI-Day1 group. Bifidobacterium , Akkermansia and Desemzia were significantly more abundant in the aNEMM-Day15, aNEMM-Day21 and aNEMM-Day28 groups, respectively. While Ligilactobacillus , unclassified_Bacilli and Lacticaseibacillus were most abundant in aSHI-Day6, aSHI-Day21 and aSHI-Day28 groups, respectively. For samples grown aerobically, significant bacterial differences were also detected, with Veillonella , Desemzia and Lacticaseibacillus were remarkably enriched in the nNEMM-Day1, nNEMM-Day9 and nNEMM-Day28 groups at the genus level, while Streptococcus exhibited higher abundance in the nSHI-Day1 group ( Figs. 11A and 11B , LDA > 4, P < 0.05).

To study the difference in microbial community abundance between two groups of samples, Metastats software was used to conduct t -test on genus richness data between groups ( Figure 12 ). At the genus level, Actinobacillus , Pseudoxanthomonas were significantly enriched in aSHI medium while Bilophila was significantly enriched in aNEMM medium ( Figure 12A ). Unclassified Clade III, NS3a marine group and unclassified Saprospiraceae were more abundant in nNEMM suspension ( Figure 12B ). The same tendency can also be seen in Fig. S4 . From the taxa tree, aNEMM was suitable for culturing Anaerolinea , Aerococcus , Lactobacillus and Negativicutes . Planococcus , Desemzia and Psychrobacter were enriched in NEMM no matter under normoxia or anaerobic condition. Geobacillus , Lactococcus , Clostridiaceae , Hyphomicrobium , Pseudolabrys , Paracoccus , Rhodobacter , Rickettsiaceae , Sphingobium , Legionella , Moraxella , Luteimonas and Lysobacter were suitable for growth in nNEMM. Bacillus , Pediococcus , Brevundimonas , Caulobacter , Phenylobacterium , Sphingomonas and Pseudoxanthomonas were suitable for growing in aSHI. Blastomonas were enriched in nSHI. Abiotrophia grew well in SHI both under normoxia and anaerobic conditions ( Fig. S4 ).

Microbiota developed in NEMM and SHI media, under anaerobic and normoxia conditions exhibited distinct community structures. According to the relative abundance ratios of genera shown on taxa tree ( Fig. S4 ), we constructed Table 2 for clearly showing the different microbial profiles cultured in NEMM medium and SHI medium. At the genus level, Anaerolinea , Planococcaceae , Acerococcus , Desemzia , Sphongopyxia , Psychrobacter and some not assigned species were significantly enriched in aNEMM. Planococcaceae , Geobacillus , Facklamia , Desemzia , Lactococcus , Blastomonas , Sphongopyxia , Legionella , Psychrobacter , Luteimonas , Lysobacter and Pseudoxanthomonas were more abundant in nNEMM conditions. The aSHI condition was suitable for Pseudograclibacillus , Sphingomonas and Pseudoxanthomonas , while Lactococcus , Blastomonas , Sphingomonas , Legionella , Luteimonas , Lysobacter and Pseudoxanthomonas were more abundant in nSHI condition.

Notes:

Note that means significantly enriched.

NEMM supported the growth of many oral species that have not been cultured from S-ECC clinical samples so far. All 128 previously uncultivated species were detected by 16S rRNA gene sequencing using SILVA as the reference database.

To study the difference in microbial community abundance between groups of samples, Metastats analysis was used to conduct t -test on species richness data between groups. Uncultured bacteria including methanogenic_archaeon , Truepera_sp. , Acidobacterium_sp. , Desulfovibrionaceae_bacterium and Nitrosomonadaceae_bacterium could only be detected in aNEMM group. Comparing relative abundance of these uncultured species between aNEMM and aSHI media, uncultured bacteria including forest_soil_bacterium , Clostridiales_bacterium , Desulfuromonadales_bacterium , Firmicutes_bacterium and Alphaproteobacteria_bacterium were observed higher in aNEMM group ( Table S4 ). Discrepancies were also detected between the nNEMM and nSHI groups, with higher abundances of uncultured bacteria including Desulfuromonadales_bacterium , Bacteroidetes_bacterium , Clostridium_sp. , Mollicutes_bacterium , Bacteroidales_bacterium , Holophaga_sp. and Acidobacteria_bacterium in the nNEMM group. Moreover, uncultured bacteria including methanogenic_archaeon , Sphingobacteriales_bacterium , Aminicenantes_bacterium , marine_bacterium , forest_soil_bacterium , Acidobacterium_sp. , Xanthomonadaceae_bacterium , Sorangiineae_bacterium , marine_microorganism and candidate_division_SBR1093_bacterium were only observed in the nNEMM group ( P < 0.05) ( Table S5 ).

The serial diluted culturomics samples derived from above broth medium conditions were inoculated into the two media-agar plates and incubated under normoxia and anaerobic, respectively, for 1-5 days ( Fig. S5 ). The bacterial diversity gradually decreased over time. With the extension of culture time, the difference of dominant bacteria in the two media gradually decreased, and from the morphological point of view, SHI and NEMM agar plates tended to cultivate similar bacterial diversity. After this process, the slow-growing colonies were subcultured and identified to isolate previously uncultivated phylotypes. Through culturomics, uncultivated bacteria belonging to the Prevotella , Actinomyces , Pseudomonas , Granulicatella , Lactobacillus , Campylobacter , Streptococcus and Bacillus genera were identified in the biofilm growth media. The oral taxa identified under the four conditions were shown in Fig. S6 , displaying that both NEMM and SHI media could culture the uncultivated bacterium in aerobic and anaerobic conditions.

Also, it turned out that Achromobacter spanius , Lactobacillus delbrueckii subsp and Actinomyces naeslundii in aNEMM, Pseudomonas sp . and Veillonella parvula in aSHI, Brevibacillus agri in nNEMM, Leuconostoc lactis , Chryseobacterium indologenes , Cardiobacterium hominis and Kingella bonacorsii in nSHI may tend to be enriched in each growth condition.