Is that a quote?
What is the recommended way to group transcripts by "gene"? For example, rnaSPAdes returns a predicted gene identifier with the transcript in the identifier.
Is there a way to merge similar transcripts from different assemblies?
I'm sure this uses the graph files in the backend so what I'm asking may be really out of scope. I could use something like CD-HIT but I'm wondering if there was a better way. Maybe a way to use the de brujn graphs together?
2 answers
Transcripts were clustered using CD-HIT (Cluster Database at High Identity with Tolerance) package. Used to remove the shorter redundant transcripts when they were 100% covered by other transcripts with more than 90% identity. The non-redundant clustered transcripts were then designated as unigenes.
If you want a common "baseline" for several assemblies, the only way I could think of would be annotation with a common database. For instance, if your assemblies come from closely related species, you could annotate CDSs with a common reference (e.g. a closely related model species).
For individual assemblies, if you would like to go from transcript-level to gene-level (which has advantages), you could also cluster transcripts using Corset, or Grouper based on shared read support.
Another way forward would be to combine both techniques. For each assembly, you could cluster assembled transcripts using one of the two previously mentioned programs. You could then construct SuperTranscripts using Lace. You could then try and annotate supertranscripts, and compare (likely) homologous genes with one another. Hope that helps.
A suggestion. If you could provide links for the programs mentioned your answer would become more complete. Programs can have similar names and searching with the names above is likely to lead to not-useful-for-science results.
Thanks for the suggestion, but I included links in my answer. They are but a click on the name away :) This is sadly hard to see on biostars when you highlight the names of programs, while also linking to a repository...
Log in to answer this question.