禾谷镰刀菌转录组分析--准备数据

1、准备数据:测序数据和参考基因组
测序数据:NCBI中SRA数据库,Accession: PRJNA522013

ASPERA,类似下载加速器,要安装在个人目录才能使用

wget http://download.asperasoft.com/download/sw/connect/3.7.4/aspera-connect-3.7.4.147727-linux-64.tar.gz
tar zxvf aspera-connect-3.7.4.147727-linux-64.tar.gz
#安装
bash aspera-connect-3.7.4.147727-linux-64.sh
# 然后cd到根目录下看看是不是存在了.aspera文件夹,有的话表示安装成功
cd && ls -a
# 将aspera软件加入环境变量,并激活
echo 'export PATH=~/.aspera/connect/bin:$PATH' >> ~/.bashrc
source ~/.bashrc
# 最后检查ascp是不是能用了
ascp --help

下载参考文件

http://fungi.ensembl.org/Fusarium_graminearum_gca_000240135/Info/Index
基因组文件(fasta)
解释文件(gtf)
蛋白序列

image.png

genome.fa 由四条染色体单独下载cat合并,防止未拼到染色体上面的序列增加分析难度。
但是发现比对率特别低,所以重新准备整理基因组文件
NCBI下载数据
https://www.ncbi.nlm.nih.gov/genome/?term=Fusarium+graminearum

$ gffread -T -o genome.gtf GCF_000240135.3_ASM24013v3_genomic.gff

下载测序数据,
https://sra-explorer.info/?
流程:FastQ全部选中--add to collection--saved datasets--Aspera commands for downloading FastQ files

SRA下载页面

复制全部代码,存于download.sh nohup sh download.sh &

#!/usr/bin/env bash
ascp -QT -l 300m -P33001 -i $HOME/.aspera/connect/etc/asperaweb_id_dsa.openssh [email protected]:vol1/fastq/SRR856/002/SRR8568982/SRR8568982_1.fastq.gz . && mv SRR8568982_1.fastq.gz SRR8568982_Strand-specific_RNA-Seq_of_wheat_spikelets_after_72h_post-inoculation_of_Fusarium_graminearum_biological_replicate_3__1.fastq.gz
ascp -QT -l 300m -P33001 -i $HOME/.aspera/connect/etc/asperaweb_id_dsa.openssh [email protected]:vol1/fastq/SRR856/002/SRR8568982/SRR8568982_2.fastq.gz . && mv SRR8568982_2.fastq.gz SRR8568982_Strand-specific_RNA-Seq_of_wheat_spikelets_after_72h_post-inoculation_of_Fusarium_graminearum_biological_replicate_3__2.fastq.gz
ascp -QT -l 300m -P33001 -i $HOME/.aspera/connect/etc/asperaweb_id_dsa.openssh [email protected]:vol1/fastq/SRR856/004/SRR8568984/SRR8568984_1.fastq.gz . && mv SRR8568984_1.fastq.gz SRR8568984_Strand-specific_RNA-Seq_of_wheat_spikelets_after_72h_post-inoculation_of_Fusarium_graminearum_biological_replicate_2__1.fastq.gz
ascp -QT -l 300m -P33001 -i $HOME/.aspera/connect/etc/asperaweb_id_dsa.openssh [email protected]:vol1/fastq/SRR856/004/SRR8568984/SRR8568984_2.fastq.gz . && mv SRR8568984_2.fastq.gz SRR8568984_Strand-specific_RNA-Seq_of_wheat_spikelets_after_72h_post-inoculation_of_Fusarium_graminearum_biological_replicate_2__2.fastq.gz
ascp -QT -l 300m -P33001 -i $HOME/.aspera/connect/etc/asperaweb_id_dsa.openssh [email protected]:vol1/fastq/SRR856/002/SRR8569392/SRR8569392_1.fastq.gz . && mv SRR8569392_1.fastq.gz SRR8569392_Strand-specific_RNA-Seq_of_wheat_spikelets_after_0h_post-inoculation_of_Fusarium_graminearum_biological_replicate_3__1.fastq.gz
ascp -QT -l 300m -P33001 -i $HOME/.aspera/connect/etc/asperaweb_id_dsa.openssh [email protected]:vol1/fastq/SRR856/002/SRR8569392/SRR8569392_2.fastq.gz . && mv SRR8569392_2.fastq.gz SRR8569392_Strand-specific_RNA-Seq_of_wheat_spikelets_after_0h_post-inoculation_of_Fusarium_graminearum_biological_replicate_3__2.fastq.gz
ascp -QT -l 300m -P33001 -i $HOME/.aspera/connect/etc/asperaweb_id_dsa.openssh [email protected]:vol1/fastq/SRR856/001/SRR8569391/SRR8569391_1.fastq.gz . && mv SRR8569391_1.fastq.gz SRR8569391_Strand-specific_RNA-Seq_of_wheat_spikelets_after_0h_post-inoculation_of_Fusarium_graminearum_biological_replicate_2__1.fastq.gz
ascp -QT -l 300m -P33001 -i $HOME/.aspera/connect/etc/asperaweb_id_dsa.openssh [email protected]:vol1/fastq/SRR856/001/SRR8569391/SRR8569391_2.fastq.gz . && mv SRR8569391_2.fastq.gz SRR8569391_Strand-specific_RNA-Seq_of_wheat_spikelets_after_0h_post-inoculation_of_Fusarium_graminearum_biological_replicate_2__2.fastq.gz
ascp -QT -l 300m -P33001 -i $HOME/.aspera/connect/etc/asperaweb_id_dsa.openssh [email protected]:vol1/fastq/SRR856/003/SRR8568983/SRR8568983_1.fastq.gz . && mv SRR8568983_1.fastq.gz SRR8568983_Strand-specific_RNA-Seq_of_wheat_spikelets_after_72h_post-inoculation_of_Fusarium_graminearum_biological_replicate_1__1.fastq.gz
ascp -QT -l 300m -P33001 -i $HOME/.aspera/connect/etc/asperaweb_id_dsa.openssh [email protected]:vol1/fastq/SRR856/003/SRR8568983/SRR8568983_2.fastq.gz . && mv SRR8568983_2.fastq.gz SRR8568983_Strand-specific_RNA-Seq_of_wheat_spikelets_after_72h_post-inoculation_of_Fusarium_graminearum_biological_replicate_1__2.fastq.gz
ascp -QT -l 300m -P33001 -i $HOME/.aspera/connect/etc/asperaweb_id_dsa.openssh [email protected]:vol1/fastq/SRR856/003/SRR8569393/SRR8569393_1.fastq.gz . && mv SRR8569393_1.fastq.gz SRR8569393_Strand-specific_RNA-Seq_of_wheat_spikelets_after_24h_post-inoculation_of_Fusarium_graminearum_biological_replicate_1__1.fastq.gz
ascp -QT -l 300m -P33001 -i $HOME/.aspera/connect/etc/asperaweb_id_dsa.openssh [email protected]:vol1/fastq/SRR856/003/SRR8569393/SRR8569393_2.fastq.gz . && mv SRR8569393_2.fastq.gz SRR8569393_Strand-specific_RNA-Seq_of_wheat_spikelets_after_24h_post-inoculation_of_Fusarium_graminearum_biological_replicate_1__2.fastq.gz
ascp -QT -l 300m -P33001 -i $HOME/.aspera/connect/etc/asperaweb_id_dsa.openssh [email protected]:vol1/fastq/SRR856/004/SRR8569394/SRR8569394_1.fastq.gz . && mv SRR8569394_1.fastq.gz SRR8569394_Strand-specific_RNA-Seq_of_wheat_spikelets_after_48h_post-inoculation_of_Fusarium_graminearum_biological_replicate_3__1.fastq.gz
ascp -QT -l 300m -P33001 -i $HOME/.aspera/connect/etc/asperaweb_id_dsa.openssh [email protected]:vol1/fastq/SRR856/004/SRR8569394/SRR8569394_2.fastq.gz . && mv SRR8569394_2.fastq.gz SRR8569394_Strand-specific_RNA-Seq_of_wheat_spikelets_after_48h_post-inoculation_of_Fusarium_graminearum_biological_replicate_3__2.fastq.gz
ascp -QT -l 300m -P33001 -i $HOME/.aspera/connect/etc/asperaweb_id_dsa.openssh [email protected]:vol1/fastq/SRR856/000/SRR8569390/SRR8569390_1.fastq.gz . && mv SRR8569390_1.fastq.gz SRR8569390_Strand-specific_RNA-Seq_of_wheat_spikelets_after_0h_post-inoculation_of_Fusarium_graminearum_biological_replicate_1__1.fastq.gz
ascp -QT -l 300m -P33001 -i $HOME/.aspera/connect/etc/asperaweb_id_dsa.openssh [email protected]:vol1/fastq/SRR856/000/SRR8569390/SRR8569390_2.fastq.gz . && mv SRR8569390_2.fastq.gz SRR8569390_Strand-specific_RNA-Seq_of_wheat_spikelets_after_0h_post-inoculation_of_Fusarium_graminearum_biological_replicate_1__2.fastq.gz
ascp -QT -l 300m -P33001 -i $HOME/.aspera/connect/etc/asperaweb_id_dsa.openssh [email protected]:vol1/fastq/SRR856/007/SRR8569387/SRR8569387_1.fastq.gz . && mv SRR8569387_1.fastq.gz SRR8569387_Strand-specific_RNA-Seq_of_wheat_spikelets_after_24h_post-inoculation_of_Fusarium_graminearum_biological_replicate_3__1.fastq.gz
ascp -QT -l 300m -P33001 -i $HOME/.aspera/connect/etc/asperaweb_id_dsa.openssh [email protected]:vol1/fastq/SRR856/007/SRR8569387/SRR8569387_2.fastq.gz . && mv SRR8569387_2.fastq.gz SRR8569387_Strand-specific_RNA-Seq_of_wheat_spikelets_after_24h_post-inoculation_of_Fusarium_graminearum_biological_replicate_3__2.fastq.gz
ascp -QT -l 300m -P33001 -i $HOME/.aspera/connect/etc/asperaweb_id_dsa.openssh [email protected]:vol1/fastq/SRR856/009/SRR8569389/SRR8569389_1.fastq.gz . && mv SRR8569389_1.fastq.gz SRR8569389_Strand-specific_RNA-Seq_of_wheat_spikelets_after_48h_post-inoculation_of_Fusarium_graminearum_biological_replicate_2__1.fastq.gz
ascp -QT -l 300m -P33001 -i $HOME/.aspera/connect/etc/asperaweb_id_dsa.openssh [email protected]:vol1/fastq/SRR856/009/SRR8569389/SRR8569389_2.fastq.gz . && mv SRR8569389_2.fastq.gz SRR8569389_Strand-specific_RNA-Seq_of_wheat_spikelets_after_48h_post-inoculation_of_Fusarium_graminearum_biological_replicate_2__2.fastq.gz
ascp -QT -l 300m -P33001 -i $HOME/.aspera/connect/etc/asperaweb_id_dsa.openssh [email protected]:vol1/fastq/SRR856/008/SRR8569388/SRR8569388_1.fastq.gz . && mv SRR8569388_1.fastq.gz SRR8569388_Strand-specific_RNA-Seq_of_wheat_spikelets_after_48h_post-inoculation_of_Fusarium_graminearum_biological_replicate_1__1.fastq.gz
ascp -QT -l 300m -P33001 -i $HOME/.aspera/connect/etc/asperaweb_id_dsa.openssh [email protected]:vol1/fastq/SRR856/008/SRR8569388/SRR8569388_2.fastq.gz . && mv SRR8569388_2.fastq.gz SRR8569388_Strand-specific_RNA-Seq_of_wheat_spikelets_after_48h_post-inoculation_of_Fusarium_graminearum_biological_replicate_1__2.fastq.gz
ascp -QT -l 300m -P33001 -i $HOME/.aspera/connect/etc/asperaweb_id_dsa.openssh [email protected]:vol1/fastq/SRR856/006/SRR8569386/SRR8569386_1.fastq.gz . && mv SRR8569386_1.fastq.gz SRR8569386_Strand-specific_RNA-Seq_of_wheat_spikelets_after_24h_post-inoculation_of_Fusarium_graminearum_biological_replicate_2__1.fastq.gz
ascp -QT -l 300m -P33001 -i $HOME/.aspera/connect/etc/asperaweb_id_dsa.openssh [email protected]:vol1/fastq/SRR856/006/SRR8569386/SRR8569386_2.fastq.gz . && mv SRR8569386_2.fastq.gz SRR8569386_Strand-specific_RNA-Seq_of_wheat_spikelets_after_24h_post-inoculation_of_Fusarium_graminearum_biological_replicate_2__2.fastq.gz

存于download.sh nohup sh downlaod.sh & 下载到全部测序数据

使用 fastp 软件对数据进行过滤

由于双端测序,使用官网推荐方法


官网使用方法

脚本


for i in FG_0h_REP1 FG_0h_REP2 FG_0h_REP3 ; do
{
    fastp -i /ifs1/User/haozhigang/rnaseq-fg/1.Rawdate/${i}_1.fastq.gz -o /ifs1/User/haozhigang/rnaseq-fg/clean_data/${i}_1.fastq.gz -I /ifs1/User/haozhigang/rnaseq-fg/1.Rawdate/${i}_2.fastq.gz -O /ifs1/User/haozhigang/rnaseq-fg/clean_data/${i}_2.fastq.gz -h /ifs1/User/haozhigang/rnaseq-fg/clean_data/${i}.html -j /ifs1/User/haozhigang/rnaseq-fg/clean_data/${i}.json
}&
done
wait

第一部分:Mapping

step1.hisat_build.sh
hisat2-build ../ref/genome.fasta ../ref/genome 1>hisat2-build.log 2>&1

建立基因组索引

建立基因组索引
step2.run_hisat.sh
测序数据比对到基因组,生成sam文件

for i in  FG_0h_REP1 FG_0h_REP2 FG_0h_REP3 FG_24h_REP1 FG_24h_REP2 FG_24h_REP3 FG_48h_REP1 FG_48h_REP2 FG_48h_REP3 FG_72h_REP1 FG_72h_REP2 FG_72h_REP3 ; do
         {
             hisat2 --new-summary -p 10 -x /ifs1/User/haozhigang/rnaseq-fg/ref/genome -1 /ifs1/User/haozhigang/rnaseq-fg/clean_data/${i}_1.fastq.gz -2 /ifs1/User/haozhigang/rnaseq-fg/clean_data/${i}_2.fastq.gz -S /ifs1/User/haozhigang/rnaseq-fg/clean_data/${i}.sam 1>/ifs1/User/haozhigang/rnaseq-fg/clean_data/${i}.log 2>&1
                                                                                               }&
               done

step3 sam2bam.sh

for i in FG_0h_REP1  FG_0h_REP2 FG_0h_REP3 FG_24h_REP1 FG_24h_REP2 FG_24h_REP3 FG_48h_REP1 FG_48h_REP2 FG_48h_REP3 FG_72h_REP1 FG_72h_REP2 FG_72h_REP3 ; do
          {
        samtools sort -o /ifs1/User/haozhigang/rnaseq-fg/clean_data/${i}.bam /ifs1/User/haozhigang/rnaseq-fg/clean_data/${i}.sam
                                                                                                                              }&
               done


这一步由于计算资源占用太大,每一个重复一组进行数据整理

step4bamindex.sh

for i in FG_0h_REP1  FG_0h_REP2 FG_0h_REP3 FG_24h_REP1 FG_24h_REP2 FG_24h_REP3 FG_48h_REP1 FG_48h_REP2 FG_48h_REP3 FG_72h_REP1 FG_72h_REP2 FG_72h_REP3; do
          {
       /ifs1/Software/bin/samtools index /ifs1/User/haozhigang/rnaseq-fg/clean_data/${i}.bam
          }&
               done

bam文件建立索引

2020年/4/20
方法没错,但是


image.png

samtools flagstat FG_0h_REP1.bam
查看比对率:重要的质控数据

你可能感兴趣的:(禾谷镰刀菌转录组分析--准备数据)