HLA typing (Long-read RNA)

HLA typing in long-read RNA sequencing FASTQ files

tools

The following tools run by default (methods: "all"):

hlaminer, specimmune

usage

nexus run --nf-workflow hla_typing_long-read-rna.nf -params-file params.yaml

Note

Nextflow config files are available here. Use the config file that matches your installed nexus version (e.g. nexus_v0.2.0_nextflow_slurm.config).

parameters

Download params.yaml

# =============================================================================
# params.yaml — hla_typing_long-read-rna
#
# Usage:
#   nextflow run hla_typing_long-read-rna.nf -params-file params.yaml
#
# Fill in the required fields below. Optional fields can be left at their
# defaults or removed entirely.
# =============================================================================


# -----------------------------------------------------------------------------
# Required
# -----------------------------------------------------------------------------

# TSV file with columns:
#   sample_id
#   fastq_file             long-read RNA FASTQ (PacBio or Nanopore), gzipped
samples_tsv_file: ""

# Directory to which output files will be copied
output_dir: ""


# -----------------------------------------------------------------------------
# Optional — general
# -----------------------------------------------------------------------------

# Methods to run. Comma-separated list or 'all'.
# Allowed values:
#   all, hlaminer, specimmune
methods: "all"


# -----------------------------------------------------------------------------
# HLAminer (long-read RNA-seq)
# Optional arguments.
# -----------------------------------------------------------------------------
# Pipeline:
#   minimap2 (against HLAminer's HLA CDS reference) | HLAminer.pl
hlaminer:
  # Extra CLI arguments passed directly to minimap2 for the upstream
  # alignment step. Use -ax map-hifi for PacBio HiFi, -ax map-ont for ONT.
  minimap2_extra_args: "-ax map-hifi --secondary=no --MD"

  # Extra CLI arguments passed directly to HLAminer.pl. Common flags:
  #   -e <expect value threshold>
  #   -s <minimum alignment score>
  #   -q <minimum mapping quality>
  #   -i <minimum percent identity>
  hlaminer_extra_args: "-e 1 -s 500 -q 1 -i 1"


# -----------------------------------------------------------------------------
# SpecImmune
# Optional arguments.
# -----------------------------------------------------------------------------
specimmune:
  # Extra CLI arguments passed directly to SpecImmune (scripts/main.py).
  #   -i   typing target:  HLA | KIR | CYP | IG_TR | extend
  #   -y   read type:      nanopore | pacbio | pacbio-hifi
  #   --RNA_type           traditional | 2D | Direct | SIRV  (selects RNA mode)
  extra_args: "-i HLA -y pacbio-hifi --RNA_type traditional"