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

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 era-fasp@fasp.sra.ebi.ac.uk: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 era-fasp@fasp.sra.ebi.ac.uk: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 era-fasp@fasp.sra.ebi.ac.uk: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 era-fasp@fasp.sra.ebi.ac.uk: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 era-fasp@fasp.sra.ebi.ac.uk: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 era-fasp@fasp.sra.ebi.ac.uk: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 era-fasp@fasp.sra.ebi.ac.uk: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 era-fasp@fasp.sra.ebi.ac.uk: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 era-fasp@fasp.sra.ebi.ac.uk: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 era-fasp@fasp.sra.ebi.ac.uk: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 era-fasp@fasp.sra.ebi.ac.uk: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 era-fasp@fasp.sra.ebi.ac.uk: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 era-fasp@fasp.sra.ebi.ac.uk: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 era-fasp@fasp.sra.ebi.ac.uk: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 era-fasp@fasp.sra.ebi.ac.uk: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 era-fasp@fasp.sra.ebi.ac.uk: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 era-fasp@fasp.sra.ebi.ac.uk: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 era-fasp@fasp.sra.ebi.ac.uk: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 era-fasp@fasp.sra.ebi.ac.uk: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 era-fasp@fasp.sra.ebi.ac.uk: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 era-fasp@fasp.sra.ebi.ac.uk: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 era-fasp@fasp.sra.ebi.ac.uk: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 era-fasp@fasp.sra.ebi.ac.uk: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 era-fasp@fasp.sra.ebi.ac.uk: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
查看比对率：重要的质控数据

最后编辑于：2021.05.19 23:50:29

人面猴
序言：七十年代末，一起剥皮案震惊了整个滨河市，随后出现的几起案子，更是在滨河造成了极大的恐慌，老刑警刘岩，带你破解...
沈念sama阅读 216,125评论 6赞 498
死咒
序言：滨河连续发生了三起死亡事件，死亡现场离奇诡异，居然都是意外死亡，警方通过查阅死者的电脑和手机，发现死者居然都...
沈念sama阅读 92,293评论 3赞 392
救了他两次的神仙让他今天三更去死
文/潘晓璐我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
开封第一讲书人阅读 162,054评论 0赞 351
道士缉凶录：失踪的卖姜人
文/不坏的土叔我叫张陵，是天一观的道长。经常有香客问我，道长，这世上最难降的妖魔是什么？我笑而不...
开封第一讲书人阅读 58,077评论 1赞 291
港岛之恋（遗憾婚礼）
正文为了忘掉前任，我火速办了婚礼，结果婚礼上，老公的妹妹穿的比我还像新娘。我一直安慰自己，他们只是感情好，可当我...
茶点故事阅读 67,096评论 6赞 388
恶毒庶女顶嫁案：这布局不是一般人想出来的
文/花漫我一把揭开白布。她就那样静静地躺着，像睡着了一般。火红的嫁衣衬着肌肤如雪。梳的纹丝不乱的头发上，一...
开封第一讲书人阅读 51,062评论 1赞 295
城市分裂传说
那天，我揣着相机与录音，去河边找鬼。笑死，一个胖子当着我的面吹牛，可吹牛的内容都是我干的。我是一名探鬼主播，决...
沈念sama阅读 39,988评论 3赞 417
双鸳鸯连环套：你想象不到人心有多黑
文/苍兰香墨我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
开封第一讲书人阅读 38,817评论 0赞 273
万荣杀人案实录
序言：老挝万荣一对情侣失踪，失踪者是张志新（化名）和其女友刘颖，没想到半个月后，有当地人在树林里发现了一具尸体，经...
沈念sama阅读 45,266评论 1赞 310
护林员之死
正文独居荒郊野岭守林人离奇死亡，尸身上长有42处带血的脓包…… 初始之章·张勋以下内容为张勋视角年9月15日...
茶点故事阅读 37,486评论 2赞 331
白月光启示录
正文我和宋清朗相恋三年，在试婚纱的时候发现自己被绿了。大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
茶点故事阅读 39,646评论 1赞 347
活死人
序言：一个原本活蹦乱跳的男人离奇死亡，死状恐怖，灵堂内的尸体忽然破棺而出，到底是诈尸还是另有隐情，我是刑警宁泽，带...
沈念sama阅读 35,375评论 5赞 342
日本核电站爆炸内幕
正文年R本政府宣布，位于F岛的核电站，受9级特大地震影响，放射性物质发生泄漏。R本人自食恶果不足惜，却给世界环境...
茶点故事阅读 40,974评论 3赞 325
男人毒药：我在死后第九天来索命
文/蒙蒙一、第九天我趴在偏房一处隐蔽的房顶上张望。院中可真热闹，春花似锦、人声如沸。这庄子的主人今日做“春日...
开封第一讲书人阅读 31,621评论 0赞 21
一桩弑父案，背后竟有这般阴谋
文/苍兰香墨我抬头看了看天上的太阳。三九已至，却和暖如春，着一层夹袄步出监牢的瞬间，已是汗流浃背。一阵脚步声响...
开封第一讲书人阅读 32,796评论 1赞 268
情欲美人皮
我被黑心中介骗来泰国打工，没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留，地道东北人。一个月前我还...
沈念sama阅读 47,642评论 2赞 368
代替公主和亲
正文我出身青楼，却偏偏与公主长得像，于是被迫代替她去往敌国和亲。传闻我的和亲对象是个残疾皇子，可洞房花烛夜当晚...
茶点故事阅读 44,538评论 2赞 352

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

下载参考文件

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

第一部分：Mapping

推荐阅读更多精彩内容