シングルセルRNAシーケンシング

山形方人
Harvard University
DOI：10.14931/bsd.8038　原稿受付日：2020年10月22日　原稿完成日：2020年12月23日
外部査読委員：京都大学メディカルイノベーションセンター 渡辺　亮
理化学研究所 生命機能科学研究センターバイオインフォマティクス研究開発チーム/東京医科歯科大学 難治疾患研究所 ゲノム応用医学部門 ゲノム機能情報分野 二階堂 愛

英：single-cell RNA sequencing, scRNA-seq

　シングルセルRNAシーケンシング(single-cell RNA sequencing, 以下scRNA-seq)は、次世代シーケンサー（next generation sequencer、以下NGS）を用いることで、個々の細胞が保持しているmRNA全体を質的、量的に網羅的に調べる方法である。次元圧縮などの数理的な解析と組み合わせることで、遺伝子発現の状態に基づいた細胞の分類を行うことが可能であり、従来の組織学的、あるいは細胞生物学的手法では知られなかった新規の細胞種の同定や細胞状態の推定を行うことが可能になった。また、遺伝子発現プロファイルの変化に基づく擬時系列解析(pseudotime analysis)によって、刺激や発生に伴う細胞状態の遷移の描写ができる。神経系では、この方法により、神経細胞や非神経細胞の分類や状態についての知見が深まり、新しい神経細胞タイプ、細胞マーカー、病態の理解、更に機能的な遺伝子の同定などが系統的かつ網羅的に行われるようになった。scRNA-seqに、空間的情報、エピゲノム情報、タンパク質情報などの複数モダリティを取り入れた統合解析（multimodal single-cell omics）も行われている。

背景

トランスクリプトーム

　トランスクリプトーム（transcriptome）は、細胞中に存在する全ての転写産物（タンパク質をコードするmRNA、タンパク質をコードしないノンコーディングRNA、マイクロRNAなど）の総体である。トランスクリプトームは、ゲノムとは異なり、同一の個体でも、組織ごとに、更には発生段階や細胞外環境や刺激によって変化する。トランスクリプトームは、同質あるいは異質の多数の細胞集団（組織、培養細胞）からRNA抽出後、cDNAに変換し、それを1990年代に出現したDNAマイクロアレイのように数多くの既知mRNAを識別する技術によって解析されるようになった。その後、NGSの利用により、希少mRNAやノンコーディングRNAを含めた未知の転写産物の高感度検出が可能になるとともに、スプライシングで成熟していく過程のmRNAなど、転写産物の種類だけでなく、転写産物の構造的差異（スプライシングバリアント、SNPs、変異など）の解析もできるようになった。加えて、ヒトやモデル実験生物（マウス、ゼブラフィッシュ、ショウジョウバエ、線虫など）だけでなく、多種多様な生物のトランスクリプトームの把握も可能になった。従来から行われてきた組織全体などの多数の細胞を対象としたRNA-seq（バルクRNA-seq）では、複数の細胞における転写産物の平均を観察しているが、本項目では個々の細胞における転写産物を解析するscRNA-seqの原理とその応用について概説する。

開発史

　1つの細胞の持つ生体物質を解明し、定量しようとする試みは古くからあった。1960年代になると、フローサイトメトリーを利用した蛍光活性化セルソーティング（Fluorescence-activated cell sorting, FACS）が発明され、標識抗体などのプローブと組み合わせることで、多数の細胞集団の中で１つの細胞が保持している生体分子の種類や量についての断片的な研究が可能になり、この方法は現在でも汎用されている。その後、免疫組織化学やin situ hybridizationnなどにより、タンパク質やmRNAの種類や量が観察できるようになり、組織中に存在するそれぞれの細胞の同定などに活用されてきている。最近では、それぞれの細胞が持つ抗原分子を、異なった金属イオンで標識した抗体とフローサイトメトリーを組み合わせた方法で検出するマスサイトメトリー（CyTOFなど）も開発されてきている。

　細胞種にもよるが、1つの細胞内にある全RNA（ribosomal RNAを含む）は細胞種にもよるが1-50pgである。そのうち、mRNAの占める割合は1-5%程度である。この微量のmRNAをcDNAに変換してから大幅に増幅できる方法が発明されたことで、1つの細胞が発現するmRNAを高感度で検出できるようになった 。例えば、1991年、Linda BuckとRichard Axelは、嗅覚受容体がGタンパク質であると仮定し、個々の嗅覚細胞で特異的に観察されるGタンパク質mRNAを比較することで、嗅覚受容体の同定に成功した。1995年になると、Catherine DulacとRichard Axelは、異なる鋤鼻神経細胞で特異的に発現する遺伝子を１つの細胞から作製したcDNAライブラリーを比較するディファレンシャル・スクリーニングを行うことで、フェロモン受容体を同定した。同じ手法で異なる種類の神経細胞で発現している遺伝子も同定され、１つの細胞の持つトランスクリプトームを比較するアプローチが神経細胞で特徴的に発現している遺伝子の同定に効果的なことが示された。

　一方で多くの種類のmRNAを１細胞レベルで一挙に観察するための技術には感度やスループット、そしてコストの観点からブレークスルーが待たれた。1つの問題は多種類のcDNAを簡便に識別することを可能にする方法の開発であった。これを可能にしたのが、PCRなどのcDNA増幅法の改良とマイクロアレイの利用であった。しかしながら、細胞ごとに高価なマイクロアレイを使用することは、多数の細胞のトランスクリプトームの観察には限界があった。2009年になると、これらの問題を解決できる可能性として、NGSを利用するscRNA-seqプロトコールがAzim Suraniのグループによって報告された。しかしながら、多数のマイクロアレイでなく1回のNGS使用で済ませることができるものの、この報告でもわずか8個の細胞の解析に留まっており、1つの細胞ごとに処理を行うという操作が必要で、多数の細胞についてのトランスクリプームを一挙に理解することはできなかった。また、塩基配列の違うcDNAごとにPCR効率に差がある結果生じる増幅バイアス、また3’末端側が選択的に補足されることなどの課題があった。

現状

分子生物学的反応

　その後、5’末端側の領域まで効率よく増幅するscRNA-seqのプロトコールが考案された。特に、SMART-seq（Switching mechanism at the 5' End of RNA Templates）およびその改良されたプロトコールであるSMART-seq2 の使用例が多い（既に、SMART-seq3という改良プロトコールもあるが、以下SMART-seqと呼ぶ）。また、類似法としてSTRT（single-cell tagged reverse transcription）などがある。

　一方、CEL-seq（Cell Expression by Linear amplification and Sequencing）、CEL-seq2 、MARS-seq（Massively parallel single-cell RNA-seq）では、T7 RNAポリメラーゼによるin vitro転写を用いることにより、ＰＣＲによる増幅で見られるバイアスを低減させようとしている。

　また、Quartz-SeqやQuartz-Seq2ではPCR用のアダプターを付加する反応にポリAテーリングを利用することで、他の手法と比較して1.5-5倍程度の遺伝子を検出できる。

バーコード技術

　増幅バイアス除去のアプローチとして特に重要なのは、2011年に発表された核酸配列バーコードを利用した方法で、分子識別子（unique molecular identifiers: UMI）を持つcDNAを増幅させ、NGS後の情報処理を用いるものであると考えられる。この方法では逆転写反応の際、ランダム塩基配列から構成されるUMIをcDNA末端に付加した後、増幅反応、NGSを行い、cDNA配列とUMI配列の両方を読む。cDNAにはRNA1分子に1つのUMIが付加されるので、同一のUMIを持っていれば、逆転写時に同一のcDNA由来とカウントする。UMIをカウントすることで、増幅前のmRNAのコピー数を知ることができる 。

多様なプラットフォーム

　細胞を分別するプラットフォームには、マイクロピペットによる捕獲、セルソーター、レーザー捕獲などを用いるマルチウェル法、あるいは半導体集積回路様の製作技術で作った流体回路を利用するFluidigm C1の装置（C1 Single- Cell Auto Prep ）[4] 、更にドロップレット使用（下記）などがある。これらは、SMART-seqと組み合わせて利用されることが多い。SMART-seqプロトコールの特徴は、全長mRNAのトランスクリプトーム情報を得ることができることであり、mRNAのスプライシングバリアントなどのアイソフォーム、SNPsの情報を利用したアリル特異的発現、変異の検出にも利用できる。また、それぞれ細胞ごとの反応を独立した場所で行うため、反応中に別の細胞の反応と混じる可能性が低い。小型のナノウェルを用いるSeq-Wellも同様に反応自体が混じる可能性が低い。これらの点が、次に説明するドロップレットを使用して3’末端のみを標的にしたscRNA-seqと比べた場合の長所であるが、その高コスト（1細胞あたり数十ドル）と処理可能な細胞数の少なさが短所である。

　これらとは別に、ハイスループットで安価な方法として、それぞれの細胞を独立に標識するのではなく、プールされた細胞を異なるウェルにランダムに振り分け、ウェル固有のバーコードで転写物を標識していく操作を複数回繰り返していくことで細胞を区別するSplit-seqやsci-RNA-seq3などの方法も用いられている。

ドロップレット使用の3’エンドリード法

　scRNA-seqのプラットフォームと方法について重要と考えられる進歩は、2015年、Harvard Medical Schoolの独立した２つのグループが、inDropsそしてDrop-seqという類似した2つのハイスループットな方法を開発したことであろう（inDropsはT7RNAポリメラーゼ、Drop-seqはPCRで増幅）。これらの方法では、マイクロ流体力学 (Microfluidics) 、 UMI（上述）と細胞ごとのバーコード（Cell Barcode）という2種類のDNAバーコーディング、そしてNGSと情報解析法を利用している。そして、多く細胞のサンプル調製の自動化と容易さから、1つの細胞あたりに要するコストを大幅に低下させることに成功した（Drop-seqは発表時で、１細胞あたり約5セント）。つまり、細胞1つずつをマイクロ流体力学によるエマルジョン作製技術を利用した装置に流入させ、その１細胞を1つのドロップレットに自動的に閉じ込める。そのドロップレット中には、ドロップレットごとにCell barcode/UMIとしてユニークなDNAバーコードを持つゲルビーズ(Gel Beads in Emulsion, GEMs)が入っており、それを足場に3’末端のみを標的にしたcDNA合成反応を実施することで、同じ細胞に含まれていた1分子のmRNAが同じCell barcodeを持つcDNAとして合成され、そのmRNA/cDNAが由来した細胞を識別できるということを利用している（図１）。

図1．ドロップレット使用の3’エンドリード法
組織から解離させた細胞それぞれを、マイクロ流体力学を利用した装置で、バーコードプライマーが結合したゲルビーズとともにドロップレットに封じ込める。ドロップレット中には、ドロップレットごとにCell barcode/UMIとしてユニークなDNA配列を持つゲルビーズ(GEMs)が入っており、それを足場にcDNA合成反応を実施することで、同じ細胞に含まれていたmRNAが同じCell barcodeを持つDNAとして合成され、それを増幅する。

　DropSeqはコストが低いが、細胞の取得率と検出感度が低い弱点がある。inDropsはDropSeqより細胞取得率が高く、パラメータを調整することにより、低レベルで発現される遺伝子の検出にも有利であるとされる。DropSeqのセットアップはDolomite Bio ([5])、inDropは1 Cellbio社から販売されている[6]。しかし、その後、10x Genomics社が同様の原理を用いたシングルセル遺伝子発現解析システムを市販することで、多くの研究者が利用できるようになっている[7]。Svenssonらによる最近のデータベース [8]では、scRNA-seqを用いた論文で用いられた方法について調査しており、この数年、10x Genomics社のプラットフォームを用いた論文が飛躍的に増加していることがわかる。10x Genomics社のプラットフォームは市販であるので導入が容易であり、DropSeqやinDropsに比べ多くの転写産物の検出が可能であるが、それらよりランニングコストは高価である。

　なお、3’エンドリード法だけでなく、抗体やT 細胞レセプターのN末端側に位置する可変領域の配列決定が可能である5'末端のシーケンシングには5’エンドリード法が利用されることがある。

実際

　ここでは主流になっている10x Genomics社のChromium controllerなどのドロップレットを用いた方法とSMART-seqなどを用いた他のプラットフォームに共通する方法の実際について概説する。scRNA-seqの利用には、4つのステップがある（図２）。

図２．scRNA-seqの実際のステップ
細胞の単離、ライブラリ作製とNGS、データの前処理から次元圧縮、データ解析。図の一部は2016 DBCLS TogoTV、あるいはSeuratを用いて10x Genomics社のPBMCデータ（[1]から執筆者が作製。

1. 個体や組織を採集し、そこから細胞あるいは細胞核を個別に解離された状態にすること。
2. ドロップレット法やSMART-seq対応のプラットフォームなどによる個々の細胞からのライブラリーの作製とNGS。
3. 得られた配列情報の前処理（preprocessing）。
4. データ解析。これらのうち、2.の段階については、上に記述したように市販の機器や試薬を利用する機会が多くなっているので、詳細は説明しない。

組織からの細胞、細胞核の分離

　浮遊細胞（血液細胞など）ではない場合、物理的あるいは酵素処理などによって解離することで、生組織から状態の良い個々に分散した細胞を調製する必要がある。神経系組織の酵素処理には、パパインを用いる方法が広く用いられている。ここで、しばしば問題となるのが、酵素処理による短時間加温や機械的刺激で、発現量が変化する遺伝子が存在することである。特に、脳のミクログリアの解析には、低温下で組織をホモゲナイズするなどの工夫が必要であった。また、このような現象を抑制するために、酵素処理時に転写阻害剤であるアクチノマイシンで処理したり、ヒマラヤ氷河から得られた細菌Bacillus licheniformisから得られた低温プロテアーゼを用いる方法も報告されている。また、細胞解離後に、メタノールで固定しscRNA-seqに使用したり、クロスリンカーを用いる方法もある。

　単離した細胞は、そのまま10x Genomicsのシングルセル遺伝子発現解析のプラットフォームに導入することができるが、細胞表面分子マーカーに対する抗体や蛍光タンパク質レポーターなどを用いたFACS、パニング、MACS(磁気ビーズカラム)などによって、細胞の選択的濃縮や除去を行う場合もある。更に、抗体に抗体表示バーコードDNAをカップリングさせるCITE-seq(Cellular Indexing of Transcriptomes and Epitopes by Sequencing) については、下記の「統合解析」でも述べる。

　なお、ヒト組織や希少生物などから生細胞を得ることは困難なことが多い。この場合、scRNA-seqの変法として、凍結した組織から、各細胞由来の核を調製し、核内のmRNAを分析するsnRNA-seq (single-nucleus RNA-seq)が利用されている。ただ、snRNA-seqでは、FACSなどによる特定細胞集団の分離が困難であることが多い。また、細胞質を持つ生細胞を利用したscRNA-seqとは違って、スプライシングの途上にある未成熟な核内転写産物を検出すること、更に検出できる遺伝子数も少なく、同等な結果が必ずしも得られない 。一方で、snRNA-seqでは、組織をそのまま凍結することから開始することが可能であるので、上述したscRNA-seqの問題である細胞解離酵素による処理などを避けることができる。更に、核を用いることで、大きな細胞体はマイクロ流体力学の流路で詰まりやすいなど、特に神経細胞で顕著である細胞の形状の多様性に伴うバイアスを減らすことができるといったメリットもある。こうしたプロトコールの一部は、protocols.ioのHuman Cell Atlasのグループ[9]で公開されている。 通常のscRNA-seqは、ポリアデニル化されたmRNAを検出しているが、MATQ-seq（multiple annealing and dC-tailing-based quantitative single-cell RNA-seq）、RamDA-seqなどを用いると、ポリアデニル化されていないRNAの検出も可能である [10]。

　更に、RNAを分析するscRNA-seqではないが、遺伝子発現状態との関係が想定されるオープンクロマチン領域をトランスポゾンを用いることで個々の細胞レベルで選択的に検出するsingle cell ATAC-seq (Assay for Transposase-Accessible Chromatin) , , single cell THS-seq (transposome hypersensitive-site) や DNAメチル化領域を観察するsnmC-seq、RRBSのような方法も利用されている。

データ処理の流れ

総論

　Illumina社に代表されるNGSを用いて得られた結果は、ベースコールや細胞バーコードを用いたdemultiplexingなどの基礎解析を行うことで、各細胞における遺伝子の発現量のマトリックスを出力する。例えば、10XGenomics社のChromiumプラットフォームを用いた場合、10XGenomics社が提供するCell Rangerのmkrefコマンド(Linux上で作動, [11])などにより、各生物種ごとのレファレンス配列リスト（[12]、マウスやヒトでは既製のものを利用できる）などを参考にしながら、細胞と転写産物量の対応マトリックスを作製する。その後のデータの処理についても、10x Genomics社がソフトウェアLoupeを提供している。しかしながら、その後のデータ解析を考慮して、R, Python, MATLABなどのデータ解析のための汎用プログラミング言語やコードで扱えるオブジェクトに変換するのが通常である。

　scRNA-seq解析のためには、数多くのツールが公開されている。これらのツールは、バージョンが更新されたり、新しいものに置き換えられることがあるので、実際に利用する場合は最新の動向に注意を払う必要がある。scRNA-seqの解析に必要なツールは、scRNA-tools [13], Awesome single cell [14], Bioconductor[15]などで紹介されており、ほとんどがダウンロード可能である。また、bioRxivなどの査読前のプレプリントサーバで公開されて、随時試用、評価されていくものが多く、scRNA-seqのデータ（下記参考）とともに、オープンサイエンス実践の好例となっている。

Seurat

　ここでは、scRNA-seqデータ解析のために最もよく利用されているRを用いたパッケージ「Seurat」 を中心に紹介しておきたい。なお、一部の解析操作は、University of WashingtonのCole Trapnell研究室で開発されてきた軌道推定（下記参考）によく使用されるMonocle3でも可能である（[16]）。Pythonを利用したものでは、ドイツ・ミュンヘンInstitute of Computational Biologyの Fabian Theisらが開発しているScanpyが有名である。

　New York UniversityのRahul Satija研究室が開発しているSeurat（画家スーラに由来）は、scRNA-seqデータ解析のために広く利用されているRパッケージであり、2020年秋現在、Seurat4のβバージョンが公開されている。論文の正式発表前から、サポート情報提供やコード修正なども頻繁に行っており、Satija研究室のウェッブサイト（ [17]）、Github([18])、更にTwitterアカウント（@satijalab）などで最新情報を得ることできる。

品質の検討事項

　最初に行うのは、scRNA-seqデータの品質管理である。ここでは、質の低い細胞のデータ（例えば、壊れた細胞では、転写産物の種類が少なくミトコンドリア由来の転写産物が多い）を取り除く。また、複数の試料を組み合わせる場合には、バッチごとの違いについて検討する 。現実には、実験ごとのバッチの違いによる影響（Batch effect）がscRNA-seqの最大の問題であると示されてきており、試料の処理を同時に行うなど実験デザインを工夫する必要がある。

　また、ドロップレットを使用するscRNA-seqでしばしば問題になるのが、ドロップレットに２つ以上の細胞が封じ込められ、それらが同一のCell barcodeを持ってしまうアーティファクトである。通常「Doublet」「Multilet」と呼ばれるこの問題はデータ解析を混乱させるので、細胞単離の段階から注意する必要があるが、scRNA-seqデータ取得後にもデータ処理で解決できる可能性もある [19][20][21]。なお、この手法を利用することで、バッチ効果を抑えるために、異なるバーコードを持つ複数の試料を混ぜて一つの試料として扱い、計算機的に再び分離、解析する手法が注目されている。

　scRNA-seqデータの次のノイズは、ある遺伝子の発現が低いために、本来同じタイプの細胞であっても、その遺伝子の発現が全く見られない「Dropout」と呼ばれる現象であり解析に影響を与えるので、これについても検討が必要である。

次元圧縮

　このような品質管理、ノーマライゼーションの過程を経て、scRNA-seqのデータ解析において、最初に行うのが、次元圧縮 (dimensionality reduction)である。PCA (Principal component analysis, 主成分分析)、更に発展させたUMAP（Uniform Manifold Approximation and Projection, 均一マニフォールド近似と投影）、Diffusion maps, tSNE(t-distributed Stochastic Neighbor Embedding , t分布型確率的近傍埋込み)などの手法が用いられる。 特に、tSNE[22]とUMAP[23]は、高次元データを低次元の点の集合として可視化することで、それぞれの細胞の持つ遺伝子発現状態の類似度についての直観的な表示が可能でありしばしば用いられる（図３）。tSNEよりUMAPの方が迅速に類似集団間の関係が明確になるので、最近はUMAPを利用することが多くなってきている。次に、Louvainアルゴリズムなどでクラスタリング（コミュニティ分割）を行いグラフ上に表示できる（図３の色分け）。こうして、異なる転写状態を示す細胞の集合が別のクラスターとして表示され、同定可能になる。

図３． tSNEとUMAPによる同じデータの可視化
網膜（ニワトリ）の視細胞のデータを用いて執筆者が作製[2]。

データ解析

細胞クラスターの解釈とマーカー遺伝子候補の発見

　scRNA-seqデータから得られる生物学的知見には、内在的に存在する細胞の種類、外部刺激や環境で変化した細胞の状態、そして種類や変化により特徴的に発現するマーカー遺伝子候補の発見がある。クラスタリングにより、異なった細胞集団の存在が認識されると、それぞれのクラスター（群）に特徴的に発現している具体的な遺伝子を探索し、細胞集団の持つバイオマーカーによって、そのクラスター（群）の同定が可能になる。例えば、既に神経細胞とグリア細胞に特異的に発現する典型的マーカーはよく知られており、それぞれのクラスターの識別は容易である。更に、神経細胞のタイプ（下記参考）を区別できるマーカーや、外部刺激によって遺伝子発現状態が変化した神経細胞の状態は、In situ hybridizationや免疫組織化学などにより確認できる。このようなクラスターごとに発現が異なる遺伝子（差次的発現遺伝子）を見つけるためには（Differential expression analysis, DE analysis）、SeuratのFindMarkersコマンド中でも利用可能であるコード（MAST、DESeq2など）を用いることができる。細胞ごとの差次的発現遺伝子のVisualization（表示可視化）には、ドットプロット（dot plot）、ヴァイオリンプロット（violin plot）、リッジプロット（Ridge plot, joy plot）、UMAPなどの次元圧縮図上に転写物量をプロットするFeatureプロット（feature plot）などが、目的に応じて頻繁に用いられる（図４）。

図４．scRNA-seqデータの可視化の例
A. ドットプロット。B.ヴァイオリンプロット。C. リッジプロット。D. UMAP（灰色）に転写物量（青）をプロットした Featureプロット。網膜の視細胞のデータを用いて執筆者が作製[3]。

擬時系列解析

　実験的なノイズとは別に生物学的に意味のある遺伝子発現の変動には、位置情報、細胞周期、概日リズム、発現変動が大きい破裂型プロモーターの作動などの理由で 変動が見られるものもある。特に、刺激・薬剤処理やさまざまな病態の進行や治療に伴う細胞の変化、発生途上にある細胞系譜や細胞分化といった細胞の遷移状態の解析（軌道推定（Trajectory inference）；擬時系列解析（擬似時系列解析）、Pseudo-time analysis ）には、scRNA-seqデータを用いることが効果的である。しばしば用いられるMonocle3 [24]など、多くのコードを収集、比較しているサイトがある [25][26]。RNA velocityといったスプライシングされていく転写産物の量から細胞の分化状態を推定する方法もある。しかし、これらの方法は、あくまで細胞系譜や細胞分化の推定に過ぎない。細胞系譜を更に確実に観察しつつ、scRNA-seqを行うことで、細胞タイプの系統関係を調べる方法として、CRISPR-Cas9を用いたゲノム編集による痕跡追跡記録法を導入したscGESTALT、ScarTrace 、LINNAEUS、あるいはアデノシンデアミナーゼでRNA編集を行いタイムスタンプを入れる方法がある。1塩基バリアント（Single-nucleotide variants: SNV）の系統的解析は、細胞の不均一性や系統的な関係を明らかにするための最も有望なアプローチの一つである。

遺伝子制御ネットワーク、パスウェイ解析など

　また細胞分化や刺激などによる変動に伴う特徴的な遺伝子発現状態をscRNA-seqで観察することは、遺伝子制御ネットワーク（例、SCENIC, [27]）、代謝経路やシグナル伝達系のためのパスウェイ解析(例、Metascape, [28]、Gene Ontolgoy[29])を理解するシステム生物学的な研究として有用である。更に、scRNA-seqで得られた結果をもとに、細胞間相互作用の理解を深めるのを目的とするCellPhoneDB[30]、NicheNet,　SVCAなどがある。特に、Perturb-seq やその変法は、CRISPRライブラリーによるゲノム編集を施した細胞をscRNA-seqで解析することで、ゲノム編集で破壊された遺伝子の機能や遺伝子間の相互作用の理解を可能にしている後述する複数モダリティ情報を取り込んだscRNA-seqの１つであり、注目されている。

神経科学研究への適用

神経系細胞ビッグデータとしてのscRNA-seq

　様々な神経・精神疾患について理解しその診断や治療に役立てるためには、神経細胞、グリア細胞を中心にした神経系にある細胞の種類や状態を識別し、それぞれの細胞における分子的な変化を観察することが重要である 。本項目で解説してきたscRNA-seq技術は、神経系に見られるそれぞれの細胞のトランスクリプトームについてビッグデータを提供することで、この細胞の種類や状態の識別に新たな判断材料を与えつつある。近年、中枢神経系のアストロサイト、オリゴデンドロサイト、ミクログリアといったグリア細胞も均一ではなく、内在的な多様性や外部因子による状態の変動が報告されてきている。神経細胞は、著しく多様であり、この多様性が神経系の多彩で複雑な機能発現の基盤となっている。従来の神経科学では、神経細胞の多様性は、それぞれの神経細胞の解剖学的な位置、発現している分子、電気生理学、結合性、形態、神経伝達物質、神経伝達物質受容体とシグナル伝達によって識別されてきている。こうした神経細胞の多様性を便宜的に記述するのに、タイプ（type）、クラス（class）、サブクラス（subclass）、サブタイプ（subtype） というような用語が用いられてきた。しかし、ここでは混乱を防ぐため、Masland（2004）が提唱し、広く受けいれられている「クラス」と「タイプ」という単語を用いることとする。タイプは、これ以上分類することができないとされる階層であり、共通性を持つタイプの集団がクラスである。例えば、大脳皮質の錐体細胞、網膜神経節細胞といった大雑把な区分はクラスである。大脳皮質の錐体細胞というクラスは、層や領野によって異なるタイプ、網膜神経節細胞には視覚情報に対して応答が異なるタイプが存在する。scRNA-seqは、「タイプ」の理解に新たな視点を提供している。

神経系へのscRNA-seqの適用

　scRNA-seqの神経系での利用については、次々と新しい論文やプレプリントが発表されており、ここではscRNA-seqで得られてきた情報の典型例を紹介することにとどめる。

　大脳皮質には、錐体細胞や非錐体細胞などの神経細胞や様々なグリア細胞などが見られ、古くから神経細胞タイプの識別が行われてきた。初期のFluidigm C1を用いたscRNA-seq技術でも、マウス皮質の小規模な細胞数を分類した研究で、これまで知られていた主要な神経細胞タイプとは違うタイプが見つかりscRNA-seqの有効性が示された。その後のドロップレット使用の3’エンドリード法を利用した多数の細胞数の解析で、更に多数の神経細胞のタイプが見つかっている[31][32] [33]。特に、GABA作動性介在神経細胞タイプの多様性とその発生についての、これまでの組織化学的な研究からは得られていなかった多くの情報は重要であろう。また、初期の発生過程、老化の理解が、scRNA-seq技術を利用することで進んでいる。更に、神経活動や臨界期に伴い変化するmRNAも細胞ごとに調査され興味深い 。 ヒトを含めた霊長類の大脳についても発達段階を含めてscRNA-seqが適用されてきている [34][35][36]。ヒトや霊長類に特徴的とされる島のvon Economo神経細胞（紡錘細胞）のような希少な神経細胞のscRNA-seqにも成功している。

　海馬では、これまでの研究で記載されてきた神経細胞のタイプの存在が確認され、更に新規のタイプが見つかった。中枢神経系では、その他、外側膝状体、大脳基底核(足底核) 、視床下部 、線条体、中脳 、手綱、発生中の間脳 、さらに小脳などの結果が報告されてきている。例えば、構成する細胞についての情報が詳細に研究されてきたと思われていたマウスの小脳においても、分子層にこれまでの星状細胞、バスケット細胞というカテゴリーとは違ったギャップジャンクションに特徴を持つ２種類の神経細胞があることが示唆されている[37]。

　脳の外部では、運動神経[38]、感覚神経、らせん神経節 、臭覚神経、腸神経系 [39]　、網膜[40][41][42] [43]でのscRNA-seqデータがある。 またiPS細胞やES細胞由来の神経組織オルガノイドに含まれる神経細胞タイプを知る上でも利用されている。このようなアプローチは、ネアンデルタール人型の遺伝子を持つ脳オルガノイドの解析やSARS-CoV-2に感染する脳オルガノイド中の細胞の同定など、新たな応用例が発表されてきており興味深い。

神経細胞以外の細胞

　上衣細胞は、神経幹細胞としての役割が示唆されてきたが、scRNA-seqによる解析ではその可能性が支持されなかった。グリア細胞では、ラジアルグリア、アストロサイト[44]に多様性があることが示唆されてきている。また、オリゴデンドロサイト[45]については、これまで細胞生物学的に研究されてきた分化の過程がscRNA-seqにより詳細に解析されている。ミクログリアは、神経系の発達、老化、損傷などに伴う重要な遺伝子発現状態の変化がscRNA-seqにより詳細に明らかになった 。また、CNS境界関連マクロファージ(BAM) や脳血管系のscRNA-seqも実施されている。

疾患

　scRNA-seqは、疾患の理解にも有用である。scRNA-seqでは、疾患に伴う遺伝子発現状態の変化を細胞タイプごとに観察することができるので、バルクRNA-seqでは埋もれていた遺伝子発現状態の変化や細胞ごとの変化を検出できるという長所がある。 例えば、筋萎縮性側索硬化症、多発性硬化症 、アルツハイマー病やそのモデル動物 [46]、統合失調症、てんかん、自閉症やレット症候群、シャルコー・マリー・トゥース病、ダウン症[47]、パーキンソン病 、ハンチントン病、がんなどに適用されている。最近、Perturb-seqにより、自閉症に関わる遺伝子の欠損に伴う細胞状態の変化などもscRNA-seqで報告されており、疾患の理解のための新たな実験系の開発も始まりつつある。

scRNA-seqの展望

神経系の多様性と進化

　NGSを用いることで、どんな生物種にも適用可能なscRNA-seqは、既に多様な生物の神経系の細胞の理解、更には種間の相同性や差異の研究に利用されており、神経系の進化を細胞レベルで考察するのに有用であろう（例、線虫、ショウジョウバエ 、カタユウレイボヤCiona intestinalis 、ゼブラフィッシュ、アカミミガメTrachemys scripta elegans、トカゲPogona vitticeps, pv、ニワトリ[48]、霊長類[49]）。ただ、遺伝子やトランスクリプトームの研究が進んでいる生物種では比較的容易であるが、遺伝子のアノテーションが十分でない生物種を用いる場合、scRNA-seqのデータ解析は困難を伴う。また種を超えた細胞タイプの相同性の理解には様々な工夫が必要である[50]。

データベースと統合

　獲得されたscRNA-seqのデータは様々な目的で利用できるので、データベース化し利用できるようにする必要がある。神経系のトランスクリプトーム一般のデータベースが多数公開されており、scRNA-seqのデータも基本的にNCBIのGene Expression Omnibus[51]に登録されている。また、オープンサイエンス推進のためにcommon coordinate framework (CCF) やcentral annotation platform (CAP)という概念のもと、特にscRNA-seqを意識したものとして、米国のBRAIN Initiative Cell Census Consortium 、Human Cell Atlas ProjectのHuman Cell Atlas Data Portal [52]、そのマウス版であるTabula Muris [53]やSten Linnarssonラボのマウス脳発生データベース[54]、アレン脳研究所のAllen Brain Atlas [55]、ブロード研究所のSingle Cell Portal [56]などのデータベースが稼働している。また、異なった方法や実験で得られたscRNA-seqのデータや後述の複数モダリティのシングルセルオミクスのデータを体系的に比較することも重要であり、CCA (Canonical correlation analysis), Seurat 3.0以降に組み込まれたMMN (Mutual Nearest Neighbors)、LIGER 、Harmony 、MetaNeighber、Conos[57]のようなアルゴリズムが開発され、後述の複数モダリティのシングルセルオミクスを組み込んだ統合サイトもでき始めている（[58]）。またデータベースを利用して発現類似性検索も研究されている。

空間トランスクリプトミクス

　多数の細胞を扱うscRNA-seqの弱点は、組織から細胞や細胞核を解離する必要があるので、その細胞が存在していた解剖学的あるいは空間的な位置の情報を消去してしまうということである。組織切片におけるタンパク質などの分布は免疫組織化学、mRNAの分布はin situ hybridizationで検出することができるが、数多くのmRNAの分布を情報処理技術と組み合わせ一気に同定する方法がscRNA-seqと同様に開発されてきている（Slide-seq、HDST、Expansion sequencing[59]など, ）、更に10x Genomics社が市販するVisium（現時点ではシングルセルレベルではない）などがある。現状では、大きな組織の空間トランスクリプトミクスは、空間解像度は限定されており、技術普及の観点からも課題が多い。しかし、そのデータを解析するためのアルゴリズム [60][61] [62]、更にMerFish 、corrFISHのように、subcellularレベルで多数のmRNAを検出する方法が多数開発されてきており（文献： osmFISH 30377364、STARmap (spatially-resolved transcript amplicon readout mapping) 、seqFISH+、pciSeq(probabilistic cell typing by in situ sequencing)[63]、DSP(Digital Spatial Profiling) 、scRNA-seqと組み合わせることで、その弱点を補う空間トランスクリプトミクスにも利用され始め[64]、今後の発展が期待される分野である。

統合解析

　同一の細胞からscRNA-seqの情報だけでなく、ゲノム配列、ATAC-seqなどによるエピゲノム解析、少数のタンパク質、あるいはプロテオームなど複数のモダリティを同時に観察するオミクス（Single-cell multimodal omics）が注目されている。2019年には、Nature Methodsの「Methods of the Year」に選ばれており、現状については、その特集号などを参考にされたい。例えば、細胞表面に提示されているマーカー分子に対する抗体にDNAを付加することで、そのマーカーを発現する細胞のトランスクリプトームを解析するCITE-seq、 REAP-seqも、既知の細胞表面マーカーの発現とscRNA-seqが同時に観察できる方法である。

　複数モダリティのシングルセルオミクスとして、神経科学分野で注目されるのは、scRNA-seqをパッチクランプによる電気生理学的情報と組み合わせたPatch-seq である。また、ゲノムDNAとscRNA-seqを同時に観察することによって、近年、精神疾患の観点から注目されている発生途中で生じる遺伝子変異を研究するPRDD-seqは今後の展開が注目される。最後に、BARseq (barcoded anatomy resolved by sequencing) 、CONNECTID[65]、Epi-Retro-seq[66]のような方法は、コネクトーム（神経細胞の結合性）と遺伝子発現状態を記録できるオミクスの新たな方向として興味深い。

関連項目

• ゲノムワイド関連解析
• ディファレンシャルディスプレイ
• In situハイブリダイゼーション法
• 免疫組織化学法
• エピジェネティクス
• コネクトーム