シングルセルRNAシーケンシング
山形方人
Harvard University
DOI：10.14931/bsd.8038 原稿受付日：年月日 原稿完成日：年月日
担当編集委員：
英：single cell RNA sequencing, scRNA-seq
シングルセルRNAシーケンシング(scRNA-seq)は、次世代シーケンシング （Next Generation Sequencing、NGS）技術を使用して個々の細胞が発現しているmRNA全体、つまりトランスクリプトームを質的、量的に網羅的に調べ、細胞ごとの違いを高解像度で検出、分類することで、細胞の分類を行うことができる技術である。また、刺激、発生など細胞の状況に応じて、個々の細胞のトランスクリプトームの情報を得ることで、病態や細胞系譜などの解析も可能である。特に多様なニューロンが存在する神経系では、この方法により、神経細胞の分類や状態について、深い理解が進んできている。
トランスクリプトーム
トランスクリプトーム（transcriptome）は、細胞中に存在する全ての転写産物（タンパク質をコードするmRNA、タンパク質をコードしないノンコーディングRNA、マイクロRNAなど）の総体である。トランスクリプトームは、ゲノムとは異なり、同一の個体でも、組織ごとに、更には発生段階や細胞外からの影響によって固有のものである。このようなトランスクリプトームは、同質あるいは異質の多数の細胞集団（例、培養細胞株、組織）からRNAを抽出し、1990年代に開発されたDNAマイクロアレイのように数多くの既知のmRNAを一気に識別する技術によって解析されるようになった。その後、次世代シーケンシング（NGS）の利用により、希少mRNAやノンコーディングRNAを含めた未知の転写産物の高感度検出も可能になるとともに、スプライシングを経て成熟していく過程のmRNAの構造など、転写産物の種類だけでなく、転写産物の構造の理解も進むことになった。また、NGSは、ヒトやモデル実験生物（マウス、センチュウ、ショウジョウバエなど）として広く利用される生物だけでなく、多様な生物のトランスクリプトームの理解も可能にした。
シングルセルトランスクリプトーム研究史の概観
1つの細胞の持つ生体物質を定量しようとする試みは古くからあった。1960年代になると、Fluorescence-activated cell sorting (FACS)が発明され、標識抗体などのプローブと組み合わせることで、多くの細胞の中で１つの細胞が持っている分子の種類や量についての断片的な研究が可能になり、この方法は現在でも利用されている。その後、免疫組織化学やin situ hybridizationなどにより、タンパク質やmRNAの種類や量が観察できるようになり、組織中の1つの細胞の同定などに活用されてきている。 一つの細胞内にある全RNAは細胞種によるが1-50pgである。そのうち、mRNAの占める割合は1-5%程度である。この微量のmRNAをcDNAに変換してから大幅に増幅できるPCRが発明されることで、1つの細胞が発現するmRNAを高感度で検出できるようになった。例えば、1991年、Linda BuckとRichard Axelは、嗅覚受容体がGタンパク質であると仮定し、個々の嗅覚細胞で特異的に観察されるGタンパク質mRNAを比較することで、嗅覚受容体候補の同定に成功した（年、ノーベル生理学・医学賞）。1995年になると、Catherine DulacとRichard Axelは、異なる鋤鼻神経細胞で特異的に発現する遺伝子を単細胞cDNAライブラリーのディファレンシャル・スクリーニングという方法で、フェロモン受容体候補を同定した。同様な手法で異なる種類の網膜神経節細胞で発現している遺伝子も同定されており、このようなアプローチが生理的に重要な機能を持つ遺伝子の発見に効果的であることを示した。 一方で多くの種類のmRNAを１細胞レベルで観察する単細胞トランスクリプトームには技術的なブレークスルーが待たれた。１つ大きな問題はPCRなどの増幅に伴うバイアスなどのアーティファクトが頻繁に観察されること、そしてもう一つの課題は多くの種類のcDNAを簡便に観察することを可能にする方法の開発であった（PMID: 16547197）。これを可能にしたのが、増幅法の改良とマイクロアレイの利用であった。しかしながら、増幅に伴うアーティファクトの解決は依然として不十分で、また１つの細胞ごとに高価なマイクロアレイを利用することは、多数の細胞のトランスクリプトームを観察するのには限界があった。2009年に、これらの問題を解決できる可能性として、High-throughput sequencing (HTS)を利用するscRNA-seqプロトコルがAzim Suraniのグループによって報告されたPMID:19349980。しかしながら、この方法でも一つの細胞ごとに処理を行うという操作が必要で、この論文でもたった8個の細胞の解析に留まっており、非常に多くの細胞についてのトランスクリプームを一挙に理解することはできなかった。
シングルセルRNAシーケンシング(scRNA-seq)の現状
以来、完全長cDNAまたは分子識別子（unique molecular identifiers: UMI）を持つcDNAを生成するためにmRNA転写産物を増幅する方法が異なるscRNA-seqが考案されてきた。2013年には、このような１細胞のシーケンシング技術が、Nature Methods誌のMethod of the Year に選ばれた。たとえば、SMART-seq（Switch mechanism at the 5' End of RNA Templates）（ 18 ）およびその改良されたプロトコルであるSMART-seq2（ 19、20 ）は、完全長cDNA合成のためのプロトコルである。また、MARS-seq（並列RNA単一細胞配列決定）（ 21 ）、STRT（単一細胞タグ付き逆転写）（ 22、23 ）、CEL-seq（線形増幅および配列決定による細胞発現）（ 24 ）、CEL-seq2（ 25 ）などが報告されてきた。特にSMART-seq（SMART-seq2）は、ピペット、限界希釈、レーザー捕獲法などを用いる多穴プレート法、更に半導体集積回路製作技術で作った流体集積回路を利用するFluidigm C1のシステム（https://jp.fluidigm.com）と組み合わせることで利用される機会が多い。このプロトコールの特徴は、全長のトランスクリプトームを得ることができることであり、mRNAのスプライシングバリアントなどのアイソフォーム、SNPs、変異の検出にも利用できる点で次に説明するUMIを用いる方法に比べて利点があるが、そのコストと処理できる細胞数の点で極めて不利である。
最も重要なscRNA-seqの方法論についての進歩は、2015年、Harvard Medical Schoolの独立した２つのグループから、inDropそしてDrop-seqという類似した2つの方法が発表されたことであろう。マイクロ流体力学 (Microfluidics) 、 UMIとしてDNAバーコーディング (DNA barcoding) 、そしてNGSを利用することで、自動化とサンプル調製の容易さから、一つの細胞あたりに要するコストを大幅に低下させることに成功した（Drop-seqはその発表時で、6セント/細胞）。これらの方法では、細胞1つずつをマイクロ流体力学によるエマルジョン技術を利用した装置に流入させ、その１細胞を試薬を封入した1つのDroplet（油滴）に自動的に閉じ込める。そのDroplet中には、DropletごとにUMIとして異なったDNAバーコードを持つゲルビーズが入っており、そこからcDNA合成反応を行うことで、それぞれの同じ細胞に含まれていたmRNAが同じUMIを持つcDNAとして合成され、Dropletを破壊した後も、そのcDNAが由来した細胞が区別できるということを利用している（図１）。このようにして増幅したUMI付きのｃDNAをNGSで配列決定することによりscRNA-seqが可能になる。なお、DropSeqはコストが低いが、細胞の取得率と検出感度が低い弱点がある。inDropはDropSeqより細胞取得率が高く、パラメータを調整することで低レベルで発現される遺伝子の検出に有利である。 inDropの方法は、1 Cellbio社（https://1cell-bio.com）から販売されているが、特に重要なのは同様の原理を用いた10xGenomics社（https://www.10xgenomics.com/jp/）がChromiumと命名された市販機器と試薬を発売することで、多くの研究者に利用できることになったことである。Svenssonらのデータベース（www.nxn.se/single-cell-studies/gui）では、scRNA-seqを用いた論文で用いられた方法について調査しているが、この数年、10xGenomics社のChromiumを用いた方法が飛躍的に増加し、ほぼ寡占状態になりつつあることがわかる（現在、10XGenomics社とBioRad社の間で関連特許をめぐる係争がある。）。10X Genomics Chromiumは市販であるので導入が容易であり、inDropやDropSeqに比べ最大数の転写産物の検出に敏感であるが、コストが高い。
シングルセルRNAシーケンシング(scRNA-seq)の実際
ここでは主流になっている10xGenomics社のChromiumを用いた方法とSMART-seqなどを用いた方法に共通する方法の実際について議論する。シングルセルRNAシーケンシングの利用には、4つのステップがある。１）細胞をバラバラに単離すること。２）ライブラリーの作製とNGSシーケンシング。３）前処理（preprocessing、得られた配列の整理）。４）ダウンストリーム分析（生物学的な情報を得る）。これらのうち、２）の段階については、上に記述したように市販の機器や試薬を利用することが多くなっているので、各社のマニュアル等を参考にするのが現実的である。
組織からの細胞の分離=
血液細胞のように浮遊した細胞ではない場合、物理的あるいは酵素処理などによって、生組織から状態の良い細胞をdissociationする必要がある。神経系組織の酵素処理には、パパインを用いる方法が広く用いられている。ただ、しばしば問題となるのが、酵素処理のため短時間加温することで、発現が変化する遺伝子が存在することである。例えば、脳のミクログリアの解析には、低温下で組織をホモゲナイズするなどの工夫が必要であった（Hammond et al., 2019）。また、酵素処理時に転写阻害剤であるアクチノマイシンで処理することで、このような現象を抑制できる（Wu et al., 2017）。更に、ヒマラヤ氷河から得られた細菌Bacillus licheniformisから得られた低温プロテアーゼを用いる方法も報告されている（Adam et al., 2017）。 単離した細胞は、そのまま10xGenomicsのChromiumのプラットフォームに導入することができるが、抗体などを用いたFACS、パニング、磁気ビーズカラムなどによる細胞の単離を行う場合もある。 なお、ヒト組織などから生細胞を得ることは困難なことが多い。この場合、凍結した組織から、核を調製し、これを分析する方法がある（）。更に、RNAを分析するscRNA-seqではないが、シングルセルの遺伝子発現を推定する方法として、トランスポゾンを用いることでゲノムのオープンクロマチン領域を選択的に検出し、ライブラリーを作製しシーケンスするATAC-seq (Assay for Transposase-Accessible Chromatin)がある。
scRNA-seqデータの前処理
Seurat, Scanpyなどのソフト。 Transcriptomeとの照合。質のチェック。 視覚化（Visualization。tSNE。
ダウンストリーム分析
Dimensionality Reductionとクラスタリング。 マーカー遺伝子とクラスタリングの同定。 DE遺伝子、の検出、MAST。 組成解析。 Trajectory interference 発生。発現の動態。
神経科学への応用
ニューロンのクラスとタイプ
様々な神経・精神疾患について理解しその診断や治療に役立てるためには、ニューロン、グリア細胞を中心にした神経系にある細胞の「タイプ」を識別し、それぞれの細胞における分子的な変化を観察することが重要である。近年、中枢神経系のグリア細胞にも、多様なアストロサイト、オリゴデンドロサイト、ミクログリアの存在が報告されてきている。一方で、ニューロンは著しく多様であり、このニューロンの多様性こそが、神経系を特徴づけており、その多彩で複雑な機能の発現に必須であることは疑う余地がない。 解剖学的な視点から言えば、すべてのニューロンの存在する位置は異なるので、すべてのニューロンは異なるという見方もできる。しかし、これは極論であり、従来の神経科学では、ニューロンの多様性は、それぞれのニューロンの解剖学的な位置、発現している分子、電気生理学、結合性、形態、神経伝達物質、神経伝達物質受容体とシグナル伝達によって識別されてきた。こうしたニューロンの多様性を便宜的に記述するのに、タイプ（type）、クラス（class）、サブクラス（subclass）、サブタイプ（subtype） というような用語が用いられてきた。しかし、本稿では混乱を防ぐため、Masland（2004）の提唱に従い、「クラス」と「タイプ」という単語を用いることとする。タイプは、これ以上分類することができないとされる階層である。例えば、大脳皮質の錐体細胞、網膜神経節細胞といった大雑把な識別は「クラス」と呼ぶ。大脳皮質の錐体細胞というクラスは、層や領野によって「タイプ」が異なるし、網膜神経節細胞には視覚情報によって応答が異なる「タイプ」が存在する。この分類は、免疫組織化学、形態、電気生理学などの技術により識別可能である暫定的なものに過ぎない。本稿で解説するscRNA-seqの技術は、その網羅性からそれぞれのニューロンについてこれまでにないビッグデータを提供することで、このニューロンのタイプの理解に確実な根拠を与えつつある。
大脳
その他のCNS
疾患
アルツハイマー、Autism
網膜
展望
データベース Human Cell Atlas Human Brain Transcriptome project Single cell portal Allen Brain Atlas
空間トランスクリプミクス
関連項目
- [[]]
- 大脳皮質
- 樹状突起
- 軸索
- 大脳皮質の局所神経回路