This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Removing UniVec sequences in Assembly

My first time uploading my de novo transcriptome assembly to TSA. I got flagged with Code(VECTOR_MATCH). I tried to trim the vector using bbduk and cutadapt, but I couldn't seem to make it work. For example, this adapter: 

>gnl|uv|NGB00150.1:1-46 Ambion FirstChoice RLM-RACE 3' RACE adapter
GCGAGCACAGAATTAATACGACTCACTATAGGTTTTTTTTTTTTVN

It matches to my example sequence (BOLD) according to VecScreen. 

> sample
GCAAAGAAGCATTTTGGCAAAAAATTGCGTAATATTCTGCCGTATGTTACTGCAATGTACACGTTTATAA
TTATTGTAATAAGAATGTCTCATATTGCCTGCTTGATGTGGCAGGGTCACTTGTCAAGTGAGGAAAAGTC
ACAGTGTGAGGACTGTCTATAAAAATTTAGGCATCATATTAAAATGTGTGGATGCCTTATTGTATAGAAT
ATTTCAAATTTTGCAAAATTTGAACAAAGCATATAAAATAAAAGGAACGAAATTGAAAAAAAAAAAAAAA
A**GTCGTATTAATTCTGTGCTCG**

My problem is: sequence trimmers could not recognize the adapter? I had to manually reverse complement the adapter and remove the "extra" sequences just to have that exact match on my sequence. Had it been 5 sequences, I can manually remove them; but, more than a thousand sequences were flagged. I am not sure what am I doing wrong. This pandemic is making me too exhausted to read more bioinformatics...

assembly rna-seq tsa cutadapt bbduk
• 2,037 views • link
updated by GenoMax GenoMax • written by jfo • jfo

2 answers

I tried to trim the vector using bbduk and cutadapt, but I couldn't seem to make it work.

Using bbduk.sh that sequence is definitely detected and trimmed. Can you let us know how you used bbduk.sh?

$ more test.fa
>test
GCAAAGAAGCATTTTGGCAAAAAATTGCGTAATATTCTGCCGTATGTTACTGCAATGTACACGTTTATAATTATTGTAATAAGAATGTCTCATATTGCCTGCTTGATGTGGCAGGGTCACTTGTCAAGTGAGGAAAAGTCACAGTGTGAGGACTGTCTATAAAAATTTAGGCATCATATTAAAATGTGTGGATGCCTTATTGTATAGAATATTTCAAATTTTGCAAAATTTGAACAAAGCATATAAAATAAAAGGAACGAAATTGAAAAAAAAAAAAAAAAGTCGTATTAATTCTGTGCTCG

$ bbduk.sh in=test.fa literal=GCGAGCACAGAATTAATACGACTCACTATAGGT ktrim=r out=new.fa k=10

$ more new.fa
>test
GCAAAGAAGCATTTTGGCAAAAAATTGCGTAATATTCTGCCGTATGTTACTGCAATGTACACGTTTATAA
TTATTGTAATAAGAATGTCTCATATTGCCTGCTTGATGTGGCAGGGTCACTTGTCAAGTGAGGAAAAGTC
ACAGTGTGAGGACTGTCTATAAAAATTTAGGCATCATATTAAAATGTGTGGATGCCTTATTGTATAGAAT
ATTTCAAATTTTGCAAAATTTGAACAAAGCATATAAAATAAAAGGAACGAAATTGAAAAAAAAAAAAAAA

If you used the full adapter sequence then more sequence is removed.

$ bbduk.sh in=test.fa literal=GCGAGCACAGAATTAATACGACTCACTATAGGTTTTTTTTTTTTVN overwrite=t ktrim=r out=new.fa k=10

$ more new.fa
>test
GCAAAGAAGCATTTTGGCAAAAAATTGCGTAATATTCTGCCGTATGTTACTGCAATGTACACGTTTATAA
TTATTGTAATAAGAATGTCTCATATTGCCTGCTTGATGTGGCAGGGTCACTTGTCAAGTGAGGAAAAGTC
ACAGTGTGAGGACTGTCTATAAAAATTTAGGCATCATATTAAAATGTGTGGATGCCTTATTGTATAGAAT
ATTTCAAATTTTGCAAAATTTGAACAAAGCATA
• 0 views • link
written by GenoMax GenoMax

This worked. I used k=21. I am not sure why. But, thanks....

• 0 views • link
written by jfo • jfo

Please consider accepting this answer (green checkmark) to provide closure to this thread.

• 0 views • link
written by GenoMax GenoMax

We find adapters in genomes all the time which is terrible for metagenomics.

This simple tool might help you to at least diagnose your problem and or replace the adapters with NNNs https://github.com/colindaven/blacklister

Try additional adapter trimmers like -trimmomatic -fastp

and after adapter trimming really check the results using FASTQC and multiqc.

• 0 views • link
written by colindaven colindaven

Log in to answer this question.

More posts like this