数据科学工厂 发表于 2024-6-29 14:44:06

ATAC-seq分析:教程简介(1)

## 简介

[本课程](https://rockefelleruniversity.github.io/RU_ATACseq/ "Source")介绍 `Bioconductor` 中的 `ATACseq` 分析。

该课程由 2 个部分组成。这将引导您完成正常 `ATACseq` 分析工作流程的每个步骤。它涵盖比对、`QC`、`peak calling`、基因组富集测试、基序富集和差异可及性测试。

## 环境准备

### IGV

IGV 可以从 `BROAD` 网站安装。 》 https://www.broadinstitute.org/igv/

### MACS2

(https://macs3-project.github.io/MACS/ "MACS") 没有 R 包,但 MACS2 可在适用于 Linux 或 MacOS 的 Anaconda 包存储库中找到。安装 MACS2 的最简单方法是使用 R 包 Herper。 Herper 允许您从 R 中管理和安装 Anaconda 包。

```R
BiocManager::install("Herper")
library(Herper)
```

安装 Herper 后,您可以使用 install_CondaTools 函数安装 MACS2。在幕后,Herper 将安装最小版本的 conda(称为 miniconda),然后创建一个新环境来安装 MACS2。当您运行该函数时,它会打印出 MACS2 的安装位置。

env 参数是您要为创建的环境指定的名称。 pathToMiniConda 指定您要安装 Miniconda 的位置,以及所有 conda 工具(如 MACS2)。

```R
install_CondaTools(tools="macs2", env="PeakCalling_analysis", pathToMiniConda="/path/to/install")
```

### R



### RStudio



### 包

- 课程包

```R
install.packages('BiocManager')
BiocManager::install('RockefellerUniversity/RU_ATACseq',subdir='atacseq')
```

- 来自 CRAN 和 Bioconductor

```R
install.packages('BiocManager')
BiocManager::install('methods')
BiocManager::install('ggplot2')
BiocManager::install('rmarkdown')
BiocManager::install('ShortRead')
BiocManager::install('ashr')
BiocManager::install('ChIPQC')
BiocManager::install('DiffBind')
BiocManager::install('BSgenome.Hsapiens.UCSC.hg19')
BiocManager::install('Rsubread')
BiocManager::install('Rbowtie2')
BiocManager::install('R.utils')
BiocManager::install('Rsamtools')
BiocManager::install('BSgenome.Hsapiens.UCSC.hg38')
BiocManager::install('rtracklayer')
BiocManager::install('ChIPseeker')
BiocManager::install('soGGi')
BiocManager::install('GenomicAlignments')
BiocManager::install('TxDb.Hsapiens.UCSC.hg19.knownGene')
BiocManager::install('DESeq2')
BiocManager::install('BSgenome.Mmusculus.UCSC.mm10')
BiocManager::install('TxDb.Hsapiens.UCSC.hg38.knownGene')
BiocManager::install('tracktables')
BiocManager::install('clusterProfiler')
BiocManager::install('TxDb.Mmusculus.UCSC.mm10.knownGene')
BiocManager::install('devtools')
BiocManager::install('tidyr')
BiocManager::install('DT')
BiocManager::install('dplyr')
BiocManager::install('rGREAT')
BiocManager::install('MotifDb')
BiocManager::install('Biostrings')
BiocManager::install('GenomicRanges')
BiocManager::install('pheatmap')
BiocManager::install('universalmotif')
BiocManager::install('seqLogo')
BiocManager::install('org.Mm.eg.db')
BiocManager::install('ATACseqQC')
BiocManager::install('JASPAR2020')
BiocManager::install('motifmatchr')
BiocManager::install('chromVAR')
BiocManager::install('ggseqlogo')
BiocManager::install('TFBSTools')
BiocManager::install('motifStack')
BiocManager::install('knitr')
BiocManager::install('testthat')
BiocManager::install('yaml')
```

## 内容

### Part_1

本节介绍Bioconductor Session部分对ATACseq数据的分析:

- 在 R 中预处理 ATACseq 数据
- 数据比对
- 为可视化创建 bigWig

!(data/attachment/forum/plugin_zhanmishu_markdown/202406/ba886d45efd2983072106418562231b3_1719643443_6398.png)

### Part_2

本节演示如何使用 `ATACseq` 数据评估可访问性的全局变化。会话部分:

- 在 R 中注释 ATACseq 数据
- 绘制无核小体和单核小体信号
- 绘制 DNA 结合蛋白周围的切割位点

!(data/attachment/forum/plugin_zhanmishu_markdown/202406/8c99e51c4ccd676e52a6099e183ac0ab_1719643443_2285.png)

### Part_3

本节演示如何评估 ATAC-seq 数据中的基序。

- 从数据库中检索 `motifs`。
- 绘制 `motifs`。
- 识别已知的 `motifs`。

!(data/attachment/forum/plugin_zhanmishu_markdown/202406/761aa72a3eec2e589ca9db5794d657fd_1719643443_5000.png)
页: [1]
查看完整版本: ATAC-seq分析:教程简介(1)