Skip to content

Marker Search

Summary

RNA virus marker protein search - using pre-made/user-supplied DBs.

Description

Most pre-made DBs are based on RdRp domain (except for geNomad). Input can be nucleotide contigs or amino acid seqs. If nucleotide, by default all contigs will be translated to six end-to-end frames (with stops replaced by X), or into ORFs called by pyrodigal (meta) or callgenes.sh

Pre-compiled options are:

• NeoRdRp2.1

    GitHub: https://github.com/shoichisakaguchi/NeoRdRp  | Paper: https://doi.org/10.1264/jsme2.ME22001

• RVMT

    GitHub: https://github.com/UriNeri/RVMT  | Zenodo: https://zenodo.org/record/7368133  |  Paper: https://doi.org/10.1016/j.cell.2022.08.023

• RdRp-Scan

    GitHub: https://github.com/JustineCharon/RdRp-scan  |  Paper: https://doi.org/10.1093/ve/veac082

        ⤷ (which IIRC incorporated PALMdb, GitHub: https://github.com/rcedgar/palmdb, Paper: https://doi.org/10.7717/peerj.14055

# • TSA_Olendraite (TSA_2018)

#     Data: https://drive.google.com/drive/folders/1liPyP9Qt_qh0Y2MBvBPZQS6Jrh9X0gyZ?usp=drive_link  |  Paper: https://doi.org/10.1093/molbev/msad060

#     Thesis: https://www.repository.cam.ac.uk/items/1fabebd2-429b-45c9-b6eb-41d27d0a90c2
• Pfam_RTs_RdRp

    RdRps and RT profiles from PFAM_A v.37 --- PF04197.17,PF04196.17,PF22212.1,PF22152.1,PF22260.1,PF05183.17,PF00680.25,PF00978.26,PF00998.28,PF02123.21,PF07925.16,PF00078.32,PF07727.19,PF13456.11
    Data: https://ftp.ebi.ac.uk/pub/databases/Pfam/releases/Pfam37.0/Pfam-A.hmm.gz | Paper https://doi.org/10.1093/nar/gkaa913
• geNomad

    RNA virus marker genes from geNomad v1.9 --- https://zenodo.org/records/14886553
For custom path, either an .hmm file, a directory with .hmm files, or a folder with MSA files (which would be used to build an HMM DB).
Please cite accordingly based on the DBs you select.

Usage

rolypoly marker-search [OPTIONS]

Options

  • -i, --input: Input fasta file. Preferably nucleotide contigs, but you can provide amino acid input too (the script would skip 6 frame translation) (type: TEXT; required; default: Sentinel.UNSET)
  • -o, --output: Path to output directory. Note - if multiple DBs are used and the resolve-mode is none, multiple outputs are made (DB name appended as suffix). (type: TEXT; default: /clusterfs/jgi/scratch/science/metagen/neri/code/rolypoly/marker_search_out)
  • -rm, --resolve-mode: How to deal with regions in your query that match multiple profiles? - merge: all overlapping hits are merged into one range - one_per_range: one hit per range (ali_from-ali_to) is reported - one_per_query: one hit per query sequence is reported - split: each overlapping domain is split into a new row - drop_contained: hits that are contained within (i.e. enveloped by) other hits are dropped. - none: no resolution of overlapping hits is performed. NOTE - EXPECT A POTENTIALLY LARGE OUTPUT - simple: heuristic/personal observation based - chains drop_contained output with split mode. (type: CHOICE; default: simple)
  • -mo, --min-overlap-positions: Minimal number of overlapping positions between two intersecting ranges before they are considered as overlapping (used in some resolve_mode(s) (type: INTEGER; default: 10)
  • -ie, --inc-evalue: Maximal e-value for including a domain match in the results (type: FLOAT; default: 0.05)
  • -s, --score: Minimal score for including a domain match in the results (type: INTEGER; default: 20)
  • -am, --aa-method: Method to translate nucleotide sequences into amino acids. Options: six frame translation using seqkit, pyrodigal-rv uses pyrodigal-meta with additional genetic codes, bbmap callgenes.sh (quick but less accurate for metagenomic data) (type: CHOICE; default: six_frame)
  • -db, --database: comma separated list of databases to search against (or all), or path to a custom database. options: NeoRdRp_v2.1, RdRp-scan, RVMT, TSA_2018, Pfam_RTs_RdRp, genomad, all, For custom path, either an .hmm file, a directory with .hmm files, or a folder with MSA files (which would be used to build an HMM DB) (type: TEXT; default: NeoRdRp_v2.1,genomad)
  • -t, --threads: Number of threads to use for searching (type: INTEGER; default: 1)
  • -g, --log-file: Absolute path to logfile (type: TEXT; default: ./marker_search_logfile.txt)
  • -ow, --overwrite: Do not overwrite the output directory if it already exists (type: BOOLEAN; default: False)
  • -td, -tempdir, --temp-dir: Path to temporary directory (type: TEXT; default: ./marker_search_tmp/)
  • --write-matched-regions, --no-write-matched-regions: Write matched query regions to FASTA (enabled by default; disable with --no-write-matched-regions) (type: BOOLEAN; default: True)
  • -mro, --matched-regions-output: Output FASTA path for matched regions (default: /marker_search_matched_regions.faa) (type: TEXT)
  • --include-aligned-region, --no-include-aligned-region: Include aligned query region sequence in marker_search_results.tsv (enabled by default) (type: BOOLEAN; default: True)
  • --include-alignment-string, --no-include-alignment-string: Include alignment identity string in marker_search_results.tsv (disabled by default) (type: BOOLEAN; default: False)