Demiulge: DEMIx Universal Lineage and Genotype with amplicon sequencing.

Demiulge is a tool for deconvolution the strain abundance of universal species in mixed samples (e.g. wastewater). The old version of this tool is NextVpower.

"Deconvolution" means "demixing", analyzing the relative abundance of lineage, genotype and phylogeny of a certain taxonomy in one sample.

What can Demiulge demix:

• Human Adenovirus (HAdV)
• Human Astrovirus (HAstV)
• Norovirus
• Rotavirus
• Pathogens in NextClade, using barcodes extracted by BarcodeExtracter.py
• ...

Uusing BarcodeBuilder.py, you can locally build barcodes of your concerned lineages of species with provided sequences. Then you can demix your high-throughput targeted-sequenced samples using Demiulge.

Requirements:

For demixing, in python3:

• numpy
• pandas
• cvxpy

For building phylogenetic trees:

• biopython (in python3)
• mafft
• trimal
• iqtree

For processing NGS data:

• bwa
• minimap2 (for long reads)
• samtools
• freebayes
• bcftools

Usage:

1. Demix from input sample mutation table file (sites_rate.tsv), and save result to output (demix_result.tsv):

python Demiulge.py -i sites_rate.tsv -o demix_result.tsv

2. Demix from input *.vcf files under a folder (vcfs/), filter mutation sites with mutation rate lower than 0.1, filter mutation sites with depth lower than 10, and save result to output (demix_result_all.tsv):

python Demiulge.py -i vcfs/ -r 0.1 -d 10 -o demix_result_all.tsv

3. Add annotation to sample mutation table (sites_rate.tsv, converted from vcf files) according to variation annotation table (varAnno.tsv) and demix:

python Demiulge.py -i vcfs/ -o demix_result_vcf_all.tsv --ann_file varAnno.tsv --ann_vcsample annotated_sites_rate_all.tsv

4. Demix and add meta information of lineages (lineage_metadata.tsv) to demix results:

python Demiulge.py -i vcfs/ -o demix_result_vcf_all.tsv --meta_file lineage_metadata.tsv

5. Demix, save result in result.tsv, and save middle data to files/folder:

python Demiulge.py -i vcfs/ -o demix_result_vcf_all.tsv --vcsample sites_rate_all.tsv --fbarcodes filtered_barcode_matrix/ --fsample filtered_sites_rate_all.tsv --potentials potential_sites_all.tsv

Please see detailed usage by running python Demiulge.py -h or in source code

An example pipeline for demixing HAdV-F is shown in analysis pipeline

Publications

This project was not published yet, but you can still have a try on your amplicon sequencing data of target species.

License

This project is licensed under the MIT License - see the LICENSE file for details.