Characterization of the Esophageal Microbiota and Prediction of the Metabolic Pathways Involved in Esophageal Cancer

Esophageal microbiota plays important roles in esophageal cancer. Esophagectomy, as the most important therapeutic way, contributes to changes of esophageal microbiome. However, there are few studies examining the esophageal microbiome and the metabolic changes before and after esophagectomy. The present study characterized the esophageal microbiome of 17 patients with esophageal squamous cell carcinoma (ESCC), 11 patients with esophagogastric junction (EGJ) cancer, 15 patients at 9–12 months after radical esophagectomy and 16 healthy controls (HC). 16S ribosomal RNA gene sequencing was used to evaluate the microbiome and predict the metabolic pathways. Our results showed that the microbial diversity was significantly lower in ESCC, EGJ and post-ESCC groups than that in the HC group. The abundance of Fusobacteria was higher (7.01 vs. 1.12%, P = 0.039) and the abundance of Actinobacteria (1.61 vs. 4.04%) was lower in the ESCC group than that in the HC group. We found significant differences in the abundance of Bacteroidetes (20.45 vs. 9.86%, P = 0.026), Fusobacteria (7.01 vs. 1.66%, P = 0.030) between ESCC and post-ESCC groups. The results of microbial composition analysis and PICRUSt demonstrated significant differences between ESCC and HC groups. The β diversity and PICRUSt suggested that the microbial composition and metabolic pathways were similar to HC group after esophagectomy. The monitoring of the esophagus microbiota may be an essential method to predict the recurrence of tumor.

Introduction

Esophageal cancer is a rapidly growing concern worldwide, with ~572,000 new cases annually and 509,000 deaths in 2018 (Bray et al., 2018). Esophageal squamous cell carcinoma (ESCC) is the most common histological subtype of esophageal cancer, particularly in areas of eastern Asia and eastern and southern Africa. Approximately 90% of esophageal cancer cases are ESCC in China, where the disease is a major public health problem (Arnold et al., 2015; Fitzmaurice et al., 2015). The incidence of adenocarcinoma of the esophagogastric junction (EGJ) also increased rapidly over the last few decades, accounting for one-third of all esophagogastric adenocarcinoma cases (Wu et al., 2009).

Esophageal cancer is highly invasive, and it has a poor prognosis, with a five-year survival rate of ~30% in China (Bray et al., 2018). The reported risk factors of esophageal cancer included drinking, smoking, ingestion of hot food, obesity, gastroesophageal reflux disease, and Barrett's esophagus (BE) (Bray et al., 2018; Ferlay et al., 2018). The gut microbiota has been shown to play an essential role in several cancers, including gastric cancer (Brawner et al., 2014), colon cancer (Castellarin et al., 2012; Yang et al., 2017), and pancreatic cancer (Riquelme et al., 2019). As an important part of the upper gastrointestinal tract, the esophageal mucosa is also colonized by microbes. Previous studies have reported that an imbalance of the microbiome may promote esophageal cancer development and progression. Some studies have shown lower microbial richness in the upper digestive tract to be associated with esophageal squamous dysplasia, regarded as a precursor lesion of ESCC (Yu et al., 2014). The Fusobacterium nucleatum and Porphyromonas gingivalis were associated with shorter survival and might contribute to aggressive tumor behavior through the activation of chemokines in ESCC patients (Bao et al., 2014; Yamamura et al., 2019). The microbiota of the oral cavity (Peters et al., 2017) and gastric cancer (Nasrollahzadeh et al., 2015) from patients with ESCC also revealed a decrease in overall oral microbial diversity and enrichment in Clostridiales and Erysipelotrichales in the gastric corpus of patients with ESCC. Microbiota dysbiosis, including the presence of Veillonella, Prevotella, Haemophilus, Neisseria, Campylobacter, and Fusobacterium, has also been reported in association with gastroesophageal reflux disease (GRED) and is hypothesized to contribute to the evolution toward BE and adenocarcinoma at the esophagogastric junction (EGJ) (Di Pilato et al., 2016). Although a few studies have suggested the effect of esophageal microbiota on ESCC and EGJ, further data that can reveal the microbial composition in ESCC and EGJ are needed.

The treatment of ESCC includes many choices, such as endoscopic treatment, oncological treatment, surgery, and chemoradiotherapy (Lagergren et al., 2017). Esophagectomy, as one of the most radical therapeutic methods of esophageal cancer, can achieve the goal of a radical cure for early-stage patients and extend overall survival time for middle- and advanced-stage patients. Gastroesophageal reflux has been a long-standing complication after esophagectomy due to the resection of cardia. Previous studies have shown that bile reflux after subtotal gastrectomy was associated with the presence of Streptococcus and Veillonella in gastric aspirates and Escherichia, Klebsiella, and Clostridium in the intestine (Tseng et al., 2016). However, there was no study reporting a microbial characterization of ESCC patients receiving esophagectomy. We plan to investigate the composition of the esophageal microbiota in healthy tissues, tumor tissues, and after esophagectomy in the present study.

Materials and Methods

Study Population and Sample Collection

Normal esophagus, esophageal squamous cell carcinoma (ESCC), esophagogastric junction cancer (EGJ, in accordance with the 8th Edition of the AJCC TNM Classification), and postoperative esophageal squamous cell carcinoma (post-ESCC) specimens were collected under electronic gastroscopy from January 2018 to June 2019 at Renmin Hospital, Wuhan University (Wuhan, China). Healthy controls (HC) were normal esophageal (upper, middle, and lower segments) specimens obtained from 16 healthy volunteers with no digestive symptoms and esophagogastric mucosal lesions, as confirmed by electronic gastroscopy. ESCC (n = 17) and EGJ (n = 11) specimens were harvested from 28 patients with primary tumors, and all EGJ cases were adenocarcinoma. Post-ESCC (n = 15) specimens were obtained from the esophageal stump (1.5 cm above the anastomosis site of the esophageal stump and stomach) at 9–12 months after radical esophagectomy without chemoradiotherapy. All patients were diagnosed by pathological examinations.

The exclusion criteria were as follows: a systemic infectious disease, other coexisting malignant tumors, preoperative neoadjuvant chemoradiotherapy, biotherapy, or a history of gastrointestinal surgery. Patients using antibiotics and microecological preparations (probiotics, prebiotics, or symbiotics) within 2 months were also excluded. Demographic data and general clinical information related to the microecological analysis were acquired, and written informed consent was obtained from each participant. The experiments were approved by the Institutional Review Board of the Renmin Hospital, Wuhan University.

DNA Extraction, Amplification, and Sequencing

Microbial genomic DNA was extracted from normal esophagus, ESCC, EGJ, and post-ESCC specimens using the E.Z.N.A.^® Soil DNA Kit (Omega Bio-Tek, Norcross, GA, USA). The quantity and purity of DNA were determined with a NanoDrop 2000 UV-Vis Spectrophotometer (Thermo Scientific, Wilmington, MA, USA). The hypervariable region V3-V4 directionally targeted by the bacterial 16S rRNA gene was amplified with primers (338F: 5′– ACTCCTACGGGAGGCAGCAG−3′; 806R: 5′–GGACTACHVGGGTWTCTAAT−3′) in a GeneAmp 9700 PCR System (Applied Biosystems, Foster City, CA, USA) (Edgar, 2013). The PCR products were purified using the AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, USA) and quantified using QuantiFluor™-ST (Promega, USA) according to the manufacturers' instructions. Purified amplicons were pooled in equimolar concentrations, and paired-end high-throughput sequencing was performed using a 2 × 300 kit on the Illumina MiSeq Platform (Illumina, San Diego, CA, USA). The protocol of Majorbio Bio-Pharm Technology Co, Ltd. (Shanghai, China) was followed.

Data Processing and Statistical Analysis

Raw Fastq files were merged in FLASH software (Magoc and Salzberg, 2011) and quality-filtered with Trimmomatic (Bolger et al., 2014) using the following criteria (Schirmer et al., 2015): (i) reads were truncated at any site receiving an average quality score <20 over a 50 bps sliding window, which retained sequences with overlaps >10 bps and mismatches of no more than 2 bp, and (ii) sequence data were demultiplexed and assigned to samples based on barcodes (exact matches) and primers (two nucleotide mismatches were allowed). Reads containing ambiguous bases were removed.

Operational taxonomic units (OTUs) were clustered with a 97% sequence similarity cutoff using UPARSE (version 7.1 https://drive5.com/uparse/) (McDonald et al., 2012; Edgar, 2013) and a novel “greedy” algorithm that simultaneously performed chimera filtering and OTU clustering. The taxonomic assignment of each sequence was carried out using the Ribosomal Database Project (RDP) Classifier algorithm (http://rdp.cme.msu.edu/) against the Silva (SSU123) 16S rRNA database with a confidence threshold of 70% (Liu et al., 2018).

Statistical analysis was carried out using GraphPad Prism (version 8.0.2) and the Majorbio I-Sanger Cloud Platform (http://www.i-sanger.com). The Kruskal–Wallis test, Fisher's exact test, and paired t-test were used to compare demographic characteristics. Alpha diversity indexes were calculated using MOTHUR (version 1.30.1) (Schloss et al., 2009). Principal coordinates analysis (PCoA) based on the weighted and unweighted UniFrac distance and analysis of similarities (ANOSIM) was carried out to compare the global microbial composition at the operational taxonomic unit (OTU) level (Kageyama et al., 2019). Welch's t-test and the Wilcoxon rank-sum test were used to identify species at phylum and genus levels. The detection of discriminant bacterial species was performed using the linear discriminant analysis effect size (LEfSe) (Segata et al., 2011). The linear discriminant analysis score (LDA score) indicated the effect size of each OTU, and OTUs with an LDA score >3.0 were defined as differentially abundant OTUs (Yang et al., 2018; Kageyama et al., 2019). Phylogenetic investigation of communities by reconstruction of unobserved states (PICRUSt) (Langille et al., 2013) analysis was used to identify Kyoto Encyclopedia of Genes and Genomes (KEGG) biochemical pathways, and the results were visualized as a heatmap with the Multiple Experiment Viewer (version 4.9.0).

Results

Demographic Characteristics of All Individuals

The demographic characteristics of all individuals (n = 59) are shown in Table 1. No differences were observed in age, sex, alcohol intake, tobacco smoking, and TNM stage. The percentages of II and III patients in the EGJ, ESCC, and post-ESCC groups were 90.9, 88.2, and 86.7%, respectively. It is of note that gastroesophageal refluxing, which is the most common complication after esophagectomy, existed in 11 (73.3%) patients of the post-ESCC group but was not observed in other groups.

TABLE 1

Table 1. Demographic characteristics of all individuals.

Decreased Microbial Diversity and Richness in the ESCC, EGJ, and Post-ESCC Groups Compared With the HC Group

For community coverage, the value of the Good's coverage for all samples was >99% (Figure S1A), indicating that the sequencing results reflected the bacterial composition of the samples. We also evaluated the flat refraction curve, and the results suggested that the quantity of sequencing data was sufficient (Figure S1B). We examined estimators of community richness (Sobs index), diversity, and evenness (Shannon index). Significant differences were observed in the Sobs and Shannon indexes between the ESCC and HC groups (Sobs, 324.59 vs. 565.69, P = 0.008; Shannon, 2.95 vs. 3.64, P = 0.017), the EGJ and HC groups (Sobs, 331.36 vs. 565.69, P = 0.044; Shannon, 2.48 vs. 3.64, P = 0.039), and the post-ESCC and HC groups (Sobs, 318.47 vs. 565.69, P = 0.021; Shannon, 2.48 vs. 3.64, P = 0.014, Figures 1A,B), indicating that microbial diversity was significantly lower in the ESCC, EGJ, and post-ESCC groups than that in the HC group. In addition, no significant differences were observed in Sobs and Shannon indexes between ESCC and EGJ groups, whose α diversity was not significantly different (Figure S1C).

FIGURE 1

Figure 1. Comparison of α diversity and the relative abundance of taxa in the healthy patients and those of the EGJ, ESCC, and post-ESCC groups. Estimators of community richness [Sobs index, (A)] and diversity [Shannon index, (B)] in OTU levels. (C) Average relative abundance of taxa at the phylum level. (D) Average relative abundance of taxa at the genus level. OTU, operational taxonomic unit.

Microbial Differences Between Tumor and Healthy Tissues

We analyzed the microbiomes of the ESCC, EGJ, and HC groups. Firmicutes was the predominant phylum in ESCC and EGJ tumor specimens, whereas Proteobacteria was the predominant phylum in healthy specimens (Figure 1C). We also examined the microbial community at the genus level. ESCC specimens were rich in Streptococcus, Lactobacillus, Prevotella, and Fusobacterium, while EGJ specimens were abundant in Streptococcus, Staphylococcus, and Pseudomonas. By contrast, HC specimens were rich in Streptococcus, Ralstonia, and Burkholderia–Caballeronia–Paraburkholderia (Figure 1D). We carried out tests of significance between tumor and healthy specimens. The relative abundance of Fusobacteria was higher (7.01 vs. 1.12%, P = 0.039) and that of Actinobacteria was lower (1.61 vs. 4.04%, P = 0.002) in the ESCC group than in the HC group (Figure 2A), while the EGJ group showed an increased abundance of Firmicutes (61.24 vs. 26.66%, P = 0.001) but a decreased abundance of Proteobacteria compared with the HC group (16.59 vs. 39.15%, P = 0.002, Figure 2B). A difference in the relative abundance of only one phylum, Firmicutes, was observed between the ESCC and EGJ groups (38.38 vs. 61.24%, P = 0.025, Figure S2A). Interestingly, we identified an unclassified bacterium whose relative abundance was higher in the HC group than in the ESCC and EGJ groups (13.69 vs. 1.50%, P = 0.005; 13.69 vs. 0.82%, P = 0.003, Figures 2A,B). Compared with the HC group, the relative proportion of Pseudomonas increased (5.33 vs. 1.46%, P = 0.003) and those of Ralstonia and Burkholderia–Caballeronia–Paraburkholderia decreased (0.21 vs. 11.38%, P < 0.001; 0.77 vs. 9.28%, P < 0.001, Table S1) in the ESCC group.

FIGURE 2

Figure 2. Composition of the esophageal microbiomes of healthy patients and those with ESCC, and β diversity for the EGJ, ESCC, and HC groups. (A) Significant differences were observed between the microbiomes of the ESCC and HC groups. (B) Significant differences were observed between the microbiomes of the EGJ and HC groups. (C) PCoA based on unweighted UniFrac distances between the ESCC and HC groups. P-values were calculated by analysis of similarities (ANOSIM) *, 0.01 < p < 0.05; **, p < 0.01.

To further evaluate differences in the esophageal microbiome, we examined the β diversity and represented the results of principal coordinate analysis (PCoA) as a plot. The results revealed significant differences between the esophageal microbiomes of the ESCC and HC groups (Figure 2C, Figure S2B). We also examined the LDA coupled with LEfSe. We found that Clostridiales, Pseudomonas, and Selenomonadales were the key taxa contributing to the changes in the microbiome of ESCC patients, whereas Burkholderiaceae, Ralstonia, and Burkholderia–Caballeronia–Paraburkholderia were the key taxa contributing to the changes in the microbiome of HC patients (Figures 3A,B). The results of LDA coupled with LEfSe between the EGJ and HC groups are also presented (Figures S2C, S3A).

FIGURE 3

Figure 3. Results of linear discriminant analysis (LDA) and effect size measurements (LEfSe) between the ESCC and HC groups. (A) Bar plot shows taxa with LDA score >3.0 from the order to the genus level. (B) LEfSe analysis shows the most abundant taxa from the phylum to the genus level between the ESCC and HC groups.

Differences Between the Microbial Communities of ESCC and Post-ESCC Patients

To examine microbial differences caused by esophagectomy, we examined the α and β diversities of the ESCC and post-ESCC groups. There are no significant differences in their Sobs and Shannon indexes (Figure S3B), indicating that the community richness, diversity, and evenness were similar between the two groups. However, the results of PCoA demonstrated that the microbial composition was significantly different (Figure 4A). We found significant differences in the proportions of Bacteroidetes (20.45 vs. 9.86%, P = 0.026) and Fusobacteria (7.01 vs. 1.66%, P = 0.030) at the phylum level (Figure 4B) and in the proportions of Pseudomonas (5.33 vs. 13.83%, P = 0.008), Fusobacterium (5.77 vs. 1.25%, P = 0.045), and Prevotella (4.42 vs. 1.02%, P = 0.001) at the genus level between the ESCC and post-ESCC groups (Table S2). The β diversity results suggested that the microbial community of the post-ESCC group was similar to that of the HC group (Figure S4). The results of LDA and LEfSe confirmed that Bacteroidetes and Pseudomonas were the key taxa contributing to the changes in the microbiome of the ESCC and post-ESCC groups (Figures 5A,B).

FIGURE 4

Figure 4. Composition of the esophageal microbiomes and β diversity of the ESCC and post-ESCC groups. (A) PCoA based on weighted UniFrac distances between the ESCC and post-ESCC groups. P-values were calculated by analysis of similarities (ANOSIM). (B) Significant differences were observed between the microbiomes of the ESCC and post-ESCC groups. *P < 0.05.

FIGURE 5

Figure 5. Results of linear discriminant analysis (LDA) and effect size measurements (LEfSe) between the ESCC and HC groups from the phylum to the genus level. (A) Bar plot shows results of taxa with LDA score >3.0. (B) LEfSe analysis shows the most differentially abundant taxa between the ESCC and post-ESCC groups.

Differences in Microbial Communities Lead to Differences Between the Metabolic Pathways in ESCC and HC Patients

Microbial imbalances can induce systemic metabolic alterations (Nieuwdorp et al., 2014) and vice versa (Cani, 2017). We carried out PICRUSt to predict the metagenomes from the 16S data, and these were used to identify the KEGG pathways involved in ESCC, HC, and post-ESCC specimens. The results are presented as a heatmap. We found that the HC group was enriched in pathways related to the metabolism of fatty acids, short-chain fatty acids (SCFAs), butanoate, propanoate, tryptophan, and beta-alanine, as well as in pathways related to the degradation of benzoate, lysine, geraniol, aminobenzoate, limonene, and pinene. The ESCC group was enriched in pathways related to metabolism of cysteine, methionine, fructose, galactose, and starch, as well as in pathways related to DNA repair and recombination, protein translation, chromosomal dynamics, and peptidase activity (Figure 6). The PICRUSt taxonomic functional relationships suggested that the microbial composition determined which metabolic pathways were involved. Interestingly, the pathways associated with the post-ESCC group were moderately similar to those of the HC group (Figure 6).

FIGURE 6

Figure 6. Hierarchical clustering analysis to identify different metabolic pathways. Heatmap of PICRUSt analysis demonstrating the enrichment or loss of metabolic pathways in the ESCC, HC, and post-ESCC groups.

Discussion

Microbes colonize the esophageal mucosa, and increasing evidence suggests that the esophageal microbiome plays important roles in esophageal diseases. Most studies on the esophageal microbiome and esophageal diseases address BE and EAC, but there are few studies concerning ESCC. This is the first study to profile the esophageal microbiome in cancer patients before and after esophagectomy and to predict systemic metabolic alterations. We found the microbial richness and diversity (Sobs and Shannon indexes) to be significantly lower in esophageal cancer and postoperative patients than in healthy controls, which was consistent with the results of a previous study (Castaño-Rodríguez et al., 2017). The α diversity and β diversity were not statistically different between the ESCC and EGJ groups, which is possibly because that the tumor microenvironment creates a niche for similar microbial communities to colonize and thrive even in unique tissue types (Shao et al., 2019). Proteobacteria and Streptococcus were the predominant taxa in normal esophagus, which was in accordance with earlier studies. In addition, we detected Ralstonia and Burkholderia–Caballeronia–Paraburkholderia, which were members of the Betaproteobacteria family, in healthy controls. Ralstonia was the most abundant genus in normal breast tissue and gastric cancer (Tseng et al., 2016; Costantini et al., 2018). Our results showed that Ralstonia comprised the microbiome of the normal esophagus, and further studies are needed to verify these results.

The esophagus of patients with ESCC was enriched in Streptococcus, Prevotella, Fusobacterium, Veillonella, and Lactobacillus, which is in conformity with a previous study (Shao et al., 2019). In Japanese patients with esophageal cancer, Fusobacterium nucleatum was enriched in ESCC tissue compared with matched normal esophagus tissue, and is quantity was negatively correlated with survival (Yamamura et al., 2016, 2019). F. nucleatum has potential prognostic value for ESCC and had been detected in colon cancer. Other studies have reported a robust association between F. nucleatum and colorectal cancer (Castellarin et al., 2012; Kostic et al., 2013; Rubinstein et al., 2013; Yu et al., 2017). The F. nucleatum presented in these tissues had been verified by isolating F. nucleatum strains directly from biopsy samples and from patient-derived xenografts passaged in mice (Castellarin et al., 2012; Bullman et al., 2017). F. nucleatum can induce expression of the pro-inflammatory cytokines in epithelial cells, including IL-6 and IL-8 (Ahn et al., 2017), which can contribute to the dynamic cross-talk between tumor cells and cancer-associated fibroblasts (CAF) in the TME for ESCC (Karakasheva et al., 2018). F. nucleatum infection can induce expression of the antimicrobial peptide β-defensin 2 and high-mobility group box 1 protein (HMGB1) (Bui et al., 2016), which can result in increased proliferation and modulated autophagy in ESCC cell lines (Di et al., 2019). In the ESCC group, the relative abundances of Clostridiales and Pseudomonas also significantly increased compared to the HC group. Clostridiales are obligate anaerobes that thrive in hypoxic environments such as those created by tumors. A previous study reported that Clostridiales in the gastric corpus could contribute to esophageal squamous dysplasia (Nasrollahzadeh et al., 2015). Transplanting Clostridium to germ-free nutrition-deficient mice can increase the acylcarnitine level in the gut and decrease protein synthesis and amino acid oxidation in the liver (Blanton et al., 2016), thereby promoting the development of cancer. Clostridiumalos difficile, as a member of the Clostridiales class, also plays a pivotal role in regulating Clostridium difficile infections (CDI) (Farowski et al., 2019). The Gram-negative Pseudomonas is over-represented in cases of oral squamous cell carcinoma (OSCC) (Perera et al., 2018), a kind of tumor that is similar to ESCC. Pseudomonas can induce Toll-like receptors to activate cytokines, chemokines, and COX-2 and recruit cells of the innate and adaptive immune system (Markou and Apidianakis, 2014). Infection with Pseudomonas could activate the c-Jun N-terminal kinase (JNK) pathway to induce enterocyte apoptosis, and intestinal stem cell proliferation (Apidianakis et al., 2009), intestinal innate immune responses, and stem cells may initiate the development of tumors and metastasis (Schwitalla et al., 2013). In addition, no significant change was observed in the relative abundance of Porphyromonas gingivalis in our study, which was evidenced in esophageal cancer and dysplasia tissues but was rarely found in non-cancerous and normal tissues (Gao et al., 2016).

We also noticed differences in the esophageal microbiome after surgery. At the genus level, the abundance of Fusobacterium and Prevotella decreased after surgery. In another study, the relative proportions of Firmicutes and Bacteroidetes increased, while those of Proteobacteria and Actinobacteria decreased after gastrectomy (Tseng et al., 2016). Bacteroidetes are Gram-negative anaerobes and micro-aerophiles that reside in the diseased esophagus (Yang et al., 2009). Compared to the HC group, Lactobacillus significantly increased in the post-ESCC group, which was associated with the acidic microenvironment caused by the increased gastric acid reflux after esophagectomy. Another study also indicated that a low pH is needed for the proliferation of Lactobacillus (Elliott et al., 2017). The gastric microenvironment is populated by microbial communities mainly comprised of the Lactobacillus, Streptococcus, and Propionibacterium genera (Nardone et al., 2017).

The results of PICRUSt showed that the pathways related to the metabolism of cysteine and methionine, monosaccharides, starch were upregulated and the metabolic levels of fatty acids, short-chain fatty acids (SCFAs), tryptophan, and beta-alanine were low in the ESCC group. Meanwhile, pathways related to other cellular functions, including DNA repair and recombination, protein translation, chromosomal dynamics, and peptidase activity were also altered. Methionine uptake and breakdown associate with various cellular processes such as methylation reactions, polyamine synthesis, and redox maintenance. SCFAs are produced by colonic microbial fermentation of undigested or partially digested dietary fibers and have broad effects on host immune system development and function, which can enhance epithelial barrier function and immune tolerance (Rooks and Garrett, 2016). Tryptophan has been shown to play a crucial role in the balance between intestinal immune tolerance and gut microbiota maintenance (Lee and Lee, 2010). The metabolic heterogeneity between the ESCC and HC groups was essential for host immune function and tumorigenesis. After radical esophagectomy, the microbial metabolism became similar to that of healthy control, suggesting that tumors can alter the microbial microenvironment and metabolic profile. The microbial composition of the esophagus may tend to become similar to normal status after surgery. It is suggested that the microbiota composition and metabolism profile are associated with the physiological or pathological status of the esophagus and that the emergence of some microbes, which dominated in the tumor patients, may be predictors of the recurrence of the tumor.

Of course, there were several limitations to this study. First, the sample size of each group was limited, and larger studies should be performed to validate our findings. Second, the microbiome that we studied was likely to be contaminated by saliva during endoscopic sampling. Finally, we only examined the microbiome at 9–12 months after surgery. Additional time points may better describe the changes in the microbial microenvironment for postoperative cases. In addition, investigations of the biological mechanisms of the microbiota are needed to elucidate the association between complex microbial environments and esophageal cancer.

In conclusion, we compared the microbial composition of esophagus in different statuses and found some microbiota associated with the tumorigenesis. We also discovered the microbial metabolism of postoperative esophagus to be similar to that of healthy tissue. Monitoring of the esophagus microbiota may be an essential method for predicting the recurrence of tumor.

Data Availability Statement

All sequencing data associated with this study have been uploaded to the NCBI SRA database (SRA accession: PRJNA628659).

Ethics Statement

The studies involving human participants were reviewed and approved by the Institutional Review Board of the Renmin Hospital of Wuhan University. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

DL contributed to study concept and design, data collection, analysis, and interpretation, and review of the manuscript. RH was responsible for the data collection, analysis, and interpretation and drafting of the manuscript. TF worked on the interpretation and review of the manuscript. GH helped with data collection, interpretation, and review of the manuscript. WM, LC, LZ, WJ, WW, ZL, and HF was involved in data collection and review of the manuscript. QG contributed to study concept and design, data collection, analysis, and interpretation, and critical review of the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This study was supported by the National Natural Science Foundation of China (Grant No. 81770095).

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.

Acknowledgments

We are grateful to TopEdit (www.topeditsci.com) for English language editing of this manuscript and endoscopists at Renmin Hospital, Wuhan University for specimen collection. We also thank pathologists at Renmin Hospital, Wuhan University for pathological diagnosis.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2020.00268/full#supplementary-material

References

Ahn, S.-H., Chun, S.-M., Park, C., Lee, J.-H., Lee, S.-W., and Lee, T.-H. (2017). Transcriptome profiling analysis of senescent gingival fibroblasts in response to Fusobacterium nucleatum infection. PLoS ONE 12:e0188755. doi: 10.1371/journal.pone.0188755

PubMed Abstract | CrossRef Full Text | Google Scholar

Apidianakis, Y., Pitsouli, C., Perrimon, N., and Rahme, L. (2009). Synergy between bacterial infection and genetic predisposition in intestinal dysplasia. Proc. Natl. Acad. Sci. U.S.A. 106, 20883–20888. doi: 10.1073/pnas.0911797106

PubMed Abstract | CrossRef Full Text | Google Scholar

Arnold, M., Soerjomataram, I., Ferlay, J., and Forman, D. (2015). Global incidence of oesophageal cancer by histological subtype in 2012. Gut 64, 381–387. doi: 10.1136/gutjnl-2014-308124

PubMed Abstract | CrossRef Full Text | Google Scholar

Bao, K., Belibasakis, G. N., Thurnheer, T., Aduse-Opoku, J., Curtis, M. A., and Bostanci, N. (2014). Role of Porphyromonas gingivalis gingipains in multi-species biofilm formation. BMC Microbiol. 14:258. doi: 10.1186/s12866-014-0258-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Blanton, L. V., Charbonneau, M. R., Salih, T., Barratt, M. J., Venkatesh, S., Ilkaveya, O., et al. (2016). Gut bacteria that prevent growth impairments transmitted by microbiota from malnourished children. Science 351:aad3311. doi: 10.1126/science.aad3311

PubMed Abstract | CrossRef Full Text | Google Scholar

Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. doi: 10.1093/bioinformatics/btu170

PubMed Abstract | CrossRef Full Text | Google Scholar

Brawner, K. M., Morrow, C. D., and Smith, P. D. (2014). Gastric microbiome and gastric cancer. Cancer J. 20, 211–216. doi: 10.1097/PPO.0000000000000043

PubMed Abstract | CrossRef Full Text | Google Scholar

Bray, F., Ferlay, J., Soerjomataram, I., Siegel, R. L., Torre, L. A., and Jemal, A. (2018). Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 68, 394–424. doi: 10.3322/caac.21492

PubMed Abstract | CrossRef Full Text | Google Scholar

Bui, F. Q., Johnson, L., Roberts, J., Hung, S.-C., Lee, J., Atanasova, K. R., et al. (2016). Fusobacterium nucleatum infection of gingival epithelial cells leads to NLRP3 inflammasome-dependent secretion of IL-1β and the danger signals ASC and HMGB1: inflammasome activation by Fusobacterium nucleatum. Cell. Microbiol. 18, 970–981. doi: 10.1111/cmi.12560

PubMed Abstract | CrossRef Full Text | Google Scholar

Bullman, S., Pedamallu, C. S., Sicinska, E., Clancy, T. E., Zhang, X., Cai, D., et al. (2017). Analysis of Fusobacterium persistence and antibiotic response in colorectal cancer. Science 358, 1443–1448. doi: 10.1126/science.aal5240

PubMed Abstract | CrossRef Full Text | Google Scholar

Cani, P. D. (2017). Gut cell metabolism shapes the microbiome. Science 357, 548–549. doi: 10.1126/science.aao2202

PubMed Abstract | CrossRef Full Text | Google Scholar

Castaño-Rodríguez, N., Goh, K.-L., Fock, K. M., Mitchell, H. M., and Kaakoush, N. O. (2017). Dysbiosis of the microbiome in gastric carcinogenesis. Sci. Rep. 7:15957. doi: 10.1038/s41598-017-16289-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Castellarin, M., Warren, R. L., Freeman, J. D., Dreolini, L., Krzywinski, M., Strauss, J., et al. (2012). Fusobacterium nucleatum infection is prevalent in human colorectal carcinoma. Genome Res. 22, 299–306. doi: 10.1101/gr.126516.111

PubMed Abstract | CrossRef Full Text | Google Scholar

Costantini, L., Magno, S., Albanese, D., Donati, C., Molinari, R., Filippone, A., et al. (2018). Characterization of human breast tissue microbiota from core needle biopsies through the analysis of multi hypervariable 16S-rRNA gene regions. Sci. Rep. 8:16893. doi: 10.1038/s41598-018-35329-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Di Pilato, V., Freschi, G., Ringressi, M. N., Pallecchi, L., Rossolini, G. M., and Bechi, P. (2016). The esophageal microbiota in health and disease: Esophageal microbiota and disease. Ann. N. Y. Acad. Sci. 1381, 21–33. doi: 10.1111/nyas.13127

PubMed Abstract | CrossRef Full Text | Google Scholar

Di, X., He, G., Chen, H., Zhu, C., Qin, Q., Yan, J., et al. (2019). High-mobility group box 1 protein modulated proliferation and radioresistance in esophageal squamous cell carcinoma. J. Gastroenterol. Hepatol. 34, 728–735. doi: 10.1111/jgh.14371

PubMed Abstract | CrossRef Full Text | Google Scholar

Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi: 10.1038/nmeth.2604

PubMed Abstract | CrossRef Full Text | Google Scholar

Elliott, D. R. F., Walker, A. W., O'Donovan, M., Parkhill, J., and Fitzgerald, R. C. (2017). A non-endoscopic device to sample the oesophageal microbiota: a case-control study. Lancet Gastroenterol. Hepatol. 2, 32–42. doi: 10.1016/S2468-1253(16)30086-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Farowski, F., Els, G., Tsakmaklis, A., Higgins, P. G., Kahlert, C. R., Stein-Thoeringer, C. K., et al. (2019). Assessment of urinary 3-indoxyl sulfate as a marker for gut microbiota diversity and abundance of Clostridiales. Gut Microbes 10, 133–141. doi: 10.1080/19490976.2018.1502536

PubMed Abstract | CrossRef Full Text | Google Scholar

Ferlay, J., Colombet, M., Soerjomataram, I., Mathers, C., Parkin, D. M., Piñeros, M., et al. (2018). Estimating the global cancer incidence and mortality in 2018: GLOBOCAN sources and methods. Int. J. Cancer 144, 1941–1953. doi: 10.1002/ijc.31937

PubMed Abstract | CrossRef Full Text | Google Scholar

Fitzmaurice, C., Dicker, D., Pain, A., Hamavid, H., Moradi-Lakeh, M., MacIntyre, M. F., et al. (2015). The global burden of cancer 2013. JAMA Oncol. 1:505. doi: 10.1001/jamaoncol.2015.0735

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, S., Li, S., Ma, Z., Liang, S., Shan, T., Zhang, M., et al. (2016). Presence of Porphyromonas gingivalis in esophagus and its association with the clinicopathological characteristics and survival in patients with esophageal cancer. Infect. Agents Cancer 11:3. doi: 10.1186/s13027-016-0049-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kageyama, S., Takeshita, T., Takeuchi, K., Asakawa, M., Matsumi, R., Furuta, M., et al. (2019). Characteristics of the salivary microbiota in patients with various digestive tract cancers. Front. Microbiol. 10:1780. doi: 10.3389/fmicb.2019.01780

PubMed Abstract | CrossRef Full Text | Google Scholar

Karakasheva, T. A., Lin, E. W., Tang, Q., Qiao, E., Waldron, T. J., Soni, M., et al. (2018). IL-6 Mediates cross-talk between tumor cells and activated fibroblasts in the tumor microenvironment. Cancer Res. 78, 4957–4970. doi: 10.1158/0008-5472.CAN-17-2268

PubMed Abstract | CrossRef Full Text | Google Scholar

Kostic, A. D., Chun, E., Robertson, L., Glickman, J. N., Gallini, C. A., Michaud, M., et al. (2013). Fusobacterium nucleatum potentiates intestinal tumorigenesis and modulates the tumor-immune microenvironment. Cell Host Microbe 14, 207–215. doi: 10.1016/j.chom.2013.07.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Lagergren, J., Smyth, E., Cunningham, D., and Lagergren, P. (2017). Oesophageal cancer. Lancet 390, 2383–2396. doi: 10.1016/S0140-6736(17)31462-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Langille, M. G. I., Zaneveld, J., Caporaso, J. G., McDonald, D., Knights, D., Reyes, J. A., et al. (2013). Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat. Biotechnol. 31, 814–821. doi: 10.1038/nbt.2676

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, J.-H., and Lee, J. (2010). Indole as an intercellular signal in microbial communities. FEMS Microbiol. Rev. 34, 426–444. doi: 10.1111/j.1574-6976.2009.00204.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, H., Li, J., Zhao, Y., Xie, K., Tang, X., Wang, S., et al. (2018). Ammonia oxidizers and nitrite-oxidizing bacteria respond differently to long-term manure application in four paddy soils of south of China. Sci. Total Environ. 633, 641–648. doi: 10.1016/j.scitotenv.2018.03.108

PubMed Abstract | CrossRef Full Text | Google Scholar

Magoc, T., and Salzberg, S. L. (2011). FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963. doi: 10.1093/bioinformatics/btr507

PubMed Abstract | CrossRef Full Text

Markou, P., and Apidianakis, Y. (2014). Pathogenesis of intestinal Pseudomonas aeruginosa infection in patients with cancer. Front. Cell. Infect. Microbiol. 3:115. doi: 10.3389/fcimb.2013.00115

PubMed Abstract | CrossRef Full Text | Google Scholar

McDonald, D., Price, M. N., Goodrich, J., Nawrocki, E. P., DeSantis, T. Z., Probst, A., et al. (2012). An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 6, 610–618. doi: 10.1038/ismej.2011.139

PubMed Abstract | CrossRef Full Text | Google Scholar

Nardone, G., Compare, D., and Rocco, A. (2017). A microbiota-centric view of diseases of the upper gastrointestinal tract. Lancet Gastroenterol. Hepatol. 2, 298–312. doi: 10.1016/S2468-1253(16)30108-X

PubMed Abstract | CrossRef Full Text | Google Scholar

Nasrollahzadeh, D., Malekzadeh, R., Ploner, A., Shakeri, R., Sotoudeh, M., Fahimi, S., et al. (2015). Variations of gastric corpus microbiota are associated with early esophageal squamous cell carcinoma and squamous dysplasia. Sci. Rep. 5:8820. doi: 10.1038/srep08820

PubMed Abstract | CrossRef Full Text | Google Scholar

Nieuwdorp, M., Gilijamse, P. W., Pai, N., and Kaplan, L. M. (2014). Role of the microbiome in energy regulation and metabolism. Gastroenterology 146, 1525–1533. doi: 10.1053/j.gastro.2014.02.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Perera, M., Al-hebshi, N. N., Perera, I., Ipe, D., Ulett, G. C., Speicher, D. J., et al. (2018). Inflammatory bacteriome and oral squamous cell carcinoma. J. Dent. Res. 97, 725–732. doi: 10.1177/0022034518767118

PubMed Abstract | CrossRef Full Text | Google Scholar

Peters, B. A., Wu, J., Pei, Z., Yang, L., Purdue, M. P., Freedman, N. D., et al. (2017). Oral microbiome composition reflects prospective risk for esophageal cancers. Cancer Res. 77, 6777–6787. doi: 10.1158/0008-5472.CAN-17-1296

PubMed Abstract | CrossRef Full Text | Google Scholar

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.e12. doi: 10.1016/j.cell.2019.07.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Rooks, M. G., and Garrett, W. S. (2016). Gut microbiota, metabolites and host immunity. Nat. Rev. Immunol. 16, 341–352. doi: 10.1038/nri.2016.42

PubMed Abstract | CrossRef Full Text | Google Scholar

Rubinstein, M. R., Wang, X., Liu, W., Hao, Y., Cai, G., and Han, Y. W. (2013). Fusobacterium nucleatum promotes colorectal carcinogenesis by modulating E-Cadherin/β-Catenin signaling via its FadA adhesin. Cell Host Microbe 14, 195–206. doi: 10.1016/j.chom.2013.07.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Schirmer, M., Ijaz, U. Z., D'Amore, R., Hall, N., Sloan, W. T., and Quince, C. (2015). Insight into biases and sequencing errors for amplicon sequencing with the Illumina MiSeq platform. Nucleic Acids Res. 43:e37. doi: 10.1093/nar/gku1341

PubMed Abstract | CrossRef Full Text | Google Scholar

Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/AEM.01541-09

PubMed Abstract | CrossRef Full Text | Google Scholar

Schwitalla, S., Fingerle, A. A., Cammareri, P., Nebelsiek, T., Göktuna, S. I., Ziegler, P. K., et al. (2013). Intestinal tumorigenesis initiated by dedifferentiation and acquisition of stem-cell-like properties. Cell 152, 25–38. doi: 10.1016/j.cell.2012.12.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol. 12:R60. doi: 10.1186/gb-2011-12-6-r60

PubMed Abstract | CrossRef Full Text

Shao, D., Vogtmann, E., Liu, A., Qin, J., Chen, W., Abnet, C. C., et al. (2019). Microbial characterization of esophageal squamous cell carcinoma and gastric cardia adenocarcinoma from a high-risk region of China. Cancer. 125:3993–4002. doi: 10.1002/cncr.32403

PubMed Abstract | CrossRef Full Text | Google Scholar

Tseng, C.-H., Lin, J.-T., Ho, H. J., Lai, Z.-L., Wang, C.-B., Tang, S.-L., et al. (2016). Gastric microbiota and predicted gene functions are altered after subtotal gastrectomy in patients with gastric cancer. Sci. Rep. 6:20701. doi: 10.1038/srep20701

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, H., Rusiecki, J. A., Zhu, K., Potter, J., and Devesa, S. S. (2009). Stomach carcinoma incidence patterns in the United States by histologic type and anatomic site. Cancer Epidemiol. Biomarkers. Prev. 18, 1945–1952. doi: 10.1158/1055-9965.EPI-09-0250

PubMed Abstract | CrossRef Full Text | Google Scholar

Yamamura, K., Baba, Y., Nakagawa, S., Mima, K., Miyake, K., Nakamura, K., et al. (2016). Human Microbiome Fusobacterium Nucleatum in esophageal cancer tissue is associated with prognosis. Clin. Cancer Res. 22, 5574–5581. doi: 10.1158/1078-0432.CCR-16-1786

PubMed Abstract | CrossRef Full Text | Google Scholar

Yamamura, K., Izumi, D., Kandimalla, R., Sonohara, F., Baba, Y., Yoshida, N., et al. (2019). Intratumoral Fusobacterium nucleatum levels predict therapeutic response to neoadjuvant chemotherapy in esophageal squamous cell carcinoma. Clin. Cancer Res. 25, 6170–6179. doi: 10.1158/1078-0432.CCR-19-0318

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J., Mu, X., Wang, Y., Zhu, D., Zhang, J., Liang, C., et al. (2018). Dysbiosis of the salivary microbiome is associated with non-smoking female lung cancer and correlated with immunocytochemistry markers. Front. Oncol. 8:520. doi: 10.3389/fonc.2018.00520

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, L., Lu, X., Nossa, C. W., Francois, F., Peek, R. M., and Pei, Z. (2009). Inflammation and intestinal metaplasia of the distal esophagus are associated with alterations in the microbiome. Gastroenterology 137, 588–597. doi: 10.1053/j.gastro.2009.04.046

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Y., Weng, W., Peng, J., Hong, L., Yang, L., Toiyama, Y., et al. (2017). Fusobacterium nucleatum increases proliferation of colorectal cancer cells and tumor development in mice by activating toll-like receptor 4 signaling to nuclear factor–κB, and up-regulating expression of microRNA-21. Gastroenterology 152, 851–866.e24. doi: 10.1053/j.gastro.2016.11.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, G., Gail, M. H., Shi, J., Klepac-Ceraj, V., Paster, B. J., Dye, B. A., et al. (2014). Association between upper digestive tract microbiota and cancer-predisposing states in the esophagus and stomach. Cancer Epidemiol. Biomarkers Prevention 23, 735–741. doi: 10.1158/1055-9965.EPI-13-0855

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, T., Guo, F., Yu, Y., Sun, T., Ma, D., Han, J., et al. (2017). Fusobacterium nucleatum promotes chemoresistance to colorectal cancer by modulating autophagy. Cell 170, 548–563.e16. doi: 10.1016/j.cell.2017.07.008

PubMed Abstract | CrossRef Full Text | Google Scholar