Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Would it work with Wolbachia (or other symbionts) ? #33

Open
sam217pa opened this issue May 18, 2022 · 6 comments
Open

Would it work with Wolbachia (or other symbionts) ? #33

sam217pa opened this issue May 18, 2022 · 6 comments

Comments

@sam217pa
Copy link
Contributor

Hi Remi,

I was wondering how mitochondrial-genomes specific is MitoFinder ?
Do you think it would work on Wolbachia for instance ? To me it looks like it should, given a set of Wolbachia genes instead of a set of mitochondrial genes.
We are looking for integrated and "free-living" Wolbachia, the parallel with NUMTS / free-living mitochondria made me thought we could use mitofinder for this task.

Thanks !

@RemiAllio
Copy link
Owner

Hi,

Thank you for your message.

Indeed, when I developed MitoFinder I thought about he fact that it could be also used for Wolbachia genomes. That's why I added the option "--new-genes". The specificity of MitoFinder is indeed mainly due to the fact that it is able to read/understand mitochondrial gene names (e.g. COX1=COI) and thus easily create a reference database for further mitochondrial search. The other reason is that mitochondrial genes are (often) monoexonic and thus easier to annotate correctly (I have also added an option for allowing introns but you need a closer reference genome to get accurate results). However, I encourage you to test it by giving a set of annotated Wolbachia genomes as reference and specify the appropriate genetic code and the option "--new-genes" to see if it works. It should at least succeed in catching the contig(s) corresponding to Wolbachia in my opinion. The only constraint is that the annotated gene names of the different reference Wolbachia (if you give more than one) must be identical (MitoFinder is not designed to find synonymy on these genes).

Please let me know if it works and if you are satisfied with the results. I haven't had a chance to try it myself.
Best,
Rémi

@sam217pa
Copy link
Contributor Author

Nice ! I'll let you know, probably in 2-3 weeks 🤞

@RemiAllio
Copy link
Owner

Nice 👍
I look forward to hearing from you!
Cheers,
Rémi

@Andy-B-123
Copy link

Oooh I would be very interested in your results @sam217pa!

And just a note, thank you @RemiAllio for this program, loving the versatility and annotation output of MitoFinder and working with the data this produces has been a genuine joy!

@RemiAllio
Copy link
Owner

Thank you @Andy-B-123 for such a positive feedback!
Glad you liked it 😃

@sam217pa
Copy link
Contributor Author

sam217pa commented Jun 8, 2022

Oops shouldn't have said 3 weeks it's been 21 days already 😬😬😬
we're putting up some controls and a set of "universal" wolbachia core-genome genes to try it out rn,
I'll keep you posted @Andy-B-123

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants