Illumina Metagenomics data

[ssharma]

Hi All,
I am new member to this forum.
Earlier i used to work with 454 data, now i am switching to illumina.
I am getting around 300 million reads (100bp) and its a metagenomic sample. So i am really confused about how to start my analysis.
Earlier i used approaches like blastx but now i think this is not a good option.
So i was just wondering if anyone had done something like this or have some idea on this.

I would really appreciate your help.

Thanks
SS

[rwenang]

Is it a 16S or WGS sample?

[ssharma]

Its a metagenomics environmental sample (no 16s).

[themerlin]

The real question is....what's the question? Are you looking for specific genes, or want to take an inventory of all genes?

Have you tried assembling the reads yet? That's always a little sketchy with mixed communities, but it might be a good place to start.

[ssharma]

thanks for your input themerlin,
Actually mainly its going to be a community study (in nut shell i need to annotate all of the sequences)
Yes i tried assembly but it doesn't look good, but yes i will try again with different programs.

[kmewis]

Is this just sequence data from DNA straight from the environment, or did you clone it into vectors first?

I handle metagenomics data, we do it in fosmids though, so it's easy to assemble contigs from one fosmid (phred/phrap). Trying to do the whole environment at once will likely be tougher. Once I have contigs, we use blastx to looks for homology and tools like fgenesb to find ORFs.

[Dilipmohana]

hi i am new to this site can anyone tell me about effective working in schrodinger plz pass useful video tutorials if possible,

[greigite]

Take a look at MG-RAST for annotation of your data http://metagenomics.nmpdr.org.

Hi All,
I am new member to this forum.
Earlier i used to work with 454 data, now i am switching to illumina.
I am getting around 300 million reads (100bp) and its a metagenomic sample. So i am really confused about how to start my analysis.
Earlier i used approaches like blastx but now i think this is not a good option.
So i was just wondering if anyone had done something like this or have some idea on this.

I would really appreciate your help.

Thanks
SS

[Eric]

Hi,

What do you mean by "annotate" ? Are you looking at "who is there" or "what are the functions" ?
Do you have reference genomes at hand, or genomes of organisms close to the ones in your sample ? Do you have an idea of the complexity of the population ? Is it eukaryote or microbes, or both ?
You can consider first trying to have an idea of the composition of your population, looking at some marker genes (eg : trying to find 16S or 18S reads in your dataset by mapping against reference databases)
If you have known reference genomes, you can also map reads against them, to evaluate the complexity/diversity
For a first glimpse at functions, you can try UniRef50 or KEGG genes (or any other functionally classified reference protein set) as a proxy.

[gridbird]

You can try WebMGA: http://weizhong-lab.ucsd.edu/metagenomic-analysis/

[colindaven]

You can try approaches like
-de novo assembly (metaVelvet, Abyss etc)
-fast clustering - (CD-Hit, RAMMCAP)
-reference based alignment (Genometa)

[faozhi]

I would trim the reads (based on qual and remove adapters), then start assembling.
If you would like to know who are there, you could use MG-RAST or just blastn your trimmed reads against greengenes or SILVA 16S databases.

[seb567]

Hi All,
I am new member to this forum.
Earlier i used to work with 454 data, now i am switching to illumina.
I am getting around 300 million reads (100bp) and its a metagenomic sample. So i am really confused about how to start my analysis.
Earlier i used approaches like blastx but now i think this is not a good option.
So i was just wondering if anyone had done something like this or have some idea on this.

I would really appreciate your help.

Thanks
SS

You can reduce the volume of data by doing a de novo assembly.

mpiexec -n 64 Ray \
 -k \
 31 \
 -p \
 Sample/ERR011142_1.fastq.gz \
 Sample/ERR011142_2.fastq.gz \
 -p \
 Sample/ERR011143_1.fastq.gz \
 Sample/ERR011143_2.fastq.gz \
 -o \
 Assembly

Sébastien Boisvert