BLASTmap: A shiny-based application to visualize BLAST results as interactive heat maps and a tool to design gene-specific baits for bespoke target enrichment sequencing

Abstract

Numerous genes that determine the outcome of plant–pathogen interactions are currently being discovered and include, for example, immune receptors, susceptibility factors and pathogen effectors and their host targets. Target enrichment sequencing provides a means to preferentially resequence these genes of interest without the need to first generate a genotype-specific genome assembly. The Basic Local Alignment Search Tool (BLAST), in combination with the here developed BLASTmap, can be used to design probes that specifically target such gene(s), either by using the target species or the closest related genus as a reference. BLAST is a ubiquitous tool in biological sequence analysis and a multitude of programs are available for the visualization of BLAST alignments. However, there are currently no dedicated programs for visual comparison of large-scale BLAST output attributes such as bit score. The need to quickly and efficiently compare many thousands of BLAST results led to the development of BLASTmap, an interactive web application created using the Shiny R package, customized for clustering and viewing BLAST results as an interactive heat map. Here we show an example of how BLASTmap was successfully applied to analyze custom DNA/RNA probe sequences and to visually determine that four probes are sufficient for the specific yet inclusive enrichment of the potato R2 disease resistance gene family.