This KBase annotation App (Annotate Microbial Genome) uses components from the RAST (Rapid Annotations using Subsystems Technology) toolkit [1,2,3] to annotate a prokaryotic genome or to update the annotations of a genome.

The release versions of the RASTtk component services used in this app are:

• kb_seed: tag 20200922
• kmer_annotation_figfam: tag 20200922
• genome_annotation: tag 20200922

The Annotate Microbial Genome App takes a Genome-typed object as input and allows users to annotate or re-annotate the genome. This will make the annotations consistent with other KBase genomes and prepare the genome for further analysis by other KBase Apps, especially the Metabolic Modeling Apps.

A Genome object can be imported or generated with one of the following annotation Apps or their multi-object versions:

• Annotate Microbial Assembly
• Annotate Microbial Genome
• Annotate Assembly and Re-annotate Genomes with Prokka

The Default Annotation Pipeline
Clicking "Run" will run the default pipeline. For a typical 2-5 MBp genome, the default annotation pipeline should take about 5 minutes. It is assumed that Genomes already have some annotation. As a result, the default behavior of this App is to use SEED to re-annotate just the protein-encoding genes. The default pipeline for this App consists of the following steps:

• DNA/RNA-based predictions
  1. Call rRNAs (default = off)
     Predict rRNAs in the genome. This is a custom BLAST-based tool for finding rRNAs.
  2. Call tRNAs with tRNAscan (default = off)
     Predict tRNAs in the genome with tRNAscan-SE [6].
  3. Call CRISPRs (default = off)
     This is a custom tool that uses a perl regular expression-based search to find CRISPR elements.
  4. Find prophage elements with phispy (default = off)
     This will use the phispy program to find prophage elements [13].
• Gene predictions
  1. Call protein-encoding genes with both Prodigal [11] and Glimmer3 [12] (default = off)
     These options will delete all existing genes in the genome object and replace them with the selected predictions.
  2. Call selenoproteins and pyrrolysylproteins [7] (default = off)
     These are custom BLAST-based tools.
• Repeats
  1. Call SEED large repeat regions (default = off)
     This is a BLASTn search within the genome for regions greater than 95% nucleotide similarity greater than or equal to 100bp in length.
  2. Find Streptococcus repeat regions [9, 10] (default = off)
     This is a command should only be implemented if the genus is Streptococcus.
• Add SEED Functions/Annotation to protein-encoding genes (k-mers needed for Metabolic Modeling)
  1. Annotate protein-encoding genes with k-mers (version 2; default = on)
     This is a set of signature k-mers (amino acid 8-mers) built from the annotations in the CoreSEED. The CoreSEED is a database of ~1,000 diverse microbial genomes and is currently the main focus of the RAST manual annotation efforts. Annotating using this k-mer set provides the user with our most stable and best estimate of the core gene functions.
  2. Annotate remaining hypothetical proteins with k-mers (version 1; default = on)
     This set of k-mers is built from the FigFam collection [8] in the PubSEED, which is the publically annotated version of the SEED database that consists of ~12,000 microbial genomes. The "classic" version of RAST on the RAST website (http://rast.nmpdr.org) uses the FigFam-based k-mers (hence the version 1 designation).
  3. Annotate remaining hypothetical proteins by protein similarity (default = on)
     We have several non-redundant databases for the most common genera. If the genus name of your organism matches one of these, a search will be performed against the remaining hypothetical proteins to attempt to find a function. The search uses a combination of BLAT [4] and BLAST [5].
• Other
  1. Perform a basic gene overlap removal (default = off)
     Using multiple gene calling algorithms can result in overlapping gene calls. This program is a custom tool that attempts to minimize overlaps and gaps to provide a set of calls that has a smaller number of gene calling errors. We do not recommend using overlap removal if you are attempting to annotate phage.
  2. Retain old annotations for hypotheticals (default = off)
     In instances where the pipeline fails to find an annotation for a gene, this will retain the original annotation from the input Genome-typed object.

Advanced Annotation Options
If you wish to customize the features in your annotation, click the "show advanced options" link. This will display the full set of available annotation options.

The Results

• The Objects section has a table of all the data objects that were created by the App. Click on the name of the data object to open a data viewer cell (below the currently selected cell).
• The Summary section give details about the coding and noncoding features that were created and the average protein length.

GUI Output
The GUI output currently consists of three tabs. The "Overview" tab provides basic information on the annotation job, the "Browse Features" tab allows the user to scroll through the features that were called, and the "Browse Contigs" tab provides information on the contigs in the genome. Users can sort on the various types of features. Note that some features will overlap (e.g., "prophage" and "CDS").

Additional Information
For more information on the steps of the default RAStk pipeline please refer to our publication on this (publication forthcoming). For more detailed tutorial information and to explore the additional functionality of RASTtk not currently available in the Narrative interface, please refer to http://tutorial.theseed.org.

Team members who developed & deployed algorithm in KBase: Thomas Brettin, James Davis, Terry Disz, Robert Edwards, Chris Henry, Gary Olsen, Robert Olson, Ross Overbeek, Bruce Parrello, Gordon Pusch, Roman Sutormin, and Fangfang Xia. For questions, please contact us.

The authors of RAST request that if you use the results of this annotation in your work, please cite the first three listed publications: