I am working with a non-model species for which another group has published a genome assembly. The assembly is still at the level of scaffolds (N50 = 3244), and no genes have been annotated. For their study, this group used around 100,000 de novo assembled transcripts (from paired end RNA-seq) to assist in the scaffolding.
I'm working with a significantly larger data set of RNA-seq reads (also paired end, 150 bp Illumina reads) that I've assembled into a de novo transcriptome. The transcriptome itself is, needless to say, also larger and more complete (according to BUSCO) than the one the genome study assembled.
I am wondering if I can now use this data to improve the existing genome assembly. Or would this be a fool's errand? If that is possible, could someone perhaps suggest a good pipeline or set of tools for this?
1 answer
Did you try it with MOSGA? It can handle FASTQ RNA-seq data with BRAKER to predict protein-coding genes. Additionally, you can provide an existing genome annotation as GBFF file and it will merge both results.
Select "BRAKER," click on "More" and select "RNA evidence-based". You should there define in the "Settings" if you are providing a hint file or FASTQ file.
You can alternatively host MOSGA on your computer via docker or a Linux (virtual) machine.
Log in to answer this question.
-
Using a different reference transcriptome assembly from different RNA seq libraries (but same speci…
written by ZuelTech •
-
Trimmomatic before or after or nope for genome guided assembly?
written by slin023 •
-
Should I map RNA reads to genome or transcript ?
written by Takuma •
-
Need guidance for RNA_seq analysis of Ecoli with expression vector
written by Morgan S. •
-
Denovo assembly or Reference based assembly for identifying novel lncRNAs
written by newbie
-
Genome assembly contig ordering
written by mittu1602
-
PolyA selection for transcriptome sequencing
written by Mehmet
-
How I can increase the length of contig in denovo assembly
written by ebrahimiet •
-
Any experience with scaffolding using SSPACE tool
written by seta
-
strand-specific transcriptome Oases vs. CLC
written by wd •
AFAIK, transcripts are not usually used for genome assembly. You can use them for gene annotation though. In any case, an assembly with N50 of 3k may be very challenging to work with, so you may consider improving it using more common ways: more sequencing data, long reads, Hi-C/optical/genetic maps etc.
Thank you for the feedback. Do you have any recommendations for tools I could use to annotate the genome as you suggested?
PASA can be a good tool to start with.
Thank you, I'll take a look at
PASAas you suggested.