Objective
The objective for this lab was to apply the skills we learned in the previous labs to analyze a data set.
Purpose
The purpose for this lab was to exercise our knowledge about Qiime2 and analyze a set of data in preparation for the analysis of our own.
Procedure/Data
As always, we bagan by installing Miniconda and Qiime2 on our computers. Then, we downloaded the folder containing all the data we would need. After changing the directory to that folder and activating Qiime2, we proceeded to import the sequences as a Qiime2 artifact using the code
qiime tools import \
–type EMPPairedEndSequences \
–input-path emp-paired-end-sequences \
–output-path emp-paired-end-sequences.qza
Once this was accomplished, we demultiplexed the sequences using the follwing codes.
qiime demux emp-paired \
–m-barcodes-file sample-metadata.tsv \
–m-barcodes-column BarcodeSequence \
–i-seqs emp-paired-end-sequences.qza \
–o-per-sample-sequences demux.qza \
qiime demux summarize \
–i-data demux.qza \
–o-visualization demux.qzv
We then used DADA2 to denoise the data and ensure our results are reliable.
qiime dada2 denoise-paired \
–i-demultiplexed-seqs demux.qza \
–p-trunc-len-f 220 \
–p-trunc-len-r 220 \
–o-table table.qza \
–o-representative-sequences rep-seqs.qza \
–o-denoising-stats denoising-stats.qza
After this, we created a feature table of our data and converted it into a .qzv file so it could be viewed.
qiime feature-table summarize \
–i-table table.qza \
–o-visualization table.qzv \
–m-sample-metadata-file sample-metadata.tsv
qiime feature-table tabulate-seqs \
–i-data rep-seqs.qza \
–o-visualization rep-seqs.qzv
qiime metadata tabulate \
–m-input-file denoising-stats.qza \
–o-visualization denoising-stats.qzv
Then, we created a phylogenetic tree that connected our DNA sequences to their respective classifications.
qiime phylogeny align-to-tree-mafft-fasttree \
–i-sequences rep-seqs.qza \
–o-alignment aligned-rep-seqs.qza \
–o-masked-alignment masked-aligned-rep-seqs.qza \
–o-tree unrooted-tree.qza \
–o-rooted-tree rooted-tree.qza
Finally, we created a taxa bar plot to display all this data in an interactive way that could better show all the specifics of the data.
qiime feature-classifier classify-sklearn –i-classifier silva-132-99-515-806-nb-classifier.qza –i-reads rep-seqs.qza –o-classification taxonomy.qza
qiime metadata tabulate \
–m-input-file taxonomy.qza \
–o-visualization taxonomy.qzv
qiime taxa barplot \
–i-table table.qza \
–i-taxonomy taxonomy.qza \
–m-metadata-file sample-metadata.tsv \
–o-visualization taxa-bar-plots.qzv
Conclusion
Overall, this lab was much like last lab where we practiced using the Moving Pictures tutorial on the Qiime2 website. Practicing with Qiime2 is much needed, and I feel as though I am getting more comfortable using it. Additionally, my understanding of each step and why it is important is growing. I am excited to perform these anaylses on our own data in the future.