「シングルセルRNAシーケンシング」の版間の差分

提供:脳科学辞典
ナビゲーションに移動 検索に移動
編集の要約なし
編集の要約なし
9行目: 9行目:


{{box|text=
{{box|text=
シングルセルRNAシーケンシング(single cell RNA sequencing, 以下scRNA-seq)は、次世代シーケンシング (next generation sequencing、以下NGS)技術を使用して個々の細胞が発現しているmRNA全体、つまりトランスクリプトームを質的、量的に網羅的に調べ、細胞ごとの違いを高解像度で検出、分類することで、細胞の分類を行うことができる分子生物学的、コンピュータ生物学的技術である。また、刺激、発生など細胞の状況に応じて、個々の細胞のトランスクリプトームの情報を得ることで、病態や細胞系譜などの解析も可能である。特に多様なニューロンが存在する神経系では、この方法により、神経細胞の分類や状態についての知見が深まり、更に新しいバイオマーカー(biomarker)の同定などが網羅的に行われるようになった。
シングルセルRNAシーケンシング(single cell RNA sequencing, 以下scRNA-seq)は、次世代シーケンシング (next generation sequencing、以下NGS)技術を使用して個々の細胞が発現しているmRNA全体、つまりトランスクリプトームを質的、量的に網羅的に調べ、細胞ごとの違いを高解像度で検出、分類することで、細胞の分類を行うことができる分子生物学的、コンピュータ生物学的技術である。また、刺激、発生など細胞の状況に応じて、個々の細胞のトランスクリプトームの情報を得ることで、病態や細胞系譜などの解析も可能である。特に多様なニューロンが存在する神経系では、この方法により、神経細胞の分類や状態についての知見が深まり、更に新しいバイオマーカー(biomarker)の発見などが網羅的に行われるようになった。。
。}}
。}}


==トランスクリプトーム==
==scRNA-seqとその開発史==
===トランスクリプトーム===
トランスクリプトーム(transcriptome)は、細胞中に存在する全ての転写産物(タンパク質をコードするmRNA、タンパク質をコードしないノンコーディングRNA、マイクロRNAなど)の総体である<ref><pubmed>19015660</pubmed></ref><ref><pubmed>31341269 </pubmed></ref>。トランスクリプトームは、ゲノムとは異なり、同一の個体でも、組織ごとに、更には発生段階や細胞外環境や刺激によって変化する。トランスクリプトームは、同質あるいは異質の多数の細胞集団(組織、培養細胞)からRNA抽出後、cDNAに変換し、それを1990年代に出現したDNAマイクロアレイのように数多くの既知mRNAを識別する技術によって解析されるようになった。その後、NGSの利用により、希少mRNAやノンコーディングRNAを含めた未知の転写産物の高感度検出が可能になるとともに、スプライシングで成熟していく過程のmRNAなど、転写産物の種類だけでなく、転写産物の構造的な違い(スプライシングバリアント、SNPs、変異など)の解析もできるようになった。また、NGSは、ヒトやモデル実験生物(マウス、ゼブラフィッシュ、ショウジョウバエ、センチュウなど)だけでなく、多種多様な生物のトランスクリプトームの把握も可能になった。本稿では、このような多数の細胞集団、つまり個体や特定の組織全体ではなく、細胞1つの持つトランスクリプトームを解析する方法(scRNA-seq)とそのscRNA-seqデータを利用することで得られる知見について概説する。
トランスクリプトーム(transcriptome)は、細胞中に存在する全ての転写産物(タンパク質をコードするmRNA、タンパク質をコードしないノンコーディングRNA、マイクロRNAなど)の総体である<ref><pubmed>19015660</pubmed></ref><ref><pubmed>31341269 </pubmed></ref>。トランスクリプトームは、ゲノムとは異なり、同一の個体でも、組織ごとに、更には発生段階や細胞外環境や刺激によって変化する。トランスクリプトームは、同質あるいは異質の多数の細胞集団(組織、培養細胞)からRNA抽出後、cDNAに変換し、それを1990年代に出現したDNAマイクロアレイのように数多くの既知mRNAを識別する技術によって解析されるようになった。その後、NGSの利用により、希少mRNAやノンコーディングRNAを含めた未知の転写産物の高感度検出が可能になるとともに、スプライシングで成熟していく過程のmRNAなど、転写産物の種類だけでなく、転写産物の構造的な違い(スプライシングバリアント、SNPs、変異など)の解析もできるようになった。また、NGSは、ヒトやモデル実験生物(マウス、ゼブラフィッシュ、ショウジョウバエ、センチュウなど)だけでなく、多種多様な生物のトランスクリプトームの把握も可能になった。本稿では、このような多数の細胞集団、つまり個体や特定の組織全体ではなく、細胞1つの持つトランスクリプトームを解析する方法(scRNA-seq)とそのscRNA-seqデータを利用することで得られる知見について概説する。


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


1つの細胞内にある全RNA(ribosomal RNAを含む)は細胞種にもよるが1-50pgである。そのうち、mRNAの占める割合は1-5%程度である<ref><pubmed>15239941</pubmed></ref>。この微量のmRNAをcDNAに変換してから大幅に増幅できるPCRが発明されることで、1つの細胞が発現するmRNAを高感度で検出できるようになった。例えば、1991年、Linda BuckとRichard Axelは、嗅覚受容体がGタンパク質であると仮定し、個々の嗅覚細胞で特異的に観察されるGタンパク質mRNAを比較することで、嗅覚受容体の同定に成功した<ref><pubmed>1840504</pubmed></ref>(2004年、ノーベル生理学・医学賞)。1995年になると、Catherine DulacとRichard Axelは、異なる鋤鼻神経細胞で特異的に発現する遺伝子を1つの細胞から作製したcDNAライブラリーを比較するディファレンシャル・スクリーニングにより、フェロモン受容体を同定した<ref><pubmed>7585937</pubmed></ref>。同じ手法で異なる種類の神経細胞で発現している遺伝子も同定されており<ref><pubmed>9778248</pubmed></ref><ref><pubmed>12230981</pubmed></ref>、1つの細胞の持つトランスクリプトームを比較するアプローチが神経系で特徴的に発現している遺伝子の同定に原理的に効果的であることを示した。
1つの細胞内にある全RNA(ribosomal RNAを含む)は細胞種にもよるが1-50pgである。そのうち、mRNAの占める割合は1-5%程度である<ref><pubmed>15239941</pubmed></ref>。この微量のmRNAをcDNAに変換してから大幅に増幅できるPCRが発明されることで、1つの細胞が発現するmRNAを高感度で検出できるようになった。例えば、1991年、Linda BuckとRichard Axelは、嗅覚受容体がGタンパク質であると仮定し、個々の嗅覚細胞で特異的に観察されるGタンパク質mRNAを比較することで、嗅覚受容体の同定に成功した<ref><pubmed>1840504</pubmed></ref>(2004年、ノーベル生理学・医学賞)。1995年になると、Catherine DulacとRichard Axelは、異なる鋤鼻神経細胞で特異的に発現する遺伝子を1つの細胞から作製したcDNAライブラリーを比較するディファレンシャル・スクリーニングにより、フェロモン受容体を同定した<ref><pubmed>7585937</pubmed></ref>。同じ手法で異なる種類の神経細胞で発現している遺伝子も同定されており<ref><pubmed>9778248</pubmed></ref><ref><pubmed>12230981</pubmed></ref>、1つの細胞の持つトランスクリプトームを比較するアプローチが神経系で特徴的に発現している遺伝子の同定に原理的に効果的であることを示した。
22行目: 23行目:
一方で多くの種類のmRNAを1細胞レベルで観察する単細胞トランスクリプトームには技術的なブレークスルーが待たれた。第一の問題はPCRなどの増幅に伴うcDNAごとのバイアスなどのアーティファクトが頻繁に観察されること、そしてもう一つの問題は多種類のcDNAを簡便に識別することを可能にする方法の開発であった。これを可能にしたのが、cDNA増幅法の改良とマイクロアレイの利用であった<ref><pubmed>16547197</pubmed></ref>。しかしながら、増幅に伴うアーティファクトの解決は依然として不十分で、また1つの細胞ごとに高価なマイクロアレイを使用することは、多数の細胞のトランスクリプトームを観察するのには限界があった。2009年になると、これらの問題を解決できる可能性として、High-throughput sequencing (HTS)を利用するscRNA-seqプロトコルがAzim Suraniのグループによって報告された<ref><pubmed>19349980</pubmed></ref>。しかしながら、この論文でもたった8個の細胞の解析に留まっており、この方法でも一つの細胞ごとに処理を行うという操作が必要で、多数の細胞についてのトランスクリプームを一挙に理解することはできなかった。
一方で多くの種類のmRNAを1細胞レベルで観察する単細胞トランスクリプトームには技術的なブレークスルーが待たれた。第一の問題はPCRなどの増幅に伴うcDNAごとのバイアスなどのアーティファクトが頻繁に観察されること、そしてもう一つの問題は多種類のcDNAを簡便に識別することを可能にする方法の開発であった。これを可能にしたのが、cDNA増幅法の改良とマイクロアレイの利用であった<ref><pubmed>16547197</pubmed></ref>。しかしながら、増幅に伴うアーティファクトの解決は依然として不十分で、また1つの細胞ごとに高価なマイクロアレイを使用することは、多数の細胞のトランスクリプトームを観察するのには限界があった。2009年になると、これらの問題を解決できる可能性として、High-throughput sequencing (HTS)を利用するscRNA-seqプロトコルがAzim Suraniのグループによって報告された<ref><pubmed>19349980</pubmed></ref>。しかしながら、この論文でもたった8個の細胞の解析に留まっており、この方法でも一つの細胞ごとに処理を行うという操作が必要で、多数の細胞についてのトランスクリプームを一挙に理解することはできなかった。


==scRNA-seqの現在==
===scRNA-seqの現状===
それ以来、完全長cDNAを増幅したり、細胞ごとに異なる分子識別子(unique molecular identifiers: UMI)を持つcDNAを増幅させるscRNA-seqが考案され始め、2013年には、このような1細胞のシーケンシング技術が、Nature Methods誌のMethod of the Year に選ばれた https://www.nature.com/collections/mysbdwgfll。たとえば、SMART-seq(Switch mechanism at the 5' End of RNA Templates)<ref><pubmed>22820318</pubmed></ref>およびその改良されたプロトコルであるSMART-seq2 <ref><pubmed> 24056875 </pubmed></ref> <ref><pubmed>24385147</pubmed></ref>は、完全長cDNA合成のためのプロトコルである。また、MARS-seq(Massively parallel single-cell RNA-seq)<ref><pubmed> 24531970 </pubmed></ref>、STRT(single-cell tagged reverse transcription)<ref><pubmed>21543516</pubmed></ref> <ref><pubmed>24363023</pubmed></ref>、CEL-seq(Cell Expression by Linear amplification and Sequencing)<ref><pubmed>22939981</pubmed></ref>、CEL-seq2<ref><pubmed> 27121950 </pubmed></ref>などが報告されてきた。特にSMART-seq(SMART-seq2)は、微小管によるマニュアル捕獲、限界希釈、セルソーター、レーザー捕獲法などを用いる多穴プレート法、更に半導体集積回路製作技術で作った流体集積回路を利用するFluidigm C1の装置https://jp.fluidigm.comと組み合わせることで利用される機会が多い。このプロトコールの特徴は、全長トランスクリプトームを得ることができることであり、mRNAのスプライシングバリアントなどのアイソフォーム、SNPs、変異の検出にも利用できる点で次に説明するUMIを用いる方法に比べて利点があるが、そのコストと処理可能な細胞数の点で弱点がある。
それ以来、完全長cDNAを増幅したり、細胞ごとに異なる分子識別子(unique molecular identifiers: UMI)を持つcDNAを増幅させるscRNA-seqが考案され始め、2013年には、このような1細胞のシーケンシング技術が、Nature Methods誌のMethod of the Year に選ばれた[https://www.nature.com/collections/mysbdwgfll]。たとえば、SMART-seq(Switch mechanism at the 5' End of RNA Templates)<ref><pubmed>22820318</pubmed></ref>およびその改良されたプロトコルであるSMART-seq2 <ref><pubmed> 24056875 </pubmed></ref> <ref><pubmed>24385147</pubmed></ref>は、完全長cDNA合成のためのプロトコルである。また、MARS-seq(Massively parallel single-cell RNA-seq)<ref><pubmed> 24531970 </pubmed></ref>、STRT(single-cell tagged reverse transcription)<ref><pubmed>21543516</pubmed></ref> <ref><pubmed>24363023</pubmed></ref>、CEL-seq(Cell Expression by Linear amplification and Sequencing)<ref><pubmed>22939981</pubmed></ref>、CEL-seq2<ref><pubmed> 27121950 </pubmed></ref>、Seq-Well <ref><pubmed> 28192419</pubmed></ref>
などが報告されてきた。最近になって、sci-RNA-seq (single-cell combinatorial indexing RNA sequencing) <ref><pubmed> 28818938 </pubmed></ref>, SPLiT-seq(split-pool ligation-based transcriptome sequencing)<ref><pubmed>29545511</pubmed></ref>のように特殊な機器を利用せずに細胞特異的UMIを保持するcDNAを作製する方法も報告されている。


しかしながら、更に、重要なscRNA-seqの方法論についての進歩は、2015年、Harvard Medical Schoolの独立した2つのグループが、inDropそしてDrop-seqという類似した2つの高スループットな方法を開発したことであろう<ref><pubmed>26000487</pubmed></ref> <ref><pubmed>26000488</pubmed></ref>。これらの方法では、マイクロ流体力学 (Microfluidics) 、 UMIとしてDNAバーコーディング (DNA barcoding) 、そしてNGSを利用することで、自動化とサンプル調製の容易さから、1つの細胞あたりに要するコストを大幅に低下させることに成功した(Drop-seqは発表時で、1細胞あたり6セント)。つまり、細胞1つずつをマイクロ流体力学によるエマルジョン技術を利用した装置に流入させ、その1細胞を試薬を封入した1つのDroplet(油滴)に自動的に閉じ込める。そのDroplet中には、DropletごとにUMIとして異なったDNAバーコードを持つゲルビーズが入っており、それを足場にcDNA合成反応を実施することで、同じ細胞に含まれていたmRNAが同じUMIを持つcDNAとして合成され、Dropletを破壊した後も、そのmRNA/cDNAが由来した細胞を識別できるということを利用している(図1)。このようにして増幅したUMI付きcDNAをNGSで配列決定することによりscRNA-seqが可能になる。なお、DropSeqはコストが低いが、細胞の取得率と検出感度が低い弱点がある。inDropはDropSeqより細胞取得率が高く、パラメータを調整することで低レベルで発現される遺伝子の検出に有利である。
これらの方法のうち、SMART-seq、その改良法であるSMART-seq2は、微小ピペットによるマニュアル捕獲、セルソーター、レーザー捕獲法などを用いる多穴プレート法、更に半導体集積回路製作技術で作った流体集積回路を利用するFluidigm C1の装置[https://jp.fluidigm.com]と組み合わせることで利用される機会が多い。このSMART-seq2プロトコールの特徴は、mRNAの全領域を読むことで、全長トランスクリプトームを得ることができることであり、mRNAのスプライシングバリアントなどのアイソフォーム、アリルごとの発現情報が得られるSNPs、変異の検出にも利用できる。また、それぞれ細胞ごとの反応を独立した場所で行うため、別の細胞の反応と混じる可能性がない。これらの点が、次に説明するDropletを使用して3’末端のみを標的にしたscRNA-seqに比べた場合の長所であるが、その高コストと処理可能な細胞数の少なさが短所である。


inDropの方法は、1 Cellbio社から販売されているhttps://1cell-bio.com。しかし、特に重要なのは10xGenomics社が同様の原理を用いた「Chromium」と命名された機器と試薬のシステムを市販することで、多くの研究者に利用できることになったことであるhttps://www.10xgenomics.com/jp/。Svenssonらによる最近のデータベースhttps://www.biorxiv.org/content/10.1101/742304v2, https://www.nxn.se/single-cell-studies/guiでは、scRNA-seqを用いた論文で用いられた方法について網羅的に調査しているが、この数年、10xGenomics社のChromiumを用いた方法が飛躍的に増加し、ほぼ寡占状態になりつつあることがわかる(現在、10XGenomics社とBioRad社の間で関連特許をめぐる係争がある。)。このシステムは市販であるので導入が容易であり、inDropやDropSeqに比べ、多くの転写産物の高感度検出が可能であるが、ランニングコストが高い。
===Droplet使用の3’エンドリード法===
しかしながら、もっとも重要なscRNA-seqの方法論についての進歩は、2015年、Harvard Medical Schoolの独立した2つのグループが、inDropそしてDrop-seqという類似した2つの高スループットな方法を開発したことであろう<ref><pubmed>26000487</pubmed></ref> <ref><pubmed>26000488 </pubmed></ref>。これらの方法では、マイクロ流体力学 (Microfluidics) 、 Cell BarcodeとUMIとしてDNAバーコーディング (DNA barcoding) 、そしてNGSを利用することで、自動化とサンプル調製の容易さから、1つの細胞あたりに要するコストを大幅に低下させることに成功した(Drop-seqは発表時で、1細胞あたり6セント)。つまり、細胞1つずつをマイクロ流体力学によるエマルジョン技術を利用した装置に流入させ、その1細胞を1つのDroplet(油中水滴)に自動的に閉じ込める。そのDroplet中には、DropletごとにCell barcode/UMIとして異なったDNAバーコードを持つゲルビーズ(Gel Beads in Emulsion, GEMs)が入っており、それを足場に3’末端のみを標的にしたcDNA合成反応を実施することで、同じ細胞に含まれていたmRNAが同じUMIを持つcDNAとして合成され、そのmRNA/cDNAが由来した細胞を識別できるということを利用している(図1)。このようにして3’末端のみを増幅したバーコード付きcDNAをNGSで配列決定することによりscRNA-seqが可能になる。なお、DropSeqはコストが低いが、細胞の取得率と検出感度が低い弱点がある。inDropはDropSeqより細胞取得率が高く、パラメータを調整することで低レベルで発現される遺伝子の検出にも有利である<ref><pubmed>30472192</pubmed></ref>。
DropSeqのセットアップはDolomite Bio ([https://www.dolomite-bio.com])、inDropは1 Cellbio社から販売されている[https://1cell-bio.com]。しかし、特に重要なのは10xGenomics社が同様の原理を用いた「Chromium」と命名された機器と試薬のシステムを市販することで、多くの研究者に利用できることになったことである[https://www.10xgenomics.com/jp/]。Svenssonらによる最近のデータベース[https://www.biorxiv.org/content/10.1101/742304v2], [http://www.nxn.se/single-cell-studies/gui]では、scRNA-seqを用いた論文で用いられた方法について調査しているが、この数年、10xGenomics社Chromiumを用いた論文が飛躍的に増加し、scRNA-seqの方法として、最も一般的になりつつあることがわかる(現在、10XGenomics社とBioRad社の間で関連特許をめぐる係争がある。)。このシステムは市販であるので導入が容易であり、DropSeqやinDropに比べ、多くの転写産物の高感度検出が可能であるが、ランニングコストは高い<ref><pubmed>30472192</pubmed></ref>。


==scRNA-seqの実際==
==scRNA-seqの実際==
ここでは主流になっている10xGenomics社のChromiumを用いた方法とSMART-seqなどを用いた方法に共通する方法の実際について議論する。シングルセルRNAシーケンシングの利用には、4つのステップがある。1)細胞をバラバラに単離すること。2)ライブラリーの作製とNGSシーケンシング。3)前処理(preprocessing、得られた配列の整理)。4)ダウンストリーム分析(生物学的な情報を得る)。これらのうち、2)の段階については、上に記述したように市販の機器や試薬を利用することが多くなっているので、各社のマニュアル等を参考にするのが現実的である。
ここでは主流になっている10xGenomics社のChromiumを用いた方法とSMART-seqなどを用いた方法に共通する方法の実際について俯瞰する。scRNA-seqの利用には、4つのステップがある(図2)<ref><pubmed>30089861</pubmed></ref>。1)個体や組織を採集し、そこから細胞あるいは細胞核を個別にすること。2)ChromiumやSMART-seq2などによる個々の細胞からのライブラリーの作製とNGSシーケンシング。3)前処理(preprocessing、得られた配列の整理)。4)ダウンストリーム解析(生物学的な情報を得る)。これらのうち、2)の段階については、上に記述したように市販の機器や試薬を利用する機会が多くなっているので、そのためのマニュアル等を参考にするのが現実的である。


===組織からの細胞の分離===
===組織からの細胞、細胞核の分離===
血液細胞のように浮遊した細胞ではない場合、物理的あるいは酵素処理などによって、生組織から状態の良い細胞をdissociationする必要がある。神経系組織の酵素処理には、パパインを用いる方法が広く用いられている。ただ、しばしば問題となるのが、酵素処理のため短時間加温することで、発現量が変化する遺伝子が存在することである。例えば、脳のミクログリアの解析には、低温下で組織をホモゲナイズするなどの工夫が必要であった<ref><pubmed>30471926</pubmed></ref>。また、酵素処理時に転写阻害剤であるアクチノマイシンで処理することで、このような現象を抑制できる<ref><pubmed>29024657</pubmed></ref>。更に、ヒマラヤ氷河から得られた細菌Bacillus licheniformisから得られた低温プロテアーゼを用いる方法も報告されている<ref><pubmed>28851704</pubmed></ref>。なお、ヒト組織などから生細胞を得ることは困難なことが多い。この場合、凍結した組織から、各細胞由来の核を調製し、核内のmRNAを分析するアプローチもあり、細胞質を持つ生細胞を利用した場合より感度は劣るが、目的によっては利用可能である<ref><pubmed>30586455</pubmed></ref>。
血液細胞のように浮遊した細胞ではない場合、物理的あるいは酵素処理などによって、生組織から状態の良い細胞をdissociationする必要がある。神経系組織の酵素処理には、パパインを用いる方法が広く用いられている。ただ、しばしば問題となるのが、酵素処理のため短時間加温することで、発現量が変化する遺伝子が存在することである<ref><pubmed>27090946</pubmed></ref>。例えば、脳のミクログリアの解析には、低温下で組織をホモゲナイズするなどの工夫が必要であった<ref><pubmed>30471926</pubmed></ref>。また、このような現象を抑制するために、酵素処理時に転写阻害剤であるアクチノマイシンで処理したり<ref><pubmed>29024657</pubmed></ref>、ヒマラヤ氷河から得られた細菌Bacillus licheniformisから得られた低温プロテアーゼを用いる方法も報告されている<ref><pubmed>28851704</pubmed></ref>。また、細胞解離後に、メタノールで固定しscRNA-seqに使用することも可能である<ref><pubmed>28526029</pubmed></ref>。
単離した細胞は、そのまま10xGenomicsのChromiumのプラットフォームに導入することができるが、抗体などを用いたFACS、パニング、磁気ビーズカラムなどによる特定のマーカーを細胞表面などに発現する細胞の単離を行う場合もある。
なお、ヒト組織などから生きた細胞を得ることは困難なことが多い。この場合、凍結した組織から、各細胞由来の核を調製し、核内のmRNAを分析する方法 snRNA-seq (single-nucleus RNA-seq)があり、細胞質を持つ生細胞を利用した場合より感度は劣るが、細胞の同定などの目的には十分使用に耐える<ref><pubmed>27471252</pubmed></ref><<<ref><pubmed>28846088</pubmed></ref><<ref><pubmed>29220646</pubmed></ref><ref><pubmed>28846088</pubmed></ref><ref><pubmed>30586455</pubmed></ref><ref><pubmed>28729663</pubmed></ref>[https://www.biorxiv.org/content/10.1101/630087v1]。snRNA-seqでは、組織をそのまま凍結することから開始するので、上述したscRNA-seqの問題である酵素処理や加温などを避けることができる。こうしたプロトコールの一部は、protocols.ioのHuman Cell Atlasのグループ[https://www.protocols.io/groups/hca]で公開されている。
更に、RNAを分析するscRNA-seqではないが、シングルセルgenome-seqの変法として、シングルセルの遺伝子発現を推定する方法として、トランスポゾンを用いることでゲノムのオープンクロマチン領域を選択的に検出し、ライブラリーを作製しシーケンスするsingle cell ATAC-seq (Assay for Transposase-Accessible Chromatin)がある<ref><pubmed>26083756</pubmed></ref>


単離した細胞は、そのまま10xGenomicsのChromiumのプラットフォームに導入することができるが、抗体などを用いたFACS、パニング、MACS(磁気ビーズカラム)などによる特定のマーカーを細胞表面などに発現する細胞の単離を行う場合もある。


更に、RNAを分析するscRNA-seqではないが、single cell genome-seqの変法として、シングルセルの遺伝子発現を推定する方法として、トランスポゾンを用いることでゲノムのオープンクロマチン領域を選択的に検出し、ライブラリーを作製しシーケンスするsingle cell ATAC-seq (Assay for Transposase-Accessible Chromatin)<ref><pubmed>26083756</pubmed></ref>, <ref><pubmed>29434377</pubmed></ref><ref><pubmed>25953818</pubmed></ref>, single cell THS-seq (transposome hypersensitive-site) <ref><pubmed>29227469</pubmed></ref>がある。またDNAメチル化をシングルセルレベルで観察する方法も報告されている<ref><pubmed>28798132</pubmed></ref>。
===scRNA-seqデータの前処理===
===scRNA-seqデータの前処理===
血液細胞のように浮遊した細胞ではない場合、物理的あるいは酵素処理などによって、生組織から状態の良い細胞をdissociationする必要がある。神経系組織の酵素処理には、パパインを用いる方法が広く用いられている。ただ、しばしば問題となるのが、酵素処理のため短時間加温することで、発現量が変化する遺伝子が存在することである。例えば、脳のミクログリアの解析には、低温下で組織をホモゲナイズするなどの工夫が必要であった<ref><pubmed>30471926</pubmed></ref>。また、酵素処理時に転写阻害剤であるアクチノマイシンで処理することで、このような現象を抑制できる<ref><pubmed>29024657</pubmed></ref>。更に、ヒマラヤ氷河から得られた細菌Bacillus licheniformisから得られた低温プロテアーゼを用いる方法も報告されている<ref><pubmed>28851704</pubmed></ref>。なお、ヒト組織などから生細胞を得ることは困難なことが多い。この場合、凍結した組織から、各細胞由来の核を調製し、核内のmRNAを分析するアプローチもあり、細胞質を持つ生細胞を利用した場合より感度は劣るが、目的によっては利用可能である<ref><pubmed>30586455</pubmed></ref>
10x Genomics社のChromium、IlluminaのNGSを利用した場合、Cell Ranger(Linux上で作動)を用いて、各生物種ごとのレファレンス配列リスト(https://www.ncbi.nlm.nih.gov/grc)やEggNOG ([http://eggnogdb.embl.de])などを利用し、細胞とトランスクリプトーム(各遺伝子の発現)の対応マトリックスを作製する。その後のデータの処理についても、10xGenomics社がソフトウェアLoupeを提供している。しかしながら、その後のダウンストリーム解析を考慮して、R言語, Python, MATLABなどのデータ解析のための汎用プログラミング言語やコードで扱えるオブジェクトに変換するのが通常である。ここでは、scRNA-seqデータ解析のために最もよく利用されているR言語を用いたパッケージ「Seurat」<ref><pubmed> 29608179 </pubmed></ref> <ref><pubmed> 31178118 </pubmed></ref>を中心に紹介したい。なお、Pythonを利用したものでは、ドイツ・ミュンヘンInstitute of Computational Biologyの Fabian J. Theisらが開発しているScanpyが有名である<ref><pubmed> 29409532</pubmed></ref>。
単離した細胞は、そのまま10xGenomicsのChromiumのプラットフォームに導入することができるが、抗体などを用いたFACS、パニング、磁気ビーズカラムなどによる特定のマーカーを細胞表面などに発現する細胞の単離を行う場合もある。
更に、RNAを分析するscRNA-seqではないが、シングルセルgenome-seqの変法として、シングルセルの遺伝子発現を推定する方法として、トランスポゾンを用いることでゲノムのオープンクロマチン領域を選択的に検出し、ライブラリーを作製しシーケンスするsingle cell ATAC-seq (Assay for Transposase-Accessible Chromatin)がある<ref><pubmed>26083756</pubmed></ref>。


New York UniversityのRahul Satija研究室が開発しているSeuratは、scRNA-seq解析のために広く利用されているR言語のパッケージであり、2019年8月現在、その最新バージョンはSeurat 3.1である。論文の正式発表前から、サポート情報提供やコード修正なども頻繁に行っており、Satija研究室のWeb site(satijalab.org/Seurat, [http://satijalab.org/Seurat])、Github([https://github.com/satijalab/Seurat])、更にTwitterアカウント(@satijalab)などで最新情報を得ることできる。


===ダウンストリーム解析===
最初に行うのは、scRNA-seqデータの品質管理である。ここでは、質の低い細胞のデータ(転写産物の種類が少ない、ミトコンドリア由来の転写産物が多い)を取り除く。また、複数の試料を組み合わせる場合には、バッチごとの違いについて検討する。特にDropletを使用するscRNA-seqの多くのケースで問題になるのが、Dropletに2つ以上の細胞が封じ込められ、それらが同じCell barcodeを持つために生じるアーティファクトである。通常Doubletsと呼ばれるこの問題はダウンストリーム解析を混乱させるので、細胞単離の段階から注意する必要があるが、明確なマーカー遺伝子が知られていればscRNA-seqデータ取得後にある程度のデータ処理で検討することは可能である。
クラスタリングにより、異なった細胞集団の存在が認識されると、それぞれのクラスターに特徴的に発現している遺伝子を具体的に探索し、細胞集団の持つバイオマーカーによって、そのクラスターの同定が可能になる。例えば、既にニューロンとグリア細胞に特異的に発現する典型的マーカーはよく知られており、それぞれのクラスターの識別は容易である。更に、ニューロンのタイプごとに区別されるマーカーは、In situ hybridizationや免疫組織化学などにより確認できる。このようなクラスターごとに発現が異なる遺伝子(差次的発現遺伝子)を見つけるためには、Seuratのコードや他のダウンストリーム解析のためのコード(MAST <ref><pubmed>26653891</pubmed></ref>など)を用いることができる。scRNA-seqの解析に必要な様々なコードは、scRNA-tools <ref>https://www.scrna-tools.org</ref>やBioconductor<ref> https://www.bioconductor.org</ref>で紹介されており、ダウンロードも可能である。また、最新の情報については、bioRxivなどのプレプリントサーバで公開されていることが多い。細胞ごとの差次的発現遺伝子の可視化には、ドットプロットやヴァイオリンプロットなどが頻繁に用いられる(図4)。
 
また、scRNA-seqでしばしば得られる情報には、細胞の分類とそれらの類縁関係、状態、それらのバイオマーカー遺伝子だけでなく、病態や発生途上の細胞系譜などの動態がある。これらの分析のためには、成分分析(compositional analysis)や系譜干渉(Trajectory interference)の解析手法が用いられる。
このようなノーマライゼーションの過程を経て、scRNA-seqのデータ解析において、最初に行うのが、次元圧縮 (dimensionality reduction)である<ref><pubmed> 30617341</pubmed></ref>
https://github.com/dynverse/dynmethods
。PCA (Principal component analysis, 主成分分析)、UMAP(Uniform Manifold Approximation and Projection, 均一マニフォールド近似と投影)、Diffusion maps,  t-SNE(t-distributed Stochastic Neighbor Embedding , t分布型確率的近傍埋込み)などの手法が用いられる。 特に、t-SNE[http://www.jmlr.org/papers/v9/vandermaaten08a.html](ティースニーと読むのが通常)は、高次元データを低次元の点の集合として可視化することで、それぞれの細胞の持つトランスクリプトームの類似度についての直観的な表示が可能でありしばしば用いられる(図3)。次に、Louvainアルゴリズムなどでクラスタリング(コミュニティ分割)を行い、tSNEグラフ上に表示できる。こうして、違ったタイプの細胞の集合が別のクラスターとして表示される。
また発現している遺伝子側から見た場合、制御ネットワークやパスウェイ解析といったシステム生物学で用いられてきた手法も用いられる。
 
==ダウンストリーム解析==
===細胞クラスターの解釈とマーカー遺伝子候補の発見===
scRNA-seqデータから得られる生物学的知見には、内在的に存在する細胞の種類、外部刺激や環境で変化した細胞の状態、そして種類や変化により特徴的に発現するマーカー遺伝子候補の発見がある。クラスタリングにより、異なった細胞集団の存在が認識されると、それぞれのクラスターに特徴的に発現している遺伝子を具体的に探索し、細胞集団の持つバイオマーカーによって、そのクラスターの同定が可能になる。例えば、既にニューロンとグリア細胞に特異的に発現する典型的マーカーはよく知られており、それぞれのクラスターの識別は容易である。更に、ニューロンのタイプごとに区別されるマーカーや神経活動により変化したニューロンの状態は、In situ hybridizationや免疫組織化学などにより確認できる。このようなクラスターごとに発現が異なる遺伝子(差次的発現遺伝子)を見つけるためには(Differential expression analysis, DE analysis)、SeuratのFindMarkersコマンドでも利用可能である目的別の解析のための専用コード(MAST <ref><pubmed>26653891</pubmed></ref>、DESeq2 <ref><pubmed>25516281</pubmed></ref>
など)を用いることができる。scRNA-seqの解析に必要なコードは、scRNA-tools [https://www.scrna-tools.org], Awesome single cell [https://github.com/seandavi/awesome-single-cell], Bioconductor[https://www.bioconductor.org]で紹介されており、ほとんどがダウンロード可能である。また、最新の情報については、bioRxivなどのプレプリントサーバで公開されていることが多く、scRNA-seqのデータ(下記参考)とともに、オープンサイエンス実践の好例となっている。細胞ごとの差次的発現遺伝子の可視化には、ドットプロットやヴァイオリンプロットなどが頻繁に用いられる(図4)。
 
===偽時系列解析、制御ネットワーク、パスウェイ解析===
実験的なノイズとは別に生物学的に意味のある遺伝子発現の変動には、位置情報、細胞周期、概日リズム、破裂型プロモーターの作動などの理由で変動が見られるものもある<ref><pubmed> 31217225 </pubmed></ref><ref><pubmed> 26000846</pubmed></ref>。特に、刺激・薬剤処理やさまざまな病態の進行や治療に伴う細胞の変化、発生途上の細胞系譜や細胞分化といった細胞の遷移状態の解析(偽時系列解析Pseudo-time analysis )には、scRNA-seqデータを用いることが極めて効果的である<ref><pubmed>29576429</pubmed></ref><ref><pubmed>28813177</pubmed></ref><ref><pubmed>29565398</pubmed></ref>
。これらの分析のためには軌道推定(Trajectory inference)の解析手法が用いられる。しばしば用いられるMonocle3 <ref><pubmed>30787437</pubmed></ref>など、多くのコードを収集しているGithubのサイトがある [https://github.com/dynverse/dynmethods][https://github.com/agitter/single-cell-pseudotime]。RNA velocityといった、転写産物のスプライシングの状態から細胞の分化状態を推定する方法もある<ref><pubmed>30089906</pubmed></ref>。しかし、これらの方法は、あくまで発生途上の細胞系譜や細胞分化の推定に過ぎない。細胞系譜を更に確実に観察しつつ、scRNA-seqを行うことで、細胞タイプの系統関係を調べる方法として、CRISPR-Cas9を用いたゲノム編集による記録法を導入したscGESTALT<ref><pubmed>29608178</pubmed></ref>、ScarTrace<ref><pubmed>29590089</pubmed></ref> 、LINNAEUS<ref><pubmed>29644996</pubmed></ref>がある。
 
また変化している遺伝子発現を解釈するために、制御ネットワーク(例、SCENIC<ref><pubmed>28991892</pubmed></ref>, [https://github.com/aertslab/SCENIC])やパスウェイ解析(例、Metascape<ref><pubmed>30944313</pubmed></ref>, [http://metascape.org])といったシステム生物学で用いられてきた手法も適用できる。


==神経科学への応用==
==神経科学への応用==
===ニューロンのクラスとタイプ===
===ニューロンのクラスとタイプ===
様々な神経・精神疾患について理解しその診断や治療に役立てるためには、ニューロン、グリア細胞を中心にした神経系にある細胞の「タイプ」を識別し、それぞれの細胞における分子的な変化を観察することが重要である。近年、中枢神経系のグリア細胞にも、多様なアストロサイト、オリゴデンドロサイト、ミクログリアの存在が報告されてきている。一方で、ニューロンは著しく多様であり、このニューロンの多様性こそが、神経系を特徴づけており、その多彩で複雑な機能の発現に必須であることは疑う余地がない。
様々な神経・精神疾患について理解しその診断や治療に役立てるためには、ニューロン、グリア細胞を中心にした神経系にある細胞の「タイプ」を識別し、それぞれの細胞における分子的な変化を観察することが重要である <ref><pubmed>28775344</pubmed></ref><ref><pubmed>29738987</pubmed></ref>。近年、中枢神経系のグリア細胞にも、多様なアストロサイト、オリゴデンドロサイト、ミクログリアの存在が報告されてきている。一方で、ニューロンは著しく多様であり、このニューロンの多様性こそが、神経系を特徴づけており、その多彩で複雑な機能の発現に必須であることは疑う余地がない。
解剖学的な視点から言えば、すべてのニューロンの存在する位置は異なるので、すべてのニューロンは異なるという見方もできる。しかし、これは極論であり、従来の神経科学では、ニューロンの多様性は、それぞれのニューロンの解剖学的な位置、発現している分子、電気生理学、結合性、形態、神経伝達物質、神経伝達物質受容体とシグナル伝達によって識別されてきた。こうしたニューロンの多様性を便宜的に記述するのに、タイプ(type)、クラス(class)、サブクラス(subclass)、サブタイプ(subtype) というような用語が用いられてきた。しかし、本稿では混乱を防ぐため、Masland(2004)の提唱に従い、「クラス」と「タイプ」という単語を用いることとする。タイプは、これ以上分類することができないとされる階層である。例えば、大脳皮質の錐体細胞、網膜神経節細胞といった大雑把な識別は「クラス」と呼ぶ。大脳皮質の錐体細胞というクラスは、層や領野によって「タイプ」が異なるし、網膜神経節細胞には視覚情報によって応答が異なる「タイプ」が存在する。この分類は、免疫組織化学、形態、電気生理学などの技術により識別可能である暫定的なものに過ぎない。本稿で解説するscRNA-seqの技術は、その網羅性からそれぞれのニューロンについてこれまでにないビッグデータを提供することで、このニューロンのタイプの理解に確実な根拠を与えつつある。
解剖学的な視点から言えば、すべてのニューロンの存在する位置は異なるので、すべてのニューロンは異なるという見方もできる。しかし、これは極論であり、従来の神経科学では、ニューロンの多様性は、それぞれのニューロンの解剖学的な位置、発現している分子、電気生理学、結合性、形態、神経伝達物質、神経伝達物質受容体とシグナル伝達によって識別されてきた。こうしたニューロンの多様性を便宜的に記述するのに、タイプ(type)、クラス(class)、サブクラス(subclass)、サブタイプ(subtype) というような用語が用いられてきた。しかし、本稿では混乱を防ぐため、Masland(2004)<ref><pubmed>15242626</pubmed></ref>が提唱し、広く受けいれられている「クラス」と「タイプ」という単語を用いることとする。タイプは、これ以上分類することができないとされる階層である。例えば、大脳皮質の錐体細胞、網膜神経節細胞といった大雑把な識別は「クラス」と呼ぶ。大脳皮質の錐体細胞というクラスは、層や領野によって「タイプ」が異なるし、網膜神経節細胞には視覚情報によって応答が異なる「タイプ」が存在する。この分類は、免疫組織化学、形態、電気生理学などの技術により識別可能である暫定的なものに過ぎない。本稿で解説するscRNA-seqの技術は、その網羅性からそれぞれのニューロンについてこれまでにないビッグデータを提供することで、このニューロンのタイプの理解に確実な根拠を与えつつある。


===大脳===
===大脳===
63行目: 78行目:
===網膜===
===網膜===


===展望===
==展望==
データベース
===進化と他の動物種===
===データベース===
Human Cell Atlas
Human Cell Atlas
Human Brain Transcriptome project
Human Brain Transcriptome project
Single cell portal
Single cell portal
Allen Brain Atlas
Allen Brain Atlas
統合 LIGER, MetaNeighbor
===空間トランスクリプミクス===


空間トランスクリプミクス
===マルチモーダルなオミクス===


== 関連項目 ==
== 関連項目 ==

2020年1月8日 (水) 07:26時点における版

山形方人
Harvard University
DOI:10.14931/bsd.8038 原稿受付日:年月日 原稿完成日:年月日
担当編集委員:

英:single cell RNA sequencing, scRNA-seq

シングルセルRNAシーケンシング(single cell RNA sequencing, 以下scRNA-seq)は、次世代シーケンシング (next generation sequencing、以下NGS)技術を使用して個々の細胞が発現しているmRNA全体、つまりトランスクリプトームを質的、量的に網羅的に調べ、細胞ごとの違いを高解像度で検出、分類することで、細胞の分類を行うことができる分子生物学的、コンピュータ生物学的技術である。また、刺激、発生など細胞の状況に応じて、個々の細胞のトランスクリプトームの情報を得ることで、病態や細胞系譜などの解析も可能である。特に多様なニューロンが存在する神経系では、この方法により、神経細胞の分類や状態についての知見が深まり、更に新しいバイオマーカー(biomarker)の発見などが網羅的に行われるようになった。。 。

scRNA-seqとその開発史

トランスクリプトーム

トランスクリプトーム(transcriptome)は、細胞中に存在する全ての転写産物(タンパク質をコードするmRNA、タンパク質をコードしないノンコーディングRNA、マイクロRNAなど)の総体である[1][2]。トランスクリプトームは、ゲノムとは異なり、同一の個体でも、組織ごとに、更には発生段階や細胞外環境や刺激によって変化する。トランスクリプトームは、同質あるいは異質の多数の細胞集団(組織、培養細胞)からRNA抽出後、cDNAに変換し、それを1990年代に出現したDNAマイクロアレイのように数多くの既知mRNAを識別する技術によって解析されるようになった。その後、NGSの利用により、希少mRNAやノンコーディングRNAを含めた未知の転写産物の高感度検出が可能になるとともに、スプライシングで成熟していく過程のmRNAなど、転写産物の種類だけでなく、転写産物の構造的な違い(スプライシングバリアント、SNPs、変異など)の解析もできるようになった。また、NGSは、ヒトやモデル実験生物(マウス、ゼブラフィッシュ、ショウジョウバエ、センチュウなど)だけでなく、多種多様な生物のトランスクリプトームの把握も可能になった。本稿では、このような多数の細胞集団、つまり個体や特定の組織全体ではなく、細胞1つの持つトランスクリプトームを解析する方法(scRNA-seq)とそのscRNA-seqデータを利用することで得られる知見について概説する。

scRNA-seqの背景

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

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

一方で多くの種類のmRNAを1細胞レベルで観察する単細胞トランスクリプトームには技術的なブレークスルーが待たれた。第一の問題はPCRなどの増幅に伴うcDNAごとのバイアスなどのアーティファクトが頻繁に観察されること、そしてもう一つの問題は多種類のcDNAを簡便に識別することを可能にする方法の開発であった。これを可能にしたのが、cDNA増幅法の改良とマイクロアレイの利用であった[10]。しかしながら、増幅に伴うアーティファクトの解決は依然として不十分で、また1つの細胞ごとに高価なマイクロアレイを使用することは、多数の細胞のトランスクリプトームを観察するのには限界があった。2009年になると、これらの問題を解決できる可能性として、High-throughput sequencing (HTS)を利用するscRNA-seqプロトコルがAzim Suraniのグループによって報告された[11]。しかしながら、この論文でもたった8個の細胞の解析に留まっており、この方法でも一つの細胞ごとに処理を行うという操作が必要で、多数の細胞についてのトランスクリプームを一挙に理解することはできなかった。

scRNA-seqの現状

それ以来、完全長cDNAを増幅したり、細胞ごとに異なる分子識別子(unique molecular identifiers: UMI)を持つcDNAを増幅させるscRNA-seqが考案され始め、2013年には、このような1細胞のシーケンシング技術が、Nature Methods誌のMethod of the Year に選ばれた[1]。たとえば、SMART-seq(Switch mechanism at the 5' End of RNA Templates)[12]およびその改良されたプロトコルであるSMART-seq2 [13] [14]は、完全長cDNA合成のためのプロトコルである。また、MARS-seq(Massively parallel single-cell RNA-seq)[15]、STRT(single-cell tagged reverse transcription)[16] [17]、CEL-seq(Cell Expression by Linear amplification and Sequencing)[18]、CEL-seq2[19]、Seq-Well [20] などが報告されてきた。最近になって、sci-RNA-seq (single-cell combinatorial indexing RNA sequencing) [21], SPLiT-seq(split-pool ligation-based transcriptome sequencing)[22]のように特殊な機器を利用せずに細胞特異的UMIを保持するcDNAを作製する方法も報告されている。

これらの方法のうち、SMART-seq、その改良法であるSMART-seq2は、微小ピペットによるマニュアル捕獲、セルソーター、レーザー捕獲法などを用いる多穴プレート法、更に半導体集積回路製作技術で作った流体集積回路を利用するFluidigm C1の装置[2]と組み合わせることで利用される機会が多い。このSMART-seq2プロトコールの特徴は、mRNAの全領域を読むことで、全長トランスクリプトームを得ることができることであり、mRNAのスプライシングバリアントなどのアイソフォーム、アリルごとの発現情報が得られるSNPs、変異の検出にも利用できる。また、それぞれ細胞ごとの反応を独立した場所で行うため、別の細胞の反応と混じる可能性がない。これらの点が、次に説明するDropletを使用して3’末端のみを標的にしたscRNA-seqに比べた場合の長所であるが、その高コストと処理可能な細胞数の少なさが短所である。

Droplet使用の3’エンドリード法

しかしながら、もっとも重要なscRNA-seqの方法論についての進歩は、2015年、Harvard Medical Schoolの独立した2つのグループが、inDropそしてDrop-seqという類似した2つの高スループットな方法を開発したことであろう[23] [24]。これらの方法では、マイクロ流体力学 (Microfluidics) 、 Cell BarcodeとUMIとしてDNAバーコーディング (DNA barcoding) 、そしてNGSを利用することで、自動化とサンプル調製の容易さから、1つの細胞あたりに要するコストを大幅に低下させることに成功した(Drop-seqは発表時で、1細胞あたり6セント)。つまり、細胞1つずつをマイクロ流体力学によるエマルジョン技術を利用した装置に流入させ、その1細胞を1つのDroplet(油中水滴)に自動的に閉じ込める。そのDroplet中には、DropletごとにCell barcode/UMIとして異なったDNAバーコードを持つゲルビーズ(Gel Beads in Emulsion, GEMs)が入っており、それを足場に3’末端のみを標的にしたcDNA合成反応を実施することで、同じ細胞に含まれていたmRNAが同じUMIを持つcDNAとして合成され、そのmRNA/cDNAが由来した細胞を識別できるということを利用している(図1)。このようにして3’末端のみを増幅したバーコード付きcDNAをNGSで配列決定することによりscRNA-seqが可能になる。なお、DropSeqはコストが低いが、細胞の取得率と検出感度が低い弱点がある。inDropはDropSeqより細胞取得率が高く、パラメータを調整することで低レベルで発現される遺伝子の検出にも有利である[25]。 DropSeqのセットアップはDolomite Bio ([3])、inDropは1 Cellbio社から販売されている[4]。しかし、特に重要なのは10xGenomics社が同様の原理を用いた「Chromium」と命名された機器と試薬のシステムを市販することで、多くの研究者に利用できることになったことである[5]。Svenssonらによる最近のデータベース[6], [7]では、scRNA-seqを用いた論文で用いられた方法について調査しているが、この数年、10xGenomics社Chromiumを用いた論文が飛躍的に増加し、scRNA-seqの方法として、最も一般的になりつつあることがわかる(現在、10XGenomics社とBioRad社の間で関連特許をめぐる係争がある。)。このシステムは市販であるので導入が容易であり、DropSeqやinDropに比べ、多くの転写産物の高感度検出が可能であるが、ランニングコストは高い[26]

scRNA-seqの実際

ここでは主流になっている10xGenomics社のChromiumを用いた方法とSMART-seqなどを用いた方法に共通する方法の実際について俯瞰する。scRNA-seqの利用には、4つのステップがある(図2)[27]。1)個体や組織を採集し、そこから細胞あるいは細胞核を個別にすること。2)ChromiumやSMART-seq2などによる個々の細胞からのライブラリーの作製とNGSシーケンシング。3)前処理(preprocessing、得られた配列の整理)。4)ダウンストリーム解析(生物学的な情報を得る)。これらのうち、2)の段階については、上に記述したように市販の機器や試薬を利用する機会が多くなっているので、そのためのマニュアル等を参考にするのが現実的である。

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

血液細胞のように浮遊した細胞ではない場合、物理的あるいは酵素処理などによって、生組織から状態の良い細胞をdissociationする必要がある。神経系組織の酵素処理には、パパインを用いる方法が広く用いられている。ただ、しばしば問題となるのが、酵素処理のため短時間加温することで、発現量が変化する遺伝子が存在することである[28]。例えば、脳のミクログリアの解析には、低温下で組織をホモゲナイズするなどの工夫が必要であった[29]。また、このような現象を抑制するために、酵素処理時に転写阻害剤であるアクチノマイシンで処理したり[30]、ヒマラヤ氷河から得られた細菌Bacillus licheniformisから得られた低温プロテアーゼを用いる方法も報告されている[31]。また、細胞解離後に、メタノールで固定しscRNA-seqに使用することも可能である[32]。 なお、ヒト組織などから生きた細胞を得ることは困難なことが多い。この場合、凍結した組織から、各細胞由来の核を調製し、核内のmRNAを分析する方法 snRNA-seq (single-nucleus RNA-seq)があり、細胞質を持つ生細胞を利用した場合より感度は劣るが、細胞の同定などの目的には十分使用に耐える[33]<<[34]<[35][36][37][38][8]。snRNA-seqでは、組織をそのまま凍結することから開始するので、上述したscRNA-seqの問題である酵素処理や加温などを避けることができる。こうしたプロトコールの一部は、protocols.ioのHuman Cell Atlasのグループ[9]で公開されている。

単離した細胞は、そのまま10xGenomicsのChromiumのプラットフォームに導入することができるが、抗体などを用いたFACS、パニング、MACS(磁気ビーズカラム)などによる特定のマーカーを細胞表面などに発現する細胞の単離を行う場合もある。

更に、RNAを分析するscRNA-seqではないが、single cell genome-seqの変法として、シングルセルの遺伝子発現を推定する方法として、トランスポゾンを用いることでゲノムのオープンクロマチン領域を選択的に検出し、ライブラリーを作製しシーケンスするsingle cell ATAC-seq (Assay for Transposase-Accessible Chromatin)[39], [40][41], single cell THS-seq (transposome hypersensitive-site) [42]がある。またDNAメチル化をシングルセルレベルで観察する方法も報告されている[43]

scRNA-seqデータの前処理

10x Genomics社のChromium、IlluminaのNGSを利用した場合、Cell Ranger(Linux上で作動)を用いて、各生物種ごとのレファレンス配列リスト(https://www.ncbi.nlm.nih.gov/grc)やEggNOG ([10])などを利用し、細胞とトランスクリプトーム(各遺伝子の発現)の対応マトリックスを作製する。その後のデータの処理についても、10xGenomics社がソフトウェアLoupeを提供している。しかしながら、その後のダウンストリーム解析を考慮して、R言語, Python, MATLABなどのデータ解析のための汎用プログラミング言語やコードで扱えるオブジェクトに変換するのが通常である。ここでは、scRNA-seqデータ解析のために最もよく利用されているR言語を用いたパッケージ「Seurat」[44] [45]を中心に紹介したい。なお、Pythonを利用したものでは、ドイツ・ミュンヘンInstitute of Computational Biologyの Fabian J. Theisらが開発しているScanpyが有名である[46]

New York UniversityのRahul Satija研究室が開発しているSeuratは、scRNA-seq解析のために広く利用されているR言語のパッケージであり、2019年8月現在、その最新バージョンはSeurat 3.1である。論文の正式発表前から、サポート情報提供やコード修正なども頻繁に行っており、Satija研究室のWeb site(satijalab.org/Seurat, [11])、Github([12])、更にTwitterアカウント(@satijalab)などで最新情報を得ることできる。

最初に行うのは、scRNA-seqデータの品質管理である。ここでは、質の低い細胞のデータ(転写産物の種類が少ない、ミトコンドリア由来の転写産物が多い)を取り除く。また、複数の試料を組み合わせる場合には、バッチごとの違いについて検討する。特にDropletを使用するscRNA-seqの多くのケースで問題になるのが、Dropletに2つ以上の細胞が封じ込められ、それらが同じCell barcodeを持つために生じるアーティファクトである。通常Doubletsと呼ばれるこの問題はダウンストリーム解析を混乱させるので、細胞単離の段階から注意する必要があるが、明確なマーカー遺伝子が知られていればscRNA-seqデータ取得後にある程度のデータ処理で検討することは可能である。

このようなノーマライゼーションの過程を経て、scRNA-seqのデータ解析において、最初に行うのが、次元圧縮 (dimensionality reduction)である[47] 。PCA (Principal component analysis, 主成分分析)、UMAP(Uniform Manifold Approximation and Projection, 均一マニフォールド近似と投影)、Diffusion maps, t-SNE(t-distributed Stochastic Neighbor Embedding , t分布型確率的近傍埋込み)などの手法が用いられる。 特に、t-SNE[13](ティースニーと読むのが通常)は、高次元データを低次元の点の集合として可視化することで、それぞれの細胞の持つトランスクリプトームの類似度についての直観的な表示が可能でありしばしば用いられる(図3)。次に、Louvainアルゴリズムなどでクラスタリング(コミュニティ分割)を行い、tSNEグラフ上に表示できる。こうして、違ったタイプの細胞の集合が別のクラスターとして表示される。

ダウンストリーム解析

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

scRNA-seqデータから得られる生物学的知見には、内在的に存在する細胞の種類、外部刺激や環境で変化した細胞の状態、そして種類や変化により特徴的に発現するマーカー遺伝子候補の発見がある。クラスタリングにより、異なった細胞集団の存在が認識されると、それぞれのクラスターに特徴的に発現している遺伝子を具体的に探索し、細胞集団の持つバイオマーカーによって、そのクラスターの同定が可能になる。例えば、既にニューロンとグリア細胞に特異的に発現する典型的マーカーはよく知られており、それぞれのクラスターの識別は容易である。更に、ニューロンのタイプごとに区別されるマーカーや神経活動により変化したニューロンの状態は、In situ hybridizationや免疫組織化学などにより確認できる。このようなクラスターごとに発現が異なる遺伝子(差次的発現遺伝子)を見つけるためには(Differential expression analysis, DE analysis)、SeuratのFindMarkersコマンドでも利用可能である目的別の解析のための専用コード(MAST [48]、DESeq2 [49] など)を用いることができる。scRNA-seqの解析に必要なコードは、scRNA-tools [14], Awesome single cell [15], Bioconductor[16]で紹介されており、ほとんどがダウンロード可能である。また、最新の情報については、bioRxivなどのプレプリントサーバで公開されていることが多く、scRNA-seqのデータ(下記参考)とともに、オープンサイエンス実践の好例となっている。細胞ごとの差次的発現遺伝子の可視化には、ドットプロットやヴァイオリンプロットなどが頻繁に用いられる(図4)。

偽時系列解析、制御ネットワーク、パスウェイ解析

実験的なノイズとは別に生物学的に意味のある遺伝子発現の変動には、位置情報、細胞周期、概日リズム、破裂型プロモーターの作動などの理由で変動が見られるものもある[50][51]。特に、刺激・薬剤処理やさまざまな病態の進行や治療に伴う細胞の変化、発生途上の細胞系譜や細胞分化といった細胞の遷移状態の解析(偽時系列解析Pseudo-time analysis )には、scRNA-seqデータを用いることが極めて効果的である[52][53][54] 。これらの分析のためには軌道推定(Trajectory inference)の解析手法が用いられる。しばしば用いられるMonocle3 [55]など、多くのコードを収集しているGithubのサイトがある [17][18]。RNA velocityといった、転写産物のスプライシングの状態から細胞の分化状態を推定する方法もある[56]。しかし、これらの方法は、あくまで発生途上の細胞系譜や細胞分化の推定に過ぎない。細胞系譜を更に確実に観察しつつ、scRNA-seqを行うことで、細胞タイプの系統関係を調べる方法として、CRISPR-Cas9を用いたゲノム編集による記録法を導入したscGESTALT[57]、ScarTrace[58] 、LINNAEUS[59]がある。

また変化している遺伝子発現を解釈するために、制御ネットワーク(例、SCENIC[60], [19])やパスウェイ解析(例、Metascape[61], [20])といったシステム生物学で用いられてきた手法も適用できる。

神経科学への応用

ニューロンのクラスとタイプ

様々な神経・精神疾患について理解しその診断や治療に役立てるためには、ニューロン、グリア細胞を中心にした神経系にある細胞の「タイプ」を識別し、それぞれの細胞における分子的な変化を観察することが重要である [62][63]。近年、中枢神経系のグリア細胞にも、多様なアストロサイト、オリゴデンドロサイト、ミクログリアの存在が報告されてきている。一方で、ニューロンは著しく多様であり、このニューロンの多様性こそが、神経系を特徴づけており、その多彩で複雑な機能の発現に必須であることは疑う余地がない。 解剖学的な視点から言えば、すべてのニューロンの存在する位置は異なるので、すべてのニューロンは異なるという見方もできる。しかし、これは極論であり、従来の神経科学では、ニューロンの多様性は、それぞれのニューロンの解剖学的な位置、発現している分子、電気生理学、結合性、形態、神経伝達物質、神経伝達物質受容体とシグナル伝達によって識別されてきた。こうしたニューロンの多様性を便宜的に記述するのに、タイプ(type)、クラス(class)、サブクラス(subclass)、サブタイプ(subtype) というような用語が用いられてきた。しかし、本稿では混乱を防ぐため、Masland(2004)[64]が提唱し、広く受けいれられている「クラス」と「タイプ」という単語を用いることとする。タイプは、これ以上分類することができないとされる階層である。例えば、大脳皮質の錐体細胞、網膜神経節細胞といった大雑把な識別は「クラス」と呼ぶ。大脳皮質の錐体細胞というクラスは、層や領野によって「タイプ」が異なるし、網膜神経節細胞には視覚情報によって応答が異なる「タイプ」が存在する。この分類は、免疫組織化学、形態、電気生理学などの技術により識別可能である暫定的なものに過ぎない。本稿で解説するscRNA-seqの技術は、その網羅性からそれぞれのニューロンについてこれまでにないビッグデータを提供することで、このニューロンのタイプの理解に確実な根拠を与えつつある。

大脳

その他のCNS

疾患

アルツハイマー、Autism

網膜

展望

進化と他の動物種

データベース

Human Cell Atlas Human Brain Transcriptome project Single cell portal Allen Brain Atlas 統合 LIGER, MetaNeighbor

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

マルチモーダルなオミクス

関連項目

参考文献 

  1. Wang, Z., Gerstein, M., & Snyder, M. (2009).
    RNA-Seq: a revolutionary tool for transcriptomics. Nature reviews. Genetics, 10(1), 57-63. [PubMed:19015660] [PMC] [WorldCat] [DOI]
  2. Stark, R., Grzelak, M., & Hadfield, J. (2019).
    RNA sequencing: the teenage years. Nature reviews. Genetics, 20(11), 631-656. [PubMed:31341269] [WorldCat] [DOI]
  3. Picot, J., Guerin, C.L., Le Van Kim, C., & Boulanger, C.M. (2012).
    Flow cytometry: retrospective, fundamentals and recent instrumentation. Cytotechnology, 64(2), 109-30. [PubMed:22271369] [PMC] [WorldCat] [DOI]
  4. Spitzer, M.H., & Nolan, G.P. (2016).
    Mass Cytometry: Single Cells, Many Features. Cell, 165(4), 780-91. [PubMed:27153492] [PMC] [WorldCat] [DOI]
  5. Livesey, F.J. (2003).
    Strategies for microarray analysis of limiting amounts of RNA. Briefings in functional genomics & proteomics, 2(1), 31-6. [PubMed:15239941] [WorldCat] [DOI]
  6. Buck, L., & Axel, R. (1991).
    A novel multigene family may encode odorant receptors: a molecular basis for odor recognition. Cell, 65(1), 175-87. [PubMed:1840504] [WorldCat] [DOI]
  7. Dulac, C., & Axel, R. (1995).
    A novel family of genes encoding putative pheromone receptors in mammals. Cell, 83(2), 195-206. [PubMed:7585937] [WorldCat] [DOI]
  8. Tanabe, Y., William, C., & Jessell, T.M. (1998).
    Specification of motor neuron identity by the MNR2 homeodomain protein. Cell, 95(1), 67-80. [PubMed:9778248] [WorldCat] [DOI]
  9. Yamagata, M., Weiner, J.A., & Sanes, J.R. (2002).
    Sidekicks: synaptic adhesion molecules that promote lamina-specific connectivity in the retina. Cell, 110(5), 649-60. [PubMed:12230981] [WorldCat] [DOI]
  10. Kurimoto, K., Yabuta, Y., Ohinata, Y., Ono, Y., Uno, K.D., Yamada, R.G., ..., & Saitou, M. (2006).
    An improved single-cell cDNA amplification method for efficient high-density oligonucleotide microarray analysis. Nucleic acids research, 34(5), e42. [PubMed:16547197] [PMC] [WorldCat] [DOI]
  11. Tang, F., Barbacioru, C., Wang, Y., Nordman, E., Lee, C., Xu, N., ..., & Surani, M.A. (2009).
    mRNA-Seq whole-transcriptome analysis of a single cell. Nature methods, 6(5), 377-82. [PubMed:19349980] [WorldCat] [DOI]
  12. Ramsköld, D., Luo, S., Wang, Y.C., Li, R., Deng, Q., Faridani, O.R., ..., & Sandberg, R. (2012).
    Full-length mRNA-Seq from single-cell levels of RNA and individual circulating tumor cells. Nature biotechnology, 30(8), 777-82. [PubMed:22820318] [PMC] [WorldCat] [DOI]
  13. Picelli, S., Björklund, Å.K., Faridani, O.R., Sagasser, S., Winberg, G., & Sandberg, R. (2013).
    Smart-seq2 for sensitive full-length transcriptome profiling in single cells. Nature methods, 10(11), 1096-8. [PubMed:24056875] [WorldCat] [DOI]
  14. Picelli, S., Faridani, O.R., Björklund, A.K., Winberg, G., Sagasser, S., & Sandberg, R. (2014).
    Full-length RNA-seq from single cells using Smart-seq2. Nature protocols, 9(1), 171-81. [PubMed:24385147] [WorldCat] [DOI]
  15. Jaitin, D.A., Kenigsberg, E., Keren-Shaul, H., Elefant, N., Paul, F., Zaretsky, I., ..., & Amit, I. (2014).
    Massively parallel single-cell RNA-seq for marker-free decomposition of tissues into cell types. Science (New York, N.Y.), 343(6172), 776-9. [PubMed:24531970] [PMC] [WorldCat] [DOI]
  16. Islam, S., Kjällquist, U., Moliner, A., Zajac, P., Fan, J.B., Lönnerberg, P., & Linnarsson, S. (2011).
    Characterization of the single-cell transcriptional landscape by highly multiplex RNA-seq. Genome research, 21(7), 1160-7. [PubMed:21543516] [PMC] [WorldCat] [DOI]
  17. Islam, S., Zeisel, A., Joost, S., La Manno, G., Zajac, P., Kasper, M., ..., & Linnarsson, S. (2014).
    Quantitative single-cell RNA-seq with unique molecular identifiers. Nature methods, 11(2), 163-6. [PubMed:24363023] [WorldCat] [DOI]
  18. Hashimshony, T., Wagner, F., Sher, N., & Yanai, I. (2012).
    CEL-Seq: single-cell RNA-Seq by multiplexed linear amplification. Cell reports, 2(3), 666-73. [PubMed:22939981] [WorldCat] [DOI]
  19. Hashimshony, T., Senderovich, N., Avital, G., Klochendler, A., de Leeuw, Y., Anavy, L., ..., & Yanai, I. (2016).
    CEL-Seq2: sensitive highly-multiplexed single-cell RNA-Seq. Genome biology, 17, 77. [PubMed:27121950] [PMC] [WorldCat] [DOI]
  20. Gierahn, T.M., Wadsworth, M.H., Hughes, T.K., Bryson, B.D., Butler, A., Satija, R., ..., & Shalek, A.K. (2017).
    Seq-Well: portable, low-cost RNA sequencing of single cells at high throughput. Nature methods, 14(4), 395-398. [PubMed:28192419] [PMC] [WorldCat] [DOI]
  21. Cao, J., Packer, J.S., Ramani, V., Cusanovich, D.A., Huynh, C., Daza, R., ..., & Shendure, J. (2017).
    Comprehensive single-cell transcriptional profiling of a multicellular organism. Science (New York, N.Y.), 357(6352), 661-667. [PubMed:28818938] [PMC] [WorldCat] [DOI]
  22. Rosenberg, A.B., Roco, C.M., Muscat, R.A., Kuchina, A., Sample, P., Yao, Z., ..., & Seelig, G. (2018).
    Single-cell profiling of the developing mouse brain and spinal cord with split-pool barcoding. Science (New York, N.Y.), 360(6385), 176-182. [PubMed:29545511] [WorldCat] [DOI]
  23. Klein, A.M., Mazutis, L., Akartuna, I., Tallapragada, N., Veres, A., Li, V., ..., & Kirschner, M.W. (2015).
    Droplet barcoding for single-cell transcriptomics applied to embryonic stem cells. Cell, 161(5), 1187-1201. [PubMed:26000487] [PMC] [WorldCat] [DOI]
  24. Macosko, E.Z., Basu, A., Satija, R., Nemesh, J., Shekhar, K., Goldman, M., ..., & McCarroll, S.A. (2015).
    Highly Parallel Genome-wide Expression Profiling of Individual Cells Using Nanoliter Droplets. Cell, 161(5), 1202-1214. [PubMed:26000488] [PMC] [WorldCat] [DOI]
  25. Zhang, X., Li, T., Liu, F., Chen, Y., Yao, J., Li, Z., ..., & Wang, J. (2019).
    Comparative Analysis of Droplet-Based Ultra-High-Throughput Single-Cell RNA-Seq Systems. Molecular cell, 73(1), 130-142.e5. [PubMed:30472192] [WorldCat] [DOI]
  26. Zhang, X., Li, T., Liu, F., Chen, Y., Yao, J., Li, Z., ..., & Wang, J. (2019).
    Comparative Analysis of Droplet-Based Ultra-High-Throughput Single-Cell RNA-Seq Systems. Molecular cell, 73(1), 130-142.e5. [PubMed:30472192] [WorldCat] [DOI]
  27. Hwang, B., Lee, J.H., & Bang, D. (2018).
    Single-cell RNA sequencing technologies and bioinformatics pipelines. Experimental & molecular medicine, 50(8), 96. [PubMed:30089861] [PMC] [WorldCat] [DOI]
  28. Lacar, B., Linker, S.B., Jaeger, B.N., Krishnaswami, S.R., Barron, J.J., Kelder, M.J.E., ..., & Gage, F.H. (2016).
    Nuclear RNA-seq of single neurons reveals molecular signatures of activation. Nature communications, 7, 11022. [PubMed:27090946] [PMC] [WorldCat] [DOI]
  29. Hammond, T.R., Dufort, C., Dissing-Olesen, L., Giera, S., Young, A., Wysoker, A., ..., & Stevens, B. (2019).
    Single-Cell RNA Sequencing of Microglia throughout the Mouse Lifespan and in the Injured Brain Reveals Complex Cell-State Changes. Immunity, 50(1), 253-271.e6. [PubMed:30471926] [PMC] [WorldCat] [DOI]
  30. Wu, Y.E., Pan, L., Zuo, Y., Li, X., & Hong, W. (2017).
    Detecting Activated Cell Populations Using Single-Cell RNA-Seq. Neuron, 96(2), 313-329.e6. [PubMed:29024657] [WorldCat] [DOI]
  31. Adam, M., Potter, A.S., & Potter, S.S. (2017).
    Psychrophilic proteases dramatically reduce single-cell RNA-seq artifacts: a molecular atlas of kidney development. Development (Cambridge, England), 144(19), 3625-3632. [PubMed:28851704] [PMC] [WorldCat] [DOI]
  32. Alles, J., Karaiskos, N., Praktiknjo, S.D., Grosswendt, S., Wahle, P., Ruffault, P.L., ..., & Rajewsky, N. (2017).
    Cell fixation and preservation for droplet-based single-cell transcriptomics. BMC biology, 15(1), 44. [PubMed:28526029] [PMC] [WorldCat] [DOI]
  33. Habib, N., Li, Y., Heidenreich, M., Swiech, L., Avraham-Davidi, I., Trombetta, J.J., ..., & Regev, A. (2016).
    Div-Seq: Single-nucleus RNA-Seq reveals dynamics of rare adult newborn neurons. Science (New York, N.Y.), 353(6302), 925-8. [PubMed:27471252] [PMC] [WorldCat] [DOI]
  34. Habib, N., Avraham-Davidi, I., Basu, A., Burks, T., Shekhar, K., Hofree, M., ..., & Regev, A. (2017).
    Massively parallel single-nucleus RNA-seq with DroNc-seq. Nature methods, 14(10), 955-958. [PubMed:28846088] [PMC] [WorldCat] [DOI]
  35. Hu, P., Fabyanic, E., Kwon, D.Y., Tang, S., Zhou, Z., & Wu, H. (2017).
    Dissecting Cell-Type Composition and Activity-Dependent Transcriptional State in Mammalian Brains by Massively Parallel Single-Nucleus RNA-Seq. Molecular cell, 68(5), 1006-1015.e7. [PubMed:29220646] [PMC] [WorldCat] [DOI]
  36. Habib, N., Avraham-Davidi, I., Basu, A., Burks, T., Shekhar, K., Hofree, M., ..., & Regev, A. (2017).
    Massively parallel single-nucleus RNA-seq with DroNc-seq. Nature methods, 14(10), 955-958. [PubMed:28846088] [PMC] [WorldCat] [DOI]
  37. Bakken, T.E., Hodge, R.D., Miller, J.A., Yao, Z., Nguyen, T.N., Aevermann, B., ..., & Tasic, B. (2018).
    Single-nucleus and single-cell transcriptomes compared in matched cortical cell types. PloS one, 13(12), e0209648. [PubMed:30586455] [PMC] [WorldCat] [DOI]
  38. Lake, B.B., Codeluppi, S., Yung, Y.C., Gao, D., Chun, J., Kharchenko, P.V., ..., & Zhang, K. (2017).
    A comparative strategy for single-nucleus and single-cell transcriptomes confirms accuracy in predicted cell-type expression from nuclear RNA. Scientific reports, 7(1), 6031. [PubMed:28729663] [PMC] [WorldCat] [DOI]
  39. Buenrostro, J.D., Wu, B., Litzenburger, U.M., Ruff, D., Gonzales, M.L., Snyder, M.P., ..., & Greenleaf, W.J. (2015).
    Single-cell chromatin accessibility reveals principles of regulatory variation. Nature, 523(7561), 486-90. [PubMed:26083756] [PMC] [WorldCat] [DOI]
  40. Preissl, S., Fang, R., Huang, H., Zhao, Y., Raviram, R., Gorkin, D.U., ..., & Ren, B. (2018).
    Single-nucleus analysis of accessible chromatin in developing mouse forebrain reveals cell-type-specific transcriptional regulation. Nature neuroscience, 21(3), 432-439. [PubMed:29434377] [PMC] [WorldCat] [DOI]
  41. Cusanovich, D.A., Daza, R., Adey, A., Pliner, H.A., Christiansen, L., Gunderson, K.L., ..., & Shendure, J. (2015).
    Multiplex single cell profiling of chromatin accessibility by combinatorial cellular indexing. Science (New York, N.Y.), 348(6237), 910-4. [PubMed:25953818] [PMC] [WorldCat] [DOI]
  42. Lake, B.B., Chen, S., Sos, B.C., Fan, J., Kaeser, G.E., Yung, Y.C., ..., & Zhang, K. (2018).
    Integrative single-cell analysis of transcriptional and epigenetic states in the human adult brain. Nature biotechnology, 36(1), 70-80. [PubMed:29227469] [WorldCat] [DOI]
  43. Luo, C., Keown, C.L., Kurihara, L., Zhou, J., He, Y., Li, J., ..., & Ecker, J.R. (2017).
    Single-cell methylomes identify neuronal subtypes and regulatory elements in mammalian cortex. Science (New York, N.Y.), 357(6351), 600-604. [PubMed:28798132] [PMC] [WorldCat] [DOI]
  44. Butler, A., Hoffman, P., Smibert, P., Papalexi, E., & Satija, R. (2018).
    Integrating single-cell transcriptomic data across different conditions, technologies, and species. Nature biotechnology, 36(5), 411-420. [PubMed:29608179] [PMC] [WorldCat] [DOI]
  45. Stuart, T., Butler, A., Hoffman, P., Hafemeister, C., Papalexi, E., Mauck, W.M., ..., & Satija, R. (2019).
    Comprehensive Integration of Single-Cell Data. Cell, 177(7), 1888-1902.e21. [PubMed:31178118] [PMC] [WorldCat] [DOI]
  46. Wolf, F.A., Angerer, P., & Theis, F.J. (2018).
    SCANPY: large-scale single-cell gene expression data analysis. Genome biology, 19(1), 15. [PubMed:29409532] [PMC] [WorldCat] [DOI]
  47. Kiselev, V.Y., Andrews, T.S., & Hemberg, M. (2019).
    Challenges in unsupervised clustering of single-cell RNA-seq data. Nature reviews. Genetics, 20(5), 273-282. [PubMed:30617341] [WorldCat] [DOI]
  48. Finak, G., McDavid, A., Yajima, M., Deng, J., Gersuk, V., Shalek, A.K., ..., & Gottardo, R. (2015).
    MAST: a flexible statistical framework for assessing transcriptional changes and characterizing heterogeneity in single-cell RNA sequencing data. Genome biology, 16, 278. [PubMed:26653891] [PMC] [WorldCat] [DOI]
  49. Love, M.I., Huber, W., & Anders, S. (2014).
    Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome biology, 15(12), 550. [PubMed:25516281] [PMC] [WorldCat] [DOI]
  50. Luecken, M.D., & Theis, F.J. (2019).
    Current best practices in single-cell RNA-seq analysis: a tutorial. Molecular systems biology, 15(6), e8746. [PubMed:31217225] [PMC] [WorldCat] [DOI]
  51. Kolodziejczyk, A.A., Kim, J.K., Svensson, V., Marioni, J.C., & Teichmann, S.A. (2015).
    The technology and biology of single-cell RNA sequencing. Molecular cell, 58(4), 610-20. [PubMed:26000846] [WorldCat] [DOI]
  52. Cembrowski, M.S., & Menon, V. (2018).
    Continuous Variation within Cell Types of the Nervous System. Trends in neurosciences, 41(6), 337-348. [PubMed:29576429] [WorldCat] [DOI]
  53. Marioni, J.C., & Arendt, D. (2017).
    How Single-Cell Genomics Is Changing Evolutionary and Developmental Biology. Annual review of cell and developmental biology, 33, 537-553. [PubMed:28813177] [WorldCat] [DOI]
  54. Telley, L., & Jabaudon, D. (2018).
    A mixed model of neuronal diversity. Nature, 555(7697), 452-454. [PubMed:29565398] [WorldCat] [DOI]
  55. Cao, J., Spielmann, M., Qiu, X., Huang, X., Ibrahim, D.M., Hill, A.J., ..., & Shendure, J. (2019).
    The single-cell transcriptional landscape of mammalian organogenesis. Nature, 566(7745), 496-502. [PubMed:30787437] [PMC] [WorldCat] [DOI]
  56. La Manno, G., Soldatov, R., Zeisel, A., Braun, E., Hochgerner, H., Petukhov, V., ..., & Kharchenko, P.V. (2018).
    RNA velocity of single cells. Nature, 560(7719), 494-498. [PubMed:30089906] [PMC] [WorldCat] [DOI]
  57. Raj, B., Wagner, D.E., McKenna, A., Pandey, S., Klein, A.M., Shendure, J., ..., & Schier, A.F. (2018).
    Simultaneous single-cell profiling of lineages and cell types in the vertebrate brain. Nature biotechnology, 36(5), 442-450. [PubMed:29608178] [PMC] [WorldCat] [DOI]
  58. Alemany, A., Florescu, M., Baron, C.S., Peterson-Maduro, J., & van Oudenaarden, A. (2018).
    Whole-organism clone tracing using single-cell sequencing. Nature, 556(7699), 108-112. [PubMed:29590089] [WorldCat] [DOI]
  59. Spanjaard, B., Hu, B., Mitic, N., Olivares-Chauvet, P., Janjuha, S., Ninov, N., & Junker, J.P. (2018).
    Simultaneous lineage tracing and cell-type identification using CRISPR-Cas9-induced genetic scars. Nature biotechnology, 36(5), 469-473. [PubMed:29644996] [PMC] [WorldCat] [DOI]
  60. Aibar, S., González-Blas, C.B., Moerman, T., Huynh-Thu, V.A., Imrichova, H., Hulselmans, G., ..., & Aerts, S. (2017).
    SCENIC: single-cell regulatory network inference and clustering. Nature methods, 14(11), 1083-1086. [PubMed:28991892] [PMC] [WorldCat] [DOI]
  61. Zhou, Y., Zhou, B., Pache, L., Chang, M., Khodabakhshi, A.H., Tanaseichuk, O., ..., & Chanda, S.K. (2019).
    Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nature communications, 10(1), 1523. [PubMed:30944313] [PMC] [WorldCat] [DOI]
  62. Zeng, H., & Sanes, J.R. (2017).
    Neuronal cell-type classification: challenges, opportunities and the path forward. Nature reviews. Neuroscience, 18(9), 530-546. [PubMed:28775344] [WorldCat] [DOI]
  63. Tasic, B. (2018).
    Single cell transcriptomics in neuroscience: cell classification and beyond. Current opinion in neurobiology, 50, 242-249. [PubMed:29738987] [WorldCat] [DOI]
  64. Masland, R.H. (2004).
    Neuronal cell types. Current biology : CB, 14(13), R497-500. [PubMed:15242626] [WorldCat] [DOI]