HybPiper: Extracting coding sequence and introns for phylogenetics from high‐throughput sequencing reads using target enrichment

Abstract

Premise of the study: Using sequence data generated via target enrichment for phylogenetics requires reassembly of highthroughput sequence reads into loci, presenting a number of bioinformatics challenges. We developed HybPiper as a userfriendly platform for assembly of gene regions, extraction of exon and intron sequences, and identification of paralogous gene copies. We test HybPiper using baits designed to target 333 phylogenetic markers and 125 genes of functional significance in Artocarpus (Moraceae). Methods and Results: HybPiper implements parallel execution of sequence assembly in three phases: read mapping, contig assembly, and target sequence extraction. The pipeline was able to recover nearly complete gene sequences for all genes in 22 species of Artocarpus. HybPiper also recovered more than 500 bp of nontargeted intron sequence in over half of the phylogenetic markers and identified paralogous gene copies in Artocarpus. Conclusions: HybPiper was designed for Linux and Mac OS X and is freely available at https://github.com/mossmatters/HybPiper.