Finding stories in noise: Mitochondrial portraits from RAD data

Mitochondrial DNA (mtDNA) has formed the backbone of phylogeographic research for many years; however, recent trends focus on genome-wide analyses. One method proposed for calibrating inferences from noisy next-generation data, such as RAD sequencing, is to compare these results with analyses of mitochondrial sequences. Most researchers using this approach appear to be unaware that many single nucleotide polymorphisms (SNPs) identified from genome-wide sequence data are themselves mitochondrial, or assume that these are too few to bias analyses. Here, we demonstrate two methods for mining mitochondrial markers using RAD sequence data from three South African species of yellowfish, Labeobarbus. First, we use a rigorous SNP discovery pipeline using the program stacks, to identify variant sites in mtDNA, which we then combine into haplotypes. Second, we directly map sequence reads against a mitochondrial genome reference. This method allowed us to reconstruct up to 98% of the Labeobarbus mitogenome. We validated these mitogenome reconstructions through blast database searches and by comparison with cytochrome b gene sequences obtained through Sanger sequencing. Finally, we investigate the organismal consequences of these data including ancient genetic exchange and a recent translocation among populations of L. natalensis, as well as interspecific hybridization between L. aeneus and L. kimberleyensis.