Comments (7)
@wensulin93 this seems to be a generic task, you may try agat: https://www.biostars.org/p/9465973/
from tesorter.
You may just use the concatenated domain CDS sequences for each LTR, or the region start at the first domain and end at the last domain (CM030788.1 10392696 - 10396473 ), or the entire LTR (CM030788.1:10391873..10397114) sequence.
from tesorter.
Excellent advice, lucky you can see this post! @oushujun
from tesorter.
@wensulin93 this seems to be a generic task, you may try agat: https://www.biostars.org/p/9465973/
Hello, I thought about this question carefully yesterday. It might help you optimize your software.
Because the genome of the selected species is incomplete (including much information missing from the genomic GFF3 file), it seems that extracting CDS is not an easy task.
I looked at the genetic structure of one LTR. It is found that all the sequences generated by the corresponding.dom. Gff3 file are not complete CDS sequences after splicing together, but they cover most of them. The main parts not covered are the beginning (including ATG) or the end or whatever. However, in the LTR sequence I observed, the EXACT division of CDS sequence was not observed in the GFF3 file of the genome, or even directly skipped that section (not shown in the GFF file). This represents an incomplete GFF file for the genome.
Therefore, if there is no new progress, I will splicing the CDS sequences of all domains together and perform synteny analysis with MCscanX. The ultimate goal is to perform synteny analysis. I think stitching together CDS sequences of all domains will also give me the ideal synteny data. @zhangrengang
from tesorter.
I tried the concatenate_domains.py script, but I got an error message like this:
Traceback (most recent call last):
File "/home/manager/miniconda3/bin/concatenate_domains.py", line 33, in
sys.exit(load_entry_point('TEsorter==1.4.1', 'console_scripts', 'concatenate_domains.py')())
File "/home/manager/miniconda3/bin/concatenate_domains.py", line 25, in importlib_load_entry_point
return next(matches).load()
StopIteration
from tesorter.
@wensulin93
Yes, TEsorter do not define the EXACT start, end, or division positions of domains. But these are not neccesary for synteny analyses, so previously I gave you three solutions: the first contains only domains, the second contians almost the full GAG-POL except that the start and end regions may be not incomplete, and the second contains full CDS but also incudes non-coding regions.
Regarding the issue of concatenate_domains.py
, how do you use the script? It should be used like:
concatenate_domains.py rice6.9.5.liban.rexdb.cls.pep RH RT INT > rice6.9.5.liban.rexdb.cls.pep_RT_RH_INT.aln
In this example, the RH, RT and INT domains are aligned seperately and then concatennated together. The LTR-RTs that do not contain all the three domains will be excluded.
from tesorter.
@wensulin93 Yes, TEsorter do not define the EXACT start, end, or division positions of domains. But these are not neccesary for synteny analyses, so previously I gave you three solutions: the first contains only domains, the second contians almost the full GAG-POL except that the start and end regions may be not incomplete, and the second contains full CDS but also incudes non-coding regions.
Regarding the issue of
concatenate_domains.py
, how do you use the script? It should be used like:concatenate_domains.py rice6.9.5.liban.rexdb.cls.pep RH RT INT > rice6.9.5.liban.rexdb.cls.pep_RT_RH_INT.aln
In this example, the RH, RT and INT domains are aligned seperately and then concatennated together. The LTR-RTs that do not contain all the three domains will be excluded.
Hello,
I think I might solve a problem. Tesorter was previously installed with Conda (1.3.0, "Conda Insall Tesorter"). Then I tried to find the "concatenate_domains.py" script. The result is in a "python-scripts" folder in the conda folder. So I typed the command:
~/miniconda3/pkgs/tesorter-1.3.0-py_0/python-scripts/concatenate_domains.py 2_oute.rexdb-plant.cls.pep RT > rt2.aln
Results hint:
"ZSH: / home/manager/miniconda3 / PKGS/tesorter - 1.3.0 - py_0 / python scripts/concatenate_domains.py: bad interpreter: /opt/conda/condabld/tesorter_1604435924377/_h_env_placehold_placehold_placehold_placehold_placehold_placehold_placehold _place: no such file or directory"
And then I realized that the first line seems to be wrong. So I changed it to
"#! /usr/bin/env python"
Then came the reminder:
"Traceback (most recent call last):
File "/home/manager/miniconda3/pkgs/tesorter-1.3.0-py_0/python-scripts/concatenate_domains.py", line 7, in
from .RunCmdsMP import run_cmd
ImportError: attempted relative import with no known parent package"
I didn't find the script ".RunCmdsMP". Finally I found the script in the "modules" folder. I also found that "concatenate_domains.py" in the modules folder only works. Then the instructions were given:
~/miniconda3/pkgs/tesorter-1.3.0-py_0/site-packages/TEsorter/modules/concatenate_domains.py 2_oute.rexdb-plant.cls.pep RT > rt2.aln
The result file was successfully obtained. The conda installation cauesed the problem. Hope it can be reference! Thank you anyway!
from tesorter.
Related Issues (20)
- error HOT 3
- error to get phylogeny tree using LTR_tree.R script HOT 9
- Does TEsorter results only contain positive strain? HOT 5
- keyError HOT 5
- Hi~ TEsorter can be worked in animal genoem ? HOT 2
- ValueError: invalid literal for int() with base 10: '0.5' HOT 8
- How to analyze the effect of transposons on plant traits in de novo transcriptome assembly. HOT 4
- ModuleNotFoundError: No module named 'RunCmdsMP' HOT 5
- Problem with concatenate_domains.py HOT 1
- How to identify the homology (synteny) of LTRs? HOT 10
- Insertion time calculation HOT 4
- Crash when special characters occurred in sequence IDs HOT 3
- Can TEsorter classify Class II elements(DNA transposons) into clade-level? HOT 2
- Assistance with custom installation directory HOT 5
- How to obtain the set of distances between LTRs and their adjacent genes? HOT 3
- TEsorter genome.fasta -genome -p 20 -prob 0.9 HOT 1
- Allocation into lineages for metazoan LTR-RTs HOT 4
- How to merge the TEsorter repeat libraires HOT 18
- get_full_seqs in LTR_retriever.py generate some empty sequences which should be generated. HOT 4
- Exploring the transposition profile of specific LTRs HOT 1
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from tesorter.