<p>I want to recreate the UGENE workflow proposed in <a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5348765/" rel="nofollow noreferrer"><em>Convergence of retrotransposons in oomycetes and plants</em></a>. After being able to understand the data and hmmsearch <a href="https://bioinformatics.stackexchange.com/q/3039/2026">Here</a> and <a href="https://bioinformatics.stackexchange.com/q/3036/2026">Here</a> and getting most of the result I want to try to recreate the workflow in UGENE.</p>

<p>I manged to implement the workflow for the search <a href="https://i.sstatic.net/5VTqt.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/5VTqt.png" alt="Ugene hmm search Workflow"></a>
However in the latter parts I struggle to implement the steps of the algorithm</p>

<blockquote>
  <p>sequences surrounding the regions of significant similarity to the aRNH profile were expanded, when possible, to 10,000 bp in both directions</p>
</blockquote>

<ul>
<li>How can sequences expanded in UGENE?</li>
<li>How can implement the next steps also? for step three </li>
</ul>

<blockquote>
  <p>Third, the enlarged sequences were screened for the presence of significant similarity to RT domains of non-LTR-RTs and LTR-RTs HMM profiles using hmmsearch.</p>
</blockquote>

<p>I assume it is very simmilar to the first step in building the right profile and running hmmsearch but the following steps are still unclear to me. </p>

CC BY-SA 3.0
2017-12-14T08:10:14.477Z
Expanding means looking outside the original sequence (in this case by 10kb), so you need to find the position of the sequences and look in the genome for a larger range. Could you clarify what problems do you have with the step three? To me it means to use hmmsearch with specific profiles to search in those sequences
llrs
CC BY-SA 3.0
2017-12-14T09:58:32.257Z
Well I thought that the search result of step one returns the loction of aRNH  then from that location expand in the old data by 10k bp and hmmsearh with the new LTR-RT and non-LTR-RT profiles. but how is the exapnsion from the search result done that was unclear to me.
A.Dumas
<sequence-analysis><hidden-markov-models>
How to rebuild the workflow for mining in oomynete genomes in UGENE
# How to rebuild the workflow for mining in oomynete genomes in UGENE

(Asked by: username_0 on 2017-12-12)

I want to recreate the UGENE workflow proposed in [*Convergence of retrotransposons in oomycetes and plants*](https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5348765/). After being able to understand the data and hmmsearch [Here](https://bioinformatics.stackexchange.com/q/3039/2026) and [Here](https://bioinformatics.stackexchange.com/q/3036/2026) and getting most of the result I want to try to recreate the workflow in UGENE.

I manged to implement the workflow for the search [![Ugene hmm search Workflow](https://i.sstatic.net/5VTqt.png)](https://i.sstatic.net/5VTqt.png)
However in the latter parts I struggle to implement the steps of the algorithm

> sequences surrounding the regions of significant similarity to the aRNH profile were expanded, when possible, to 10,000 bp in both directions

* How can sequences expanded in UGENE?
* How can implement the next steps also? for step three

> Third, the enlarged sequences were screened for the presence of significant similarity to RT domains of non-LTR-RTs and LTR-RTs HMM profiles using hmmsearch.

I assume it is very simmilar to the first step in building the right profile and running hmmsearch but the following steps are still unclear to me.

**Tags:** sequence-analysis, hidden-markov-models
> **Comments:**
> **username_1**: Expanding means looking outside the original sequence (in this case by 10kb), so you need to find the position of the sequences and look in the genome for a larger range. Could you clarify what problems do you have with the step three? To me it means to use hmmsearch with specific profiles to search in those sequences  
> **username_0**: Well I thought that the search result of step one returns the loction of aRNH  then from that location expand in the old data by 10k bp and hmmsearh with the new LTR-RT and non-LTR-RT profiles. but how is the exapnsion from the search result done that was unclear to me.  

---

## 0 answers

