Build a core-genome phylogeny

Overview

This tutorial describes step by step the reconstruction of a core genome phylogeny for a given clade of bacteria.

Setup

1. clone eutils repository

   1 # download repo from metagenlab
   2 git clone https://github.com/metagenlab/eutils.git
   3 # move to eutils directory
   4 cd eutils
   5 # add path to pythonpath and path
   6 echo "export PATH=:$PWD"':$PATH' >> ~/.bashrc
   7 echo "export PYTHONPATH=:$PWD"':$PYTHONPATH' >> ~/.bashrc
   8 # reload bashrc
   9 source ~/.bashrc
   

2. install miniconda
3. create a dedicated environment and install orthofinder

   1# create a new environment
   2conda create -n orthofinder
   3source activate orthofinder
   4# install orthofinder
   5conda install orthofinder
   6# install biopython
   7conda install biopython
   

Get table of genome assemblies

The list of available genome assemblies can be downloaded with the script get_ncbi_genome_assemblies.py.

 1usage: get_ncbi_genome_assemblies.py [-h] [-t TAXON_ID] [-c] [-rf] [-rp] [-rs]
 2                                     [-gb] [-md] [-e] [-d] [-ra RANK]
 3                                     [-o OUTPUT]
 4
 5optional arguments:
 6  -h, --help            show this help message and exit
 7  -t TAXON_ID, --taxon_id TAXON_ID
 8                        taxon_idd
 9  -c, --complete        Filter complete genomes (default: False)
10  -rf, --reference      Filter reference genomes (default: False)
11  -rp, --representative
12                        Filter representative genomes (default: False)
13  -rs, --refseq         Filter refseq assemblies (default: False)
14  -gb, --genbank        Filter genbank assemblies (default: False)
15  -md, --exclude_metagenome_derived
16                        Exclude metagenome derived assemblies (default: False)
17  -e, --exclude_unclassified
18                        Exclude species labelled as "unclassified xxx" or as
19                        "xxxx sp. yyy" (e.g. Blautia sp. TF11-31AT)
20                        (default=False)
21  -d, --complete_list   Report complete list of assemblies (no taxonomy
22                        clustering) (default: False)
23  -ra RANK, --rank RANK
24                        Keep one prepresentative/clade of specified taxonomic
25                        <rank>
26  -o OUTPUT, --output OUTPUT
27                        Output name (default: assemblies_<taxid>.tab )

Download gbff files

1df -h