Six representative fungal species and their whole genome assemblies are prepared as a sample dataset.
|Cryptococcus neoformans JEC21||Basidiomycota||
|Saccharomyces cerevisiae S288C||Ascomycota||
|Aspergillus nidulans FGSC A4||Ascomycota||
|Spizellomyces punctatus DAOM BR117||Chytridiomycota||
|Trametes versicolor FP-101664 SS1||Basidiomycota||
|Rhizopus oryzae 99-892||Mucoromycota||
After running following command, you will find all the required files from sample/ directory.
The directory contains following components.
By following this guide, you will
Launch profile module in interactive mode by entering following command on your terminal.
The program will tell you to enter the input file or directory. We want to use the entire directory as an input, therefore you should enter the path where the genomes are lying.
Enter the output directory in the same way.
We want to extract protein markers from the sequences. Type 'PRO' to continue.
We do not want to maintain temporary files. Type 'n'.
Now the pipeline will automatically solve dependency and inform you to enter the path of the binaries remaining unsolved.
After then, the program will ask you to decide the number of CPU threads to use. Allocate appropriate amount of threads considering the hardware specification and your affordable time. In here, 10 threads are allocated as an example.
Then the pipeline will ask for the metadata. In this tutorial, we will use the prepared table with complete information.
Finally, after refusing to configurate the advanced options, the pipeline will display the single-line command to run, and ask for your final confirmation. Enter 'y' to launch the pipeline.
Tip. If you copy and save the command, you can skip the entire interactive process next time.
By launching the pipeline, the process will begin by reading the input files in and preparing queries to process.
After the preparation finishes, you will see the visualized progress bar such as following.
Alphabets in the progress bar, which dynamically change during the extraction process, indicate the current status of the genes from the UFCG set as following.
As the extraction finishes, the program will report the result by the count of genes based on their existence, and write the .ucg file on the output folder.
The pipeline will automatically terminate as it finishes extracting the gene profiles from all six queries. You can check the result files lying in the output directory by entering the following command.
If you have followed the guide through here, you have successfully extracted the core gene profiles from the genome assemblies. You are now ready to proceed.
By following this guide, you will
Launch the module by entering the directory with profiles as an input, labeling them with unique IDs. Since we do not want to perform precise phylogeny here, we will use
The pipeline will automatically detect .ucg files from sample/out/ directory, perform multiple alignment on each core gene, and reconstruct the phylogenetic tree along with the single gene trees on sample/tree/tutorial/ directory.
Run following command to check the contents of the result directory.
Reconstructed phylogenetic tree is written in concatenated.nwk. Printing the Newick tree in a plain text will display its content such as following (unique ID number may differ).
Oops! Since we labelled the leaves with their unique IDs, this tree does not inform us the actual taxonomic relationship in a glance. We have to replace their label with proper metadata to solve this problem.
Launch prune module with following command, indicating it to replace the labels of UFCG tree with taxon names and strain names.
This process will create the replaced UFCG tree as replaced.UFCG.nwk in current directory. Move it to our tutorial directory and check its content.
If you have followed the guide through here, you have successfully reconstructed the phylogenetic tree with core genes labelled with proper taxonomic information. You are now ready to proceed.
MEGA, provided by MEGA software, is used as a visualization toolkit in this tutorial. Use this external link to download the software.
Launch the MEGA application and load up the replaced tree with Analysis - User Tree - Display Newick Trees menu. The tree will be visualized such as following image.
By rearranging the subtrees and labeling them with phyla, the tree can be modified like following.
Considering Spizellomyces punctatus as a root, you can easily see that the classification of remaining species are intact, corresponding to the known phylogenetic relationship of fungal phyla.
This is the end of the tutorial. You are now ready to reconstruct the phylogenetic relationship with your own fungal genome assemblies. Feel free to email us if you have any additional questions to ask, want some technical difficulties to be solved, or wish to give us some kind suggestions for improvement.