#Attempt 4 at 16S rRNA metagonomic samples #Erinne Stirling - 07/09/2018 #################### START #################### #In terminal source activate qiime2-2018.6 source tab-qiime #Change to working directory cd ~/Documents/BioInf/Divstats/ITS_workspace/ #################### FILES #################### #Validate metadata file used Keemei: cloud-based validation of tabular bioinformatics file formats in Google Sheets. Rideout JR, Chase JH, Bolyen E, Ackermann G, González A, Knight R, Caporaso JG. GigaScience. 2016;5:27. http://dx.doi.org/10.1186/s13742-016-0133-6 #Metadata file: ~/divstats-meta.tsv #Manifest file: /Documents/BioInf/Divstats/files/AGRF_CAGRF17715_BWV5G/manifest.txt #################### IMPORT #################### cd ~/Documents/BioInf/Divstats/files/AGRF_CAGRF17715_BWV5G qiime tools import \ --type 'SampleData[PairedEndSequencesWithQuality]' \ --input-path manifest.txt \ --source-format PairedEndFastqManifestPhred33 \ --output-path ~/Documents/BioInf/Divstats/ITS_workspace/demux.qza cd ~/Documents/BioInf/Divstats/ITS_workspace/ #Visualise the data qiime demux summarize \ --i-data demux.qza \ --o-visualization demux.qzv qiime tools view demux.qzv #Add metadata qiime metadata tabulate \ --m-input-file divstats-meta.tsv \ --o-visualization tabulated_meta.qzv qiime tools view tabulated_meta.qzv #################### DADA #################### #Denoise qiime dada2 denoise-paired \ --i-demultiplexed-seqs demux.qza \ --p-trim-left-f 0 \ --p-trunc-len-f 250 \ --p-trim-left-r 0 \ --p-trunc-len-r 250 \ --o-representative-sequences rep-seqs.qza \ --o-table table.qza \ --o-denoising-stats stats.qza qiime metadata tabulate \ --m-input-file stats.qza \ --o-visualization stats.qzv qiime tools view stats.qzv #################### ANALYSIS PREP #################### #Sequence number per sample qiime feature-table summarize \ --i-table table.qza \ --o-visualization table.qzv \ --m-sample-metadata-file divstats-meta.tsv qiime tools view table.qzv #------------------------------------------------------Sample depth = 30046 qiime feature-table tabulate-seqs \ --i-data rep-seqs.qza \ --o-visualization rep-seqs.qzv qiime tools view rep-seqs.qzv #################### TREES #################### #Alignment qiime alignment mafft \ --i-sequences rep-seqs.qza \ --o-alignment aligned-rep-seqs.qza #Mask the highly variable bits qiime alignment mask \ --i-alignment aligned-rep-seqs.qza \ --o-masked-alignment masked-aligned-rep-seqs.qza #Fasttree qiime phylogeny fasttree \ --i-alignment masked-aligned-rep-seqs.qza \ --o-tree unrooted-tree.qza #Root it qiime phylogeny midpoint-root \ --i-tree unrooted-tree.qza \ --o-rooted-tree rooted-tree.qza #################### DIVERSITY - ALPHA #################### qiime diversity core-metrics-phylogenetic \ --i-phylogeny rooted-tree.qza \ --i-table table.qza \ --p-sampling-depth 5900 \ --m-metadata-file divstats-meta.tsv \ --output-dir core-metrics-results #***FAITH***# qiime diversity alpha-group-significance \ --i-alpha-diversity core-metrics-results/faith_pd_vector.qza \ --m-metadata-file divstats-meta.tsv \ --o-visualization core-metrics-results/faith-pd-group-significance.qzv qiime tools view core-metrics-results/faith-pd-group-significance.qzv #------------------------------------------------------Kruskall-Wallace #------------------------------------------------------H 17.9003 #------------------------------------------------------P 0.0013 #------------------------------------------------------PAIRWISE #------------------------------------------------------P<0.05 40.C-> control #------------------------------------------------------ 260 #------------------------------------------------------ 320 #------------------------------------------------------ 150.C-> control #------------------------------------------------------ 260 #------------------------------------------------------ 320 #------------------------------------------------------ 260.C-> 40 #------------------------------------------------------ 150 #------------------------------------------------------ 320.C-> 40 #------------------------------------------------------ 150 qiime diversity alpha-group-significance \ --i-alpha-diversity core-metrics-results/evenness_vector.qza \ --m-metadata-file divstats-meta.tsv \ --o-visualization core-metrics-results/evenness-group-significance.qzv qiime tools view core-metrics-results/evenness-group-significance.qzv #------------------------------------------------------Kruskall-Wallace #------------------------------------------------------H 10.8775 #------------------------------------------------------P 0.0280 #------------------------------------------------------PAIRWISE #------------------------------------------------------P<0.05 40.C-> 150 #------------------------------------------------------ 150.C-> 260 #------------------------------------------------------ 320 #------------------------------------------------------ 260.C-> ns #------------------------------------------------------ 320.C-> ns #------------------------------------------------------ qiime diversity alpha-correlation \ --i-alpha-diversity core-metrics-results/faith_pd_vector.qza \ --m-metadata-file divstats-meta.tsv \ --o-visualization core-metrics-results/faith-alpha-correlation.qzv qiime tools view core-metrics-results/faith-alpha-correlation.qzv #------------------------------------------------------Respiration #------------------------------------------------------Spearman #------------------------------------------------------Test-stat -0.7838 #------------------------------------------------------P 0.0000 #***SHANNON***# qiime diversity alpha \ --i-table table.qza \ --p-metric shannon \ --o-alpha-diversity shannon.qza qiime diversity alpha-group-significance \ --i-alpha-diversity shannon.qza \ --m-metadata-file divstats-meta.tsv \ --o-visualization shannon.qzv qiime tools view Shannon.qzv #################### DIVERSITY - BETA #################### qiime diversity beta-group-significance \ --i-distance-matrix core-metrics-results/unweighted_unifrac_distance_matrix.qza \ --m-metadata-file divstats-meta.tsv \ --m-metadata-column Temp \ --o-visualization core-metrics-results/unweighted-unifrac-temp-significance.qzv \ --p-pairwise qiime tools view core-metrics-results/unweighted-unifrac-temp-significance.qzv #------------------------------------------------------PERMANOVA #------------------------------------------------------pseudo-F #------------------------------------------------------Test-stat 3.10431 #------------------------------------------------------P 0.001 #------------------------------------------------------PAIRWISE #------------------------------------------------------P<0.05 40.C-> control #------------------------------------------------------ 260 #------------------------------------------------------ 320 #------------------------------------------------------ 150.C-> control #------------------------------------------------------ 260 #------------------------------------------------------ 320 #------------------------------------------------------ 260.C-> ns #------------------------------------------------------ 320.C-> control #------------------------------------------------------ ns qiime diversity beta-group-significance \ --i-distance-matrix core-metrics-results/weighted_unifrac_distance_matrix.qza \ --m-metadata-file divstats-meta.tsv \ --m-metadata-column Temp \ --o-visualization core-metrics-results/weighted-unifrac-temp-significance.qzv \ --p-pairwise qiime tools view core-metrics-results/weighted-unifrac-temp-significance.qzv #################### TRAIN A CLASSIFIER #################### #Couldn't download the files; used a premed classifier instead. mv unite-ver7-99-classifier-01.12.2017.qza classifier.qza #################### TAXONOMY #################### qiime feature-classifier classify-sklearn \ --i-classifier classifier.qza \ --i-reads rep-seqs.qza \ --o-classification taxonomy.qza qiime metadata tabulate \ --m-input-file taxonomy.qza \ --o-visualization taxonomy.qzv qiime tools view taxonomy.qzv qiime taxa barplot \ --i-table table.qza \ --i-taxonomy taxonomy.qza \ --m-metadata-file divstats-meta.tsv \ --o-visualization taxa-bar-plots.qzv qiime tools view taxa-bar-plots.qzv #################### q2-composition #################### qiime composition add-pseudocount \ --i-table table.qza \ --o-composition-table comp-table.qza qiime composition ancom \ --i-table comp-table.qza \ --m-metadata-file divstats-meta.tsv \ --m-metadata-column Temp \ --o-visualization ancom-Temp.qzv ### qiime taxa collapse \ --i-table table.qza \ --i-taxonomy taxonomy.qza \ --p-level 4 \ --o-collapsed-table table-l4.qza qiime composition add-pseudocount \ --i-table table-l4.qza \ --o-composition-table comp-table-l4.qza qiime composition ancom \ --i-table comp-table-l4.qza \ --m-metadata-file divstats-meta.tsv \ --m-metadata-column Temp \ --o-visualization ancom-Temp-l4.qzv