search

jtlovell / GENESPACE

Other
191 stars 27 forks source link

Synteny / homoeologs in polyploids #144

Closed tallnuttrbgv closed 5 months ago

tallnuttrbgv commented 9 months ago

Hi,

unfortunately I could not work out how to do this from the docs and after several attempts need to ask for help.

I have a genome that is a putative triploid. BUSCOs on the purged assembly are 55% duplicated and the species is known to form triploids. I wanted to use gs to identify homoeolog contigs (the assembly is hifiasm and mostly complete chromosomes). When I ran gs on the assembly to itself the riparian plot only showed single syntenitc groups. I set onlyOgAnchors=F, and nSecondHits=3. What is your recommended method to do this?

Thanks.

image

jtlovell commented 7 months ago

Sorry about the slow response ... were you able to resolve this? We go into detail about how to handle genomes w/ WGDs in the paper. In short, use an outgroup that predates the WGD.

tallnuttrbgv commented 7 months ago

Hi,

I don't think that will achieve the plot I was after - I need to see more than one syntenic hit. The current plots only show the best one (I assume). I have started using other methods now (pairwise dotplots) but it would still be nice to be able to show a riparian plot. Another option is perhaps circos, but it is rather difficult to use.

Thanks.

jtlovell commented 7 months ago

The problem isn't plotting, but your definition of orthology. Without an outgroup, you'll never see the WGD from orthogroups.