aucome 0.5

Automatic Comparison of Metabolism

AuCoMe: Automatic Comparison of Metabolism

Workflow to reconstruct multiple metabolic networks in order to compare them.

Table of contents

• Installation

  • Dependencies

  • Docker

  • pip

• Usage

  • Initialization

  • Check command

  • Renconstruction command

  • Orthology command

  • Draft command

  • Workflow command

  • Analysis command

  • Compare command

Installation

Dependencies

This tool needs:

• Orthofinder (which needs Diamond and FastME)

• Pathway Tools (which needs Blast)

And some python package:

• biopython

• docopt

• eventlet

• intervene (which needs R package UpSetR)

• lxml

• matplotlib

• mpwt

• pandas

• padmet

• padmet-utils

• requests

• seaborn

To run annotation based reconstruction, you need to install Pathway Tools. This tool is available at the Pathway Tools website. A command in the package install the tools:

aucome --installPWT=path/to/pathway/tools/installer

You can also provide an option to this commande: –ptools=ptools_path

This option let you choose the path where the ptools-local folder will be installed. PGDBs created by Pathway-Tools are stored in this folder.

Docker

From git repository:

git clone https://github.com/AuReMe/aucome.git

cd aucome

docker build -t "my_image".

pip

If you have installed all the dependencies, you can just install acuome with:

pip install aucome

Usage

Initialization

You have to create the working folder for AuCoMe, with the –init argument:

aucome --init=run_ID [-v]

This command will create a folder name “run_ID” inside the working folder. In this “run_ID” folder, the command will create all the folders used during the analysis.

run_ID
├── analysis
        ├── group_template.tsv
        ├──
├── annotation_based
        ├── PADMETs
                ├──
        ├── PGDBs
                ├──
        ├── SBMLs
                ├──
├── config.txt
├── logs
        ├──
├── model_organisms
        ├──
├── networks
        ├── PADMETs
                ├──
        ├── SBMLs
                ├──
├── orthology_based
        ├── Orthofinder_WD
                ├──
├── studied_organisms
        ├──

analysis will store the result of padmet analysis.

annotation_based contains three sub-folders. The folder PGDBs will contain all the results from Pathway-Tools (in dat format). These results will be stored in padmet and sbml inside PADMETs and SBMLs.

config.txt contains numerous paths used by the script.

model_organisms contains the model organisms you want to use for the orthology. In this folder you put a new folder with the name of the species and in this folder you put the proteome and the sbml of the metabolic network of your species. Proteome and metabolic network names must be the same than the name of the folder.

├── model_organisms
        ├── A_thaliana
                ├── A_thaliana.fasta
                ├── A_thaliana.sbml

networks will contain all the metabolic network created by aucome in padmet format.

orthology_based contains one folder Orthofinder_WD. This folder will contain all the run of Orthofinder.

studied_organisms: you put all the species that you want to studies in this folder. For each species you create a folder and in this folder you put the genbank file of this species. Like for model_organisms, file and folder must have the same name. And the genbank file must end with a ‘.gbk’.

├── studied_organisms
        ├── species_1
                ├── species_1.gbk
        ├── species_2
                ├── species_2.gbk

Once you have put your species in the studied_organisms folder and the model in model_organisms, a check must be done on the data using:

Check command

aucome check --run=run_ID [--cpu=INT] [-v]

This command will check if there is no character that will cause trouble. It will also create the proteome fasta file from the genbank.

Also, this command will fill the ‘all’ row of analysis/group_template.tsv, with all the species from the studied_organisms folder.

And for the annotation_based folder, if PGDBs contains folder, it will create the padmet and the sbml corresponding to these draft in PADMETs and SBMLs.

Renconstruction command

A run of Pathway-Tools can be launched using the command:

aucome reconstruction --run=run_ID [--cpu=INT] [-v]

├── annotation_based
        ├── PADMETs
                ├── output_pathwaytools_species_1.padmet
                ├── output_pathwaytools_species_2.padmet
        ├── PGDBs
                ├── species_1
                        ├── PGDB dat files
                        ├── ...
                ├── species_2
                        ├── PGDB dat files
                        ├── ...
        ├── SBMLs
                ├── output_pathwaytools_species_1.sbml
                ├── output_pathwaytools_species_2.sbml
├── logs
        ├── log_error.txt
        ├── resume_inference.tsv

Using the package mpwt, it will create the input file for Pathway-Tools inside studied_organisms and if there is no error, it will create for each species inside this folder a folder inside PGDBs containing all the dat files ofthe draft metabolic network.

Orthology command

Orthofinder can be launched using:

aucome orthology --run=run_ID [-S=STR] [--orthogroups] [--cpu=INT] [-v]

├── orthology_based
        ├── Orthofinder_WD
                ├── species_1
                        ├── output_orthofinder_from_species_2.sbml
                ├── species_2
                        ├── output_orthofinder_from_species_1.sbml
                ├── Orthofinder_WD
                        ├── species_1.faa
                        ├── species_2.faa
                        ├── OrthoFinder
                                ├── Results_MonthDay
                                        ├── Orthogroups
                                        ├── Orthologues
                                        ├── ..

Then the proteome from the studied organisms and from the models will be moved to the Orthofinder_WD folder and orthofinder will be launch on them. Orthofinder result will be in this folder and in orthology_based, there will be all the metabolic network reconstructed from orthology.

Draft command

Then you can merge all the metabolic network with:

aucome draft --run=run_ID [--cpu=INT] [-v]

├── networks
        ├── PADMETs
                ├── species_1.padmet
                ├── species_2.padmet
        ├── SBMLs
                ├── species_1.sbml
                ├── species_2.sbml

This will output the result inside the networks folder.

Workflow command

You can launch the all workflow with the command:

aucome workflow --run=ID [-S=STR] [--orthogroups] [--cpu=INT] [-v]

Analysis command

You can launch group analysis with the command:

aucome analysis --run=ID [--cpu=INT] [-v]

You must write the groups of species that you want to analyze in the analysis/group_template.tsv file: The first line of the file contains ‘all’ (it will launch the analysis on all the species).

When you create the repository with –init, the file will only contain ‘all’ row:

all

After the check (with check or workflow command), it will add all the species that you have in your studied_organisms folder:

all

species_1

species_2

species_3

species_4

Then you can create a new row to add another group. The name of the group is in the first column. Then for each species you add a column with the species name. You must at least give 2 species.

Example:

all	species_1	species_2	species_3	species_4
group_1	species_1	species_2
group_2	species_1	species_2	species_4

This script will create one folder for each group:

├── analysis
        ├── group_template.tsv
        ├── all
                ├──
        ├── group_1
                ├──
        ├── group_2
                ├──

Compare command

You can launch group analysis with the command:

aucome analysis --run=ID [--cpu=INT] [-v]

This script will read the group_template.tsv file and create a folder containing an upset graph comparing the group that you selected:

├── analysis
        ├── group_template.tsv
        ├── upgset_graph
                ├── genes.csv
                ├── Intervene_upset.R
                ├── Intervene_upset.svg
                ├── Intervene_upset_combinations.txt
                ├── metabolites.csv
                ├── pathways.csv
                ├── reactions.csv
                ├── tmp_data
                        ├──