Stabilizing selection and mitochondrial heteroplasmy in the Canada lynx (<i>Lynx canadensis)</i>

Mitochondrial DNA is commonly used in population genetic studies to investigate spatial structure, intraspecific variation, and phylogenetic relationships. The control region is the most rapidly evolving and largest non-coding region, but its analysis can be complicated by heteroplasmic signals of genome duplication in many mammals, including felids. Here, we describe the presence of heteroplasmy in the control region of Canada lynx (<i>Lynx canadensis</i>) through intra-individual sequence variation. Our results demonstrate multiple haplotypes of varying length in each lynx, resulting from different copy numbers of the repetitive sequence RS-2 and suggest possible heteroplasmic single nucleotide polymorphisms in both repetitive sequences RS-2 and RS-3. Intra-individual variation was only observed in the repetitive sequences while inter-individual variation was detected in the flanking regions outside of the repetitive sequences, indicating that heteroplasmic mutations are restricted to these repeat regions. Although each lynx displayed multiple haplotypes of varying length, we found the most common variant contained three complete copies of the RS-2 repeat unit, suggesting copy number is regulated by stabilizing selection. While genome duplication offers potential for increased diversity, heteroplasmy may lead to a selective advantage or detriment in the face of mitochondrial function and disease, which could have significant implications for wildlife populations experiencing decline (e.g., bottlenecks) as a result of habitat modification or climate change.