![]() |
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
De-novo assembly of bacteria genomes - which tools? | maasha | Bioinformatics | 8 | 08-12-2009 02:41 PM |
DGE on procaryotic organisms | huguesparri | Illumina/Solexa | 1 | 05-26-2009 06:47 AM |
Comparing SNP/indel detection tools on 454 | dePhi | 454 Pyrosequencing | 2 | 03-07-2009 05:41 AM |
comparing results from two different reference genomes | BAJ | Bioinformatics | 2 | 02-24-2009 08:38 AM |
Genomic enrichment on non-model organisms... | MattB | General | 2 | 08-08-2008 02:07 AM |
![]() |
|
Thread Tools |
![]() |
#1 |
Junior Member
Location: Australia Join Date: May 2012
Posts: 1
|
![]()
Hi Bioinformatics gurus,
Could you suggest me any tool(s) for comparing genomes of two (or more) different organisms to have at the end a list of genes of organism A which are lacking in organism B? And also, a list of common genes in both organisms which are drastically different. Thanks |
![]() |
![]() |
![]() |
#2 | |
Senior Member
Location: Cambridge, UK Join Date: May 2010
Posts: 311
|
![]()
Maybe InParanoid http://www.ncbi.nlm.nih.gov/pubmed/19892828 can help you? From the documentation http://inparanoid.sbc.su.se/download...t/relnotes.txt
Quote:
|
|
![]() |
![]() |
![]() |
#3 |
Member
Location: UK Join Date: Jan 2009
Posts: 15
|
![]()
Hi,
I'm assuming they've been annotated and are bacteria/virus. You can use MAUVE, which can do the alignment of 2 or more genomes and there is an option output orthologues based on your cutoffs. The output is a table (which can be opened in excel), where across the top (X axis) are the organism name, below which are the genes present. Comparing the columns, you can pull out whats shared and whats not. Alternatively, you could you Blast to compare and pull out similar sequences between the genomes, again using what ever cutoff you like - that can be the more labour intensive way and gets complicated with 3 or more genomes. Cheers, Raj |
![]() |
![]() |
![]() |
#4 |
Member
Location: Raleigh, NC Join Date: Nov 2008
Posts: 51
|
![]()
A clutering tool like cd-hit would work.
|
![]() |
![]() |
![]() |
#5 |
Member
Location: RTP, NC, USA Join Date: Jun 2009
Posts: 28
|
![]()
No second dariober's suggestion to use InParanoid. Its summary information and definitions of orthologs vs paralogs is quite useful. The key is that you must have annotated organisms to start with. And keep in mind the most likely differences are going to arise from genome assembly errors. For example my first pass filter on gene loss is to run peptide BLAT/BLAST and look for hits on the *_random chromosomes. I'd often find fragments of the gene at the edges of genome gaps. Sometimes half the gene would be on chr6, the other half on chr6_random, but you could see pretty much anything, depending upon the quality of the assembly and the annotation. But InParanoid should in theory help you figure out what to follow up with a deeper analysis, ie what holes there are to fill.
|
![]() |
![]() |
![]() |
#6 |
Member
Location: Montreal Join Date: Jan 2014
Posts: 21
|
![]()
Hi, I'm facing the same kind of situation. Anybody here used InParanoid successfully?
|
![]() |
![]() |
![]() |
#8 |
Member
Location: Montreal Join Date: Jan 2014
Posts: 21
|
![]()
Thanks GenoMax, I'm now trying Mauve. As for InParanoid, I'm not able to make it work on my transcriptomes. It's working with a small subset of sequences (a few thousands sequences) but it freezes forever when I try with the whole transcriptomes (hundreds of thousands sequences).
|
![]() |
![]() |
![]() |
Tags |
genes diversity |
Thread Tools | |
|
|