Complete mitochondrial genome sequences of three bats species and whole genome mitochondrial analyses reveal patterns of codon bias and lend support to a basal split in Chiroptera

Comparative mitochondrial genomics of endemic Mexican vesper yellow bats genus Rhogeessa (Chiroptera: Vespertilionidae) and insights into internal relationships in the family Vespertilionidae

In the species-rich family Vespertilionidae, vesper yellow bats in the genus Rhogeessa include eleven species, three of them endemic to Mexico. These insectivorous bats provide important ecosystem services, including pest control. Even though some aspects of their biology are well- known, only a few genomic resources are available for these species, which limits our understanding of their biology. In this study, we assembled and annotated the mitochondrial genome of four species: R. aenea, R. genowaysi, R. mira, and R. parvula. We generated a phylomitogenomic hypothesis based on translated protein-coding genes for a total of 52 species in the family Vespertilionidae and examined the phylogenetic position of the genus Rhogeessa and species within the family. The AT-rich mitogenomes of R. aenea, R. genowaysi, R. mira, and R. parvula are 16,763, 16,781, 16,807, and 16,794 pb in length, respectively. Each studied mitogenome encodes 13 Protein Coding Genes (PCGs), 22 transfer RNA genes, and 2 rRNA genes, and contains a putative control region (CR). All tRNAs exhibit a 'cloverleaf' secondary structure, except tRNA-Serine-1 that lacked the DHU arm in all studied mitogenomes. Selective pressure analyses indicated that all protein-coding genes are exposed to purifying selection. The phylomitogenomic analysis supported the monophyletic status of the family Vespertilionidae, confirmed the placement of Rhogeessa within the tribe Antrozoini, and clarified phylogenetic relationships within and among subfamilies and tribes in this family. Our results indicate that phylomitogenomics are useful to explore the evolutionary history of vesper bats. The assembly and comprehensive analysis of mitochondrial genomes offer the potential to generate molecular references and resources beneficial for genetic analyses aimed at understanding the ecology and evolution of these remarkable bats.

Analysis of the complete mitochondrial genome sequence of Hipposideros pratti

In order to explore the characteristics of the mitochondrial genome sequence of Pratt's leaf-nosed bat (Hipposideros pratti Thomas 1891) and understand their phylogenetic status in Chiroptera, this study determined the mitochondrial genome sequences of H. pratti from five regions in China using high-throughput sequencing technology, sequence assembly, and genome annotation. The results showed that these sequences contained 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and 1 non-coding region, all exhibiting a significant AT bias. Based on the phylogenetic tree constructed using 13 protein-coding genes from 15 Chiroptera species, the study found that H. pratti from the five regions clustered together, and then clustered with H. lylei into a single clade. Meanwhile, H. pratti from Jiangxi, Fujian, and Guangdong regions of China showed closer genetic relationships, while H. pratti from Yunnan and Henan regions of China exhibited closer genetic relationships. This study not only supplemented the mitochondrial genome database of H. pratti but also laid a foundation for genetic variation, molecular classification, and evolutionary studies of H. pratti.

Characterization of the mitochondrial genomes of the Mexican endemic bats Corynorhinus mexicanus and Corynorhinus leonpaniaguae (Chiroptera: Vespertilionidae)

BACKGROUND

The genus Corynorhinus is composed of four recognized species: C. rafinesquii, C. townsendii, C. mexicanus, and C. leonpaniaguae, the latter two being endemic to Mexico. According to the IUCN, C. mexicanus is considered "Near Threatened", as its populations are dwindling and habitats are affected by anthropogenic disturbance. Corynorhinus leonpaniaguae has not been assigned to an IUCN Red List risk category due to its recent description.

METHODS AND RESULTS

In this study, the mitochondrial genomes of C. mexicanus and C. leonpaniaguae were assembled and characterized in detail. The mitochondrial genomes (mtDNA) of C. mexicanus and C. leonpaniaguae have lengths of 16,470 and 16,581 bp respectively, with a predominant nucleotide usage of adenine (31.670% and 31.729%, respectively) and thymine (26.15% and 26.18%, respectively). The mtDNA of C. mexicanus and C. leonpaniaguae is composed of 37 coding and non-coding elements: 22 transfer RNAs (tRNA), 13 protein-coding genes (PCGs), two ribosomal RNAs and a non-coding region, the control region, which has a length of 933 bp and 1,149 bp, respectively. All tRNAs exhibited a cloverleaf secondary structure, with the exception of trn-Ser1 which showed a deletion of the dihydrouridine arm in the two species. All PCGs are subjected to purifying selection, with atp8 being the gene showing the highest Ka/Ks value.

CONCLUSIONS

These are the first whole mitogenomic resources developed for C. mexicanus and C. leonpaniaguae and enhance our knowledge of the ecology of these species and aid in their conservation.

Mitochondrial Genomic Evidence of Selective Constraints in Small-Bodied Terrestrial Cetartiodactyla

Body size may drive the molecular evolution of mitochondrial genes in response to changes in energy requirements across species of different sizes. In this study, we perform selection pressure analysis and phylogenetic independent contrasts (PIC) to investigate the association between molecular evolution of mitochondrial genome protein-coding genes (mtDNA PCGs) and body size in terrestrial Cetartiodactyla. Employing selection pressure analysis, we observe that the average non-synonymous/synonymous substitution rate ratio (ω) of mtDNA PCGs is significantly reduced in small-bodied species relative to their medium and large counterparts. PIC analysis further confirms that ω values are positively correlated with body size (R2 = 0.162, p = 0.0016). Our results suggest that mtDNA PCGs of small-bodied species experience much stronger purifying selection as they need to maintain a heightened metabolic rate. On the other hand, larger-bodied species may face less stringent selective pressures on their mtDNA PCGs, potentially due to reduced relative energy expenditure per unit mass. Furthermore, we identify several genes that undergo positive selection, possibly linked to species adaptation to specific environments. Therefore, despite purifying selection being the predominant force in the evolution of mtDNA PCGs, positive selection can also occur during the process of adaptive evolution.

Reference genome of Townsend's big-eared bat, Corynorhinus townsendii

Townsend's big-eared bat, Corynorhinus townsendii, is a cave- and mine-roosting species found largely in western North America. Considered a species of conservation concern throughout much of its range, protection efforts would greatly benefit from understanding patterns of population structure, genetic diversity, and local adaptation. To facilitate such research, we present the first de novo genome assembly of C. townsendii as part of the California Conservation Genomics Project (CCGP). Pacific Biosciences HiFi long reads and Omni-C chromatin-proximity sequencing technologies were used to produce a de novo genome assembly, consistent with the standard CCGP reference genome protocol. This assembly comprises 391 scaffolds spanning 2.1 Gb, represented by a scaffold N50 of 174.6 Mb, a contig N50 of 23.4 Mb, and a benchmarking universal single-copy ortholog (BUSCO) completeness score of 96.6%. This high-quality genome will be a key tool for informed conservation and management of this vulnerable species in California and across its range.

Comparative mitochondrial genomics of American nectar-feeding long-nosed bats Leptonycteris spp. with insights into the phylogeny of the family Phyllostomidae

Among leaf-nosed bats (family Phyllostomidae), the genus Leptonycteris (subfamily Glossophaginae), contains three migratory and obligate nectar-dwelling species of great ecological and economic importance; the greater long-nosed bat L. nivalis, the lesser long-nosed bat L. yerbabuenae, and the southern long-nosed bat L. curasoae. According to the IUCN, the three species are categorized as 'vulnerable', 'endangered', and 'near threatened', respectively. In this study, we assembled and characterized in detail the mitochondrial genome of Leptonycteris spp. and examined the phylogenetic position of this genus in the family Phyllostomidae based on protein coding genes (PCGs). The mitogenomes of L. nivalis, L. curasoae, and L. yerbabuenae are 16,708, 16,758, and 16,729 bp in length and each encode 13 PCGs, 22 transfer RNA genes, 2 rRNA genes, and a putative control region (CR). Mitochondrial gene order is identical to that reported before for the family Phyllostomidae. All tRNAs exhibit a 'cloverleaf' secondary structure, except tRNA-Serine-1 that is missing the DHU arm in the three species. All PCGs are exposed to purifying selection with atp8 experiencing the most relaxed purifying selection as the ω ratio was higher for this gene compared to the other PCGs in each species. The CR of each species contains three functional domains: extended termination associated sequence (ETAS), Central, and conserved sequence block domain (CSB). A phylomitogenomic analysis revealed that Leptonycteris is monophyletic and most closely related to the genus Glossophaga. The analysis also supported the monophyly of the family Glossophaginae in the speciose family Phyllostomidae. The mitochondria characterization of these species provides relevant information to develop molecular markers for conservation purposes.

The complete mitochondrial genomes of the Leuciscus baicalensis and Rutilus rutilus: a detailed genomic comparison among closely related species of the Leuciscinae subfamily

Cyprinidae is the largest family in the order of freshwater fish Cypriniformes. Increased subfamily members of Cyprinidae have been suggested to be re-classified for decades. In this study, we sequenced the mitochondrial genomes (mitogenomes) of Leuciscus baicalensis and Rutilus rutilus collected from northwest China and compared with other closely related species to determine their associated family or subfamily. We used Illumina NovaSeq to sequence the entire mitochondrial genomes of Leuciscus baicalensis and Rutilus rutilus and characterized the mitogenomes by the gene structure, gene order, and the secondary structures of the 22 tRNA genes. We compared mitogenome features of Leuciscinae with other subfamilies in Cyprinidae. We used the analytic Bayesian Information and Maximum Likelihood methods to determine phylogenetic trees of 13 PCGs. The mitogenomes of Leuciscus baicalensis and Rutilus rutilus were 16,607 bp and 16,606 bp, respectively. Organization and location of these genes were consistent with already studied Leuciscinae fishes. Synonymous codon usage was conservative in Leuciscinae as compared with other subfamilies in Cyprinidae. Phylogenetic analysis indicated that Leuciscinae was a monophyletic group, and genus Leuciscus was a paraphyletic group. Our approach, for the first time, of studying comparative mitochondrial genomics and phylogenetics together provided a supportive platform to the analysis of population genetics and phylogeny for Leuciscinae. Our results indicated a promising potential of comparative mitochondrial genomics in the manifestation of phylogenetic relationships between fishes, leading us to a suggestion that mitogenomes should be routinely considered in clarifying phylogenetics of family and subfamily members of fish.

The mitochondrial genomes of big-eared bats, Macrotus waterhousii and Macrotus californicus (Chiroptera: Phyllostomidae: Macrotinae)

In the species-rich family Phyllostomidae, the genus Macrotus ('big eared' bats) contains only two species; Macrotus waterhousii, distributed in western, central, and southern Mexico, Guatemala and some Caribbean Islands, and Macrotus californicus, distributed in the southwestern USA, and in the Baja California peninsula and the state of Sonora in Mexico. In this study, we sequenced and assembled the mitochondrial genome of Macrotus waterhousii and characterized in detail this genome and that of the congeneric M. californicus. Then, we examined the phylogenetic position of Macrotus in the family Phyllostomidae based on protein coding genes (PCGs). The AT-rich mitochondrial genomes of M. waterhousii and M. californicus are 16,792 and 16,691 bp long, respectively, and each encode 13 PCGs, 22 tRNA genes, 2 rRNA genes, and a putative non-coding control region 1,336 and 1,232 bp long, respectively. Mitochondrial synteny in Macrotus is identical to that reported before for all other cofamilial species. In the two studied species, all tRNAs exhibit a 'typical' cloverleaf secondary structure with the exception of trnS1, which lacks the D arm. A selective pressure analysis demonstrated that all PCGs are under purifying selection. The CR of the two species feature three domains previously reported in other mammals, including bats: extended terminal associated sequences (ETAS), central (CD), and conserved sequence block (CSB). A phylogenetic analysis based on the 13 mitochondrial PCGs demonstrated that Macrotus is monophyletic and the subfamily Macrotinae is a sister group of all remaining phyllostomids in our analysis, except Micronycterinae. The assembly and detailed analysis of these mitochondrial genomes represents a step further to continue improving the understanding of phylogenetic relationships within the species-rich family Phyllostomidae.

The metazoan landscape of mitochondrial DNA gene order and content is shaped by selection and affects mitochondrial transcription

Mitochondrial DNA (mtDNA) harbors essential genes in most metazoans, yet the regulatory impact of the multiple evolutionary mtDNA rearrangements has been overlooked. Here, by analyzing mtDNAs from ~8000 metazoans we found high gene content conservation (especially of protein and rRNA genes), and codon preferences for mtDNA-encoded tRNAs across most metazoans. In contrast, mtDNA gene order (MGO) was selectively constrained within but not between phyla, yet certain gene stretches (ATP8-ATP6, ND4-ND4L) were highly conserved across metazoans. Since certain metazoans with different MGOs diverge in mtDNA transcription, we hypothesized that evolutionary mtDNA rearrangements affected mtDNA transcriptional patterns. As a first step to test this hypothesis, we analyzed available RNA-seq data from 53 metazoans. Since polycistron mtDNA transcripts constitute a small fraction of the steady-state RNA, we enriched for polycistronic boundaries by calculating RNA-seq read densities across junctions between gene couples encoded either by the same strand (SSJ) or by different strands (DSJ). We found that organisms whose mtDNA is organized in alternating reverse-strand/forward-strand gene blocks (mostly arthropods), displayed significantly reduced DSJ read counts, in contrast to organisms whose mtDNA genes are preferentially encoded by one strand (all chordates). Our findings suggest that mtDNA rearrangements are selectively constrained and likely impact mtDNA regulation.

Revised phylogeny from complete mitochondrial genomes of phyllostomid bats resolves subfamilial classification

Classically, molecular phylogenetic trees of Phyllostomidae have been inferred using a combination of a few mitochondrial and nuclear markers. However, there is still uncertainty in the relationships, especially among deep clades within the family. In this study, we provide newly sequenced complete mitochondrial genomes from 26 bat species, including genomes of 23 species reported here for the first time. By carefully analysing these genomes using maximum likelihood and Bayesian methods and different ingroup and outgroup samples, partition schemes and data types, we investigated the robustness and sensitivity of our phylogenetic results. The optimal topologies were those inferred from the complete data matrix of nucleotides, with complex and highly parameterized substitution models and partition schemes. Our results show a statistically robust picture of the evolutionary relationships between phyllostomid subfamilies and clarify hitherto uncertain relationships of Lonchorhininae and Macrotinae.

Complete mitochondrial genome of the greater round-eared bat, Tonatia bidens (Chiroptera: Phyllostomidae) from Brazil and phylogenetic relationships within the family

The greater round-eared bat, Tonatia bidens, is a locally rare species belonging to the highly diverse family Phyllostomidae. In this study, the complete mitogenome of T. bidens was sequenced using optimized protocols of DNA extraction from fixed cells originally prepared for cytogenetic studies. Here we present the complete mitogenome and place our results in a phylogenetic context with other data generated for the family Phyllostomidae. The circular genome had 16,717 bp in size, comprising 37 genes and GC content of 42.24%. Furthermore, the phylogenetic tree indicated a well-supported relationship between the representatives of Tonatia into the subfamily Phyllostominae.

A genomic portrait of Sturnira parvidens: mitochondrial chromosome, repetitive elements, and microsatellite discovery

The yellow-shouldered bat, Sturnira parvidens (Phyllostomidae), is an abundant and widespread species in southern North America and Mesoamerica. Despite its important ecological role, no genomic resources exist for this species. Using low-coverage short Illumina 150 bp pair-end reads sequencing, this study reports the mitochondrial chromosome and nuclear repetitive elements, including microsatellites, in S. parvidens. The mitochondrial genome of S. parvidens is 16,612 bp in length and is comprised of 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. Repetitive elements constituted ~67% of the nuclear genome while ~33% of the genome represented single- or low-copy sequences. A moderate proportion of repetitive sequences (31% putative families) could not be assigned to known repeat element families. Considering only annotated repetitive elements, the most ubiquitous repetitive elements belonged to Class I-LINE and Satellite DNA, which were considerably more abundant than Class I-LTR elements and Class II-DNA transposons (TcMar-Mariner and hAT-Charlie). A total of 193 microsatellites were identified.

Complete mitochondrial genomes reveal robust phylogenetic signals and evidence of positive selection in horseshoe bats

BACKGROUND

In genus Rhinolophus, species in the Rhinolophus philippinensis and R. macrotis groups are unique because the horseshoe bats in these group have relatively low echolocation frequencies and flight speeds compared with other horseshoe bats with similar body size. The different characteristics among bat species suggest particular evolutionary processes may have occurred in this genus. To study the adaptive evidence in the mitochondrial genomes (mitogenomes) of rhinolophids, especially the mitogenomes of the species with low echolocation frequencies, we sequenced eight mitogenomes and used them for comparative studies of molecular phylogeny and adaptive evolution.

RESULTS

Phylogenetic analysis using whole mitogenome sequences produced robust results and provided phylogenetic signals that were better than those obtained using single genes. The results supported the recent establishment of the separate macrotis group. The signals of adaptive evolution discovered in the Rhinolophus species were tested for some of the codons in two genes (ND2 and ND6) that encode NADH dehydrogenases in oxidative phosphorylation system complex I. These genes have a background of widespread purifying selection. Signals of relaxed purifying selection and positive selection were found in ND2 and ND6, respectively, based on codon models and physicochemical profiles of amino acid replacements. However, no pronounced overlap was found for non-synonymous sites in the mitogenomes of all the species with low echolocation frequencies. A signal of positive selection for ND5 was found in the branch-site model when R. philippinensis was set as the foreground branch.

CONCLUSIONS

The mitogenomes provided robust phylogenetic signals that were much more informative than the signals obtained using single mitochondrial genes. Two mitochondrial genes that encoding proteins in the oxidative phosphorylation system showed some evidence of adaptive evolution in genus Rhinolophus and the positive selection signals were tested for ND5 in R. philippinensis. These results indicate that mitochondrial protein-coding genes were targets of adaptive evolution during the evolution of Rhinolophus species, which might have contributed to a diverse range of acoustic adaptations in this genus.

Geographic Variation of Phyllodiaptomus tunguidus Mitogenomes: Genetic Differentiation and Phylogeny

Phyllodiaptomus tunguidus (Copepoda: Calanoida) is largely endemic to and widespread in freshwater in southern China, where it inhabits a complex landscape from lowland to highland across an elevation gradient of 2000m. A deep genetic differentiation can be expected between its most distant geographic populations. Here, we sequenced nine mitogenomes from diverse populations. All mitogenomes contained 37 genes, including 13 protein-coding genes (PCG), two rRNA genes, 22 tRNA genes and one control region. Their base composition, genetic distance and tRNA structure indeed revealed a wide differentiation between mitogenomes. Two P. tunguidus from Guangxi near Vietnam differed from the other seven by up to 10.1%. Their tRNA-Arg had a complete clover-leaf structure, whereas that of the others did not contain an entire dihydrouridine arm. The nine mitogenomes also differed in the length of rRNA. NJ, ML, and Bayesian analyses all split them into two clades, viz. the two P. tunguidus from Guangxi (Clade 1), and the other seven (Clade 2). Both the structure and phylogeny of the mitogenomes suggest that P. tunguidus has complex geographic origin, and its populations in Clade 1 have long lived in isolation from those in Clade 2. They currently reach the level of subspecies or cryptic species. An extensive phylogenetic analysis of Copepoda further verified that Diaptomidae is the most recently diverging family in Calanoida and that P. tunguidus is at the evolutionary apex of the family.

Assessing Confidence in Root Placement on Phylogenies: An Empirical Study Using Non-Reversible Models for Mammals

Using time-reversible Markov models is a very common practice in phylogenetic analysis, because although we expect many of their assumptions to be violated by empirical data, they provide high computational efficiency. However, these models lack the ability to infer the root placement of the estimated phylogeny. In order to compensate for the inability of these models to root the tree, many researchers use external information such as using outgroup taxa or additional assumptions such as molecular clocks. In this study, we investigate the utility of nonreversible models to root empirical phylogenies and introduce a new bootstrap measure, the rootstrap, which provides information on the statistical support for any given root position. [Bootstrap; nonreversible models; phylogenetic inference; root estimation.]

The complete mitochondrial genome of the Honduran white bat Ectophylla alba (Allen 1982) (Chiroptera: Phyllostomidae)

The Honduran white bat, Ectophylla alba (Allen 1982), is one of eight species belonging to the family Phyllostomidae that exclusively roosts in tents. Due to its restricted distribution, habitat specificity, and diet requirements, E. alba has been strongly affected by habitat loss and fragmentation during the last decade. In this study, we developed the first genomic resource for this species; we assembled and analyzed in detail the complete mitochondrial genome of E. alba. The mitogenome of E. alba is 16,664 bp in length and is comprised of 13 protein coding genes (PCGs), 2 ribosomal RNA genes, 22 transfer RNA genes (tRNAs), and a putative Control Region (CR) 1,232 bp in length. Gene arrangement in the mitochondrial chromosome of E. alba is identical to that reported before in other species of co-familiar bats. All PCGs are under purifying selection, with atp8 experiencing the least selective pressure. In all PCGs, codons ending with adenine are preferred over others ending in thymine and cytosine. Except tRNA-Serine 1, all tRNAs exhibit a cloverleaf secondary structure. The CR of E. alba exhibits three domains commonly described in other mammals, including bats; extended terminal associated sequences (ETAS), central, and conserved sequence block (CSB). A ML phylogenetic reconstruction of the family Phyllostomidae based on all 13 mitochondrial PCGs confirms the monophyletic status of the subfamily Sternodermatinae and indicates the close relationship between E. alba and the genus Artibeus. This is the first genomic resource developed for E. alba and represents the first step to improving our understanding of the genomic underpinnings involved in the evolution of specialization as well as acclimatization and adaptation to local and global change of specialist bats.

Sequencing and analysis of the complete mitochondrial genome of the lesser bandicoot rat (Bandicota bengalensis) from China and its phylogenetic analysis

The complete mitogenome sequence of the lesser bandicoot rat (Bandicota bengalensis Gray and Hardwicke, 1833) was determined using long PCR. The genome was 16,327 bp in length and contained 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, 1 origin of L strand replication and 1 control region. The overall base composition of the heavy strand is A (34.2%), C (24.9%), T (28.5%) and G (12.4%). The base compositions present clearly the A-T skew, which is most obviously in the control region and protein-coding genes. Mitochondrial genome analyses based on MP, ML, NJ and Bayesian analyses yielded identical phylogenetic trees. This study verifies the evolutionary status of Bandicota bengalensis in Muridae at the molecular level. The mitochondrial genome would be a significant supplement for the Bandicota bengalensis genetic background. The two Bandicota species formed a monophyletic group with the high bootstrap value (100%) in all examinations.

Sequencing and analysis of the complete mitochondrial genome of Ochotona hyperborea from China and its phylogenetic analysis

The complete mitogenome sequence of Ochotona hyperborea was determined using long PCR. The genome was 17,063 bp in length and contained 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, one origin of L strand replication, and one control region. The overall base composition of the heavy strand is A (31.1%), C (28.7%), T (26.3%), and G (13.9%). The base compositions present clearly the A-T skew, which is most obvious in the control region and protein-coding genes. Mitochondrial genome analyses based on MP, ML, NJ, and Bayesian analyses yielded identical phylogenetic trees. This study verifies the evolutionary status of Ochotona hyperborea in Ochotonidae at the molecular level. The mitochondrial genome would be a significant supplement for the Ochotona hyperborea genetic background. The eight Ochotona species formed a monophyletic group with the high bootstrap value (100%) in all examinations.

Sequencing and analysis of the complete mitochondrial genome of Micromys erythrotis from China and its phylogenetic analysis

The complete mitogenome sequence of Micromys erythrotis was determined using long PCR. The genome was 16,238 bp in length and contained 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, 1 origin of L strand replication and 1 control region. The overall base composition of the heavy strand is A (33.7%), C (24.8%), T (29.1%) and G (12.4%). The base compositions present clearly the A-T skew, which is most obviously in the control region and protein-coding genes. Mitochondrial genome analyses based on MP, ML, NJ and Bayesian analyses yielded identical phylogenetic trees. This study verifies the evolutionary status of Micromys erythrotis in Muridae at the molecular level. The mitochondrial genome would be a significant supplement for the Micromys erythrotis genetic background.

Sequence comparison of the mitochondrial genomes in two species of the genus Nerita (Gastropoda: Neritimorpha: Neritidae): phylogenetic implications and divergence time estimation for Neritimorpha

Many Nerita species live in warm-water environments, and they are some of the few organisms from the intertidal zone that can live in both freshwater and seawater. Previous comparative studies of the mitogenomes of Nerita species suggest that the genome rearrangements are very conservative. Generally, the species possess a set of similar mitochondrial gene arrangements, but nucleotide sequences can be used to elucidate phylogenetic relationships at various levels of divergence. Here, the mitogenomes of Nerita undata and Nerita balteata were sequenced and found to be 15,583 bp and 15,571 bp, respectively. The mitogenomes of both species contain 13 protein-coding genes, 22 tRNA genes, and two rRNA genes. The nucleotides of the two mitogenomes are highly similar, with the same gene composition and genomic organization as those present in other Nerita species. The tRNA secondary structures were different from those of other gastropods: trnS2 is not folded into typical secondary structures, and the dihydrouridine (DHU) arm simply forms a loop. The phylogenetic analysis showed that Neritimorpha is a sister group of Vetigastropoda and Caenogastropoda. Nerita balteata is a sister group of Nerita versicolor and Nerita undata, and all three species belong to Neritimorpha. This study contributes towards the comparative mitogenomic analysis of Neritidae and phylogenetic considerations among Neritimorpha species. The estimation of divergence time revealed that the two Nerita species were differentiated in the late Paleogene of the Cenozoic Era, and their evolution may be related to environmental changes.

The published complete mitochondrial genome of Eptesicus serotinus is a chimera of Vespertilio sinensis and Hypsugo alaschanicus (Mammalia: Chiroptera)

The mitogenome of Eptesicus serotinus (Serotine bat) was published in 2013 with GenBank accession number KF111725 and NCBI Reference Sequence number NC_022474. This sequence was placed with Vespertilio sinensis (Asian parti-colored bat) in a COI gene tree but with Hypsugo alashanicus (Alashanian pipistrelle) in a cytochrome b gene tree. Direct comparison of mitogenomes showed that 92.4% of this mitogenome is similar to Vespertilio sinensis, 5.9% to Hypsugo alaschanicus, and that 1.6% of the mitogenome could not be attributed to either species, or any other species. This mitogenome has been re-used in at least 17 phylogenies. Our findings suggest that mitogenomes are best verified with multiple gene trees, followed by direct comparison of sequences. We conclude that greater vigilance is warranted to ensure that problematic sequences do not enter the scientific record, and are not re-used in subsequent studies.

Sequencing and analysis of the complete mitochondrial genome of Blarinella griselda from China and its phylogenetic analysis

The complete mitogenome sequence of Blarinella griselda was determined using long PCR. The genome was 16,947 bp in length and contained 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, 1 origin of L strand replication and 1 control region. The overall base composition of the heavy strand is A (33.1%), C (22.6%), T (31.6%) and G (12.7%). The base compositions present clearly the A-T skew, which is most obviously in the control region and protein-coding genes. Mitochondrial genome analyses based on MP, ML, NJ and Bayesian analyses yielded identical phylogenetic trees. This study verifies the evolutionary status of Blarinella griselda in Soricidae at the molecular level. The mitochondrial genome would be a significant supplement for the Blarinella griselda genetic background. The three Blarinella species formed a monophyletic group with the high bootstrap value (100%) in all examinations.

Sequencing and analysis of the complete mitochondrial genome of Crocidura tanakae from China and its phylogenetic analysis

The complete mitogenome sequence of Crocidura tanakae was determined using long PCR. The genome was 16,969 bp in length and contained 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, 1 origin of L strand replication and 1 control region. The overall base composition of the heavy strand is A (32.5%), C (22.3%), T (31.9%), and G (13.3%). The base compositions present clearly the A-T skew, which is most obviously in the control region and protein-coding genes. Mitochondrial genome analyses based on MP, ML, NJ, and Bayesian analyses yielded identical phylogenetic trees. The five Crocidura species formed a monophyletic group with the high bootstrap value (100%) in all examinations. This study verifies the evolutionary status of C. tanakae in Soricidae at the molecular level. The mitochondrial genome would be a significant supplement for the C. tanakae genetic background.

The complete mitochondrial genome of the striped hamster (Cricetulus barabesis) China and its phylogenetic analysis

The complete mitogenome sequence of the striped hamster was determined using long PCR. The genome was 16,282 bp in length and contained 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, one origin of L strand replication, and one control region. The overall base composition of the heavy strand is A (33.7%), C (22.8%), T (30.5%), and G (13.0%). The base compositions present clearly the A-T skew, which is most obviously in the control region and protein-coding genes. Mitochondrial genome analyses based on MP, ML, NJ, and Bayesian analyses yielded identical phylogenetic trees. Results of phylogenetic analysis showed that Cricetulus had close relationship with Meriones. This study verifies the evolutionary status of the striped hamster in Cricetulus at the molecular level. The mitochondrial genome would be a significant supplement for the striped hamster genetic background. Results of phylogenetic analysis showed that the striped hamster had close relationship with C. griseus in Cricetulus.

Sequencing and analysis of the complete mitochondrial genome of Chodsigoa hoffmanni from China and its phylogenetic analysis

The complete mitogenome sequence of Chodsigoa hoffmanni was determined using long PCR. The genome was 17,138 bp in length and contained 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, one origin of L strand replication, and one control region. The overall base composition of the heavy strand is A (32.8%), C (24.4%), T (29.8%), and G (13.0%). The base compositions present clearly the A-T skew, which is most obviously in the control region and protein-coding genes. Mitochondrial genome analyses based on MP, ML, NJ, and Bayesian analyses yielded identical phylogenetic trees. Chodsigoa hoffmanni is the first species to have been reported on the mitochondrial genome in Chodsigoa genus. This study verifies the evolutionary status of C. hoffmanni in Soricidae at the molecular level. The mitochondrial genome would be a significant supplement for the C. hoffmanni genetic background.

Sequencing and analysis of the complete mitochondrial genome of the taiga shrew (Sorex isodon) from China

The complete mitogenome sequence of the taiga shrew (Sorex isodon) was determined using long PCR. The genome was 17,008 bp in length and contained 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, one origin of L strand replication and one control region. The overall base composition of the heavy strand is A (32.5%), C (24.5%), T (28.5%), and G (13.5%). The base compositions present clearly the A-T skew, which is most obviously in the control region and protein-coding genes. The extended termination-associated sequence domain, the central conserved domain and the conserved sequence block domain are defined in the mitochondrial genome control region of the taiga shrew. Mitochondrial genome analyses based on MP, ML, NJ, and Bayesian analyses yielded identical phylogenetic trees. The eight Sorex species formed a monophyletic group with the high bootstrap value (100%) in all examinations.

Mitochondrial genomes of African pangolins and insights into evolutionary patterns and phylogeny of the family Manidae

BACKGROUND

This study used next generation sequencing to generate the mitogenomes of four African pangolin species; Temminck's ground pangolin (Smutsia temminckii), giant ground pangolin (S. gigantea), white-bellied pangolin (Phataginus tricuspis) and black-bellied pangolin (P. tetradactyla).

RESULTS

The results indicate that the mitogenomes of the African pangolins are 16,558 bp for S. temminckii, 16,540 bp for S. gigantea, 16,649 bp for P. tetradactyla and 16,565 bp for P. tricuspis. Phylogenetic comparisons of the African pangolins indicated two lineages with high posterior probabilities providing evidence to support the classification of two genera; Smutsia and Phataginus. The total GC content between African pangolins was observed to be similar between species (36.5% - 37.3%). The most frequent codon was found to be A or C at the 3rd codon position. Significant variations in GC-content and codon usage were observed for several regions between African and Asian pangolin species which may be attributed to mutation pressure and/or natural selection. Lastly, a total of two insertions of 80 bp and 28 bp in size respectively was observed in the control region of the black-bellied pangolin which were absent in the other African pangolin species.

CONCLUSIONS

The current study presents reference mitogenomes of all four African pangolin species and thus expands on the current set of reference genomes available for six of the eight extant pangolin species globally and represents the first phylogenetic analysis with six pangolin species using full mitochondrial genomes. Knowledge of full mitochondrial DNA genomes will assist in providing a better understanding on the evolution of pangolins which will be essential for conservation genetic studies.

Sequencing and analysis of the complete mitochondrial genome of the Ussuri shrew (Sorex mirabilis) from China

The complete mitogenome sequence of the Ussuri shrew (Sorex mirabilis) was determined using long PCR. The genome was 17,315 bp in length and contained 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, one origin of L strand replication, and one control region. The overall base composition of the heavy strand is A (32.6%), C (25.2%), T (28.8%), and G (13.4%). The base compositions present clearly the A–T skew, which is most obviously in the control region and protein-coding genes. The extended termination-associated sequence domain, the central conserved domain and the conserved sequence block domain are defined in the mitochondrial genome control region of the Ussuri shrew. Mitochondrial genome analyses based on MP, ML, NJ, and Bayesian analyses yielded identical phylogenetic trees. The five Sorex species formed a monophyletic group with the high bootstrap value (100%) in all examinations.

Sequencing and analysis of the complete mitochondrial genome of the slender shrew (Sorex gracillimus) from China

The complete mitogenome sequence of the slender shrew (Sorex gracillimus) was determined using long PCR. The genome was 17,002 bp in length and contained 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, one origin of L strand replication, and one control region. The overall base composition of the heavy strand is A (32.5%), C (25.5%), T (28.5%), and G (13.5%). The base compositions present clearly the A–T skew, which is most obviously in the control region and protein-coding genes. The extended termination-associated sequence domain, the central conserved domain and the conserved sequence block domain are defined in the mitochondrial genome control region of the slender shrew. Mitochondrial genome analyses based on MP, ML, NJ, and Bayesian analyses yielded identical phylogenetic trees. The five Sorex species formed a monophyletic group with the high bootstrap value (100%) in all examinations.

Sequencing and analysis of the complete mitochondrial genome of the masked shrew (Sorex caecutiens) from China

The complete mitogenome sequence of the masked shrew (Sorex caecutiens) was determined using long PCR. The genome was 17,096 bp in length and contained 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, one origin of L strand replication, and one control region. The overall base composition of the heavy strand is A (32.9%), C (24.5%), T (29.3%), and G (13.3%). The base compositions present clearly the A–T skew, which is most obviously in the control region and protein-coding genes. The extended termination-associated sequence domain, the central conserved domain and the conserved sequence block domain are defined in the mitochondrial genome control region of the masked shrew. Mitochondrial genome analyses based on MP, ML, NJ, and Bayesian analyses yielded identical phylogenetic trees. The five Sorex species formed a monophyletic group with the high bootstrap value (100%) in all examinations.

Sequencing and analysis of the complete mitochondrial genome of flat-skulled shrew (Sorex roboratus) from China

The complete mitogenome sequence of flat-skulled shrew (Sorex roboratus) was determined using long PCR. The genome was 17,153 bp in length and contained 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, 1 origin of L strand replication and 1 control region. The overall base composition of the heavy strand is A (33.1%), C (24.4%), T (29.4%), and G (13.1%). The base compositions present clearly the A–T skew, which is most obviously in the control region and protein-coding genes. The extended termination-associated sequence domain, the central conserved domain and the conserved sequence block domain are defined in the mitochondrial genome control region of flat-skulled shrew. Mitochondrial genome analyses based on MP, ML, NJ and Bayesian analyses yielded identical phylogenetic trees. The five Sorex species formed a monophyletic group with the high bootstrap value (100%) in all examinations.

Mitochondrial genome of Pteronotus personatus (Chiroptera: Mormoopidae): comparison with selected bats and phylogenetic considerations

We described the complete mitochondrial genome (mitogenome) of the Wagner's mustached bat, Pteronotus personatus, a species belonging to the family Mormoopidae, and compared it with other published mitogenomes of bats (Chiroptera). The mitogenome of P. personatus was 16,570 bp long and contained a typically conserved structure including 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and one control region (D-loop). Most of the genes were encoded on the H-strand, except for eight tRNA and the ND6 genes. The order of protein-coding and rRNA genes was highly conserved in all mitogenomes. All protein-coding genes started with an ATG codon, except for ND2, ND3, and ND5, which initiated with ATA, and terminated with the typical stop codon TAA/TAG or the codon AGA. Phylogenetic trees constructed using Maximum Parsimony, Maximum Likelihood, and Bayesian inference methods showed an identical topology and indicated the monophyly of different families of bats (Mormoopidae, Phyllostomidae, Vespertilionidae, Rhinolophidae, and Pteropopidae) and the existence of two major clades corresponding to the suborders Yangochiroptera and Yinpterochiroptera. The mitogenome sequence provided here will be useful for further phylogenetic analyses and population genetic studies in mormoopid bats.

Comparative Mitogenomics of the Genus Odontobutis (Perciformes: Gobioidei: Odontobutidae) Revealed Conserved Gene Rearrangement and High Sequence Variations

To understand the molecular evolution of mitochondrial genomes (mitogenomes) in the genus Odontobutis, the mitogenome of Odontobutis yaluensis was sequenced and compared with those of another four Odontobutis species. Our results displayed similar mitogenome features among species in genome organization, base composition, codon usage, and gene rearrangement. The identical gene rearrangement of trnS-trnL-trnH tRNA cluster observed in mitogenomes of these five closely related freshwater sleepers suggests that this unique gene order is conserved within Odontobutis. Additionally, the present gene order and the positions of associated intergenic spacers of these Odontobutis mitogenomes indicate that this unusual gene rearrangement results from tandem duplication and random loss of large-scale gene regions. Moreover, these mitogenomes exhibit a high level of sequence variation, mainly due to the differences of corresponding intergenic sequences in gene rearrangement regions and the heterogeneity of tandem repeats in the control regions. Phylogenetic analyses support Odontobutis species with shared gene rearrangement forming a monophyletic group, and the interspecific phylogenetic relationships are associated with structural differences among their mitogenomes. The present study contributes to understanding the evolutionary patterns of Odontobutidae species.

The complete mitochondrial genome of Eurasian water shrew (Neomys fodiens)

The complete mitogenome sequence of Eurasian water shrew (Neomys fodiens) was determined using long PCR. The genome was 17,260 bp in length and contained 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes and 1 control region. The overall base composition of the heavy strand is A (33.0%), C (23.2%), T (30.9%) and G (12.9%). The base compositions present clearly the A-T skew, which is most obviously in the control region and protein-coding genes. The extended termination-associated sequence domain, the central conserved domain and the conserved sequence block domain are defined in the mitochondrial genome control region of Eurasian water shrew. Mitochondrial genomes analyses based on MP, ML, NJ and Bayesian analyses yielded identical phylogenetic trees. Neomys elegans, Soriculus fumidus, and N. fodiens formed a monophyletic group with the high bootstrap value (100%) in all examinations. This clade with the Anourosorex squamipes as the sister taxon to Sorex.

Are mutagenic non D-loop direct repeat motifs in mitochondrial DNA under a negative selection pressure?

Non D-loop direct repeats (DRs) in mitochondrial DNA (mtDNA) have been commonly implicated in the mutagenesis of mtDNA deletions associated with neuromuscular disease and ageing. Further, these DRs have been hypothesized to put a constraint on the lifespan of mammals and are under a negative selection pressure. Using a compendium of 294 mammalian mtDNA, we re-examined the relationship between species lifespan and the mutagenicity of such DRs. Contradicting the prevailing hypotheses, we found no significant evidence that long-lived mammals possess fewer mutagenic DRs than short-lived mammals. By comparing DR counts in human mtDNA with those in selectively randomized sequences, we also showed that the number of DRs in human mtDNA is primarily determined by global mtDNA properties, such as the bias in synonymous codon usage (SCU) and nucleotide composition. We found that SCU bias in mtDNA positively correlates with DR counts, where repeated usage of a subset of codons leads to more frequent DR occurrences. While bias in SCU and nucleotide composition has been attributed to nucleotide mutational bias, mammalian mtDNA still exhibit higher SCU bias and DR counts than expected from such mutational bias, suggesting a lack of negative selection against non D-loop DRs.

Sequencing and analysis of the complete mitochondrial genome of tundra shrew (Sorex tundrensis) from China

The complete mitogenome sequence of tundra shrew (Sorex tundrensis) was determined using long PCR. The genome was 17,444 bp in length and contained 13 protein-coding genes, 2 ribosomal RNA genes, 23 transfer RNA genes, 1 origin of L strand replication and 1 control region. The overall base composition of the heavy strand is A (32.9%), C (24.8%), T (29.0%), and G (13.3%). The base compositions present clearly the A-T skew, which is most obviously in the control region and protein-coding genes. The extended termination-associated sequence domain, the central conserved domain and the conserved sequence block domain are defined in the mitochondrial genome control region of tundra shrew. Mitochondrial genome analyses based on MP, ML, NJ and Bayesian analyses yielded identical phylogenetic trees. The three Sorex species formed a monophyletic group with the high bootstrap value (100 %) in all examinations.

The complete mitochondrial genome of Vespertilio sinensis from China

The complete mitogenome sequence of Vespertilio sinensis was determined using long PCR. The genome was 17,146 bp in length and contained 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes and 1 control region. The overall base composition of the heavy strand is A (32.73%), C (23.64%), T (29.72%), and G (13.91%). The base compositions present clearly the A-T skew, which is most obviously in the control region and protein-coding genes. The extended termination-associated sequence domain, the central conserved domain and the conserved sequence block domain are defined in the mitochondrial genome control region of V. sinensis.

Complete mitochondrial genome of Cynopterus sphinx (Pteropodidae: Cynopterus)

We have characterized the complete mitochondrial genome of Cynopterus sphinx (Pteropodidae: Cynopterus) and described its organization in this study. The total length of C. sphinx complete mitochondrial genome was 16,895 bp with the base composition of 32.54% A, 14.05% G, 25.82% T and 27.59% C. The complete mitochondrial genome included 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes (12S rRNA and 16S rRNA) and 1 control region (D-loop). The control region was 1435 bp long with the sequence CATACG repeat 64 times. Three protein-coding genes (ND1, COI and ND4) were ended with incomplete stop codon TA or T.

Next-generation sequencing and phylogenetic signal of complete mitochondrial genomes for resolving the evolutionary history of leaf-nosed bats (Phyllostomidae)

Leaf-nosed bats (Phyllostomidae) are one of the most studied groups within the order Chiroptera mainly because of their outstanding species richness and diversity in morphological and ecological traits. Rapid diversification and multiple homoplasies have made the phylogeny of the family difficult to solve using morphological characters. Molecular data have contributed to shed light on the evolutionary history of phyllostomid bats, yet several relationships remain unresolved at the intra-familial level. Complete mitochondrial genomes have proven useful to deal with this kind of situation in other groups of mammals by providing access to a large number of molecular characters. At present, there are only two mitogenomes available for phyllostomid bats hinting at the need for further exploration of the mitogenomic approach in this group. We used both standard Sanger sequencing of PCR products and next-generation sequencing (NGS) of shotgun genomic DNA to obtain new complete mitochondrial genomes from 10 species of phyllostomid bats, including representatives of major subfamilies, plus one outgroup belonging to the closely-related mormoopids. We then evaluated the contribution of mitogenomics to the resolution of the phylogeny of leaf-nosed bats and compared the results to those based on mitochondrial genes and the RAG2 and VWF nuclear makers. Our results demonstrate the advantages of the Illumina NGS approach to efficiently obtain mitogenomes of phyllostomid bats. The phylogenetic signal provided by entire mitogenomes is highly comparable to the one of a concatenation of individual mitochondrial and nuclear markers, and allows increasing both resolution and statistical support for several clades. This enhanced phylogenetic signal is the result of combining markers with heterogeneous evolutionary rates representing a large number of nucleotide sites. Our results illustrate the potential of the NGS mitogenomic approach for resolving the evolutionary history of phyllostomid bats based on a denser species sampling.