Skip to content

treangenlab/wasterwater_arg_metagenomics

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

5 Commits
 
 
 
 
 
 

Repository files navigation

wasterwater_arg_metagenomics

antibiotic resistance genes (ARGs) discovery and comparison among metagenomics sequenced data pipeline for paper "Profiling the fate of antibiotic resistance across a WWTP using metagenomics and epicPCR for risk estimation"

Part1 Sample ARG Identification

Required libraries, databases and tools:

megahit https://github.com/voutcn/megahit 
rgi https://github.com/arpcard/rgi
card-databse https://card.mcmaster.ca/download 
mob-suite https://github.com/phac-nml/mob-suite
bowtie2 https://github.com/BenLangmead/bowtie2
blastn and blastx 
MGE-database https://bench.cs.vt.edu/ftp/data/databases/MGEs90.fasta 
centrifuge https://github.com/DaehwanKimLab/centrifuge 
PlasFlow https://github.com/smaegol/PlasFlow
seqkit https://github.com/shenwei356/seqkit
seqtk https://github.com/lh3/seqtk 
os, subprocess, collections, pathlib, numpy, pandas, BioPython

Usage

  1. Install all the tools, libraries and their related databases.
  2. Replace the second block in the jupyter notebook with custom path.
  3. Run the code, but need to switch environment for PlasFlow.

Part2 EpicPCR Analysis

Methods

Step1 Trim and filter out low qulity reads

a) Trim reads using porechop

porechop -t 15 -i [temp]_input.fastq --discard_middle -o trimmed_[temp].fastq

b) Remove reads with quality score below 7

NanoFilt -q 7 trimmed_[temp].fastq > nanofilt_[temp].fastq

Step2 Blast agaisnt linker primer sequences

blastn -db reference_[temp]_bridge -query [temp]_input_seq.fasta -outfmt 6 -num_threads 10 -out arg_mapping/[temp]_bridge_blast

Step3 Split fusion reads using a customized script

python split_read_blast.py

Step4 Taxonomic Classification

emu abundance [temp]_input.fasta --output-dir output_dir --output-basename [temp]_emu --threads 10 --db dbs

About

No description, website, or topics provided.

Resources

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published