转录组自动化分析流程搭建及使用

白墨石

发布于 2021-03-24 14:06:23

9180

发布于 2021-03-24 14:06:23

文章被收录于专栏：生信情报站生信情报站

这次分析流程搭建使用基于Nextflow 的 nf-core，该工具可以实现自动化的转录组上游分析。

官网：https://nf-co.re/rnaseq

GitHub：https://github.com/nf-core/rnaseq

安装 nf-core rnaseq

可以使用Git clone，也可以下载好解压到流程目录

安装Nextflow

curl -s https://get.nextflow.io | bash

检测版本是否符合nf-core使用，可以升级

nextflow self-update

安装aws

curl "https://awscli.amazonaws.com/awscli-exe-linux-x86_64.zip" -o "awscliv2.zip"
unzip awscliv2.zip
sudo ./aws/install

下载参考基因组

aws s3 --no-sign-request --region eu-west-1 sync s3://ngi-igenomes/igenomes/Homo_sapiens/Ensembl/GRCh37/Annotation/Genes/ ./references/Homo_sapiens/Ensembl/GRCh37/Annotation/Genes/ --exclude "*" --include "genes.gtf"
aws s3 --no-sign-request --region eu-west-1 sync s3://ngi-igenomes/igenomes/Homo_sapiens/Ensembl/GRCh37/Sequence/WholeGenomeFasta/ ./references/Homo_sapiens/Ensembl/GRCh37/Sequence/WholeGenomeFasta/
aws s3 --no-sign-request --region eu-west-1 sync s3://ngi-igenomes/igenomes/Homo_sapiens/Ensembl/GRCh37/Sequence/STARIndex/ ./references/Homo_sapiens/Ensembl/GRCh37/Sequence/STARIndex/
aws s3 --no-sign-request --region eu-west-1 sync s3://ngi-igenomes/igenomes/Homo_sapiens/Ensembl/GRCh37/Sequence/BWAIndex/ ./references/Homo_sapiens/Ensembl/GRCh37/Sequence/BWAIndex/
aws s3 --no-sign-request --region eu-west-1 sync s3://ngi-igenomes/igenomes/Homo_sapiens/Ensembl/GRCh37/Sequence/Bowtie2Index/ ./references/Homo_sapiens/Ensembl/GRCh37/Sequence/Bowtie2Index/
aws s3 --no-sign-request --region eu-west-1 sync s3://ngi-igenomes/igenomes/Homo_sapiens/Ensembl/GRCh37/Annotation/Genes/ ./references/Homo_sapiens/Ensembl/GRCh37/Annotation/Genes/ --exclude "*" --include "genes.bed"

https://ewels.github.io/AWS-iGenomes/

测试数据

数据来源GSE101571

构建测试数据信息表，rnaseq-test.csv

group,replicate,fastq_1,fastq_2,strandedness
2cell,1,/data/baimoc/data/rnaseq-test/SRR5837392_1.fastq.gz,/data/baimoc/data/rnaseq-test/SRR5837392_2.fastq.gz,unstranded
2cell,1,/data/baimoc/data/rnaseq-test/SRR5837393_1.fastq.gz,/data/baimoc/data/rnaseq-test/SRR5837393_2.fastq.gz,unstranded
8cell,1,/data/baimoc/data/rnaseq-test/SRR5837402_1.fastq.gz,/data/baimoc/data/rnaseq-test/SRR5837402_2.fastq.gz,unstranded
8cell,1,/data/baimoc/data/rnaseq-test/SRR5837403_1.fastq.gz,/data/baimoc/data/rnaseq-test/SRR5837403_2.fastq.gz,unstranded

项目目录

启动流程

../../nextflow run ../../rnaseq --input /data/baimoc/data/rnaseq-test/rnaseq-test.csv --genome GRCh37 --igenomes_base /data/baimoc/references/ -profile docker

本文参与腾讯云自媒体分享计划，分享自作者个人站点/博客。

原始发表：2021-03-22 ，如有侵权请联系 cloudcommunity@tencent.com 删除

https