多発性骨髄腫（MM）は、骨髄内において単一クローンの形質細胞が異常増殖する造血器腫瘍です。多くの患者において、無症状の前駆病態である意義不明の単クローン性ガンマグロブリン血症（MGUS）やくすぶり型多発性骨髄腫（SMM）を経て、治療を要する初発多発性骨髄腫（NDMM）、さらには再発・難治性多発性骨髄腫（RRMM）へと段階的に進行します。これまで病態解明の研究は主に腫瘍細胞自体の遺伝子変異に焦点を当てて行われてきましたが、疾患の進行や治療抵抗性の獲得には、腫瘍細胞を取り巻く腫瘍微小環境（TME）の改変が不可欠な要素として深く関与していることが明らかになってきました。
本研究では、疾患の全スペクトラムにわたる235例の患者から採取された骨髄サンプルを用い、単一細胞RNAシーケンス（scRNA-seq）、T細胞受容体シーケンス（scTCR-seq）、B細胞受容体シーケンス（scBCR-seq）を統合した大容量の単一細胞マルチオミクス解析を実施しました。計968,681個の細胞（TME細胞563,957個、形質細胞107,310個）を詳細に解析した結果、骨髄微小環境に存在する免疫細胞および間質細胞の多様な集団動態と、機能障害に至る分子機構が網羅的に解明されました。
キーポイント
- 5つの再現可能なTME「エコタイプ」の確立: 非負値行列因子分解（NMF）を用いた細胞構成の解析により、骨髄微小環境内の免疫・間質細胞が形成する協調的な組織化構造（E1〜E5）を導き出しました。
- 疾患病期を超えた微小環境の構造化と分子ネットワーク: エコタイプ分類は従来の臨床病期（MGUS、SMM、NDMM、RRMM）の枠組みを超えた微小環境の不均一性を捉えており、特定の細胞間コミュニケーションや転写因子ネットワーク、共通遺伝子モジュールと密接に連動しています。
- ストレス応答に伴うエネルギー代謝不全と免疫機能障害: 疾患の進行に伴い出現する「ストレス関連T/NK細胞」は、AP-1やNF-κBの亢進、小胞体ストレス応答を特徴とし、酸化的リン酸化（ATP合成）の低下という代謝不全と直接結合した機能低下状態にあることが見出されました。
- 臨床予後および高度免疫療法への不応答性の予測価値: エコタイプは初発患者の生存期間（PFS/OS）と関連するとともに、T細胞の機能的整合性（フィットネス）の欠如を示すストレス型エコタイプ（E5）は、CAR-T療法やT細胞誘導抗体（TCE）療法における治療不応答の主要な要因となることが実証されました。
5つのTME「エコタイプ」による骨髄微小環境の体系的分類
腫瘍微小環境（TME）を体系的に理解するためには、単一の免疫細胞サブセットの増減を個別に評価するだけでは不十分です。多発性骨髄腫の骨髄内では、多様な細胞成分が密接に相互作用し合いながら、ひとつの協調的な細胞社会（エコタイプ）を形成しています。細胞成分の構造化された集団としてTMEを定義することは、病態の解明および患者ごとの微小環境の不均一性を解き明かす上で重要な視点となります。
本論文では、患者ごとの細胞サブセット構成比率に対して非負値行列因子分解（NMF）を適用し、再現性の高い5つのTMEエコタイプ（E1〜E5）を識別しました。各エコタイプは独自の細胞構成とシグナル網を有しています。
- E1（前駆細胞富裕型）: 造血幹細胞・前駆細胞や未熟B細胞が豊富で、成熟免疫細胞の浸潤が抑制された微小環境です。マクロファージ等から腫瘍細胞や前駆細胞へのCXCL12-CXCR4 シグナルに加え、ICAM、VCAM、PECAM、BTLA、CD200、SIRPAなどの接着・免疫チェックポイント分子の発現が維持されています。
- E2（恒常的免疫監視型）: ナイーブ/メモリT細胞、制御性T細胞（Treg）、メモリB細胞、NK細胞を中心とする恒常的状態です。Midkineシグナルの選択的な更新が認められ、非炎症性の環境下でT細胞の生存が維持されています。
- E3（エフェクターT細胞優位型）: 細胞毒性を有するエフェクターCD8+ T細胞が集積した状態です。補体経路やncWNT経路とともに、CCL3/CCL5-CCR1を介した自然免疫系と獲得免疫系間のコミュニケーションが活発に行われています。
- E4（炎症性モノサイト優位型）: 炎症性モノサイトが主成分を占める状態です。NECTINやIL1B-IL1R2を中心とする炎症性シグナルが更新し、好中球やNK細胞、γδT細胞、MAIT細胞などを炎症部位へ動員するネットワークが構築されています。
- E5（ストレス関連・免疫不全型）: ストレス応答を示すT細胞やNK細胞、および好中球が豊富な免疫不全状態です。IFN-II、TNF、TGF-βといった強固な炎症性・微小環境ストレスシグナルに支配されています。SCENIC解析により、NF-κBファミリー（BCL3, NFKB1, RELB）、AP-1（FOS, JUN）、ERストレス応答（ATF4, DDIT3）、および酸化・代謝ストレス応答（NFE2L2, FOXO1/3）といった転写因子群が共通して活性化していることが確認されました。
今回NMFにより同定された5つのエコタイプは、病期（Stage）という従来の分類枠組みを超えた、組織化された細胞構造のバリエーションを捉えており、細胞間コミュニケーション、サイトカインシグナル伝達、転写因子プログラム、共有遺伝子モジュールと密接に関連しています。
解析の結果、同一のエコタイプに分類される患者サンプルは、臨床診断名（病期）が異なっていても酷似した細胞構成およびシグナル特性を示すことが示されました。
病勢進行に伴う免疫機能の低下とストレス応答プログラムの獲得
骨髄微小環境における細胞群の再構成と機能変化を解明するため、本研究では「疾患進行に伴う微小環境全体の機能再構成およびストレス応答の分子機構」という課題に対し、単一細胞レベルでの網羅的解析に基づく検証を試みました。
リンパ球系から骨髄系へのシフトと二相性動態
疾患がMGUS/SMMからNDMM、RRMMへと進行するにつれて、CD4+ T細胞、NK細胞、γδT細胞などのリンパ球系細胞の割合が著しく減少する一方、CD14+モノサイトなどの骨髄系細胞が増加する顕著な集団シフトが確認されました。成熟B細胞は一貫して低下し、CD8+ T細胞の全体的な存在比率は比較的安定して推移します。 一方で、樹状細胞（DC）および好中球は、前駆病態から進行期にかけて単純に増減するのではなく、二相性（Biphasic）の動態を示すことが明らかとなりました。また、CD56+ 未熟NK細胞（NK_GZMK_SELL）はMGUSからNDMMにかけて減少した後にRRMMで再び上昇するという、V字回復的な推移を見せました。さらに、腫瘍細胞自体において新規のNK細胞チェックポイント分子であるIGSF8の発現が亢進していることも確認されました。
骨髄系細胞による抗原提示能低下と免疫抑制環境の構築
モノサイトや好中球においては、疾患進行に伴いCD14+モノサイトにおけるMHC class-II分子（HLA-DRA、CD74など）の発現が段階的に減退し、抗原提示能の低下が明瞭となりました。また、骨髄由来免疫抑制細胞（MDSC）様の遺伝子プログラムや、免疫チェックポイント分子（CSF1R、PILRA、VSIR、LILRB2、TIM-3/HAVCR2など）の高発現が観察され、多角的な免疫抑制機能が構築されていることが示されました。
ストレス関連T/NK細胞の出現と代謝不全の解明
進行期（NDMM/RRMM）で増加する「ストレス関連T細胞/NK細胞」の解析では、FOS/JUNをはじめとするAP-1転写因子群およびNF-κBシグナルの増強、小胞体（ER）ストレスや低酸素応答の亢進が認められました。Monocle2を用いたPseudotime軌跡解析により、これらの細胞は従来の疲弊（Exhaustion）状態や熱ショックタンパク（HSP）応答とは異なる独自の分化経路を辿っていることが示されました。
さらに、Spectra解析を用いた無方向性遺伝子モジュール（GM）の解析により、以下の重要な統合的知見が得られました。
- GM1（IFN-I応答）: MX1, ISG15, IFI44L, IFI6などを特徴とし、型Iインターフェロン応答を反映します。
- GM3（酸化的リン酸化・ATP合成）: 正常なミトコンドリアエネルギー代謝を担うモジュールです。
- GM5（ストレス/AP-1/ケモカイン）: NR4A2, DUSP2, JUN, TNFAIP3を含み、微小環境ストレスに応答します。
多変量解析の結果、GM5（ストレスモジュール）はGM3（酸化的リン酸化モジュール）と有意な負の相関（r < 0）を示し、GM1（炎症モジュール）とは相関を示さないことが明らかとなりました。これは、微小環境におけるストレス応答プログラムが慢性炎症の副産物ではなく、エネルギー代謝の破綻（酸化的リン酸化およびATP合成機能の低下）と直接的に結びついた機能不全状態であることを示しています。
こうした細胞状態および代謝プログラムの変容は、病期のみによって一義的に決まるものではなく、腫瘍自体の量的な広がりや過去の治療内容によっても修飾されます。
腫瘍負荷および治療曝露が及ぼす動的なTMEへの影響
骨髄微小環境の改変は、腫瘍自体の量的な広がり（腫瘍負荷：Tumor Burden）および患者が受けた過去の治療薬剤という、2つの動的な要素によって強く修飾されます。
腫瘍負荷（Tumor Burden）と染色体異常の影響
RRMMを除くコホートでの解析により、腫瘍負荷の増大はCD4+ T細胞、NK細胞、γδT細胞の減少と相関し、前駆細胞富裕型であるE1シグナルの増加と正の相関を示しました。一方で、腫瘍細胞の増殖能や遺伝的不均一性といった質的な特徴は、免疫構成全体に対する直接的な影響が限定的でした。
染色体異常との関連においては、17p欠損（del(17p)）を有する患者において、CD4+ T細胞の減少およびCD14+/CD16+ モノサイトの増加とともに、E4（炎症性モノサイト型）シグナルが有意に高値を示す（特にSMMおよびRRMMにおいて）ことが明らかとなりました。
既往治療が及ぼす長期的・短期的影響
既往治療、特に再発・難治期における各種薬剤の曝露は、微小環境に異なる時間軸での変化をもたらします。
- 抗CD38抗体薬: 治療終了直後はNK細胞やCD8+ T細胞に一時的な変動が生じますが、時間経過とともに回復します。一方、治療終了から6ヶ月以上経過したサンプルでは、CD14+モノサイトの増加ならびにCD4+ T細胞や成熟B細胞の減少という持続的な免疫改変が確認されました。
- レナリドミド: MAIT細胞の割合は、レナリドミド既治療・抵抗性サンプルと比較して、レナリドミド未治療（Naïve）または感受性（Sensitive）サンプルにおいて有意に低値であることが示されました。
- 造血幹細胞移植（ASCT）およびCAR-T療法: ASCT既往やレナリドミド抵抗性患者ではモノサイトやリンパ球サブセットに変化が見られる一方、過去のCAR-T療法曝露による長期的な微小環境への影響は比較的小規模であり、CD14+モノサイトの軽度な増加にとどまりました。
臨床的アウトカムおよび免疫療法応答性の予測価値
TMEのエコタイプ分類は、従来の臨床病期（ISSなど）を補完し、患者の予後や高度免疫療法に対する治療応答性を予測するうえで重要な役割を持っています。
NDMMコホート（MMRF Immune Atlas）における検証と機能的二面性
MMRF Immune Atlasコホート（初発多発性骨髄腫263例）を用いた検証において、エコタイプは無増悪生存期間（PFS）および全生存期間（OS）の独立した予後予測因子であることが示されました（多変量Cox解析にてPFS: P = 1.40 × 10⁻⁷）。
- E1（前駆細胞富裕型）および E3（エフェクターT細胞優位型）: PFSおよびOSが有意に不良であり、急速進行（Rapid Progressor）群に多く観察されました。E3が高レベルのエフェクターT細胞を集積しているにもかかわらず予後不良を示す背景には、E3を呈する患者群が有意に高齢であり、T細胞の加齢変化・細胞老化（Senescence）や機能的消耗が進行していることが関与しています。
- E2（恒常的免疫監視型）および E5（ストレス関連・免疫不全型）: PFSおよびOSが良好であり、非進行（Non-Progressor）群に富裕していました。初発（NDMM）における標準的な化学療法下においてE5が良好な生存を示すのは、このストレスプログラムが深刻な細胞死や即時の機能脱失ではなく、微小環境のストレスに対する「適応的応答」を反映しているためと考察されます。
免疫療法（CAR-T / TCE）における不応答性と機能的フィットネス
一方で、CAR-T療法やT細胞誘導抗体（TCE）療法が直後に投与された患者群の解析では、治療不応答者および短期奏効者（奏効期間6ヶ月未満）の64.3%がE5（ストレス関連・免疫不全型）に分類されました。これは、標準治療下での適応的応答とは対照的に、T細胞の再導向（Redirection）を主軸とする高度な免疫療法においては、E5に認められるミトコンドリア代謝障害やストレス応答に伴うT細胞の「機能的整合性（フィットネス）」の低下が、治療抵抗性を引き起こす主要因となることを示しています。
他コホートでの検証と臨床解釈上の注意点
本研究の妥当性を検証するため、遺伝子発現予測モデルを用いて外部コホートであるUAMSコホート（GSE136400、N=434、バルクマイクロアレイデータ）への適用を行いました。 その結果、UAMSコホートにおいてはE1/E3が比較的良好なPFSを示し、E2が不良なPFSを示すなど、MMRFコホートとは一部相違する傾向が観察されました。この乖離の背景には、解析技術の違い（バルク組織マイクロアレイ vs 単一細胞RNA-seq）や投与された治療レジメンの異質性が存在するため、エコタイプと生存期間の関連を解釈する際には、解析プラットフォームや治療文脈に応じた慎重な評価が必要となります。
さいごに
本研究は、多発性骨髄腫の病態スペクトラム全体にわたる大規模単一細胞マルチオミクス解析により、骨髄微小環境の体系的なアトラスを構築し、5つの再現可能な「エコタイプ」概念を確立しました。この分類は、疾患の進行、治療曝露、そして患者の予後や免疫療法への応答性を統合的に理解するための基盤となると考えられます。
実臨床への応用に向けては、マルチパラメーターフローサイトメトリーや多重免疫染色などの既存技術を用いて、主要な免疫浸潤指標やT細胞の機能状態（ストレス・老化マーカー）、骨髄系細胞の分布を計測し、エコタイプを簡略的に同定する判定体系の開発が期待されます。一方で、本研究には縦断的（経時的）追跡サンプルの不足や、再発・難治期における追跡期間の制限といった限界が存在するため、今後は前向きコホートによるさらなる検証が望まれます。
