A bioinformatics tool designed to enable rapid recovery of high-quality draft genomes of microbial eukaryotes from various environmental metagenomic samples.
- reference-independent and cultivation-free
- separate reads or contigs into five groups: Bacteria, Archaeal, Eukaryotic, Viral, and Unknown.
- generate a fasta file with Euk and unknown contigs for binning
Schematic representation of eukfinder workflows. eukfinder is a taxonomic classification-based bioinformatics approach to retrieve microbial eukaryotic nuclear and mitochondrial genomes from WGS metagenomic sequencing data. eukfinder has two different workflows based on the input files:
- (a) using Illumina short reads, it first classifies Illumina reads into 5 distinct taxonomic categories (Archaeal, Bacterial, Viral, Eukaryotic, and Unknown), after assembling the Eukaryotic and Unknown reads together, second round of classification and supervised binning, it will output Eukaryotic nuclear and mitochondrial genomes.
- (b) using MAG assembled contigs or long-read sequencing data generated by Nanopore or Pacbio platforms. It goes through one round of classification to select Eukaryotic and Unknown contigs, and after supervised binning generates Eukaryotic nuclear and mitochondrial genomes.
Anaconda or miniconda required*
conda create -n eukfinder -c bioconda eukfinderDefault reference databases can be downloaded from Eukfinder Databases
- Plast Database
- Centrifuge Database
- acc2tax Database
- Human Genome for read decontamination
- Read Adapters for Illumina sequencing
./download_db.shUsers can flexibly customize the reference data (see here)
source activate eukfinderSee Wiki for detailed description
eukfinder read_prep
Run Trimmomatic to remove low-quality reads, and adaptors Run Bowtie2 to remove host reads Run Centrifuge for the first round of classification
eukfinder read_prep [-h] --r1 R1 --r2 R2 -n THREADS -i ILLUMINA_CLIP
--hcrop HCROP -l LEADING_TRIM -t TRAIL_TRIM --wsize
WSIZE --qscore QSCORE --mlen MLEN --hg HG -o
OUT_NAME --cdb CDB
eukfinder short_seqs
eukfinder short_seqs [-h] --r1 R1 --r2 R2 --un UN -o OUT_NAME -n
NUMBER_OF_THREADS -z NUMBER_OF_CHUNKS -t
TAXONOMY_UPDATE -p PLAST_DATABASE -m PLAST_ID_MAP
[-p2 ANCILLARY_PLAST_DATABASE]
[-m2 ANCILLARY_PLAST_ID_MAP]
[--force-pdb FORCE_PDB] -a ACC2TAX_DATABASE --cdb
CDB -e E_VALUE --pid PID --cov COV --max_m MAX_M
--mhlen MHLEN --pclass PCLASS --uclass UCLASS
eukfinder long_seqs
eukfinder long_seqs [-h] -l LONG_SEQS -o OUT_NAME --mhlen MHLEN --cdb
CDB -n NUMBER_OF_THREADS -z NUMBER_OF_CHUNKS -t
TAXONOMY_UPDATE -p PLAST_DATABASE -m PLAST_ID_MAP
-a ACC2TAX_DATABASE -e E_VALUE --pid PID --cov COV
Contributions are what makes the open-source community such an amazing place to learn, inspire, and create. Any contributions you make are greatly appreciated.
If you have a suggestion that would make this better, please fork the repo and create a pull request. You can also simply open an issue with the tag "enhancement". Don't forget to give the project a star! Thanks again!
Zhao, D., Salas-Leiva, D.E., Williams, S.K., Dunn, K.A. and Roger, A.J., 2023. Eukfinder: a pipeline to retrieve microbial eukaryote genomes from metagenomic sequencing data. bioRxiv, pp.2023-12.