▶️DRAGEN Microbial Amplicon App Documentation

Overview

DRAGEN Microbial Amplicon is a software application designed to analyze sequencing data from amplicon library preps (both DNA and RNA) on microbiological samples, with an emphasis on viruses. Illumina sequencing reads are processed to generate consensus sequences that represent a best estimate of the population of viral sequences in each sample. Where appropriate, these consensus sequences are further analyzed by the phylogenetic analysis tools Nextclade and/or Pangolin to provide an identification of the clade or lineage of each sequence.

Input

Data can be provided in one of the following ways:

• Samples / biosamples with FASTQ datasets (see details in library preparation documents)

• A project containing one or more samples / biosamples with FASTQ datasets

  • All samples / biosamples in the selected project will be analyzed

Supported amplicon primer schemes

• Chikungunya

  • Grubaugh Lab

  • Illumina

• Dengue

  • Serotype 1 - Illumina

  • All serotypes - DengueSeq from Grubaugh Lab

• Influenza A/B - Universal

• Mpox

  • Pan-clade - ARTIC

  • Clade I - Illumina

  • Clade II - Grubaugh Lab

• RSV

  • CDC

  • WCCRRI

• SARS-CoV-2 - ARTIC

  • v5.4.2

  • v5.3.2, v4.1, v4, v3

• Zika - Grubaugh Lab

Custom genome and primer sets

Users can upload custom files to provide user-defined reference genome set and primer definitions. Multiplexed amplicon panels targeting multiple organisms in the same reaction are supported.

🧬Custom reference

Pipeline steps

1. Trim and filter reads using Trimmomatic

2. Remove off-target reads using DRAGEN v4.3.6 kmer classifier (for custom reference, remove human reads using a modified version of the SRA Human Read Scrubber tool v2.2.1)

3. For organisms with one default reference genome, skip this step. For organisms with multiple candidates, trim primer sequences in reads using Trimmomatic, perform assembly using MEGAHIT, cluster contigs using CD-HIT-EST, map contigs to candidate reference genomes using minimap2, then select reference genomes based on the mapping

4. Align reads to the default reference genome or selected reference genomes using DRAGEN v4.3.6

5. Trim primer sequences in aligned reads based on coordinates

6. Filter out samples with insufficient amplicon coverage

7. Call sequence variants from the alignments using DRAGEN Somatic v4.3.6 and apply them to the corresponding reference genomes to create consensus sequences

8. If applicable, run Nextclade/Pangolin on the consensus sequences

💠Pipeline Logic

Output

• Consensus sequences representing a best estimate of targeted sequences

• Tables and plots reporting read counts, coverage, and Nextclade/Pangolin results

📂Output files

Currently supported platforms

• BaseSpace Sequence Hub

Important Notes

• The sequences are labeled according to the best match in the reference database, which is not exhaustive and the labels should not be taken as definitive for strain-typing. If strain typing is needed, the built-in Nextclade and/or Pangolin tools can be used for supported organisms. Alternatively, a BLAST or similar search of nucleotide databases may provide a more detailed match.

• Because of sequence homology, it is possible that organisms with very few reads will result in the generation of a sequence not present (false positive). Although the de novo assembly step of this software largely mitigates such instances, sequences with very low horizontal coverage (< 5%) should be treated with caution.