Abstract
Cancer immunotherapy including immune checkpoint inhibitors (ICI) has revolutionized non-small cell lung cancer (NSCLC) therapy. Recently, the microbiota status “before” initiation of ICI therapy has been emphasized as a predictive biomarker in patients undergoing ICI therapy. However, the microbiota diversity and composition “during” ICI therapy is unknown. This multicenter, prospective observational study analyzed both saliva and feces from 28 patients with NSCLC. We performed 16S ribosomal RNA gene sequencing, then analyzed associations of oral and gut microbiota diversity or composition with ICI response. At the genus level, the alpha diversity of the gut microbiota was significantly greater in responders (n = 17) than in non-responders (n = 11) (Chao 1, p = 0.0174; PD whole tree, p = 0.0219; observed species, p = 0.0238; Shannon, p = 0.0362), while the beta diversity of the gut microbiota was significantly different (principal coordinates analysis, p = 0.035). Compositional differences in the gut microbiota were observed between the two groups; in particular, g_Blautia was enriched in responders, whereas o_RF32 order unclassified was enriched in non-responders. The progression-free survival (PFS) of patients enriched gut microbiota of g_Blautia was significantly longer [median survival time (MST): not reached vs. 549 days, p = 0.0480] and the PFS of patients with gut microbiota of o_RF32 unclassified was significantly shorter (MST: 49 vs. 757 days, p = 0.0205). There were no significant differences between groups in the oral microbiota. This study revealed a strong association between gut microbiota diversity and ICI response in NSCLC patients. Moreover, specific gut microbiota compositions may influence the ICI response. These findings might be useful in identifying biomarkers to predict ICI response.
Introduction
Immunotherapy with immune checkpoint inhibitors (ICI) is widely used to treat various malignancies, including non-small cell lung cancer (NSCLC); it has revolutionized therapeutic approaches to cancer. Programmed death-ligand-1 (PD-L1) is an immune checkpoint protein expressed on tumor cells and tumor-infiltrating immune cells, which can mediate anticancer immunosuppression (Ribas and Wolchok, 2018). Anti-PD-1 antibodies (e.g., nivolumab and pembrolizumab) and anti-PD-L1 antibodies (e.g., atezolizumab and durvalumab) enable T-cell activation and immune system recognition.
Although tumorous PD-L1 expression is a potential biomarker of the ICI therapeutic response, there is no widely accepted optimal biomarker to predict the efficacy of ICI, because ICI response and survival outcomes show heterogeneity in NSCLC patients receiving ICI therapy, regardless of PD-L1 expression level (Rittmeyer et al., 2017; Gandhi et al., 2018; Paz-Ares et al., 2018; Socinski et al., 2018).
We recently reported that the pretreatment host immune-nutritional condition was a prognostic marker for NSCLC patients receiving ICI therapy (Shoji et al., 2019). Host immunity is clearly associated with the ICI response. The internal microbiome is regarded as a controlling factor in host immunity. In particular, the gut microbiome can modulate the host immune response (e.g., anti-tumor immunity) and optimize both innate and adaptive immune responses (Littman, 2012). Recently, preclinical studies have shown that the gut microbiome composition and its modification in murine models could influence the efficacy of ICI (Corrales et al., 2015; Vétizou M, 2015). Therefore, the microbiome has been emphasized as a predictive biomarker of ICI therapy, mainly in studies from the United States or Europe. Additionally, the gut microbiota diversity or abundance of specific gut microbiome components has been correlated with the efficacy of anti-PD-1 antibody in melanoma patients (Fessler et al., 2018). Moreover, fecal microbiota transplantation (FMT) in murine models might restore the ICI response (Gopalakrishnan et al., 2018; Routy et al., 2018). In a recent study, FMT from ICI responders to ICI non-responders produced ICI efficacy in melanoma patients (Davar et al., 2021).
Furthermore, the oral microbiota has been associated with several diseases (e.g., inflammatory bowel disease and allergic diseases) through its influence on the gut microbiota (Pietrantoni et al., 2003; Yamasaki et al., 2013; Atarashi K, 2017). A recent study revealed that variation in the oral microbiota was associated with a risk of lung cancer (Hosgood et al., 2021). However, samples were collected prior to ICI therapy in most previous studies, and thus minimal information has been available regarding the microbiota status during ICI therapy. Accordingly, FMT or biotics therapy approaches are needed to investigate changes in the microbiota during ICI therapy. Notably, there are definite differences in microbiota composition among ethnicities (Nishijima et al., 2016); to the best of our knowledge, few reports have been published regarding Japanese NSCLC patients (Hakozaki et al., 2020). In addition, the present study is based on the clinical question how is the condition of host microbiome in NSCLC patients during ICI therapy. Therefore, the present study might be meaningful as one of the pioneer studies highlighted the microbiome status during ICI therapy.
Here, we performed a prospective study to clarify the microbiota diversity and composition in Japanese NSCLC patients by analyzing samples collected during ICI therapy.
Materials and methods
Study design and participants
This prospective observational study was conducted at multiple centers: Department of Thoracic Surgery and Department of Respiratory Medicine, Clinical Research Institute, National Hospital Organization (NHO) Kyushu Medical Center and Department of Thoracic Oncology, NHO Kyushu Cancer Center. Eligibility criteria were as follows: pathologically or cytologically confirmed diagnosis of locally advanced/unresectable or postoperative recurrent NSCLC; receipt of ICI monotherapy including nivolumab (Opdivo, Bristol-Myers Squibb), pembrolizumab (Keytruda, Merck), atezolizumab (Tecentriq, Genentech), and durvalumab (Imfinzi, Astra Zeneca) or platinum-based therapy combined with these ICI agents. Patients were also enrolled if they had discontinued these therapies but had not received any additional therapies. Patients with the presence of ongoing antibiotics therapy for infectious diseases before/during ICI therapy, were excluded in the present study.
From July 2019 to December 2020, 34 NSCLC patients were eligible and enrolled. Of those 34 patients, 28 had both saliva and feces samples available for this study (Supplementary Figure S1). All enrolled patients had at least one measurable target lesion based on the Response Evaluation Criteria in Solid Tumors (RECIST), version 1.1 (Eisenhauer et al., 2009). Clinical/pathological stage was based on the Tumor Node Metastasis (TNM) classification established by the International Union Against Cancer (Goldstraw et al., 2016). For TNM staging, all patients underwent computed tomography (CT) of the thorax and upper abdomen, as well as bone scintigrams, brain CT scans, magnetic resonance imaging (MRI), or fluorodeoxyglucose-positron emission tomography (FDG-PET). Postoperative local or distant recurrence was defined as described previously (Varlotto et al., 2009). ICI therapy was continued until radiographic progression. PD-L1 protein expression was evaluated using antibody clone 22C3 (Dako, Agilent Technologies, Santa Clara, CA, United States).
Sample collection
Salivary and fecal samples were collected in sterile containers and immediately placed at 4°C, then frozen at −80°C. Individual periods of sample collection are shown in Figure 1. The mean numbers of days between the initiation of ICI therapy and the day of sample collection were 307 (29–945) days in ICI responders and 117 (23–491) days in ICI non-responders.
FIGURE 1
DNA extraction, gene amplification, sequencing, and data analysis procedures
Preliminary treatment of fecal samples was conducted in accordance with a previously described method (Takahashi et al., 2014); DNA was then extracted using an automated DNA isolation system (Gene Prep Star PI-480, Kurabo, Japan). DNA was extracted from saliva using the Mora-Extract kit (Kyokuto Pharmaceutical, Japan). The V3–V4 regions of bacterial 16S rRNA genes were amplified using the Pro341F/Pro805R primers (Takahashi et al., 2014) and dual-index method (Hisada et al., 2015) under hemi-nested PCR conditions (Hell et al., 2013). Barcoded amplicons were paired-end sequenced on a 2 × 284-bp cycle using the MiSeq system with MiSeq Reagent Kit chemistry, version 3 (600 Cycle). Paired-end sequencing reads were merged using the fastq-join program with default settings (Aronesty, 2013). Only reads with quality value (QV) scores of ≥ 33 were extracted with split_libraries_fastq.py command in QIIME, version 1.8.0 (Caporaso et al., 2010). Chimeric sequences were removed using USEARCH61 (Edgar et al., 2011) with the identify_chimeric_seqs.py command in QIIME (Caporaso et al., 2010). Operational taxonomic units (OTUs) were aligned using the pick_open_reference_otus.py command in QIIME (Caporaso et al., 2010). OTUs with 97% similarity were identified with the Greengenes database, version 13.8 (Aitchison, 1986). Alpha diversity indices (e.g., observed species, Chao-1, Shannon, and PD_whole_tree) and beta diversity indices [e.g., principal coordinates analysis (PCoA)] were analyzed using the alpha_rarefaction.py and beta_diversity.py commands in QIIME, respectively (Edgar et al., 2011). Sampling depth in alpha and beta diversity was 50,895, which was the minimum number of read counts among samples. To account for compositional artifacts, we transformed relative abundances using the Centered Log-Ratio (CLR) transformation (Aitchison, 1986). The Chao-1 index was used to determine community richness and the Shannon index was used to determine community diversity. The PD_whole_tree index was used to compute Faith’s phylogenetic diversity. PCoA was used to show differences between the two groups. Unweighted UniFrac metrics were used for beta diversity (DeSantis et al., 2006). In order to identify the distinct gut microbiota between ICI responders and non-responders, the Linear Discriminative Analysis Effect Size (LEfSe) method was used to compare the composition of the gut microbiota in genus level using an online tool (http://www.ehbio.com/ImageGP/index.php/Home/Index/LEFSe.html).
Statistical analysis
Categorical variables were analyzed using Fisher’s exact test. Continuous variables were compared using the chi-squared test. The Mann-Whitney U test was used to determine significant differences among the different groups using alpha diversity, which showed the diversity in each individual sample. Logistic regression analysis to calculate odds ratios for ICI response with respect to clinic-pathological characteristics was used. Kaplan-Meier statistics and log-rank testing to evaluate progression-free survival (PFS) was applied.
Statistical analyses were performed using JMP software, version 14.0 (SAS Institute, Inc., Cary, NC, United States). The adonis function in the vegan package of R software, version 3.6.1, was used to conduct permutational multivariate analysis of variance (PERMANOVA) with respect to microbiome composition. p-values < 0.05 were considered statistically significant.
Results
Patient characteristics
The results were determined in follow-up examinations over a mean duration of 598 days (range, 81–1,225 days) after initial ICI therapy. Patient characteristics are shown in Supplementary Table S1. The study group included seven women and 21 men, with a mean age at diagnosis of 71 years (range, 56–88 years). Fifteen patients (53.6%) had ECOG performance status (PS) 0 and 13 (46.4%) had ECOG-PS 1. Seven patients (25.0%) had never smoked, and the remaining 21 patients were current or former smokers. The histological types were adenocarcinoma in 16 patients (57.1%) and squamous cell carcinoma in 12 patients (42.9%). Of the 28 included patients, one (3.6%) had stage IIA, seven (25.0%) had unresectable stage III (two with IIIA, three with IIIB, and two with IIIC), 11 (39.3%) had stage IV (seven with IVA and four with IVB), and nine (32.1%) had postoperative recurrence. Seven patients (25.0%) had mutant epithelial growth factor receptor (EGFR) and 21 patients (75.0%) had wild-type EGFR or no data regarding EGFR status. ICI was first-line therapy in 11 patients (39.3%), second-line therapy in 10 patients (35.7%), third-line therapy in five patients (17.8%), and fourth-line or later therapy in two patients (7.2%). Twenty-four patients (85.7%) received ICI monotherapy (nivolumab: 3 mg/kg or 240 mg/body intravenously at 2-week intervals; pembrolizumab: 200 mg/body intravenously at 3-week intervals; atezolizumab: 1,200 mg/body intravenously at 3-week intervals; or durvalumab: 10 mg/kg intravenously at 2-week intervals), while the remaining four patients received ICI therapy combined with platinum-doublet chemotherapy [pembrolizumab: 200 mg/body plus carboplatin, area under the curve (AUC) for concentration-time: 5 mg/ml/minute; plus pemetrexed 500 mg/m2 or nab-paclitaxel 200 mg/m2 intravenously at 3-week intervals]. Ten patients (35.7%) had more than 50% tumorous PD-L1 expression, 10 patients (35.7%) had 1%–49% PD-L1 expression, and eight patients (28.6%) had no PD-L1 expression or no data regarding PD-L1 status.
Immune checkpoint inhibitors response in non-small cell lung cancer patients
No patient experienced a complete response (CR) (0%), 17 patients had a partial response (PR) (60.7%), two patients had stable disease (SD) (7.2%), and nine patients (32.1%) had progressive disease (PD). Therefore, 17 patients (PR) were regarded as ICI responders; the remaining 11 patients (SD or PD) were regarded as ICI non-responders (Supplementary Figure S1). Logistic regression analysis to calculate odds ratios for ICI response with respect to clinic-pathological characteristics such as gender, smoking status, number of prior systemic therapy, regimen, histology and tumor proportional score, was performed (Supplementary Table S2). As a result, there was no significant difference between ICI response and any clinic-pathological features in multivariate analysis.
Relative abundances in oral and gut microbiomes
We analyzed the relative abundances of oral bacteria at the phylum and genus levels. At the phylum level, p_Firmicutes, p_Bacteroides, p_Proteobacteria, p_Actinobacteria, p_Fusobacteria, and p_TM 7 were the main taxa; these taxa comprised more than 99% in all groups. At the genus level, 169 genera were detected in the ICI responder saliva microbiota, while 152 species were detected in the ICI non-responder saliva microbiota (Figures 2A,B). Among the genus taxa with < 1% relative abundance, 157 genera (92.9%) were identified in the ICI responder saliva microbiota, while 138 genera (90.8%) were identified in the ICI non-responder saliva microbiota. Genera in both groups mainly included g_Streptococcus, g_Veillonella, g_Prevotella, g_Haemophilus, and g_Neisseria (Figures 2C,D). Regarding the gut microbiota, at the phylum level, p_Firmicutes, p_Bacteroides, p_Actinobacteria, p_Proteobacteria, p_Fusobacteria, and p_Verrucomicrobia were the main taxa (comprising more than 99%) in ICI responders. p_Firmicutes, p_Bacteroides, p_Actinobacteria, p_Proteobacteria and p_Fusobacteria were the main taxa (comprising more than 99%) in ICI non-responders. At the genus level, 180 genera were detected in the ICI responder gut microbiota, while 136 genera were detected in the ICI non-responder gut microbiota (Figures 2E,F). Among the genus taxa with < 1% relative abundance, 167 genera (92.8%) were identified in the ICI responder gut microbiota, while 122 genera (89.7%) were identified in the ICI non-responder gut microbiota. Genera in ICI responders mainly included g_Bacteroides, f_Ruminococcaceae; g_, f_Lachnospiraceae; g_, g_Streptococcus, and g_Blautia; genera in ICI non-responders mainly included g_Bacteroides, f_Ruminococcaceae; g_, f_Lachnospiraceae; g_, g_Prevotella, and g_Bifidobacterium (Figures 2G,H).
FIGURE 2
Diversity metrics in oral and gut microbiomes
We used alpha and beta diversity indices to evaluate the intersample and intrasample relationships in the oral and gut microbiota. Figures 3, 4 show the alpha diversity metrics of oral and gut microbiota at the genus level. In the oral microbiota, there were no significant differences between the two groups (Chao 1, p = 0.632; Observed species, p = 0.523; PD_whole_tree, p = 0.471; and Shannon, p = 0.474) (Figures 3A–D). In the gut microbiota, all alpha diversity metrics were significantly higher in ICI responders than in ICI non-responders (Chao 1, p = 0.017; Observed species, p = 0.024; PD_whole_tree, p = 0.022; and Shannon, p = 0.036) (Figures 3E–H). PCoA assessment of beta diversity was conducted based on weighted UniFrac distance (Figure 4). In both oral and gut microbiota, similar patterns were evident in the two groups. The greatest variations in the oral and gut microbiota of the two groups were 4.39% (PC1) and 4.32% (PC2), and 4.51% (PC1) and 4.22% (PC2), respectively. PERMANOVA based on unweighted UniFrac distance confirmed significant differences between the two groups in the gut microbiota alone (Figure 4A, oral: p = 0.904 and Figure 4C, gut: p = 0.035). Additionally, The Mann-Whitney U test showed significant differences in the gut microbiota alone (Figure 4B, oral: p = 0.760 and Figure 4D, gut: p = 0.005).
FIGURE 3
FIGURE 4
Linear discriminative analysis effect size results
We used LEfSe to perform high-dimensional genus comparisons regarding oral and gut microbiota between ICI responders and non-responders. Figure 5 shows that the ICI responder gut microbiota was significantly enriched for g_Blautia, compared with the ICI non-responder gut microbiota. In contrast, the ICI non-responder gut microbiota was significantly enriched for o_RF32 unclassified, compared with the ICI responder gut microbiota. There were no significant differences in the oral microbiota between the two groups.
FIGURE 5
Analyses for progression-free survival and correlation between the microbiome and immune adverse events
During a median follow-up of 598 days, 9 of 28 patients had immune-related adverse events (irAE). Six patients had hypothyroidism, 2 interstitial pneumonitis, 1 thyroiditis, 1 colitis and 1 hypoadrenalism. We applied Kaplan-Meier statistics and log-rank testing to evaluate PFS by ICI response, irAE, relative abundance of g. Blautia and o. RF32 unclassified (Figure 6). Median PFS was not reached in ICI responders, but median PFS was 92 days for non-responders (p < 0.0001). Median PFS was not reached in patients with irAE, but median PFS was 549 days for patients without irAE (p = 0.0581). Median PFS was not reached in patients enriched gut microbiota of g. Blautia, but median PFS was 549 days for patients reduced gut microbiota of g. Blautia (p = 0.0480). Median PFS was 49 days in patients with gut microbiota of o. RF32 unclassified, but median PFS was 757 days for patients without gut microbiota of o. RF32 unclassified (p = 0.0205).
FIGURE 6
Discussion
Similar studies regarding ICI therapy have focused on the microbiome in patients who have not yet received ICI therapy. Those studies revealed that the microbiota diversity and composition before ICI therapy was a predictive biomarker for ICI response. Although variations in gut microbiota composition were observed in the previous studies, there has been minimal information regarding gut microbiota status during ICI therapy. This information is important for efforts to enhance ICI therapy through biotics therapy (e.g., pre-, pro-, and synbiotics) and/or FMT. This analysis of oral and gut microbiota profiles in Japanese NSCLC patients during ICI therapy produced several novel findings.
The gut microbiota might have an important role in the ICI response in NSCLC patients, although the oral microbiota conveyed information distinct from the gut microbiota. Greater numbers of gut microbiota species were observed in ICI responders than in ICI non-responders. Additionally, those microbiota mainly consisted of minor species (<1%) at the genus level. Moreover, the fourth and fifth majority of gut microbiota species in ICI responder were differed from that in ICI non-responder, while first to fifth majority of oral microbiota were same between two groups.
An important tool for objective evaluation of the above data involves analysis of microbiota diversity: the numbers or abundances of microorganisms colonizing the gut. Greater alpha diversity indicates larger numbers of species in the gut, which implies a distinct gut microbiota composition. Several studies have reported that the ICI response was influenced by the alpha diversity of the gut microbiota before ICI therapy (Hamady and Knight, 2007; Gopalakrishnan et al., 2018). In the present study, high alpha diversity was observed in the ICI responder gut microbiota. Thus, our results indicated that ICI responders had more abundant gut microbiota, compared with ICI non-responders, during ICI therapy.
Our study also revealed a significant association between beta diversity and ICI response. Higher beta diversity indicates a significant difference in gut microbiota composition between two samples. In this study, the intersample distance was significantly shorter in ICI responders than in ICI non-responders. Therefore, the ICI responder gut microbiota had significant similarity, compared with the ICI non-responder gut microbiota, which implies a simple approach to control the gut microbiota by adding biotics therapy to ICI therapy.
Additionally, we identified specific gut microbiota species in Japanese NSCLC patients receiving ICI therapy by using LEfSe. The results indicated that g_Blautia was enriched in ICI responders, whereas o_RF32 unclassified was enriched in ICI non-responders. Several studies have shown significant associations between various gut microbiota components [e.g., Bifidobacterium longum (Fessler et al., 2018), Collinsella aerofaciens (Fessler et al., 2018), Enterococcus faecium (Fessler et al., 2018), Akkermansia muciniphila (Routy et al., 2018), Ruminococcaceae (Gopalakrishnan et al., 2018), Faecalibacterium genii (Jin et al., 2019) and Firmicutes (Froehlich et al., 2017)] and clinical response to ICI therapy. These results differed among individual studies, including the present study, presumably due to factors such as patient disease (e.g., melanoma, renal cell carcinoma, and NSCLC) and ICI regimen. Indeed, Chaput et al. (2017) indicated that the gut microbiota composition depended on the ICI regimen. Moreover, these differences might be caused by ethnicities. Nishijima et al. (2016) compared the gut microbiota between the Japanese population and individuals from 11 other nations. Notably, the Japanese gut microbiota was considerably different from the microbiota of other populations; it was characterized by the highest abundances of Blautia, Bifidobacterium, Collinsella, Streptococcus, and an unclassified Clostridiales genus, compared with the microbiota from the remaining 11 countries.
In addition to its status as a species characteristic of the gut microbiota in Japanese individuals, Blautia coccoides is regarded as an effective probiotic species. B. coccoides is an anaerobe and Gram-positive species found in human fecal samples. Reduced numbers of B. coccoides are associated with several benign diseases (e.g., hepatic cirrhosis and encephalopathy, irritable bowel syndrome, acute diarrhea, idiopathic inflammatory bowel disease, intestinal inflammation, and diabetes mellitus) and some malignancies (e.g., colorectal and breast cancers) (Ghanem et al., 2012; Liu and Tong, 2012). Furthermore, some reports have revealed that increased numbers of B. coccoides might be beneficial for human health. B. coccoides is increased among individuals with diets high in resistant starch and arabinoxylan (Nielsen et al., 2014); moreover, it can reduce NF-κB activity in human colon cancer cells (Tap et al., 2011). Additionally, Blautia obeum is a gut microbiota component involved in the transformation of carcinogenic heterocyclic amines, and reduced abundance of this species may increase heterocyclic amine-induced colorectal cancer risk (Ocvirk et al., 2019). Myles et al. (2014) reported that high omega-3 intake altered colonic inflammation and increased Blautia abundance in a murine model. Recently, Martini et al. (2022) reported that Blautia species could be a potential biomarker of outcome in metastatic colorectal cancer and non-small cell lung cancer patients treated with combined ICI therapy. Therefore, Blautia may be a key gut microbiota component involved in the ICI response in Japanese NSCLC patients. Conversely, the relative abundance of RF32 unclassified has been positively correlated with colonic damage and inflammation (Castro-Mejía et al., 2016), which are presumably negative influences on immunity in NSCLC patients.
The findings in this study indicate that controlling both gut microbiome diversity and the abundances of specific gut microbiome species during ICI therapy might lead to ICI response enhancement in Japanese NSCLC patients. Future therapies targeting the gut microbiome by means of pre-, pro-, or synbiotics to enhance the ICI response might be considered from our findings.
This study had the following limitations. First, only a single sample collection was performed. Thus, it was unclear how the microbiota diversity and composition might have changed before and after ICI therapy. Second, the sample size was small, which may have interfered with meaningful conclusions. Total 28 patients consisted of 18 with advanced stage and 10 with postoperative recurrence. Therefore, a statistical power was insufficient to analyze patients with each stages separately. Third, the delay between microbiota analysis and start of ICI therapy was almost 3 times longer in responders versus non responders, which might have any modifying effect on microbiota alpha and beta-diversity. Lastly, this study allowed various ICI regimens with or without combination chemotherapy. In the future, we plan to perform a large, multicenter, prospective observational study to evaluate the association between ICI response and changes in gut microbiota by collecting samples at multiple points (before and during ICI therapy) for NSCLC patients receiving a specific ICI regimen.
In conclusion, this study revealed a strong association between gut microbiota diversity and ICI response in Japanese NSCLC patients. Moreover, specific gut microbiota compositions may influence the ICI response. These findings might be useful in identifying biomarkers to predict ICI response, as well as in developing biotic therapies to enhance the ICI response.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
The studies involving human participants were reviewed and approved by the ethics committee of Kyushu Medical Center (approval number, 19C031) and the ethics committee of Kyushu Cancer Center (approval number, 2019-44). The patients/participants provided their written informed consent to participate in this study.
Author contributions
FS, MY, MO, ST, KY, TO, and YM: conception, design of the study protocol, data curation, methodology and writing original draft. FS, MY, KY, and TO: statistical analysis and interpretation of data. MY, MO, ST, KY, TO, and YM: supervision and writing original draft revision. All the authors read and approved the final paper.
Funding
This study was supported by JSPS KAKENHI Grant Number JP20K09188 and the Japanese Foundation for Multidisciplinary Treatment of Cancer (JFMC).
Acknowledgments
We thank Ryan Chastain-Gross, from Edanz (https://jp.edanz.com/ac) for editing a draft of this manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2022.1040424/full#supplementary-material
SUPPLEMENTARY FIGURE S1Study flow chart.
Abbreviations
CR, complete response; ECOG, European Clinical Oncology Group; EGFR, epidermal growth factor receptor; FMT, fecal microbiome transplantation; ICI, immune-check point inhibitors; irAE, immune-related adverse events; MST, median survival time; NSCLC, non-small cell lung cancer; PD, progressive disease; PD-L1, programmed death-ligand 1; PFS, progression free survival; PS, performance status; PR, partial response; SD, stable disease.
References
1
AitchisonJ. (1986). The statistical analysis of compositional data. London; New York: Chapman & Hall, 416.Google Scholar
2
AronestyE. (2013). Comparison of sequencing utility programs. Open Bioinforma. J.7, 1–8. 10.2174/1875036201307010001CrossRef Full Text | Google Scholar
3
AtarashiK.SudaW.LuoC.et alEctopic colonization of oral bacteria in the intestine drives TH1 cell induction and inflammation. Science. 2017; 358: 359–365. PubMed Abstract | CrossRef Full Text | Google Scholar
4
CaporasoJ. G.KuczynskiJ.StombaughJ.BittingerK.BushmanF. D.CostelloE. K.et al (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods7, 335–336. 10.1038/nmeth.f.303PubMed Abstract | CrossRef Full Text | Google Scholar
5
Castro-MejíaJ.JakesevicM.KrychŁ.et al (2016). Treatment with a Monoclonal Anti-IL-12p40 Antibody Induces Substantial Gut Microbiota Changes in an Experimental Colitis Model. Gastroenterol Res Pract2016, 4953120. 10.1155/2016/4953120PubMed Abstract | CrossRef Full Text | Google Scholar
6
ChaputN.LepageP.CoutzacC.SoularuEE.Le RouxK.MonotC.et al (2017). Baseline gut microbiota predicts clinical response and colitis in metastatic melanoma patients treated with ipilimumab. Ann. Oncol28, 1368–1379. 10.1093/annonc/mdx108PubMed Abstract | CrossRef Full Text | Google Scholar
7
CorralesL.WilliamsJ. B.Aquino-MichaelsK.EarleyZ. M.et al (2015). Commensal Bifidobacterium promotes antitumor immunity and facilitates anti-PD-L1 efficacy. Science350, 1084–1089. 10.1126/science.aac4255PubMed Abstract | CrossRef Full Text | Google Scholar
8
DavarD.DzutsevA. K.McCullochJ. A.RodriguesR. R.ChauvinJ. M.MorrisonR. M.et al (2021). Fecal microbiota transplant overcomes resistance to anti-PD-1 therapy in melanoma patients. Science371, 595–602. 10.1126/science.abf3363PubMed Abstract | CrossRef Full Text | Google Scholar
9
DeSantisT. Z.HugenholtzP.LarsenN.RojasM.BrodieE. L.KellerK.et al (2006). Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl. Environ. Microbiol72, 5069–5072. 10.1128/AEM.03006-05PubMed Abstract | CrossRef Full Text | Google Scholar
10
EdgarR. C.HaasB. J.ClementeJ. C.QuinceC.KnightR. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics27, 2194–2200. 10.1093/bioinformatics/btr381PubMed Abstract | CrossRef Full Text | Google Scholar
11
EisenhauerE. A.TherasseP.BogaertsJ.SchwartzL. H.SargentD.FoRdR.et al (2009). New response evaluation criteria in solid tumours: Revised RECIST guideline (version 1.1). Eur. J. Cancer45, 228–247. 10.1016/j.ejca.2008.10.026PubMed Abstract | CrossRef Full Text | Google Scholar
12
FesslerJ.ChongsuwatT.ZhaY.AlegreM. L.et al (2018). The commensal microbiome is associated with anti-PD-1 efficacy in metastatic melanoma patients. Science359, 104–108. 10.1126/science.aao3290PubMed Abstract | CrossRef Full Text | Google Scholar
13
FroehlichT. W.XieY.FrenkelE. P.et al (2017). Metagenomic shotgun sequencing and unbiased metabolomic profiling identify specific human gut microbiota and metabolites associated with immune checkpoint therapy efficacy in melanoma patients. Neoplasia19, 848–855. 10.1016/j.neo.2017.08.004PubMed Abstract | CrossRef Full Text | Google Scholar
14
GandhiL.Rodríguez-AbreuD.GadgeelS.EstebanE.FelipE.De AngelisF.et al (2018). Pembrolizumab plus chemotherapy in metastatic non-small-cell lung cancer. N. Engl. J. Med378, 2078–2092. 10.1056/NEJMoa1801005PubMed Abstract | CrossRef Full Text | Google Scholar
15
GhanemK. Z.BlautM.ClavelT.LohG. (2012). Lignan transformation by gut bacteria lowers tumor burden in a gnotobiotic rat model of breast cancer. Carcinogenesis33, 203–208. 10.1093/carcin/bgr256PubMed Abstract | CrossRef Full Text | Google Scholar
16
GoldstrawP.ChanskyK.CrowleyJ.Rami-PortaR.AsamuraH.EberhardtW. E. E.et al (2016). The IASLC lung cancer staging project: Proposals for revision of the TNM stage groupings in the forthcoming (eighth) edition of the TNM classification for lung cancer. J. Thorac. Oncol11, 39–51. 10.1016/j.jtho.2015.09.009PubMed Abstract | CrossRef Full Text | Google Scholar
17
GopalakrishnanV.SpencerC. N.NeziL.ReubenA.AndrewsM. C.KarpinetsT. V.et al (2018). Gut microbiome modulates response to anti-PD-1 immunotherapy in melanoma patients. Science359, 97–103. 10.1126/science.aan4236PubMed Abstract | CrossRef Full Text | Google Scholar
18
HakozakiT.RichardC.ElkriefA.HosomiY.BenlaifaouiM.MimpenI.et al (2020). The gut microbiome associates with immune checkpoint inhibition outcomes in patients with advanced non-small cell lung cancer. Cancer Immunol. Res8, 1243–1250. 10.1158/2326-6066.CIR-20-0196PubMed Abstract | CrossRef Full Text | Google Scholar
19
HamadyM.KnightR. (2007). Quantitative and qualitative beta diversity measures lead to different insights into factors that structure microbial communities. Appl. Environ. Microbiol73, 1576–1585. 10.1128/AEM.01996-06PubMed Abstract | CrossRef Full Text | Google Scholar
20
HellK.EdwardsA.ZarskyJ.PodmirsegS. M.GirdwoodS.PachebatJ. A.et al (2013). The dynamic bacterial communities of a melting High Arctic glacier snowpack. ISME J7, 1814–1826. 10.1038/ismej.2013.51PubMed Abstract | CrossRef Full Text | Google Scholar
21
HisadaT.EndohK.KurikiK. (2015). Inter-and intra-individual variations in seasonal and daily stabilities of the human gut microbiota in Japanese. Arch. Microbiol.197, 919–934. 10.1007/s00203-015-1125-0PubMed Abstract | CrossRef Full Text | Google Scholar
22
HosgoodH. D.CaiQ.HuaX.LongJ.ShiJ.WanY.et al (2021). Variation in oral microbiome is associated with future risk of lung cancer among never-smokers. Thorax76, 256–263. 10.1136/thoraxjnl-2020-215542PubMed Abstract | CrossRef Full Text | Google Scholar
23
JinY.DongH.XiaL.YangY.ZhuY.ShenY.et al (2019). The diversity of gut microbiome is associated with favorable responses to anti-programmed death 1 immunotherapy in Chinese patients with NSCLC. J. Thorac. Oncol14, 1378–1389. 10.1016/j.jtho.2019.04.007PubMed Abstract | CrossRef Full Text | Google Scholar
24
LittmanD. R. (2012). Interactions between the microbiota and the immune system. Science336, 1268–1273. 10.1126/science.1223490PubMed Abstract | CrossRef Full Text | Google Scholar
25
LiuF.TongX. (2012). Human intestinal lumen and mucosa-associated microbiota in patients with colorectal cancer. PLoS One7, e39743. 10.1371/journal.pone.0039743PubMed Abstract | CrossRef Full Text | Google Scholar
26
MartiniG.CiardielloD.DallioM.FamigliettiV.EspositoL.CorteC. M. D.et al (2022). Gut microbiota correlates with antitumor activity in patients with mCRC and NSCLC treated with cetuximab plus avelumab. Int. J. Cancer151, 473–480. 10.1002/ijc.34033PubMed Abstract | CrossRef Full Text | Google Scholar
27
MylesI. A.PincusN. B.FontecillaN. M.DattaS. K. (2014). Effects of parental omega-3 fatty acid intake on offspring microbiome and immunity. PLoS One9, e87181. 10.1371/journal.pone.0087181PubMed Abstract | CrossRef Full Text | Google Scholar
28
NielsenT. S.LærkeH. N.TheilP. K.SorensenJ. F.SaarinenM.ForsstenS.et al (2014). Diets high in resistant starch and arabinoxylan modulate digestion processes and SCFA pool size in the large intestine and faecal microbial composition in pigs. Br. J. Nutr112, 1837–1849. 10.1017/S000711451400302XPubMed Abstract | CrossRef Full Text | Google Scholar
29
NishijimaS.SudaW.OshimaK.KimS. W.HiroseY.MoritaH.et al (2016). The gut microbiome of healthy Japanese and its microbial and functional uniqueness. DNA Res23, 125–133. 10.1093/dnares/dsw002PubMed Abstract | CrossRef Full Text | Google Scholar
30
OcvirkS.LiJ.MahachiC. B.SoniD.et al (2019). Differences in fecal gut microbiota, short-chain fatty acids and bile acids link colorectal cancer risk to dietary changes associated with urbanization among Zimbabweans. Nutr. Cancer71, 1313–1324. 10.1080/01635581.2019.1602659PubMed Abstract | CrossRef Full Text | Google Scholar
31
Paz-AresL.LuftA.VicenteD.TafreshiA.GumusM.MazieresJ.et al (2018). Pembrolizumab plus chemotherapy for squamous non-small-cell lung cancer. N. Engl. J. Med.379, 2040–2051. 10.1056/NEJMoa1810865PubMed Abstract | CrossRef Full Text | Google Scholar
32
PietrantoniC.AquilinaA. T.OkadaM.GroverV.et al (2003). Microbiology of severe aspiration pneumonia in institutionalized elderly. Am. J. Respir. Crit. Care Med.167, 1650–1654. 10.1164/rccm.200212-1543OCPubMed Abstract | CrossRef Full Text | Google Scholar
33
RibasA.WolchokJ. D. (2018). Cancer immunotherapy using checkpoint blockade. Science359, 1350–1355. 10.1126/science.aar4060PubMed Abstract | CrossRef Full Text | Google Scholar
34
RittmeyerA.BarlesiF.WaterkampD.ParkK.CiardielloF.von PawelJ.et al (2017). Atezolizumab versus docetaxel in patients with previously treated non-small-cell lung cancer (OAK): A phase 3, open-label, multicentre randomised controlled trial. Lancet389, 255–265. 10.1016/S0140-6736(16)32517-XPubMed Abstract | CrossRef Full Text | Google Scholar
35
RoutyB.Le ChatelierE.DerosaL.DuongC. P. M.AlouM. T.DaillereR.et al (2018). Gut microbiome influences efficacy of PD-1-based immunotherapy against epithelial tumors. Science359, 91–97. 10.1126/science.aan3706PubMed Abstract | CrossRef Full Text | Google Scholar
36
ShojiF.TakeokaH.KozumaY.ToyokawaG.YamazakiK.IchikiM.et al (2019). Pretreatment prognostic nutritional index as a novel biomarker in non-small cell lung cancer patients treated with immune checkpoint inhibitors. Lung Cancer136, 45–51. 10.1016/j.lungcan.2019.08.006PubMed Abstract | CrossRef Full Text | Google Scholar
37
SocinskiM. A.JotteR. M.CappuzzoF.OrlandiF.StroyakovskiyD.NogamiN.et al (2018). Atezolizumab for first-line treatment of metastatic nonsquamous NSCLC. N. Engl. J. Med378, 2288–2301. 10.1056/NEJMoa1716948PubMed Abstract | CrossRef Full Text | Google Scholar
38
TakahashiS.TomitaJ.NishiokaK.HisadaT.NishijimaM. (2014). Development of a prokaryotic universal primer for simultaneous analysis of Bacteria and Archaea using next-generation sequencing. PLoS One9, e105592. 10.1371/journal.pone.0105592PubMed Abstract | CrossRef Full Text | Google Scholar
39
TapJ.Beguet-CrespelF.Le RouxK.de WoutersT.CultroneA.et al (2011). Identification of NF-κB modulation capabilities within human intestinal commensal bacteria. J. Biomed. Biotechnol2011, 282356. 10.1155/2011/282356PubMed Abstract | CrossRef Full Text | Google Scholar
40
VarlottoJ. M.RechtA.FlickingerJ. C.Medford-DavisL. N.DyerA. M.DeCampM. M. (2009). Factors associated with local and distant recurrence and survival in patients with resected nonsmall cell lung cancer. Cancer115, 1059–1069. 10.1002/cncr.24133PubMed Abstract | CrossRef Full Text | Google Scholar
41
VétizouM.PittJ. M.DaillèreR.et alAnticancer immunotherapy by CTLA-4 blockade relies on the gut microbiota, Science. 2015; 350: 1079–1084. 10.1126/science.aad1329PubMed Abstract | CrossRef Full Text | Google Scholar
42
YamasakiK.KawanamiT.YateraK.FukudaK.NoguchiS.NagataS.et al (2013). Significance of anaerobes and oral bacteria in community-acquired pneumonia. PLoS One8, e63103. 10.1371/journal.pone.0063103PubMed Abstract | CrossRef Full Text | Google Scholar
Summary
Keywords
oral and gut microbiota, diversity, specific composition, non-small cell lung cancer, immune checkpoint inhibitors response
Citation
Shoji F, Yamaguchi M, Okamoto M, Takamori S, Yamazaki K, Okamoto T and Maehara Y (2022) Gut microbiota diversity and specific composition during immunotherapy in responders with non-small cell lung cancer. Front. Mol. Biosci. 9:1040424. doi: 10.3389/fmolb.2022.1040424
Received
09 September 2022
Accepted
11 October 2022
Published
24 October 2022
Volume
9 - 2022
Edited by
Libo Wang, First Affiliated Hospital of Zhengzhou University, China
Reviewed by
Mark Uhlenbruch, Kaiserswerther Diakonie, Germany
Thomas Kaiser, University of Minnesota Twin Cities, United States
Updates
Copyright
© 2022 Shoji, Yamaguchi, Okamoto, Takamori, Yamazaki, Okamoto and Maehara.
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: Fumihiro Shoji, fumshojifumshoji@gmail.com
This article was submitted to Molecular Diagnostics and Therapeutics, a section of the journal Frontiers in Molecular Biosciences
Disclaimer
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.