The complete nucleotide sequence of a snake (Dinodon semicarinatus) mitochondrial genome with two identical control regions

Complete mitochondrial genome and phylogenetic analysis of Mancinella alouina

BACKGROUND

The Muricidae family in the Class Gastropoda comprises numerous species with a vast range of morphological features and a worldwide presence. The phylogeny of the Muricidae has been analyzed in previous studies; however, the evolutionary relationships among the main branches of the Muricidae remain unknown.

METHODS AND RESULTS

In the present study, the mitochondrial genome of Mancinella alouina was sequenced. The mitochondrial genome was found to be 16,671 bp in length and made up of 37 genes (13 protein-coding genes, 22 transfer RNA and 2 ribosomal RNA genes). The genome has an A-T-rich region (66.5% A + T content) and all of the PCGs use the ATN start codon and the TAG or TAA stop codons. The mitochondrial gene arrangement of Mancinella alouina is similar to that of other Muricidae, except for Ocinebrellus inornatus and Ceratostoma burnetti. On the basis of a flexible molecular clock model, time-calibrated phylogenetic results indicate that the genus Mancinella diverged roughly 18.09 Mya, and that the family Muricidae emerged in the Late Cretaceous.

CONCLUSIONS

This study reveals the structural and sequence information features of the mitochondrial genome of Mancinella alouina. This study provides evidence for the relationships within the family Muricidae at the molecular level, and infer the divergence time. The results of phylogenetic analyses strongly support the current classification.

Comparative Mitogenomic Analysis of Two Snake Eels Reveals Irregular Gene Rearrangement and Phylogenetic Implications of Ophichthidae

The family Ophichthidae has the largest number and the most various species (about 359 valid species) in the order Anguilliformes worldwide. Both morphological and molecular characteristics have been used to assess their taxonomic status. However, due to the ambiguous morphological features, molecular data such as mitochondrial DNA sequences have been implemented for the correct identification and classification of these fishes. In this study, the gene arrangement and structure characteristics of two Ophichthidae mitochondrial genomes were investigated for the first time. The total mitogenome lengths of O. evermanni and O. erabo were 17,759 bp and 17,856 bp, respectively. Comparing with the ancestral mitochondrial gene order, the irregular gene rearrangement happened between ND6 and tRNA-Pro (P) genes with another similar control region emerging between tRNA-Thr (T) and ND6 genes, which could be explained by the tandem duplication and random loss (TDRL) model appropriately. ML phylogenetic tree demonstrated that the family Ophichthidae was monophyletic origin, but genus Ophichthus might be polyphyletic because of the confused cluster relationships among different species.

A comparison of gene organisations and phylogenetic relationships of all 22 squamate species listed in South Korea using complete mitochondrial DNA

Studies using complete mitochondrial genome data have the potential to increase our understanding on gene organisations and evolutionary species relationships. In this study, we compared complete mitochondrial genomes between all 22 squamate species listed in South Korea. In addition, we constructed Maximum Parsimony (MP), Maximum Likelihood (ML) and Bayesian Inference (BI) phylogenetic trees using 13 mitochondrial protein-coding genes. The mitochondrial genes for all six species in the suborder Sauria followed the same organisation as the sequenced Testudines (turtle) outgroup. In contrast, 16 snake species in the suborder Serpentes contained some gene organisational variations. For example, all snake species contained a second control region (CR2), while three species in the family Viperidae had a translocated tRNA-Pro gene region. In addition, the snake species, Elapheschrenckii, carried a tRNA-Pro pseudogene. We were also able to identify a translocation of a tRNA-Asn gene within the five tRNA (WANCY gene region) gene clusters for two true sea snake species in the subfamily Hydrophiinae. Our BI phylogenetic tree was also well fitted against currently known Korean squamate phylogenetic trees, where each family and genus unit forms monophyletic clades and the suborder Sauria is paraphyletic to the suborder Serpentes. Our results may form the basis for future northeast Asian squamate phylogenetic studies.

Complete mitochondrial genomes of Boigakraepelini and Hebiuscraspedogaster (Reptilia, Squamata, Colubridae) and their phylogenetic implications

The complete sequence of the mitochondrial genome is a powerful tool for studying phylogenetic relationships and molecular evolution in various species. In this work, the mitogenomes of Boigakraepelini and Hebiuscraspedogaster were sequenced and characterized for the first time. The lengths of the B.kraepelini and H.craspedogaster mitogenomes were 17,124 bp and 17,120 bp, respectively, and both included 13 protein-coding genes, 22 tRNAs, two rRNAs and two control regions. The arrangements of these mitochondrial genes were the same in B.kraepelini and H.craspedogaster. In addition, both genome compositions showed A+T bias (59.03%, 60.93%) and had positive AT skews (0.179, 0.117) and negative GC skews (-0.397, -0.348). The phylogenetic results illustrated a close relationship between B.kraepelini and the genus Lycodon. Moreover, H.craspedogaster was clustered with other Hebius snakes and closely related to other Natricinae species. These results will provide references for further research on the phylogeny of Colubridae.

Complete Mitogenomes of Polypedates Tree Frogs Unveil Gene Rearrangement and Concerted Evolution within Rhacophoridae

New developments in sequencing technology and nucleotide analysis have allowed us to make great advances in reconstructing anuran phylogeny. As a clade of representative amphibians that have radiated from aquatic to arboreal habitats, our understanding of the systematic status and molecular biology of rhacophorid tree frogs is still limited. We determined two new mitogenomes for the genus Polypedates (Rhacophoridae): P. impresus and P. mutus. We conducted comparative and phylogenetic analyses using our data and seven other rhacophorid mitogenomes. The mitogenomes of the genera Polypedates, Buergeria, and Zhangixalus were almost identical, except that the ATP8 gene in Polypedates had become a non-coding region; Buergeria maintained the legacy "LTPF" tRNA gene cluster compared to the novel "TLPF" order in the other two genera; and B. buergeri and Z. dennysi had no control region (CR) duplication. The resulting phylogenetic relationship supporting the above gene rearrangement pathway suggested parallel evolution of ATP8 gene loss of function (LoF) in Polypedates and CR duplication with concerted evolution of paralogous CRs in rhacophorids. Finally, conflicting topologies in the phylograms of 185 species reflected the advantages of phylogenetic analyses using multiple loci.

Evolution of the Noncoding Features of Sea Snake Mitochondrial Genomes within Elapidae

Mitochondrial genomes of four elapid snakes (three marine species [Emydocephalus ijimae, Hydrophis ornatus, and Hydrophis melanocephalus], and one terrestrial species [Sinomicrurus japonicus]) were completely sequenced by a combination of Sanger sequencing, next-generation sequencing and Nanopore sequencing. Nanopore sequencing was especially effective in accurately reading through long tandem repeats in these genomes. This led us to show that major noncoding regions in the mitochondrial genomes of those three sea snakes contain considerably long tandem duplications, unlike the mitochondrial genomes previously reported for same and other sea snake species. We also found a transposition of the light-strand replication origin within a tRNA gene cluster for the three sea snakes. This change can be explained by the Tandem Duplication—Random Loss model, which was further supported by remnant intervening sequences between tRNA genes. Mitochondrial genomes of true snakes (Alethinophidia) have been shown to contain duplicate major noncoding regions, each of which includes the control region necessary for regulating the heavy-strand replication and transcription from both strands. However, the control region completely disappeared from one of the two major noncoding regions for two Hydrophis sea snakes, posing evolutionary questions on the roles of duplicate control regions in snake mitochondrial genomes. The timing and molecular mechanisms for these changes are discussed based on the elapid phylogeny.

Divergent Specialization of Simple Venom Gene Profiles among Rear-Fanged Snake Genera (Helicops and Leptodeira, Dipsadinae, Colubridae)

Many venomous animals express toxins that show extraordinary levels of variation both within and among species. In snakes, most studies of venom variation focus on front-fanged species in the families Viperidae and Elapidae, even though rear-fanged snakes in other families vary along the same ecological axes important to venom evolution. Here we characterized venom gland transcriptomes from 19 snakes across two dipsadine rear-fanged genera (Leptodeira and Helicops, Colubridae) and two front-fanged genera (Bothrops, Viperidae; Micrurus, Elapidae). We compared patterns of composition, variation, and diversity in venom transcripts within and among all four genera. Venom gland transcriptomes of rear-fanged Helicops and Leptodeira and front-fanged Micrurus are each dominated by expression of single toxin families (C-type lectins, snake venom metalloproteinase, and phospholipase A2, respectively), unlike highly diverse front-fanged Bothrops venoms. In addition, expression patterns of congeners are much more similar to each other than they are to species from other genera. These results illustrate the repeatability of simple venom profiles in rear-fanged snakes and the potential for relatively constrained venom composition within genera.

Phylogeography of the smooth greensnake, Opheodrys vernalis (Squamata: Colubridae): divergent lineages and variable demographics in a widely distributed yet enigmatic species

Phylogeographic studies can uncover robust details about the population structure, demographics, and diversity of species. The smooth greensnake, Opheodrys vernalis, is a small, cryptic snake occupying mesic grassland and sparsely wooded habitats. Although O. vernalis has a wide geographical range, many metapopulations are patchy and some are declining. We used mitochondrial DNA and double digest restriction-site associated DNA sequencing to construct the first phylogeographic assessment of O. vernalis. Genomic analysis of 119 individuals (mitochondrial DNA) and a subset of another 45 smooth greensnakes (nuclear DNA; N = 3031 single nucleotide polymorphisms) strongly supports two longitudinally separated lineages, with admixture in the Great Lakes region. Post-Pleistocene secondary contact best explains admixture from populations advancing northwards. Overall, populations expressed low heterozygosity, variable inbreeding rates, and moderate to high differentiation. Disjunct populations in the Rocky Mountains and central Great Plains regions might be contracting relicts, whereas northerly populations in more continuous mesic habitats (e.g., Prairie Pothole region, southern Canada) had signals of population expansion. Broadly, conservation management efforts should be focused on local populations, because habitat connectivity may facilitate gene flow and genetic diversity.

Mitogenomes of Accipitriformes and Cathartiformes Were Subjected to Ancestral and Recent Duplications Followed by Gradual Degeneration

The rearrangement of 37 genes with one control region, firstly identified in Gallus gallus mitogenome, is believed to be ancestral for all Aves. However, mitogenomic sequences obtained in recent years revealed that many avian mitogenomes contain duplicated regions that were omitted in previous genomic versions. Their evolution and mechanism of duplication are still poorly understood. The order of Accipitriformes is especially interesting in this context because its representatives contain a duplicated control region in various stages of degeneration. Therefore, we applied an appropriate PCR strategy to look for duplications within the mitogenomes of the early diverged species Sagittarius serpentarius and Cathartiformes, which is a sister order to Accipitriformes. The analyses revealed the same duplicated gene order in all examined taxa and the common ancestor of these groups. The duplicated regions were subjected to gradual degeneration and homogenization during concerted evolution. The latter process occurred recently in the species of Cathartiformes as well as in the early diverged lineages of Accipitriformes, that is, Sagittarius serpentarius and Pandion haliaetus. However, in other lineages, that is, Pernis ptilorhynchus, as well as representatives of Aegypiinae, Aquilinae, and five related subfamilies of Accipitriformes (Accipitrinae, Circinae, Buteoninae, Haliaeetinae, and Milvinae), the duplications were evolving independently for at least 14-47 Myr. Different portions of control regions in Cathartiformes showed conflicting phylogenetic signals indicating that some sections of these regions were homogenized at a frequency higher than the rate of speciation, whereas others have still evolved separately.

Structural features and phylogenetic implications of Cicadellidae subfamily and two new mitogenomes leafhoppers

Complete mitochondrial genome sequences facilitate species identification and analyses of phylogenetic relationships. However, the available data are limited to the diverse and widespread insect family Cicadellidae. This study analyzes and summarizes the complete mitochondrial genome structure characteristics of 11 leafhopper subfamilies and two newly sequenced Typhlocybinae species, Empoascanara wengangensis and E. gracilis. Moreover, using 13PCGs and rRNA data to analyze the nucleotide diversity, evolution rate, and the phylogenetic relationship between the subfamilies of 56 species, verifying the taxonomic status analysis of E. wengangensis and E. gracilis. The analysis results show that the genome structures of the subfamilies and the newly sequenced two species are very similar, and the size of the CR region is significantly related to the repeat unit. However, in the entire AT-skews and CG-skews, the AT-skews of other subfamilies are all positive, and CG-skews are negative, while Empoascini of Typhlocybinae and Ledrinae are the opposite. Furthermore, among 13PCGs, the AT-skews of 13 species are all negative while CG-skews are positive, which from Empoascini in Typhlocybinae, Idiocerinae, Cicadellinae, Ledrinae, and Evacanthinae. Phylogenetic analysis shows that ML and PB analysis produce almost consistent topologies between different data sets and models, and some relationships are highly supported and remain unchanged. Mileewinae is a monophyletic group and is a sister group with Typhlocybinae, and the sister group of Evacanthinae is Ledrinae + Cicadellinae. Phylogenetic analysis grouped the two newly sequenced species with other species of Typhlocybinae, which was separated from other subfamilies, and all Erythroneurini insects gathered together. However, E. gracilis grouped into a single group, not grouped with species of the same genus (Empoascanara). This result does not match the traditional classification, and other nuclear genes or transcriptome genes may be needed to verify the result. Nucleotide diversity analysis shows that nad4 and nad5 may be evaluated as potential DNA markers defining the Cicadellidae insect species.

Comparative mitogenomes of three species in Moenkhausia: Rare irregular gene rearrangement within Characidae

Generally, a teleostean group possesses only one type or a set of similar mitochondrial gene arrangements. However, a new type of gene arrangement has been identified in the mitochondrial genomes (mitogenomes) of Moenkhausia. Here, three newly sequenced complete mitogenomes of tetras (Characidae: Moenkhausia) are presented (M. costae, M. pittieri, and M. sanctaefilomenae). The three mitogenomes had a classical circular structure, with total lengths ranging from 15,811 to 18,435 bp. Base composition analysis indicated that the sequences were biased toward adenine (A) and thymine (T), with A + T content of 54.63% in M. costae, 58.47% in M. pittieri, and 59.98% in M. sanctaefilomenae. The gene order and organization of M. sanctaefilomenae differed from those of typical teleostean mitogenomes. The genes tRNA-Ile, tRNA-Gln, and tRNA-Pro were translocated between tRNA-Trp and tRNA-Asn. One extra tRNA-Met and an extra CR were also discovered in the mitogenome. BI and ML analyses based on sequences of 38 different mitogenomes showed that M. costae and M. pittieri were classified together, and M. sanctaefilomenae was slightly further from other fish of the same genus. These results provide insight into the gene arrangement features of Characidae mitogenomes and lay the foundation for further phylogenetic studies on Characidae.

Novel gene rearrangement in the mitochondrial genome of Muraenesox cinereus and the phylogenetic relationship of Anguilliformes

The structure and gene sequence of the fish mitochondrial genome are generally considered to be conservative. However, two types of gene arrangements are found in the mitochondrial genome of Anguilliformes. In this paper, we report a complete mitogenome of Muraenesox cinereus (Anguilliformes: Muraenesocidae) with rearrangement phenomenon. The total length of the M. cinereus mitogenome was 17,673 bp, and it contained 13 protein-coding genes, two ribosomal RNAs, 22 transfer RNA genes, and two identical control regions (CRs). The mitochondrial genome of M. cinereus was obviously rearranged compared with the mitochondria of typical vertebrates. The genes ND6 and the conjoint trnE were translocated to the location between trnT and trnP, and one of the duplicated CR was translocated to the upstream of the ND6. The tandem duplication and random loss is most suitable for explaining this mitochondrial gene rearrangement. The Anguilliformes phylogenetic tree constructed based on the whole mitochondrial genome well supports Congridae non-monophyly. These results provide a basis for the future Anguilliformes mitochondrial gene arrangement characteristics and further phylogenetic research.

New view on the organization and evolution of Palaeognathae mitogenomes poses the question on the ancestral gene rearrangement in Aves

BACKGROUND

Bird mitogenomes differ from other vertebrates in gene rearrangement. The most common avian gene order, identified first in Gallus gallus, is considered ancestral for all Aves. However, other rearrangements including a duplicated control region and neighboring genes have been reported in many representatives of avian orders. The repeated regions can be easily overlooked due to inappropriate DNA amplification or genome sequencing. This raises a question about the actual prevalence of mitogenomic duplications and the validity of the current view on the avian mitogenome evolution. In this context, Palaeognathae is especially interesting because is sister to all other living birds, i.e. Neognathae. So far, a unique duplicated region has been found in one palaeognath mitogenome, that of Eudromia elegans.

RESULTS

Therefore, we applied an appropriate PCR strategy to look for omitted duplications in other palaeognaths. The analyses revealed the duplicated control regions with adjacent genes in Crypturellus, Rhea and Struthio as well as ND6 pseudogene in three moas. The copies are very similar and were subjected to concerted evolution. Mapping the presence and absence of duplication onto the Palaeognathae phylogeny indicates that the duplication was an ancestral state for this avian group. This feature was inherited by early diverged lineages and lost two times in others. Comparison of incongruent phylogenetic trees based on mitochondrial and nuclear sequences showed that two variants of mitogenomes could exist in the evolution of palaeognaths. Data collected for other avian mitogenomes revealed that the last common ancestor of all birds and early diverging lineages of Neoaves could also possess the mitogenomic duplication.

CONCLUSIONS

The duplicated control regions with adjacent genes are more common in avian mitochondrial genomes than it was previously thought. These two regions could increase effectiveness of replication and transcription as well as the number of replicating mitogenomes per organelle. In consequence, energy production by mitochondria may be also more efficient. However, further physiological and molecular analyses are necessary to assess the potential selective advantages of the mitogenome duplications.

Characteristics of the mitochondrial genome of Rana omeimontis and related species in Ranidae: Gene rearrangements and phylogenetic relationships

The Omei wood frog (Rana omeimontis), endemic to central China, belongs to the family Ranidae. In this study, we achieved detail knowledge about the mitogenome of the species. The length of the genome is 20,120 bp, including 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and a noncoding control region. Similar to other amphibians, we found that only nine genes (ND6 and eight tRNA genes) are encoded on the light strand (L) and other genes on the heavy strand (H). Totally, The base composition of the mitochondrial genome included 27.29% A, 28.85% T, 28.87% C, and 15.00% G, respectively. The control regions among the Rana species were found to exhibit rich genetic variability and A + T content. R. omeimontis was clustered together with R. chaochiaoensis in phylogenetic tree. Compared to R. amurensis and R. kunyuensi, it was more closely related to R. chaochiaoensis, and a new way of gene rearrangement (ND6-trnE-Cytb-D-loop-trnL2 (CUN)-ND5-D-loop) was also found in the mitogenome of R. amurensis and R. kunyuensi. Our results about the mitochondrial genome of R. omeimontis will contribute to the future studies on phylogenetic relationship and the taxonomic status of Rana and related Ranidae species.

Complete mitochondrial genome of a blue-tailed skink Plestiodon capito (Reptilia, Squamata, Scincidae) and comparison with other Scincidae lizards

Vertebrate mitochondrial genomes (mitogenomes) are valuable for studying phylogeny, evolutionary genetics and genomics. To date, however, compared to other vertebrate groups, our knowledge about the mitogenomes of skinks (the family Scincidae), even of reptile, has been relatively limited. In the present study, we determined the complete mitogenome of a blue-tailed skink Plestiodon capito for the first time, and compared it with other skinks available in GenBank. The circular genome is 17,344 bp long, showing a typical vertebrate pattern with 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and one control region (CR). The gene organization, nucleotide composition, and codon usage are similar to those from skinks previously published. Twelve out of 13 PCGs initiates with canonical start codon (ATG), while COX1 starts with GTG. The codon usage analysis revealed a preferential use of the Leu^CUN (Leu1), Pro, and Thr codons with the A/U ending. All tRNAs in P. capito were predicted to fold into typical clover-leaf secondary structure, except tRNA-Ser ^AGY. The secondary structures of 12S rRNA and 16S rRNA comprises 34 helices and 56 helices, respectively. The alignment of the Plesitodon species CRs exhibited high genetic variability and rich A + T content. Besides, variable types and numbers of tandem repeat units were also identified in the CR of Plestiodon. Phylogenetic analyses recovered P. capito as the sister species to P. tunganus; monophyly of the Scincidae is well supported. Our results will help to better understand structure and evolution of the mitochondrial DNA control region in reptiles as well as the evolutionary status of P. capito, and to lay foundation for further phylogenetic study of skinks in a mitogenomic framework.

Structural Variation of the Turtle Mitochondrial Control Region

The present study describes the most comprehensive comparison of turtle mtD-loop regions to date. The primary structure was compared from DNA sequences accessed from GenBank from 48 species in 13 families of extant turtles, and secondary structures of the mtD-loop region were inferred from thermal stabilities, using the program Mfold, for each superfamiliy of turtles. Both primary and secondary structures were found to be highly variable across the order. The Cryptodira showed conservation in the primary structure at conserved sequence blocks (CSBs), but the Pleurodira displayed limited conservation of primary structural characters, other than the coreTAS, a binding site for the helicase TWINKLE, which was highly conserved in the Central and Right Domains across the order. No secondary structure was associated with a TAS, but an AT-rich fold (secondary structure) near the 3' terminus of the mtD-loop region was detected in all turtle superfamilies. Mapping of character states of structural features of the mtD-loop region revealed that most character states were autapomorphies and inferred a number of homoplasies. The Left Domain of turtles, containing no highly conserved structural elements, likely does not serve a functional role; therefore, the Central Domain in turtles is likely equivalent to the Left Domain of mammals. The AT-rich secondary structural element near the 3' terminus of the mtD-loop region may be conserved across turtles because of a functional role, perhaps containing the Light Strand Promotor, or perhaps interacting with the TWINKLE-coreTAS complex in the Central and Right Domains to regulate mtDNA replication and transcription.

A new species of the genus Lycodon (Serpentes, Colubridae) from Guangxi, China

A new species of colubrid snake, Lycodoncathayasp. nov., is described based on two adult male specimens collected from Huaping Nature Reserve, Guangxi, southern China. In a phylogenetic analyses, the new species is shown to be a sister taxon to the clade composed of L.futsingensis and L.namdongensis with low statistical support, and can be distinguished from all known congeners by the significant genetic divergence in the mitochondrial cytochrome b gene fragment (p-distance ≥ 7.9%), and morphologically by the following combination of characters: (1) dorsal scales in 17–17–15 rows, smooth throughout; (2) supralabials eight, third to fifth in contact with eye, infralabials nine; (3) ventral scales 199–200 (plus two preventral scales), subcaudals 78; (4) loreal single, elongated, in contact with eye or not, not in contact with internasals; (5) a single preocular not in contact with frontal, supraocular in contact with prefrontal, two postoculars; (6) maxillary teeth 10 (4+2+2+2); (7) two anterior temporals, three posterior temporals; (8) precloacal plate entire; (9) ground color from head to tail brownish black, with 31–35 dusty rose bands on body trunk, 13–16 on tail; (10) bands in 1–2 vertebral scales broad in minimum width; (11) bands separate ground color into brownish black ellipse patches arranged in a row along the top of body and tail; (12) elliptical patches in 3–6 scales of the vertebral row in maximum width; (13) ventral surface of body with wide brownish black strip, margined with a pair of continuous narrow greyish white ventrolateral lines. With the description of the new species, 64 congeners are currently known in the genus Lycodon, with 16 species occurring in China.

Six additional mitochondrial genomes for North American nightsnakes (Dipsadidae: Hypsiglena) and a novel gene feature for advanced snakes

The North American nightsnakes in the genus Hypsiglena is composed of nine named and at least two unnamed species. Here, we provide the first mt-genome of H. affinis, an additional mt-genome for H. sp. nov. 1, and four additional mt-genomes from the widespread H. jani. These mtDNA genomes were sequenced using both Illumina and Ion Torrent sequencing technologies. The resulting genomes contained the expected 13 protein coding genes, 22 tRNA genes, 2 rRNA genes, and 2 control regions typical of colubroid snakes. Two of the H. jani samples had partial tRNA^Ile genes upstream of CR2 which has not been previously documented in colubroid snakes. A maximum likelihood gene-tree based on these data combined with previously published sequence data recovers a well-supported phylogeny and is in concordance with previous estimates of evolutionary relationships in this group.

Resolving Phylogenetic Relationships within Passeriformes Based on Mitochondrial Genes and Inferring the Evolution of Their Mitogenomes in Terms of Duplications

Mitochondrial genes are placed on one molecule, which implies that they should carry consistent phylogenetic information. Following this advantage, we present a well-supported phylogeny based on mitochondrial genomes from almost 300 representatives of Passeriformes, the most numerous and differentiated Aves order. The analyses resolved the phylogenetic position of paraphyletic Basal and Transitional Oscines. Passerida occurred divided into two groups, one containing Paroidea and Sylvioidea, whereas the other, Passeroidea and Muscicapoidea. Analyses of mitogenomes showed four types of rearrangements including a duplicated control region (CR) with adjacent genes. Mapping the presence and absence of duplications onto the phylogenetic tree revealed that the duplication was the ancestral state for passerines and was maintained in early diverged lineages. Next, the duplication could be lost and occurred independently at least four times according to the most parsimonious scenario. In some lineages, two CR copies have been inherited from an ancient duplication and highly diverged, whereas in others, the second copy became similar to the first one due to concerted evolution. The second CR copies accumulated over twice as many substitutions as the first ones. However, the second CRs were not completely eliminated and were retained for a long time, which suggests that both regions can fulfill an important role in mitogenomes. Phylogenetic analyses based on CR sequences subjected to the complex evolution can produce tree topologies inconsistent with real evolutionary relationships between species. Passerines with two CRs showed a higher metabolic rate in relation to their body mass.

The revised complete mitogenome sequence of the tree frog Polypedatesmegacephalus (Anura, Rhacophoridae) by next-generation sequencing and phylogenetic analysis

The mitochondrial genome (mitogenome) sequence of the tree frog Polypedates megacephalus (16,473 bp) was previously reported as having the unusual characteristic of lacking the ND5 gene. In this study, a new mitogenome of P. megacephalus (19,952 bp) was resequenced using the next-generation sequencing (NGS) and standard Sanger sequencing technologies. It was discovered that the ND5 gene was not lost but translocated to the control region (CR) from its canonical location between the ND4 and ND6 genes. In addition, a duplicated control region was found in the new mitogenome of this species. Conservative region identification of the ND5 gene and phylogenetic analysis confirmed that the ND5 gene was located between two control regions. The phylogenetic relationship among 20 related species of anura revealed a rearrangement of the ND5 gene during the evolutionary process. These results also highlighted the advantages of next-generation sequencing. It will not only decrease the time and cost of sequencing, but also will eliminate the errors in published mitogenome databases.

Using high-resolution annotation of insect mitochondrial DNA to decipher tandem repeats in the control region

In this study, we used a small RNA sequencing (sRNA-seq) based method to annotate the mitochondrial genome of the insect Erthesina fullo Thunberg at 1 bp resolution. The high-resolution annotations cover both entire strands of the mitochondrial genome without any gaps or overlaps. Most of the new annotations were consistent with the previous annotations which had been obtained using PacBio full-length transcripts. Two important findings were that animals transcribe both entire strands of mitochondrial genomes and the tandem repeats in the control region of the E. fullo mitochondrial genome contains the repeated Transcription Initiation Sites (TISs) of the heavy strand. In addition, we found that the copy numbers of tandem repeats showed a great diversity within an individual, suggesting that mitochondrial DNA recombination occurs in an individual. In conclusion, the sRNA-seq based method uses 5' and 3' end small RNAs to annotate nuclear non-coding and mitochondrial genes at 1 bp resolution, and can be used to identify new steady RNAs, particularly long non-coding RNAs (lncRNAs). The high-resolution annotations of mitochondrial genomes can also be used to study the molecular phylogenetics and evolution of animals or to investigate mitochondrial gene transcription, RNA processing, RNA maturation and several other related topics. The complete mitochondrial genome sequence of E. fullo with the new annotations using the sRNA-seq based method is available at the NCBI GenBank database under the accession number MK374364. We publish our theories, methods, the high quality sRNA-seq and RNA-seq data (SRA: SRP174926) for extensive use.

New Insight into Parrots' Mitogenomes Indicates That Their Ancestor Contained a Duplicated Region

Mitochondrial genomes of vertebrates are generally thought to evolve under strong selection for size reduction and gene order conservation. Therefore, a growing number of mitogenomes with duplicated regions changes our view on the genome evolution. Among Aves, order Psittaciformes (parrots) is especially noteworthy because of its large morphological, ecological, and taxonomical diversity, which offers an opportunity to study genome evolution in various aspects. Former analyses showed that tandem duplications comprising the control region with adjacent genes are restricted to several lineages in which the duplication occurred independently. However, using an appropriate polymerase chain reaction strategy, we demonstrate that early diverged parrot groups contain mitogenomes with the duplicated region. These findings together with mapping duplication data from other mitogenomes onto parrot phylogeny indicate that the duplication was an ancestral state for Psittaciformes. The state was inherited by main parrot groups and was lost several times in some lineages. The duplicated regions were subjected to concerted evolution with a frequency higher than the rate of speciation. The duplicated control regions may provide a selective advantage due to a more efficient initiation of replication or transcription and a larger number of replicating genomes per organelle, which may lead to a more effective energy production by mitochondria. The mitogenomic duplications were associated with phenotypic features and parrots with the duplicated region can live longer, show larger body mass as well as predispositions to a more active flight. The results have wider implications on the presence of duplications and their evolution in mitogenomes of other avian groups.

Complete mitochondrial genome of Benthodytes marianensis (Holothuroidea: Elasipodida: Psychropotidae): Insight into deep sea adaptation in the sea cucumber

Complete mitochondrial genomes play important roles in studying genome evolution, phylogenetic relationships, and species identification. Sea cucumbers (Holothuroidea) are ecologically important and diverse members, living from the shallow waters to the hadal trench. In this study, we present the mitochondrial genome sequence of the sea cucumber Benthodytes marianensis collected from the Mariana Trench. To our knowledge, this is the first reported mitochondrial genome from the genus Benthodytes. This complete mitochondrial genome is 17567 bp in length and consists of 13 protein-coding genes, two ribosomal RNA genes and 22 transfer RNA genes (duplication of two tRNAs: trnL and trnS). Most of these genes are coded on the positive strand except for one protein-coding gene (nad6) and five tRNA genes which are coded on the negative strand. Two putative control regions (CRs) have been found in the B. marianensis mitogenome. We compared the order of genes from the 10 available holothurian mitogenomes and found a novel gene arrangement in B. marianensis. Phylogenetic analysis revealed that B. marianensis clustered with Peniagone sp. YYH-2013, forming the deep-sea Elasipodida clade. Positive selection analysis showed that eleven residues (24 S, 45 S, 185 S, 201 G, 211 F and 313 N in nad2; 108 S, 114 S, 322 C, 400 T and 427 S in nad4) were positively selected sites with high posterior probabilities. We predict that nad2 and nad4 may be the important candidate genes for the further investigation of the adaptation of B. marianensis to the deep-sea environment.

Resolving deep-sea pelagic saccopharyngiform eel mysteries: Identification of Neocyema and Monognathidae leptocephali and establishment of a new fish family "Neocyematidae" based on larvae, adults and mitogenomic gene orders

Deep-sea midwater "saccopharyngiform" eels of the families Cyematidae, Monognathidae, Eurypharyngidae and Saccopharyngidae (order Anguilliformes) are extraordinary fishes having major skeletal reductions and modifications compared to the general anguilliform body structure. Little is known about most aspects of the systematics, phylogeny, and ecology of these families, and few of the approximately 30 species described from adult specimens have been matched with their leptotocephalus larvae. Based on mitogenomic sequence data from rare new specimens, we show that the long-speculated-about larval form referred to as "Leptocephalus holti", which was thought to possibly be the larva of the rare orange-colored eels of Neocyema (5 known specimens; speculated to belong to the Cyematidae) are actually the larvae of the one-jaw eels of the family Monognathidae. One of the 5 types of L. holti larvae that were collected in the Pacific is genetically matched with Monognathus jesperseni, but multiple species exist based on larval sequence data and the morphology of adult specimens. A rare leptocephalus from the Sargasso Sea, with unique morphological characteristics including many small orange spots on the gut, was found to be the larva of Neocyema, which is presently only known from the Atlantic Ocean. We demonstrate that Neocyema constitutes a separate family being most closely related to Eurypharyngidae and Saccopharyngidae based on mitogenomic DNA sequences and unique mitochondrial gene orders.

Multiple independent structural dynamic events in the evolution of snake mitochondrial genomes

BACKGROUND

Mitochondrial DNA sequences have long been used in phylogenetic studies. However, little attention has been paid to the changes in gene arrangement patterns in the snake's mitogenome. Here, we analyzed the complete mitogenome sequences and structures of 65 snake species from 14 families and examined their structural patterns, organization and evolution. Our purpose was to further investigate the evolutionary implications and possible rearrangement mechanisms of the mitogenome within snakes.

RESULTS

In total, eleven types of mitochondrial gene arrangement patterns were detected (Type I, II, III, III-A, III-B, III-B1, III-C, III-D, III-E, III-F, III-G), with mitochondrial genome rearrangements being a major trend in snakes, especially in Alethinophidia. In snake mitogenomes, the rearrangements mainly involved three processes, gene loss, translocation and duplication. Within Scolecophidia, the O_L was lost several times in Typhlopidae and Leptotyphlopidae, but persisted as a plesiomorphy in the Alethinophidia. Duplication of the control region and translocation of the tRNA^Leu gene are two visible features in Alethinophidian mitochondrial genomes. Independently and stochastically, the duplication of pseudo-Pro (P*) emerged in seven different lineages of unequal size in three families, indicating that the presence of P* was a polytopic event in the mitogenome.

CONCLUSIONS

The WANCY tRNA gene cluster and the control regions and their adjacent segments were hotspots for mitogenome rearrangement. Maintenance of duplicate control regions may be the source for snake mitogenome structural diversity.

Monoplacophoran mitochondrial genomes: convergent gene arrangements and little phylogenetic signal

BACKGROUND

Although recent studies have greatly advanced understanding of deep molluscan phylogeny, placement of some taxa remains uncertain as different datasets support competing class-relationships. Traditionally, morphologists have placed Monoplacophora, a group of morphologically simple, limpet-like molluscs as sister group to all other conchiferans (shelled molluscs other than Polyplacophora), a grouping that is supported by the latest large-scale phylogenomic study that includes Laevipilina. However, molecular datasets dominated by nuclear ribosomal genes support Monoplacophora + Polyplacophora (Serialia). Here, we evaluate the potential of mitochondrial genome data for resolving placement of Monoplacophora.

RESULTS

Two complete (Laevipilina antarctica and Vema ewingi) and one partial (Laevipilina hyalina) mitochondrial genomes were sequenced, assembled, and compared. All three genomes show a highly similar architecture including an unusually high number of non-coding regions. Comparison of monoplacophoran gene order shows a gene arrangement pattern not previously reported; there is an inversion of one large gene cluster. Our reanalyses of recently published polyplacophoran mitogenomes show, however, that this feature is also present in some chiton species. Maximum Likelihood and Bayesian Inference analyses of 13 mitochondrial protein-coding genes failed to robustly place Monoplacophora and hypothesis testing could not reject any of the evaluated placements of Monoplacophora.

CONCLUSIONS

Under both serialian or aculiferan-conchiferan scenarios, the observed gene cluster inversion appears to be a convergent evolution of gene arrangements in molluscs. Our phylogenetic results are inconclusive and sensitive to taxon sampling. Aculifera (Polyplacophora + Aplacophora) and Conchifera were never recovered. However, some analyses recovered Serialia (Monoplacophora + Polyplacophora), Diasoma (Bivalvia + Scaphopoda) or Pleistomollusca (Bivalvia + Gastropoda). Although we could not shed light on deep evolutionary traits of Mollusca we found unique patterns of gene arrangements that are common to monoplacophoran and chitonine polyplacophoran species but not to acanthochitonine Polyplacophora. Complete mitochondrial genome of Laevipilina antarctica.

Mitochondrial DNA Sequence and Lack of Response to Anoxia in the Annual Killifish Austrofundulus limnaeus

The annual killifish Austrofundulus limnaeus inhabits ephemeral ponds in regions of Venezuela, South America. Permanent populations of A. limnaeus are maintained by production of stress-tolerant embryos that are able to persist in the desiccated sediment. Previous work has demonstrated that A. limnaeus have a remarkable ability to tolerate extended periods of anoxia and desiccating conditions. After considering temperature, A. limnaeus embryos have the highest known tolerance to anoxia when compared to any other vertebrate yet studied. Oxygen is completely essential for the process of oxidative phosphorylation by mitochondria, the intracellular organelle responsible for the majority of adenosine triphosphate production. Thus, understanding the unique properties of A. limnaeus mitochondria is of great interest. In this work, we describe the first complete mitochondrial genome (mtgenome) sequence of a single adult A. limnaeus individual and compare both coding and non-coding regions to several other closely related fish mtgenomes. Mitochondrial features were predicted using MitoAnnotator and polyadenylation sites were predicted using RNAseq mapping. To estimate the responsiveness of A. limnaeus mitochondria to anoxia treatment, we measure relative mitochondrial DNA copy number and total citrate synthase activity in both relatively anoxia-tolerant and anoxia-sensitive embryonic stages. Our cross-species comparative approach identifies unique features of ND1, ND5, ND6, and ATPase-6 that may facilitate the unique phenotype of A. limnaeus embryos. Additionally, we do not find evidence for mitochondrial degradation or biogenesis during anoxia/reoxygenation treatment in A. limnaeus embryos, suggesting that anoxia-tolerant mitochondria do not respond to anoxia in a manner similar to anoxia-sensitive mitochondria.

Novel duplication pattern of the mitochondrial control region in Cantor's Giant softshell turtle Pelochelys cantorii

Cantor's Giant Softshell Turtle, Pelochelys cantorii has become one of the most critically endangered species in the world. When comparative analyses of the P. cantorii complete mitochondrial genome sequences were conducted, we discovered a duplication of a segment of the control region in the mitochondrial genome of P. cantorii. The duplication is characterized by two copies of conserved sequence box 2 (CSB2) and CSB3 in a single control region. In contrast to previous reports of duplications involving the control regions of other animals, this particular pattern of duplications appears to be unique to P. cantorii. Copies of the CSB2 and CSB3 show many of the conserved sequence features typically found in mitochondrial control regions, and rare differences were found between the paralogous copies. Using the primer design principle of simple sequence repeats (SSR) and the reference sequence of the duplicated CSBs, specific primers were designed to amplify the duplicated CSBs. These primers were validated among different individuals and populations of P. cantorii. This unique duplication structure suggests the two copies of the CSB2 and CSB3 may have arisen through occasional tandem duplication and subsequent concerted evolution.

Rearrangement and evolution of mitochondrial genomes in parrots

Mitochondrial genome rearrangements that result in control region duplication have been described for a variety of birds, but the mechanisms leading to their appearance and maintenance remain unclear, and their effect on sequence evolution has not been explored. A recent survey of mitochondrial genomes in the Psittaciformes (parrots) found that control region duplications have arisen independently at least six times across the order. We analyzed complete mitochondrial genome sequences from 20 parrot species, including representatives of each lineage with control region duplications, to document the gene order changes and to examine effects of genome rearrangements on patterns of sequence evolution. The gene order previously reported for Amazona parrots was found for four of the six independently derived genome rearrangements, and a previously undescribed gene order was found in Prioniturus luconensis, representing a fifth clade with rearranged genomes; the gene order resulting from the remaining rearrangement event could not be confirmed. In all rearranged genomes, two copies of the control region are present and are very similar at the sequence level, while duplicates of the other genes involved in the rearrangement show signs of degeneration or have been lost altogether. We compared rates of sequence evolution in genomes with and without control region duplications and did not find a consistent acceleration or deceleration associated with the duplications. This could be due to the fact that most of the genome rearrangement events in parrots are ancient, and additionally, to an effect of body size on evolutionary rate that we found for mitochondrial but not nuclear sequences. Base composition analyses found that relative to other birds, parrots have unusually strong compositional asymmetry (AT- and GC-skew) in their coding sequences, especially at fourfold degenerate sites. Furthermore, we found higher AT skew in species with control region duplications. One potential cause for this compositional asymmetry is that parrots have unusually slow mtDNA replication. If this is the case, then any replicative advantage provided by having a second control region could result in selection for maintenance of both control regions once duplicated.

Concerted Evolution of Duplicate Control Regions in the Mitochondria of Species of the Flatfish Family Bothidae (Teleostei: Pleuronectiformes)

Mitogenomes of flatfishes (Pleuronectiformes) exhibit the greatest diversity of gene rear-rangements in teleostean fishes. Duplicate control regions (CRs) have been found in the mito-genomes of two flatfishes, Samariscus latus (Samaridae) and Laeops lanceolata (Bothidae), which is rare in teleosts. It has been reported that duplicate CRs have evolved in a concerted fashion in fishes and other animals, however, whether concerted evo-lution exists in flatfishes remains unknown. In this study, based on five newly sequenced and six previously reported mitogenomes of lefteye flounders in the Bothidae, we explored whether duplicate CRs and concerted evolution exist in these species. Results based on the present study and previous reports show that four out of eleven bothid species examined have duplicate CRs of their mitogenomes. The core regions of the duplicate CRs of mitogenomes in the same species have identical, or nearly identical, sequences when compared to each other. This pattern fits the typical characteristics of concerted evolution. Additionally, phylogenetic and ancestral state reconstruction analysis also provided evidence to support the hypothesis that duplicate CRs evolved concertedly. The core region of concerted evolution is situated at the conserved domains of the CR of the mitogenome from the termination associated sequences (TASs) to the conserved sequence blocks (CSBs). Commonly, this region is con-sidered to regulate mitochondrial replication and transcription. Thus, we hypothesize that the cause of concerted evolution of the duplicate CRs in the mtDNAs of these four bothids may be related to some function of the conserved sequences of the CRs during mitochondrial rep-lication and transcription. We hope our results will provide fresh insight into the molecular mechanisms related to replication and evolution of mitogenomes.

Complete mitochondrial genome of Epigonichthys cultellus (Cephalochordata: Branchiostomatidae)

The systematic position of the amphioxus species with dextral gonads distributed in the southeast Pacific and Indian Oceans has remained to be clarified due to the adoption of different names by different researchers. Mitochondrial (mt) DNA is generally considered to be a powerful molecular marker in taxonomic studies. For a reliable systematics of the amphioxi collected from the South China Sea, we sequenced the complete mtDNA from a single specimen and compared it with those of the other eight amphioxus species. The present mtDNA genome contains 13 protein-coding genes, 22 tRNA genes, and two rRNA genes, with the same gene order as those in Branchiostoma and Epigonichthys, which, however, is different from that in Asymmetron. Based on our morphological data (including measurements of some characters) and the features of the mt genome, together with the distribution records of the dextral-gonad amphioxi, we conclude that the present mtDNA sequence most likely represents that of E. cultellus.

The first complete mitochondrial genome of Pygopodidae (Aprasia parapulchella Kluge)

The Pygopodidae comprise an enigmatic group of legless lizards endemic to the Australo-Papuan region. Here we present the first complete mitochondrial genome for a member of this family, Aprasia parapulchella, from Australia. The mitochondrial genome of A. parapulchella is 16 528 base pairs long and contains 13 protein-coding genes, 22 tRNA genes, two rRNA genes and the control region, conforming to the typical vertebrate gene order. The overall mitochondrial nucleotide composition is 31.7% A, 24.5% T, 30.5% C and 13.2% G. This corresponds to a total A+T content of 56.3%, which is similar to that of other squamate lizard genomes.

An evolutionary preserved intergenic spacer in gadiform mitogenomes generates a long noncoding RNA

Background

Vertebrate mitogenomes are economically organized and usually lack intergenic sequences other than the control region. Intergenic spacers located between the tRNA^Thr and tRNA^Pro genes (“T-P spacers”) have been observed in several taxa, including gadiform species, but information about their biological roles and putative functions is still lacking.

Results

Sequence characterization of the complete European hake Merluccius merluccius mitogenome identified a complex T-P spacer ranging in size from 223–532 bp. Further analyses of 32 gadiform species, representing 8 families and 28 genera, revealed the evolutionary preserved presence of T-P spacers across all taxa. Molecular complexity of the T-P spacers was found to be coherent with the phylogenetic relationships, supporting a common ancestral origin and gain of function during codfish evolution. Intraspecific variation of T-P spacer sequences was assessed in 225 Atlantic cod specimens and revealed 26 haplotypes. Pyrosequencing data representing the mito-transcriptome poly (A) fraction in Atlantic cod identified an abundant H-strand specific long noncoding RNA of about 375 nt. The T-P spacer corresponded to the 5’ part of this transcript, which terminated within the control region in a tail-to-tail configuration with the L-strand specific transcript (the 7S RNA).

Conclusions

The T-P spacer is inferred to be evolutionary preserved in gadiform mitogenomes due to gain of function through a long noncoding RNA. We suggest that the T-P spacer adds stability to the H-strand specific long noncoding RNA by forming stable hairpin structures and additional protein binding sites.

The complete mitochondrial genomes of sixteen ardeid birds revealing the evolutionary process of the gene rearrangements

BACKGROUND

The animal mitochondrial genome is generally considered to be under selection for both compactness and gene order conservation. As more mitochondrial genomes are sequenced, mitochondrial duplications and gene rearrangements have been frequently identified among diverse animal groups. Although several mechanisms of gene rearrangement have been proposed thus far, more observational evidence from major taxa is needed to validate specific mechanisms. In the current study, the complete mitochondrial DNA of sixteen bird species from the family Ardeidae was sequenced and the evolution of mitochondrial gene rearrangements was investigated. The mitochondrial genomes were then used to review the phylogenies of these ardeid birds.

RESULTS

The complete mitochondrial genome sequences of the sixteen ardeid birds exhibited four distinct mitochondrial gene orders in which two of them, named as "duplicate tRNA(Glu)-CR" and "duplicate tRNAThr-tRNA(Pro) and CR", were newly discovered. These gene rearrangements arose from an evolutionary process consistent with the tandem duplication--random loss model (TDRL). Additionally, duplications in these gene orders were near identical in nucleotide sequences within each individual, suggesting that they evolved in concert. Phylogenetic analyses of the sixteen ardeid species supported the idea that Ardea ibis, Ardea modesta and Ardea intermedia should be classified as genus Ardea, and Ixobrychus flavicollis as genus Ixobrychus, and indicated that within the subfamily Ardeinae, Nycticorax nycticorax is closely related to genus Egretta and that Ardeola bacchus and Butorides striatus are closely related to the genus Ardea.

CONCLUSIONS

The duplicate tRNAThr-CR gene order is found in most ardeid lineages, suggesting this gene order is the ancestral pattern within these birds and persisted in most lineages via concerted evolution. In two independent lineages, when the concerted evolution stopped in some subsections due to the accumulation of numerous substitutions and deletions, the duplicate tRNAThr-CR gene order was transformed into three other gene orders. The phylogenetic trees produced from concatenated rRNA and protein coding genes have high support values in most nodes, indicating that the mitochondrial genome sequences are promising markers for resolving the phylogenetic issues of ardeid birds when more taxa are added.

The complete mitochondrial genome of the striped-tailed rat-snake, Orthriophis taeniurus (Reptilia, Serpentes, Colubridae)

The complete mitochondrial genome (mitogenome) of the striped-tailed rat-snake Orthriophis taeniurus was determined in the present study. The genome is 17,183 bp in size, containing 2 ribosomal RNA (rRNA) genes, 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, and 2 control regions (CRI and CRII). The gene order and orientation in O. taeniurus mitogenome are basically identical to that of other alethinophidian snakes. Nucleotide composition is very similar with other vertebrates, showing an AT bias.

Structure and evolution of the Phasianidae mitochondrial DNA control region

The mitochondrial DNA control region is an area of the mitochondrial genome which is non-coding DNA. To infer the structural and evolutionary characteristics of Phasianidae mitochondrial DNA control region, the entire control region sequences of 34 species were analyzed. The length of the control region sequences ranged from 1144 bp (Phasianus colchicus) to 1555 bp (Coturnix japonica) and can be separated into three domains. The average genetic distances among the species within the genera varied from 1.96% (Chrysolophus) to 12.05% (Coturnix). The average genetic distances showed significantly negative correlation with ts/tv. In most genera (except Coturnix), domain I is the most variable among the three domains. However, the first 150 nucleotides apparently evolved at unusually low rates. Four conserved sequence boxes in the domain II of Phasianidae sequences were identified. The alignment of the Phasianidae four boxes and CSB-1 sequences showed considerable sequence variation.

Recombination and evolution of duplicate control regions in the mitochondrial genome of the Asian big-headed turtle, Platysternon megacephalum

Complete mitochondrial (mt) genome sequences with duplicate control regions (CRs) have been detected in various animal species. In Testudines, duplicate mtCRs have been reported in the mtDNA of the Asian big-headed turtle, Platysternon megacephalum, which has three living subspecies. However, the evolutionary pattern of these CRs remains unclear. In this study, we report the completed sequences of duplicate CRs from 20 individuals belonging to three subspecies of this turtle and discuss the micro-evolutionary analysis of the evolution of duplicate CRs. Genetic distances calculated with MEGA 4.1 using the complete duplicate CR sequences revealed that within turtle subspecies, genetic distances between orthologous copies from different individuals were 0.63% for CR1 and 1.2% for CR2app:addword:respectively, and the average distance between paralogous copies of CR1 and CR2 was 4.8%. Phylogenetic relationships were reconstructed from the CR sequences, excluding the variable number of tandem repeats (VNTRs) at the 3' end using three methods: neighbor-joining, maximum likelihood algorithm, and Bayesian inference. These data show that any two CRs within individuals were more genetically distant from orthologous genes in different individuals within the same subspecies. This suggests independent evolution of the two mtCRs within each P. megacephalum subspecies. Reconstruction of separate phylogenetic trees using different CR components (TAS, CD, CSB, and VNTRs) suggested the role of recombination in the evolution of duplicate CRs. Consequently, recombination events were detected using RDP software with break points at ≈290 bp and ≈1,080 bp. Based on these results, we hypothesize that duplicate CRs in P. megacephalum originated from heterological ancestral recombination of mtDNA. Subsequent recombination could have resulted in homogenization during independent evolutionary events, thus maintaining the functions of duplicate CRs in the mtDNA of P. megacephalum.

The largest unassigned regions of the male- and female-transmitted mitochondrial DNAs in Musculista senhousia (Bivalvia Mytilidae)

Musculista senhousia is a marine mussel with doubly uniparental inheritance (DUI) of mitochondria. In this study we analyzed the largest unassigned region (LUR) of its female- and male-transmitted mitochondrial genomes, described their fine characteristics and searched for shared features. Our results suggest that both LURs contain the control region of their respective mitochondrial genomes. The female-transmitted control region is duplicated in tandem, with the two copies evolving in concert. This makes the F-mtDNA of M. senhousia the first Bivalve mitochondrial genome with this feature. We also compared M. senhousia control regions to that of other Mytilidae, and demonstrated that signals for basic mtDNA functions are retained over evolutionary times even among the fast-evolving mitochondrial genomes of DUI species. Finally, we discussed how similarities between female and male LURs may be explained in the context of DUI evolution and if the duplicated female control region might have influenced the DUI system in this species.

Afrobatrachian mitochondrial genomes: genome reorganization, gene rearrangement mechanisms, and evolutionary trends of duplicated and rearranged genes

Background

Mitochondrial genomic (mitogenomic) reorganizations are rarely found in closely-related animals, yet drastic reorganizations have been found in the Ranoides frogs. The phylogenetic relationships of the three major ranoid taxa (Natatanura, Microhylidae, and Afrobatrachia) have been problematic, and mitogenomic information for afrobatrachians has not been available. Several molecular models for mitochondrial (mt) gene rearrangements have been proposed, but observational evidence has been insufficient to evaluate them. Furthermore, evolutionary trends in rearranged mt genes have not been well understood. To gain molecular and phylogenetic insights into these issues, we analyzed the mt genomes of four afrobatrachian species (Breviceps adspersus, Hemisus marmoratus, Hyperolius marmoratus, and Trichobatrachus robustus) and performed molecular phylogenetic analyses. Furthermore we searched for two evolutionary patterns expected in the rearranged mt genes of ranoids.

Results

Extensively reorganized mt genomes having many duplicated and rearranged genes were found in three of the four afrobatrachians analyzed. In fact, Breviceps has the largest known mt genome among vertebrates. Although the kinds of duplicated and rearranged genes differed among these species, a remarkable gene rearrangement pattern of non-tandemly copied genes situated within tandemly-copied regions was commonly found. Furthermore, the existence of concerted evolution was observed between non-neighboring copies of triplicated 12S and 16S ribosomal RNA regions.

Conclusions

Phylogenetic analyses based on mitogenomic data support a close relationship between Afrobatrachia and Microhylidae, with their estimated divergence 100 million years ago consistent with present-day endemism of afrobatrachians on the African continent. The afrobatrachian mt data supported the first tandem and second non-tandem duplication model for mt gene rearrangements and the recombination-based model for concerted evolution of duplicated mt regions. We also showed that specific nucleotide substitution and compositional patterns expected in duplicated and rearranged mt genes did not occur, suggesting no disadvantage in employing these genes for phylogenetic inference.

Mitochondrial genome sequence of Sibynophis chinensis (Squamata, Colubridae)

The complete mitochondrial (mt) genome of the Chinese many-toothed snake, Sibynophis chinensis, was sequenced and found to be 17,163 bp in length. The arrangement of 13 protein-coding genes, tRNAs and rRNAs was identical to that of other common snake mt genomes. The mt protein-coding genes of S. chinensis utilized ATA, ATG, ATA and GTG as initiation codons and AGA, AGG, TAA, TAG and T as termination codons. Among three tRNA clusters (LQM, WANCY and HSL), LQM was found instead of IQM, which is common in other vertebrates. We also identified two control regions that contained several conserved elements known as conserved sequence blocks and termination-associated sequences related to mt replication and transcription.

The complete mitochondrial genome of Oocatochus rufodorsatus (Reptilia, Serpentes, Colubridae)

The complete mitochondrial genome from the red-backed ratsnake Oocatochus rufodorsatus was determined. The following four distinctive features are observed: duplicate control regions that have nearly identical sequences at two different locations of mitogenome, translocation of tRNA-Leu (UUR), a pseudogene for tRNA-Pro and truncations of TΨC arm for most tRNA genes.

The complete mitochondrial genomes of deep-sea squid (Bathyteuthis abyssicola), bob-tail squid (Semirossia patagonica) and four giant cuttlefish (Sepia apama, S. latimanus, S. lycidas and S. pharaonis), and their application to the phylogenetic analysis of Decapodiformes

We determined the complete mitochondrial (mt) genomes of the deep-sea squid (Bathyteuthis abyssicola; supperfamily Bathyteuthoidea), the bob-tail squid (Semirossia patagonica; order Sepiolida) and four giant cuttlefish (Sepia apama, S. latimanus, S. lycidas and S. pharaonis; order Sepiida). The unique structures of the mt genomes of Bathyteuthis and Semirossia provide new information about the evolution of decapodiform mt genomes. We show that the mt genome of B. abyssicola, like those of other oegopsids studied so far, has two long duplicated regions that include seven genes (COX1-3, ATP6 and ATP8, tRNA(Asn), and either ND2 or ND3) and that one of the duplicated COX3 genes has lost its function. The mt genome of S. patagonica is unlike any other decapodiforms and, like Nautilus, its ATP6 and ATP8 genes are not adjacent to each other. The four giant cuttlefish have identical mt gene order to other cuttlefish determined to date. Molecular phylogenetic analyses using maximum likelihood and Bayesian methods suggest that traditional order Sepioidea (Sepiolida+Sepiida) is paraphyletic and Sepia (cuttlefish) has the sister-relationship with all other decapodiforms. Taking both the phylogenetic analyses and the mt gene order analyses into account, it is likely that the octopus-type mt genome is an ancestral state and that it had maintained from at least the Cephalopoda ancestor to the common ancestor of Oegopsida, Myopsida and Sepiolida.

Complete mitochondrial genome of the sea slug, Placida sp., contains unusual noncoding regions

The complete sequence (14,751 bp) of the mitochondrial DNA (mtDNA) of the opisthobranch gastropod Placida sp. was determined using long PCR and genome-walking techniques. The genome contains 13 protein-coding genes (PCG), 2 ribosomal RNA (rRNA), 22 transfer RNA (tRNA) genes that are typical of metazoan mtDNA, and 2 lengthy noncoding regions with a total length of 1441 bp. All the tRNA genes have general secondary structures without reduced T or D stems. The mitochondrial gene arrangement of Placida sp. is almost identical to that of sea slugs, whereas the tRNA gene arrangement is different from other animals outside Sacoglossa in the Euthyneura. A phylogenetic analysis based on the mitochondrial genome of Euthyneura indicated that Placida sp. and Placida dendritica are closely related to Elysia chlorotica as a sister taxon, which strongly supports their close phylogenetic affinity.

Evolution of extensively fragmented mitochondrial genomes in the lice of humans

Bilateral animals are featured by an extremely compact mitochondrial (mt) genome with 37 genes on a single circular chromosome. The human body louse, Pediculus humanus, however, has its mt genes on 20 minichromosomes. We sequenced the mt genomes of two other human lice: the head louse, P. capitis, and the pubic louse, Pthirus pubis. Comparison among the three human lice revealed the presence of fragmented mt genomes in their most recent common ancestor, which lived ∼7 Ma. The head louse has exactly the same set of mt minichromosomes as the body louse, indicating that the number of minichromosomes, and the gene content and gene arrangement in each minichromosome have remained unchanged since the body louse evolved from the head louse ∼107,000 years ago. The pubic louse has the same pattern of one protein-coding or rRNA gene per minichromosome (except one minichromosome with two protein-coding genes, atp6 and atp8) as the head louse and the body louse. This pattern is apparently ancestral to all human lice and has been stable for at least 7 Myr. Most tRNA genes of the pubic louse, however, are on different minichromosomes when compared with their counterparts in the head louse and the body louse. It is evident that rearrangement of four tRNA genes (for leucine, arginine and glycine) was due to gene-identity switch by point mutation at the third anticodon position or by homologous recombination, whereas rearrangement of other tRNA genes was by gene translocation between minichromosomes, likely caused by minichromosome split via gene degeneration and deletion.

Low metabolic rates in salamanders are correlated with weak selective constraints on mitochondrial genes

Mitochondria are the site for the citric acid cycle and oxidative phosphorylation (OXPHOS), the final steps of ATP synthesis via cellular respiration. Each mitochondrion contains its own genome; in vertebrates, this is a small, circular DNA molecule that encodes 13 subunits of the multiprotein OXPHOS electron transport complexes. Vertebrate lineages vary dramatically in metabolic rates; thus, functional constraints on mitochondrial-encoded proteins likely differ, potentially impacting mitochondrial genome evolution. Here, we examine mitochondrial genome evolution in salamanders, which have the lowest metabolic requirements among tetrapods. We show that salamanders experience weaker purifying selection on protein-coding sequences than do frogs, a comparable amphibian clade with higher metabolic rates. In contrast, we find no evidence for weaker selection against mitochondrial genome expansion in salamanders. Together, these results suggest that different aspects of mitochondrial genome evolution (i.e., nucleotide substitution, accumulation of noncoding sequences) are differently affected by metabolic variation across tetrapod lineages.