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

Interpreting a possible haploid #147

Closed
jacquelinemattos opened this issue May 15, 2024 · 5 comments
Closed

Interpreting a possible haploid #147

jacquelinemattos opened this issue May 15, 2024 · 5 comments
Labels
genomescope_included smudgeplot_included if smudgeplot was posted with the quesiton / problem

Comments

@jacquelinemattos
Copy link

Hi,
I'm having trouble understanding my smudgeplot.
We have some reasons to believe that our species' genome should be a haploid, but based on what you replied on another thread regarding the 1n coverage peak, we are no longer sure about the ploidy. Could you please help figure it out?
It's a fungal species that was sequenced using PacBio HiFi.

Our plot looks like this:

pucciniasilphii_ploidy_smudgeplot

And this is the genomescope profile for the same genome:

transformed_linear_plot

If you could please help understanding this ploidy, it would be greatly appreciated.
Thanks so much!

@KamilSJaron
Copy link
Owner

Hi Jacqueline, smudgeplot find all the k-mer pairs that are similar to each other (i.e. different by 1 nt) and generates this plot out of it while guessing the 1n coverage.

Basically what these two plots give you together, that there is no peak ~100x, all the k-mer pairs that are similar to each other have 200x coverage each - suggestion that that's the 1n peak. And given how big it is, it's quite possibly a haploid genome (obviously, there are some exceptions, very heterozygous genomes can occasionally look like this as well, but I don't think it's your case)

What's the species?

@KamilSJaron KamilSJaron added smudgeplot_included if smudgeplot was posted with the quesiton / problem genomescope_included labels May 16, 2024
@jacquelinemattos
Copy link
Author

Hi @KamilSJaron,
Thanks so much for your response! This is really helpful.
Yes, we're working with a Puccinia sp. isolate (a rust fungus). We also thought it would be a haploid, so I'm guessing we could conclude that afterwards?
Thank you!

@KamilSJaron
Copy link
Owner

Oooh, a rust fungus! Cool!

Yeah, i would be comfortable thinking that. I would try to fit the haploid model in genomescope (-p 1) so you have a meaningful genome size estimate. Then keep an eye on the size of the assembled genome. if they would be a lot different you can post here again and we can give it a think again.

@jacquelinemattos
Copy link
Author

Hi @KamilSJaron!
Oh, for sure, I'll do that. But this has already been super helpful!! We really wanted to underestand the ploidy of this species and now it seems to be much more clear!
Thank you so much for you help! I'll come back and post again if needed, for sure.

@KamilSJaron
Copy link
Owner

Great, I will close for now, feel free to reopen

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
genomescope_included smudgeplot_included if smudgeplot was posted with the quesiton / problem
Projects
None yet
Development

No branches or pull requests

2 participants