2019 5/23 tweet追記
2019 9/26 動画追加
2019 11/30 ツイート追記
2019 12/22 統合TVリンク追加
2020 2/2 ８章補足資料リンク追加
RNAシークエンシング（RNA-Seq）[1]は、ゲノムワイドな発現解析のための日常的な技術となった。ますます低コストで、ライブラリー構築およびシーケンシングはしばしば標準的なプロトコルに従って実施することができる。多くの研究者、特にバイオインフォマティクスの経験がない研究者にとって、この技術を十分に活用するためのボトルネックは、発現プロファイルを実用的な洞察に変換する方法である。典型的な分析ワークフローは多くのステップを含み、各ステップは異なるツールを必要とする。これらのツールを正しく学習し、調整し、繋ぐのは時間がかかる。もう1つの障害は、多様な種類の遺伝子IDを持つ散在したアノテーションデータベースである。これらの問題を軽減するため、研究者は、RNA-Seqデータを分析するのに必要な時間と労力を大幅に削減できるアプリケーションの開発を目指している。
RNA-Seqデータ解析は、クオリティ管理、前処理、マッピング、そしてrawシークエンスリードのsummarizingから始まる。これらのステップは、従来のTuxedo Suite [ref.2、3]またはより高速でアライメントフリーの定量方法[ref.4、5]のいずれかを使用して完了したと想定する。これらのツールは、スタンドアロンでも、GenePattern [ref.6]、Galaxy [ref.7]、CyVerse [ref.8]（紹介）のようなプラットフォームでも使用できる。
リードマッピングの後、よくあるのは、遺伝子レベルのリードカウントまたは正規化された発現レベルの行列（Fragments Per Kilobase Million、またはFPKM）を取得することである。Rは、そのDNAマイクロアレイデータのような表形式データの強力な可視化と統計分析ツールである。さらに、differentially expressed genes（DEG）および変更されたパスウェイを同定するために、多くの専用のRおよびBioconductor [ref.9]パッケージが開発されてきた。 DESeq2 [ref.10]のようないくつかのパッケージは、特にリードカウントの統計的モデリングのために開発され、広く使われている。しかし、これらのパッケージは時間がかかり、コーディング経験がないと研究者にとって手が届かないことさえある。
発現データを分析するいくつかのウェブアプリケーションが開発されている（論文 表１）。 STARTアプリ（shinnyトランスクリプトーム解析リソースツール）は、階層的クラスタリング、主成分分析（PCA）、遺伝子レベルのボックスプロット、およびDEGを実行するshinnyアプリである[ref.11]。別の同様のツールであるDegust [ref.12]はEdgeR [ref.13]またはlimma-voom [ref.14]を使って発現解析を実行し、インタラクティブに結果をプロットできる。他のツールには、Sleuth [ref.15]とShinyNGS [ref.16]がある。shinnyでないアプリケーションも開発された。これにはDEIVA [ref.17]とVisRseq [ref.18]が含まれる。いくつかのツールはパスウェイ分析のいくつかの能力を組み込んでいる。定量された発現データについては、ASAP（自動化シングルセル分析パイプライン）[ref.19]により、Gene Ontology（GO）[ref.20]およびKEGG [ref.21]データベースに基づいて正規化、フィルタリング、クラスタリング、およびエンリッチメント分析を行うことができる。 EXPath Tool [ref.22]を使用すると、ユーザーはパスウェイ探索、GOエンリッチメント解析および共発現解析を実行できる。 IRISなどの他のいくつかのShinyベースのツール[ref.23]も開発されている。過去数年間におけるこれらのツールの開発は、RNA-Seqデータの解釈を容易にした。
本研究では、（1）広範囲にわたる自動遺伝子ID変換、（2）植物と動物の両方に対する包括的な遺伝子アノテーションとパスウェイデータベース、（3）いくつかの方法、詳細なEDAおよびパスウェイ解析、（4）アプリケーションプログラミングインタフェース（API）を介したKEGG [ref.21]、STRING-db [ref.24]などのWebサービスへのアクセス、および（5）スタンドアロン解析用のRスクリプトの生成による再現性の向上、を含むwebアプリケーションを開発する。iDEPを使用して、2つのサンプルデータセットを分析し、論文表1と図1以外のすべての図と表を生成した。(以下略)
iDEP overview
https://idepsite.wordpress.com/
FAQ
https://idepsite.wordpress.com/faq/
データフォーマット
Data format – iDEP: Gain Insights from RNA-seq
Our iDEP paper is finally out. In addition to the source code in Github, all the supporting databases are available on Zenodo.
— Steven Ge (@StevenXGe) December 20, 2018
iDEP: an integrated web application for differential expression and pathway analysis of RNA-Seq data https://t.co/9D1J1eX04A #bmcbioinformatics
追記
Big update to iDEP!
— Steven Ge (@StevenXGe) May 21, 2019
Besides the 315 genomes in Ensembl, we incorporated the annotation for 115 archaeal, 1678 bacterial, and 238 eukaryotic species from STRING-db. iDEP is designed to help smaller labs making sense of transcriptomic data.
Please help spread the word! pic.twitter.com/6rw0X7tFNw
Re-analyze public RNA-seq data can be done easily with interactive tools like iDEP and uniformly processed read counts data from ARCHS4 and DEE2. See how I did it in under 12 minutes.
— Steven Ge (@StevenXGe) September 23, 2019
No more excuses!
We just need good questions and critical thinkers!https://t.co/CZB28e0U0Y
iDEPに関するツイート
使い方
local環境でも簡単に実行できるが、ここではidepサイトにアクセスして動作を確認する。
http://bioinformatics.sdstate.edu/idep/
デモデータで読み込むフォーマットを確認する。左上のClick Hereをクリック。
読み込まれた。
１列目がgene IDになる。Ensembl IDを使うが、一般的なgene IDなら自動認識してEnsembl IDに内部変換してくれる。１行目にはサンプル名を記載する。Control、 TreatmentA、TreatmentBの３条件でそれぞれbiological/technical replicatesが２つずつあるなら、Ctrl_1, Ctrl_2, TrtA_1, TrtA_2, TrtB_1, TrtB_2のように記載する。これでreplicatesとして認識される。
実際に読み込めるデータは、RNA seqのリードカウントデータ（正規化前）、またはFPKMで正規化したデータ、マイクロアレイのデータなどになる。
他の入力例（データフォーマットの解説より）
Gene expression matrix ex.1
Gene expression matrix ex.2
fold-change and P values
Public RNA-Seq and ChIP-Seq dataでは、ARCHS4でマイニングされた7000以上のSRAサンプルの定量データ（CSVファイル）（paper）を解析できる。
DownloadからARCHS4でリードカウントされたデータのCSVをダウンロードして使ってみる。1行目はこのようにした。
左上のリードカウントからCSVをアップロード。
アップロード完了。フォーマットにエラーがあると赤字のメッセージが出る。gene IDは、Ensembl IDに変換され、以後の解析が行われる。IDはEnsembl release 92に基づいており、220の生物種に対応している（link）（2018/12/29現在）。
１、Pre-process（説明）: データの正規化
replicatesとして認識されたサンプルはグラフで同じ色がアサインされている。動作は非常に俊敏で、データのアップロード完了から１０秒くらいで結果は可視化される。
左上のウィンドウからパラメータを変更できる。デフォルトでは1サンプルで0.5 counts per million (CPM) 以上の遺伝子が分析対象となる。 pseudo countとして追加される値はデフォルト4になっている。
２、Heatmap（説明）階層的クラスタリングおよびデンドログラムとヒートマップによる可視化
defaultでは全サンプル間のSD top1000が対象となる。左上のgene SD distributionで可視化。
左上のinteractive heatmapボタンからスケールをリアルタイムに変更できる。50遺伝子に絞った。
グラフはマウスで操作できる。
図と対象遺伝子のCSVは左下からダウンロードできる。
３、k-means（説明）非階層クラスタリングおよびデンドログラムとヒートマップによる可視化
defaultでは全サンプル間のSD top 2000がクラスター分析対象となる。クラスター数のdefaultは4。
t-SNEによる次元圧縮結果
エンリッチされた転写因子（TF）結合モチーフ
エンリッチされたpathway
デンドログラムで可視化（左のVisualize enrichmentボタンより）
４、PCA（説明）主成分分析
５、DEG（説明）発現変動遺伝子の検出
デフォルトではDESeq2が使われる。FDR cutoffは0.1、Min fold changeは2。３つ以上のグループがある場合、全てのペアワイズ比較が実施される。これまでと同様に、結果は左のメニューからダウンロードできる。
ベン図（Venn Diagramボタンより実行できる）
３サンプルの場合
DEGのタブは２つある。DEG2のタブでは、３群以上のデータでも、組み合わせを選び、それぞれ２群比較を実行できる。結果はヒートマップ、MA-plot、Scatter plotなどで可視化される。また、検定されてエンリッチされていると判定されたpathwayが下にp-value付きで表示される。
操作パネル左下にはShinyGO（prepirnt）へのリンクもある。
ShinyGOはユーザーが指定した遺伝子リストを元に、エンリッチされた系を可視化する。機能はiDEPとかなり重複しているが、iDEPにないツールもある。
２群間比較結果の可視化
Scatter plot
MA plot
Volcano plot
６、Pathway（説明）GOエンリッチメント解析
pathway解析はDEGで指定してフィルタリングした条件ではなく、全データからDESeq2 / limmaで出力したfold change値を使って行われることに注意する。有意水準を決めるFDR cutoffはpathway解析パネルにもあるので、こちらで厳密さは調節する。
有名な、いくつかのGOエンリッチメント解析ツールを利用できる。
GAGE（Generally Applicable Gene-set Enrichment ）（pubmed）（遅い）
GSEA (Gene Set Enrichment Analysis) (preranked fgsea)（link）（解説HP）
PGSEA: PAGE (Parametric Analysis of Gene Set Enrichment) の実装 （PDF link）２群間比較
PGSEA w/all sample: 全サンプル間比較
上記の分析は全てbuint-inのデータを使っているが、最後のReactomePA (Reactome Pathway Analysis) は、 Reactomeの遺伝子セットデータベースを使う（Reactome 統合TV解説）。
Reactome (HP) は歴史あるデータベースで、ヒト以外のモデル生物種にも対応している。ペーパーは多数出ている (Google scholar検索結果)。
たくさんのgenesetを利用できる。一覧はマニュアル参照（link）。KEGGに切り替えてエンリッチと判定されたpathwayを可視化する。
Down regulationが35で、NES（Normalised enrichment score for the given gene set）がもっとも低かったHedgehog signaling pathway
可視化するKEGG pathwayは図の上のメニューから選ぶ。
７、Genomes DGEのGenome上の位置を可視化する
８、Biclustering（説明）サンプル数の多いラージデータセットの解析で（>10）、相関関係にある遺伝子のグループを検出する
biclust R、BCCC、BCXmotifs、BCPlaid、BCSpectral、BCBimax、BCQuest、QUBICなどのRパッケージを利用できる。解説は上の説明リンク参照。
９、Network（説明）共発現解析およびネットワークの可視化。かなり巨大なデータセットで（>15）相関関係がありそうな遺伝子セットを探すために使う（小さなデータセットではクラスタリング解析を行う）。
ネットワークは巨大になるので、全ネットワークの可視化には、localでcytoscape等を使うことが推奨されている。
感想
インタラクティブにパラメータを変更して、結果を見比べながら進められる素晴らしいツールですね。salmonやkallistoなどと組み合わせれば、rawデータが手に入ってから数時間以内にラフな結果は出せるのではないでしょうか（*1）。もちろん、方法について理解していないと使いこなせない訳で、誰でもRNA seq解析できるとは言いません。しかし、RNA seq解析の敷居が下がっているのは確かで、これだけ完成度の高いアプリケーションが出てくると、リードカウントはCyverse（紹介）、DEG検出などはiDEP、と使い分けることで、コンソールが一切使えない環境でも結果を出せてしまえますね。
2019 11/30追記
羊土社さんからRNA-Seqデータ解析の専門書が発売されました。その中の１章で、iDEPの使い方を紹介しています。よろしければ手にとってご確認下さい。どの章もたいへん見応えのある内容になっています。
【新刊】『RNA-Seqデータ解析 WETラボのための鉄板レシピ』発売！ 食材調達（データ入手）から盛付（論文投稿）までシェフ直伝！「内容の更新情報が届く」という驚きの付録つき！ #RNAseqRecipe https://t.co/bPPUTmA2tt
— 羊土社 実験医学 (@Yodosha_EM) 2019年11月29日
実験医学別冊：RNA-Seqデータ解析 WETラボのための鉄板レシピ - 羊土社
８章補足資料
2019 12/22
ShinyGO
2020 4/21
統合TVでBioJupiesの使い方が紹介されています。
BioJupiesを使ってウェブブラウザ上でRNA-seqデータ解析を行う
2020 5/26
引用
iDEP: an integrated web application for differential expression and pathway analysis of RNA-Seq data
Steven Xijin Ge, Eun Wo Son, Runan Yao
BMC Bioinformatics 2018 19:534
参考
slideshare バイオインフォマティクスによる遺伝子発現解析
https://www.slideshare.net/sesejun/ss-24923282
*1
75-bpのシングルエンドシーケンシングと組み合わせれば、ラン開始から24時間で解析結果まで得られることになる。
pre-processing
mapping
read count
関連ツール