BUSCO: Assessing Genomic Data Quality and Beyond

Abstract

Evaluation of the quality of genomic “data products” such as genome assemblies or gene sets is of critical importance in order to recognize possible issues and correct them during the generation of new data. It is equally essential to guide subsequent or comparative analyses with existing data, as the correct interpretation of the results necessarily requires knowledge about the quality level and reliability of the inputs. Using datasets of near universal single‐copy orthologs derived from OrthoDB, BUSCO can estimate the completeness and redundancy of genomic data by providing biologically meaningful metrics based on expected gene content. These can complement technical metrics such as contiguity measures (e.g., number of contigs/scaffolds, and N50 values). Here, we describe the use of the BUSCO tool suite to assess different data types that can range from genome assemblies of single isolates and assembled transcriptomes and annotated gene sets to metagenome‐assembled genomes where the taxonomic origin of the species is unknown. BUSCO is the only tool capable of assessing all these types of sequences from both eukaryotic and prokaryotic species. The protocols detail the various BUSCO running modes and the novel workflows introduced in versions 4 and 5, including the batch analysis on multiple inputs, the auto‐lineage workflow to run assessments without specifying a dataset, and a workflow for the evaluation of (large) eukaryotic genomes. The protocols further cover the BUSCO setup, guidelines to interpret the results, and BUSCO “plugin” workflows for performing common operations in genomics using BUSCO results, such as building phylogenomic trees and visualizing syntenies. © 2021 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Assessing an input sequence with a BUSCO dataset specified manually Basic Protocol 2: Assessing an input sequence with a dataset automatically selected by BUSCO Basic Protocol 3: Assessing multiple inputs Alternate Protocol: Decreasing analysis runtime when assessing a large number of small genomes with BUSCO auto‐lineage workflow and Snakemake Support Protocol 1: BUSCO setup Support Protocol 2: Visualizing BUSCO results Support Protocol 3: Building phylogenomic trees

Molecular Characterization of Candida auris Isolates at a Major Tertiary Care Center in Lebanon

“…The ribosomal RNA (rRNA) identification was done by downloading C. auris rRNA sequences from the National Center for Biotechnology Information (NCBI) and then by using BLAST + v2.11.0 ( Camacho et al, 2009 ). Assemblies and annotations were assessed using BUSCO V5.2.2 ( Manni et al, 2021 ).…”

The mosaic oat genome gives insights into a uniquely healthy cereal crop

“…1a and Supplementary Table 1 ), with a BUSCO (v5.1.2; ref. 6 ) score of 98.7% (Extended Data Fig. 2a ), following the short-read strategy used for wheat 7 , barley 8 and rye 9 .…”

Odorant receptor orthologues in conifer-feeding beetles display conserved responses to ecologically relevant odors

“…The average transcript length was 824 bp with an N50 of 1,570 bp. The completeness of the sex-combined assembly was firstly assessed using the Benchmarking Universal Single-Copy Orthologs (BUSCO v5.2.2; https://busco.ezlab.org/) tool performed against the Insecta odb10 dataset, including 1,367 reference genes (Manni, et al 2021). This analysis revealed 96.6% complete (C) BUSCOs, of which 59.8% were present as single copy genes (S) and 40.2% as duplicated genes (D).…”