Hi,
I have an OTU biom file (obtained from Closed reference QIIME 1.8.0 v) contains 65 samples, I am trying to do analysis for PAN/CORE genome.
I have filtered out the taxonomy from the abundance file (with particular threshold, lets say 60 %), now I have an taxonomy column only in file from all the 65 samples (with threshold 60%), Is there is a way where I can do the functional annotation for it?
Any server/ software is there which can do that? or which do pan (complete) /core (shared) analysis
Any suggestions ?
Best!
Shashank
•
written
by
bioinformaticssrm2011
2 answers
Dear you can also try DAVID functional annotation database. For more detail plz see this link https://david.ncifcrf.gov/tools.jsp
Hopefully this may help you.
If you used GreenGenes 13_5 as a reference database, you can associate the OTUs with protein content of nearest sequenced reference genomes with PICRUST. However, in my opinion this approach is pretty much worthless. The 16S sequence of your OTU representative being relatively similar (60% is not even remotely similar, threshold should be like 99.9% for this stuff) to the 16S sequence of some reference genome does not mean that that the protein contents of these two genomes are even remotely similar..
Log in to answer this question.
-
how to copy and rename files in nextflow
written by fabiomarcelomedvet •
-
Process to merge biom files
written by dpc
-
BPGA Functional Analysis not working
written by rncrake •
-
Visualyzing an otu table
written by life99945 •
-
16S rRNA de novo otu picking analysis
written by dcron13 •
-
UniFrac in R with just the OTU table?
written by David
-
Convert biom to tabular file and parse taxonomic levels
written by naphtali.paul •
-
16S sequencing analysis: from fasta to otu to PCA&Heatmap
written by shlomsauce •
-
How to normalize OTU abundances in 16S-rDNA sequencing data by individual sample amount?
written by yuki.weber81 •
-
PCOA plot for special genus
written by agata88
What you mean with PAN and CORE genome analysis is that you want to find the complete genes/proteins and shared genes/proteins among your OTUs, right? I'm not familiar with biom file, but does it contain the genome sequence(s) of the organisms that you're analyzing? You will need the gene sequences of the whole genomes (or protein sequences of the whole proteomes) to get the pan-genome and core-genome.
Yes, Complete gene and Shared gene.
Biom file don't have the genomic sequence. It looks like-
Where number represents the OTU ID, followed by taxonomy. OTU ID represents the particular sequence associated with the particular taxonomy.
If I incorporated the gene sequence by using the OTU ID corresponding to the taxonomy, Now I have a gene sequence file, than how can I use it for further analysis ?
Cheers!
"OTU ID represents the particular sequence associated with the particular taxonomy." What particular sequence is it? From one gene only? You can't do pan- and core-genome analysis using only one gene from each species/OTU. You need the genome (or better, the proteome) from each OTU, find orthologous gene/protein among the OTUs (I used OrthoMCL for my bacteria), and there you have the core-proteome. The pan-proteome would be the core plus any other proteins of each OTU.