Comments (5)
Hi Maksym,
I guess you're talking about the Resolved_Gene_Trees/*.txt that are produced by Orthofinder?
Can you confirm?
If so, the ENSEMBL.buildProteinTrees.py script is not written for that. (This script is used to move duplication nodes that are not well supported in Ensembl trees).
I think I have something that could be used for that (transform nwk to nhx with species and duplication node annotations?) but I have to make some tests before.
best regards,
Alex
from agora.
Hi Alex,
Thanks for answer. Yes, I was talking about Resolved gene trees, which Orthofinder produces.
If you have something that could help with transformation of nwk to nhx with species and duplication node annotations and could share, this will be great. Thank you.
Kind regards,
Maksym Prylutskyi
from agora.
Hi Maksym,
sorry for the delay.
To annotate a nwk with species and duplications, you can use a tool that has been developped by a former PostDoc in the lab:
https://github.com/delehef/chainsaw
you should be able to install it with:
cargo install chainsaw
if you have any problem with chainsaw, open an issue on the chainsaw repository.
To use it with OF resolved gene trees, you will have to first transform the leaf names of the trees to get the corresponding extant species tag:
I've added a script in the dev branch og AGORA repo:
src/import/orthofinder_resolved_gene_trees/addSpeciesAnnotationToOFTrees.py
And here is the procedure that was successfull for me in an Orthofinder test directory (OrthoFinder/Results_May24>):
mkdir AgoraData/Resolved_Gene_trees
python $AGORA/src/import/orthofinder_resolved_gene_trees/addSpeciesAnnotationToOFTrees.py -of_trees=Resolved_Gene_Trees -outdir=AgoraData/Resolved_Gene_trees/
for i in AgoraData/Resolved_Gene_trees/*; do chainsaw $i annotate --species-tree Species_Tree/SpeciesTree_rooted_node_labels.txt ; done
=> you should have annotated nhx trees that could be concatenated in a single forestTree files to use AGORA on it.
I hope this helps,
regards,
Alex
from agora.
Hi Alex,
Sorry for late response. Thank you for script and explanation. You saved me. Thank you again.
Kind regards,
Maksym
from agora.
Hi Maksym,
I'd love any feedback on the method and the Resolved_Gene_trees you use.
The first tests I did on our data were disappointing (a lot of duplications), so if you can tell me what you feel about other data, I'd be grateful.
I don't now if it is due to the OF algorithm or our methodology to annotate the trees... :-(
best regards,
Alex
from agora.
Related Issues (20)
- Error when running agora-basic.py HOT 5
- Agora in living sps HOT 1
- Error Running agora-basic.py: "assert oldName not in seen" HOT 5
- Request for Custom Scripts: Rearrangement Breakpoints filtering and Interchromosomal Rearrangement Identification
- Clarification on Identifying Successive Genomes and Synthetic Blocks between Ancestral Nodes and Extant Genomes in PhylDiag HOT 2
- Creating input for analysis HOT 4
- How to locate the support scores of the AGORA adjacency graph for adjacent genes in ancGenome?
- Orthology groups correspond to HOGs? HOT 2
- Agora-generic pipeline seems to be βstuckβ for several days after adding 1 more genome HOT 4
- AssertionError in ALL.reformatGeneFamilies.py HOT 4
- drawKaryotype trouble HOT 2
- Visualization of results HOT 1
- Is it possible to get ancestral nucleotide sequence? HOT 2
- How to merge SCORPiOs HOT 2
- list indices must be integers or slices, not str HOT 10
- assert oldName not in seen HOT 2
- Too more CARs and how to determine the number of chromosomes HOT 2
- how to deal with many CARs HOT 3
- ValueError: max() arg is an empty sequence HOT 2
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 agora.