Multiple independent origins of mitochondrial gene order in birds

The mitochondrial genome of Bottapotamon fukienense (Brachiura: Potamidae) is fragmented into two chromosomes

BACKGROUND

China is the hotspot of global freshwater crab diversity, but their wild populations are facing severe pressures associated with anthropogenic factors, necessitating the need to map their taxonomic and genetic diversity and design conservation policies.

RESULTS

Herein, we sequenced the mitochondrial genome of a Chinese freshwater crab species Bottapotamon fukienense, and found that it is fragmented into two chromosomes. We confirmed that fragmentation was not limited to a single specimen or population. Chromosome 1 comprised 15,111 base pairs (bp) and there were 26 genes and one pseudogene (pseudo-nad1) encoded on it. Chromosome 2 comprised 8,173 bp and there were 12 genes and two pseudogenes (pseudo-trnL2 and pseudo-rrnL) encoded on it. Combined, they comprise the largest mitogenome (23,284 bp) among the Potamidae. Bottapotamon was the only genus in the Potamidae dataset exhibiting rearrangements of protein-coding genes. Bottapotamon fukienense exhibited average rates of sequence evolution in the dataset and did not differ in selection pressures from the remaining Potamidae.

CONCLUSIONS

This is the first experimentally confirmed fragmentation of a mitogenome in crustaceans. While the mitogenome of B. fukienense exhibited multiple signs of elevated mitogenomic architecture evolution rates, including the exceptionally large size, duplicated genes, pseudogenisation, rearrangements of protein-coding genes, and fragmentation, there is no evidence that this is matched by elevated sequence evolutionary rates or changes in selection pressures.

Phylomitogenomics elucidates the evolution of symbiosis in Thoracotremata (Decapoda: Cryptochiridae, Pinnotheridae, Varunidae)

Background

Thoracotremata belong to the large group of "true" crabs (infraorder Brachyura), and they exhibit a wide range of physiological and morphological adaptations to living in terrestrial, freshwater and marine habitats. Moreover, the clade comprises various symbiotic taxa (Aphanodactylidae, Cryptochiridae, Pinnotheridae, some Varunidae) that are specialised in living with invertebrate hosts, but the evolutionary history of these symbiotic crabs is still partially unresolved.

Methods

Here we assembled and characterised the complete mitochondrial genomes (hereafter mitogenomes) of three gall crab species (Cryptochiridae): Kroppcarcinus siderastreicola, Opecarcinus hypostegus and Troglocarcinus corallicola. A phylogenetic tree of the Thoracotremata was reconstructed using 13 protein-coding genes and two ribosomal RNA genes retrieved from three new gall crab mitogenomes and a further 72 available thoracotreme mitogenomes. Furthermore, we applied a comparative analysis to characterise mitochondrial gene order arrangement, and performed a selection analysis to test for selective pressure of the protein-coding genes in symbiotic Cryptochiridae, Pinnotheridae, and Varunidae (Asthenognathus inaequipes and Tritodynamia horvathi).

Results

The results of the phylogenetic reconstruction confirm the monophyly of Cryptochiridae, which clustered separately from the Pinnotheridae. The latter clustered at the base of the tree with robust branch values. The symbiotic varunids A. inaequipes and T. horvathi clustered together in a clade with free-living Varunidae species, highlighting that symbiosis in the Thoracotremata evolved independently on multiple occasions. Different gene orders were detected in symbionts and free-living species when compared with the ancestral brachyuran gene order. Lastly, the selective pressure analysis detected two positively selected sites in the nad6 gene of Cryptochiridae, but the evidence for positive selection in Pinnotheridae and A. inaequipes and T. horvathi was weak. Adaptive evolution of mitochondrial protein-coding genes is perhaps related to the presumably higher energetic demands of a symbiotic lifestyle.

Identification and phylogenetic analysis in Pterorhinuschinensis (Aves, Passeriformes, Leiothrichidae) based on complete mitogenome

The Black-throated Laughingthrush (Pterorhinuschinensis) is a bird belonging to the order Passeriformes and the family Leiothrichidae, and is found in Cambodia, China, Laos, Myanmar, Thailand and Vietnam. Pterorhinuschinensis was once classified as belonging to the genus Garrulax. However, recent research has reclassified it in the genus Pterorhinus. In this study, we sequenced and characterized the complete mitogenome of P.chinensis. The complete mitochondrial genome of P.chinensis is 17,827 bp in length. It consists of 13 PCGs, 22 tRNAs, two rRNAs, and two control regions. All genes are coded on the H-strand, except for one PCG (nad6) and eight tRNAs. All PCGs are initiated with ATG and stopped by five types of stop codons. Our comparative analyses show irregular gene rearrangement between trnT and trnP genes with another similar control region emerging between trnE and trnF genes compared with the ancestral mitochondrial gene order, called "duplicate CR gene order". The phylogenetic position of P.chinensis and phylogenetic relationships among members of Leiothrichidae are assessed based on complete mitogenomes. Phylogenetic relationships based on Bayesian inference and maximum likelihood methods showed that Garrulax and (Pterorhinus + Ianthocincla) formed a clade. Leiothrix and Liocichla also formed a clade. Our study provides support for the transfer of P.chinensis from Garrulax to Pterorhinus. Our results provide mitochondrial genome data to further understand the mitochondrial genome characteristics and taxonomic status of Leiothrichidae.

Characterization and Phylogenetic Analysis of the Complete Mitochondrial Genome of Aythya marila

Aythya marila is a large diving duck belonging to the family Anatidae. However, the phylogenetic relationship among these Aythya species remains unclear due to the presence of extensive interspecific hybridization events within the Aythya genus. Here, we sequenced and annotated the complete mitochondrial genome of A. marila, which contained 22 tRNAs, 13 protein-coding genes (PCGs), 2 ribosomal RNAs, and 1 D-loop, with a length of 16,617 bp. The sizes of the PCGs ranged from 297 to 1824 bp and were all, except for ND6, located on the heavy chain (H). ATG and TAA were the most common start and termination codons of the 13 PCGs, respectively. The fastest- and slowest-evolving genes were ATP8 and COI, respectively. Codon usage analysis indicated that CUA, AUC, GCC, UUC, CUC, and ACC were the six most frequent codons. The nucleotide diversity values indicated a high level of genetic diversity in A. marila. F_ST analysis suggested a widespread gene exchange between A. baeri and A. nyroca. Moreover, phylogenetic reconstructions using the mitochondrial genomes of all available Anatidae species showed that, in addition to A. marila, four major clades among the Anatidae (Dendrocygninae, Oxyurinae, Anserinae, and Anatinae) were closely related to A. fuligula. Overall, this study provides valuable information on the evolution of A. marila and new insights into the phylogeny of Anatidae.

Rapid evolution, rearrangements and whole mitogenome duplication in the Australian stingless bees Tetragonula (Hymenoptera: Apidae): A steppingstone towards understanding mitochondrial function and evolution

The extreme conservation of mitochondrial genomes in metazoans poses a significant challenge to understanding mitogenome evolution. However, the presence of variation in gene order or genome structure, found in a small number of taxa, can provide unique insights into this evolution. Previous work on two stingless bees in the genus Tetragonula (T. carbonaria and T. hockingsi) revealed highly divergent CO1 regions between them and when compared to the bees from the same tribe (Meliponini), indicating rapid evolution. Using mtDNA isolation and Illumina sequencing, we elucidated the mitogenomes of both species. In both species, there has been a duplication of the whole mitogenome to give a total genome size of 30,666 bp in T. carbonaria; and 30,662 bp in T. hockingsi. These duplicated genomes present a circular structure with two identical and mirrored copies of all 13 protein coding genes and 22 tRNAs, with the exception of a few tRNAs that are present as single copies. In addition, the mitogenomes are characterized by rearrangements of two block of genes. We believe that rapid evolution is present in the whole Indo-Malay/Australasian group of Meliponini but is extraordinarily elevated in T. carbonaria and T. hockingsi, probably due to founder effect, low effective population size and the mitogenome duplication. All these features - rapid evolution, rearrangements, and duplication - deviate significantly from the vast majority of the mitogenomes described so far, making the mitogenomes of Tetragonula unique opportunities to address fundamental questions of mitogenome function and evolution.

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.

Mitogenomic Codon Usage Patterns of Superfamily Certhioidea (Aves, Passeriformes): Insights into Asymmetrical Bias and Phylogenetic Implications

The superfamily Certhioidea currently comprises five families. Due to the rapid diversification, the phylogeny of Certhioidea is still controversial. The advent of next generation sequencing provides a unique opportunity for a mitogenome-wide study. Here, we first provided six new complete mitogenomes of Certhioidea (Certhia americana, C. familiaris, Salpornis spilonota, Cantorchilus leucotis, Pheugopedius coraya, and Pheugopedius genibarbis). We further paid attention to the genomic characteristics, codon usages, evolutionary rates, and phylogeny of the Certhioidea mitogenomes. All mitogenomes we analyzed displayed typical ancestral avian gene order with 13 protein-coding genes (PCGs), 22 tRNAs, 2 rRNAs, and one control region (CR). Our study indicated the strand-biased compositional asymmetry might shape codon usage preferences in mitochondrial genes. In addition, natural selection might be the main factor in shaping the codon usages of genes. Additionally, evolutionary rate analyses indicated all mitochondrial genes were under purifying selection. Moreover, MT-ATP8 and MT-CO1 were the most rapidly evolving gene and conserved genes, respectively. According to our mitophylogenetic analyses, the monophylies of Troglodytidae and Sittidae were strongly supported. Importantly, we suggest that Salpornis should be separated from Certhiidae and put into Salpornithidae to maintain the monophyly of Certhiidae. Our findings are useful for further evolutionary studies within Certhioidea.

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.

Caryophylliids (Anthozoa, Scleractinia) and mitochondrial gene order: insights from mitochondrial and nuclear phylogenomics

Molecularly, the family Caryophylliidae is polyphyletic and different sets of genetic data converge towards a consensus that a taxonomic review of this family is necessary. Overall, the order of genes in the mitochondrial genome (mitogenome) together with DNA sequences have been used to successfully untangle evolutionary relationships in several groups of organisms. Published mitogenomes of two caryophylliid genera (Desmophyllum and Solenosmilia) present a transposition of the gene block containing cob, nad2, and nad6, which is located between nad5 5' exon and trnW, while that of Polycyathus chaishanensis presents the same gene order as the majority of scleractinian corals. In molecular-based evolutionary reconstructions, caryophylliids that have the mitochondrial gene rearrangement were recovered as a monophyletic lineage ("true" caryophylliids), while members of the genus Polycyathus were placed in a different position. In this study, additional mitogenomes of this family were assembled and included in evolutionary reconstructions of Scleractinia in order to improve our understanding on whether the mitogenome gene rearrangement is limited to and, therefore, could be a synapomorphy of the actual members of Caryophylliidae. Specimens of Caryophyllia scobinosa, Premocyathus sp., Heterocyathus sulcatus, and Trochocyathus caryophylloides, as well as Desmophyllum pertusum and Solenosmilia variabilis from the Southwest Atlantic were sequenced using Illumina platforms. Then, mitochondrial genomes were assembled and annotated, and nuclear datasets were recovered in-silico from assembled contigs using a previously published set of baits. Evolutionary reconstructions were performed using mitochondrial and nuclear datasets and based on Maximum Likelihood and Bayesian Inference. Obtained mitogenomes are circular and range between 15,816 and 18,225 bp in size and from 30.76% to 36.63% in GC content. The gene rearrangement is only seen in C. scobinosa, D. pertusum, Premocyathus sp., and S. variabilis, which were recovered as a monophyletic clade in both mitochondrial and nuclear phylogenies. On the other hand, the "caryophylliids" with the canonical mitogenome gene order were not recovered within this clade. Differences in features of the skeleton of "true" caryophylliids in comparison to traditional members of the family were observed and offer further support that the gene rearrangement might be seen as a synapomorphy of family Caryophylliidae.

Comparative genomic analysis of vertebrate mitochondrial reveals a differential of rearrangements rate between taxonomic class

Vertebrate mitochondrial genomes have been extensively studied for genetic and evolutionary purposes, these are normally believed to be extremely conserved, however, different cases of gene rearrangements have been reported. To verify the level of rearrangement and the mitogenome evolution, we performed a comparative genomic analysis of the 2831 vertebrate mitochondrial genomes representing 12 classes available in the NCBI database. Using a combination of bioinformatics methods, we determined there is a high number of errors in the annotation of mitochondrial genes, especially in tRNAs. We determined there is a large variation in the proportion of rearrangements per gene and per taxonomic class, with higher values observed in Actinopteri, Amphibia and Reptilia. We highlight that these are results for currently available vertebrate sequences, so an increase in sequence representativeness in some groups may alter the rearrangement rates, so in a few years it would be interesting to see if these rates are maintained or altered with the new mitogenome sequences. In addition, within each vertebrate class, different patterns in rearrangement proportion with distinct hotspots in the mitochondrial genome were found. We also determined that there are eleven convergence events in gene rearrangement, nine of which are new reports to the scientific community.

Exploring mitogenome evolution in Branchiopoda (Crustacea) lineages reveals gene order rearrangements in Cladocera

The class Branchiopoda, whose origin dates back to Cambrian, includes ~ 1200 species which mainly occupy freshwater habitats. The phylogeny and systematics of the class have been debated for long time, until recent phylogenomic analyses allowed to better clarify the relationships among major clades. Based on these data, the clade Anostraca (fairy and brine shrimps) is sister to all other branchiopods, and the Notostraca (tadpole shrimps) results as sister group to Diplostraca, which includes Laevicaudata + Spinicaudata (clam shrimps) and Cladoceromorpha (water fleas + Cyclestherida). In the present analysis, thanks to an increased taxon sampling, a complex picture emerges. Most of the analyzed mitogenomes show the Pancrustacea gene order while in several other taxa they are found rearranged. These rearrangements, though, occur unevenly among taxa, most of them being found in Cladocera, and their taxonomic distribution does not agree with the phylogeny. Our data also seems to suggest the possibility of potentially homoplastic, alternative gene order within Daphniidae.

Positive selection over the mitochondrial genome and its role in the diversification of gentoo penguins in response to adaptation in isolation

Although mitochondrial DNA has been widely used in phylogeography, evidence has emerged that factors such as climate, food availability, and environmental pressures that produce high levels of stress can exert a strong influence on mitochondrial genomes, to the point of promoting the persistence of certain genotypes in order to compensate for the metabolic requirements of the local environment. As recently discovered, the gentoo penguins (Pygoscelis papua) comprise four highly divergent lineages across their distribution spanning the Antarctic and sub-Antarctic regions. Gentoo penguins therefore represent a suitable animal model to study adaptive processes across divergent environments. Based on 62 mitogenomes that we obtained from nine locations spanning all four gentoo penguin lineages, we demonstrated lineage-specific nucleotide substitutions for various genes, but only lineage-specific amino acid replacements for the ND1 and ND5 protein-coding genes. Purifying selection (dN/dS < 1) is the main driving force in the protein-coding genes that shape the diversity of mitogenomes in gentoo penguins. Positive selection (dN/dS > 1) was mostly present in codons of the Complex I (NADH genes), supported by two different codon-based methods at the ND1 and ND4 in the most divergent lineages, the eastern gentoo penguin from Crozet and Marion Islands and the southern gentoo penguin from Antarctica respectively. Additionally, ND5 and ATP6 were under selection in the branches of the phylogeny involving all gentoo penguins except the eastern lineage. Our study suggests that local adaptation of gentoo penguins has emerged as a response to environmental variability promoting the fixation of mitochondrial haplotypes in a non-random manner. Mitogenome adaptation is thus likely to have been associated with gentoo penguin diversification across the Southern Ocean and to have promoted their survival in extreme environments such as Antarctica. Such selective processes on the mitochondrial genome may also be responsible for the discordance detected between nuclear- and mitochondrial-based phylogenies of gentoo penguin lineages.

Mitochondrial genome of the critically endangered Baer's Pochard, Aythya baeri, and its phylogenetic relationship with other Anatidae species

Historically, the diving duck, Baer’s Pochard (Aythya baeri) was widely distributed in East and South Asia, but according to a recent estimate, its global population is now less than 1000 individuals. To date, the mitochondrial genome of A. baeri has not been deposited and is not available in GenBank. Therefore, we aimed to sequence the complete mitochondrial genome of this species. The genome was 16,623 bp in length, double stranded, circular in shape, and contained 13 protein-coding genes, 22 tRNA genes, two rRNA genes, and one non-coding control region. Many structural and compositional similarities were discovered between A. baeri and the other three Aythya mitochondrial genomes. Among 13 protein-coding genes of the four Aythya species, the fastest-evolving gene was ATP8 while the slowest-evolving gene was COII. Furthermore, the phylogenetic tree of Anatidae based on Bayesian inference and maximum likelihood methods showed that the relationships among 15 genera of the Anatidae family were as follows: Dendrocygna was an early diverging lineage that was fairly distant from the other ingroup taxa; Cygnus, Branta, and Anser were clustered into one branch that corresponded to the Anserinae subfamily; and Aythya, Asarcornis, Netta, Anas, Mareca, Mergus, Lophodytes, Bucephala, Tadorna, Cairina, and Aix were clustered into another branch that corresponded to the Anatinae subfamily. Our target species and three other Aythya species formed a monophyletic group. These results provide new mitogenomic information to support further phylogenetic and taxonomic studies and genetic conservation of Anatidae species.

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.

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.

The queen conch mitogenome: intra- and interspecific mitogenomic variability in Strombidae and phylogenetic considerations within the Hypsogastropoda

Aliger gigas is an economically important and vulnerable marine species. We present a new mitogenome of A. gigas from the Mexican Caribbean and use the eight publicly available Strombidae mitogenomes to analyze intra- and interspecific variation. We present the most complete phylogenomic understanding of Hypsogastropoda to date (17 superfamilies, 39 families, 85 genera, 109 species) to revisit the phylogenetic position of the Stromboidea and evaluate divergence times throughout the phylogeny. The A. gigas mitogenome comprises 15,460 bp including 13 PCGs, 22 tRNAs, and two rRNAs. Nucleotide diversity suggested divergence between the Mexican and Colombian lineages of A. gigas. Interspecific divergence showed high differentiation among Strombidae species and demonstrated a close relationship between A. gigas and Strombus pugilis, between Lambis lambis and Harpago chiragra, and among Tridentarius dentatus/Laevistrombus canarium/Ministrombus variabilis. At the intraspecific level, the gene showing the highest differentiation is ATP8 and the lowest is NAD4L, whereas at the interspecific level the NAD genes show the highest variation and the COX genes the lowest. Phylogenomic analyses confirm that Stromboidea belongs in the non-Latrogastropoda clade and includes Xenophoridea. The phylogenomic position of other superfamilies, including those of previously uncertain affiliation, is also discussed. Finally, our data indicated that Stromboidea diverged into two principal clades in the early Cretaceous while Strombidae diversified in the Paleocene, and lineage diversification within A. gigas took place in the Pleistocene.

Complete vertebrate mitogenomes reveal widespread repeats and gene duplications

BACKGROUND

Modern sequencing technologies should make the assembly of the relatively small mitochondrial genomes an easy undertaking. However, few tools exist that address mitochondrial assembly directly.

RESULTS

As part of the Vertebrate Genomes Project (VGP) we develop mitoVGP, a fully automated pipeline for similarity-based identification of mitochondrial reads and de novo assembly of mitochondrial genomes that incorporates both long (> 10 kbp, PacBio or Nanopore) and short (100-300 bp, Illumina) reads. Our pipeline leads to successful complete mitogenome assemblies of 100 vertebrate species of the VGP. We observe that tissue type and library size selection have considerable impact on mitogenome sequencing and assembly. Comparing our assemblies to purportedly complete reference mitogenomes based on short-read sequencing, we identify errors, missing sequences, and incomplete genes in those references, particularly in repetitive regions. Our assemblies also identify novel gene region duplications. The presence of repeats and duplications in over half of the species herein assembled indicates that their occurrence is a principle of mitochondrial structure rather than an exception, shedding new light on mitochondrial genome evolution and organization.

CONCLUSIONS

Our results indicate that even in the "simple" case of vertebrate mitogenomes the completeness of many currently available reference sequences can be further improved, and caution should be exercised before claiming the complete assembly of a mitogenome, particularly from short reads alone.

Complete mitogenome of endemic plum-headed parakeet Psittacula cyanocephala - characterization and phylogenetic analysis

Psittacula cyanocephala is an endemic parakeet from the Indian sub-continent that is widespread in the illegal bird trade. Previous studies on Psittacula parakeets have highlighted taxonomic ambiguities, warranting studies to resolve the issues. Since the mitochondrial genome provides useful information concerning the species evolution and phylogenetics, we sequenced the complete mitogenome of P. cyanocephala using NGS, validated 38.86% of the mitogenome using Sanger Sequencing and compared it with other available whole mitogenomes of Psittacula. The complete mitogenome of the species was 16814 bp in length with 54.08% AT composition. P. cyanocephala mitogenome comprises of 13 protein-coding genes, 2 rRNAs and 22 tRNAs. P. cyanocephala mitogenome organization was consistent with other Psittacula mitogenomes. Comparative codon usage analysis indicated the role of natural selection on Psittacula mitogenomes. Strong purifying selection pressure was observed maximum on nad1 and nad4l genes. The mitochondrial control region of all Psittacula species displayed the ancestral avian CR gene order. Phylogenetic analyses revealed the Psittacula genus as paraphyletic nature, containing at least 4 groups of species within the same genus, suggesting its taxonomic reconsideration. Our results provide useful information for developing forensic tests to control the illegal trade of the species and scientific basis for phylogenetic revision of the genus Psittacula.

Structure and comparative analysis of the mitochondrial genomes of Liolaemus lizards with different modes of reproduction and ploidy levels

Liolaemus is the most specious genus of the Squamata lizards in South America, presenting exceptional evolutionary radiation and speciation patterns. This recent diversification complicates the formal taxonomic treatment and the phylogenetic analyses of this group, causing relationships among species to remain controversial. Here we used Next-Generation Sequencing to do a comparative analysis of the structure and organization of the complete mitochondrial genomes of three differently related species of Liolaemus and with different reproductive strategies and ploidy levels. The annotated mitochondrial genomes of ca. 17 kb are the first for the Liolaemidae family. Despite the high levels of sequence similarity among the three mitochondrial genomes over most of their lengths, the comparative analyses revealed variations at the stop codons of the protein coding genes and the structure of the tRNAs among species. The presence of a non-canonical dihydrouridine loop is a novelty for the pleurodonts iguanians. But the highest level of variability was observed in two repetitive sequences of the control region, which were responsible for most of the length heterogeneity of the mitochondrial genomes. These tandem repeats may be useful markers to analyze relationships of closely related species of Liolaemus and related genera and to conduct population and phylogenetic studies.

Quantification and evolution of mitochondrial genome rearrangement in Amphibians

BACKGROUND

Rearrangement is an important topic in the research of amphibian mitochondrial genomes ("mitogenomes" hereafter), whose causes and mechanisms remain enigmatic. Globally examining mitogenome rearrangements and uncovering their characteristics can contribute to a better understanding of mitogenome evolution.

RESULTS

Here we systematically investigated mitogenome arrangements of 232 amphibians including four newly sequenced Dicroglossidae mitogenomes. The results showed that our new sequenced mitogenomes all possessed a trnM tandem duplication, which was not exclusive to Dicroglossidae. By merging the same arrangements, the mitogenomes of ~ 80% species belonged to the four major patterns, the major two of which were typical vertebrate arrangement and typical neobatrachian arrangement. Using qMGR for calculating rearrangement frequency (RF) (%), we found that the control region (CR) (RF = 45.04) and trnL2 (RF = 38.79) were the two most frequently rearranged components. Forty-seven point eight percentage of amphibians possessed rearranged mitogenomes including all neobatrachians and their distribution was significantly clustered in the phylogenetic trees (p < 0.001). In addition, we argued that the typical neobatrachian arrangement may have appeared in the Late Jurassic according to possible occurrence time estimation.

CONCLUSION

It was the first global census of amphibian mitogenome arrangements from the perspective of quantity statistics, which helped us to systematically understand the type, distribution, frequency and phylogenetic characteristics of these rearrangements.

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.

Characterization of three new mitochondrial genomes of Coraciiformes (Megaceryle lugubris, Alcedo atthis, Halcyon smyrnensis) and insights into their phylogenetics

Coraciiformes contains more than 200 species with great differences on external morphology and life-style. The evolutionary relationships within Coraciiformes and the phylogenetic placement of Coraciiformes in Aves are still questioned. Mitochondrial genome (mitogenome) sequences are popular markers in molecular phylogenetic studies of birds. This study presented the genome characteristics of three new mitogenomes in Coraciiformes and explored the phylogenetic relationships among Coraciiformes and other five related orders with mitogenome data of 30 species. The sizes of three mitogenomes were 17,383 bp (Alcedo atthis), 17,892 bp (Halcyon smyrnensis) and 17,223 bp (Megaceryle lugubris). Each mitogenome contained one control region and 37 genes that were common in vertebrate mitogenomes. The organization of three mitogenomes was identical to the putative ancestral gene order in Aves. Among 13 available Coraciiform mitogenomes, 12 protein coding genes showed indications of negative selection, while the MT-ND6 presented sign of positive selection or relaxed purifying selection. The phylogenetic results supported that Upupidae and Bucerotidae should be separated from Coraciiformes, and that Coraciiformes is more closely related to Piciformes than to Strigiformes, Trogoniformes and Cuculiformes. Our study provide valuable data for further phylogenetic investigation of Coraciiformes.

A Highly Contiguous Genome for the Golden-Fronted Woodpecker (Melanerpes aurifrons) via Hybrid Oxford Nanopore and Short Read Assembly

Woodpeckers are found in nearly every part of the world and have been important for studies of biogeography, phylogeography, and macroecology. Woodpecker hybrid zones are often studied to understand the dynamics of introgression between bird species. Notably, woodpeckers are gaining attention for their enriched levels of transposable elements (TEs) relative to most other birds. This enrichment of TEs may have substantial effects on molecular evolution. However, comparative studies of woodpecker genomes are hindered by the fact that no high-contiguity genome exists for any woodpecker species. Using hybrid assembly methods combining long-read Oxford Nanopore and short-read Illumina sequencing data, we generated a highly contiguous genome assembly for the Golden-fronted Woodpecker (Melanerpes aurifrons). The final assembly is 1.31 Gb and comprises 441 contigs plus a full mitochondrial genome. Half of the assembly is represented by 28 contigs (contig L50), each of these contigs is at least 16 Mb in size (contig N50). High recovery (92.6%) of bird-specific BUSCO genes suggests our assembly is both relatively complete and relatively accurate. Over a quarter (25.8%) of the genome consists of repetitive elements, with 287 Mb (21.9%) of those elements assignable to the CR1 superfamily of transposable elements, the highest proportion of CR1 repeats reported for any bird genome to date. Our assembly should improve comparative studies of molecular evolution and genomics in woodpeckers and allies. Additionally, the sequencing and bioinformatic resources used to generate this assembly were relatively low-cost and should provide a direction for development of high-quality genomes for studies of animal biodiversity.

The complete mitochondrial genome of the endangered Assam Roofed Turtle, Pangshura sylhetensis (Testudines: Geoemydidae): Genomic features and phylogeny

The Assam Roofed Turtle, Pangshura sylhetensis is an endangered and least studied species endemic to India and Bangladesh. The present study decodes the first complete mitochondrial genome of P. sylhetensis (16,568 bp) by using next-generation sequencing. The assembly encodes 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and one control region (CR). Most of the genes were encoded on the majority strand, except NADH dehydrogenase subunit 6 (nad6) and eight tRNAs. All PCGs start with an ATG initiation codon, except for Cytochrome oxidase subunit 1 (cox1) and NADH dehydrogenase subunit 5 (nad5), which both start with GTG codon. The study also found the typical cloverleaf secondary structures in most of the predicted tRNA structures, except for serine (trnS1) which lacks of conventional DHU arm and loop. Both Bayesian and maximum-likelihood phylogenetic inference using 13 concatenated PCGs demonstrated strong support for the monophyly of all 52 Testudines species within their respective families and revealed Batagur trivittata as the nearest neighbor of P. sylhetensis. The mitogenomic phylogeny with other amniotes is congruent with previous research, supporting the sister relationship of Testudines and Archosaurians (birds and crocodilians). Additionally, the mitochondrial Gene Order (GO) analysis indicated plesiomorphy with the typical vertebrate GO in most of the Testudines species.

Complete mitogenomic and phylogenetic characteristics of the speckled wood-pigeon (Columba hodgsonii)

The speckled wood-pigeon, Columba hodgsonii, is mainly distributed in Bhutan, China, India, Laos, Myanmar, Nepal, Pakistan, and Thailand. Although there are several studies on birds in the family Columbidae, no study has focused on C. hodgsonii, a member of this family. Therefore, this study aimed to clarify the phylogenetic status of C. hodgsonii. The complete mitochondrial genome (mitogenome) of C. hodgsonii was sequenced and characterized and compared with those of other Columba species. The C. hodgsonii mitogenome was found to be 17,477 bp in size and contained 13 PCGs, two rRNAs, 22 tRNAs, and one CR. Of the 37 genes encoded by the C. hodgsonii mitogenome, 28 were on the heavy strand and nine were on the light strand. Twelve PCGs were initiated by ATN codons and one PCG harbored an incomplete termination codon (T-). The base composition of C. hodgsonii PCGs was A = 29.44%, T = 24.37%, G = 12.43%, and C = 33.76%. For the whole mitogenome, including PCGs, rRNAs, tRNAs, and the control region, the AT-skew was positive, and the GC-skew was negative. Phylogenetic analysis based on the base sequences of 13 PCGs from 28 Columbidae species and one outgroup using maximum likelihood and Bayesian inference indicated that C. hodgsonii belongs to the genus Columba and that the family Columbidae is monophyletic.

Mitogenome of the little owl Athene noctua and phylogenetic analysis of Strigidae

New advances in molecular approaches for DNA analysis have enhanced our understanding of the phylogenetic relationship of birds. The Little Owl (Athene noctua) is of great significance for the integrated management of forest diseases and control of regional pests. Here, we sequenced and annotated the 17,772 bp complete mitogenome of A. noctua. The mitogenome encoded 37 typical mitochondrial genes: 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA, and one non-coding control region (D-loop). The organization and location of genes in the A. noctua mitogenome were consistent with those reported for other Strigidae birds. Phylogenetic relationships based on Bayesian inference and Maximum likelihood methods showed that A. noctua has close relationships with Athene brama and Glaucidium cuculoides, confirming that A. noctua belongs to the Strigidae family. The phylogenetic relationships among seven genera of the Strigidae family used in this study were: Ninox and the other six genera were far apart, Otus and the clade ((Bubo + Strix) + Asio) were clustered into one branch, and Athene and Glaucidium were clustered into one branch. This phylogenetic classification is consistent with prior taxonomic studies on the Strigidae family. Our results provide new mitogenomic data to support further phylogenetic and taxonomic studies of Strigidae.

Mitogenomics suggests a sister relationship of Relicanthus daphneae (Cnidaria: Anthozoa: Hexacorallia: incerti ordinis) with Actiniaria

Relicanthus daphneae (formerly Boloceroides daphneae) was first described in 2006 as a giant sea anemone based on morphology. In 2014, its classification was challenged based on molecular data: using five genes, Relicanthus was resolved sister to zoanthideans, but with mixed support. To better understand the evolutionary relationship of Relicanthus with other early-branching metazoans, we present 15 newly-sequenced sea anemone mitochondrial genomes and a mitogenome-based phylogeny including all major cnidarian groups, sponges, and placozoans. Our phylogenetic reconstruction reveals a moderately supported sister relationship between Relicanthus and the Actiniaria. Morphologically, the cnidae of Relicanthus has apical flaps, the only existing synapomorphy for sea anemones. Based on both molecular and morphological results, we propose a third suborder (Helenmonae) within the Actiniaria to accommodate Relicanthus. Although Relicanthus shares the same gene order and content with other available actiniarian mitogenomes, it is clearly distinct at the nucleotide level from anemones within the existing suborders. The phylogenetic position of Relicanthus could reflect its association with the periphery of isolated hydrothermal vents, which, although patchy and ephemeral, harbor unique chemosynthetic communities that provide a relatively stable food source to higher trophic levels over long evolutionary timescales. The ability to colonize the deep sea and the periphery of new vent systems may be facilitated by Relicanthus’ large and extremely yolky eggs.

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.

Homoplasy or plesiomorphy? Reconstruction of the evolutionary history of mitochondrial gene order rearrangements in the subphylum Neodermata

Recent mitogenomic studies have exposed a gene order (GO) shared by two classes, four orders and 31 species ('common GO') within the flatworm subphylum Neodermata. There are two possible hypotheses for this phenomenon: convergent evolution (homoplasy) or shared ancestry (plesiomorphy). To test those, we conducted a meta-analysis on all available mitogenomes to infer the evolutionary history of GO in Neodermata. To improve the resolution, we added a newly sequenced mitogenome that exhibited the common GO, Euryhaliotrema johni (Ancyrocephalinae), to the dataset. Phylogenetic analyses conducted on two datasets (nucleotides of all 36 genes and amino acid sequences of 12 protein coding genes) and four algorithms (MrBayes, RAxML, IQ-TREE and PhyloBayes) produced topology instability towards the tips, so ancestral GO reconstructions were conducted using TreeREx and MLGO programs using all eight obtained topologies, plus three unique topologies from previous studies. The results consistently supported the second hypothesis, resolving the common GO as a plesiomorphic ancestral GO for Neodermata, Cestoda, Monopisthocotylea, Cestoda + Trematoda and Cestoda + Trematoda + Monopisthocotylea. This allowed us to trace the evolutionary GO scenarios from each common ancestor to its descendants amongst the Monogenea and Cestoda classes, and propose that the common GO was most likely retained throughout all of the common ancestors, leading to the extant species possessing the common GO. Neodermatan phylogeny inferred from GOs was largely incongruent with all 11 topologies described above, but it did support the mitogenomic dataset in resolving Polyopisthocotylea as the earliest neodermatan branch. Although highly derived GOs might be of some use in resolving isolated taxonomic and phylogenetic uncertainties, we conclude that, due to the discontinuous nature of their evolution, they tend to produce artefactual phylogenetic relationships, which makes them unsuitable for phylogenetic reconstruction in Neodermata. Wider and denser sampling of neodermatan mitogenomic sequences will be needed to infer the evolutionary pathways leading to the observed diversity of GOs with confidence.

Complete mitochondrial genome sequence of the Himalayan Griffon, Gyps himalayensis (Accipitriformes: Accipitridae): Sequence, structure, and phylogenetic analyses

This is the first study to describe the mitochondrial genome of the Himalayan Griffon, Gyps himalayensis, which is an Old World vulture belonging to the family Accipitridae and occurring along the Himalayas and the adjoining Tibetan Plateau. Its mitogenome is a closed circular molecule 17,381 bp in size containing 13 protein-coding genes, 22 tRNA coding genes, two rRNA-coding genes, a control region (CR), and an extra pseudo-control region (CCR) that are conserved in most Accipitridae mitogenomes. The overall base composition of the G. himalayensis mitogenome is 24.55% A, 29.49% T, 31.59% C, and 14.37% G, which is typical for bird mitochondrial genomes. The alignment of the Accipitridae species control regions showed high levels of genetic variation and abundant AT content. At the 5' end of the domain I region, a long continuous poly-C sequence was found. Two tandem repeats were found in the pseudo-control regions. Phylogenetic analysis with Bayesian inference and maximum likelihood based on 13 protein-coding genes indicated that the relationships at the family level were (Falconidae + (Cathartidae + (Sagittariidae + (Accipitridae + Pandionidae))). In the Accipitridae clade, G. himalayensis is more closely related to Aegypius monachus than to Spilornis cheela. The complete mitogenome of G. himalayensis provides a potentially useful resource for further exploration of the taxonomic status and phylogenetic history of Gyps species.

Comparative mitogenomics of the Decapoda reveals evolutionary heterogeneity in architecture and composition

The emergence of cost-effective and rapid sequencing approaches has resulted in an exponential rise in the number of mitogenomes on public databases in recent years, providing greater opportunity for undertaking large-scale comparative genomic and systematic research. Nonetheless, current datasets predominately come from small and disconnected studies on a limited number of related species, introducing sampling biases and impeding research of broad taxonomic relevance. This study contributes 21 crustacean mitogenomes from several under-represented decapod infraorders including Polychelida and Stenopodidea, which are used in combination with 225 mitogenomes available on NCBI to investigate decapod mitogenome diversity and phylogeny. An overview of mitochondrial gene orders (MGOs) reveals a high level of genomic variability within the Decapoda, with a large number of MGOs deviating from the ancestral arthropod ground pattern and unevenly distributed among infraorders. Despite the substantial morphological and ecological variation among decapods, there was limited evidence for correlations between gene rearrangement events and species ecology or lineage specific nucleotide substitution rates. Within a phylogenetic context, predicted scenarios of rearrangements show some MGOs to be informative synapomorphies for some taxonomic groups providing strong independent support for phylogenetic relationships. Additional comparisons for a range of mitogenomic features including nucleotide composition, strand asymmetry, unassigned regions and codon usage indicate several clade-specific trends that are of evolutionary and ecological interest.

Unravelling population processes over the Late Pleistocene driving contemporary genetic divergence in Palearctic buzzards

Population range expansions and contractions as a response to climate and habitat change throughout the Quaternary are known to have contributed to complex phylogenetic and population genetic events. Speciation patterns and processes in Palearctic buzzards (genus Buteo) are a long-standing example of morphological and genetic data incongruence, attributed to panmixia, habitat range shifts, contact zones, and climate change. Here we assess the systematics, phylogeography and population genetic structure of three nominal species of Palearctic buzzards, Buteo buteo (including B. b. vulpinus), B. rufinus (including B. r. cirtensis) and B. hemilasius. Phylogenetic analyses inferred from mitochondrial data recover B. hemilasius as sister to the sister clades B. r. rufinus and B. buteo complex (B. b. buteo, B. b. vulpinus, but also including B. r. cirtensis). In contrast, we find an unresolved genetic delimitation inferred from four nuclear loci, suggesting an ancestral genetic pool for all species. Time-trees suggest population contractions and expansions throughout the Pleistocene, which likely reflect habitat change and contrasting ecological niche requirements between species. Microsatellite-based extended Bayesian skyline plots reveal relatively constant population sizes for B. hemilasius, B. r. rufinus, and B. b. vulpinus, in contrast to a dramatic population expansion in B. r. cirtensis within the last 3 kya. Overall, our study illustrates how complex population processes over the Late Pleistocene have shaped the patterns of genetic divergence in Palearctic buzzards, due to the joint effects of shared ancestral polymorphisms, population expansions and contractions, with hybridization at contact zones leading to admixture and introgression.

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.

Evolutionary progression of mitochondrial gene rearrangements and phylogenetic relationships in Strigidae (Strigiformes)

The bird mitogenome is generally considered to have a conservative genome size, consistent gene content, and similar gene order. As more mitogenomes are sequenced, mitochondrial (mt) gene rearrangements have been frequently identified among diverse birds. Within two genera (Bubo and Strix) of typical owls (Strigidae, Strigiformes), the rearrangement of the mt gene has been a subject of debate. In the current study, we first sequenced the whole mitogenomes of S. uralensis and B. scandiaca and resequenced the entire mitogenome of B. bubo. By combining our data with previously sequenced mitogenomes in Strigidae, we examined the mt gene rearrangements in the family and attempted to reconstruct the evolutionary progression of these rearrangements. The mitogenomes were then used to review the phylogenies of Strigidae. Most mitogenomes exhibited the ancestral gene order (A) in Strigidae. The ancestral gene order in the previously published mitogenome of B. bubo was found to be incorrect. We determined the mt gene order (the duplicate tRNA^Thr-CR, B) and discovered two additional mt gene orders (the duplicate tRNA^Glu-L-CR and CR, C and D) in the Bubo and Strix genera. Gene order B was likely derived from A by a tandem duplication of the region spanning from tRNA^Thr to CR. The other two modified gene orders, C and D, were likely derived from B by further degenerations or deletions of one copy of specific duplicated genes. We also preliminarily reconstructed the evolutionary progression of mt gene rearrangements and discussed maintenance of the duplicated CR in the genera. Additionally, the phylogenetic trees based on the mitogenomes supported the division of Strigidae into three subfamilies: Ninoxinae + (Surniinae + Striginae). Within the Striginae clade, the four genera formed a phylogenetic relationship: Otus + (Asio + (Bubo + Strix)). This suggests that Otus firstly diverges in their evolutionary history, and Bubo and Strix show a close relationship. B. bubo, B. blakistoni and B. scandiaca form a clade should be considered members of the same genus. The well-supported topology obtained in our Bayesian inference (BI) and maximum likelihood (ML) analyses of Strigid mitogenomes suggests that these genomes are informative for constructing phylogenetic relationships.

The complete mitochondrial genomes of Tarsiger cyanurus and Phoenicurus auroreus: a phylogenetic analysis of Passeriformes

Passeriformes is the largest group within aves and the phylogenetic relationships between Passeriformes have caused major disagreement in ornithology. Particularly, the phylogenetic relationships between muscicapoidea and sylvioidea are complex, and their taxonomic boundaries have not been clearly defined. Our aim was to study the status of two bird species: Tarsiger cyanurus and Phoenicurus auroreus. Furthermore, we analyzed the phylogenetic relationships of Passeriformes. Complete mitochondrial DNA (mtDNA) sequences of both species were determined and the lengths were 16,803 (T. cyanurus) and 16,772 bp (P. auroreus), respectively. Thirteen protein-coding genes, 22 tRNA genes, two rRNA genes, and one control region were identified in these mtDNAs. The contents of A and T at the base compositions was significantly higher than the content of G and C, and this AT skew was positive, while the GC skew was negative. The monophyly of Passeriformes is divided into four major clades: Corvoidea, Sylvioidea, Passeroidea, and Musicicapoidea. Paridae should be separated from the superfamily Sylvioidea and placed within the superfamily Muscicapoidea. The family Muscicapidae and Corvida were paraphyly, while Carduelis and Emberiza were grouped as a sister taxon. The relationships between some species of the order passeriformes may remain difficult to resolve despite an effort to collect additional characters for phylogenetic analysis. Current research of avian phylogeny should focus on adding characters and taxa and use both effectively to obtain a better resolution for deeper and shallow nodes.

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.

Mitogenomes of two neotropical bird species and the multiple independent origin of mitochondrial gene orders in Passeriformes

At least four mitogenome arrangements occur in Passeriformes and differences among them are derived from an initial tandem duplication involving a segment containing the control region (CR), followed by loss or reduction of some parts of this segment. However, it is still unclear how often duplication events have occurred in this bird order. In this study, the mitogenomes from two species of Neotropical passerines (Sicalis olivascens and Lepidocolaptes angustirostris) with different gene arrangements were first determined. We also estimated how often duplication events occurred in Passeriformes and if the two CR copies demonstrate a pattern of concerted evolution in Sylvioidea. One tissue sample for each species was used to obtain the mitogenomes as a byproduct using next generation sequencing. The evolutionary history of mitogenome rearrangements was reconstructed mapping these characters onto a mitogenome Bayesian phylogenetic tree of Passeriformes. Finally, we performed a Bayesian analysis for both CRs from some Sylvioidea species in order to evaluate the evolutionary process involving these two copies. Both mitogenomes described comprise 2 rRNAs, 22 tRNAs, 13 protein-codon genes and the CR. However, S. olivascens has 16,768 bp showing the ancestral avian arrangement, while L. angustirostris has 16,973 bp and the remnant CR2 arrangement. Both species showed the expected gene order compared to their closest relatives. The ancestral state reconstruction suggesting at least six independent duplication events followed by partial deletions or loss of one copy in some lineages. Our results also provide evidence that both CRs in some Sylvioidea species seem to be maintained in an apparently functional state, perhaps by concerted evolution, and that this mechanism may be important for the evolution of the bird mitogenome.

Complete mitochondrial genome sequences of the northern spotted owl (Strix occidentalis caurina) and the barred owl (Strix varia; Aves: Strigiformes: Strigidae) confirm the presence of a duplicated control region

We report here the successful assembly of the complete mitochondrial genomes of the northern spotted owl (Strix occidentalis caurina) and the barred owl (S. varia). We utilized sequence data from two sequencing methodologies, Illumina paired-end sequence data with insert lengths ranging from approximately 250 nucleotides (nt) to 9,600 nt and read lengths from 100–375 nt and Sanger-derived sequences. We employed multiple assemblers and alignment methods to generate the final assemblies. The circular genomes of S. o. caurina and S. varia are comprised of 19,948 nt and 18,975 nt, respectively. Both code for two rRNAs, twenty-two tRNAs, and thirteen polypeptides. They both have duplicated control region sequences with complex repeat structures. We were not able to assemble the control regions solely using Illumina paired-end sequence data. By fully spanning the control regions, Sanger-derived sequences enabled accurate and complete assembly of these mitochondrial genomes. These are the first complete mitochondrial genome sequences of owls (Aves: Strigiformes) possessing duplicated control regions. We searched the nuclear genome of S. o. caurina for copies of mitochondrial genes and found at least nine separate stretches of nuclear copies of gene sequences originating in the mitochondrial genome (Numts). The Numts ranged from 226–19,522 nt in length and included copies of all mitochondrial genes except tRNA^Pro, ND6, and tRNA^Glu. Strix occidentalis caurina and S. varia exhibited an average of 10.74% (8.68% uncorrected p-distance) divergence across the non-tRNA mitochondrial genes.

Northern Spotted Owl (Strix occidentalis caurina) Genome: Divergence with the Barred Owl (Strix varia) and Characterization of Light-Associated Genes

We report here the assembly of a northern spotted owl (Strix occidentalis caurina) genome. We generated Illumina paired-end sequence data at 90× coverage using nine libraries with insert lengths ranging from ∼250 to 9,600 nt and read lengths from 100 to 375 nt. The genome assembly is comprised of 8,108 scaffolds totaling 1.26 × 109 nt in length with an N50 length of 3.98 × 106 nt. We calculated the genome-wide fixation index (FST) of S. o. caurina with the closely related barred owl (Strix varia) as 0.819. We examined 19 genes that encode proteins with light-dependent functions in our genome assembly as well as in that of the barn owl (Tyto alba). We present genomic evidence for loss of three of these in S. o. caurina and four in T. alba. We suggest that most light-associated gene functions have been maintained in owls and their loss has not proceeded to the same extent as in other dim-light-adapted vertebrates.

Rapid and recent diversification patterns in Anseriformes birds: Inferred from molecular phylogeny and diversification analyses

The Anseriformes is a well-known and widely distributed bird order, with more than 150 species in the world. This paper aims to revise the classification, determine the phylogenetic relationships and diversification patterns in Anseriformes by exploring the Cyt b, ND2, COI genes and the complete mitochondrial genomes (mito-genomes). Molecular phylogeny and genetic distance analyses suggest that the Dendrocygna species should be considered as an independent family, Dendrocygnidae, rather than a member of Anatidae. Molecular timescale analyses suggests that the ancestral diversification occurred during the Early Eocene Climatic Optimum (58 ~ 50 Ma). Furthermore, diversification analyses showed that, after a long period of constant diversification, the median initial speciation rate was accelerated three times, and finally increased to approximately 0.3 sp/My. In the present study, both molecular phylogeny and diversification analyses results support that Anseriformes birds underwent rapid and recent diversification in their evolutionary history, especially in modern ducks, which show extreme diversification during the Plio-Pleistocene (~ 5.3 Ma). Therefore, our study support that the Plio-Pleistocene climate fluctuations are likely to have played a significant role in promoting the recent diversification for Anseriformes.

Complete mitochondrial genome of white-backed woodpecker Dendrocopos leucotos (Piciformes: Picidae) and its phylogenetic position

The white-backed woodpecker (Dendrocopos leucotos) is an ecologically important bird in Eurasian forest ecosystems. In this study, the complete mitochondrial genome of the species was sequenced using next-generation sequencing technology. The genome was 16,838 bp in length, consisting of 13 protein coding genes, two rRNAs, 22 tRNAs, a non-coding control region and a repeat region. Phylogenetic analysis using available complete mitochondrial genomes of the Coraciimorphae supported the monophyly of the Piciformes. The complete mitochondrial genome of D. leucotos will be a useful genetic resource for population genetics, phylogenetic analysis and conservation of the species.

Striking pseudogenization in avian phylogenetics: Numts are large and common in falcons

Nuclear copies of mitochondrial genes (numts) are a well-known feature of eukaryotic genomes and a concern in systematics, as they can mislead phylogenetic inferences when inadvertently used. Studies on avian numts initially based on the chicken genome suggest that numts may be uncommon and relatively short among birds. Here we ask how common numts are in falcons, based on recently sequenced genomes of the Saker falcon (Falco cherrug) and Peregrine falcon (F. peregrinus). We identified numts by BLASTN searches and then extracted CYTB, ND2 and COI sequences from them, which were then used for phylogeny inference along with several sequences from other species in Falconiformes. Our results indicate that avian numts may be much more frequent and longer than previously thought. Phylogenetic inferences revealed multiple independent nuclear insertions throughout the history of the Falconiformes, including cases of sequences available in public databases and wrongly identified as authentic mtDNA. New sequencing technologies and ongoing efforts for whole genome sequencing will provide exciting opportunities for avian numt research in the near future.

Intraspecific rearrangement of mitochondrial genome suggests the prevalence of the tandem duplication-random loss (TDLR) mechanism in Quasipaa boulengeri

BACKGROUND

Tandem duplication followed by random loss (TDRL) is the most frequently invoked model to explain the diversity of gene rearrangements in metazoan mitogenomes. The initial stages of gene rearrangement are difficult to observe in nature, which limits our understanding of incipient duplication events and the subsequent process of random loss. Intraspecific gene reorganizations may represent intermediate states, and if so they potentially shed light on the evolutionary dynamics of TDRL.

RESULTS

Nucleotide sequences in a hotspot of gene-rearrangement in 28 populations of a single species of frog, Quasipaa boulengeri, provide such predicted intermediate states. Gene order and phylogenetic analyses support a single tandem duplication event and a step-by-step process of random loss. Intraspecific gene rearrangements are not commonly found through comparison of all mitochondrial DNA records of amphibians and squamate reptiles in GenBank.

CONCLUSIONS

The intraspecific variation in Q. boulengeri provides insights into the rate of partial duplications and deletions within a mitogenome, and reveals that fixation and gene-distribution in mitogenomic reorganization is likely non-adaptive.

Two complete mitochondrial genomes from Praticolella mexicana Perez, 2011 (Polygyridae) and gene order evolution in Helicoidea (Mollusca, Gastropoda)

Helicoidea is a diverse group of land snails with a global distribution. While much is known regarding the relationships of helicoid taxa, comparatively little is known about the evolution of the mitochondrial genome in the superfamily. We sequenced two complete mitochondrial genomes from Praticolellamexicana Perez, 2011 representing the first such data from the helicoid family Polygyridae, and used them in an evolutionary analysis of mitogenomic gene order. We found the mitochondrial genome of Praticolellamexicana to be 14,008 bp in size, possessing the typical 37 metazoan genes. Multiple alternate stop codons are used, as are incomplete stop codons. Mitogenome size and nucleotide content is consistent with other helicoid species. Our analysis of gene order suggested that Helicoidea has undergone four mitochondrial rearrangements in the past. Two rearrangements were limited to tRNA genes only, and two involved protein coding genes.