Complete Mitogenomes of Xinjiang Hares and Their Selective Pressure Considerations

Comparative analysis based on the mitogenomes of hares in Xinjiang, China, is limited. In this study, the complete mitochondrial genomes of seven hare samples including four hare species and their hybrids from different environments were sequenced, assembled, and annotated. Subsequently, we performed base content and bias analysis, tRNA analysis, phylogenetic analysis, and amino acid sequence analysis of the annotated genes to understand their characteristics and phylogenetic relationship. Their mitogenomes are circular molecules (from 16,691 to 17,598 bp) containing 13 protein-coding genes, two rRNA genes, 22 tRNA genes, and a control region, which are similar with other Lepus spp. worldwide. The relative synonymous codon usage analysis revealed that the adaptation of Lepus yarkandensis to its unique arid and hot environment might be associated with synthesizing amino acids like alanine, leucine, serine, arginine, and isoleucine and the terminator caused by the different usage of codons. Further, we utilized the MEME model and identified two positive selection genes (ND4, ND5) in Lepus tibetanus pamirensis and one (ND5) in L. yarkandensis that might be important to their adaptation to the plateau and dry and hot basin environments, respectively. Meanwhile, Lepus tolai lehmanni and Lepus timidus may have evolved different adaptive mechanisms for the same cold environment. This study explored the evolutionary dynamics of Xinjiang hares' mitochondrial genomes, providing significant support for future research into their adaptation mechanisms in extreme environments.