Oral microbial taxa associated with risk for SARS-CoV-2 infection

Hypothesis and objective: The oral and digestive tract microbial ecosystem has sparked interest because of its impact on various systemic diseases and conditions. The oral cavity serves not only as a reservoir for many potentially virulent microbiota but also as an important entry point and portal to the human body system. This is especially significant in the transmissibility of the virulent current pandemic virus SARS-CoV-2. The oral and digestive microbiome influences the inflammatory burden and effectiveness of the immune system and serves as a marker of activity of these host processes. The host immune response plays a role in infection susceptibility, including SARS-CoV-2. The purpose of this study is to investigate the role of specific salivary oral microbiome in susceptibility to SARS-CoV-2 infection.

Methods and results: One hundred six subjects of known medical and dental history who consented to provide saliva samples between January 2017 and December 2019 were included in this study. Sixteen had become COVID-19 positive based on the PCR test by 3/01/2021. A comparison of oral microbiome bacteria taxa profiles based on 16S rRNA sequencing revealed differences between the two groups in this pilot study.

Conclusions: These bacteria taxa may be markers of increased susceptibility to SARS-CoV-2 infection in the unvaccinated population.

Introduction

The microbiome of the oral cavity and the digestive tract exists in a delicate homeostasis, and this balance of bacteria plays an important role in the body's immune and inflammatory potency levels [1, 2]. The oral and nasal mucosal membranes are the first entry point for many respiratory viruses, and the microbiome plays a significant role in the pathogenesis of these viral infections [3]. There is certainly precedent for the idea that microbiota, whether digestive (gut/oral) or respiratory, could play a role in the resistance and response to infection by viruses such as the SARS-CoV-2 virus [4], associated with 5.3 million deaths and 272 million cases worldwide by the end of 2021. Salivary microbiota, specifically bacteria, play a significant role in facilitating attachment of incoming microbes into its ecosystem and, from there, breakthrough to the mucosal membranes of the host system. Salivary microbiota have the potential to control viral uptake into cells and whether virus introduction into the host is sufficient for a productive, detectable infection. In addition, these microbiota can create molecules that alter host epithelial barriers [5].

Many studies have identified changes in the gut, oral, and nasopharyngeal bacteria that occur during and post SARS-CoV-2 infection vs. healthy controls [6–11]; however, the mechanism of these changes remains unknown. The SARS-CoV-2 virus has a profound effect on the oral cavity's function, with symptoms including loss of taste, dry mouth, dysgeusia, and opportunist infections that may alter oral microbiota [12]. We can presume that changes in other oral microbiota due to SARS-CoV-2 infection may contribute to these symptoms within the oral cavity. Additional work has attempted to study the changes in oral/gut bacteria that occur during the course of SARS-CoV-2 active infection and its correlation with severity of disease [8, 11, 13–15]. These studies shed light on bacteria that have the potential to increase the morbidity rate with SARS-CoV-2 infection at sites such as the oro-pharynx [16, 17].

The COVID-19 etiological agent, SARS-CoV-2 betacoronavirus, has high spreading capability, mediated by several vehicles, such as oral droplets [18–20]. The oral cavity might be considered a preferred route for SARS-CoV-2 entry, and there is ample evidence that spread primarily via a droplet [2]. These droplets can enter the human body through the mucosal membranes of, primarily, the nasal cavity via the ACE2 + TMPRSS2 + nasal epithelial cells [3]. While there has been more investigation of the nasal-lung axis of infection [4], infection can occur via the oral cavity [1, 2]. This is supported by high expression of both the ACE2 and TMPRSS2 receptors in the salivary glands and the oral cavity mucosa. While oral manifestation of SARS-CoV-2 infection occurs in oral tissues, it is not clear how often it is the primary site of SARS-CoV-2 initial infection [21]. After entry into the oral cavity, a successful infection requires survival of the virus in the saliva or on the mucosa until its adherence to the mucosal surface, followed by penetrance through the ECM binding to the ACE receptor for uptake, and then successful viral reproduction to the point that the level of viral particles exceed the threshold for detection. The question addressed here is whether certain salivary oral microbiota may facilitate attachment and penetration, increasing the risk of a successful infection process.

This was a retrospective study; the purpose was to investigate preexisting oral salivary bacterial diversity and its impact on susceptibility to SARS2-CoV-2 infection. Salivary samples were collected from patients of record, with known pre-pandemic medical and dental history, who were followed through the pandemic. The saliva microbiome is known to be relatively stable over time for individual adults [22, 23], and this method offered the unusual proposition that the salivary microbiome was known prior to infection. This allowed the inquiry into whether there are salivary bacteria taxa profiles associated with susceptibility to successful COVD-19 infections. Based on the location and time of the study in central Illinois, this analysis is of the association between successful infection by alpha, beta, and gamma variants of SARS-CoV-2 and levels of specific oral bacteria taxa (https://www.cdc.gov/coronavirus/2019-ncov/variants/variant-classifications.html).

Materials and methods

Study population and patient characteristics

This study was on a subset of subjects from a cross-sectional study conducted at the University of Illinois at Chicago College of Dentistry, General Practice and Prosthodontic clinics between 01/08/2017 and 06/21/2019 [24]. Approximately, 272 subjects provided written informed consent to participate in accordance with guidelines of the institutional ethics committee of the University of Illinois at Chicago, Institutional Review Board 1, which approved this study #2016-0696. This study was done in full accordance with the principles of the Declaration of Helsinki. Inclusion criteria: 18 years of age and older, available medical and dental record including a current medication list, a full periodontal exam, a visual, tactile, and radiographic caries exam, dentate or edentulous, and consent to supply a saliva sample. Exclusion criteria: presence of restored dental implants, removable partial dentures, maxillofacial defects scaling of teeth within the past 3 months; acute disease that requires urgent care, less than twenty natural teeth for the dentate subjects, antibiotic use within the past month or mouthwash in the last 12 h, vaccination of SARS-CoV-2 prior to the March 1, 2021 cutoff date, unknown COVID-19 history.

Sample collection

Stimulated saliva was collected from patients asked to chew paraffin over a 5-min period as part of an earlier study [25]. The subjects all had had full dental exams prior to the study.

Evaluation of COVID-19 status

All the subjects were contacted by telephone at least three times between April and May of 2021 before exclusion due to no response. Those who responded were subjected to a scripted telephone interview to determine presence of respiratory illness between February 1, 2020 and March 1, 2021, whether they had tested positive for COVID-19 by the PCR test, and whether they had been vaccinated for SARS-CoV-2 before March 1, 2021.

Characterization of microbial community structure

Genomic DNA was extracted from saliva samples, followed by amplification of the V1–V3 variable region of bacterial 16S ribosomal RNA (rRNA) genes using the primer pair 27F/534R and a two-stage targeted amplicon sequencing as described previously [25, 26]. The second PCR amplification used FluidigmAccess Array barcoded primers to allow unique labeling of each sample at the University of Illinois at Chicago Sequencing Core, which was followed by sequencing on an Illumina MiSeq with the V3 kit to allow 600 cycles total.

For taxa assignment and measurement, reverse sequence reads from the FASTQ files were analyzed using the software package QIIME2 (v2022.2) [27]. Sequences were trimmed if the average quality was lower than 25. As a result, the read sequences were truncated at 262nt. DADA2-plugin in QIIME2 was used to denoise the sequence and generate feature data and feature tables for the dataset of DNA sequences [28]. Taxonomy assignment was done by classify-consensus-blast function using the Blast+ consensus taxonomy classifier to determine 98% match identity of the query sequences to the Human Oral Microbiome Database (v15.22) [29]. Of the 106 samples, there were 25,136 reads per sample on average (range, 12,645 – 36,906). There were 2,664,478 million reads total.

Statistical analysis

Alpha diversity and Beta diversity analyses were performed using MicrobiomeAnalyst [30]. This analysis included calculation of the Shannon's Diversity Index of both species number and their distribution, and Chao1 indices of richness. These were compared between groups using a t-test within MicrobiomeAnalyst. Beta diversity analysis was visualized using the Bray–Curtis dissimilarity (non-phylogenetic) metric. MaAsLin2, using a Zero-inflated negative binomial regression with CSS normalization, allowed determination of taxa that were differentially represented in saliva from the subjects who became COVID positive vs. those who did not [31, 32]. This was done after internal batch correction by setting random_effects in the program to one of the two sequence runs. At least 2 reads were required per taxa measurement, and taxa that appeared in less that 10% of the total number of the samples were omitted. For statistical analysis of the study population demographics, two-tailed Fisher's Exact and Student's T-tests were run within the Kalediagraph program (Synergy Software). For analyses using MaAsLin2, all q values are determined using the Benjamini Hochberg procedure.

Results

Participants and population description

Approximately, 272 subjects who had provided saliva samples between November 2017 and December 2019 were contacted by phone during April 2021. Of these, 166 either did not know their COVID-19 status or were not reachable by phone and were excluded from the study. One hundred six subjects knew their prior COVID-19 infection status at the cut-off date of March 1, 2021; 16 had been PCR-verified COVID-19 positive at some point and were classified as positive. There were 90 subjects in the COVID-19 negative group. These people reported no respiratory illness between 1 February 2020 and 1 March 2021 and no positive SARS-CoV-2 PCR test at any point.

The salivary microbiome was compared between the positive and negative groups in this pilot study. Table 1 shows the general characteristics of the study group at the time of sample acquisition in terms of dental health, tobacco usage, number of prescribed medications, and age. The groups are similar overall, with a trend toward a higher level of tobacco usage in the COVID-19 positive group.

TABLE 1

Table 1. Demographics of the study population.

Microbial diversity analysis and differentially abundant taxa

The Alpha diversity of the two groups shows no difference between the groups (Figure 1). The Beta diversity overall does not show a significant difference between the two groups (Figure 2). This is typical for saliva samples, which are a mixture of taxa from multiple niches.

FIGURE 1

Figure 1. Phylogenetic diversity of saliva from subjects that self-report as being COVID-19 positive and those that did not. A box plot of Comparison of Shannon Diversity indices shows insignificant differences in taxa richness in the two groups, p < 0.543 and [t-test] statistic: −0.616. Comparison of taxa richness based on the CHAO1 indices reveals the same, p < 0.2510, (t-test) statistic: −1.19.

FIGURE 2

Figure 2. Principal coordinate analysis (PCoA) of saliva microbiome profiles of both groups based on Bray-Curtis Distances.

Taxonomic analysis reveals limited differences between the two groups. Usage of analysis tool MaAsLin2 in Table 2 shows bacterial taxa: Schaaliacardiffensis, g_Bergeyella, g_Bacteroidetes_[G-5] s__bacterium_HMT_511, g__Neisseria. s__elongata, g__Prevotella.s__dentalis, which differ in levels between the two groups at FDR <0.1 (Table 2) [31, 32]. Multivariable analysis with correction for levels of tobacco usage, age, and medication count is shown in Figure 3; Supplementary Table 1.

TABLE 2

Table 2. MaAsLin2-based determination of taxa that are differentially abundant in saliva of the subjects that became COVID-19 positive.

FIGURE 3

Figure 3. Significant associations (FDR < 0.25) between bacterial taxa and clinical variables detected by MaAsLin2-based multiple variable analysis. Each association such as future COVID-19 status and salivary taxa were adjusted for the remaining factors, such as subject's tobacco usage, age, and the number of medications prescribed to the subject.

Discussion

A small group of oral taxa was shown to be differentially abundant in subjects who went on to be COVID-19 positive. Two of these, Prevotella and Neisseria, are of interest in that they have been tied to COVID-19 in the past. Multivariable analysis (Figure 3) revealed these taxa may be associated in part with differences in age, tobacco usage, and/or numbers of medications required in the groups studied.

A species of the genus Prevotella was observed to be more abundant in the saliva of patients who went on to become SARS-CoV-2 positive (Table 2).

Earlier studies have noted Prevotella pallens was enriched in saliva samples from patients with COVID-19 with active disease, while Rothia mucilaginosa and some Streptoccoci species were at lower levels [13]. A second study of patients with active disease vs. healthy controls noted oral mucosa were enriched for Haemophilus parainfluenzae, Veillonella infantium, Soonwooa purpurea, Prevotella salivae, Prevotella jejuni, and Capnocytophaga gingivalis, while displaying lower levels of a number of taxa [7]. A third study did not see this higher level of Prevotella at least on the genus level in saliva samples of patients with COVID-19 [33]. Based on the enrichment of Prevotella genes in lung samples from COVID-19 patients with active disease and computational identification of possible roles for Prevotella-encoded proteins in host-viral interactions, there is speculation that this enrichment of Prevotella in some patients could play a fairly direct role in stimulating SARS-CoV2 infection and subsequent disease [34]. Comparison between studies is made difficult by the variety of 16S rRNA amplicons that are examined in the different studies. Nevertheless, the higher level of the single Prevotella species seen in the subjects who became COVID-19 positive (Figure 3; Table 2) may provide some support for that model.

Low levels of the bacteria genus Neisseria in patients with COVID-19 have been proposed to play a role in increased inflammation in some patients [35]. There are several reports of lower levels of Nesseria overall in patients with COVID-19, and Iebba et al. report lower levels, specifically of Neiserria elongata in oral rinses from patients with COVID-19 at FDR <0.11 [7, 8, 15]. As shown in Table 2, Neiserria elongata was at a lower level in the patients who went on to contract COVID-19. One might speculate that this species is lower in COVID-19-positive patients before they get the disease.

While the taxa that differ between the two groups may serve as markers of COVID-19 disease susceptibility, the mechanism is obscure. Differences in resistance to viral establishment could be tied to genetics of the host, inflammatory state, immune function, and the microbial organisms that can affect each other [36]. At the time of the initiation of this work, it was not clear if the oral cavity is a primary site of SARS-CoV-2 infection, and this is still unclear [37, 38]. For that reason, the few differences seen in salivary microbiome in the group that had been COVID-19 positive could be manifestations of lifestyle associated with increased risk anywhere in the body. This is not an unprecedented concept, in that there is ample evidence of oral microbiome differences being linked to differential disease risk throughout the body [39].

This study has several limitations: it was an observational, retrospective, single-center study. Multiple factors can cause minor differences in oral microbiota, so any relationships may be indirect, and the number of COVID-positive subjects limits sensitivity of the study [39, 40]. In addition, it is possible that a small number of COVID-negative subjects had been positive at some time, but this was not known to them [41]. However, the large number of the subjects in the COVID-negative group would minimize that as a source of error. We used V2-V3 16S rRNA analysis, which provides reproducible information on the levels of a large percentage of known oral bacteria taxa, but not all oral bacteria were measured [42]. This study has additional limitations: (1) A lack of functional analysis of microbial genes, (2) Collection of saliva in some cases had occurred over a year prior to the COVID-19 pandemic, (3) The date of COVID-19 positivity was not recorded. (4) The presence of restorations such as dental crowns was not considered, although the patients with implants were excluded. A prospective study with more frequent saliva collection would have addressed some of these limitations.

Given the March 1, 2021 cutoff for previous infection status and the location in the United States, findings may be limited to the alpha and beta SARS2-CoV-2 strains [43]. In addition, salivary properties of viscosity, pH level, and immunoglobulin and inflammatory cytokines levels were not assessed, which may play a role in the protective mechanism against microbial attachment and penetration [44].

Conclusion

These bacteria taxa present in the saliva may be associated with increased risk of early SARS-CoV-2 infection in the unvaccinated population. This study can serve as a prelude to the study of the diverse oral microbiota and its effect on the protective mechanism or facilitation of infection by respiratory viruses.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: BIoproject accession numbers: PRJNA674379 and PRJNA739492.

Ethics statement

The studies involving human participants were reviewed and approved by UIC IRB. The patients/participants provided their written informed consent to participate in this study.

Author contributions

NC, MH, TB, GA, and NK contributed to the performed data collection and analysis. NC and GA wrote the first draft of the manuscript. All the authors contributed to the editing and agree with the content of the manuscript.

Acknowledgments

The authors acknowledge the Research Open Access Publishing (ROAAP) Fund of the University of Illinois Chicago for financial support towards the open access publishing fee for this article.

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/froh.2022.886341/full#supplementary-material

References

1. Herrera D, Serrano J, Roldán S. Sanz, M. Is the oral cavity relevant in SARS-CoV-2 pandemic? Clin Oral Investig. (2020) 24:2925–30. doi: 10.1007/s00784-020-03413-2

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Huang N, Pérez P, Kato T, Mikami Y, Okuda K, Gilmore RC, et al. SARS-CoV-2 infection of the oral cavity and saliva. Nat Med. (2021) 27:892–903. doi: 10.1038/s41591-021-01296-8

PubMed Abstract | CrossRef Full Text | Google Scholar

3. World Health Organization. Modes of Transmission of Virus Causing COVID-19: Implications for IPC Precaution Recommendations (2020). Available online at: https://www.who.int/docs/defaultsource/coronaviruse/20200329-scientific-brief1-final-corrected2-revised-final.pdf

Google Scholar

4. Tada A, Senpuku H. The impact of oral health on respiratory viral infection. Dent J. (2021) 9:4–43. doi: 10.3390/dj9040043

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Lin D, Yang L, Wen L, Lu H, Chen Q, Wang Z. Crosstalk between the oral microbiota, mucosal immunity, and the epithelial barrier regulates oral mucosal disease pathogenesis. Mucosal Immunol. (2021) 14:1247–58. doi: 10.1038/s41385-021-00413-7

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Gupta A, Karyakarte R, Joshi S, Das R, Jani K, Shouche Y, et al. Nasopharyngeal microbiome reveals the prevalence of opportunistic pathogens in SARS-CoV-2 infected individuals and their association with host types. Microbes Infect. (2022) 24:104880. doi: 10.1016/j.micinf.2021.104880

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Iebba V, Zanotta N, Campisciano G, Zerbato V, Di Bella S, Cason C, et al. Profiling of oral microbiota and cytokines in COVID-19 patients. Front Microbiol. (2021) 12:671813. doi: 10.3389/fmicb.2021.671813

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Ma S, Zhang F, Zhou F, Li H, Ge W, Gan R, et al. Metagenomic analysis reveals oropharyngeal microbiota alterations in patients with COVID-19. Signal Transduct Target Ther. (2021) 6:191. doi: 10.1038/s41392-021-00614-3

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Sarkar A, Harty S, Moeller AH, Klein SL, Erdman SE, Friston KJ, et al. The gut microbiome as a biomarker of differential susceptibility to SARS-CoV-2. Trends Mol Med. (2021) 27:1115–34. doi: 10.1016/j.molmed.2021.09.009

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Wang H, Wang H, Sun Y, Ren Z, Zhu W, Li A, et al. Potential associations between microbiome and COVID-19. Front Med (Lausanne). (2021) 8:785496. doi: 10.3389/fmed.2021.785496

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Yamamoto S, Saito M, Tamura A, Prawisuda D, Mizutani T, Yotsuyanagi H. The human microbiome and COVID-19: A systematic review. PLoS ONE. (2021) 16:e0253293. doi: 10.1371/journal.pone.0253293

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Cuevas-Gonzalez MV, Espinosa-Cristóbal LF, Donohue-Cornejo A, Tovar-Carrillo KL, Saucedo-Acuña RA, García-Calderón AG, et al. COVID-19 and its manifestations in the oral cavity: a systematic review. Medicine (Baltimore). (2021) 100:e28327. doi: 10.1097/MD.0000000000028327

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Miller EH, Annavajhala MK, Chong AM, Park H, Nobel YR, Soroush A, et al. Oral microbiome alterations and SARS-CoV-2 saliva viral load in patients with COVID-19. Microbiol Spectr. (2021) 9:e0005521. doi: 10.1128/Spectrum.00055-21

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Soffritti I, D'Accolti M, Fabbri C, Passaro A, Manfredini R, Zuliani G, et al. Oral Microbiome dysbiosis is associated with symptoms severity and local immune/inflammatory response in COVID-19 patients: a cross-sectional study. Front Microbiol. (2021) 12:687513. doi: 10.3389/fmicb.2021.687513

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Wu Y, Cheng X, Jiang G, Tang H, Ming S, Tang L, et al. Altered oral and gut microbiota and its association with SARS-CoV-2 viral load in COVID-19 patients during hospitalization. NPJ Biofilms Microbiomes. (2021) 7:61. doi: 10.1038/s41522-021-00232-5

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Lansbury L, Lim B, Baskaran V, Lim WS. Co-infections in people with COVID-19: a systematic review and meta-analysis. J Infect. (2020) 81:266–75. doi: 10.1016/j.jinf.2020.05.046

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Zhang H, Zhang Y, Wu J, Li Y, Zhou X, Li X, et al. Risks and features of secondary infections in severe and critical ill COVID-19 patients. Emerg Microbes Infect. (2020) 9:1958–64. doi: 10.1080/22221751.2020.1812437

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Bao L, Zhang C, Dong J, Zhao L, Li Y, Sun J. Oral microbiome and SARS-CoV-2: beware of lung co-infection. Front Microbiol. (2020) 11:1840. doi: 10.3389/fmicb.2020.01840

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Netz RR, Eaton WA. Physics of virus transmission by speaking droplets. Proc Natl Acad Sci U S A. (2020) 117:25209–11. doi: 10.1073/pnas.2011889117

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Stadnytskyi V, Anfinrud P, Bax A. Breathing, speaking, coughing or sneezing: what drives transmission of SARS-CoV-2? J Intern Med. (2021) 290:1010–27. doi: 10.1111/joim.13326

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Brandão TB, Gueiros LA, Melo TS, Prado-Ribeiro AC, Nesrallah A, Prado GVB, et al. Oral lesions in patients with SARS-CoV-2 infection: could the oral cavity be a target organ? Oral Surg Oral Med Oral Pathol Oral Radiol. (2021) 131:e45–51. doi: 10.1016/j.oooo.2020.07.014

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Cameron SJ, Huws SA, Hegarty MJ, Smith DP, Mur LA. The human salivary microbiome exhibits temporal stability in bacterial diversity. FEMS Microbiol Ecol. (2015) 91:fiv091. doi: 10.1093/femsec/fiv091

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Yamanaka W, Takeshita T, Shibata Y, Matsuo K, Eshima N, Yokoyama T, et al. Compositional stability of a salivary bacterial population against supragingival microbiota shift following periodontal therapy. PLoS ONE. (2012) 7:e42806. doi: 10.1371/journal.pone.0042806

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Schwartz JL, Peña N, Kawar N, Zhang A, Callahan N, Robles SJ, et al. Old age and other factors associated with salivary microbiome variation. BMC Oral Health. (2021) 21:490. doi: 10.1186/s12903-021-01828-1

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Adami GR, Ang MJ, Kim EM. Comparison of microbiome in stimulated saliva in edentulous and dentate subjects. Methods Mol Biol. (2021) 2327:69–86. doi: 10.1007/978-1-0716-1518-8_5

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Naqib A, Poggi S, Wang W, Hyde M, Kunstman K, Green SJ. Making and Sequencing Heavily Multiplexed, High-Throughput 16S Ribosomal RNA Gene Amplicon Libraries Using a Flexible, Two-Stage PCR Protocol. Methods Mol Biol. (2018) 1783:149–69. doi: 10.1007/978-1-4939-7834-2_7

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. (2019) 37:852–7. doi: 10.1038/s41587-019-0209-9

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJ, Holmes SP. DADA2: High-resolution sample inference from Illumina amplicon data. Nat Methods. (2016) 13:581–3. doi: 10.1038/nmeth.3869

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Dewhirst FE, Chen T, Izard J, Paster BJ, Tanner AC, Yu WH, et al. The human oral microbiome. J Bacteriol. (2010) 192:5002–17. doi: 10.1128/JB.00542-10

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Chong J, Liu P, Zhou G, Xia J. Using microbiomeanalyst for comprehensive statistical, functional, and meta-analysis of microbiome data. Nat Protoc. (2020) 15:799–821. doi: 10.1038/s41596-019-0264-1

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Mallick H, Rahnavard A, McIver LJ, Ma S, Zhang Y, Nguyen LH, et al. Multivariable association discovery in population-scale meta-omics studies. PLoS Comput Biol. (2021) 17:e1009442. doi: 10.1371/journal.pcbi.1009442

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Zhang X, Mallick H, Yi N. Zero-inflated negative binomial regression for differential abundance testing in microbiome studies. J Bioinf Genom. (2016) 2, 2–4. doi: 10.18454/jbg.2016.2.2.1

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Rafiqul Islam SM, Foysal MJ, Hoque MN, Mehedi HMH, Rob MA, Salauddin A, et al. Dysbiosis of Oral and Gut Microbiomes in SARS-CoV-2 Infected Patients in Bangladesh: Elucidating the Role of Opportunistic Gut Microbes. Front Med (Lausanne). (2022) 9:821777. doi: 10.3389/fmed.2022.821777

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Khan AA, Khan Z. COVID-2019-associated overexpressed Prevotella proteins mediated host-pathogen interactions and their role in coronavirus outbreak. Bioinformatics. (2020) 36:4065–9. doi: 10.1093/bioinformatics/btaa285

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Demirci M. Could Neisseria in oral microbiota modulate the inflammatory response of COVID-19? Oral Dis. (2021). doi: 10.1111/odi.14082. [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Uchida H, Ovitt CE. Novel impacts of saliva with regard to oral health. J Prosthet Dent. (2022) 127:383–91. doi: 10.1016/j.prosdent.2021.05.009

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Casillas Santana MA, Dipp Velázquez FA, Sámano Valencia C, Martínez Zumarán A, Zavala Alonso NV, Martínez Rider R, et al. Saliva: what dental practitioners should know about the role of this biofluid in the transmission and diagnostic of SARS-CoV-2. Medicina (Kaunas). (2021) 57. doi: 10.3390/medicina57040349

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Sakaguchi W, Kubota N, Shimizu T, Saruta J, Fuchida S, Kawata A, et al. Existence of SARS-CoV-2 Entry Molecules in the Oral Cavity. Int J Mol Sci. (2020) 21. doi: 10.3390/ijms2117600

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Belstrøm D. The salivary microbiota in health and disease. J Oral Microbiol. (2020) 12:1723975. doi: 10.1080/20002297.2020.1723975

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Wade WG. The oral microbiome in health and disease. Pharmacol Res. (2013) 69:137–43.

Google Scholar

41. Gao Z, Xu Y, Sun C, Wang X, Guo Y, Qiu S, et al. A systematic review of asymptomatic infections with COVID-19. J Microbiol Immunol Infect. (2021) 54:12–6. doi: 10.1016/j.jmii.2020.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Lazarevic V, Whiteson K, Gaïa N, Gizard Y, Hernandez D, Farinelli L, et al. Analysis of the salivary microbiome using cultureindependent techniques. J Clin Bioinforma. (2012) 2:4. doi: 10.1186/2043-9113-2-4

PubMed Abstract | CrossRef Full Text | Google Scholar

43. SARS-CoV-2 Variant Classifications and Definitions. Centers for Disease Control and Prevention (2021).

Google Scholar

44. Wade WG. Resilience of the oral microbiome. Periodontol 2000. (2021) 86:113–22. doi: 10.1111/prd.12365

PubMed Abstract | CrossRef Full Text | Google Scholar