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


21行目: 21行目:


一方で多くの種類の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>、Seq-Well <ref><pubmed> 28192419</pubmed></ref>
それ以来、完全長cDNAを増幅したり、細胞ごとに異なる分子識別子(unique molecular identifiers: UMI)を持つcDNAを増幅させるscRNA-seqが考案され始め、2013年には、このような1細胞のシーケンシング技術が、Nature Methods誌のMethod of the Year に選ばれた<ref> https://www.nature.com/collections/mysbdwgfll<ref>。たとえば、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合成のためのプロトコルである(既に、SMART-seq3という改良プロトコールもあるhttps://doi.org/10.1101/817924)。また、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>、Microwell-seq<ref><pubmed>29474909</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を作製する方法も報告されている。
などが報告されてきた。最近になって、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を作製する方法も報告されている。
 
これらの方法のうち、SMART-seq、その改良法であるSMART-seq2は、微小ピペットによるマニュアル捕獲、セルソーター、レーザー捕獲法などを用いる多穴プレート法、更に半導体集積回路製作技術で作った流体集積回路を利用するFluidigm C1の装置[https://jp.fluidigm.com]と組み合わせることで利用される機会が多い<ref><pubmed>30405621</pubmed></ref>。このSMART-seq2プロトコールの特徴は、mRNAの全領域を読むことで、全長トランスクリプトームを得ることができることであり、mRNAのスプライシングバリアントなどのアイソフォーム、アリルごとの発現情報が得られるSNPs、変異の検出にも利用できる。また、それぞれ細胞ごとの反応を独立した場所で行うため、別の細胞の反応と混じる可能性がない。これらの点が、次に説明するDropletを使用して3’末端のみを標的にしたscRNA-seqに比べた場合の長所であるが、その高コスト(1細胞あたり数十ドル)と処理可能な細胞数の少なさが短所である。


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


===Droplet使用の3’エンドリード法===
===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>。
しかしながら、もっとも重要な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細胞あたり約5セント)。つまり、細胞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>。
http://www.youtube.com/watch?v=fHq9ewdYEWM
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>。
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などを用いた方法に共通する方法の実際について俯瞰する。scRNA-seqの利用には、4つのステップがある(図2)<ref><pubmed>30089861</pubmed></ref>。1)個体や組織を採集し、そこから細胞あるいは細胞核を個別にすること。2)ChromiumやSMART-seq2などによる個々の細胞からのライブラリーの作製とNGSシーケンシング。3)前処理(preprocessing、得られた配列の整理)。4)ダウンストリーム解析(生物学的な情報を得る)。これらのうち、2)の段階については、上に記述したように市販の機器や試薬を利用する機会が多くなっているので、そのためのマニュアル等を参考にするのが現実的である。
ここでは主流になっている10xGenomics社のChromiumを用いた方法とSMART-seq2などを用いた方法に共通する方法の実際について俯瞰する。scRNA-seqの利用には、4つのステップがある(図2)<ref><pubmed>30089861</pubmed></ref>。1)個体や組織を採集し、そこから細胞あるいは細胞核を個別にすること。2)ChromiumやSMART-seq2などによる個々の細胞からのライブラリーの作製とNGSシーケンシング。3)前処理(preprocessing、得られた配列の整理)。4)ダウンストリーム解析(生物学的な情報を得る)。これらのうち、2)の段階については、上に記述したように市販の機器や試薬を利用する機会が多くなっているので、そのためのマニュアル等を参考にするのが現実的である。


===組織からの細胞、細胞核の分離===
===組織からの細胞、細胞核の分離===
血液細胞のように浮遊した細胞ではない場合、物理的あるいは酵素処理などによって、生組織から状態の良い細胞を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>。
血液細胞のように浮遊した細胞ではない場合、物理的あるいは酵素処理などによって解離することで、生組織から状態の良いバラバラになった細胞を調製する必要がある。神経系組織の酵素処理には、パパインを用いる方法が広く用いられている。ここで、しばしば問題となるのが、酵素処理のため短時間加温することで、発現量が変化する遺伝子が存在することである<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>。
なお、ヒト組織などから生きた細胞を得ることは困難なことが多い。この場合、凍結した組織から、各細胞由来の核を調製し、核内の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]で公開されている。
 
単離した細胞は、そのまま10xGenomicsのChromiumのプラットフォームに導入することができるが、抗体や蛍光タンパク質レポーターなどを用いたFACS、パニング、MACS(磁気ビーズカラム)などによる特定のマーカーを細胞表面などに発現する細胞の選択的濃縮や除去を行う場合もある。更に、抗体にUMIをカップリングさせるCITE-seqについては、下記のマルチモーダルなオミクスの項目で述べる。
 
なお、ヒト組織や希少生物などから生細胞を得ることは困難なことが多い。この場合、scRNA-seqの変法として、凍結した組織から、各細胞由来の核を調製し、核内のmRNAを分析するsnRNA-seq (single-nucleus RNA-seq)が利用されている。ただ、この場合、FACSなどによる特定細胞集団の同定が難しく、細胞質を持つ生細胞を利用した場合と同等な結果が必ずしも得られない<ref><pubmed>24248345</pubmed></ref><ref><pubmed>26890679</pubmed></ref>  <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><ref><pubmed>31728515</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]で公開されている。


単離した細胞は、そのまま10xGenomicsのChromiumのプラットフォームに導入することができるが、抗体や蛍光タンパク質レポータ=などを用いたFACS、パニング、MACS(磁気ビーズカラム)などによる特定のマーカーを細胞表面などに発現する細胞の単離を行う場合もある。更に、抗体にUMIをカップリングさせることで、細胞を単離せずにscRNA-seqを行うCITE-seqなどの方法については、下記のマルチモーダルなオミクスの項で述べる。


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


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データ取得後にある程度のデータ処理で検討することは可能である。また、この問題を解決する新たなアプローチも試みられている<ref><pubmed>30954476</pubmed></ref> <ref><pubmed>31836005</pubmed></ref><ref><pubmed>31856883</pubmed></ref> [https://doi.org/10.1101/2019.12.17.879304][https://doi.org/10.1101/699637]


最初に行うのは、scRNA-seqデータの品質管理である。ここでは、質の低い細胞のデータ(転写産物の種類が少ない、ミトコンドリア由来の転写産物が多い)を取り除く。また、複数の試料を組み合わせる場合には、バッチごとの違いについて検討する。特にDropletを使用するscRNA-seqの多くのケースで問題になるのが、Dropletに2つ以上の細胞が封じ込められ、それらが同じCell barcodeを持つために生じるアーティファクトである。通常Doubletsと呼ばれるこの問題はダウンストリーム解析を混乱させるので、細胞単離の段階から注意する必要があるが、明確なマーカー遺伝子が知られていればscRNA-seqデータ取得後にある程度のデータ処理で検討することは可能である。
このようなノーマライゼーションの過程を経て、scRNA-seqのデータ解析において、最初に行うのが、次元圧縮 (dimensionality reduction)である<ref><pubmed>30617341</pubmed></ref><ref><pubmed>31780648</pubmed></ref>
。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グラフ上に表示できる。こうして、違ったタイプの細胞の集合が別のクラスターとして表示される。しかし、データによっては、tSNEだけでなく、他の方法でパラメータの調節を行うことで違ったデータ解釈(別の細胞クラスターの同定)ができるケースもあり、いくつかの方法を試してみることが推奨される。


このようなノーマライゼーションの過程を経て、scRNA-seqのデータ解析において、最初に行うのが、次元圧縮 (dimensionality reduction)である<ref><pubmed> 30617341</pubmed></ref>
。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データから得られる生物学的知見には、内在的に存在する細胞の種類、外部刺激や環境で変化した細胞の状態、そして種類や変化により特徴的に発現するマーカー遺伝子候補の発見がある<ref><pubmed>27824854</pubmed></ref><ref><pubmed>32033589</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)。
。クラスタリングにより、異なった細胞集団の存在が認識されると、それぞれのクラスターに特徴的に発現している遺伝子を具体的に探索し、細胞集団の持つバイオマーカーによって、そのクラスターの同定が可能になる。例えば、既にニューロンとグリア細胞に特異的に発現する典型的マーカーはよく知られており、それぞれのクラスターの識別は容易である。更に、ニューロンのタイプごとに区別されるマーカーや神経活動により変化したニューロンの状態は、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>
実験的なノイズとは別に生物学的に意味のある遺伝子発現の変動には、位置情報、細胞周期、概日リズム、破裂型プロモーターの作動などの理由で変動が見られるものもある<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>がある。
。これらの分析のためには軌道推定(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>がある。
64行目: 72行目:
また変化している遺伝子発現を解釈するために、制御ネットワーク(例、SCENIC<ref><pubmed>28991892</pubmed></ref>, [https://github.com/aertslab/SCENIC])やパスウェイ解析(例、Metascape<ref><pubmed>30944313</pubmed></ref>, [http://metascape.org])といったシステム生物学で用いられてきた手法も適用できる。
また変化している遺伝子発現を解釈するために、制御ネットワーク(例、SCENIC<ref><pubmed>28991892</pubmed></ref>, [https://github.com/aertslab/SCENIC])やパスウェイ解析(例、Metascape<ref><pubmed>30944313</pubmed></ref>, [http://metascape.org])といったシステム生物学で用いられてきた手法も適用できる。


==神経科学への応用==
==scRNA-seqの神経科学研究への適用==
===ニューロンのタイプ===
===神経細胞多様性の新たな根拠===
様々な神経・精神疾患について理解しその診断や治療に役立てるためには、ニューロン、グリア細胞を中心にした神経系にある細胞の「タイプ」を識別し、それぞれの細胞における分子的な変化を観察することが重要である <ref><pubmed>28775344</pubmed></ref><ref><pubmed>29738987</pubmed></ref>。近年、中枢神経系のグリア細胞にも、多様なアストロサイト、オリゴデンドロサイト、ミクログリアの存在が報告されてきている。一方で、ニューロンは著しく多様であり、このニューロンの多様性こそが、神経系を特徴づけており、その多彩で複雑な機能の発現に必須であることは疑う余地がない。
様々な神経・精神疾患について理解しその診断や治療に役立てるためには、神経細胞、グリア細胞を中心にした神経系にある細胞の「タイプ」を識別し、それぞれの細胞における分子的な変化を観察することが重要である <ref><pubmed>28775344</pubmed></ref><ref><pubmed>29738987</pubmed></ref>。近年、中枢神経系のグリア細胞にも、多様なアストロサイト、オリゴデンドロサイト、ミクログリアの存在が報告されてきている。一方で、神経細胞は著しく多様であり、この多様性が神経系を特徴づけており、その多彩で複雑な機能の発現に必須であることは疑う余地がない。解剖学的な視点から言えば、すべての神経細胞の存在する位置は異なるので、すべての神経細胞は異なるという見方もできる。しかし、従来の神経科学では、神経細胞の多様性は、それぞれの神経細胞の解剖学的な位置、発現している分子、電気生理学、結合性、形態、神経伝達物質、神経伝達物質受容体とシグナル伝達によって識別されてきている。こうした神経細胞の多様性を便宜的に記述するのに、タイプ(type)、クラス(class)、サブクラス(subclass)、サブタイプ(subtype) というような用語が用いられてきた。しかし、本稿では混乱を防ぐため、Masland(2004)<ref><pubmed>15242626</pubmed></ref>が提唱し、広く受けいれられている「クラス」と「タイプ」という単語を用いることとする。タイプは、これ以上分類することができないとされる階層である。例えば、大脳皮質の錐体細胞、網膜神経節細胞といった大雑把な識別は「クラス」と呼ぶ。大脳皮質の錐体細胞というクラスは、層や領野によって「タイプ」が異なるし、網膜神経節細胞には視覚情報によって応答が異なる「タイプ」が存在する。この分類は、免疫組織化学、形態、電気生理学などの技術により識別可能である暫定的なものに過ぎない。本項目で解説してきたscRNA-seqの技術は、その網羅性からそれぞれの神経細胞についてこれまでにないビッグデータを提供することで、この神経細胞のタイプの理解に新たな根拠を与えつつある。
解剖学的な視点から言えば、すべてのニューロンの存在する位置は異なるので、すべてのニューロンは異なるという見方もできる。しかし、これは極論であり、従来の神経科学では、ニューロンの多様性は、それぞれのニューロンの解剖学的な位置、発現している分子、電気生理学、結合性、形態、神経伝達物質、神経伝達物質受容体とシグナル伝達によって識別されてきた。こうしたニューロンの多様性を便宜的に記述するのに、タイプ(type)、クラス(class)、サブクラス(subclass)、サブタイプ(subtype) というような用語が用いられてきた。しかし、本稿では混乱を防ぐため、Masland(2004)<ref><pubmed>15242626</pubmed></ref>が提唱し、広く受けいれられている「クラス」と「タイプ」という単語を用いることとする。タイプは、これ以上分類することができないとされる階層である。例えば、大脳皮質の錐体細胞、網膜神経節細胞といった大雑把な識別は「クラス」と呼ぶ。大脳皮質の錐体細胞というクラスは、層や領野によって「タイプ」が異なるし、網膜神経節細胞には視覚情報によって応答が異なる「タイプ」が存在する。この分類は、免疫組織化学、形態、電気生理学などの技術により識別可能である暫定的なものに過ぎない。本稿で解説するscRNA-seqの技術は、その網羅性からそれぞれのニューロンについてこれまでにないビッグデータを提供することで、このニューロンのタイプの理解に確実な根拠を与えつつある。
 
===神経細胞へのscRNA-seqの適用===
大脳皮質には、錐体細胞や非錐体細胞などの神経細胞、更に種々のグリア細胞などが見られる。マウスを用いた初期のscRNA-seqを使った小規模な細胞数を分類した研究で、既に、これまで知られていた主要な細胞タイプを超えるタイプが見つかった<ref><pubmed>25700174</pubmed></ref>。その後のDroplet使用の3’エンドリード法を利用した多数の細胞数の解析で、更に多数の神経細胞のタイプが見つかっている<ref><pubmed>30096299</pubmed></ref><ref><pubmed>30096314</pubmed></ref><ref><pubmed>30382198</pubmed></ref><ref><pubmed>29320739</pubmed></ref><ref><pubmed>28846088</pubmed></ref>。特に、GABA作動性介在神経細胞タイプの多様性<ref><pubmed>28942923</pubmed></ref><ref><pubmed>29472441</pubmed></ref><ref><pubmed>28134272</pubmed></ref><ref><pubmed>29513653</pubmed></ref>  や発生過程<ref><pubmed>31073041</pubmed></ref><ref><pubmed>230635555</pubmed></ref><ref><pubmed>29539641</pubmed></ref><ref><pubmed>30485812</pubmed></ref><ref><pubmed>26940868</pubmed></ref>の理解が、scRNA-Seq技術を利用することで進んでいる。更に、神経活動によって変化するトランスクリプトームの変化も細胞ごとに調査され興味深い<ref><pubmed>29230054</pubmed></ref> 。
ヒトを含めた霊長類の大脳についても発達段階を含めてscRNA-seqが適用されてきている<ref><pubmed>28846088</pubmed></ref><ref><pubmed>26060301</pubmed></ref>  <ref><pubmed>27339989</pubmed></ref><ref><pubmed>29539641</pubmed></ref>     <ref><pubmed>31303374</pubmed></ref>  <ref><pubmed>29227469</pubmed></ref> <ref><pubmed>29217575</pubmed></ref><ref><pubmed>31435019</pubmed></ref> <ref><pubmed>29867213</pubmed></ref>  。https://doi.org/10.1101/709501。
島にあり、ヒトや霊長類に特徴的とされるvon Economoニューロン(紡錘細胞)のような希少な細胞のscRNA-seqにも成功している<ref><pubmed>32127543</pubmed></ref>。
 
海馬<ref><pubmed>29912866</pubmed></ref><ref><pubmed>29335606</pubmed></ref><ref><pubmed>31942070</pubmed></ref><ref><pubmed>29241552</pubmed></ref>では、これまでの研究で記載されてきた神経細胞のタイプが確認され、更に新しいタイプが見つかった。中枢神経系では、外側膝状体<ref><pubmed>29343640</pubmed></ref>、大脳基底核(足底核)<ref><pubmed>28384468</pubmed></ref> 、視床下部<ref><pubmed>28355573</pubmed></ref>  <ref><pubmed>27991900</pubmed></ref>  <ref><pubmed>31249056</pubmed></ref>  <ref><pubmed>30858605</pubmed></ref><ref><pubmed>28166221</pubmed></ref>  <ref><pubmed>30385464</pubmed></ref>  、線条体<ref><pubmed>31875543</pubmed></ref>  <ref><pubmed>27425622</pubmed></ref>  、
中脳<ref><pubmed>27716510</pubmed></ref>  <ref><pubmed>30718509</pubmed></ref>  <ref><pubmed>29499164</pubmed></ref>  、発生中の間脳<ref><pubmed>30872278</pubmed></ref> 、さらに小脳<ref><pubmed>30735127</pubmed></ref><ref><pubmed>30220501</pubmed></ref><ref><pubmed>30690467</pubmed></ref><ref><pubmed>29576475</pubmed></ref>などの結果が得られている。マウスの小脳においては、分子層にこれまでの星状細胞、バスケット細胞というカテゴリーとは違った2種類の神経細胞があることが示唆されている[https://doi.org/10.1101/2020.03.04.976407]。脳の外部では、感覚神経<ref><pubmed>25420068</pubmed></ref>  <ref><pubmed>26691752</pubmed></ref>  、らせん神経節<ref><pubmed>3007810</pubmed></ref> 、臭覚神経<ref><pubmed>26541607</pubmed></ref>  、腸神経系 <ref><pubmed>29483303</pubmed></ref>[https://doi.org/10.1101/2020.03.02.955757] 、網膜<ref><pubmed>27565351</pubmed></ref><ref><pubmed>29909983</pubmed></ref> <ref><pubmed>31260032</pubmed></ref><ref><pubmed>31128945</pubmed></ref><ref><pubmed>30018341</pubmed></ref><ref><pubmed>30712875</pubmed></ref><ref><pubmed>30548510</pubmed></ref><ref><pubmed>31075224</pubmed></ref><ref><pubmed>31399471</pubmed></ref><ref><pubmed>31848347</pubmed></ref><ref><pubmed>31673015</pubmed></ref><ref><pubmed>31653841</pubmed></ref><ref><pubmed>31784286</pubmed></ref>[https://doi.org/10.1101/2020.02.26.966093][https://www.biorxiv.org/content/10.1101/779694][https://www.biorxiv.org/content/10.1101/617555]でのscRNA-seqデータがある。またiPS細胞やES細胞由来のオルガノイドに含まれる神経細胞タイプを知る上でも利用されている<ref><pubmed>31168097</pubmed></ref><ref><pubmed>28094016</pubmed></ref><ref><pubmed>28279351</pubmed></ref><ref><pubmed>31996853</pubmed></ref><ref><pubmed>31968264</pubmed></ref>。
 
 
===神経細胞以外の細胞===
上衣細胞<ref><pubmed>29727663</pubmed></ref>は、神経幹細胞としての役割が示唆されてきたが、scRNA-seqの結果ではその可能性が支持されていない。グリア細胞では、ラジアルグリア<ref><pubmed>29217575</pubmed></ref><ref><pubmed>26406371</pubmed></ref><ref><pubmed>29539641</pubmed></ref><ref><pubmed>25734491</pubmed></ref><ref><pubmed>29281841</pubmed></ref>  、アストロサイト<ref><pubmed>32139688</pubmed></ref><ref><pubmed>32203496</pubmed></ref>に多様性があることが示唆されてきている。また、オリゴデンドロサイト<ref><pubmed>27284195</pubmed></ref>については、これまで細胞生物学的に研究されてきた分化の過程がscRNA-seqにより検出されている。


===大脳===
ミクログリアは、その活性化に伴うトランスクリプトームの変化の様子がscRNA-seqにより詳細に明らかになった
<ref><pubmed>29539641</pubmed></ref><ref><pubmed>30206190</pubmed></ref>
<ref><pubmed>27338705</pubmed></ref>
<ref><pubmed>30471926</pubmed></ref>
<ref><pubmed>31209379</pubmed></ref>
<ref><pubmed>29020624</pubmed></ref>
<ref><pubmed>31835035</pubmed></ref>。また、CNS境界関連マクロファージ(BAM) <ref><pubmed>31061494</pubmed></ref>のscRNA-seqも実施されている。


===その他の神経系===


===疾患===
===疾患===
アルツハイマー、Autism
scRNA-seqは、疾患の理解にも有用である。筋萎縮性側索硬化症ALS <ref><pubmed>30948552</pubmed></ref>、アルツハイマー病やそのモデル動物<ref><pubmed>31042697</pubmed></ref><ref><pubmed>31399126</pubmed></ref>
===網膜===
[https://doi.org/10.1101/628347]<ref><pubmed>28602351</pubmed></ref>、自閉症やレット症候<ref><pubmed>31097668</pubmed></ref><ref><pubmed>30455458</pubmed></ref>
、シャルコー・マリー・トゥース病<ref><pubmed>29888333</pubmed></ref>、ダウン症[https://www.biorxiv.org/content/10.1101/2020.01.01.892398v1]、パーキンソン病<ref><pubmed>30503143</pubmed></ref>、がん<ref><pubmed>31327527</pubmed></ref><ref><pubmed>28360267</pubmed></ref>などに適用されている。
 
==scRNA-seqの展望==
===神経系の多様性と進化===
scRNA-seqは、既に多様な生物の神経系の細胞の理解、更には種間の相同性や差異の研究に利用されており、神経系の進化を細胞レベルで考察するのに有用であろう(例、カタユウレイボヤCiona intestinalis <ref><pubmed>30069052</pubmed></ref><ref><pubmed>30228204</pubmed></ref>、ショウジョウバエ  <ref><pubmed>29909982</pubmed></ref><ref><pubmed>29149607</pubmed></ref><ref><pubmed>30703584</pubmed></ref><ref><pubmed>29909983</pubmed></ref>、ゼブラフィッシュ<ref><pubmed>31018142</pubmed></ref><ref><pubmed>30929901</pubmed></ref>、アカミミガメTrachemys scripta elegans、トカゲPogona vitticeps, pv<ref><pubmed>29724907</pubmed></ref>、ニワトリ、霊長類<ref><pubmed>30730291</pubmed></ref><ref><pubmed>31619793</pubmed></ref>[https://doi.org/10.1101/2020.03.31.016972
])。ただ、遺伝子やトランスクリプトームの研究が進んでいる生物種では比較的容易であるが、アノテーションが十分でない生物種を用いる場合、scRNA-seqのデータ解析は困難である。この場合は、NCBIのTaxonomy[https://www.ncbi.nlm.nih.gov/taxonomy]やEggNOG [http://eggnogdb.embl.de] <ref><pubmed>30418610</pubmed></ref>を利用する。


==展望==
===他の動物種と進化===
===データベースと統合===
===データベースと統合===
Human Cell Atlas
scRNA-seqのデータは様々な目的で利用できるため、データベース化し利用できるようにする必要がある。オープンサイエンスの典型として、神経系のトランスクリプトーム一般のデータベースが多数公開されており<ref><pubmed>29437890</pubmed></ref>、基本的にNCBIのGene Expression Omnibus[https://www.ncbi.nlm.nih.gov/geo/]に登録されている。また、common coordinate framework (CCF) やcentral annotation platform (CAP)という概念のもと、scRNA-seqを意識したものとして、米国のBRAIN Initiative Cell Census Consortium<ref><pubmed>29096072</pubmed></ref> 、Human Cell Atlas ProjectのHuman Cell Atlas Data Portal [https://data.humancellatlas.org]、アレン脳研究所のAllen Brain Atlas [https://portal.brain-map.org]、ブロード研究所のSingle Cell Portal, 
Human Brain Transcriptome project
[https://singlecell.broadinstitute.org/]が稼働している。
Single cell portal
 
Allen Brain Atlas
また異なった方法や実験で得られたscRNA-seqのデータを比較することも重要である。最近、LIGER<ref><pubmed>31178122</pubmed></ref> 、Harmony<ref><pubmed>31740819</pubmed></ref>  、MetaNeighber<ref><pubmed>29491377</pubmed></ref>のようなアルゴリズムが開発されている。
統合 LIGER, MetaNeighbor
 


===空間トランスクリプトミクス===
===空間トランスクリプトミクス===
scRNA-seqの弱点は、組織からそれぞれの細胞を解離する必要があるので、その細胞が存在していた空間的な位置の情報を失ってしまうということである。組織切片におけるタンパク質などの分布は免疫組織化学、mRNAの分布はin situ hybridizationで検出することができるが、数多くのmRNAの分布を一気に同定する方法がscRNA-seqと同様に開発されてきている。Slide-seq<ref><pubmed>30923225</pubmed></ref>、osmFISH<ref><pubmed>30377364</pubmed></ref>
、STARmap (spatially-resolved transcript amplicon readout mapping), <ref><pubmed>29930089</pubmed></ref>、seqFISH  <ref><pubmed>27764670</pubmed></ref>、pciSeq(probabilistic cell typing by in situ sequencing)[ https://doi.org/10.1101/431957
]、更に10xGenomics社のVisiumなどがある。現状では、組織を用いる空間トランスクリプトミクス(Spatial transcriptomics)は、空間解像度が1つずつの細胞のレベルにいたっておらず、課題も多い。しかし、そのデータを解析するためのアルゴリズム<ref><pubmed>29553578</pubmed></ref><ref><pubmed>29553579</pubmed></ref>[https://doi.org/10.1101/757096][ https://doi.org/10.1101/701680]
やMerFish  <ref><pubmed>25858977</pubmed></ref>、corrFISH  <ref><pubmed>27271198</pubmed></ref>のように、一つの細胞内でのそれぞれのmRNAを検出できる方法が開発されてきており、scRNA-seqの弱点を補う空間トランスクリプトームは今後の発展が期待される。


===マルチモーダルなオミクス===
===マルチモーダルなシングルセルオミクス===
同じ細胞からscRNA-seqの情報だけでなく、ゲノム配列、ATAC-seqなどによるエピゲノム解析、タンパク質、プロテオームなどを、同時に記録するマルチモーダルなオミクスomicsが注目されている<ref><pubmed>31907462</pubmed></ref><ref><pubmed>30696980</pubmed></ref>。Nature Methodsの2019年「Methods of the Year」に選ばれており、その特集号などを参考にされたい。神経科学分野で注目されるのは、パッチクランプの生理学的情報とscRNA-seqを組み合わせたPatch-seq <ref><pubmed>26689544</pubmed></ref> <ref><pubmed>26689543</pubmed></ref>であろう。また、細胞表面分子に対する抗体にDNAを付加することで、マーカーを発現する細胞のトランスクリプトームを観察するCITE-seq<ref><pubmed>28759029</pubmed></ref>、 REAP-seq<ref><pubmed>28854175</pubmed></ref>は細胞マーカー発現と遺伝子発現が同時に定量できるマルチモーダルなオミクスである。 また、BARseq (barcoded anatomy resolved by sequencing) <ref><pubmed>31626774</pubmed></ref>のような方法は、コネクトームと遺伝子発現を記録できるomicsとして興味深い。


== 関連項目 ==
== 関連項目 ==
*[[]]
 
*[[大脳皮質]]
*[[樹状突起]]
*[[軸索]]
*[[大脳皮質の局所神経回路]]


== 参考文献  ==
== 参考文献  ==
<references />
<references />

2020年4月3日 (金) 06:15時点における版

山形方人
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 に選ばれた引用エラー: <ref> タグに対応する </ref> タグが不足していますおよびその改良されたプロトコルであるSMART-seq2 [12] [13]は、完全長cDNA合成のためのプロトコルである(既に、SMART-seq3という改良プロトコールもあるhttps://doi.org/10.1101/817924)。また、MARS-seq(Massively parallel single-cell RNA-seq)[14]、STRT(single-cell tagged reverse transcription)[15] [16]、CEL-seq(Cell Expression by Linear amplification and Sequencing)[17]、CEL-seq2[18]、Seq-Well [19]、Microwell-seq[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の装置[1]と組み合わせることで利用される機会が多い[23]。このSMART-seq2プロトコールの特徴は、mRNAの全領域を読むことで、全長トランスクリプトームを得ることができることであり、mRNAのスプライシングバリアントなどのアイソフォーム、アリルごとの発現情報が得られるSNPs、変異の検出にも利用できる。また、それぞれ細胞ごとの反応を独立した場所で行うため、別の細胞の反応と混じる可能性がない。これらの点が、次に説明するDropletを使用して3’末端のみを標的にしたscRNA-seqに比べた場合の長所であるが、その高コスト(1細胞あたり数十ドル)と処理可能な細胞数の少なさが短所である。


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

しかしながら、もっとも重要なscRNA-seqの方法論についての進歩は、2015年、Harvard Medical Schoolの独立した2つのグループが、inDropそしてDrop-seqという類似した2つの高スループットな方法を開発したことであろう[24] [25]。これらの方法では、マイクロ流体力学 (Microfluidics) 、 Cell BarcodeとUMIとしてDNAバーコーディング (DNA barcoding) 、そしてNGSを利用することで、自動化とサンプル調製の容易さから、1つの細胞あたりに要するコストを大幅に低下させることに成功した(Drop-seqは発表時で、1細胞あたり約5セント)。つまり、細胞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より細胞取得率が高く、パラメータを調整することで低レベルで発現される遺伝子の検出にも有利である[26]http://www.youtube.com/watch?v=fHq9ewdYEWM DropSeqのセットアップはDolomite Bio ([2])、inDropは1 Cellbio社から販売されている[3]。しかし、特に重要なのは10xGenomics社が同様の原理を用いた「Chromium」と命名された機器と試薬のシステムを市販することで、多くの研究者に利用できることになったことである[4]。Svenssonらによる最近のデータベース[5], [6]では、scRNA-seqを用いた論文で用いられた方法について調査しているが、この数年、10xGenomics社Chromiumを用いた論文が飛躍的に増加し、scRNA-seqの方法として、最も一般的になりつつあることがわかる(現在、10XGenomics社とBioRad社の間で関連特許をめぐる係争がある。)。このシステムは市販であるので導入が容易であり、DropSeqやinDropに比べ、多くの転写産物の高感度検出が可能であるが、ランニングコストは高い[27]


scRNA-seqの実際

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

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

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

単離した細胞は、そのまま10xGenomicsのChromiumのプラットフォームに導入することができるが、抗体や蛍光タンパク質レポーターなどを用いたFACS、パニング、MACS(磁気ビーズカラム)などによる特定のマーカーを細胞表面などに発現する細胞の選択的濃縮や除去を行う場合もある。更に、抗体にUMIをカップリングさせるCITE-seqについては、下記のマルチモーダルなオミクスの項目で述べる。

なお、ヒト組織や希少生物などから生細胞を得ることは困難なことが多い。この場合、scRNA-seqの変法として、凍結した組織から、各細胞由来の核を調製し、核内のmRNAを分析するsnRNA-seq (single-nucleus RNA-seq)が利用されている。ただ、この場合、FACSなどによる特定細胞集団の同定が難しく、細胞質を持つ生細胞を利用した場合と同等な結果が必ずしも得られない[34][35] [36][37]<[38][39][40][41][42] [7] 。snRNA-seqでは、組織をそのまま凍結することから開始するので、上述したscRNA-seqの問題である酵素処理や加温などを避けることができる。こうしたプロトコールの一部は、protocols.ioのHuman Cell Atlasのグループ[8]で公開されている。


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


scRNA-seqデータの前処理

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

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

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

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


ダウンストリーム解析

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

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

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

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

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

scRNA-seqの神経科学研究への適用

神経細胞多様性の新たな根拠

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

神経細胞へのscRNA-seqの適用

大脳皮質には、錐体細胞や非錐体細胞などの神経細胞、更に種々のグリア細胞などが見られる。マウスを用いた初期のscRNA-seqを使った小規模な細胞数を分類した研究で、既に、これまで知られていた主要な細胞タイプを超えるタイプが見つかった[75]。その後のDroplet使用の3’エンドリード法を利用した多数の細胞数の解析で、更に多数の神経細胞のタイプが見つかっている[76][77][78][79][80]。特に、GABA作動性介在神経細胞タイプの多様性[81][82][83][84] や発生過程[85][86][87][88][89]の理解が、scRNA-Seq技術を利用することで進んでいる。更に、神経活動によって変化するトランスクリプトームの変化も細胞ごとに調査され興味深い[90] 。 ヒトを含めた霊長類の大脳についても発達段階を含めてscRNA-seqが適用されてきている[91][92] [93][94]  [95] [96] [97][98] [99]https://doi.org/10.1101/709501。 島にあり、ヒトや霊長類に特徴的とされるvon Economoニューロン(紡錘細胞)のような希少な細胞のscRNA-seqにも成功している[100]

海馬[101][102][103][104]では、これまでの研究で記載されてきた神経細胞のタイプが確認され、更に新しいタイプが見つかった。中枢神経系では、外側膝状体[105]、大脳基底核(足底核)[106] 、視床下部[107] [108] [109] [110][111] [112] 、線条体[113] [114] 、 中脳[115] [116] [117] 、発生中の間脳[118] 、さらに小脳[119][120][121][122]などの結果が得られている。マウスの小脳においては、分子層にこれまでの星状細胞、バスケット細胞というカテゴリーとは違った2種類の神経細胞があることが示唆されている[23]。脳の外部では、感覚神経[123] [124] 、らせん神経節[125] 、臭覚神経[126] 、腸神経系 [127][24] 、網膜[128][129] [130][131][132][133][134][135][136][137][138][139][140][25][26][27]でのscRNA-seqデータがある。またiPS細胞やES細胞由来のオルガノイドに含まれる神経細胞タイプを知る上でも利用されている[141][142][143][144][145]


神経細胞以外の細胞

上衣細胞[146]は、神経幹細胞としての役割が示唆されてきたが、scRNA-seqの結果ではその可能性が支持されていない。グリア細胞では、ラジアルグリア[147][148][149][150][151] 、アストロサイト[152][153]に多様性があることが示唆されてきている。また、オリゴデンドロサイト[154]については、これまで細胞生物学的に研究されてきた分化の過程がscRNA-seqにより検出されている。

ミクログリアは、その活性化に伴うトランスクリプトームの変化の様子がscRNA-seqにより詳細に明らかになった [155][156] [157] [158] [159] [160] [161]。また、CNS境界関連マクロファージ(BAM) [162]のscRNA-seqも実施されている。


疾患

scRNA-seqは、疾患の理解にも有用である。筋萎縮性側索硬化症ALS [163]、アルツハイマー病やそのモデル動物[164][165] [28][166]、自閉症やレット症候[167][168] 、シャルコー・マリー・トゥース病[169]、ダウン症[29]、パーキンソン病[170]、がん[171][172]などに適用されている。

scRNA-seqの展望

神経系の多様性と進化

scRNA-seqは、既に多様な生物の神経系の細胞の理解、更には種間の相同性や差異の研究に利用されており、神経系の進化を細胞レベルで考察するのに有用であろう(例、カタユウレイボヤCiona intestinalis [173][174]、ショウジョウバエ [175][176][177][178]、ゼブラフィッシュ[179][180]、アカミミガメTrachemys scripta elegans、トカゲPogona vitticeps, pv[181]、ニワトリ、霊長類[182][183][https://doi.org/10.1101/2020.03.31.016972 ])。ただ、遺伝子やトランスクリプトームの研究が進んでいる生物種では比較的容易であるが、アノテーションが十分でない生物種を用いる場合、scRNA-seqのデータ解析は困難である。この場合は、NCBIのTaxonomy[30]やEggNOG [31] [184]を利用する。

データベースと統合

scRNA-seqのデータは様々な目的で利用できるため、データベース化し利用できるようにする必要がある。オープンサイエンスの典型として、神経系のトランスクリプトーム一般のデータベースが多数公開されており[185]、基本的にNCBIのGene Expression Omnibus[32]に登録されている。また、common coordinate framework (CCF) やcentral annotation platform (CAP)という概念のもと、scRNA-seqを意識したものとして、米国のBRAIN Initiative Cell Census Consortium[186] 、Human Cell Atlas ProjectのHuman Cell Atlas Data Portal [33]、アレン脳研究所のAllen Brain Atlas [34]、ブロード研究所のSingle Cell Portal, [35]が稼働している。

また異なった方法や実験で得られたscRNA-seqのデータを比較することも重要である。最近、LIGER[187] 、Harmony[188] 、MetaNeighber[189]のようなアルゴリズムが開発されている。


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

scRNA-seqの弱点は、組織からそれぞれの細胞を解離する必要があるので、その細胞が存在していた空間的な位置の情報を失ってしまうということである。組織切片におけるタンパク質などの分布は免疫組織化学、mRNAの分布はin situ hybridizationで検出することができるが、数多くのmRNAの分布を一気に同定する方法がscRNA-seqと同様に開発されてきている。Slide-seq[190]、osmFISH[191] 、STARmap (spatially-resolved transcript amplicon readout mapping), [192]、seqFISH  [193]、pciSeq(probabilistic cell typing by in situ sequencing)[ https://doi.org/10.1101/431957 ]、更に10xGenomics社のVisiumなどがある。現状では、組織を用いる空間トランスクリプトミクス(Spatial transcriptomics)は、空間解像度が1つずつの細胞のレベルにいたっておらず、課題も多い。しかし、そのデータを解析するためのアルゴリズム[194][195][36][ https://doi.org/10.1101/701680] やMerFish [196]、corrFISH  [197]のように、一つの細胞内でのそれぞれのmRNAを検出できる方法が開発されてきており、scRNA-seqの弱点を補う空間トランスクリプトームは今後の発展が期待される。


マルチモーダルなシングルセルオミクス

同じ細胞からscRNA-seqの情報だけでなく、ゲノム配列、ATAC-seqなどによるエピゲノム解析、タンパク質、プロテオームなどを、同時に記録するマルチモーダルなオミクスomicsが注目されている[198][199]。Nature Methodsの2019年「Methods of the Year」に選ばれており、その特集号などを参考にされたい。神経科学分野で注目されるのは、パッチクランプの生理学的情報とscRNA-seqを組み合わせたPatch-seq [200] [201]であろう。また、細胞表面分子に対する抗体にDNAを付加することで、マーカーを発現する細胞のトランスクリプトームを観察するCITE-seq[202]、 REAP-seq[203]は細胞マーカー発現と遺伝子発現が同時に定量できるマルチモーダルなオミクスである。 また、BARseq (barcoded anatomy resolved by sequencing) [204]のような方法は、コネクトームと遺伝子発現を記録できるomicsとして興味深い。

関連項目

参考文献 

  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. 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]
  13. 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]
  14. 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]
  15. 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]
  16. 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]
  17. 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]
  18. 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]
  19. 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]
  20. Han, X., Wang, R., Zhou, Y., Fei, L., Sun, H., Lai, S., ..., & Guo, G. (2018).
    Mapping the Mouse Cell Atlas by Microwell-Seq. Cell, 172(5), 1091-1107.e17. [PubMed:29474909] [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. See, P., Lum, J., Chen, J., & Ginhoux, F. (2018).
    A Single-Cell Sequencing Guide for Immunologists. Frontiers in immunology, 9, 2425. [PubMed:30405621] [PMC] [WorldCat] [DOI]
  24. 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]
  25. 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]
  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. 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]
  28. 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]
  29. 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]
  30. 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]
  31. 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]
  32. 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]
  33. 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]
  34. Grindberg, R.V., Yee-Greenbaum, J.L., McConnell, M.J., Novotny, M., O'Shaughnessy, A.L., Lambert, G.M., ..., & Lasken, R.S. (2013).
    RNA-sequencing from single nuclei. Proceedings of the National Academy of Sciences of the United States of America, 110(49), 19802-7. [PubMed:24248345] [PMC] [WorldCat] [DOI]
  35. Krishnaswami, S.R., Grindberg, R.V., Novotny, M., Venepally, P., Lacar, B., Bhutani, K., ..., & Lasken, R.S. (2016).
    Using single nuclei for RNA-seq to capture the transcriptome of postmortem neurons. Nature protocols, 11(3), 499-524. [PubMed:26890679] [PMC] [WorldCat] [DOI]
  36. 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]
  37. 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]
  38. 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]
  39. 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]
  40. 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]
  41. 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]
  42. Amamoto, R., Zuccaro, E., Curry, N.C., Khurana, S., Chen, H.H., Cepko, C.L., & Arlotta, P. (2020).
    FIN-Seq: transcriptional profiling of specific cell types from frozen archived tissue of the human central nervous system. Nucleic acids research, 48(1), e4. [PubMed:31728515] [WorldCat] [DOI]
  43. 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]
  44. 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]
  45. 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]
  46. 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]
  47. 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]
  48. 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]
  49. 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]
  50. 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]
  51. Wolock, S.L., Lopez, R., & Klein, A.M. (2019).
    Scrublet: Computational Identification of Cell Doublets in Single-Cell Transcriptomic Data. Cell systems, 8(4), 281-291.e9. [PubMed:30954476] [PMC] [WorldCat] [DOI]
  52. Huang, Y., McCarthy, D.J., & Stegle, O. (2019).
    Vireo: Bayesian demultiplexing of pooled single-cell RNA-seq data without genotype reference. Genome biology, 20(1), 273. [PubMed:31836005] [PMC] [WorldCat] [DOI]
  53. Xu, J., Falconer, C., Nguyen, Q., Crawford, J., McKinnon, B.D., Mortlock, S., ..., & Coin, L.J.M. (2019).
    Genotype-free demultiplexing of pooled single-cell RNA-seq. Genome biology, 20(1), 290. [PubMed:31856883] [PMC] [WorldCat] [DOI]
  54. 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]
  55. Kobak, D., & Berens, P. (2019).
    The art of using t-SNE for single-cell transcriptomics. Nature communications, 10(1), 5416. [PubMed:31780648] [PMC] [WorldCat] [DOI]
  56. Wagner, A., Regev, A., & Yosef, N. (2016).
    Revealing the vectors of cellular identity with single-cell genomics. Nature biotechnology, 34(11), 1145-1160. [PubMed:27824854] [PMC] [WorldCat] [DOI]
  57. Lähnemann, D., Köster, J., Szczurek, E., McCarthy, D.J., Hicks, S.C., Robinson, M.D., ..., & Schönhuth, A. (2020).
    Eleven grand challenges in single-cell data science. Genome biology, 21(1), 31. [PubMed:32033589] [PMC] [WorldCat] [DOI]
  58. 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]
  59. 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]
  60. 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]
  61. 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]
  62. 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]
  63. 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]
  64. Telley, L., & Jabaudon, D. (2018).
    A mixed model of neuronal diversity. Nature, 555(7697), 452-454. [PubMed:29565398] [WorldCat] [DOI]
  65. 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]
  66. 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]
  67. 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]
  68. 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]
  69. 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]
  70. 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]
  71. 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]
  72. 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]
  73. Tasic, B. (2018).
    Single cell transcriptomics in neuroscience: cell classification and beyond. Current opinion in neurobiology, 50, 242-249. [PubMed:29738987] [WorldCat] [DOI]
  74. Masland, R.H. (2004).
    Neuronal cell types. Current biology : CB, 14(13), R497-500. [PubMed:15242626] [WorldCat] [DOI]
  75. Zeisel, A., Muñoz-Manchado, A.B., Codeluppi, S., Lönnerberg, P., La Manno, G., Juréus, A., ..., & Linnarsson, S. (2015).
    Brain structure. Cell types in the mouse cortex and hippocampus revealed by single-cell RNA-seq. Science (New York, N.Y.), 347(6226), 1138-42. [PubMed:25700174] [WorldCat] [DOI]
  76. Saunders, A., Macosko, E.Z., Wysoker, A., Goldman, M., Krienen, F.M., de Rivera, H., ..., & McCarroll, S.A. (2018).
    Molecular Diversity and Specializations among the Cells of the Adult Mouse Brain. Cell, 174(4), 1015-1030.e16. [PubMed:30096299] [PMC] [WorldCat] [DOI]
  77. Zeisel, A., Hochgerner, H., Lönnerberg, P., Johnsson, A., Memic, F., van der Zwan, J., ..., & Linnarsson, S. (2018).
    Molecular Architecture of the Mouse Nervous System. Cell, 174(4), 999-1014.e22. [PubMed:30096314] [PMC] [WorldCat] [DOI]
  78. Tasic, B., Yao, Z., Graybuck, L.T., Smith, K.A., Nguyen, T.N., Bertagnolli, D., ..., & Zeng, H. (2018).
    Shared and distinct transcriptomic cell types across neocortical areas. Nature, 563(7729), 72-78. [PubMed:30382198] [PMC] [WorldCat] [DOI]
  79. Chevée, M., Robertson, J.J., Cannon, G.H., Brown, S.P., & Goff, L.A. (2018).
    Variation in Activity State, Axonal Projection, and Position Define the Transcriptional Identity of Individual Neocortical Projection Neurons. Cell reports, 22(2), 441-455. [PubMed:29320739] [PMC] [WorldCat] [DOI]
  80. 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]
  81. Paul, A., Crow, M., Raudales, R., He, M., Gillis, J., & Huang, Z.J. (2017).
    Transcriptional Architecture of Synaptic Communication Delineates GABAergic Neuron Identity. Cell, 171(3), 522-539.e20. [PubMed:28942923] [PMC] [WorldCat] [DOI]
  82. Mi, D., Li, Z., Lim, L., Li, M., Moissidis, M., Yang, Y., ..., & Marín, O. (2018).
    Early emergence of cortical interneuron diversity in the mouse embryo. Science (New York, N.Y.), 360(6384), 81-85. [PubMed:29472441] [PMC] [WorldCat] [DOI]
  83. Frazer, S., Prados, J., Niquille, M., Cadilhac, C., Markopoulos, F., Gomez, L., ..., & Dayer, A. (2017).
    Transcriptomic and anatomic parcellation of 5-HT3AR expressing cortical interneuron subtypes revealed by single-cell RNA sequencing. Nature communications, 8, 14219. [PubMed:28134272] [PMC] [WorldCat] [DOI]
  84. Mayer, C., Hafemeister, C., Bandler, R.C., Machold, R., Batista Brito, R., Jaglin, X., ..., & Satija, R. (2018).
    Developmental diversification of cortical inhibitory interneurons. Nature, 555(7697), 457-462. [PubMed:29513653] [PMC] [WorldCat] [DOI]
  85. Telley, L., Agirman, G., Prados, J., Amberg, N., Fièvre, S., Oberst, P., ..., & Jabaudon, D. (2019).
    Temporal patterning of apical progenitors and their daughter neurons in the developing neocortex. Science (New York, N.Y.), 364(6440). [PubMed:31073041] [WorldCat] [DOI]
  86. Resource not found in PubMed.
  87. Zhong, S., Zhang, S., Fan, X., Wu, Q., Yan, L., Dong, J., ..., & Wang, X. (2018).
    A single-cell RNA-seq survey of the developmental landscape of the human prefrontal cortex. Nature, 555(7697), 524-528. [PubMed:29539641] [WorldCat] [DOI]
  88. Zywitza, V., Misios, A., Bunatyan, L., Willnow, T.E., & Rajewsky, N. (2018).
    Single-Cell Transcriptomics Characterizes Cell Types in the Subventricular Zone and Uncovers Molecular Defects Impairing Adult Neurogenesis. Cell reports, 25(9), 2457-2469.e8. [PubMed:30485812] [WorldCat] [DOI]
  89. Telley, L., Govindan, S., Prados, J., Stevant, I., Nef, S., Dermitzakis, E., ..., & Jabaudon, D. (2016).
    Sequential transcriptional waves direct the differentiation of newborn neurons in the mouse neocortex. Science (New York, N.Y.), 351(6280), 1443-6. [PubMed:26940868] [WorldCat] [DOI]
  90. Hrvatin, S., Hochbaum, D.R., Nagy, M.A., Cicconet, M., Robertson, K., Cheadle, L., ..., & Greenberg, M.E. (2018).
    Single-cell analysis of experience-dependent transcriptomic states in the mouse visual cortex. Nature neuroscience, 21(1), 120-129. [PubMed:29230054] [PMC] [WorldCat] [DOI]
  91. 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]
  92. Darmanis, S., Sloan, S.A., Zhang, Y., Enge, M., Caneda, C., Shuer, L.M., ..., & Quake, S.R. (2015).
    A survey of human brain transcriptome diversity at the single cell level. Proceedings of the National Academy of Sciences of the United States of America, 112(23), 7285-90. [PubMed:26060301] [PMC] [WorldCat] [DOI]
  93. Lake, B.B., Ai, R., Kaeser, G.E., Salathia, N.S., Yung, Y.C., Liu, R., ..., & Zhang, K. (2016).
    Neuronal subtypes and diversity revealed by single-nucleus RNA sequencing of the human brain. Science (New York, N.Y.), 352(6293), 1586-90. [PubMed:27339989] [PMC] [WorldCat] [DOI]
  94. Zhong, S., Zhang, S., Fan, X., Wu, Q., Yan, L., Dong, J., ..., & Wang, X. (2018).
    A single-cell RNA-seq survey of the developmental landscape of the human prefrontal cortex. Nature, 555(7697), 524-528. [PubMed:29539641] [WorldCat] [DOI]
  95. Polioudakis, D., de la Torre-Ubieta, L., Langerman, J., Elkins, A.G., Shi, X., Stein, J.L., ..., & Geschwind, D.H. (2019).
    A Single-Cell Transcriptomic Atlas of Human Neocortical Development during Mid-gestation. Neuron, 103(5), 785-801.e8. [PubMed:31303374] [PMC] [WorldCat] [DOI]
  96. 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]
  97. Nowakowski, T.J., Bhaduri, A., Pollen, A.A., Alvarado, B., Mostajo-Radji, M.A., Di Lullo, E., ..., & Kriegstein, A.R. (2017).
    Spatiotemporal gene expression trajectories reveal developmental hierarchies of the human cortex. Science (New York, N.Y.), 358(6368), 1318-1323. [PubMed:29217575] [PMC] [WorldCat] [DOI]
  98. Hodge, R.D., Bakken, T.E., Miller, J.A., Smith, K.A., Barkan, E.R., Graybuck, L.T., ..., & Lein, E.S. (2019).
    Conserved cell types with divergent features in human versus mouse cortex. Nature, 573(7772), 61-68. [PubMed:31435019] [PMC] [WorldCat] [DOI]
  99. Fan, X., Dong, J., Zhong, S., Wei, Y., Wu, Q., Yan, L., ..., & Tang, F. (2018).
    Spatial transcriptomic survey of human embryonic cerebral cortex by single-cell RNA-seq analysis. Cell research, 28(7), 730-745. [PubMed:29867213] [PMC] [WorldCat] [DOI]
  100. Hodge, R.D., Miller, J.A., Novotny, M., Kalmbach, B.E., Ting, J.T., Bakken, T.E., ..., & Lein, E.S. (2020).
    Transcriptomic evidence that von Economo neurons are regionally specialized extratelencephalic-projecting excitatory neurons. Nature communications, 11(1), 1172. [PubMed:32127543] [PMC] [WorldCat] [DOI]
  101. Harris, K.D., Hochgerner, H., Skene, N.G., Magno, L., Katona, L., Bengtsson Gonzales, C., ..., & Hjerling-Leffler, J. (2018).
    Classes and continua of hippocampal CA1 inhibitory neurons revealed by single-cell transcriptomics. PLoS biology, 16(6), e2006387. [PubMed:29912866] [PMC] [WorldCat] [DOI]
  102. Hochgerner, H., Zeisel, A., Lönnerberg, P., & Linnarsson, S. (2018).
    Conserved properties of dentate gyrus neurogenesis across postnatal development revealed by single-cell RNA sequencing. Nature neuroscience, 21(2), 290-299. [PubMed:29335606] [WorldCat] [DOI]
  103. Zhong, S., Ding, W., Sun, L., Lu, Y., Dong, H., Fan, X., ..., & Wang, X. (2020).
    Decoding the development of the human hippocampus. Nature, 577(7791), 531-536. [PubMed:31942070] [WorldCat] [DOI]
  104. Artegiani, B., Lyubimova, A., Muraro, M., van Es, J.H., van Oudenaarden, A., & Clevers, H. (2017).
    A Single-Cell RNA Sequencing Study Reveals Cellular and Molecular Dynamics of the Hippocampal Neurogenic Niche. Cell reports, 21(11), 3271-3284. [PubMed:29241552] [WorldCat] [DOI]
  105. Kalish, B.T., Cheadle, L., Hrvatin, S., Nagy, M.A., Rivera, S., Crow, M., ..., & Greenberg, M.E. (2018).
    Single-cell transcriptomics of the developing lateral geniculate nucleus reveals insights into circuit assembly and refinement. Proceedings of the National Academy of Sciences of the United States of America, 115(5), E1051-E1060. [PubMed:29343640] [PMC] [WorldCat] [DOI]
  106. Wallace, M.L., Saunders, A., Huang, K.W., Philson, A.C., Goldman, M., Macosko, E.Z., ..., & Sabatini, B.L. (2017).
    Genetically Distinct Parallel Pathways in the Entopeduncular Nucleus for Limbic and Sensorimotor Output of the Basal Ganglia. Neuron, 94(1), 138-152.e5. [PubMed:28384468] [PMC] [WorldCat] [DOI]
  107. Chen, R., Wu, X., Jiang, L., & Zhang, Y. (2017).
    Single-Cell RNA-Seq Reveals Hypothalamic Cell Diversity. Cell reports, 18(13), 3227-3241. [PubMed:28355573] [PMC] [WorldCat] [DOI]
  108. Romanov, R.A., Zeisel, A., Bakker, J., Girach, F., Hellysaz, A., Tomer, R., ..., & Harkany, T. (2017).
    Molecular interrogation of hypothalamic organization reveals distinct dopamine neuronal subtypes. Nature neuroscience, 20(2), 176-188. [PubMed:27991900] [WorldCat] [DOI]
  109. Rossi, M.A., Basiri, M.L., McHenry, J.A., Kosyk, O., Otis, J.M., van den Munkhof, H.E., ..., & Stuber, G.D. (2019).
    Obesity remodels activity and transcriptional state of a lateral hypothalamic brake on feeding. Science (New York, N.Y.), 364(6447), 1271-1274. [PubMed:31249056] [WorldCat] [DOI]
  110. Mickelsen, L.E., Bolisetty, M., Chimileski, B.R., Fujita, A., Beltrami, E.J., Costanzo, J.T., ..., & Jackson, A.C. (2019).
    Single-cell transcriptomic analysis of the lateral hypothalamic area reveals molecularly distinct populations of inhibitory and excitatory neurons. Nature neuroscience, 22(4), 642-656. [PubMed:30858605] [PMC] [WorldCat] [DOI]
  111. Campbell, J.N., Macosko, E.Z., Fenselau, H., Pers, T.H., Lyubetskaya, A., Tenen, D., ..., & Tsai, L.T. (2017).
    A molecular census of arcuate hypothalamus and median eminence cell types. Nature neuroscience, 20(3), 484-496. [PubMed:28166221] [PMC] [WorldCat] [DOI]
  112. Moffitt, J.R., Bambah-Mukku, D., Eichhorn, S.W., Vaughn, E., Shekhar, K., Perez, J.D., ..., & Zhuang, X. (2018).
    Molecular, spatial, and functional single-cell profiling of the hypothalamic preoptic region. Science (New York, N.Y.), 362(6416). [PubMed:30385464] [PMC] [WorldCat] [DOI]
  113. Märtin, A., Calvigioni, D., Tzortzi, O., Fuzik, J., Wärnberg, E., & Meletis, K. (2019).
    A Spatiomolecular Map of the Striatum. Cell reports, 29(13), 4320-4333.e5. [PubMed:31875543] [WorldCat] [DOI]
  114. Gokce, O., Stanley, G.M., Treutlein, B., Neff, N.F., Camp, J.G., Malenka, R.C., ..., & Quake, S.R. (2016).
    Cellular Taxonomy of the Mouse Striatum as Revealed by Single-Cell RNA-Seq. Cell reports, 16(4), 1126-1137. [PubMed:27425622] [PMC] [WorldCat] [DOI]
  115. La Manno, G., Gyllborg, D., Codeluppi, S., Nishimura, K., Salto, C., Zeisel, A., ..., & Linnarsson, S. (2016).
    Molecular Diversity of Midbrain Development in Mouse, Human, and Stem Cells. Cell, 167(2), 566-580.e19. [PubMed:27716510] [PMC] [WorldCat] [DOI]
  116. Tiklová, K., Björklund, Å.K., Lahti, L., Fiorenzano, A., Nolbrant, S., Gillberg, L., ..., & Perlmann, T. (2019).
    Single-cell RNA sequencing reveals midbrain dopamine neuron diversity emerging during mouse brain development. Nature communications, 10(1), 581. [PubMed:30718509] [PMC] [WorldCat] [DOI]
  117. Hook, P.W., McClymont, S.A., Cannon, G.H., Law, W.D., Morton, A.J., Goff, L.A., & McCallion, A.S. (2018).
    Single-Cell RNA-Seq of Mouse Dopaminergic Neurons Informs Candidate Gene Selection for Sporadic Parkinson Disease. American journal of human genetics, 102(3), 427-446. [PubMed:29499164] [PMC] [WorldCat] [DOI]
  118. Guo, Q., & Li, J.Y.H. (2019).
    Defining developmental diversification of diencephalon neurons through single cell gene expression profiling. Development (Cambridge, England), 146(12). [PubMed:30872278] [PMC] [WorldCat] [DOI]
  119. Wizeman, J.W., Guo, Q., Wilion, E.M., & Li, J.Y. (2019).
    Specification of diverse cell types during early neurogenesis of the mouse cerebellum. eLife, 8. [PubMed:30735127] [PMC] [WorldCat] [DOI]
  120. Carter, R.A., Bihannic, L., Rosencrance, C., Hadley, J.L., Tong, Y., Phoenix, T.N., ..., & Gawad, C. (2018).
    A Single-Cell Transcriptional Atlas of the Developing Murine Cerebellum. Current biology : CB, 28(18), 2910-2920.e2. [PubMed:30220501] [WorldCat] [DOI]
  121. Peng, J., Sheng, A.L., Xiao, Q., Shen, L., Ju, X.C., Zhang, M., ..., & Luo, Z.G. (2019).
    Single-cell transcriptomes reveal molecular specializations of neuronal cell types in the developing cerebellum. Journal of molecular cell biology, 11(8), 636-648. [PubMed:30690467] [PMC] [WorldCat] [DOI]
  122. Pandey, S., Shekhar, K., Regev, A., & Schier, A.F. (2018).
    Comprehensive Identification and Spatial Mapping of Habenular Neuronal Types Using Single-Cell RNA-Seq. Current biology : CB, 28(7), 1052-1065.e7. [PubMed:29576475] [PMC] [WorldCat] [DOI]
  123. Usoskin, D., Furlan, A., Islam, S., Abdo, H., Lönnerberg, P., Lou, D., ..., & Ernfors, P. (2015).
    Unbiased classification of sensory neuron types by large-scale single-cell RNA sequencing. Nature neuroscience, 18(1), 145-53. [PubMed:25420068] [WorldCat] [DOI]
  124. Li, C.L., Li, K.C., Wu, D., Chen, Y., Luo, H., Zhao, J.R., ..., & Zhang, X. (2016).
    Somatosensory neuron types identified by high-coverage single-cell RNA-sequencing and functional heterogeneity. Cell research, 26(1), 83-102. [PubMed:26691752] [PMC] [WorldCat] [DOI]
  125. Hirata, Y. (1986).
    [Protein and polypeptide growth factors--epidermal growth factor]. Nihon rinsho. Japanese journal of clinical medicine, 44(1), 7-15. [PubMed:3007810] [WorldCat]
  126. Hanchate, N.K., Kondoh, K., Lu, Z., Kuang, D., Ye, X., Qiu, X., ..., & Buck, L.B. (2015).
    Single-cell transcriptomics reveals receptor transformations during olfactory neurogenesis. Science (New York, N.Y.), 350(6265), 1251-5. [PubMed:26541607] [PMC] [WorldCat] [DOI]
  127. Hockley, J.R.F., Taylor, T.S., Callejo, G., Wilbrey, A.L., Gutteridge, A., Bach, K., ..., & Smith, E.S.J. (2019).
    Single-cell RNAseq reveals seven classes of colonic sensory neuron. Gut, 68(4), 633-644. [PubMed:29483303] [PMC] [WorldCat] [DOI]
  128. Shekhar, K., Lapan, S.W., Whitney, I.E., Tran, N.M., Macosko, E.Z., Kowalczyk, M., ..., & Sanes, J.R. (2016).
    Comprehensive Classification of Retinal Bipolar Neurons by Single-Cell Transcriptomics. Cell, 166(5), 1308-1323.e30. [PubMed:27565351] [PMC] [WorldCat] [DOI]
  129. Konstantinides, N., Kapuralin, K., Fadil, C., Barboza, L., Satija, R., & Desplan, C. (2018).
    Phenotypic Convergence: Distinct Transcription Factors Regulate Common Terminal Features. Cell, 174(3), 622-635.e13. [PubMed:29909983] [PMC] [WorldCat] [DOI]
  130. Buenaventura, D.F., Corseri, A., & Emerson, M.M. (2019).
    Identification of Genes With Enriched Expression in Early Developing Mouse Cone Photoreceptors. Investigative ophthalmology & visual science, 60(8), 2787-2799. [PubMed:31260032] [PMC] [WorldCat] [DOI]
  131. Clark, B.S., Stein-O'Brien, G.L., Shiau, F., Cannon, G.H., Davis-Marcisak, E., Sherman, T., ..., & Blackshaw, S. (2019).
    Single-Cell RNA-Seq Analysis of Retinal Development Identifies NFI Factors as Regulating Mitotic Exit and Late-Born Cell Specification. Neuron, 102(6), 1111-1126.e5. [PubMed:31128945] [PMC] [WorldCat] [DOI]
  132. Rheaume, B.A., Jereen, A., Bolisetty, M., Sajid, M.S., Yang, Y., Renna, K., ..., & Trakhtenberg, E.F. (2018).
    Single cell transcriptome profiling of retinal ganglion cells identifies cellular subtypes. Nature communications, 9(1), 2759. [PubMed:30018341] [PMC] [WorldCat] [DOI]
  133. Peng, Y.R., Shekhar, K., Yan, W., Herrmann, D., Sappington, A., Bryman, G.S., ..., & Sanes, J.R. (2019).
    Molecular Classification and Comparative Taxonomics of Foveal and Peripheral Cells in Primate Retina. Cell, 176(5), 1222-1237.e22. [PubMed:30712875] [PMC] [WorldCat] [DOI]
  134. Collin, J., Queen, R., Zerti, D., Dorgau, B., Hussain, R., Coxhead, J., ..., & Lako, M. (2019).
    Deconstructing Retinal Organoids: Single Cell RNA-Seq Reveals the Cellular Components of Human Pluripotent Stem Cell-Derived Retina. Stem cells (Dayton, Ohio), 37(5), 593-598. [PubMed:30548510] [PMC] [WorldCat] [DOI]
  135. Voigt, A.P., Whitmore, S.S., Flamme-Wiese, M.J., Riker, M.J., Wiley, L.A., Tucker, B.A., ..., & Scheetz, T.E. (2019).
    Molecular characterization of foveal versus peripheral human retina by single-cell RNA sequencing. Experimental eye research, 184, 234-242. [PubMed:31075224] [PMC] [WorldCat] [DOI]
  136. Lo Giudice, Q., Leleu, M., La Manno, G., & Fabre, P.J. (2019).
    Single-cell transcriptional logic of cell-fate specification and axon guidance in early-born retinal neurons. Development (Cambridge, England), 146(17). [PubMed:31399471] [WorldCat] [DOI]
  137. Liang, Q., Dharmat, R., Owen, L., Shakoor, A., Li, Y., Kim, S., ..., & Chen, R. (2019).
    Single-nuclei RNA-seq on human retinal tissue provides improved transcriptome profiling. Nature communications, 10(1), 5743. [PubMed:31848347] [PMC] [WorldCat] [DOI]
  138. Laboissonniere, L.A., Goetz, J.J., Martin, G.M., Bi, R., Lund, T.J.S., Ellson, L., ..., & Trimarchi, J.M. (2019).
    Molecular signatures of retinal ganglion cells revealed through single cell profiling. Scientific reports, 9(1), 15778. [PubMed:31673015] [PMC] [WorldCat] [DOI]
  139. Menon, M., Mohammadi, S., Davila-Velderrain, J., Goods, B.A., Cadwell, T.D., Xing, Y., ..., & Hafler, B.P. (2019).
    Single-cell transcriptomic atlas of the human retina identifies cell types associated with age-related macular degeneration. Nature communications, 10(1), 4902. [PubMed:31653841] [PMC] [WorldCat] [DOI]
  140. Tran, N.M., Shekhar, K., Whitney, I.E., Jacobi, A., Benhar, I., Hong, G., ..., & Sanes, J.R. (2019).
    Single-Cell Profiles of Retinal Ganglion Cells Differing in Resilience to Injury Reveal Neuroprotective Genes. Neuron, 104(6), 1039-1055.e12. [PubMed:31784286] [PMC] [WorldCat] [DOI]
  141. Velasco, S., Kedaigle, A.J., Simmons, S.K., Nash, A., Rocha, M., Quadrato, G., ..., & Arlotta, P. (2019).
    Individual brain organoids reproducibly form cell diversity of the human cerebral cortex. Nature, 570(7762), 523-527. [PubMed:31168097] [PMC] [WorldCat] [DOI]
  142. Yao, Z., Mich, J.K., Ku, S., Menon, V., Krostag, A.R., Martinez, R.A., ..., & Ramanathan, S. (2017).
    A Single-Cell Roadmap of Lineage Bifurcation in Human ESC Models of Embryonic Brain Development. Cell stem cell, 20(1), 120-134. [PubMed:28094016] [PMC] [WorldCat] [DOI]
  143. Close, J.L., Yao, Z., Levi, B.P., Miller, J.A., Bakken, T.E., Menon, V., ..., & Lein, E. (2017).
    Single-Cell Profiling of an In Vitro Model of Human Interneuron Development Reveals Temporal Dynamics of Cell Type Production and Maturation. Neuron, 93(5), 1035-1048.e5. [PubMed:28279351] [PMC] [WorldCat] [DOI]
  144. Bhaduri, A., Andrews, M.G., Mancia Leon, W., Jung, D., Shin, D., Allen, D., ..., & Kriegstein, A.R. (2020).
    Cell stress in cortical organoids impairs molecular subtype specification. Nature, 578(7793), 142-148. [PubMed:31996853] [WorldCat] [DOI]
  145. Nickolls, A.R., Lee, M.M., Espinoza, D.F., Szczot, M., Lam, R.M., Wang, Q., ..., & Bönnemann, C.G. (2020).
    Transcriptional Programming of Human Mechanosensory Neuron Subtypes from Pluripotent Stem Cells. Cell reports, 30(3), 932-946.e7. [PubMed:31968264] [PMC] [WorldCat] [DOI]
  146. Shah, P.T., Stratton, J.A., Stykel, M.G., Abbasi, S., Sharma, S., Mayr, K.A., ..., & Biernaskie, J. (2018).
    Single-Cell Transcriptomics and Fate Mapping of Ependymal Cells Reveals an Absence of Neural Stem Cell Function. Cell, 173(4), 1045-1057.e9. [PubMed:29727663] [WorldCat] [DOI]
  147. Nowakowski, T.J., Bhaduri, A., Pollen, A.A., Alvarado, B., Mostajo-Radji, M.A., Di Lullo, E., ..., & Kriegstein, A.R. (2017).
    Spatiotemporal gene expression trajectories reveal developmental hierarchies of the human cortex. Science (New York, N.Y.), 358(6368), 1318-1323. [PubMed:29217575] [PMC] [WorldCat] [DOI]
  148. Pollen, A.A., Nowakowski, T.J., Chen, J., Retallack, H., Sandoval-Espinosa, C., Nicholas, C.R., ..., & Kriegstein, A.R. (2015).
    Molecular identity of human outer radial glia during cortical development. Cell, 163(1), 55-67. [PubMed:26406371] [PMC] [WorldCat] [DOI]
  149. Zhong, S., Zhang, S., Fan, X., Wu, Q., Yan, L., Dong, J., ..., & Wang, X. (2018).
    A single-cell RNA-seq survey of the developmental landscape of the human prefrontal cortex. Nature, 555(7697), 524-528. [PubMed:29539641] [WorldCat] [DOI]
  150. Johnson, M.B., Wang, P.P., Atabay, K.D., Murphy, E.A., Doan, R.N., Hecht, J.L., & Walsh, C.A. (2015).
    Single-cell analysis reveals transcriptional heterogeneity of neural progenitors in human cortex. Nature neuroscience, 18(5), 637-46. [PubMed:25734491] [PMC] [WorldCat] [DOI]
  151. Yuzwa, S.A., Borrett, M.J., Innes, B.T., Voronova, A., Ketela, T., Kaplan, D.R., ..., & Miller, F.D. (2017).
    Developmental Emergence of Adult Neural Stem Cells as Revealed by Single-Cell Transcriptional Profiling. Cell reports, 21(13), 3970-3986. [PubMed:29281841] [WorldCat] [DOI]
  152. Batiuk, M.Y., Martirosyan, A., Wahis, J., de Vin, F., Marneffe, C., Kusserow, C., ..., & Holt, M.G. (2020).
    Identification of region-specific astrocyte subtypes at single cell resolution. Nature communications, 11(1), 1220. [PubMed:32139688] [PMC] [WorldCat] [DOI]
  153. Bayraktar, O.A., Bartels, T., Holmqvist, S., Kleshchevnikov, V., Martirosyan, A., Polioudakis, D., ..., & Rowitch, D.H. (2020).
    Astrocyte layers in the mammalian cerebral cortex revealed by a single-cell in situ transcriptomic map. Nature neuroscience. [PubMed:32203496] [WorldCat] [DOI]
  154. Marques, S., Zeisel, A., Codeluppi, S., van Bruggen, D., Mendanha Falcão, A., Xiao, L., ..., & Castelo-Branco, G. (2016).
    Oligodendrocyte heterogeneity in the mouse juvenile and adult central nervous system. Science (New York, N.Y.), 352(6291), 1326-1329. [PubMed:27284195] [PMC] [WorldCat] [DOI]
  155. Zhong, S., Zhang, S., Fan, X., Wu, Q., Yan, L., Dong, J., ..., & Wang, X. (2018).
    A single-cell RNA-seq survey of the developmental landscape of the human prefrontal cortex. Nature, 555(7697), 524-528. [PubMed:29539641] [WorldCat] [DOI]
  156. Sousa, C., Golebiewska, A., Poovathingal, S.K., Kaoma, T., Pires-Afonso, Y., Martina, S., ..., & Michelucci, A. (2018).
    Single-cell transcriptomics reveals distinct inflammation-induced microglia signatures. EMBO reports, 19(11). [PubMed:30206190] [PMC] [WorldCat] [DOI]
  157. Matcovitch-Natan, O., Winter, D.R., Giladi, A., Vargas Aguilar, S., Spinrad, A., Sarrazin, S., ..., & Amit, I. (2016).
    Microglia development follows a stepwise program to regulate brain homeostasis. Science (New York, N.Y.), 353(6301), aad8670. [PubMed:27338705] [WorldCat] [DOI]
  158. 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]
  159. Gunner, G., Cheadle, L., Johnson, K.M., Ayata, P., Badimon, A., Mondo, E., ..., & Schafer, D.P. (2019).
    Sensory lesioning induces microglial synapse elimination via ADAM10 and fractalkine signaling. Nature neuroscience, 22(7), 1075-1088. [PubMed:31209379] [PMC] [WorldCat] [DOI]
  160. Mathys, H., Adaikkan, C., Gao, F., Young, J.Z., Manet, E., Hemberg, M., ..., & Tsai, L.H. (2017).
    Temporal Tracking of Microglia Activation in Neurodegeneration at Single-Cell Resolution. Cell reports, 21(2), 366-380. [PubMed:29020624] [PMC] [WorldCat] [DOI]
  161. Geirsdottir, L., David, E., Keren-Shaul, H., Weiner, A., Bohlen, S.C., Neuber, J., ..., & Prinz, M. (2019).
    Cross-Species Single-Cell Analysis Reveals Divergence of the Primate Microglia Program. Cell, 179(7), 1609-1622.e16. [PubMed:31835035] [WorldCat] [DOI]
  162. Van Hove, H., Martens, L., Scheyltjens, I., De Vlaminck, K., Pombo Antunes, A.R., De Prijck, S., ..., & Movahedi, K. (2019).
    A single-cell atlas of mouse brain macrophages reveals unique transcriptional identities shaped by ontogeny and tissue environment. Nature neuroscience, 22(6), 1021-1035. [PubMed:31061494] [WorldCat] [DOI]
  163. Maniatis, S., Äijö, T., Vickovic, S., Braine, C., Kang, K., Mollbrink, A., ..., & Phatnani, H. (2019).
    Spatiotemporal dynamics of molecular pathology in amyotrophic lateral sclerosis. Science (New York, N.Y.), 364(6435), 89-93. [PubMed:30948552] [WorldCat] [DOI]
  164. Mathys, H., Davila-Velderrain, J., Peng, Z., Gao, F., Mohammadi, S., Young, J.Z., ..., & Tsai, L.H. (2019).
    Single-cell transcriptomic analysis of Alzheimer's disease. Nature, 570(7761), 332-337. [PubMed:31042697] [PMC] [WorldCat] [DOI]
  165. Del-Aguila, J.L., Li, Z., Dube, U., Mihindukulasuriya, K.A., Budde, J.P., Fernandez, M.V., ..., & Harari, O. (2019).
    A single-nuclei RNA sequencing study of Mendelian and sporadic AD in the human brain. Alzheimer's research & therapy, 11(1), 71. [PubMed:31399126] [PMC] [WorldCat] [DOI]
  166. Keren-Shaul, H., Spinrad, A., Weiner, A., Matcovitch-Natan, O., Dvir-Szternfeld, R., Ulland, T.K., ..., & Amit, I. (2017).
    A Unique Microglia Type Associated with Restricting Development of Alzheimer's Disease. Cell, 169(7), 1276-1290.e17. [PubMed:28602351] [WorldCat] [DOI]
  167. Velmeshev, D., Schirmer, L., Jung, D., Haeussler, M., Perez, Y., Mayer, S., ..., & Kriegstein, A.R. (2019).
    Single-cell genomics identifies cell type-specific molecular changes in autism. Science (New York, N.Y.), 364(6441), 685-689. [PubMed:31097668] [WorldCat] [DOI]
  168. Renthal, W., Boxer, L.D., Hrvatin, S., Li, E., Silberfeld, A., Nagy, M.A., ..., & Greenberg, M.E. (2018).
    Characterization of human mosaic Rett syndrome brain tissue by single-nucleus RNA sequencing. Nature neuroscience, 21(12), 1670-1679. [PubMed:30455458] [PMC] [WorldCat] [DOI]
  169. Sainio, M.T., Ylikallio, E., Mäenpää, L., Lahtela, J., Mattila, P., Auranen, M., ..., & Tyynismaa, H. (2018).
    Absence of NEFL in patient-specific neurons in early-onset Charcot-Marie-Tooth neuropathy. Neurology. Genetics, 4(3), e244. [PubMed:29888333] [PMC] [WorldCat] [DOI]
  170. Lang, C., Campbell, K.R., Ryan, B.J., Carling, P., Attar, M., Vowles, J., ..., & Wade-Martins, R. (2019).
    Single-Cell Sequencing of iPSC-Dopamine Neurons Reconstructs Disease Progression and Identifies HDAC4 as a Regulator of Parkinson Cell Phenotypes. Cell stem cell, 24(1), 93-106.e6. [PubMed:30503143] [PMC] [WorldCat] [DOI]
  171. Neftel, C., Laffy, J., Filbin, M.G., Hara, T., Shore, M.E., Rahme, G.J., ..., & Suvà, M.L. (2019).
    An Integrative Model of Cellular States, Plasticity, and Genetics for Glioblastoma. Cell, 178(4), 835-849.e21. [PubMed:31327527] [PMC] [WorldCat] [DOI]
  172. Venteicher, A.S., Tirosh, I., Hebert, C., Yizhak, K., Neftel, C., Filbin, M.G., ..., & Suvà, M.L. (2017).
    Decoupling genetics, lineages, and microenvironment in IDH-mutant gliomas by single-cell RNA-seq. Science (New York, N.Y.), 355(6332). [PubMed:28360267] [PMC] [WorldCat] [DOI]
  173. Horie, R., Hazbun, A., Chen, K., Cao, C., Levine, M., & Horie, T. (2018).
    Shared evolutionary origin of vertebrate neural crest and cranial placodes. Nature, 560(7717), 228-232. [PubMed:30069052] [PMC] [WorldCat] [DOI]
  174. Horie, T., Horie, R., Chen, K., Cao, C., Nakagawa, M., Kusakabe, T.G., ..., & Levine, M. (2018).
    Regulatory cocktail for dopaminergic neurons in a protovertebrate identified by whole-embryo single-cell transcriptomics. Genes & development, 32(19-20), 1297-1302. [PubMed:30228204] [PMC] [WorldCat] [DOI]
  175. Davie, K., Janssens, J., Koldere, D., De Waegeneer, M., Pech, U., Kreft, Ł., ..., & Aerts, S. (2018).
    A Single-Cell Transcriptome Atlas of the Aging Drosophila Brain. Cell, 174(4), 982-998.e20. [PubMed:29909982] [PMC] [WorldCat] [DOI]
  176. Li, H., Horns, F., Wu, B., Xie, Q., Li, J., Li, T., ..., & Luo, L. (2017).
    Classifying Drosophila Olfactory Projection Neuron Subtypes by Single-Cell RNA Sequencing. Cell, 171(5), 1206-1220.e22. [PubMed:29149607] [PMC] [WorldCat] [DOI]
  177. Bates, A.S., Janssens, J., Jefferis, G.S., & Aerts, S. (2019).
    Neuronal cell types in the fly: single-cell anatomy meets single-cell genomics. Current opinion in neurobiology, 56, 125-134. [PubMed:30703584] [WorldCat] [DOI]
  178. Konstantinides, N., Kapuralin, K., Fadil, C., Barboza, L., Satija, R., & Desplan, C. (2018).
    Phenotypic Convergence: Distinct Transcription Factors Regulate Common Terminal Features. Cell, 174(3), 622-635.e13. [PubMed:29909983] [PMC] [WorldCat] [DOI]
  179. Cosacak, M.I., Bhattarai, P., Reinhardt, S., Petzold, A., Dahl, A., Zhang, Y., & Kizil, C. (2019).
    Single-Cell Transcriptomics Analyses of Neural Stem Cell Heterogeneity and Contextual Plasticity in a Zebrafish Brain Model of Amyloid Toxicity. Cell reports, 27(4), 1307-1318.e3. [PubMed:31018142] [WorldCat] [DOI]
  180. Thyme, S.B., Pieper, L.M., Li, E.H., Pandey, S., Wang, Y., Morris, N.S., ..., & Schier, A.F. (2019).
    Phenotypic Landscape of Schizophrenia-Associated Genes Defines Candidates and Their Shared Functions. Cell, 177(2), 478-491.e20. [PubMed:30929901] [PMC] [WorldCat] [DOI]
  181. Tosches, M.A., Yamawaki, T.M., Naumann, R.K., Jacobi, A.A., Tushev, G., & Laurent, G. (2018).
    Evolution of pallium, hippocampus, and cortical cell types revealed by single-cell transcriptomics in reptiles. Science (New York, N.Y.), 360(6391), 881-888. [PubMed:29724907] [WorldCat] [DOI]
  182. Marchetto, M.C., Hrvoj-Mihic, B., Kerman, B.E., Yu, D.X., Vadodaria, K.C., Linker, S.B., ..., & Gage, F.H. (2019).
    Species-specific maturation profiles of human, chimpanzee and bonobo neural cells. eLife, 8. [PubMed:30730291] [PMC] [WorldCat] [DOI]
  183. Kanton, S., Boyle, M.J., He, Z., Santel, M., Weigert, A., Sanchís-Calleja, F., ..., & Camp, J.G. (2019).
    Organoid single-cell genomic atlas uncovers human-specific features of brain development. Nature, 574(7778), 418-422. [PubMed:31619793] [WorldCat] [DOI]
  184. Huerta-Cepas, J., Szklarczyk, D., Heller, D., Hernández-Plaza, A., Forslund, S.K., Cook, H., ..., & Bork, P. (2019).
    eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic acids research, 47(D1), D309-D314. [PubMed:30418610] [PMC] [WorldCat] [DOI]
  185. Keil, J.M., Qalieh, A., & Kwan, K.Y. (2018).
    Brain Transcriptome Databases: A User's Guide. The Journal of neuroscience : the official journal of the Society for Neuroscience, 38(10), 2399-2412. [PubMed:29437890] [PMC] [WorldCat] [DOI]
  186. Ecker, J.R., Geschwind, D.H., Kriegstein, A.R., Ngai, J., Osten, P., Polioudakis, D., ..., & Zeng, H. (2017).
    The BRAIN Initiative Cell Census Consortium: Lessons Learned toward Generating a Comprehensive Brain Cell Atlas. Neuron, 96(3), 542-557. [PubMed:29096072] [PMC] [WorldCat] [DOI]
  187. Welch, J.D., Kozareva, V., Ferreira, A., Vanderburg, C., Martin, C., & Macosko, E.Z. (2019).
    Single-Cell Multi-omic Integration Compares and Contrasts Features of Brain Cell Identity. Cell, 177(7), 1873-1887.e17. [PubMed:31178122] [PMC] [WorldCat] [DOI]
  188. Korsunsky, I., Millard, N., Fan, J., Slowikowski, K., Zhang, F., Wei, K., ..., & Raychaudhuri, S. (2019).
    Fast, sensitive and accurate integration of single-cell data with Harmony. Nature methods, 16(12), 1289-1296. [PubMed:31740819] [PMC] [WorldCat] [DOI]
  189. Crow, M., Paul, A., Ballouz, S., Huang, Z.J., & Gillis, J. (2018).
    Characterizing the replicability of cell types defined by single cell RNA-sequencing data using MetaNeighbor. Nature communications, 9(1), 884. [PubMed:29491377] [PMC] [WorldCat] [DOI]
  190. Rodriques, S.G., Stickels, R.R., Goeva, A., Martin, C.A., Murray, E., Vanderburg, C.R., ..., & Macosko, E.Z. (2019).
    Slide-seq: A scalable technology for measuring genome-wide expression at high spatial resolution. Science (New York, N.Y.), 363(6434), 1463-1467. [PubMed:30923225] [PMC] [WorldCat] [DOI]
  191. Codeluppi, S., Borm, L.E., Zeisel, A., La Manno, G., van Lunteren, J.A., Svensson, C.I., & Linnarsson, S. (2018).
    Spatial organization of the somatosensory cortex revealed by osmFISH. Nature methods, 15(11), 932-935. [PubMed:30377364] [WorldCat] [DOI]
  192. Wang, X., Allen, W.E., Wright, M.A., Sylwestrak, E.L., Samusik, N., Vesuna, S., ..., & Deisseroth, K. (2018).
    Three-dimensional intact-tissue sequencing of single-cell transcriptional states. Science (New York, N.Y.), 361(6400). [PubMed:29930089] [PMC] [WorldCat] [DOI]
  193. Shah, S., Lubeck, E., Zhou, W., & Cai, L. (2016).
    In Situ Transcription Profiling of Single Cells Reveals Spatial Organization of Cells in the Mouse Hippocampus. Neuron, 92(2), 342-357. [PubMed:27764670] [PMC] [WorldCat] [DOI]
  194. Edsgärd, D., Johnsson, P., & Sandberg, R. (2018).
    Identification of spatial expression trends in single-cell gene expression data. Nature methods, 15(5), 339-342. [PubMed:29553578] [PMC] [WorldCat] [DOI]
  195. Svensson, V., Teichmann, S.A., & Stegle, O. (2018).
    SpatialDE: identification of spatially variable genes. Nature methods, 15(5), 343-346. [PubMed:29553579] [PMC] [WorldCat] [DOI]
  196. Chen, K.H., Boettiger, A.N., Moffitt, J.R., Wang, S., & Zhuang, X. (2015).
    RNA imaging. Spatially resolved, highly multiplexed RNA profiling in single cells. Science (New York, N.Y.), 348(6233), aaa6090. [PubMed:25858977] [PMC] [WorldCat] [DOI]
  197. Coskun, A.F., & Cai, L. (2016).
    Dense transcript profiling in single cells by image correlation decoding. Nature methods, 13(8), 657-60. [PubMed:27271198] [PMC] [WorldCat] [DOI]
  198. Zhu, C., Preissl, S., & Ren, B. (2020).
    Single-cell multimodal omics: the power of many. Nature methods, 17(1), 11-14. [PubMed:31907462] [WorldCat] [DOI]
  199. Stuart, T., & Satija, R. (2019).
    Integrative single-cell analysis. Nature reviews. Genetics, 20(5), 257-272. [PubMed:30696980] [WorldCat] [DOI]
  200. Fuzik, J., Zeisel, A., Máté, Z., Calvigioni, D., Yanagawa, Y., Szabó, G., ..., & Harkany, T. (2016).
    Integration of electrophysiological recordings with single-cell RNA-seq data identifies neuronal subtypes. Nature biotechnology, 34(2), 175-183. [PubMed:26689544] [PMC] [WorldCat] [DOI]
  201. Cadwell, C.R., Palasantza, A., Jiang, X., Berens, P., Deng, Q., Yilmaz, M., ..., & Tolias, A.S. (2016).
    Electrophysiological, transcriptomic and morphologic profiling of single neurons using Patch-seq. Nature biotechnology, 34(2), 199-203. [PubMed:26689543] [PMC] [WorldCat] [DOI]
  202. Stoeckius, M., Hafemeister, C., Stephenson, W., Houck-Loomis, B., Chattopadhyay, P.K., Swerdlow, H., ..., & Smibert, P. (2017).
    Simultaneous epitope and transcriptome measurement in single cells. Nature methods, 14(9), 865-868. [PubMed:28759029] [PMC] [WorldCat] [DOI]
  203. Peterson, V.M., Zhang, K.X., Kumar, N., Wong, J., Li, L., Wilson, D.C., ..., & Klappenbach, J.A. (2017).
    Multiplexed quantification of proteins and transcripts in single cells. Nature biotechnology, 35(10), 936-939. [PubMed:28854175] [WorldCat] [DOI]
  204. Chen, X., Sun, Y.C., Zhan, H., Kebschull, J.M., Fischer, S., Matho, K., ..., & Zador, A.M. (2019).
    High-Throughput Mapping of Long-Range Neuronal Projection Using In Situ Sequencing. Cell, 179(3), 772-786.e19. [PubMed:31626774] [WorldCat] [DOI]