__ __ _ __
____ / /_ __ __/ /___ _____(_)________ _____ / /_____ _____
/ __ \/ __ \/ / / / / __ \/ ___/ / ___/ __ `/ __ \/ __/ __ \/ ___/
/ /_/ / / / / /_/ / / /_/ / /__/ / / / /_/ / /_/ / /_/ /_/ / /
/ .___/_/ /_/\__, /_/\____/\___/_/_/ \__,_/ .___/\__/\____/_/
/_/ /____/ /_/
Phylociraptor is a computational framework calculating phylogenomic trees for a specified set of species using different alignment, trimming and tree reconstruction methods. It is very scalable and runs on Linux/Unix machines and servers as well as HPC clusters. Phylociraptor automatically downloads genomes available on NCBI and combines them with additional specified genomes provided by the user. It uses BUSCO to identify single-copy orthologs which are filtered, aligned, trimmed and subjected to phylogenetic inference.
If you use phylociraptor please cite our preprint, which is available here.
Resl Philipp & Hahn Christoph (2023) Phylociraptor - A unified computational framework for reproducible phylogenomic inference (Preprint) bioRxiv doi:10.1101/2023.09.22.558970
Phylociraptor was designed in such a way that it can run on desktop computers (although this is discouraged), solitary linux servers or large HPC clusters. Depending on the system setup, requirements are different:
Local computer or solitary server:
or
On a HPC cluster:
Orthology inference:
Alignment:
Trimming:
Tree inference:
$ conda install -n base -c conda-forge mamba
$ mamba create -c conda-forge -c bioconda -n snakemake snakemake=6.0.2
$ conda activate snakemake
Additional information on how to install snakemake can be found here.
Now we can proceed with installing phylociraptor.
(snakemake) $ git clone --recursive https://github.com/reslp/phylociraptor.git
(snakemake) $ cd phylociraptor
(snakemake) $ ./phylociraptor
Welcome to
__ __ _ __
____ / /_ __ __/ /___ _____(_)________ _____ / /_____ _____
/ __ \/ __ \/ / / / / __ \/ ___/ / ___/ __ `/ __ \/ __/ __ \/ ___/
/ /_/ / / / / /_/ / / /_/ / /__/ / / / /_/ / /_/ / /_/ /_/ / /
/ .___/_/ /_/\__, /_/\____/\___/_/_/ \__,_/ .___/\__/\____/_/
/_/ /____/ /_/
the rapid phylogenomic tree calculator, ver.1.0.0
Usage: phylociraptor <command> <arguments>
Commands:
setup Setup pipeline
orthology Infer orthologs in a set of genomes
filter-orthology Filter orthology results
align Create alignments for orthologous genes
filter-align Trim and filter alignments
modeltest Calculate gene trees and perform model testing
mltree Calculate Maximum-Likelihood phylogenomic trees
speciestree Calculate species trees
njtree Calculate Neighbor-Joining trees
report Create a HTML report of the run
check Quickly check status of the run
util Utilities for a posteriori analyses of trees
-v, --version Print version
-h, --help Display help
Examples:
To see options for the setup step:
./phylociraptor setup -h
To run orthology inferrence for a set of genomes on a SLURM cluster:
./phylociraptor orthology -t slurm -c data/cluster-config-SLURM.yaml
For more information including a short tutorial please refer to the documentation.
A typical run of phylociraptor would look like this:
$ ./phylociraptor setup -t local
$ ./phylociraptor orthology -t local
config.yaml
file:$ ./phylociraptor filter-orthology -t local
$ ./phylociraptor align -t local
config.yaml
file:$ ./phylociraptor filter-align -t local
$ ./phylociraptor modeltest -t local
$ ./phylociraptor njtree -t local
$ ./phylociraptor mltree -t local
$ ./phylociraptor speciestree -t local
$ ./phylociraptor report
To customize the behavior of the pipeline to fit your needs you can edit the config.yaml
file in the data/
folder. Two things are mandatory:
species: data/data.csv
busco:
set: "fungi_odb9"
ausgustus_species: anidulans
You will also need to provide a list of genomes which should be used in your analysis. To do this, edit your data.csv file
The data.csv file should look something like this:
$ cat data.csv
species,web_local,mode
Salpingoeca rosetta,web=GCA_000188695.1,
Coccidioides posadasii,web,
Sclerophora sanguinea,data/assemblies/Sclerophora_sanguinea_Sclsan1_AssemblyScaffolds_Repeatmasked.fasta,
Capsaspora owczarzaki,web=GCA_000151315.2,
Dictyostelium lacteum,web=GCA_001606155.1,
Paraphelidium tribonemae,data/assemblies/EP00158_Paraphelidium_tribonemae.fasta,protein
Synchytrium microbalum,web=GCA_006535985.1,
Nuclearia sp,data/assemblies/Nuclearia_sp_trinity_cdhit-0.95.fasta,transcriptome
Aspergillus nidulans,data/assemblies/Asp_nid_full_proteome.fasta,protein
Stereomyxa ramosa,data/assemblies/Stereomyxa_ramosa_trinity_cdhit-0.95.fasta,transcriptome
The basis of this file can be a CSV file directly downloaded from the NCBI Genome Browser. Just mind the changed header and additional column in the example above. The mandatory columns are the species
and the web_local
column. The first is the species name and the second specifies whether the genome is provided locally (in which case you should specify the path to the assembly) or not (in which case you should specify web).
It is important that the species names correspond exactly to the names under which a genome is deposited at NCBI. Therefore it makes sense to use a downloaded file from the NCBI Genome Browser and add local species to them. However, you can also run the pipeline with only your own assemblies without downloading anything. It is also possible to use transcriptome assemblies or sets of proteins. Please refer to the documentation for more information.
Phylociraptor provides several utilities to investigate tree similarity and plot trees. Please refer to the documentation for additional details.
It is also possible to run phylociraptor inside a docker container. This container has singularity and conda installed so the only requirement is that Docker is properly set up. This is still an experimental feature and we have not tested this extensively. We still recommend using phylociraptor as it is described above, especially when you work on a HPC cluster.
In case you would like to try phylociraptor in Docker you just have to subsitute the phylociraptor
command with phylociraptor-docker
.
$ git clone --recursive https://github.com/reslp/phylociraptor.git
$ cd ./phylociraptor
$ ./phylociraptor-docker