Biodiversity genomics for species at risk: patterns of DNA sequence variation within and among complete mitochondrial genomes of three species of wolffish (Anarhichas spp.)

Abstract

The first marine fish species to be listed under the Canadian Species At Risk Act as Threatened with extinction are the spotted wolffish ( Anarhichas minor Olafsen, 1774) and the broadhead wolffish ( Anarhichas denticulatus Krøyer, 1844); a third species, the striped wolffish ( Anarhichas lupus L., 1758), is listed as a species of special concern. As part of the recovery plan for wolffish, we determined the complete mitochondrial DNA (mtDNA) genome sequences of all three species to identify the most variable gene regions for population analysis. Anarhichas genomes comprise either 16 519 or 16 520 base pairs (bp), among which there are 449 single nucleotide polymorphisms (SNPs). The most variable protein-coding loci are ND4, CYTB, and ND2, with 4.40, 4.22, and 4.19 SNPs / 100 bp, respectively. Comparisons of rates of synonymous and nonsynonymous nucleotide substitutions indicate no evidence of selection. The control region, characterized in many species as hypervariable, was less variable than 9 of 13 protein-coding loci (2.45 SNPs / 100 bp). Phylogenetic analysis shows that A. lupus and A. minor are more closely related to each other than either is to A. denticulatus. Amplification and sequence analysis of a contiguous block of 6392 bp that spans the ND4, ND5, ND6, and CYTB loci is an efficient strategy for evaluating patterns of intraspecific mtDNA variability.