hAMRonization is a project aiming at the harmonization of output file formats of antimicrobial resistance detection tools. This is a workflow acting as a proof of concept test-case for the hAMRonization parsers.
Specifically, this runs a set of AMR gene detection tools against a set of contigs/reads and uses hAMRonization
to collate the results in a single unified report.
The following tools are currently included:
- abricate
- AMRFinderPlus
- ariba
- Groot
- RGI (for complete and draft genomes)
- RGI BWT (for metagenomes)
- staramr
- resfams
- staramr
- Resfinder
- sraX
- DeepARG (requires singularity)
- CSSTAR
- AMRplusplus
- SRST2
- KmerResistance
Excluded tools:
- mykrobe (needs variant specification to be parseable)
- pointfinder (needs variant specification to be parseable)
- SEAR, ARG-ANNOT (no longer downloadable)
- RAST/PATRIC (not easily runnable on CLI)
- Single organism/or resistance tools (e.g. Kleborate, LREfinder, SSCmec Finder, U-CARE, ARGO)
- shortBRED, ARGS-OAP (rely on usearch which isn't open-source)
First clone this repository:
git clone https://github.com/pha4ge/hAMRonization_workflow
This pipeline depends on snakemake, conda, build-essentials, git, zlib-dev, and unzip. If you have conda installed, please run:
conda env create -n hamronization_workflow --file envs/hamronization_workflow.yaml
and
conda activate hamronization_workflow
All further dependencies will be installed via conda on execution.
If you want to run DeepARG
you need to have a working singularity
install on your system and invoke --use-singularity --singularity-args "-B $PWD:/data"
when running snakemake (otherwise comment out this input to the cleanup rule in the Snakefile
).
To execute the pipeline, navigate to the cloned repository, edit the config (config/config.yaml
) and input details (config/isolate_list.txt
) for your purposes.
Execute the following substitution the value for --jobs
as needed:
snakemake --configfile config/config.yaml --use-conda --conda-frontend mamba --jobs 2 --use-singularity --singularity-args "-B $PWD:/data"
To test the pipeline follow the above installation instructions and execute on the test data set:
snakemake --configfile config/test_config.yaml --use-conda --conda-frontend mamba --jobs 1 --use-singularity --singularity-args "-B $PWD:/data"
Alternatively, the workflow can be run using docker. Given the collective quirks of the bundled tools this will probably be easier for most users.
Unfortunately, deeparg is only really runnable as a container, and snakemake uses singularity, the docker version has to be run in a privileged manner i.e. docker run --privileged
.
If you are unable to run docker in privileged mode then you can just comment out the deeparg target in the main Snakefile
(expand("results/{sample}/deeparg/output.mapping.ARG", sample=samples.index),
).
First get the docker container:
docker pull finlaymaguire/hamronization:1.0.1
You can execute it in a couple of ways but the easiest is to just mount the folder containing your reads and running it interactively:
docker run -it --privileged -v $HOST_FOLDER_CONTAINING_ISOLATES:/data finlaymaguire/hamronization:1.0.1 /bin/bash
If our isolate data is in ~/isolates
the command to interactively run this container and get a bash terminal would be:
docker run -it --privileged -v ~/isolates:/data finlaymaguire/hamronization:1.0.1 /bin/bash
Then point your sample_table.tsv
to that folder, entries for this example would be:
species biosample assembly read1 read2
Mycobacterium tuberculosis SAMN02599008 /data/SAMN02599008/GCF_000662585.1.fna /data/SAMN02599008/SRR1180160_R1.fastq.gz /data/SAMN02599008/SRR1180160_R2.fastq.gz
Mycobacterium tuberculosis SAMN02599009 /data/SAMN02599009/GCF_000662586.1.fna /data/SAMN02599009/SRR1180161_R1.fastq.gz /data/SAMN02599009/SRR1180161_R2.fastq.gz
Then specify your config.yaml
to use this sample_table.tsv
and execute the pipeline from bash in the container by activating the top-level environment:
conda activate hamronization_workflow
Then the workflow:
snakemake --configfile config/config.yaml --use-conda --cores 6 --use-singularity --singularity-args "-B $PWD:/data"
WARNING You will have to extract your results folder (e.g. cp results /data
for the example mounted volume) from the container if you wish to use them elsewhere.
Note: kma/kmerresistance fails without explanation in the container (possibly zlib related, although adding the zlib headers didn't solve this). It is commented out for now.
Following datasets are currently used for result file generation:
organism Biosample Assembly Run
Salmonella enterica SAMN13012778 GCA_009009245.1 SRR10258315
Salmonella enterica SAMN13064234 GCA_009239915.1 SRR10313698
Salmonella enterica SAMN10872197 GCA_007657735.1 SRR8528923
Salmonella enterica SAMN13064249 GCA_009239785.1 SRR10313716
Salmonella enterica SAMN07255713 GCA_009439415.1 SRR5921214
Salmonella enterica SAMN03098832 GCA_006629605.1 SRR1616829
Klebsiella pneumoniae SAMN02927805 GCA_004302785.1 SRR1561295
Salmonella enterica SAMEA6058467 GCA_009625195.1 ERR3581801
E. coli SAMN05980528 GCA_004268245.1 SRR4897319
Mycobacterium tuberculosis SAMN02599008 GCA_000662585.1 SRR1182980 SRR1180160
Mycobacterium tuberculosis SAMN02599179 GCA_000665745.1 SRR1172848 SRR1172873
Mycobacterium tuberculosis SAMN02599095 GCA_000706105.1 SRR1173728 SRR1173217
Mycobacterium tuberculosis SAMN02599061 GCA_000663625.1 SRR1175151 SRR1172938
Mycobacterium tuberculosis SAMN02598983 GCA_000654735.1 SRR1174279 SRR1173257
Links to data and corresponding metadata need to be stored in a tab separated sample sheet with the following columns:
species biosample assembly reads read1 read2
The results generated on the aforementioned datasets can be retrieved here.
Please consult the PHA4GE project website for questions.
For technical questions, please feel free to consult:
- Finlay Maguire <finlaymaguire (at) gmail.com>
- Simon H. Tausch <Simon.Tausch (at) bfr.bund.de>