Study on Tongue Coating Microbiota in Patients with Atrophic Gastritis ()
1. Introduction
The human is a complex and diverse ecosystem consisting of trillions of microbes (i.e. bacteria, archaea, fungi, protists and viruses) [1], the majority of which are essential in human development, and the complex interactions between the human and the microbes confer benefits to the host in many key aspects of life [2]. The microbes are thought to be the second genome of the human body [3], as the human gut microbes alone have up to 3 million genes, 150 times the number of human genes, whose function is not just helpful in physiologic activities like digestion, absorption, protection, growth but also plays an important role in human health and disease through participating in various biochemical pathways [4] [5]. So far, the microbes have been implicated to be closely related to many chronic diseases, such as obesity [6], malignancies [7], diabetes [8] and so on.
Atrophic gastritis (AG) is a relatively common disease globally, the incidence rate of which ranged from 0% to 10.9% [9], characterized by the replacement of the original gastric glands by intestinal metaplasia (IM), pseudopyloric metaplasia, and/or fibrosis [10]. As the human model of gastric carcinogenesis with the following sequential stages: chronic gastritis → atrophy → intestinal metaplasia → dysplasia is increasingly accepted [11], the early diagnosis and treatment of atrophic gastritis are of great significance. It has been confirmed that the bacterium is the cause of gastritis in an overwhelming majority of the cases because of the discovery of the Helicobacter pylori (Hp) by Warren and Marshall in 1982 [12] and from current knowledge, the pathogenic mechanisms between bacterium and gastritis may involve in exacerbating gastric mucosal inflammation, altering the host immune system, and promoting the transformation of dietary nitrate carcinogens [13] [14] [15]. However, most of the current studies focus on the gastrointestinal flora, and the research on the relationship between tongue coating bacteria and atrophic gastritis is still in its infancy.
As it is determined that factors including age, gender can significantly alter the composition of the microbiota [1] [16], this study was conducted in patients with atrophic gastritis and non-atrophic gastritis matched with their gender and age. In the present study, we characterized the tongue coating microbiota to evaluate the differences of tongue coating microbiota in patients with atrophic gastritis and non-atrophic gastritis by using the 16S rRNA gene sequencing technology.
2. Materials and Methods
2.1. Subject
The 29 atrophic gastritis patients confirmed by endoscopy and pathology at the Department of Gastroenterology of Wenzhou People’s Hospital in Zhejiang, from May 2019 to November 2019 were recruited as cases. In the same period, we included 29 non-atrophic gastritis patients matched with age (±3) and gender of atrophic gastritis patients as controls. All participants were briefed on the purpose, procedures, potential risks and benefits of the study, and informed consent was signed and those who had taken antibiotics or microecological modulators within 2 weeks before were also ruled out. The study was approved by the Ethical Committee of Wenzhou People’s Hospital.
2.2. Sample Collections
Participants were not allowed to eat or drink before this procedure. A sterile swab was used to gently scrape the tongue coating on the patient for 5 - 6 times with standard operations. After collection, the specimen was placed in a sterile tube, transported to the laboratory with liquid nitrogen and stored at −80˚C freezing storage until further processing.
2.3. 16S rRNA Gene Sequence Analysis
The 16S rRNA sequencing [17] analysis was used to determine the tongue coating microbiota. The microbial DNA was extracted using the MiSeq Reagent Kit v3 (Illumina, USA). Agarose gel electrophoresis and Nanodrop 2000 (Thermo Fisher Scientific, USA) were used to determine the integrity of the genomic DNA and PCR amplification was carried out in an ABI 2720 Thermal Cycler (Thermo Fisher Scientific, USA). We used isolated DNA as the template for amplification of the V3-V4 hypervariable region of the 16S rRNA gene in three replicate reactions with the primer pair F (Illumina adapter sequence 1 + CCTACGGGNGGCWGCAG) and R (Illumina adapter sequence 2 + GACTACHVGGGTATCTAATCC). The replicate PCR products were pooled and purified with Agencourt AMPure XP magnetic beads (Beckman Coulter, USA).
2.4. cDNA Library Construction and Sequencing
According to the preliminary quantitative results of agarose gel electrophoresis, the Qubit was used to accurately quantify the library. Insert length and integrity of the libraries were assessed by Agilent 2100 bioanalyzer (Agilent Technologies, USA) to check the quality of the library. Later, the library was sequenced by two-terminal sequencing strategy of 2 × 250 bp on Miseq platform according to standard protocols [18].
2.5. Statistical Analysis
All statistics were analyzed by the SPSS v25.0 (SPSS Inc., Chicago, IL) and R v3.0 (The R Foundation for Statistical Computing, Vienna, Austria) and all measurement data were presented as mean ± standard deviation. Paired t-test was used to analyze the participant characteristics between the two groups. The relative abundances in each sample were calculated on the basis of the sum of 16S rRNA sequences. Furthermore, the Alpha-and beta-diversity were estimated between the two groups using the R software. A value of p < 0.05 was considered to be significant in this work.
3. Results
In this study, 29 patients with atrophic gastritis including 16 males and 13 females and 29 non-atrophic gastritis matched according to gender and age (±3), were included in the study and the average age was 46.91 ± 8.41 years old. Differences in variables, such as gender, age, smoking, alcohol drinking, BMI and the presence Hp did not achieve significance (P > 0.05) between the two groups (Table 1).
3.1. Diversity of the Microbial Community between Two Groups
A total of 1000 operational taxonomic units (OTUs) were quantified in this study. The Venn diagram of the diversity showed that 695 OTUs were shared among the two groups, 106 OTUs were detected only in atrophic gastritis group and 199 OTUs were only quantified in the controls (Figure 1).
Chao1, abundance-based coverage estimator (ACE), Shannon and Simpson are the comprehensive indexes to reflect the alpha diversity of the microbial community in the tongue coatings. Chao1 and ACE reflect community richness, the trend of which was similar to that of OTUs. Shannon and Simpson indexes reflect community diversity. As listed (Table 2), there was no significant difference (P > 0.05) in microbial community diversity between the two groups.
Table 1. Characteristics of the two groups in the study.
Table 2. The alpha diversity of two groups.
Figure 1. A Venn diagram showing the comparison of microbiota diversity between the two groups by at the OTUs level. A: Atrophic gastritis group; B: Control group.
3.2. Relative Abundances of Microbes
The relative abundances of bacteria at the class level in each of the two groups was produced (Figure 2). Dominant bacteria of tongue coating (5% ≥ relative abundance) in the two groups were similar and the greatest relative abundances were Bacteroidia, Bacilli, Betaproteobacteria, Fusobacteriia, Negativicutes, Actinobacteria, Gammaproteobacteria, Clostridia. Moreover, patients with atrophic gastritis had lower relative abundances of Betaproteobacteria and Spirochaetia than the control group (P < 0.05) (Figure 3).
At the gene level, the dominant bacteria in the two groups were basically the same and relative abundance was slightly different. The abundance of atrophic
Figure 2. The relative abundances of bacteria at the class level in two groups.
Figure 3. Relative abundance of Betaproteobacteria and Spirochaetia in two groups. (*= P < 0.05).
gastritis group from high to low were: Prevotella (23.95%), Streptococcus (8.16%), Alloprevotella (6.97%), Fusobacterium (6.71%), Neisseria (6.53%). The abundance of non-atrophic gastritis group was followed by the Prevotella (20.39%), Neisseria (9.7%), Alloprevotella (8.69%), Streptococcus (7.28%), Fusobacterium (6.9%). Further analysis revealed that the relative abundance of Neisseria and Aggregatibacter in atrophic gastritis group were significantly lower than that in control group (P < 0.05) (Figure 4).
3.3. Functionally Predicted Metabolic Pathways Underlying the Two Groups
A total of 24 KEGG metabolic pathways were identified to show significantly different microbiota abundances between the two groups (Figure 5). The metabolic pathways significantly upregulated in atrophic gastritis group at level 3 were enriched in: taurine and hyposulfonic acid metabolism, ribonucleic acid polymerase, fructose and mannose metabolism, starch and sucrose metabolism, glycolysis/gluconeogenesis, etc. (P < 0.05). Other pathways including glucose biosynthesis, metabolism, α-linoleic acid metabolism, glyoxylate and dicarboxylate metabolism (P < 0.05), had more abundance in the control group than in the atrophic gastritis group.
4. Discussion
The oral cavity is connected with the respiratory system, the digestive system and the external environment, which provides suitable environment and conditions for the growth, reproduction and settlement of various microorganisms in the oral cavity, thus forming a complex oral microecological environment. The coating on the tongue is a thin, white and moist coating on the back of the tongue, mainly composed of desquamated epithelial cells, blood cells, metabolites, nutrients and bacteria, where is the main distribution center of oral flora [19]. In recent years, with the development of genomics, new achievements have been emerging in the research on the relationship between tongue coating microbiota and disease.
It has been found that the tongue microbiota is of great significance for oral
Figure 4. Relative abundance of Neisseria and Aggregatibacter in two groups.
Figure 5. Metabolic function analysis of tongue coating flora in two groups.
health, and the imbalance of oral microecology may lead to caries [20], halitosis [21], periodontitis [22] and other oral diseases. In addition, the imbalance of tongue microbiota may directly involve in the occurrence of diseases beyond the oral. Study has demonstrated a significant association between the Porphyromonas gingivalis and rheumatoid arthritis disease, which suggests that improving oral flora may be a new direction to alleviate the symptoms of rheumatoid arthritis and delay the disease progression [23]. Moreover, a large number of candida albicans were found in tongue coating cultures in patients with chronic kidney disease, implying that tongue coating may be one of the portals for disseminated infection [24]. The relative abundance of Neisseria, Haemophilus, Fusobacterium and Porphyromonas in the healthy people were higher than that in the cancers, suggested that tongue diagnosis may provide potential screening and early diagnosis method for cancer [25].
In this study, we used the 16S rRNA gene sequencing technology to analyze the composition and diversity of the tongue coating microbiota and compared them between atrophic gastritis group and non-atrophic gastritis group. Comparing bacterial taxa at each level with the control group, Betaproteobacteria and Spirochaetia were enriched in the non-atrophic gastritis group at the class level. At the genus level, the atrophic gastritis group had lower proportion of Neisseria and Aggregatibacter compared with the control group. Its mechanism may be involved in the local microenvironment inflammation, which in turn indirectly promotes the occurrence and development of gastritis. In addition, we used PICRUSt software to predict that these changes of tongue coating flora can cause functional changes in 24 pathways, suggesting that these functional changes may be the cause of the symptoms of atrophic gastritis.
Li et al. [26] found that the abundance of Campylobacter gradually increased during the progression of gastric mucosa from normal to superficial, atrophic gastritis, and intestinal metaplasia, which suggested that tongue coating flora may be expected to become a non-invasive new biomarker suitable for long-term monitoring of gastritis progression. Another study [27] found that the relative abundance of Proteobacteria, Neisseria and Haemophilus in tongue fur of gastric cancer patients were significantly lower than that of healthy patients, indicating the potential of oral microbial detection as a means of early diagnosis or screening of gastric cancer in the future. At present, the diagnosis of atrophic gastritis is mainly invasive endoscopic examination, which causes great pain and high cost for patients. If the change of the microbial community of tongue coating is used as a screening method for patients, it will bring important changes to the early diagnosis of atrophic gastritis.
There were still limitations in this study. The current sample size of this study is still small, and dietary habits, regional environment and other factors were not taken into account. Although the association can be derived from the results of the study, the cross-sectional design does not finally prove causality. Moreover, in this paper, the differences between the bacterial flora of the tongue coating were discussed, but the metabolites and pathogenic mechanisms of the specific bacterial flora were not further analyzed. The interrelation between functional roles of identified bacterial species and atrophic gastritis should be further elaborated in the follow-up studies.
5. Conclusion
We investigated that the tongue coating microbiota differed in atrophic gastritis patients compared with non-atrophic gastritis controls by 16S rRNA sequencing technology. The composition and diversity of the tongue coating microbiota differ at the class and genus between the two groups and functionally predicted metabolic pathways between the two groups were also discussed in our finding. We provide novel evidence that tongue coating microbiota may be a biomarker for characterizing patient with atrophic gastritis, but its mechanism needs to be further elaborated.
Funding Information
This work is supported by the Wenzhou traditional Chinese medicine construction project (Grant No. 2018008).