Commit f16887b8 authored by peguerin's avatar peguerin
Browse files

readme update

parent be9d281c
# edna--smk_rapidrun_obitools
snakemake obitools workflow to process rapidrun eDNA data
## eDNA-seq Metabarcoding annotation pipeline
**eDNA-seq Metabarcoding annotation** is a bioinformatics pipeline built using Snakemake, a workflow tool to run and manage tasks in any execution environment.
## Introduction
**eDNA-seq Metabarcoding annotation** is specifically used for the analysis of environmental DNA metabarcoding NGS data, demultiplexing, filtering, dereplicating and taxonomic assignment (NCBI taxonomy and custom taxonomy).
This pipeline has been initially tested with marine environmental DNA samples, using molecular markers such as Vert01 ([Riaz et al. 2011](, Teleo01 ([Valentini et al. 2016](, Chond01 or Mamm01 ([Taberlet et al. 2018](10.1093/oso/9780198767220.001.0001)). The workflow should work with any organisms and environment. It is proven for large-scale data analysis.
## Method
This workflow performs step by step analysis of metabarcoding data using [obitools](
The workflow assembles paired-end reads together (illuminapairedend), applies demultiplexing (ngsfilter), dereplicates sequences (obiuniq), detect and removes chimera and PCR clones (obiclean; obigrep) and assigns taxonomy to each representative sequences (NCBI taxonomy and custom taxonomy; ecotag). Ultimately, abundance distribution among samples of representative sequences is generated as a table.
## Workflow
![workflow obitools](
## Quick Start
See [Installion section]( for installation instructions.
Download example data with:
curl -JLO; tar zfxz test_rapidrun_data.tar.gz -C resources/test/
* `./resources/test/test_rapidrun_data/`: this folder contains a reference database for 4 markers (Teleo01; Mamm01; Vert01; Chond01), NGS metabarcoding raw data, required metadata to handle demultiplexing on RAPIDRUN format.
If you have installed **eDNA-seq Metabarcoding annotation**, you can run the example data with:
snakemake --configfile config/config_test_rapidrun.yaml --cores 4 --use-conda
This will generate outputs into `./results` folder.
Final tables are available as:
results/projet1_chond_run1_ecotag_ncbi.csv results/projet1_teleo_run5_ecotag_ncbi.csv results/projet1_chond_run2_ecotag_ncbi.csv
results/projet1_teleo_run2_ecotag_ncbi.csv results/projet1_vert_run3_ecotag_ncbi.csv results/projet1_mamm_run4_ecotag_ncbi.csv
## Next steps
Now that you've gotten the example to work, see the [wiki]( to navigate to the more detailed descriptions and instructions for exploring your own data.
## Software Requirements
* python3
* python3 required modules included: pandas, biopython, datrie, ConfigArgParse, appdirs, gitdb2, snakemake
## Credits
**eDNA-seq Metabarcoding annotation** was coded and written by Virginie Marques and Pierre-Edouard Guerin.
We thank the following people for their help in the development of this pipeline: Agnes Duhamet, Alice Valentini, Apolline Gorry, Bastien Mace, David Mouillot, Emilie Boulanger, Laetitia Mathon, Laura Benestan, Stephanie Manel, Tony Dejean.
## Contributions and Support
:bug: If you are sure you have found a bug, please submit a bug report. You can submit your bug reports on Gitlab [here](
Visit project [wiki]( for documentation.
\ No newline at end of file
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment