#fastq #bam #complex #typing #fasta #tuberculosis #mtbc

app fastlin

an ultra-fast program for MTBC lineage typing

4 releases (2 breaking)

0.4.1 Sep 22, 2024
0.4.0 Sep 22, 2024
0.2.1 Jul 17, 2023
0.2.0 Jul 16, 2023
0.1.0 Jun 26, 2023

#61 in Biology

Download history 225/week @ 2024-09-16 115/week @ 2024-09-23 22/week @ 2024-09-30 8/week @ 2024-10-07

370 downloads per month

MIT/Apache

38KB
513 lines

Anaconda-Server Badge Crates.io GitHub release (latest SemVer)

fastlin

Overview

Fastlin is an ultra-fast program to perform lineage typing of Mycobacterium tuberculosis complex (MTBC) FASTQ read data, BAM files and FASTA assemblies. Using the split-kmer approach, it can accuratly predict MTBC lineages and strain mixtures in seconds.

Reference: fastlin: an ultra-fast program for Mycobacterium tuberculosis complex lineage typing.

Main updates since publication:

  • 0.2.3 : FASTA files as input (also using seq_io)
  • 0.3.0 : multi-threading available
  • 0.4.0 : BAM files as input (using rust-htslib)

There is no planned updates at the moment. Please open an issue if you have any suggestion or request.

Installation

To install fastlin via cargo (you must have the rust toolchain installed):

cargo install fastlin

Or you can download the latest release from this repository and compile it using cargo:

cargo install --path directory_release

Alternatively, you can install precompiled binaries using Conda:

conda install -c bioconda fastlin

You will also need a barcode file (see Input files below).

Running fastlin

The default command line is:

fastlin -d your_directory -b your_barcodes.txt

If your dataset consists of FASTQ files that are not BAM-derived, then you can apply a maximum kmer coverage threshold to reduce runtimes:

fastlin -d your_directory -b your_barcodes.txt -x 80

Input files

Fastlin takes as input the path of the directory containing FASTQ, BAM and/or FASTA files. The directory can contain a mix of FASTA geome assemblies, BAM alignment files, paired-end and single-end FASTQ files. FASTQ and FASTA files should be gzipped, with the following extensions:

  • .fastq.gz or .fq.gz for FASTQ read data. The names of paired-end files should be in the form name_1.fq.gz and name_2.fq.gz (or equivalent with fastq.gz)
  • .bam, or .BAM for BAM files. Here the maximum kmer coverage is ignored since BAM files can be sorted.
  • .fas.gz, .fasta.gz or .fna.gz for FASTA genome assemblies. Here, minimum occurence is set to 1 and the maximum kmer coverage is ignored.

Please note that BAM files are analyzed in the same way as FASTQ files, by scanning reads without considering quality scores. You may find faster scripts or programs that focus solely on specific genomic positions.

The MTBC barcode file can be downloaded from https://www.github.com/rderelle/barcodes-fastlin. Alternatively, you can build and test your own kmer barcodes using the Python scripts available in that directory.

Manual

A full description of fastlin parameters can be found here.

Output file

Fastlin output consists of a tab-delimited file with the following fields:

  • sample: sample name
  • data type: 'assembly', 'BAM', 'single' (reads) or 'paired' (-end reads)
  • k_cov: theoretical kmer coverage of the fastq files(s) based on the number of extracted kmers
  • mixture: pure ('no') or mixed ('yes') sample
  • lineages: detected lineages (median kmer occurences within paratheses)
  • log_barcodes: kmer barcodes passing the minimum occurence threshold, indicated by their kmer occurence and grouped by lineages

Here is a simple example:

#sample    data type    k_cov    mixture    lineages    log_barcodes    log_errors
ERRxxxxx    paired    118    no    2 (45)    2 (42, 48, 39, 43, 54, 47, 45), 4.1 (4)

The sample ERRxxxxx contains a single strain belonging to lineage 2. This typing is supported by 7 kmer barcodes, with a median number of occurences of 45. Since the abundance of the strain is far below the theoretical kmer coverage (equal here to 118), we can conclude that the sample is likely to contain high level of contaminations or sequencing errors.

Multi-threading

By default, fastlin runs on 1 thread. The number of threads can be increased using the '-t' parameter, which will split the sample set among all threads (for a single sample, increasing the number of threads will have no impact on runtime).

Here are some examples of runtimes (in seconds) using real-world Mtb genomic data on a M2 Macbook Air:

data 1 thread 4 threads
12 paired FASTQ 66.9 19.3
4 BAM files 26.9 9.4
190 genomes FASTA 6.7 1.8

Error handling

When fastlin cannot read a fastq file (e.g., faulty record within the fastq file, corrupt gzip file), it stops scanning it, re-initialises all values to 0 and reports the error message in the last column of the output file. Here is an example of output with 3 different errors:

#sample    data type    k_cov    mixture    lineages    log_barcodes    log_errors
dummy1   single   0   no       Error in file "reads/dummy1.fastq.gz": FASTQ parse error: sequence length is 150, but quality length is 50 (record 'ERR551806.5' at line 17).
dummy2   single   0   no       Error in file "reads/dummy2.fastq.gz": invalid gzip header
dummy3   single   0   no       Error in file "reads/dummy3.fastq.gz": corrupt deflate stream

Dependencies

~17–27MB
~433K SLoC