Mitochondrial genomes of the bird genus Piranga: rates of sequence evolution, and discordance between mitochondrial and nuclear markers

We report the characteristics of the mitochondrial genomes of 22 individuals in the bird genus Piranga, including all currently recognized species in the genus (n = 11). Elements follow the standard avian mitogenome series, including two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, 13 protein coding genes, and the mitochondrial control region. Excluding tRNA sequences, sequence divergence rate was lowest in rRNA genes and highest in genes encoding NADH (specifically ND1, ND2, ND3) and the control region. Gene trees assembled from 16 elements (non-tRNAs) varied greatly in topological concordance compared to the recognized species tree (based on thousands of nuclear loci), with no one gene tree precisely recovering all relationships in the genus. We also investigated patterns of concordance between the mitogenome tree and the nuclear species tree and found some discrepancies. Across non-tRNA gene trees (n = 16), the species tree topology was recovered by as few as three elements at a particular node and complete concordance (i.e. 16/16 gene trees matched the species tree topology) was recovered at only one node. We found mitochondrial gene regions that are often used in vertebrate systematics (e.g. CytB, ND2) recovered nearly the exact same topology as the nuclear species tree topology.