- 1Department of Traditional Chinese Medicine, Changhai Hospital, Naval Military Medical University, Shanghai, China
- 2Department of Gastroenterology, Shanghai Jiading Hospital of Traditional Chinese Medicine, Shanghai, China
Tongue diagnosis is a unique aspect of traditional Chinese medicine for diagnosing diseases before determining proper means of treatment, but it also has the disadvantage of relying on the subjective experience of medical practitioners and lack objective basis. The purpose of this article is to elucidate tongue-coating microbiota and metabolic differences in primary liver cancer (PLC) patients with thick or greasy tongue coatings. Tongue-coating samples were analyzed in 60 PLC patients (30 PLC with thick or greasy tongue-coating patients and 30 PLC with tongue-coating neither thick nor greasy) and 25 healthy controls (HC) using 16S rRNA gene sequencing technology. As compared to healthy individuals, tongue coatings of patients with PLC had elevated levels of Firmicutes and Actinobacteria. The abundance of Fusobacteria, SR1_Absconditabacteria_, and Spirochaete were higher in tongue coatings of healthy controls compared to samples in patients with PLC. In addition to site-specific differences, higher abundances of Fusobacteria and Actinobacteria were observed in thick or greasy tongue-coating patients as compared to non-thick and greasy tongue-coating patients. The inferred metagenomic pathways enriched in the PLC tongue-coating patients were mainly those involved in replication, recombination, and repair of protein. We also identify a tongue-coating microbiome signature to discriminate HC and PLC, including 15 variables on genus level. The prediction performance of the signature showed well in the training and validation cohorts. This research illustrates specific clinical features and bacterial structures in PLC patients with different tongue coatings, which facilitates understanding of the traditional tongue diagnosis.
Introduction
Primary liver cancer (PLC) is one of the most prevalent malignant tumors and the third leading cause of cancer-related death in the world (Sung et al., 2021), ranks fourth in the incidence of malignant tumors, and third in mortality in China (Chen et al., 2016). Histologically, PLC includes three main subtypes: hepatocellular carcinoma (HCC), intrahepatic cholangiocarcinoma (ICC), and mixed hepatocellular cholangiocarcinoma. Despite recent advances in treatments of PLC, including targeted and immune therapy, the length or quality of life of patients with HCC is still poor. In China, traditional Chinese medicine (TCM) is a wildly used therapy for the prevention and treatment of liver cancer (Xi and Minuk, 2018), and it has contributed to prolonging the survival and improving the quality of life of patients with PLC (Liu et al., 2019).
Tongue diagnosis is a novel diagnostic method in Chinese medicine to determine the health status of patients (Solos and Liang, 2018; Sun et al., 2018). TCM practitioners provide a basis for TCM diagnosis by observing the changes in the characteristics of the dorsum of the tongue, including the shape, size, texture, thickness, color, and so on (Anastasi et al., 2009). It will reveal the functional status of human organs and disease outcomes. Many features are used in the discrimination of TCM Syndromes, but the tongue-coating appearance is one of the major factors. For example, the greasy tongue is an important variable affecting discrimination in Shi-Re (dampness-heat) syndrome. The tongue diagnosis relies on the personal experience of Chinese medicine practitioners and the surrounding environment. However, it lacks certain reproducibility and is difficult to be recognized by the majority of modern doctors (Kamarudin et al., 2017). It is needed to provide a corresponding biological basis for tongue diagnosis in patients.
The oral cavity contains highly diverse microbiota, which includes bacteria, fungi, viruses, protozoa, and archaea. With the in-depth study of microorganisms, the influence of flora on human health and disease has been increasingly concerned by researchers. Many studies have found that the flora is closely related to colorectal cancer (Wong and Yu, 2019), pancreatic cancer (Riquelme et al., 2019), liver disease (Jia et al., 2019; Schwimmer et al., 2019), and other chronic diseases (Nishida et al., 2018; Dabke et al., 2019). The oral microbiome is the second largest microbial system in the human body, smaller than the intestinal microbial system (Human Microbiome Project Consortium, 2012). The oral cavity provides a suitable environment for microbial colonization due to its unique structure (Zhou and Zhang, 2015). The researchers found that it has a complex relationship with oral diseases (Kanasi et al., 2010; Shi et al., 2014) and systemic diseases (Pushalkar et al., 2012; Zhou et al., 2013). Recently, it was suggested that the microbiota dysbiosis of the tongue coat could facilitate the prognosis of several cancers (Ali Mohammed et al., 2021), such as pancreatic head carcinoma (Lu et al., 2019) and gastric cancer (Xu et al., 2021). But, the characteristics of the thick or greasy tongue-coating microbiota in patients with hepatocellular carcinoma are unclear.
Our findings would provide insight into the association between the oral microbiota and PLC, and the distinction between bacterial flora of thick or greasy tongue coating of PLC. It would provide the development of a new precautionary or diagnostic biomarker for patients with PLC.
Materials and methods
Study participants
The patients were recruited from the Changhai Hospital, China. From May 2017 to June 2018, after inclusion and exclusion criteria, 60 patients were diagnosed with PLC as well as 25 healthy subjects.
The inclusion criteria were as follows: (1) female or male, more than 18 years old; (2) histopathological or clinical confirmed PLC; (3) the following drugs are not used within 2 months: (a) antibiotics; (b) microbial preparations; and (c) other drugs that affect oral microbes; (4) all eligible individuals voluntarily agreed to participate and gave written informed consent; and (5) history of smoking cessation of more than 1 year. The exclusion criteria were as follows: (1) the healthy control group has a history of the gastric disease; (2) those who have used antibiotics, microbial preparations, or other drugs that affect oral microbes within 2 months; (3) receiving radiotherapy, radiofrequency, and percutaneous ethanol injection within 2 months; (4) serious infection; (5) those who have an oral disease; and (6) those who have other malignant tumors. The experiment was approved by the Ethics Committee of Changhai Hospital, and all participants signed informed consent.
Sampling collection
All tongue-coating samples were obtained and photographed in the morning before patients' food consumption to avoid the interference of food debris. Two traditional Chinese physicians with more than 20 years of clinical experience diagnosed the type of tongue coatings separately, and the DS01-A tongue diagnostic information acquisition system (DAOSH Co., Shanghai, China) was also used to photograph and analyze the images of the tongue coating. The patients were only recruited when the diagnosis of two physicians and the analysis of DS01-A were consistent. All participants were required to rinse their mouths by gargling sterile saline two times before sampling. Each tongue-coating sample was collected from the middle section of the tongue dorsum using a fresh one-off toothbrush and put into the test tube with saline. The tubes were centrifuged for 10 min at 3,000 r/min, and the precipitates were collected. Samples were stored at – 80°C.
DNA extraction and PCR amplification
Total DNA was extracted using the E.Z.N.A.® Soil DNA Kit (Omega Bio-tek, Norcross, GA, U.S.) according to the manufacturer's protocols. The integrity of the genomic DNA was assessed by electrophoresis (1% agarose gel). Polymerase chain reaction (PCR) amplification of the V3–V4 region of the bacterial 16S rRNA genes was performed using forward primers 338F (5'-ACTCCTACGGGAGGCAGCAG-3') and 806R (5'-GGACTACHVGGGTWTCTAAT-3'). The PCR reactions were performed in 20 μl PrimerStar HS Premix (AP221-02; TransGen, Beijing, China) that contained 4 μl of 5×FastPfu Buffer, 2 μl of 2.5 mM dNTPs, 0.4 μl of forward primer, 0.4 μl of reverse primer, 0.4 μl of FsatPfu Polymerase, and 10 ng DNA. All samples were amplified on an ABI GeneAmp 9700 (ABI, USA) using the following parameters: 95°C for 2 min, followed by 25 cycles of 95°C for 30 s, 55°C for 30 s, and 72°C for 30 s, and a final extension at 72°C for 5 min. The obtained PCR products were run on 2% agarose gel and purified after size selection repeated three times for each sample and quantified using the QuantiFluor™ ST system.
Sequence analysis
The PCR product was recovered using a 2% agarose gel, purified using AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, USA), eluted with Tris-HCl, and detected by 2% agarose electrophoresis. Detection quantification was performed using QuantiFluorTM-ST (Promega, USA). The purified amplified fragment was constructed into a library of PE 2*300 according to the standard operating protocol of the Illumina MiSeq platform (Illumina, San Diego, USA) and Shanghai Majorbio Bio-pharm Technology Co., Ltd.
Bioinformatics analysis
Raw fastq files were demultiplexed, quality-filtered by Trimmomatic, and merged by Fast Length Adjustment of SHort reads (FLASH) v1.2.10 software with the following criteria: (1) The reads were truncated at any site receiving an average quality score of <20 over a 50-bp sliding window; (2) Primers were exactly matched allowing two nucleotide mismatching, and reads containing ambiguous bases were removed; And (3) Sequences whose overlap is longer than 10 bp were merged according to their overlap sequence.
Operational taxonomic units (OTUs) were clustered with a 97% similarity cutoff using UPARSE (version 7.1 http://drive5.com/uparse/), and chimeric sequences were identified and removed using UCHIME. The taxonomy of each 16S rRNA gene sequence was analyzed by the RDP Classifier algorithm (http://rdp.cme.msu.edu/) against the Silva (SSU123) 16S rRNA database using a confidence threshold of 70%.
Statistical analysis
The alpha diversity, the Relative abundances, and the linear discriminant analysis (LDA) of the effect size (LEfSe) (http://huttenhower.sph.harvard.edu/lefse/) were performed by R software packages. The Wilcoxon rank sum test was used for the differential analysis of continuous variables between the two groups. One-way ANOVA was used to assess the differences between the three groups. The Chi-square test was used for counting data tests. The Chi-square test and t-test were used to compare the common characteristics of participants and the non-parametric test was used to compare the indices of alpha diversity and the relative abundances between different types of tongue coating. Correlate analysis was used to establish the commensal relationship. Receiver operating characteristic (ROC) curves were used to analyze the distinguishing ability of the taxa. P-value < 0.05 was statistically significant. Statistical analysis was performed using R software 3.6.0.
Results
Participant characteristics
A total of 60 PLC patients and 25 unrelated control subjects were assigned to three groups: (1) PLC with thick or greasy coating (LTG) group (n = 30) was comprised of 27 men and 3 women with an average of 57.7 ± 10 years old; (2) PLC with non-thick or greasy coating (LNTG) group (n = 30) had 26 men and 4 women with an average of 53.67 ± 11 years old; and (3) The healthy control (HC) group (n = 25) had 23 men and 2 women with an average of 56.68 ± 10.359 years (Table 1 and Figure 1). Age, gender, Body mass index (BMI), diabetes, and high blood pressure were comparable between the three groups (P > 0.05).
Figure 1. Different tongue-coating images of the three groups. (A) Thick or greasy coating (LTG) group (abnormal). (B) Non-thick or greasy coating (LNTG) group (abnormal). (C) Healthy control (HC) group (normal).
Tongue-coating microbial dysbiosis
The Venn diagrams were performed to present the species diversity among the three groups and between the two groups, respectively (Figures 2A,B). There were 389 species in HC and PLC, including 7 species unique to HC and 65 species unique to PLC (Figure 2A). There are 373 species in the three groups; 7 species are not detected in the LTG and LNTG groups, 3 species are unique to the LTG group, and 41 species are unique to the LNTG group (Figure 2B). It was suggested that the species composition among the groups is different and the LNTG group has unique species. Moreover, the species of tongue coating in patients with liver cancer are more than those in the HC group.
Figure 2. The tongue-coating microbial dysbiosis. (A) Venn diagram of the shared OTUs between HC and PLC groups in tongue-coating contents. (B) Venn diagram of the shared OTUs among HC, LTG, and LNTG groups in tongue-coating contents. (C,D) The relative abundance of tongue-coating microbiota at the phylum level in the two different groups (C) and the three different groups (D). (E,F) The heatmap of the bacterial community composition at genus level across the two different groups (E) and the three different groups (F). (PLC represents primary liver cancer, LTG represents liver cancer with thick or greasy coating, LNTG represents liver cancer with non-thick or greasy coating, and HC represents healthy control).
The differential abundance of OTUs at the phylum, and genus levels of bacterial taxonomic classification reveal differences in the relative abundance of several taxa among the three groups and between the two groups, respectively (Figures 2C–F). At the phylum level, the community composition of PLC and HC groups includes Bacteroidetes, Firmicutes, Proteobacteria, Fusobacteria, Actinobacteria, Saccharibacteria, and so on. PLC participants show a slightly higher abundance of Firmicutes and Actinobacteria and a slightly lower abundance of phyla Proteobacteria and Fusobacteria (Figure 2C). LTG participants show a slightly higher abundance of Bacteroidetes, Fusobacteria, and Actinobacteria, and a slightly lower abundance of phyla Firmicutes and Proteobacteria (Figure 2D). At the genus level, PLC participants showed a slightly higher abundance of Prevotella_7, Streptococcus, Veillonella, and Actinomyces, and a slightly lower abundance of Neisseria and Fusobacterium (Figure 2E). LTG participants showed a slightly higher abundance of Prevotella_7, Actinobacteria, and Veillonella than those of the HC group (Figure 2F). The major bacteria were principally consistent among these groups, but different relative abundances could be observed (Figures 2E,F).
The tongue-coating microbial diversities
Similarities of the bacterial communities among the three groups and between the two groups were compared by NMDS and ANOSIM based on Bray–Curtis. NMDS ordination plots of the tongue-coating data showed a clear distinction between the two groups (Figure 3A; NMDS, Stress = 0.181) and among the three groups (Figure 3B; NMDS, Stress = 0.181). ANOSIM revealed significant differences in the structure of tongue-coating microbiota between the two groups (Figure 3C; ANOSIM, r = 0.191, P < 0.05) and among the three groups (Figure 3D; ANOSIM, r = 0.0756, P < 0.05). At the phylum level, the Wilcoxon rank-sum test was performed among the three groups and between the two groups. PLC participants showed a higher abundance of Firmicutes and Actinobacteria and a slightly lower abundance of phylum Fusobacteria and SR1_Absconditabacteria_, Spirochaete (all P < 0.05; Figure 3E). LTG participants show a slightly higher abundance of Fusobacteria and Actinobacteria and a slightly lower abundance of phylum Firmicutes (all P < 0.05; Figure 3F).
Figure 3. The tongue-coating microbial diversities. (A,B) Beta diversity of the tongue-coating microbiota was analyzed using non-metric multidimensional scaling (NMDS) on a Brays–Curtis distance matrix between the two groups (A) and among the three groups (B). (C,D) Box diagram of the tongue-coating microbiota was analyzed using the ANOSIM test on a Bray–Curtis distance matrix between the two groups (C) and among the three groups (D). (E,F) Comparison of the microbiome on phylum level of the two groups (E) and among the three groups (F).
Differences in tongue-coating bacterial community composition
The differentially abundant bacterial taxa between the two groups and among the three groups were identified in the tongue-coating microbiota by the linear discriminant analysis effect size (LEfSe) (Figure 4). Between the two groups on genus level, the differentially enriched taxa in PLC groups were mainly from Firmicutes, Prevotellaceae, c_Negativicutes, o_selenomonadales, Veillonellaceae, and Veillonella, whereas the differentially enriched taxa in HC groups were mainly from Fusobacteriaceae, Fusobacterium, Fusobacteriales, Fusobacteria, and porphyromonas (Figures 4A,B). Among the three groups on genus level, the differentially enriched taxa in LNTG groups were mainly from Firmicutes and the differentially enriched taxa in LTG groups were mainly from c_Negativicutes, o_selenomonadales, Veillonellaceae, and Prevotella_7 (Figures 4C,D).
Figure 4. Differentially abundant bacterial taxa in the tongue-coating microbiota were identified by LDA effect size (LEfSe) analysis. (A,B) Taxonomic cladogram (A) of significant differences and histogram (B) of the LDA scores for PLC (blue) and HC (red) groups. (C,D) Taxonomic cladogram (C) of significant differences and histogram (D) of the LDA scores between HC (red), LTG (green), and LNTG (blue) groups. (PLC represents primary liver cancer, LTG represents liver cancer with thick or greasy coating, LNTG represents liver cancer with non-thick or greasy coating, and HC represents healthy control).
Identification of tongue-coating microbiome signature to discriminate HC and PLC
To explore the potential ability of the gut microbiome to classify HC and PLC, we constructed a tongue-coating microbiome signature by a random forest model. The two-dimensional scatter graph represents the similarity of samples in different groups (Figure 5A). The number of top important variables in the training set was 15 calculated by a genus-level random forest model with a decision tree for 500 times (Figure 5B). The signature consisted of top 15 variable importance on genus level, including Fusobacterium, Bergeyella, norank_c__JG37_AG_4, Megasphaera, unclassified_k__norank, Prevotella_6, norank_f__Actinomycetaceae, Moraxella, Atopobium, Mogibacterium, norank_f__Leptotrichiaceae, Veillonella, Peptococcus, Solobacterium, norank_f__Lachnospiraceae, and p__Bacteroidetes (Figure 5C). The area under the ROC curve (AUG) in the training (15 HC and 25 PLC individuals; Figure 5D) and validating (10 HC and 10 PLC individuals; Figure 5E) set were 0.63 (95% confidential interval (CI): 0.45–0.8) and 0.72 (95% CI: 0.47–0.97) based on the 15 most important variables. The 15-microbiome signature showed great diagnostic power to discriminate PLC patients from HC individuals.
Figure 5. Identification of a tongue-coating microbiome signature to discriminate HC and PLC by random forest analysis. (A) The two-dimensional scatter graph represents the similarity of samples in different groups. (B) Signature evaluation by using top 10 important features on genus level by random forest analysis. (C) Dot chart of top 15 variable importance on genus level. (D) ROC of the signature is based on the 15 most important variables in the training set. (E) ROC of the signature based on the 15 most important variables in the validation set. (PLC represents primary liver cancer, HC represents healthy control).
Functional prediction
The functional prediction was analyzed to reveal potential function annotation of the Cluster of orthologous groups (COG) family and KEGG Ortholog (KO) using Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt). The top five most abundant functions in HC, LTG, and LNTG groups were: translation, ribosomal structure, and biogenesis; general function prediction; Replication, recombination, and repair; cell wall/membrane/envelope biogenesis; and amino acid transport and metabolism (Figure 6A). The same results of COG function prediction were found between HC and PLC groups (Figure 6B).
Figure 6. Functional predictions. (A) The compositions of COG function in the HC and PLC groups. (B) The compositions of COG function in the HC, LTG, and LNTG groups. (PLC represents primary liver cancer, LTG represents liver cancer with thick or greasy coating, LNTG represents liver cancer with non-thick or greasy coating, and HC represents healthy control).
Discussion
Nowadays, PLC still represents a heavy burden for patients and healthcare systems in China, due to high incidence and high mortality. TCM is widely used in all stages of PLC that plays an active role in the prevention and treatment of PLC in China (Hamasaki, 2017; Shi et al., 2017). The TCM Syndrome, known as “Zheng,” is an integral and essential part of TCM theory. It is always used to identify the pattern of patients and guide individualized therapy (Chen and Wang, 2012; Ji et al., 2016). Many features are used in the discrimination of syndromes, but the tongue-coating appearance is one of the major factors. Shi-Re (dampness-heat) syndrome is one of the fundamental syndromes of malignancies and a greasy or thick tongue coating usually indicates dampness-heat syndrome (Sui et al., 2019). In summary, to investigate the underlying mechanisms of greasy or thick tongue coating, we can explore the nature of dampness-heat syndrome and obtain new clues and novel ideas for the objective studies of TCM syndromes.
The main purpose of this part of the study is to try to find evidence of tongue microbiota for greasy or thick tongue coating of PLC and to provide a reference to distinguish PLC from healthy individuals. It was revealed that the relative abundance of Actinobacteria, Veillonella, Actinomyces, and Prevotella_7 in the LTG group was significantly higher than that in HC and LNTG groups. For patients with liver cancer, we found that the relative abundance of Firmicutes, Actinobacteria, Veillonella, and Actinomyces in the PLC group was significantly higher than that in the HC group. However, Li et al. found that the dominant bacteria on the tongue coating of patients with liver cancer were Fusobacterium (Lu et al., 2016), while we found that Fusobacterium was highly enriched in patients with HC, which was contrary to this conclusion. We found that Firmicutes, Actinomyces, and Actinomyces were more abundant in tongue coating in the PLC group, which was consistent with Li et al.'s conclusion. It still needs further in-depth mechanism research to elaborate the results.
Through the discriminant analysis of LEfSe multistage species difference, we found that Firmicutes, Prevotellaceae, Veillonellaceae, and Veillonella in the PLC group had a great influence on the differential effect. However, Veillonellaceae and Prevotella_7 had a great influence on the differential effect in the LTG group. Our function prediction of tongue coating showed that the functions enriched in the HC group mainly included replication, recombination, and repair, while the function in the PLC group was on the low side. The PLC group is mainly enriched in the pentose phosphate pathway and lysine biosynthesis. It was found that the main function of this group may be to induce the risk of the inflammatory reaction to more inflammatory factors.
Using the 16S rRNA gene, the core function of tongue-coating microbiome could be predicted by the PICRUSt software. Pathways encoding for carbohydrate metabolism and oxidative phosphorylation metabolism were detected. The oral cavity is a major gateway to the human body. Microorganisms colonizing the oral cavity have a significant probability of spreading to the stomach, lung, and intestinal tract. The metagenome prediction validated the role of tongue coating.
Recently, microbiome signatures are presented to discriminate against different diseases, including colorectal cancer (Mangifesta et al., 2018) and gastric cancer (Xu et al., 2019; Wang et al., 2020). These results showed that the identified microbiome signatures could be non-invasive and accurate biomarkers for early diagnosis of cancers. No study to date, however, has investigated such microbiome signature in PLC in the literature. In this study, the microbiome signature to discriminate HC and PLC was identified using the random forest analysis. The signature consisted of a set of 15 bacterial genera capable of PLC with considerably high accuracy. We found that the signatures classified PLC from healthy controls with an area under the curve (AUC) values all >0.6 in the training and validating sets. The results indicate that the constructed signature is powerful to predict the diagnosis of PLC.
There are still some defects in the study. First, the limited sample size of our study might have led to a certain degree of bias in the results. Second, we only analyzed the specificity of tongue-coating microbiome in patients with the damp-heat syndrome of liver cancer and patients with liver cancer, and predict the differences in microbial function by function analysis. Third, we lacked an additional set of tongue-coating microbiome signatures for external validation. There is no further study on how differential microbial communities cause the damp-heat syndrome of liver cancer. High-throughput sequencing is particularly important in the tongue-coating flora.
In summary, we provided biological evidence of tongue diagnosis and dampness-heat syndrome in TCM. We identified a tongue-coating microbiome signature that was capable of accurately distinguishing PLC from healthy controls. The signature showed characteristics of varying abundance, high degree of bacteria interaction, and carcinogenic potentials. It could serve as a potential non-invasive biomarker, which might be suitable for long-term monitoring of PLC.
Data availability statement
The raw sequence data reported have been deposited in the Genome Sequence Archive (Genomics, Protemics & Bioinformatics 2021) in National Genomics Data Center (Nucleic Acids Res 2022), China National Center for Bioinformation/Beijing Institute of Genomics, Chinese Academy of Sciences (GSA: CRA009244) that are publicly accessible at https://ngdc.cncb.ac.cn/gsa.
Ethics statement
The studies involving human participants were reviewed and approved by the Ethics Committee of Changhai Hospital. The patients/participants provided their written informed consent to participate in this study.
Author contributions
YZ and XZ conceived and designed the experiments, edited, and reviewed the manuscript. YG, YM, SY, and BP performed the experiments. HZ wrote the original draft. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Clinical innovation project of Changhai Hospital (2019YXK029).
Acknowledgments
We thank the Clinical innovation project of Changhai Hospital (2019YXK029). We also thank Majorbio Bio-Pharm Technology Co. Ltd. (Shanghai, China); analysis of amplicon data was performed through the online Majorbio Cloud Platform (www.majorbio.com).
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.
References
Ali Mohammed, M. M., Al Kawas, S., and Al-Qadhi, G. (2021). Tongue-coating microbiome as a cancer predictor: a scoping review. Archives Oral Biol. 132, 105271. doi: 10.1016/j.archoralbio.2021.105271
Anastasi, J. K., Currie, L. M., and Kim, G. H. (2009). Understanding diagnostic reasoning in TCM practice: tongue diagnosis. Alternat. Therapies Health Med. 15, 18–28.
Chen, W., Zheng, R., Baade, P. D., Zhang, S., Zeng, H., Bray, F., et al. (2016). Cancer statistics in China, 2015. CA Cancer J. Clin. 66, 115–132. doi: 10.3322/caac.21338
Chen, Z., and Wang, P. (2012). Clinical distribution and molecular basis of traditional Chinese medicine Zheng in cancer. Evidence Based Complement. Alternat. Med. eCAM 2012, 783923. doi: 10.1155/2012/783923
Dabke, K., Hendrick, G., and Devkota, S. (2019). The gut microbiome and metabolic syndrome. J. Clin. Invest. 129, 4050–4057. doi: 10.1172/JCI129194
Hamasaki, H. (2017). Exercise and gut microbiota: clinical implications for the feasibility of Tai Chi. J. Integrat. Med. 15, 270–281. doi: 10.1016/S2095-4964(17)60342-X
Human Microbiome Project Consortium. (2012). A framework for human microbiome research. Nature 486, 215–221. doi: 10.1038/nature11209
Ji, Q., Luo, Y. Q., Wang, W. H., Liu, X., Li, Q., Su, S. B., et al. (2016). Research advances in traditional Chinese medicine syndromes in cancer patients. J. Integrat. Med. 14, 12–21. doi: 10.1016/S2095-4964(16)60237-6
Jia, X., Lu, S., Zeng, Z., Liu, Q., Dong, Z., Chen, Y., et al. (2019). Characterization of gut microbiota, bile acid metabolism, and cytokines in intrahepatic cholangiocarcinoma. Hepatology 71, 893–906. doi: 10.1002/hep.30852
Kamarudin, N. D., Ooi, C. Y., Kawanabe, T., Odaguchi, H., and Kobayashi, F. (2017). A fast SVM-based tongue's colour classification aided by k-means clustering identifiers and colour attributes as computer-assisted tool for tongue diagnosis. J. Healthcare Eng. 2017, 7460168. doi: 10.1155/2017/7460168
Kanasi, E., Dewhirst, F. E., Chalmers, N. I., Kent, R., Moore, C. V., Hughes, N., et al. (2010). Clonal analysis of the microbiota of severe early childhood caries. Caries Res. 44, 485–497. doi: 10.1159/000320158
Liu, X., Li, M., Wang, X., Dang, Z., Yu, L., Wang, X., et al. (2019). Effects of adjuvant traditional Chinese medicine therapy on long-term survival in patients with hepatocellular carcinoma. Phytomedicine Int. J. Phytotherapy Phytopharmacol. 62, 152930. doi: 10.1016/j.phymed.2019.152930
Lu, H., Ren, Z., Li, A., Li, J., Xu, S., Zhang, H., et al. (2019). Tongue coating microbiome data distinguish patients with pancreatic head cancer from healthy controls. J. Oral Microbiol. 11, 1563409. doi: 10.1080/20002297.2018.1563409
Lu, H., Ren, Z., Li, A., Zhang, H., Jiang, J., Xu, S., et al. (2016). Deep sequencing reveals microbiota dysbiosis of tongue coat in patients with liver carcinoma. Sci. Rep. 6, 33142. doi: 10.1038/srep33142
Mangifesta, M., Mancabelli, L., Milani, C., Gaiani, F., de'Angelis, N., de'Angelis, G. L., et al. (2018). Mucosal microbiota of intestinal polyps reveals putative biomarkers of colorectal cancer. Sci. Rep. 8, 13974. doi: 10.1038/s41598-018-32413-2
Nishida, A., Inoue, R., Inatomi, O., Bamba, S., Naito, Y., Andoh, A., et al. (2018). Gut microbiota in the pathogenesis of inflammatory bowel disease. Clin. J. Gastroenterol. 11, 1–10. doi: 10.1007/s12328-017-0813-5
Pushalkar, S., Ji, X., Li, Y., Estilo, C., Yegnanarayana, R., Singh, B., et al. (2012). Comparison of oral microbiota in tumor and non-tumor tissues of patients with oral squamous cell carcinoma. BMC Microbiol. 12, 144. doi: 10.1186/1471-2180-12-144
Riquelme, E., Zhang, Y., Zhang, L., Montiel, M., Zoltan, M., Dong, W., et al. (2019). Tumor microbiome diversity and composition influence pancreatic cancer outcomes. Cell 178, 795–806.e712. doi: 10.1016/j.cell.2019.07.008
Schwimmer, J. B., Johnson, J. S., Angeles, J. E., Behling, C., Belt, P. H., Borecki, I., et al. (2019). Microbiome signatures associated with steatohepatitis and moderate to severe fibrosis in children with nonalcoholic fatty liver disease. Gastroenterology 157, 1109–1122. doi: 10.1053/j.gastro.2019.06.028
Shi, H., Rui, L., Xiaowei, Z., Tao, H., Helen, Z., Alice, C., et al. (2014). Predictive modeling of gingivitis severity and susceptibility via oral microbiota. ISME J. 8, 1768. doi: 10.1038/ismej.2014.32
Shi, Z., Song, T., Wan, Y., Xie, J., Yan, Y., Shi, K., et al. (2017). A systematic review and meta-analysis of traditional insect Chinese medicines combined chemotherapy for non-surgical hepatocellular carcinoma therapy. Sci. Rep. 7, 4355. doi: 10.1038/s41598-017-04351-y
Solos, I., and Liang, Y. (2018). A historical evaluation of Chinese tongue diagnosis in the treatment of septicemic plague in the pre-antibiotic era, and as a new direction for revolutionary clinical research applications. J. Integrat. Med. 16, 141–146. doi: 10.1016/j.joim.2018.04.001
Sui, X., Guo, Y., Ni, W., Jin, H., Lin, H., Xie, T., et al. (2019). Molecular profiling analysis for colorectal cancer patients with Pi-Xu or Shi-Re syndrome. Integrat. Med. Res. 8, 21–25. doi: 10.1016/j.imr.2018.11.001
Sun, S., Wei, H., Zhu, R., Pang, B., Jia, S., Liu, G., et al. (2018). Biology of the tongue coating and its value in disease diagnosis. Complement. Med. Res. 25, 191–197. doi: 10.1159/000479024
Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global cancer statistics 2020: globocan estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clinicians 71, 209–249. doi: 10.3322/caac.21660
Wang, L., Xin, Y., Zhou, J., Tian, Z., Liu, C., Yu, X., et al. (2020). Gastric Mucosa-associated microbial signatures of early gastric cancer. Front. Microbiol. 11, 1548. doi: 10.3389/fmicb.2020.01548
Wong, S. H., and Yu, J. (2019). Gut microbiota in colorectal cancer: mechanisms of action and clinical applications. Nat. Rev. Gastroenterol. Hepatol. 16, 690–704. doi: 10.1038/s41575-019-0209-8
Xi, S. Y., and Minuk, G. Y. (2018). Role of traditional Chinese medicine in the management of patients with hepatocellular carcinoma. World J. Hepatol. 10, 799–806. doi: 10.4254/wjh.v10.i11.799
Xu, J., Xiang, C., Zhang, C., Xu, B., Wu, J., Wang, R., et al. (2019). Microbial biomarkers of common tongue coatings in patients with gastric cancer. Microbial Pathogenesis 127, 97–105. doi: 10.1016/j.micpath.2018.11.051
Xu, S., Xiang, C., Wu, J., Teng, Y., Wu, Z., Wang, R., et al. (2021). Tongue coating bacteria as a potential stable biomarker for gastric cancer independent of lifestyle. Digestive Dis. Sci. 66, 2964–2980. doi: 10.1007/s10620-020-06637-0
Zhou, M., Rong, R., Munro, D., Zhu, C., Gao, X., Zhang, Q., et al. (2013). Investigation of the effect of type 2 diabetes mellitus on subgingival plaque microbiota by high-throughput 16S rDNA pyrosequencing. PLoS ONE 8, e61516. doi: 10.1371/journal.pone.0061516
Keywords: traditional Chinese medicine (TCM), 16S rRNA gene sequencing, tongue coating microbiota, primary liver cancer, syndrome
Citation: Zhang Y, Zhao H, Guo Y, Meng Y, Yu S, Pan B and Zhai X (2022) Relationship between thick or greasy tongue-coating microbiota and tongue diagnosis in patients with primary liver cancer. Front. Microbiol. 13:903616. doi: 10.3389/fmicb.2022.903616
Received: 25 March 2022; Accepted: 20 September 2022;
Published: 29 December 2022.
Edited by:
Ana Isabel Alvarez- Mercado, University of Granada, SpainReviewed by:
Wensheng Lan, Technology Center for Animal and Plant Inspection and Quarantine, Shenzhen Customs, ChinaJulio Plaza-Diaz, Children's Hospital of Eastern Ontario (CHEO), Canada
Copyright © 2022 Zhang, Zhao, Guo, Meng, Yu, Pan and Zhai. 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: Xiaofeng Zhai, emhhaXhmY2gmI3gwMDA0MDsxNjMuY29t
†These authors have contributed equally to this work and share first authorship