Method Article
Cytofastは、クラスタリングからの出力を分析するために使用される視覚化ツールです。Cytofastを使用して、FlowSOMとCytosploreの 2 つのクラスタリング方法を比較できます。Cytofastは、質量サイトメトリーデータの定量的および定性的な概要を迅速に生成し、異なるクラスタリングアルゴリズム間の主な違いを強調することができます。
大量のサイトメトリーによって生成されるデータの複雑さは、分析結果を迅速に可視化するための新しいツールを必要としています。細胞球団やFlowSOMなどのクラスタリング方法は、細胞クラスターの可視化と同定に使用されます。下流分析のために、新しく開発されたRパッケージCytofastは、クラスタリング法から結果を迅速に可視化することができる。Cytofastは、細胞クラスターの表色特性を考慮に入れ、細胞クラスターの存在量を計算し、次いでグループを定量的に比較する。このプロトコルは、腫瘍の課題に続いて腫瘍の課題に対する免疫系の変調(すなわち、ナチュラルキラー[NK]細胞応答)に基づく質量サイトメトリーデータの使用に対する細胞化剤の適用を説明する(PD-L1遮断)。FlowSOMおよびサイトスプロアによるサイトファストの有用性の実証を示す。細胞化は、グループ関連の免疫細胞クラスターの視覚的表現と免疫系組成との相関を迅速に生成します。クラスタリング分析では違いが見られますが、グループ間の分離は両方のクラスタリング方法で表示されます。Cytofastは、活性化NK細胞サブセットのより高い存在量を含むPD-L1処理によって誘導されるパターンを視覚的に示し、活性化マーカーのより高い強度(すなわち、CD54またはCD11c)を発現する。
質量サイトメトリー(飛行時間によるサイトメトリー、またはCyTOF)は、何百万もの単一細胞における幅広い細胞内または細胞外バイオマーカーの検出を可能にします。質量サイトメトリーデータの高次元性は、SPADE1、FlowMaps2、FlowSOM 3、フェノグラフ4、VorteX 5、および足場マップ6のような細胞クラスタリング技術などの特定の解析ツールを必要とする。また、様々な次元低減ベースの技術が開発されている(すなわち、主成分解析[PCA]7、t分散確率的近傍埋め込み[t-SNE]8、階層確率論的隣人埋め込み[HSNE]9、均一多様体近似図10、及び拡散マップ11)を開発し、高速化、解釈及び可視化を向上させる。
高次元流量と質量細胞測定データの下流分析では、クラスター周波数と臨床結果とのリンクに関する統計的検査を実行する自動プロセスが欠けていることがよくあります。以前は、CytosploreまたはFlowSOMによるクラスタリング技術の視覚的および定量的な下流分析を可能にする、Cytofast 12と呼ばれるRベースのワークフローを開発しました。
ここで説明するプロトコルは、RでのCytofastの使用を明確にし、定量的および定性的なヒートマップとグラフを生成する方法を示しています。さらに、観察された免疫表現型と臨床結果との間の接続の決定を促進する。このレポートでは、FlowSOMとCytosploreという 2 つの異なるクラスタリング手順を使用した特定の質量サイトメトリー データセットの分析についても説明します。両方のクラスタリング方法でCytofastを用いることで、NK細胞の活性化表現型がPD-L1免疫チェックポイント遮断の影響を受けることを示す。
すべての動物実験はLUMC動物実験委員会によって承認され、オランダおよび欧州委員会のガイドラインに従ってLUMCの動物実験ガイドラインに従って実行された。
注:実験的セットアップのために、C57BL/6マウスを、リン酸緩衝生理食塩水(PBS)の濃度0.3 x106細胞/200μLのマウス大腸腫瘍MC38で右脇腹に皮下接種した。10日後、腫瘍が触知可能な場合、マウスをPD-L1遮断抗体(クローンMIH-5、200μg/マウス、腹腔内注射)で治療するか、または模擬処理した。腫瘍をPD-L1注射の3日後に切除し、ex vivoを処理し、38マーカー13を用いてCyTOF質量サイトメトリーにより分析した。
1. データ分析のための機器とソフトウェア
メモ: 2.4 GHz 以上のコンピュータ (Windows 7 以降) とプロセッサ I5、インストールされているメモリ RAM 6 GB、および 10 GB のハード ドライブの空き容量を使用します。RパッケージCytofastは、主にflowCore、pheatmap、およびggplotの既存の機能を使用しています。 R で実行されるコマンド ラインは、プロトコルに含まれます。R 命令のリソースはhttps://education.rstudio.com/にあります。
2. クラスターの作成
メモ:サイトファストを用いてサイトスプロアとFlowSOMの2つのクラスタリング方法を紹介するために、PD-L1治療の3日後に腫瘍マイクロ環境におけるNK細胞(CD161+)を分析する。
3. 可視化:後処理クラスタリング分析
メモ: この手順は、両方のクラスタリング方法に共通する方法です。したがって、FlowSOMまたはサイトスプロアを使用してクラスタリングした後に実行できます。
ワークフロー Cytofast(図1)は、分析ソフトウェア(FlowSOMまたはCytosplore)によって最初にクラスタ化されたデータの定量的かつ定性的な概要を提供することを意図しています。Cytofastは、解析で識別されたすべてのクラスターのヒートマップを含む、いくつかの可能な出力を実行し、マーカー式に基づいています (図 2および図 3)。上部の樹形図は、識別されたクラスター間の階層的な類似性を表します。上部パネルには、各サンプル内の対応するサブセットの相対量を示す別のヒートマップが表示されます。右側の樹形図は、サンプル間の類似性を示し、サンプル間のユークリッド距離で実行される階層クラスタリングに基づいています。組み合わせたヒートマップはFlowSOMに続き、図2ではCytofast、続いて図3にCytosploreが続きます。 図 4および図 5に示すように、Cytofastを使用してデータを定量的に表示し、結果を箱ひげ図 (cytoBoxplots 関数を使用) で表示することもできます。
2つの異なる方法の間で同様のクラスターが見出され(例えば、Cytosploreからのクラスター8はFlowSOMからのクラスター10に対応)、PD-1およびLAG-3のようないくつかの阻害マーカーの共発現は両方の方法でまだ見られた。どちらのクラスタリング方法でも、PD-L1との間での判別が可能でした。PBSはマウスを治療した。対照的に、両方の方法の違いがいくつか強調表示されます。FlowSOMは 2 つのクラスタ (MHC-II+) を識別し、Cytosploreは 1 つのクラスタ (MHC-II+dim)のみを示します。これは、NK細胞がCD161+細胞上で手動でゲートされ、その後FlowSOMによってさらに処理された最初のゲーティング戦略によるものです。しかし、Cytosploreは最初のHSNEレベルでCD45+母集団から自動的に細胞をゲートし、その後、より高い階層レベルでクラスタリングされました。したがって、Cytosploreは、手動ゲーティングがCD161に焦点を当てた方法よりも正確にNK細胞サブセットを定義した。それにもかかわらず、サンプルの階層クラスタリングは、右側のデンドログラムに示すように保存され、2つのグループ間の分離(PD-L1とPBS)が選択されたクラスタリング方法に依存していないことを示す。
クラスターの数は、両方の方法を使用して手動で定義できます。Cytofastを使用すると、ユーザーはデータの異質性を評価し、データを分割するクラスターの数を選択する方法に関する洞察を提供できます。msiPlot 関数 (ステップ 3.4.2) などの他の機能は、グループごとの各マーカーの中央値信号強度 (MSI) プロットを示します (図 6および図7)。この機能は、PD-L1処理群のNK細胞におけるCD54またはCD11cの発現の増加などの世界的な変化を検出することを可能にする。オプション機能は、棒グラフやその他のデータ表現方法でのデータの表示など、Cytofastパッケージに組み込むことができます。後者は、Rによって生成することができるggplotツールを追加する必要があります。
図1:サイトファストパッケージのワークフローこのデータは、免疫療法による治療の3日後に腫瘍から集団サイトメトリーによって生成されたか、または未治療のままにした。2 つの異なるクラスタリング技術を比較しました:サイトスプロールとFlowSOM.サイトファストは、2つの技術の違いを可視化するために使用された。この図の大きなバージョンを表示するには、ここをクリックしてください。
図 2: Cytosploreに続くCytofastによって分析されたグループごとのクラスターの概要とクラスターの存在量 。全NK細胞クラスターのヒートマップ(CD161+細胞は細胞スプロールによって自動的に定義される)、免疫療法の3日後に同定した(PD-L1)。示されたデータは、Cytosploreクラスタリングに基づいており、未処理およびPD-L1処理されたグループからプールされる。ArcSinh5 変換式マーカーのレベルは、虹のスケールで表示されます。下部パネルでは、各サンプルの相対的な存在量は緑色から紫色のスケールで表されます。右側の樹形図は、サブセット周波数に基づくサンプル間の類似性を表します。周波数スケールは平均の分散を表します。低周波または高周波は、それぞれ緑色または紫色で表されます。この図の大きなバージョンを表示するには、ここをクリックしてください。
図 3: Cytofastフォロー FlowSOMによって分析されたグループごとのクラスターの概要とクラスターの存在量免疫療法(PD-L1)の3日後に同定された全NK細胞クラスター(CD161+イベントにプリゲート)のヒートマップ。表示されるデータはFlowSOMクラスタリングに基づいており、未処理のグループと PD-L1 処理されたグループからプールされます。ArcSinh5 変換式マーカーのレベルは、虹のスケールで表示されます。下部パネルでは、各サンプルの相対的な存在量は緑色から紫色のスケールで表されます。右側の樹形図は、サブセット周波数に基づくサンプル間の類似性を表します。周波数スケールは平均の分散を表します。低周波または高周波は、それぞれ緑色または紫色で表されます。この図の大きなバージョンを表示するには、ここをクリックしてください。
図 4: Cytosploreによって定義されたクラスターの箱ひげ図を使用したサイトファスト表現 。各クラスターの頻度は箱ひげ図で表され、2つのグループ(PBSとPD-L1)に分けられます。1 つの個々のドットは、1 つのマウスに対応します。この図の大きなバージョンを表示するには、ここをクリックしてください。
図 5: FlowSOMで定義されたクラスターの箱ひげ図を使用したサイトファスト表現 。各クラスターの頻度は箱ひげ図で表され、2つのグループ(PBSとPD-L1)に分けられます。1 つの個々のドットは、1 つのマウスに対応します。この図の大きなバージョンを表示するには、ここをクリックしてください。
図6:Cytosploreによって自動的にゲートされたNK細胞からの信号強度プロットの分布。信号強度の分布は、CD45、CD11c、および CD54 の 3 つの特定のマーカーのヒストグラムに示されています。この図の大きなバージョンを表示するには、ここをクリックしてください。
図7:FlowSOMによって自動的にゲートされたNKセルからの信号強度プロットの分布。信号強度の分布は、CD45、CD11c、およびCD54の3つの特定のマーカーのヒストグラムに示され、グループPBSおよびPD-L1によって分離される。この図の大きなバージョンを表示するには、ここをクリックしてください。
補足ファイル 1.1–1.8.このファイルを表示するには、ここをクリックしてください (右クリックしてダウンロードしてください)。
補足ファイル 2.1– 2.10.このファイルを表示するには、ここをクリックしてください (右クリックしてダウンロードしてください)。
補足ファイル 3.このファイルを表示するには、ここをクリックしてください (右クリックしてダウンロードしてください)。
補足ファイル 4.1–4.8.このファイルを表示するには、ここをクリックしてください (右クリックしてダウンロードしてください)。
Cytofastは、治療固有の細胞サブセットを強調表示および定量することにより、細胞測定データの迅速かつグローバルな探索を提供する迅速な計算ツールです。説明したプロトコルは、CytosploreまたはFlowSOMを使用してクラスタリング分析をさらに処理することを目的としています。他のクラスタリング解析ツールはCytofastに適していますが、各セルをサブセットに割り当てるにはCytofastを使用する必要があります。ただし、Cytofastはクラスタリング方法ではないため、使用する前にクラスタリング手順が必要です。
ここで行った分析は、腫瘍微小環境における特定のCD161+NK細胞サブセットがPD-L1遮断に敏感であることを示した。これは、クラスタリング法としてCytosploreとFlowSOMの両方を使用して観察された表現型と豊富さの変化によって示されました。どちらの方法も、わずかに異なる周波数(Cytosploreの場合は15%-20%、FlowSOMの場合は30%~40%)でメインNK細胞クラスター(CD11b + NKG2A+)を区別しました。図2と図3の右側のパネルに表示される両方のデンドログラムが同様の結果を示したため、豊富さとこの近似の違いはグローバルパターンに影響を与えませんでした。Cytofastを使用することにより、NK細胞クラスター表現型および豊富さの分析に基づいてPD-L1処理および未処理マウスを分離することが可能である(選択されたクラスタリング方法とは無関係)。
記録されたパラメータに応じて、プロトコルを変更する必要があります。具体的には、クラスタリング分析の実行中に、時間や背景などの特定のパラメーターを削除する必要があります。さらに、各セルをサブセットに割り当てることが重要です。cfData 関数は、サンプルごとのクラスターごとの生のセル数を cfList に追加するだけです。このステップから、セクション 3 で説明したようにサイトヒートマップを構築できます。
Cytofastは、異なるクラスタリング方法13を比較するための可視化および定量ツールとして正常に使用されている。このRパッケージは、臨床変数を使用してクラスタのグループ間の関連をテストできるグローバルテスト14などの高度な機能とも互換性があります。将来的には、グローバルテストツールやその他のアルゴリズムをCytofastと統合して、より詳細な視覚化と定量化を実現できます。
著者たちは何も開示する必要はない。
我々は、提案番号675743(ISPIC)の下でH2020 MSCA賞の欧州委員会からの資金提供を認めます。私たちは、プロトコルをテストするためのテチェ・ファン・デル・スルイスとアイリス・パルディックに感謝します。
Name | Company | Catalog Number | Comments |
Computer | Dell | NA | NA |
このJoVE論文のテキスト又は図を再利用するための許可を申請します
許可を申請This article has been published
Video Coming Soon
Copyright © 2023 MyJoVE Corporation. All rights reserved