Hub genes and pathways related to caries-free dental biofilm: clinical metatranscriptomic study

Objective This study aimed to evaluate the microbial functional profile of biofilms related to caries-free (CF, n = 6) and caries-arrested (CI, n = 3) compared to caries-active (CA, n = 5) individuals. Materials and methods A metatranscriptomic was performed in supragingival biofilm from different clinical conditions related to caries or health. Total RNA was extracted and cDNAs were obtained and sequenced (Illumina HiSeq3000). Trimmed data ( SortMeRNA ) were submitted to the SqueezeMeta pipeline in the co-assembly mode for functional analysis and further differential gene expression analysis (DESeq2) and weighted gene co-expression network analysis (WCGNA) to explore and identify gene modules related to these clinical conditions. Results A total of 5303 genes were found in the metatranscriptomic analysis. A co-expression network identified the most relevant modules strongly related to specific caries status. Correlation coefficients were calculated between the eigengene modules and the clinical conditions (CA, CI, and CF) discriminating multiple modules. CA and CI showed weak correlation coefficient strength across the modules, while the CF condition presented a very strong positive correlation coefficient ( r = 0.9, p value = 4 × 10 −9 ). Pearson’s test was applied to further analyze the module membership and gene significance in CF conditions, and the most relevant were HSPA1s -K03283, Epr - K13277, and SLC1A -K05613. Gene Ontology (GO) shows important bioprocesses, such as two-component system, fructose and mannose metabolism, pentose and glucuronate interconversions, and flagellar assembly ( p -adjust < 0.05). The ability to use different carbohydrates, integrate multiple signals, swarm, and bacteriocin production are significant metabolic advantages in the oral environment related to CF. Conclusions A distinct functional health profile could be found in CF, where co-occurring genes can act in different pathways at the same time. Genes HSPA1s, Epr , and SLC1A may be appointed as potential biomarkers for caries-free biofilms. Clinical relevance Potential biomarkers for caries-free biofilms could contribute to the knowledge of caries prevention and control.