- 著者: Po-Yu Chien, Wen-Chien Cheng, Chin-Chuan Hung, Chih-Yen Tu, Te-Chun Hsia, Der-Yang Cho, Po-Ren Hsueh, Zi-Lun Lai, Yi-Cheng Shen, Yu-Chao Lin
- Corresponding author: Yi-Cheng Shen, Yu-Chao Lin (Division of Pulmonary and Critical Care, China Medical University Hospital, Taichung, Taiwan)
- 雑誌: Cancer Immunology, Immunotherapy
- 発行年: 2026
- Epub日: N/A
- Article種別: Original Article
- DOI: 10.1007/s00262-026-04463-3
背景
非小細胞肺癌 (NSCLC) は全肺癌の約85%を占め、世界的に癌関連死亡率の主要因として認識されている。免疫チェックポイント阻害薬 (ICI) は抗腫瘍免疫監視機構を回復させることで一部の患者に持続的な臨床利益をもたらしたが、奏効は個人差が極めて大きく、一次耐性・獲得耐性・免疫関連有害事象が依然として課題である。
喫煙歴はNSCLCの主要な発癌リスク因子であるが、paradoxicalに喫煙者はICIによってより良好な臨床転帰を示すことが繰り返し報告されている。これは腫瘍変異量 (tumor mutational burden, TMB) の高さやネオアンチゲン量の増加による免疫原性の亢進に起因すると考えられている。さらに、本来ICIへの反応性が低いEGFR変異NSCLCでも、喫煙者ではより顕著な免疫活性化シグネチャーを示すことが報告されており、喫煙誘発性の腫瘍免疫微小環境の変化が関与する可能性が示唆されている。
腸内細菌叢はICIの有効性を規定する重要な宿主因子として確立されており、樹状細胞の成熟・抗原提示・T細胞プライミングを介した全身免疫調節を担う。Akkermansia muciniphilaやBifidobacterium属、Faecalibacterium属などの特定菌種がPD-1遮断の効果を増強することが示されており (Hakozaki et al. CancerImmunolRes 2020)、腸-肺軸 (gut-lung axis) がNSCLCの免疫療法において重要な役割を担う概念が確立しつつある (Zhu et al. CancerMed 2026)。
一方、喫煙はニコチン・多環芳香族炭化水素・揮発性有機化合物などを介して腸粘膜の完全性・腔内pH・抗菌圧を変化させ、Bifidobacterium・Faecalibacterium の減少やPrevotella・Escherichia-Shigella の増加など、腸内細菌叢の恒常性を障害することが知られている。しかしながら、喫煙誘発性の腸内細菌叢変化がICIアウトカムにどう寄与するかは未解明であり、喫煙という交絡因子を考慮した系統的な解析が不足していた。特にBifidobacterium longumの役割については、Jin et al. JThoracOncol 2019 がNSCLC PD-1治療の奏効者でB. longumの濃縮を報告する一方、Routy et al. (Science 2018) はNSCLC患者においてB. longum低量群でPFS延長を報告するなど矛盾した知見が残されており、特に喫煙・EGFR変異・腸内細菌叢・ICI奏効の四者関連はこれまで報告されていない。本研究はこのギャップを埋めるべく実施された。
目的
NSCLCにおける喫煙状況・腸内細菌叢特徴・ICI治療後の臨床転帰の三者関連を前向きコホートで解明すること。喫煙関連の腸内細菌叢変化が差次的なICI奏効に寄与する微生物・機能的メカニズムを明らかにし、腸内細菌叢を標的とした新たな治療戦略の開発に向けた基盤知識を提供する。
結果
患者背景と登録コホートの特性： 271名がスクリーニングされ、ベースライン便検体欠損や臨床データ不完全の46名を除外した225名が最終解析に含まれた (Fig. 1)。患者は喫煙状況により非喫煙者 (n=121) と喫煙者 (n=104) に分類された。コホート全体の中央値年齢は67歳 (範囲28-93歳) で、男性が56.0%を占めた (Table 1)。病期IV疾患が68.4%と多く、組織型は腺癌が80.9%と最も多く、続いて扁平上皮癌が14.2%、腺扁平上皮癌が2.2%を占めた。平均BMIは23.46 (範囲14.89-36.59) であった。167例でEGFR変異検査が実施され、75例 (33.3%) に変異を認め、L858Rが57.3%、Del-19 (exon 19 deletion、exon 19欠失) が32.0%、稀少変異が10.7%を占め、野生型EGFRは40.9% (92例) であった。PD-L1 tumor proportion score (TPS) は≥50%が32.4%、1-49%が37.8%、<1%が29.7%であった。ICI投与を受けたのは52名 (Table 2) で、ペムブロリズマブが86.5%と最も多く使用され、アテゾリズマブ5.8%、ニボルマブ3.8%が続いた。ICI治療コホートは男性優位 (75.0%)、喫煙者優位 (73.1%) であり、腺癌55.8%・扁平上皮癌36.5%が主要組織型であった。
喫煙状況による腸内細菌叢組成の変化： 全検体にわたり26,879のASVs (amplicon sequence variants、増幅配列変異体) が同定され、うち5,627が両群間で共有され、11,367が非喫煙者固有、9,885が喫煙者固有であった (Fig. 2A)。属レベルでは、Bacteroides・Prevotella・Sutterella・Faecalibacterium・Escherichia-Shigella・Bifidobacterium・Akkermansia などの優占分類群の相対量に両群間で統計学的有意差は認められなかった (Fig. 2B)。しかし喫煙者では、dysbiosis・炎症促進状態と関連するPrevotella・Sutterella・Escherichia-Shigella の増加傾向とともに、Bifidobacterium・Faecalibacterium・Akkermansia の減少傾向が一貫して観察された。Alpha多様性はShannon entropy (p=0.755) ・observed OTUs (operational taxonomic units、操作的分類単位; p=0.901) ともに非喫煙者 (n=121) vs 喫煙者 (n=104) の間で有意差なく (Fig. 2C-D)、Bray-Curtis非類似度に基づくBeta多様性もPERMANOVA (permutational multivariate analysis of variance、置換多変量分散分析) でp=0.175と有意な群間差は認められなかった (Fig. 2E)。以上より喫煙は全体的な微生物多様性・群集構造を大きく変化させないことが示された。LEfSe (linear discriminant analysis effect size、線形判別分析効果量; LDAスコア≥2) による選択的taxa差異解析では、非喫煙者ではClostridia・Peptostreptococcales 内のtaxaが濃縮され、喫煙者ではCoriobacteriales Incertae Sedis・Streptococcaceae が高頻度であった (Fig. 2F)。
Del-19/喫煙交差解析によるB. longumの同定と機能経路の予測： exon 19欠失 (Del-19、ICIへの反応性が低いEGFR変異サブタイプ) 腫瘍で豊富に存在し、かつ喫煙者群で枯渇するtaxaを同定するoverlap解析を実施した (Fig. 3A-B)。この交差解析から、Bifidobacterium longumがDel-19変異で増加し喫煙者で有意に減少する主要taxonとして同定された (Fig. 3C)。次にPICRUSt2を用いた予測メタゲノム経路解析を行ったところ、喫煙者ではLPS生合成・Kdo₂-lipid A生合成・lipid IVa生合成など複数のLPS関連経路が顕著に濃縮された (Fig. 3D)。LPS産生を担う主体がGram陰性菌であるため、B. longum量とGram陰性菌割合の関係を解析したところ逆相関が確認され (Spearman相関解析、Fig. 3E)、B. longum減少がGram陰性菌優位な細菌叢配置と連動することが示された。さらに、Gram陰性菌割合とlpxM (hexa-acylated lipid A生成に関わる脂質Aミリストイル転移酵素をコード) の予測量の間には正相関が認められた (Fig. 3F)。lpxMは主としてGammaproteobacteria クラスのメンバーがコードするため、GammaproteobacteriaとICI奏効の関連をJonckheere-Terpstra検定で評価したところ、Gammaproteobacteria 高量がより良好なICI奏効 (PD→SD→PR) と有意に単調増加関係を示した (Fig. 3G)。
B. longum量とICI奏効・PFSの多変数解析： ICI治療コホート (n=52) において、年齢・BMI (body mass index)・性別・PD-L1発現・喫煙歴・食事習慣 (野菜摂取・果物摂取・プロバイオティクス飲料摂取・揚げ物摂取) を共変量とした多変数cumulative-logit比例オッズモデルにてB. longum量とICI奏効の関連を評価した。B. longum高量は交絡因子調整後もICI奏効不良と独立して関連した (OR 0.44、95% CI 0.22-0.89、p=0.022；Fig. 4A、Table 3)。すなわちB. longum低量がICI奏効良好の独立した微生物予測因子であることが示された。その他の有意な独立予測因子として年齢 (OR 1.18、95% CI 1.03-1.36、p=0.021) およびBMI (OR 0.58、95% CI 0.38-0.89、p=0.012) が同定され、果物摂取量も有意な正の関連を示した (OR 76.53、95% CI 1.43-4105.77、p=0.033)。喫煙変数自体はOR 13.53 (95% CI 0.47-392.42、p=0.129) と有意でなく、本コホートでは喫煙状況単独ではICI奏効改善と関連しなかった。
PFS解析はペムブロリズマブを受けたStage IVのNSCLC患者に限定して実施した。B. longum低量群 (n=19、進行イベント9件) vs 高量群 (n=16、進行イベント13件) の比較で、低量群のPFSが有意に延長した (log-rank p=0.044；Fig. 4B)。進行イベント発生率は低量群9/19 (47.4%) vs 高量群13/16 (81.3%) と低量群で顕著に低く、B. longumはNSCLCにおけるICI治療、特にペムブロリズマブ投与患者において機能的な負の予測バイオマーカーとして位置づけられる可能性が示唆された。
考察/結論
本研究はNSCLC患者の前向きコホートにおいて、喫煙関連腸内細菌叢の変容をICIアウトカムと結びつける統合的な証拠を提供した。喫煙は全体的な腸内微生物多様性を大きく変化させないものの、B. longum相対量の選択的減少とGram陰性菌 (特にGammaproteobacteria) の濃縮をもたらし、lpxMを介したhexa-acylated LPS産生増強と関連することが示された。
① 先行研究との違い: B. longumの役割については先行研究で矛盾した報告が存在する。Jin et al. JThoracOncol 2019 はPD-1阻害薬奏効者でB. longumの濃縮を報告し、Komatsu et al. (Sci Rep 2025) も奏効者でBifidobacteriaceae高量を報告している。しかし本研究はこれらと対照的に、B. longum低量がICI奏効良好と関連するという結果を示した。この相違は方法論的な差異によって説明できる。Jinらは部分奏効 (PR) と病勢安定 (SD) を統合した二項「奏効者」として分類したのに対し、本研究ではRECIST基準に基づく段階的順序アウトカム (PD < SD < PR) を用いた多変数比例オッズモデルを採用した。またKomatsuらが科 (family) レベルで評価したのに異なり、本研究では全長16S rRNA解析により種レベル精度でB. longumを特定した。Routy et al. (Science 2018) のNSCLC患者でのB. longum低量群PFS延長という報告は本研究の所見と整合的であり、本研究はより精密な分類学的・統計学的手法でこの関係を確認した形となる。
② 新規性: 本研究で初めて、喫煙→B. longum減少→Gram陰性菌 (Gammaproteobacteria) 優位な細菌叢配置→lpxM高発現→hexa-acylated LPS産生増加→TLR4を介した自然免疫活性化→ICI増感という一貫した微生物軸を前向きコホートで新規に実証した。喫煙がICIを増強するメカニズムとして従来はTMBの高さが主要因と考えられてきたが、本研究は腸内細菌叢が喫煙と免疫療法奏効を繋ぐ機能的中間因子として新規の役割を担うことを示した。さらに、B. longumとEGFR Del-19変異が同じ微生物シグネチャー (B. longum高量) を共有するという観察は、腫瘍ドライバー変異ごとに異なる腸内細菌叢-免疫軸が存在する可能性という新規の研究概念を開拓するものである。喫煙関連微生物変化が免疫療法奏効に対して有利な細菌叢配置を誘導しうるというパラダイム転換的示唆も、本研究の重要な新規性として位置づけられる。
③ 臨床応用: 本研究の知見は複数の臨床応用可能性を示唆する。第一に、B. longumは既存のPD-L1 TPS単独では不十分なICI奏効予測を補完する微生物バイオマーカーとしての臨床的有用性が期待される。第二に、BMIや食事習慣 (果物摂取量) もICI奏効と関連することが示されており、代謝・栄養状態と腸内細菌叢情報を統合した多変数バイオマーカーパネルの開発が臨床に資する可能性がある。第三に、喫煙の有害作用なしに有益な微生物変化—B. longum枯渇・lpxM/Gammaproteobacteria濃縮—を再現する腸内細菌叢調節戦略 (特定菌種の選択的制御・lipid A生合成経路の標的化・糞便微生物移植) がICI奏効改善のための臨床応用として期待される。特に、Gammaproteobacteria由来のhexa-acylated LPSがICI奏効を増強するという最近の実験的知見 (Sardar et al. Nat Microbiol 2025) を踏まえれば、lipid A acylation patternの制御が具体的な治療介入点となりうる。
④ 残された課題: 今後の課題として以下が重要である。第一に、全長16S rRNA解析は株レベルの多様性を捉えることができず、PICRUSt2による機能予測は実際の遺伝子発現・代謝産物産生・宿主免疫応答の直接的証拠ではない。今後はショットガンメタゲノミクス・メタボロミクスによる実測データが必要である。第二に、食事習慣の詳細な栄養素組成・総カロリー・縦断的変動が評価されておらず、食事由来の残余交絡は排除できない。第三に、ICI治療サブコホートはn=52と小規模であり、PFS解析はペムブロリズマブ治療例のみに限定されているため、外的妥当性の確認に独立コホートでの検証が必須である。limitation として単施設からのサンプルであること、食事介入前の詳細な生活習慣データが限られること、喫煙変数自体はICI奏効と有意に関連しなかった点 (異質性・小サンプルの影響) なども挙げられ、今後の多施設研究での克服が期待される。
方法
研究デザインと対象: 台湾・China Medical University HospitalにおいてNSCLC患者を前向き登録するコホート研究 (IRB承認取得; NCBI SRA: PRJNA1388452)。271名をスクリーニングし、ベースライン便検体欠損・臨床データ不完全な46名を除外した225名を最終解析対象とした (ヘルシンキ宣言準拠、全例文書同意取得)。ICI治療サブコホートはn=52で、ペムブロリズマブが86.5%を占めた。PFS解析はペムブロリズマブ治療Stage IV患者に限定。
マイクロバイオーム解析: 便検体はICI開始前に採取し、QIAamp PowerFecal Pro DNA Kit (QIAcube HT自動抽出系) でゲノムDNAを抽出 (最低収量500 ng/sample)。PacBio barcoded primer 27F/1492R による全長16S rRNA遺伝子増幅後、PacBio Sequel IIシーケンサでシーケンス実施。DADA2によるASVs生成・キメラ除去を行い、VSEARCH/SILVA v138データベースでtaxonomic分類した。Alpha多様性 (Shannon指数・observed OTUs、Kruskal-Wallis検定)、Beta多様性 (Bray-Curtis非類似度、PERMANOVA) を算出し、差次的豊富性はLEfSe (LDAスコア≥2) で評価した。機能的メタゲノム経路予測はPICRUSt2で実施し、STAMP で可視化・統計比較した。多重比較補正はBenjamini-Hochberg FDR法を使用した。
統計解析: 喫煙者vs非喫煙者の微生物量比較はWhiteの非パラメトリックt検定 (STAMP)、EGFR亜型間の比較はANOVAを使用した。ICI治療コホートでは臨床奏効 (PD < SD < PR) を順序アウトカムとして、cumulative-logit比例オッズモデル (clm関数、R ordinalパッケージ) による多変数解析を実施し、年齢・性別・BMI・喫煙状況・PD-L1発現・食事習慣を共変量として調整した。PFS解析はKaplan-Meier法とlog-rank検定を適用した (制限: Stage IV・ペムブロリズマブ治療患者)。全統計はR 4.4.1 (tidyverse, dplyr, ggplot2, survival, survminer) およびGraphPad Prism v5を使用し、両側p<0.05を有意水準とした。NCBIシーケンスリードアーカイブ登録番号: PRJNA1388452。