Seuratを用いたscATAC-seqの各クラスターのピークを重ね合わせる方法
【目的】
Single-cell ATAC-seq (もしくはsingle-cell RNA-seq)の各クラスターのピークを1つのラインに重ね合わせた図を作りたいです。
【環境】
Windows 10
scATAC-seq data set: 10x Genomicsのサービスを利用
Viewer: Loupe browser 5.0.0
R: version 3.5.1
Seurat: version 4.0
【質問】
バイオインフォマティクス初心者です。Single-cell ATAC-seqのデータから論文用にFigureを作ろうとしています。
特定の遺伝子領域上のピークを各クラスターごとでLoupe browserにて表示することはできるのですが、それらクラスターのピークを重ね合わせて1つの図にする方法を教えて頂けませんでしょうか。
一応、Loupe browserで表示したピークを各クラスターごとでスクリーンショットを撮り、Powerpointを使って手動で重ね合わせ、透明度をいじることで同じような図を作ることはできます(添付図)。しかし、Loupe browserで表示したピーク図の背景が透明ではなく白色なのもあり、ピークの輪郭・色が薄くなってしまい、あまり綺麗な図になりません。別の方法が無いか問い合わせたところ、そういった3次解析はViewerのLoupe browserではなく、SeuratなどのRパッケージを使うのが一般的だとお聞きしました。
しかし、自分なりに調べても、ピークの図を重ねる方法を見つけることができず、ここに質問することにしました。
素人質問で申し訳ありませんが、お答え願えますと大変幸いです。
【追記】
Seurat (version 4.0)のパッケージはSatija labのウェブページの案内に沿って、既にRにインストール済みです。
https://satijalab.org/seurat/install.html