Nanopore测序在宏基因组学中的应用

2022-10-27 本文已影响0人生信斑马鱼

视频推荐

请先在B站观看视频，毕竟看视频比看书学习快多了

1.B站检索（第一天：illumina、PacBio、Nanopore测序原理）

视频推荐.png

2.B站检索（Nanopore平台宏基因组测序的优点）

视频推荐2.png

标准流程

本流程复现自Nature communication的文章“Short- and long-read metagenomics expand individualized structural variations in gut microbiomes”

数据来源 https://www.ncbi.nlm.nih.gov/bioproject/PRJNA820119
脚本来源 https://github.com/chen318liang/Gut-Metagenome-Pipeline-Based-on-Nanopore-Sequencing

数据下载

按图示下载数据

downloading1.png

downloading2.png

downloading3.png

downloading4.png

downloading5.png

也可采用sra-toolkit批量下载方法，具体见文
https://www.jianshu.com/p/6f8f9a6d1a56

SRA数据转fastq

使用fasterq-dump转换数据

fasterq-dump  SRR5318040 -t . -e 24;
###-e 24 是24线程

脚本安装

本分析涉及到的分析软件很多，不同的软件安装方法差异很大，以下是安装所需软件的下载源或安装源。
提前熟悉linux的基本命令，awk, sed，grep等等

git 安装

OPERA-MS v0.9.0 mumandco mcaller guppy Prokka CD-hit iqtree2

conda 安装

IGV seqkit dRep gtdbtk samtools MAFFT CAT blast Nanopolish Canu Flye minimap2 MetaWRAP prokka salmon

pip install drep clustalw Emboss ViennaRNA

网页安装

rnammer
https://services.healthtech.dtu.dk/software.php

网页或UI软件可视化

ont-tombo
nanodisco
crsprdetect
prophagehunter
itol

单脚本

CRISPRDetect (https://github.com/davidchyou/CRISPRDetect_2.4)
crisprdetectparser.py (https://github.com/hwalinga/crisprdetect-parser)
~/miniconda3/envs/metawrap/bin/metawrap-scripts/split_salmon_out_into_bins.py

分析流程

测序质量评估

注意设置好对应变量，注意提前测试所需软件及脚本是否可用

#!/bin/bash

threads="20"
assembly="/data/chenliang/Zymo_Community_Standards_data_assembly/spades_hybird_assembly/Zymo-GridION-EVEN_spades_hybird_assembly/scaffolds.fasta"
out="quant_bins"
reads_1="/data/chenliang/MC.Hiseq/MC_1.fastq"
reads_2="/data/chenliang/MC.Hiseq/MC_2.fastq"
tmp=${reads_1##*/}
sample=${tmp%_*}
bin_folder="bin_refinement/metawrap_70_10_bins/"


salmon index -p $threads -t $assembly -i ${out}/assembly_index
salmon quant -i ${out}/assembly_index --libType IU -1 $reads_1 -2 $reads_2 -o ${out}/alignment_files/${sample}.quant --meta -p $threads

home=$(pwd)
cd ${out}/alignment_files/
/software_users/chenliang/miniconda3/envs/metawrap/bin/metawrap-scripts/summarize_salmon_files.py
cd $home
mkdir ${out}/quant_files
#for f in $(ls ${out}/alignment_files/ | grep .quant.counts); do mv ${out}/alignment_files/$f ${out}/quant_files/; done

#n=$(ls ${out}/quant_files/ | grep counts | wc -l)
/software_users/chenliang/miniconda3/envs/metawrap/bin/metawrap-scripts/split_salmon_out_into_bins.py ${out}/quant_files/ $bin_folder $assembly > ${out}/bin_abundance_table.tab

基因组组装

使用的5种策略,Canu, Flye, OPERA-MS, Spades, MetaSPAdes.

使用fastANI计算ANI时需要用到参考基因组，请提前下载
https://zenodo.org/record/3935737#.Y1rh37ZBxD8

详细脚本见
https://github.com/chen318liang/Gut-Metagenome-Pipeline-Based-on-Nanopore-Sequencing/blob/main/scripts/0_Zymo_community_standard_data_assembly.sh

突变分析

这一步没啥特别的，注意提前测试软件和脚本就行。观察脚本，提前将文件置于对应位置的话能省不少事情。
minimap2 samtools make_IGV_snapshots.py emapper.py prokka

详细脚本见
https://github.com/chen318liang/Gut-Metagenome-Pipeline-Based-on-Nanopore-Sequencing/blob/main/scripts/2_Gut_metagenome_SV_analysis.sh

预测原噬菌体、去除冗余、基因注释和构建进化树

这步用到了CD-hit需要编译安装，而iqtree2可在github下载，用于化系统发育树。

#####################
Part III. Prediction of prophage, redundancy removal, gene annotation, and construction of evolutionary trees
###################
# prophage dientification using the machine-learning-based tool ProphageHunter
# Webserver address
https://pro-hunter.genomics.cn/

# Candidate propages clustering by CDhit
CDhit-est -i Candidate_propahge.fasta -c 0.95 -o Candidate_propahge_clusted.fasta

# Viral genomic CDS prediction using multiphate2 (https://github.com/carolzhou/multiPhATE2)
#Set up the configuration file as required and then run cmd as follow
python3 multiPhate.py multiPhate.config


# Create phylogenetic tree using iqtree2 (https://github.com/iqtree/iqtree2)
#Join_MCP_TLS.faa is joined protein sequence with major capsid protein and terminal large subunit 
iqtree2 -s Join_MCP_TLS.faa -m MFP -B 1000 --bnni -T 40

从所有样品的contigs预测CRISPR区

CRISPRDetect.pl安装比较繁琐，现在github下载安装，之后要安装依赖环境，根据对应报错提示安装就行，我花了不少时间。crisprdetectparser.py直接在github下载就行。

#!/usr/bin/bash
## CRISPRDetect (https://github.com/davidchyou/CRISPRDetect_2.4)
## crisprdetectparser.py (https://github.com/hwalinga/crisprdetect-parser)

path="/data/nano_meta_ref/result/spades/hybird_temp/spades"

perl CRISPRDetect.pl -f $path/contigs.fasta -o CRISPRDetect_result \
-check_direction 0 -array_quality_score_cutoff 3 -T 20

python crisprdetectparser.py --spacers-directory spacer_dir --spacers-extension fna CRISPRDetect_result > metadata.tsv