# Complete mitochondrial genomes of three vulnerable cave bat species and their phylogenetic relationships within the order Chiroptera

## TL;DR

This study sequenced the complete mitochondrial genomes of three vulnerable Brazilian bat species and analyzed their evolutionary relationships within the order Chiroptera.

## Contribution

The first complete mitochondrial genomes of Furipterus horrens, Lonchorhina aurita, and Natalus macrourus are reported, along with their phylogenetic placement.

## Key findings

- The mitogenomes of the three species are similar in length and gene arrangement to other bat species.
- Phylogenetic analysis placed F. horrens as sister to Noctilio leporinus and N. macrourus as sister to Molossidae and Vespertilionidae.
- All protein-coding genes are under purifying selection, with atp8 showing the most relaxed selection.

## Abstract

The IUCN Red List of Threatened Species contains 175 Brazilian bat species that are threatened by extinction in some degree. From this perspective, it is essential to expand the knowledge about the genetic diversity of vulnerable bats. Genomic sequencing can be useful to generate robust and informative genetic references, increasing resolution when analyzing relationships among populations, species, or higher taxonomic levels. In this study, we sequenced and characterized in detail the first complete mitochondrial genomes of Furipterus horrens, Lonchorhina aurita, and Natalus macrourus, and investigated their phylogenetic position based on amino acid sequences of protein-coding genes (PCGs). The mitogenomes of these species are 16,516, 16,697, and 16,668 bp in length, respectively, and each comprises 13 PCGs, 22 tRNA genes, two rRNA genes, and a putative control region (CR). In the three species, genes were arranged similarly to all other previously described bat mitogenomes, and nucleotide composition was also consistent with the reported range. The length and arrangement of rrnS and rrnL were also consistent with those of other bat species, showing a positive AT-skew and a negative GC-skew. Except for trnS1, for which we did not observe the DHU arm, all other tRNAs showed the cloverleaf secondary structure in the three species. In addition, the mitogenomes showed minor differences in start and stop codons, and in all PCGs, codons ending in adenine were more common compared to those ending in guanine. We found that PCGs of the three species use multiple codons to encode each amino acid, following the previously documented pattern. Furthermore, all PCGs are under purifying selection, with atp8 experiencing the most relaxed purifying selection. Considering the phylogenetic reconstruction, F. horrens was recovered as sister to Noctilio leporinus, L. aurita and Tonatia bidens shared a node within Phyllostomidae, and N. macrourus appeared as sister to Molossidae and Vespertilionidae.