方法論記事

ヒト筋肉間脂肪組織からの核の単離と下流の単一核RNAシーケンシング

1.6K 回視聴

DOI:

10.3791/66784

2024年5月3日

この記事について

サマリー

筋肉間脂肪組織(IMAT)の生物学は、ヒト組織のアクセス性が限られているため、ほとんど未踏です。ここでは、このユニークな脂肪貯蔵庫の細胞組成を特定するための単一核RNAシーケンシングのための凍結ヒトIMATの核単離とライブラリ調製のための詳細なプロトコルを紹介します。

要約

筋肉間脂肪組織(IMAT)は、筋繊維の間に位置する比較的研究されていない脂肪貯蔵庫です。IMAT含有量は、年齢とBMIとともに増加し、代謝性疾患および筋肉変性疾患に関連しています。しかし、IMATの生物学的特性と周囲の筋線維との相互作用についての理解は著しく不足しています。近年、シングルセルおよび核RNAシーケンシングにより、いくつかのヒト組織の細胞タイプ特異的なアトラスが得られました。しかし、ヒトIMATの細胞組成は、ヒトの生検収集からのアクセス可能性に固有の課題のために、ほとんど調査されていないままです。ヒトIMATは、採取される組織の量が限られていることに加え、骨格筋組織や筋膜に近いため、処理が複雑になります。脂肪細胞の脂質を多く含む性質により、脂肪細胞は単一細胞の単離と互換性がありません。したがって、単一核RNAシーケンシングは、単一細胞分解能で高次元トランスクリプトミクスを得るのに最適であり、IMATの正確な細胞組成を含むこのデポの生物学を明らかにする可能性を提供します。ここでは、単一核RNAシーケンシングのための凍結ヒトIMATの核単離とライブラリ調製のための詳細なプロトコルを紹介します。このプロトコルでは、液滴ベースのアプローチを使用して数千の核をプロファイリングできるため、希少で存在量が少ない細胞タイプを検出する能力が得られます。

概要

筋肉間脂肪組織(IMAT)は、筋線維1の間および周囲に存在する異所性脂肪貯蔵庫です。Goodpasterらによる最近のレビューで詳細に説明されているように、IMATは、高解像度のコンピュータ断層撮影法(CT)および磁気共鳴画像法(MRI)(図1AB)を使用して検出でき、全身の筋線維の周囲および内部に見られます1。IMATの量は個人によって大きく異なり、BMI、年齢、性別、人種、座りがちな状況に影響されます2,3,4。さらに、IMAT沈着は、筋肉の変性5に関連する病的状態によく見られ、多くの研究で、肥満、2型糖尿病、メタボリックシンドローム、およびインスリン抵抗性6,7,8,9の個人におけるIMAT量の増加が記録されています。それにもかかわらず、IMATの細胞的および生物学的特性は解明され始めたばかりです。アクセス性が限られていること、IMATの位置と内容物が全身に異なることが、このユニークな脂肪貯蔵所2からのサンプル収集を困難にしています。さらに、サンプルは収集時に骨格筋(SM)で容易に「汚染」されるため、異なる組織からの生物学的寄与の分離を解読するのが困難になります(図1C)。そのため、過去10年間で大きな注目を集めている一核RNAシーケンシング(snRNA-seq)は、IMATおよびSM由来の遺伝子発現パターンをシングルセルの分解能で分離するための理想的な方法論として機能します。さらに、核の単離は、脂質を多く含む脂肪細胞が大きいため、脂肪組織に特に有用であり、細胞の完全性を損なうことなく単一細胞懸濁液に解離することは不可能です。最後に、この技術は、IMAT特異的脂肪細胞の新規マーカーを発見し、異なる前駆細胞集団の組成と存在を明らかにするだけでなく、病理学的および正常な状態での細胞組成の変動を研究する可能性を秘めています。

figure-introduction-1
図1:IMATの画像。 (A)中年の痩せた女性と(B)肥満の中年男性からのIMATの代表的な磁気共鳴(MRI)画像。赤:皮下脂肪組織、黄:筋肉間脂肪組織、緑:骨格筋、青:骨。画像提供:Heather Cornnell氏、AdventHealth Translational Research Institute。(C)IMATを使用した新鮮な組織サンプル(黒の破線で囲まれています)。画像提供:AdventHealth Translational Research Institute の Meghan Hopf 氏とコロラド大学の Bryan Bergman 氏。この図は、Goodpaster et al.1 の許可を得て変更されています。 この図の拡大版を表示するには、ここをクリックしてください。

豚、鶏、牛の肉の霜降り(特にIMAT)をシングルセル(sc)とsnRNA-seq10を用いて調査した畜産業界から多くの研究が発表されています。これらの研究により、脂肪細胞のいくつかの亜集団とIMATの潜在的な前駆細胞のマーカーが特定されました11,12,13;しかし、これらの細胞組成物がヒトIMATに翻訳されるかどうかは不明です。私たちの知る限り、変形性股関節症の男性患者から得られた、脂肪浸潤を伴うヒト筋肉の細胞の不均一性をsnRNA-seq14を使用して調べた研究は1つだけです。研究者らは、小さな脂肪細胞集団といくつかの線維脂肪形成前駆細胞(FAP)亜集団を報告しました 筋核14の大規模な集団内。私たちの研究は、snRNA-seqを使用して、ヒトの筋肉から手動で解剖したIMATの細胞組成を直接調べる方法を開発した最初の研究です。

重要なことは、snRNA-seqのプロトコールは、利用可能な組織の量と特定の組織の物理的特性によって最適な処理ステップが決まるため、研究対象の特定の組織に合わせてカスタマイズする必要があるということです。IMATの組織収量は通常小さく、超音波ガイド下生検を行う場合でも、50 mgを超えないことがよくあります。したがって、この希少な組織の適切な処理が不可欠です。このプロトコルは、ヒトIMATを研究する研究者にとって貴重なリソースとなると考えています。

プロトコル

このプロトコルに使用されたサンプルは、Western IRB-Copernicus Group(WCG)の治験審査委員会によって承認されヘルシンキ宣言に従って実施されたStudy of Muscle Mobility, and Aging(SOMMA)15の一部でした。参加者は、研究への参加について書面によるインフォームドコンセントを提供しました。

注:このプロトコルは、ナノウェルベースのプラットフォーム16上で100mgのヒト腹部皮下脂肪組織を使用した以前のプロトコルから適応されています。現在のプロトコールは、50 mgのヒトIMATおよび液滴ベースのプラットフォームを使用したライブラリ調製に最適化されています。非ヒトIMATまたは他の脂肪デポーからの核単離のためのこのプロトコルのさらなる最適化が必要となる場合があります。

1. バッファーと試薬の調製(表1、表2)

注:実験当日にバッファーを新たに調製し、再利用しないでください。

  1. 遠心分離機を4°Cに予冷します。
  2. 均質化バッファーと核分離培地を調製します。
    1. 氷のバケツ2つを用意し、2 x 15 mLコニカルチューブを予冷します。
    2. ホモジナイズバッファー(HB)用のすべての試薬を、 表1に記載されている順序で15 mLのコニカルチューブに混合します。氷の上を保ってください。ボルテックスで混ぜます。
    3. 核単離培地(NIM)用のすべての試薬を、 表2の順序付きリストにある15 mLのコニカルチューブに混合します。氷の上を保ってください。ボルテックスで混ぜます。
    4. 100 μL の Triton X-100 を 900 μL のヌクレアーゼフリー水に加えて、10% Triton-X を調製します。適切な混合を確保するための渦巻き。室温(RT)で保管してください。

試薬容量 (μL)最終濃度(mM)
1倍2倍
1 M MgCl210205
1 M トリスバッファー、pH 8.0204010
2 M KCl255025
1.5 Mショ糖(-4°C)334668250
1mM DTT240.001 (~1 μM)
100xプロテアーゼ阻害剤20401倍
スーパーアシン 20 U/μL40800.4 U/μL
ヌクレアーゼフリー水15493098-
総ボリューム20004000-

表1:均質化バッファー(HB)。 氷の上を保ってください。ボルテックスで混ぜます。

試薬容量 (μL)最終濃度(mM)
1倍2倍
EDTAの0.40.80.1
リボロックRNAse阻害剤(40U/μL)40800.8 U/μL
1%BSA-PBS(-/-)1959.63919.2-
総ボリューム20004000-

表2:核分離培地(NIM)。 氷の上を保ってください。ボルテックスで混ぜます。

2. 凍結組織の粉砕(図2A)

  1. 均質化のためにワークステーションをセットアップします。
    1. キャニスターに液体窒素(LN2)を入れます。
      注意: LN 2で作業するときは 常にゴーグルとクライオグローブを着用してください。
    2. 自動ダンサー用に、2x乳鉢、1x乳棒、1xマイクロスコップヘラ、1xガラスドンス、および1xステンレス鋼乳棒を入手します。
    3. 自動ダンサーを設定します。
    4. ビーカーに氷を入れ、グラスダンスを予冷します。
  2. 2つの乳鉢(乳棒とへらを含む)にLN2 を充填して、器具を冷やします。LN2 を蒸発させて繰り返します。
  3. 機器が冷えている間に、グラスドンスにHB1mLを追加します。.
  4. 最後にもう一度両方のモルタルにLN2 を充填し、50 mgのIMATサンプルをモルタルの1つに注ぎます。
  5. 乳棒を使用してIMATをティッシュピースにそっと押し付けて粉砕し、細かく砕きます。すべてのピースが粉砕されていることを確認してください。
  6. LN2は、組織を粉砕しながらゆっくりと蒸発します。組織が適切に粉砕され、LN2のモルタルの1/4〜1/2がまだ残っている場合は、モルタルをモルタルのリップに向かって傾けて、粉砕された組織をリップで収集します。LN2を完全に蒸発させます。
  7. 最後のLN2 が蒸発した直後に、粉砕した組織を1mLのHBが入ったガラス製菓にすくい取ります。

3. 粉砕組織の均質化

  1. 自動ダウンサーを使用して粉砕された組織を均質化します。ガラスをステンレス製の乳棒を上下に動かし、順方向に10ストローク、続いて逆方向に10ストロークします。
  2. 均質化後の溶液が濁り、目に見える組織片が含まれていないことを確認してください。淡いピンク色は、筋肉組織による汚染のためにしばしば予想されます。
  3. ホモジネートを氷上で予冷した1.7 mL低結合チューブに移します。
  4. 400μLのHBを使用してダウンスをすすぎ、すべての材料が転写されていることを確認し、チューブに追加します。
    注:一度に2つのサンプルを処理できます。これを行うには、HB と NIM の量を 2 倍にします。1つの組織サンプルを粉砕して均質化し、その直後に2番目の組織サンプルを粉砕して均質化することで、単離とクリーンアップのステップを並行して実行できるようにします。

4. 原子核の単離とクリーンアップ(図2B)

  1. ホモジネートに14 μLのTriton-X(10%)を加えて、0.1%濃度にします。
  2. チューブを氷の上と暗闇に10〜15分間置き、3分ごとにボルテックスします。
  3. 100 μmセルストレーナー1個と40 μmセルストレーナー1個(サンプルあたり)を、それぞれ100 μLのRT DPBSで50 mLコニカルチューブにプレウェットします。
  4. 100 μmのセルストレーナーでホモジネートをろ過します。
  5. 1.7 mLチューブを400 μLのHBですすぎ、100 μmセルストレーナーでろ過します。
  6. 次に、40 μmセルストレーナーで溶液をろ過します。
  7. 各チューブ内の~900 μLに対応する2本の予冷済み1.7 mL低結合チューブに同量の溶液を移します。
  8. チューブを2700 x g 、4°Cで10分間遠心分離します。 遠心分離後に小さなペレットが見える必要があります。
  9. 上部の脂質層と残りの上清を取り出して廃棄し、最初のチューブから~50 μLの溶液を残します。
  10. 2番目のチューブについても繰り返します。
  11. 20倍に静かにピペッティングしてペレットを最初のチューブに完全に再懸濁し、新しい1.7 mL低結合チューブに移します。バブルの作成は避けてください。
  12. 2番目のチューブについてこの手順を繰り返し、再懸濁溶液を同じチューブに移します。
  13. NIM500μLを加え、ピペッティングで混合します。
  14. バランスを1000 x g にしてチューブを4°Cで10分間遠心分離します。
  15. 上清を取り出し、~50 μL を残し、ペレットが再懸濁されるまでゆっくりと上下にピペットで動かします。オプションで、チューブの側面に脂質が残っている場合は、再懸濁したペレットを新しいきれいなチューブに移します。
  16. NIM200μLを加え、ピペッティングで混合します。

5. 核の染色とカウント(図2C および 図3)

注:カウントを容易にするために、明視野とDAPIチャネルの調整がカウントに大きく影響する可能性があるため、自動セルカウンターで「核カウント」プロトコルを設定します。破片ではなく原子核のみが捕捉されるようにチャネルを調整します。明視野チャンネルが、DAPI染色も持っている「オブジェクト」のみをマークするようにしてください。

  1. 生細胞染色液を1滴加え、氷の上、暗所に15分間放置します。
  2. 30 μmセルストレーナーで溶液をろ過します。
  3. ピペッティングで核溶液を混合し、10 μLの溶液をセルカウンティングチャンバースライドに加えます。
  4. 自動セルカウンターを使用して核をカウントします。
    注:最適濃度は、1.0 x 106/mLに対応する1000 nuclei/μLです。
    1. 核の塊が存在しないことを確認すると、チップを詰まらせて単一核の液滴を生成する可能性があります(図3)。
    2. 核濃度が十分に高くない場合は、溶液を1000 x g で4°Cで10分間スピンダウンしてペレットを取得し、上清を除去し、少量で再懸濁します。
    3. 溶液中の破片の程度が高い場合は、核溶液をより大容量のNIM(すなわち、1 mL)に再懸濁し、30 μmの細胞ストレーナーで再度ろ過します。次に、1000 x g で4°Cで10分間スピンダウンし、核の濃度に対して適切な量で再懸濁します。
  5. 核濃度を取得したら、ライブラリ調製の最初のステップに直接進みます。

figure-protocol-1
図3:単離された核の染色。 NucBlue/DAPIで染色した核のセルカウンターからの画像(左の画像)と対応する明視野の画像(右の画像)。微量の破片の存在は、明視野画像で明らかです。ここで使用する自動セルカウンターには、スケールバーを含めるオプションはありません。 この図の拡大版を表示するには、ここをクリックしてください。

6. ライブラリ調製およびシーケンシングパラメータ

  1. プロバイダーのWebページ17で入手可能な液滴ベースの単一核アプローチを使用したライブラリ調製のための詳細なプロトコルを参照してください。
    1. 目標核回収率10,000を目標にしてください。ただし、高レベルの破片や壊れやすい核を持つサンプルの場合、回収される核の数は少なくなると予想されます。
    2. ステップ2.3の後、サンプルを4°Cで最大72時間保存します。ライブラリ調製プロトコルでは、より多くのサンプルの処理を並行して組み合わせます。これを行うには、ステップ 2.3 までの 2 つのサンプルを 2 日間連続して処理し、3 日目にはステップ 3 の 4 つのサンプルを一緒に処理し、ライブラリ調製プロトコルで進めます。
  2. シーケンシングパラメータ:核あたり50,000ペアエンドリードを目指すシーケンシングプラットフォーム上でのシーケンシング。
    注:このプロトコルで示されたデータは、NovaSeq 6000プラットフォームでシーケンシングされ、核あたり50,000のペアエンドリードを目指しました。

figure-protocol-2
図2:プロトコルワークフロー。 プロトコルの(A)ステップ2および3、(B)ステップ4、および(C)ステップ5のワークフローの概略図。フィギュアは BioRender.com で作成しました。 この図の拡大版を表示するには、ここをクリックしてください。

7. データ処理と分析

注:このプロトコルでは、結果として得られるシーケンシングデータの処理に使用される推奨ソフトウェアとRパッケージの一部が、初期前処理後のステップに焦点を当てて簡単に紹介されています(表3)。この研究では、一般的な品質管理(QC)メトリクスと、 図4の均一多様体近似および投影(UMAP)の例を提供します。ただし、バイオインフォマティクス分析の詳細な説明は、このプロトコルの範囲外です。したがって、読者はHeumosらによるシングルセル解析のベストプラクティスに関する最近のレビューを参照することができます18

  1. シーケンシングデータの前処理
    1. 読み取った単一核をヒト参照ゲノムGRCh38にマッピングします。
    2. イントロンリードをカウントに含めます。
  2. Seurat Rパッケージ19を使用してデータのQCとフィルタリングを実行します。
    1. 細胞複雑性スコアは、検出された遺伝子のlog(10)数を検出されたlog(10)リード数で割ることにより、細胞複雑性スコアを計算します。
    2. ヒストグラムまたはバイオリンプロットを使用して、核ごとに検出された遺伝子の数、ミトコンドリアの読み取り率、細胞の複雑さのスコアなど、最も重要なQC指標をプロットします。
    3. 核あたりの遺伝子数が200未満または10,000個以上、ミトコンドリアの読み取りが10%を超え、複雑さのスコアが0.8未満の核をフィルタリングします。
  3. データを正規化し、次元削減を実行します。
    1. Seurat SCTransform 関数を使用して、2000 個の変数機能を使用してデータを正規化します。
    2. Seurat R パッケージの RunPCAFindNeighborsFindClustersRunUMAP の関数を使用してデータをクラスター化します。
    3. UMAP をプロットして、データのクラスタリングを可視化します。
  4. DoubletFinder R パッケージ20 を使用して予測されたダブレットを除外し、データを再クラスタリングします。
  5. 組織内に存在すると予想される細胞タイプの既知の遺伝子マーカーを使用してクラスターにアノテーションを付けるか(教師ありアプローチ)、またはクラスター間で発現差のある上位5つの遺伝子に基づいて(教師なしアプローチ)。
  6. decontX21を使用して、アンビエントRNAの汚染度を測定し、アンビエントRNAの遺伝子発現マトリックスを調整します。
    1. 生の遺伝子マトリックスをバックグラウンドとして含めます。
  7. Seurat オブジェクトを保存して、後でデータを探索することができます。
    注: QC およびクラスタリング分析のコードは、 補足ファイル 1 にあります。
データワークフローで使用されるソフトウェア/Rパッケージ代替ソフトウェア/パッケージ処理ステップ
セルレンジャースターソロ、カリストトリミング、整列、マッピング
スーラSingleCellExperiment、セルレンジャーQC、分析、データ探索
ダブレットファインダーscds、scdblFinder、スクラブレットダブレット検出
デコントXSoupX、セルベンダーアンビエントRNA調整

表3:データワークフロー用のソフトウェア/ツール。

結果

このワークフローは、凍結ヒトIMATサンプルの処理をガイドして、単一核の分解能で遺伝子発現プロファイルを取得し、細胞タイプの同定を可能にするように設計されています。ここでは、SOMMA試験の参加者から代表的なIMATサンプルを1つ紹介します。

snRNA-seqデータの解析の最初のステップは、データの品質を評価して、データセットから削除する必要がある可能性のある低品質の核を特定することです。重要なことに、フィルタリングステップと閾値は、手元にあるサンプルとデータセットの特定のタイプについて決定されるべきであり、一般的に評価される指標は組織や細胞タイプによって異なる可能性があるためです22,23図4Aは、生成されたsnRNA-seqデータの品質を評価するために使用されるいくつかの主要な指標を視覚的に示しています。核ごとに検出される遺伝子の数は、シーケンシングの深さと細胞の種類に依存しますが、良質な核18,23では200を超えると予想されます。このプロトコルを使用して生成されたデータは、合計 4662 個の核から、核あたり 1134 個の遺伝子の中央値で、予想される範囲内にあることがわかりました。

ミトコンドリアリードの割合は、損傷した核や核に付着した周囲のRNAから高度のミトコンドリア汚染が発生する可能性があり、核の品質が悪いことを示しているため、評価されます。ここで紹介したデータセットでは、ミトコンドリアの読み取り率の中央値が2.65であり、これは文献24,25,26で一般的に使用される5%-20%の閾値をはるかに下回っています。リボソームリードの割合は、細胞の種類や組織によって異なります。ただし、リボソーム遺伝子の大部分がデータのクラスタリングに影響を与える可能性があるため、クラスタリングする前にリボソーム読み取り率を確認し、リボソーム遺伝子またはリボソーム遺伝子のレベルが高い核をデータセットから削除することをお勧めします。このプロトコルで生成されたデータは、中央値が2.46%、最大が16.5%とリボソームのリードが低かったため、この指標に基づいてフィルタリングしませんでした。最後に、細胞複雑性スコアは、検出された遺伝子のlog(10)数を検出されたリードのlog(10)数で割った値に基づいて計算されました。良質な核は0.8以上であると予想され、この研究で使用されたサンプルでは中央値0.92が得られました。これらのQC指標に基づいて、データセットからどの核を除外するかを決定できます。解析では、核あたりの遺伝子数が200未満または10,000個以上、ミトコンドリアの読み取りが10%を超え、複雑さのスコアが0.8未満の核をフィルタリングすることを選択しました。

最初の品質評価とフィルタリングのステップに続いて、UMAP を生成して核のクラスタリングを視覚化できます。クラスタリングは、SCT形質転換を使用して、最も変動する2000の遺伝子に基づいて実行されました。最初のクラスタリングステップは、QC特徴のいずれか(例えば、ミトコンドリアリードが高い核)がクラスター化しているかどうかを確認するために使用できます。さらに、このプロトコルで使用されていたDoubletFinder20を含む一部のダブレット検出方法では、クラスタリング情報が必要です。ダブレットファインダーは、ドロップレットベースのプラットフォームのプロバイダーが提案しているように、予想されるマルチプレットレートを4.8%に設定して使用されました。ダブレット除去後、周囲のRNA汚染レベルが推定されましたが、これは特に単一核調製で一般的であり、RNAは細胞溶解時に細胞質から放出され、ゲルビーズインエマルジョン(GEM)に分注され、次のライブラリ調製ステップで増幅されます。そのため、アンビエントRNAコンタミネーションの固有の問題を修正するために、いくつかのツールが開発されています(表3を参照)。RパッケージのdecontX21を使用し生のバックグラウンドマトリックス(空の液滴のみを含む)を使用して遺伝子発現マトリックスを調整し、実際の遺伝子発現シグネチャーを強化しました。

クラスタリングと存在量が少ない細胞タイプを検出する能力は、核の数に依存します。この研究では、QCフィルタリング、ダブレット除去、およびアンビエントRNA調整後、合計3817個の核からIMAT(図4B)で予想されるすべての主要な細胞タイプを検出しました。これらには、幹細胞、線維脂肪形成前駆細胞(FAP)、成熟脂肪細胞、周皮細胞、平滑筋細胞、免疫細胞、筋前駆細胞、骨格筋細胞汚染による筋核細胞が含まれていました。

全体として、このプロトコルは、IMATの生物学と細胞起源を解明するために重要な細胞型アノテーションの検出を可能にする高解像度の単一核データを生成することを実証しました。

figure-results-1
(A) サンプルの評価とシーケンシング性能の評価に不可欠な指標のバイオリンプロット。これには、核ごとに検出された遺伝子の数、ミトコンドリアのリードの割合、リボソームのリードの割合、検出された遺伝子のログ数(10)をログ(10)リード数で割った値で測定された細胞の複雑さが含まれます。各メトリックの中央値は、閉じたボックス内に記載されています。原子核の総数:4662個。(B)個々の核のクラスタリングと、フィルタリング後の各クラスターの細胞型マーカー遺伝子の相対的な遺伝子発現を示す対応するDotPlotを表示するUMAP。核の数:3817。この図の拡大版を表示するには、ここをクリックしてください。

補足ファイル1:QCおよびクラスタリング分析のコード。このファイルをダウンロードするには、ここをクリックしてください。

ディスカッション

IMATとの連携には、いくつかの固有の課題があります。その限られたアクセス性に加えて、サンプル材料の収量は非常に少ないことが多く、骨格筋の「汚染」を避けることはほとんど不可能です。最高品質のサンプルを得るためには、生検針を挿入する際に筋筋膜を貫通し(皮下脂肪組織を採取しないようにするため)、採取後すぐに顕微鏡でサンプルを解剖してできるだけ多くの筋肉組織を切除し、その後の組織学、スナップ凍結などの処理を行う必要があります。これは、骨格筋から前向きに収集および解剖されたIMATサンプルを利用する最初の方法です。一方、IMATの以前のsnRNA-seqデータは、IMATだけでなく、既知の脂肪浸潤を伴う骨格筋生検の混合物の核分離によるものでした。

IMATは皮下脂肪組織に比べて密度が高いため、組織構造から核を放出するためには、核分離時に組織を適切に粉砕することが不可欠です。不適切な粉砕は核収率の低下につながる可能性があり、これは限られた量の組織では好ましくありません。さらに、調製した核溶液中の核の正しい数を計算することは、GEM生成のためにチップにロードするのに十分な(ただし多すぎない)核を確保するための重要なステップです。核が多すぎる、または大きな核クラスターや破片が存在すると、詰まりが発生し、最終的なcDNA濃度、ひいてはシーケンシング出力に影響を与える可能性があります。単一核調製物中の核凝集物や破片を回避することは困難であり、サンプルの品質に影響されます。核の凝集および/または破片の存在は、30 μmセルストレーナに追加のフィルタリングステップを含めることで減らすことができます。

RNA分解のリスクを減らすためには、サンプルを迅速に処理することが重要です。核単離のすべての可能なステップで溶液を氷上に保つか冷却し、十分な量のRNAse阻害剤を添加することで、分解を最小限に抑えることができます。各サンプルの処理時間が重要であるため、一度に最大 2 つのサンプルを実行することをお勧めします。さらに、原子核溶液は、原子核が壊れやすくなりすぎないように注意して取り扱う必要があります。壊れやすい核は、GEM生成時にバーストを起こしやすく、その結果、サンプル中に大量の周囲RNAが生じる可能性があります。このプロトコルでは、シリンジ(STARプロトコル16に記載されている)を使用したフィルタリングステップを除外して、核の完全性を維持しました。

ライブラリー調製方法の選択は、研究課題によって異なります。私たちは、液滴ベースの単一核アプローチとナノウェルベースの単一核アプローチの両方を用いて研究してきました。ナノウェルベースのアプローチでは、1200〜1600個の核と、16個の核あたり約3000〜6000個の検出遺伝子をプロファイルします。さらに、このアプローチは全長トランスクリプトミクスをサポートしており、これにより構造的変異、偽遺伝子、およびスプライスバリアントの研究が可能になる27。これに対し、液滴ベースのアプローチでは、10,000個以上の核をプロファイリングすることができますが、その代償として、核ごとに検出される遺伝子の数が少なくなります(2〜3倍少なくなります)。このプロトコルでは、プロファイリングされた核の数が多いほど、存在量が少ない希少な細胞タイプを検出する能力が高まるため、液滴ベースのプラットフォームの1つを使用することを選択しました28

snRNA-seqの主な制限は、細胞質含量の喪失であり、発現解析が核転写産物に限定されます。核は、細胞全体のそれよりも10〜100倍少ないmRNAを含んでおり、これは組織内のすべての細胞タイプを明らかにするために懸念される可能性がある29。しかし、Guptaらによる最近の研究では、核トランスクリプトミクスは、全細胞トランスクリプトミクスと比較して29、培養前駆脂肪細胞と成熟脂肪細胞の生物学的に関連する遺伝子の類似した細胞集団と発現プロファイルを検出したことがわかりました。同様の結果は、他の細胞型でも見出されている30,31

単一核単離は、scRNA-seqと不適合な脆弱細胞の単一細胞(sc)分解能の遺伝子発現パターンを提供するという利点があります。さらに、凍結バイオバンクサンプルを使用することができ、核の単離は組織の酵素分解を必要とせず、細胞内のストレス応答を誘発し、トランスクリプトーム32に影響を与える可能性がある。snRNA-seqを用いることで、脂肪細胞が脂質を多く含み、液滴ベースのアプローチと相性が悪いため、scRNA-seqでは不可能であったIMATに存在する脂肪細胞やその他の細胞種を転写的にプロファイリングすることができました。

IMATのsnRNA-Seqの重要な応用は、さまざまな代謝状態におけるIMATの細胞組成と転写の違いを解読することです。さらに、snRNA-seqデータセットを使用して、IMAT上でバルクRNA-seqをデコンボリューションし、細胞組成の違いをより大規模に評価することができる33。さらに、snRNA-seqデータの擬似時間解析は、成熟脂肪細胞の細胞起源がまだ議論されているIMATなどの組織で特に興味深い細胞成熟の軌跡を特定するのに役立ちます。

開示事項

著者は何も開示していません。

謝辞

著者らは、コロラド大学のBryan Bergman博士が、MoTrIMAT研究(R01AG077956)の図1CのIMAT生検の画像を提供してくれたことに感謝の意を表します。私たちは、代表的な結果セクションにデータが示されているIMATサンプルを提供した筋肉、可動性、老化の研究に感謝しています。国立老化研究所(NIA)は、筋肉、可動性、老化の研究(SOMMA;R01AG059416)およびその補助的な研究SOMMA AT(R01AG066474)およびSOMMA Knee OA(R01AG070647)。研究インフラストラクチャのサポートは、ピッツバーグ大学(P30AG024827)とウェイクフォレスト大学(P30AG021332)のNIAクロード・D・ペッパー高齢者アメリカ独立センターと、ウェイクフォレスト大学の国立トランスレーショナルサイエンス推進センター(UL1 0TR001420)が資金提供した臨床およびトランスレーショナルサイエンス研究所によって部分的に資金提供されました。

材料

この記事で使用された材料の一覧
名前会社カタログ番号コメント
0.2 &マイクロ;m corningシリンジフィルター Millipore Sigma SigmaCLS431229
1.7 mL DNA LoBindチューブEppendorf22431021低結合チューブ
10% Tween 20Bio-Rad1662404
100x protease inhibitorサーモフィッシャーサイエンティフィック78437
10X マグネティックセパレーター10X Genomics230003
10X Vortex Adapter10X Genomics
15 mL 標準チューブSarstedt6,25,54,502
2100 Bioanalyzer AgilentG2939BA
50 mL コニカルチューブSarstedt6,25,47,254
CellRangerGenomicsN/A
Chromium iX アクセサリーキット10X GenomicsPN1000323
Chromium iX Controller10X GenomicsPN1000326
Chromium Next GEM Chip G シングルセルキット 10X GenomicsPN1000127
Chromium Next GEM シングルセル 3' ゲルビーズキット v3.1 10X GenomicsPN1000129
Chromium Next GEM シングルセル GEM キット v3.110X GenomicsPN1000130
Countess 3 自動セルカウンターサーモフィッシャーサイエンティフィックサイエン自動セルカウンター
カウンテスセルカウンティングチャンバースライドサーモフィッシャーサイエンティフィックC10228
DoubletFinderN/A
DPBS (カルシウム、マグネシウムなし)Thermo Fisher Scientific14190144
DTTThermo Fisher ScientificR0861
Dual Index Kit TT Set A, 96 rxns10X GenomicsPN1000215
Dynabeads MyOne SILANE 10X GenomicsPN2000048
Falcon 100 & マイクロ;mセルストレーナーCorning Life Science352360
Falcon 40 µm Cell strainerCorning Life Science352340
グリセリン (グリセロール), 50% (v/v) 水溶液Ricca Chemical Company3290-32
KCLThermo Fisher ScientificAM9640G
Library Construction Kit v3.110X GenomicsPN1000196
MACS SmartStrainers (30µm)Miltenyi Biotec130-098-458
Mastercycler Nexus Gradient Thermal cyclerEppendorf6331000017
MgCl2AmbionAM9530G
Mortar and pestelヘルスケア物流 14075
NucBlue Live Ready Probes 試薬サーモフィッシャーサイエンティフィックR37605
ヌクレアーゼフリー水 (DEPC 未処理)サーモフィッシャーサイエンティフィックAM9930
プロブミン ウシ血清アルブミン 脂肪酸フリー、粉末Sigma-Aldrich820024
Qiagen バッファー EBQiagen19086
リボロック RNAse 阻害剤サーモフィッシャーサイエンティフィックEO0382
スーラN/A
ショ糖Sigma-AldrichS0389
SUPERasin 20 U/µLサーモフィッシャーサイエンティフィックAM2695
サーモミキサー Cエッペンドルフ 
組織ホモジナイザーGlass-Col099C K54
Tris バッファー pH 8.0サーモフィッシャーサイエンティフィックAM9855G
Triton X-100サーモフィッシャーサイエンティフィAC327372500
UltraPure 0.5M EDTA pH 8.0Gibco15575020
330002ティフィック AMQAX2000 5382000015ック

参考文献

  1. Goodpaster, B. H., Bergman, B. C., Brennan, A. M., Sparks, L. M. Intermuscular adipose tissue in metabolic disease. Nat Rev Endocrinol. 19 (5), 285-298 (2023).
  2. Sparks, L. M., Goodpaster, B. H., Bergman, B. C. The metabolic significance of intermuscular adipose tissue: Is IMAT a friend or a foe to metabolic health. Diabetes. 70 (11), 2457-2467 (2021).
  3. Gallagher, D., et al. Adipose tissue in muscle: A novel depot similar in size to visceral adipose tissue. Am J Clin Nutr. 81 (4), 903-910 (2005).
  4. Manini, T. M., et al. Reduced physical activity increases intermuscular adipose tissue in healthy young adults. Am J Clin Nutr. 85 (2), 377-384 (2007).
  5. Addison, O., Marcus, R. L., LaStayo, P. C., Ryan, A. S. Intermuscular fat: A review of the consequences and causes. Int J Endocrinol. 2014, 309570(2014).
  6. Goodpaster, B. H., et al. Obesity, regional body fat distribution, and the metabolic syndrome in older men and women. Arch Intern Med. 165 (7), 777-783 (2005).
  7. Goodpaster, B. H., Thaete, F. L., Kelley, D. E. Thigh adipose tissue distribution is associated with insulin resistance in obesity and in type 2 diabetes mellitus. Am J Clin Nutr. 71 (4), 885-892 (2000).
  8. Goodpaster, B. H., et al. Association between regional adipose tissue distribution and both type 2 diabetes and impaired glucose tolerance in elderly men. Diabetes Care. 26 (2), 372-379 (2003).
  9. Sachs, S., et al. Intermuscular adipose tissue directly modulates skeletal muscle insulin sensitivity in humans. Am J Physiol Endocrinol Metab. 316 (5), E866-E879 (2019).
  10. Ford, H., Liu, Q., Fu, X., Strieder-Barboza, C. White adipose tissue heterogeneity in the single-cell era: From mice and humans to cattle. Biology (Basel). 12 (10), 1289(2023).
  11. Wang, L., et al. Single-nucleus and bulk RNA sequencing reveal cellular and transcriptional mechanisms underlying lipid dynamics in high marbled pork. NPJ Sci Food. 7 (1), 23(2023).
  12. Li, J., et al. Identification of diverse cell populations in skeletal muscles and biomarkers for intramuscular fat of chicken by single-cell RNA sequencing. BMC Genomics. 21 (1), 752(2020).
  13. Lyu, P., Qi, Y., Tu, Z. J., Jiang, H. Single-cell RNA sequencing reveals heterogeneity of cultured bovine satellite cells. Front Genet. 12, 742077(2021).
  14. Fitzgerald, G., et al. MME+ fibro-adipogenic progenitors are the dominant adipogenic population during fatty infiltration in human skeletal muscle. Commun Biol. 6 (1), 111(2023).
  15. Cummings, S. R., et al. The study of muscle, mobility and aging (SOMMA): A unique cohort study about the cellular biology of aging and age-related loss of mobility. J Gerontol A Biol Sci Med Sci. 78 (11), 2083-2093 (2023).
  16. Whytock, K. L., et al. Isolation of nuclei from frozen human subcutaneous adipose tissue for full-length single-nuclei transcriptional profiling. STAR Protoc. 4 (1), 102054(2023).
  17. 10x Genomics. Chromium Single Cell 3' Reagent Kits User Guide (v3.1 Chemistry Dual Index), Document Number CG000315 RevE. , Available from: https://cdn.10xgenomics.com/image/upload/v1668017706/support-documents/CG000315_ChromiumNextGEMSingleCell3-_GeneExpression_v3.1_DualIndex__RevE.pdf (2022).
  18. Heumos, L., et al. Best practices for single-cell analysis across modalities. Nat Rev Genet. 24 (1), 550-572 (2023).
  19. Hao, Y., et al. Dictionary learning for integrative, multimodal and scalable single-cell analysis. Nat Biotechnol. 42 (2), 293-304 (2023).
  20. McGinnis, C. S., Murrow, L. M., Gartner, Z. J. DoubletFinder: Doublet detection in single-cell RNA sequencing data using artificial nearest neighbors. Cell Syst. 8 (4), 329-337 (2019).
  21. Yang, S., et al. Decontamination of ambient RNA in single-cell RNA-seq with DecontX. Genome Biol. 21 (2), 57(2020).
  22. Common considerations for quality control filters for single cell RNA-seq data. 10X Genomics. , Available from: https://www.10xgenomics.com/analysis-guides/common-considerations-for-quality-control-filters-for-single-cell-rna-seq-data (2022).
  23. Luecken, M. D., Theis, F. J. Current best practices in single-cell RNA-seq analysis: a tutorial. Mol Syst Biol. 15 (6), e8746(2019).
  24. Emont, M. P., et al. A single-cell atlas of human and mouse white adipose tissue. Nature. 603 (7903), 926-933 (2022).
  25. Hildreth, A. D., et al. Single-cell sequencing of human white adipose tissue identifies new cell states in health and obesity. Nat Immunol. 22 (5), 639-653 (2021).
  26. Whytock, K. L., et al. Single cell full-length transcriptome of human subcutaneous adipose tissue reveals unique and heterogeneous cell populations. iScience. 25 (8), 104772(2022).
  27. Probst, V., et al. Benchmarking full-length transcript single cell mRNA sequencing protocols. BMC Genomics. 23 (1), 860(2022).
  28. CG000148 Rev A Technical Note - Resolving cell types as a function of read depth and cell number. Technical note. 10X Genomics. , Available from: https://assets.ctfassets.net/an68im79xiti/6gDArDPBTOg4IIkYEO2Sis/803be2286bb a5ca67f353e6baf68d276/CG000148_10x_Technical _Note_Resolving_Cell_Types_as_Function_of_ Read_Depth_Cell_Number_RevA.pdf (2018).
  29. Gupta, A., et al. Characterization of transcript enrichment and detection bias in single-nucleus RNA-seq for mapping of distinct human adipocyte lineages. Genome Res. 32 (2), 242-257 (2022).
  30. Bakken, T. E., et al. Single-nucleus and single-cell transcriptomes compared in matched cortical cell types. PLoS One. 13 (12), e0209648(2018).
  31. Wu, H., Kirita, Y., Donnelly, E. L., Humphreys, B. D. Advantages of single-nucleus over single-cell RNA sequencing of adult kidney: Rare cell types and novel cell states revealed in fibrosis. J Am Soc Nephrol. 30 (1), 23-32 (2019).
  32. Kim, N., Kang, H., Jo, A., Yoo, S. -A., Lee, H. -O. Perspectives on single-nucleus RNA sequencing in different cell types and tissues. J Pathol Transl Med. 57 (1), 52-59 (2023).
  33. Avila Cobos, F., Alquicira-Hernandez, J., Powell, J. E., Mestdagh, P., De Preter, K. Benchmarking of cell type deconvolution pipelines for transcriptomics data. Nat Commun. 11 (1), 5650(2020).

再版と許可

タグ

核分離ヒト脂肪組織細胞タイププロファイリングドロップレットベースシーケンシング組織粉砕セルストレーナーろ過自動細胞計数器代謝性疾患研究