<p>This should be very easy to do. Here are some options:</p>

<ol>
<li><p>Use a tool like <a href="https://www.ebi.ac.uk/about/vertebrate-genomics/software/exonerate" rel="nofollow noreferrer">Exonerate</a> or <a href="https://www.ebi.ac.uk/Tools/psa/genewise/" rel="nofollow noreferrer">GeneWise</a> that can align protein sequences to genomic DNA while attempting to model splice sites etc. </p></li>
<li><p>As you said, blast. There's nothing special in what you describe, really, this is what tBLASTn is for. Just use your protein as the query sequence and tBLASTn as the blast flavor which will take a protein query and match it against all 6 possible reading frames of a DNA database. Set up your assembly as the blast DB, and then blast as you would normally. </p></li>
</ol>

<p>Of the two, using the first approach is much better since it will try to build a valid gene model for you instead of simply finding the regions of high homology as blast will do. </p>

<p>Finally, just for completeness, you can go old school and run a full de-novo gene prediction tool on your assembly. Something like <a href="http://genome.crg.es/software/geneid/" rel="nofollow noreferrer">GeneID</a> for instance. But that's certainly overkill if you're just looking for a single gene.</p>

CC BY-SA 3.0
2017-11-13T12:18:18.610Z
2017-11-13T12:18:18.610Z
terdon
<p>I have two genome assemblies of the same non-model species, call them Assembly 1 (generated from Illumina data) and Assembly 2 (generated from PacBio data).</p>

<p>For Assembly 1, I also have predicted proteome data, generated with <a href="https://evidencemodeler.github.io/" rel="nofollow noreferrer">EVM</a>. Say there is a protein, call it Protein X, for which I have the fasta sequence from predictions, but suspect an error in Assembly 1 (because the Protein X has a long segment of auto-identity). For this reason, I would like to verify the assembly of the locus coding for Protein X in Assembly 2 and/or verify the presence of Protein X in the Assembly 2. </p>

<p>Data in my possession:</p>

<ul>
<li>both genome assemblies (fasta)</li>
<li>raw rna-seq data </li>
<li>reference-based transcriptome assembly for each genome assembly (fasta/gtf)</li>
<li>set of predicted protein sequence based on Assembly 1 (fasta)</li>
</ul>

<p>Things I was thinking of: </p>

<ul>
<li>build the EVM predictions on the Assembly 2 and verify the identity of Protein X by simple fasta identity check </li>
<li>look for the protein sequence in the Assembly 2. I could use something similar to <a href="http://www.bioinformatics.org/sms2/rev_trans.html" rel="nofollow noreferrer">Reverse Translate</a> to get the most likely nucleotide sequence from the input amino-acid sequence of Protein X, then align it to the Assembly 2 to check the coordinates of the locus it originates from - but I am not very sure of the relevance and/or correctness of this procedure</li>
<li>something allowing 'blast' of the protein sequence in the Assembly 2 (but I am not aware of any software that could do that)</li>
</ul>

<p>Any ideas on how I could do that?</p>

<proteins><assembly><validation>
Verify a predicted protein in one genome in a different genome of the same species
# Verify a predicted protein in one genome in a different genome of the same species

(Asked by: username_0 on 2017-11-13)

I have two genome assemblies of the same non-model species, call them Assembly 1 (generated from Illumina data) and Assembly 2 (generated from PacBio data).

For Assembly 1, I also have predicted proteome data, generated with [EVM](https://evidencemodeler.github.io/). Say there is a protein, call it Protein X, for which I have the fasta sequence from predictions, but suspect an error in Assembly 1 (because the Protein X has a long segment of auto-identity). For this reason, I would like to verify the assembly of the locus coding for Protein X in Assembly 2 and/or verify the presence of Protein X in the Assembly 2.

Data in my possession:

* both genome assemblies (fasta)
* raw rna-seq data
* reference-based transcriptome assembly for each genome assembly (fasta/gtf)
* set of predicted protein sequence based on Assembly 1 (fasta)

Things I was thinking of:

* build the EVM predictions on the Assembly 2 and verify the identity of Protein X by simple fasta identity check
* look for the protein sequence in the Assembly 2. I could use something similar to [Reverse Translate](http://www.bioinformatics.org/sms2/rev_trans.html) to get the most likely nucleotide sequence from the input amino-acid sequence of Protein X, then align it to the Assembly 2 to check the coordinates of the locus it originates from - but I am not very sure of the relevance and/or correctness of this procedure
* something allowing 'blast' of the protein sequence in the Assembly 2 (but I am not aware of any software that could do that)

Any ideas on how I could do that?

**Tags:** proteins, assembly, validation

---

## 1 answers

### Answer by username_1 on 2017-11-13. Score: 3 | Accepted answer

This should be very easy to do. Here are some options:

1. Use a tool like [Exonerate](https://www.ebi.ac.uk/about/vertebrate-genomics/software/exonerate) or [GeneWise](https://www.ebi.ac.uk/Tools/psa/genewise/) that can align protein sequences to genomic DNA while attempting to model splice sites etc.
2. As you said, blast. There's nothing special in what you describe, really, this is what tBLASTn is for. Just use your protein as the query sequence and tBLASTn as the blast flavor which will take a protein query and match it against all 6 possible reading frames of a DNA database. Set up your assembly as the blast DB, and then blast as you would normally.

Of the two, using the first approach is much better since it will try to build a valid gene model for you instead of simply finding the regions of high homology as blast will do.

Finally, just for completeness, you can go old school and run a full de-novo gene prediction tool on your assembly. Something like [GeneID](http://genome.crg.es/software/geneid/) for instance. But that's certainly overkill if you're just looking for a single gene.


---

